Contemporary Engineering Sciences, Vol. 11, 2018, no. 5, 237-247 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.712204 Hydrodynamic Control of the Growth of the Two Dimensional Boundary Layers around a Flat Plate Placed in the Center of Convergent: Concept of Uniform Accessibility from the Point of View of the Diffusion Matter Kabouche Nabil University Batna 1, Laboratory (LEMPAU) of the University of Batna 1, Algeria Copyright 2018 Kabouche Nabil. This article is distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium provided the original work is properly cited. Abstract The authors indicated numerically, that a plane surface placed in the center of convergent and crossed by a Newtonian fluid, in laminar and permanent flow is uniformly accessible from certain distance of the entering edge when the velocity at the outer boundary of the boundary layer that develops on its increases linearly with the longitudinal component x. They explained the concept of hydrodynamic and mass uniform accessibility and set its validity area according to the Reynolds number. Such study can assist in the understanding of vapor phase deposition phenomena of semiconductors. Keywords: forced convection, vapor deposition, uniformly accessible surface, mass transfer, semiconductors
238 Kabouche Nabil 1. Introduction The control of the boundary layer finds numerous, different applications in the industrial sector. One of the applications of this process is the technique of organometallic vapor phase epitaxy (OMVPE) [1, 2]. In fact, in this technique, it is sought to standarize on a fairly large extent the deposition thicknesses of metal compounds heated to high temperature. For example, let us consider the deposit of the metallic compounds obtained from gaseous mixture consisting essentially of hydrogen acting as a carrier gas with in low concentrations of trimethylgallium, trimethylalimunium and trimethyllantimony. The gas mixture, which should be extracted at ordinary temperature of 25, is pushed inside the channel. If the conditions are well chosen, in immediate contact with the plane plate carried out in a temperature of 800 C, the organometallic gases break down and the metal atoms are deposited on the plate to form the desired solid compound. In reality, the reaction mechanism is poorly known and it is assumed that above the plane plate the following chemical reaction occurs [3, 4] σ Al + (1 σ) Ga + 1 4 Sb 4 Ga 1 σ Al σ Sb (1) Through this study, we ask whether it is possible to control the flux of the material diffusion on the plane plate by imposing a velocity on the external border of the boundary layer which develops on the latter, a law of variation with the x coordinate in the following form: U e (x) = C x. In other words we ask whether it is possible to control the thicknesses of the metallic compounds deposited on the substrate by imposing a variation law of the velocity at the external border of the boundary layer which develops on the substrate. The theoretical works of [5] indicated that it is possible to control the material diffusion flux on the plate provided that it is placed in the center of a channel with a hyperbolic walls, whose wall profile is given by : r 2 sin[2θ] = Cte which results in linear increase of the longitudinal component of the velocity according to direction of the Ox. In this case, the velocity profile of the border in the boundary layer that develops on the plane plate is given by [6,7] U e (x) = C x. The results presented by the cited authors [5] are insufficient to explain the concept of the uniform accessibility. Indeed, they do not specify the conditions of the uniform accessibility according to all parameters of the studied geometric, namely: (i) the link between the uniform accessibility and the geometric parameters of the studied system, (ii) the relationship between the uniform accessibility of the plane surface and the Reynolds number.
Hydrodynamic control of the growth 239 y H 2 Z O x Figure 1 : Schematic representation of the duct and of the co-ordinates 2. Mathematical formulation of the problem The geometry of the studied problem is presented in figure 1. It is a flat plate placed in the center of a channel and a Cartesian reference (O, x, y, z). Forcing a Newtonian fluid in the case of hydrogen to circulate in laminar and permanent flow inside the channel follows the Ox direction. Adding the following simplifying assumptions: (i) the viscous friction and the radiation transfers are negligible, (ii) the mass fractions are constant and uniform, (iii) the diffusion coefficients of the reagents are equals. 2.1 Continuity equation: U x + V y = 0 (2) 2.2 The conservation equation of the movement quantity: (U U U + V x y ) = U e du e dx 2.3 Diffusion equation of matter: + 1 Re 2 U (3) y2 (U C k x + C V k y ) = 1 2 C k (4) Pem y2 2.4 Boundary conditions: - On the surface of the substrate (y = 0) U = 0; V = 0; C k = 1
240 Kabouche Nabil 3. Numerical processing - At the external border of the boundary layer U U e ; C k 0 The transport equations, after being discretized, can all be put in the following unique form: k+1 k i,j = (1 β) i,j (5) +β [ (a k i+1,j i+1,j k + a i,j+1 i,j+1 + a i 1,j k+1 i 1,j + a i,j 1 k+1 i,j 1 ) + ( x i y S ) ] (a i+1,j + a i 1,j + a i,j+1 + a i,j 1 ) The are the dimensionless unknowns composed of quantities (C k, U, V ), the a i,j are given as functions of the network steps, the dimensions of the controlled volume and the physical constants of the problem. They are calculated using the numerical scheme PLDS (Power low Differencing Scheme) of Patankar [8, 9]. 4. Results and discussion 4.1 validation of the calculation code After validating our calculation code, we impose a profile of velocity U e (x ) = 1 at the border of the boundary layer which develops on the flat plate. Our calculation code gives us the results shown in figure 2; it is well known in the literature because it is the Blasius flow [10]. 4.2 Flow study: Definition of hydrodynamic boundary layer and hydrodynamic entry distance: Concept of uniform accessibility from hydrodynamic point of view Figure 3 shows the variation of the longitudinal component of velocity U in terms of the dimensionless ordinate y for different values of the dimensionless abscissa x and the Reynold number Re. In general, when y increases, the component U differ from 0 value on the flate plate to value U max independent of y. It is also possible to define, on the flat plate, a hydrodynamic boundary layer with thickness δ H equal to the value of y from which U = U max.
Hydrodynamic control of the growth 241 In figure 4 which illustrates the variations of δ H in terms of x for Re = 1500, we see in particular that from certain hydrodynamic distance x eh, the thickness δ H does not depending on x anymore and takes a maximum value δ Hmax : then we can say that the flat plate is uniformly accessible from the hydrodynamic point of view. Table 1 summarizes the values of δ Hmax and x eh for different Reynolds number values. Re δ Hmax x eh 100 0,053 0,222 300 0,041 0,395 500 0,038 0,405 800 0,030 0,438 1500 0,023 0,401 2000 0,016 0,533 Table 1: Values of x eh and δ Hmax as a function of Reynolds The local values variations of the coefficient of friction represented in figure 5 reaffirm that for x > x eh the longitudinal component of the velocity increases linearly with the longitudinal component x. 4.3 Study of mass transfer: Definition of the mass boundary layer and the mass entry distance: Concept of uniform accessibility from the point of view of matter diffusion The curves represent the variations of C (k) depending on y, for different values of x and the Reynolds number (figure 6). In general terms, when y increases the mass fraction C (k) decreases, it goes from a value equal to the unit away from the plate to a zero value on the surface of the flat plate. There is a critical value of x, defining a mass input distance x > x ec, from which C (k) depends more on x. This result is confirmed by the figure which represents the variations of Sherwood number according to x. Thus, for x > x ec, it can be concluded that the flat plate is uniformly accessible from the matter diffusion point of view (the thickness of the metal compound which is deposited on the substrate is constant) The entry distance from which uniform accessibility is observed depends on the Reynolds number and it gets even smaller when the Reynolds number decreases. Table 2 gives the values of x ec according to the Reynolds number.
242 Kabouche Nabil 5. Conclusion Re x ec 100 0,231 300 0,287 500 0,325 800 0,377 1500 0,403 2000 0,581 Table 2: Values of x ec, as a function of Reynolds. At the end of this numerical study we can say that a flat surface is uniformly accessible from hydrodynamic and mass point of view from certain distance of the input edge when the velocity at the border of the boundary layer that develops on it is linear increases with the dimensionless abscissax. Such case is obtained when the plane surface is placed in the center of the channel whose wall profile is an equilateral hyperbole with equation r 2 sin[2θ] = Cte. The entry distance from which the uniform accessibility is observed depends on the Reynolds number and it gets even smaller when the Reynolds number decreases. 1,0 0,8 U* 0,6 0,4 Reynolds=800 Reynolds=1500 0,2 0,0 0,00 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40 y* Figure 2: Dimensionless velocity U *versus the dimensionless normal coordinate y for different values of Reynolds
Hydrodynamic control of the growth 243 U* 2,2 2,0 1,8 1,6 1,4 1,2 1,0 0,8 0,6 0,4 0,2 Reynolds=500 x*=0,01 x*=0,8 x*=0,5 0,0 0,00 0,01 0,02 0,03 0,04 0,05 y* U* 2,2 2,0 1,8 1,6 1,4 1,2 1,0 0,8 0,6 0,4 0,2 Reynolds=1500 x*=0,8 x*=0,5 x*=0,01 0,0 0,00 0,01 0,02 0,03 0,04 y* Figure. 3. Variations of the dimensionless longitudinal velocity U versus the dimensionless normal co-ordinate y for different values of x and Reynolds Thickness of the boundary layer 0,023 0,022 0,021 0,020 0,019 0,018 0,017 0,016 Re=1500 0,015 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 X eh x* Figure. 4. Variations of δ H versus x for Re =1500
244 Kabouche Nabil 1,0 Reynolds=500 1,0 Reynolds=1500 0,8 0,8 C*(k) 0,6 0,4 x*=0,5; 0,7; 1 C*(k) 0,6 0,4 x*=0,5; 0,7; 1 0,2 x*=0,01 0,2 x*=0,01 0,0 0,00 0,01 0,02 0,03 0,04 y* 0,0 0,00 0,01 0,02 0,03 0,04 y* Figure. 5. Variations of the concentration of the diffusing substance C (k) versus the dimensionless normal co-ordinate y for different values of x and Reynolds 1200 1000 800 Cf 600 RE=1500 Re=800 Re=500 400 200 0 0,0 0,2 0,4 0,6 0,8 1,0 x* Figure. 6. Variations of the local sheer stress versus the dimensionless normal coordinate x for different values of Reynolds
Hydrodynamic control of the growth 245 Figure. 7. Variations of the Sherwood numbers versus the dimensionless normal co-ordinate x for different values of Reynolds Nomenclature. Sher 1000 900 800 700 600 500 400 300 200 100 Re=1500 Re=800 Re=500 0 0,0 0,2 0,4 0,6 0,8 1,0 Al Chemical symbol of aluminum. Cf : Coefficient of local friction.(cf = ( U / y )) Ga : Chemical symbol of Gallium L : Length of channel. Pem : Number of Peclet mass (Pem = Re. Sc) Sb : Chemical symbol of antimoine. Sher : Number of local Sherwood. Sher = [ C n y ] y =0 Sc : Number of Schmidt.(Sc = θ/d n ) Re : Number of Reynolds. C k : Mass fraction of metallic gas.. C k : Dimensionless mass fraction of metallic gas. U e : Speed of the gas mixture at the inlet of the duct. U : Component of the following speed x. U : Dimensionless component of the following velocity x. (U = U/U e ) V : Component of the following velocity y V : Dimensionless component of the following velocity y.(v = V/U e ) x : Cartesian coordinate defined in Fig. x : Dimensionless coordinate following x (x = x/l) x ec : Dimensionless mass input distance. y : Cartesian coordinate defined in figure 1 y : Dimensionless coordinate following y (y = y/l) x*
246 Kabouche Nabil Greek letters v: Kinematic viscosity of hydrogen β :Coefficient of over-relaxation. : Unknown dimensionless (C k, U, V ) σ : Constant that express the number of moles in chemical reaction References [1] Gerald B. Stringfellow, Organometallic Vapor - Phase Epitaxy: Theory and Practice, 2st Edition, Academic Press 1999. [2] K. Y. Ma, Z. M. Fang, R. M. Cohen and G. B. Stringfellow, Organometallic vapor- phase epitaxy growth and characterization of Bi- containing III/V alloys, Journal of Applied Physics, 68 (1990), 4586-4591. https://doi.org/10.1063/1.346166 [3] A. Ali Chérif, CH.R.R. Raminosoa, M. Rakotamalala, A. Daif, M. Daguenet, Contrôle hydrodynamique des couches limites thermiques en convection mixte autour d'ellipsoïdes aplatis axisymétriques, International Journal of Heat Mass Transfer, 40 (1997), no. 3, 657-663. https://doi.org/10.1016/0017-9310(96)00077-4 [4] A. Ali Chérif, M. Rakotamalala, A. Daif and M. Daguenet, Contro le hydrodynamique en convection mixte de l'épaisseur de dépo ts en phase gazeuse de semi-conducteurs sur des corps à symétrie de révolution, Can. J. Chem. Eng., 73 (1995), no. 6, 908-917. https://doi.org/10.1002/cjce.5450730614 [5] N. Kabouche, A. Nahoui, Semi analytical study of a flat surface uniformly accessible in a permanent laminar flow and a Newtonian fluid, Adv. Theor. Appl. Mech., 1 (2008), no. 3, 173-186. [6] A. Pantokratoras, The Falkner-Skan flow with constant wall temperature and variable viscosity, International Journal of Thermal Sciences, 45 (2006), no. 4, 378-389. https://doi.org/10.1016/j.ijthermalsci.2005.06.004 [7] N. Kabouche, A. Benkhaled, C.M bow, M. Daguenet, Etude Numérique d une surface plane uniformément accessible dans un écoulement permanent et laminaire d un fluide newtonien, International Journal of Thermal Sciences, 41 (2002), no. 5, 455-469. https://doi.org/10.1016/s1290-0729(02)01338-8
Hydrodynamic control of the growth 247 [8] S.V. Patankar, Numerical Heat Transfer and Fluid Flow, McGraw-Hill, New York, 1980. [9] H.K. Versteeg and W. Malalasekera, An Introduction to Computational Fluid Dynamics, The finite volume method, Pearson Prentice Hall, 1995. [10] Couche limite laminaire et Turbulence et couche limite, J. Cousteix- Cepadues - Editions, Toulouse (1989). Received: December 23, 2017; Published: February 27, 2018