Modelling and Simulation in Engineering
Volume 2016 (2016), Article ID 7619746, 9 pages
http://dx.doi.org/10.1155/2016/7619746

Research Article

Extended Macroscopic Study of Dilute Gas Flow within a Microcavity

Mohamed Hssikou, Jamal Baliti, and Mohammed Alaoui

Département de physique, Université Moulay Ismaïl, Meknès, Morocco

Received 27 August 2016; Accepted 24 October 2016

Academic Editor: Ricardo Perera

Copyright © 2016 Mohamed Hssikou et al. This is an open access article 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 behaviour of monatomic and dilute gas is studied in the slip and early transition regimes using the extended macroscopic theory. The gas is confined within a two-dimensional microcavity where the longitudinal sides are in the opposite motion with constant velocity . The microcavity walls are kept at the uniform and reference temperature . Thus, the gas flow is transported only by the shear stress induced by the motion of upper and lower walls. From the macroscopic point of view, the regularized 13-moment equations of Grad, R13, are solved numerically. The macroscopic gas proprieties are studied for different values of the so-called Knudsen number (Kn), which gives the gas-rarefaction degree. The results are compared with those obtained using the classical continuum theory of Navier-Stokes and Fourier (NSF).

1. Introduction

Recently, the technology of the Microelectromechanical Systems (MEMS) has greatly developed and they have wide areas of application [1–3]. This fast growth of MEMS use is not followed enough by the physical understanding of rarefied gas flows in these microdevices. For this purpose, several studies have been recently focused on for more understanding of the physical phenomena involved in these small devices [4]. In fact, the performances of MEMS often defy the predictions made using the scaling laws developed for large systems. In fact, the gas flows inside the MEMS, under the standard conditions, are usually characterized by a mean free path comparable to the system characteristic-length . Thus, the so-called Knudsen number of gas flow is in the slip-transition regimes range; that is, . In this case, the conventional computational fluid dynamics (CFD) scheme, based on the classical Navier-Stokes and Fourier (NSF) equations, becomes inappropriate to describe the gas flow behaviour in MEMS devices. Therefore, the Knudsen number, in MEMS, is not sufficiently small to guarantee the validity of the NSF equations and the processes in MEMS need to be modelled with more accurate transport models. Similar rarefaction effects can be found in the problems of gas flows under low pressure and atmospheric conditions [5]. For gas flows outside the hydrodynamic regime ( [6], many interesting rarefaction effects such as velocity-slip and temperature-jump at the walls [7–10], Knudsen paradox, Knudsen layers [11], transpiration flow [12, 13], thermal stress [14], and heat flux without temperature gradients can take place [15]. Hence, there is a pressing need to develop the more accurate methods allowing a good description of gas-dynamic processes into these microsystems. The direct simulation Monte Carlo (DSMC) is the largely kinetic method used to simulate a rarefied gas flow where the behaviour is mainly described by the Boltzmann equation [16]. The accuracy of this method is proved by many previous studies especially with the actual computers capabilities. But, the computational cost and fluctuations noises, especially in the low-signal flows, remain the major inconveniences of this kinetic method [17]. Indeed, many macroscopic approaches are proposed such as the Chapman-Enskog (CE) expansion and the Grad moments theory. At the first order of CE both approaches lead to the famous laws of Navier-Stokes and Fourier. However, on one hand, the instability of Burnett equations obtained at second order of CE expansion is the main problem of this approach. On the other hand, the Grad 13-moment equations are hyperbolic in nature, yielding finite wave speeds, and discontinuous subshock structures when the Mach number lies above [18]. We also note that the Grad 13-moment equations, for nonlinear problems, lack suitable boundary conditions. But, when the rarefaction degree becomes more intense, the Knudsen number value is in the range of , and a purely kinetic approach is needed for describing the gas flows [19].

Keeping the benefits of both approaches and to avoid their failures, Struchtrup has adopted recently the combination of the above approaches [19]. This leads to the set of regularized 13-moments equations (R13) used as higher order of continuum solution, that is, , to capture the rarefaction effects, described by the nonlinear terms. The main goal of this paper is to investigate the behaviour of a dilute gas flow inducing only the longitudinal shear stress using the classical theory of NSF, with slip and jump boundary conditions, and the regularized 13-moment equations of Grad approaches. In this study, the rarefaction effects are evaluated in the slip and early transition regimes range. We assume that the gas flow is induced with no synergetic contributions from external force fields.

2. Statement of Problem

A monatomic and Maxwell-molecules gas, where the collisions rate is independent of the collision-patterns velocities, is confined within a square microcavity. In this gas, the particles are interacting via a potential with being the interparticles distance [16]. The orthogonal cross section and the origin of the coordinate system are shown in Figure 1. The upper and lower sides are in the opposite motion with a constant velocity . The microcavity walls are kept at uniform and environmental temperature . The macroscopic proprieties of the gas are evaluated for different values of the Knudsen number in the slip and early transition regime, .

Figure 1: Microcavity geometry.

In the kinetic point of view, the macroscopic gas proprieties can be estimated by solving the Boltzmann equation given by

In the steady state, the main unknown of this equation is the velocity distribution function in which is the position vector and is the molecular velocity vector. The distribution function is defined such that the quantity is the number of particles in the phase-volume element near the point . The analytical solution of this kinetic equation is difficult for the multidimensional systems. This is because of the collisions integral which must be computed at each physical-space and velocity-space points. In many of the previous studies, the so-called Bhatnagar-Gross-Krook kinetic (BGK) model is the most used approximation to simplify the collisional operator . In this case, the collisions are modelled by the relaxation of towards the Maxwellian equilibrium distribution that is a function of the particle velocity vector and fluid macroscopic quantities, that is, mass density , hydrodynamic velocity , and temperature . All macroscopic characteristics of the gas flow such as mass density, hydrodynamic velocity, hydrostatic pressure, pressure tensor, and heat flux vector can be related in the steady state to the distribution function , respectively, as follows:

where , and are, respectively, the molecular mass, number density, and peculiar velocity.

3. Continuum Description

The evolution equations for density, velocity, and internal energy given by the conservation laws of mass, momentum, and energy can be written in the general form, if one neglects the external forces, as follows:

where is the hydrostatic pressure in which corresponds to the temperature in energy units. To close the above system, one needs to define the constitutive equations including the stress tensor and heat flux vector with is the Kronecker delta.

The Grad’s moment method is the best-known alternative of the CE expansion [20]. Multiplication Boltzmann equation by functions and subsequent integration over velocity space yields the 13-moment equations of Grad. To solve the closure problem of this set of equations, Grad constructed a phase density by an expansion of the Maxwellian into the Hermite polynomials. However, the discontinuous subshock structures problem is the most failure of this method. Recently, Struchtrup and his coauthors show [21–23] that the Grad’s approach accuracy can be further improved up in the CE expansion by the regularization of the thirteen Grad’s moments. Note that (3a), (3b), and (3c) can be found easily by setting and , respectively, and using the invariance property of the collisions operator and right side of (1). The two other values oflead to the nonlinear constitutive equations of stress tensor and heat flux vector given, respectively, in the so-called R13 set of equations by [19]

The indices in the angular brackets denote the symmetric trace-free parts of tensors. The condensed notation, with indices, used here can be found with more details in the appendix of [19]. The above equations contain the additional, nonlinear, quantities , , and given by considering the appropriate terms of higher order moment equations as (see [24])

If one neglects these quantities, the R13 set of equations is reduced to the well-known set of 13-moment equations of Grad. At the small values of (), the stress tensor and heat flux vector can be written in the CE expansion of Knudsen number powers as

Note that the zeroth order of this expansion corresponds to Euler equations. An analysis of the order of Knudsen number in (4)shows that only the underlined terms are of first order in Knudsen number, while all others are of higher order [25]. Thus, the underlined terms are the NSF contributions which lead to [26]

To solve the R13-moment equations one needs a set of boundary conditions for the normal components of the fluxes. Torrilhon and Struchtrup [27] show that , , , , , and are the only normal components that can be prescribed. The subscripts and denote the tangential and normal, pointing towards the gas, components of the tensors. The Maxwell accommodation model [28] for the boundary conditions in kinetic theory was used to derive boundary conditions for the R13 system, which link the moments of the gas in front of the wall to the tangential wall velocity and wall temperature . The details of the construction procedure for the boundary conditions (BCs) for the R13 equations can be found in [27]; the result reads

where and are the velocity-slip and temperature-jump at the vicinity of the wall, respectively. The parameter is the Maxwell accommodation coefficient where the value is for the full diffuse reflection and for the specular one. The boundary conditions for NSF solution, first order of velocity-slip and temperature-jump, can be derived from the R13 ones by ignoring the higher order terms, , , and :

In this case reduces to .

4. Numerical Scheme

The present scheme is based mainly on the algorithm developed recently by Rana et al. [24]. To solve the differential equations (3a)–(4) in the steady state, it is more convenient to rewrite the system in the matrix and dimensionless form as in [24] to facilitate their numerical resolution:

where , , is the vector of field variables and and are the coefficient matrices in and directions, respectively. The matrix is the production matrix. These matrices are given in Appendix  A of [24]. The matrix equation (10) can be resolved using the finite difference approach in which the BCs (8a)–(8f), in the matrix form, are implemented. Thus, the cavity domain is discretised on a network,, of nodes points. The spatial coordinates of each node, including the boundary ones, are given by and with and . Using the indices notation, (9a), (9b), and (9c) can be written for each node as

where and the matrix are functions of .

The boundary conditions (8a)–(8f) are incorporated in matrix form, and linear extrapolation is used for boundary nodes at walls and ones, that is, wall normal pointing in - or -direction:

where the superscript signs denote the positive/negative directions of the wall normal. The matrices and correspond to the coefficient matrix of field variable in (8a)–(8f)but the vectors and are related to surface inhomogeneity such as wall temperature and velocity . The expressions of these matrices are given in the appendix of [24]. By replacing the derivatives in (10) using a central difference scheme we obtain

Thus, by combining this equation with the boundary conditions given by (12a), (12b), (12c), and (12d) and setting the indices and one gets the equations for different cavity walls. The solution convergence criterion is chosen to be which corresponds in our simulation to .

For more convenience, the results are shown with the dimensionless variables using the following normalisations: ; ; ; ; ; ; , ; , where is the reference mass density related the Knudsen number and reference viscosity by the relation and is the reference pressure.

5. Results and Discussion

In the present paper, we are interested in the study of the viscous and thermal behaviours of a dilute gas flow within the microcavity. The algorithm of this new approach is tested for this standard geometry for two crucial reasons. In the first hand, the cavity is the classical geometry encountered at the lot of devices in engineering area. The second reason is to compare with the previous study especially those of the one-dimension problem such as the classical Couette flow between two parallel plates. The rarefaction effects are discussed in the slip and early transition regimes. The streamlines of gas flow obtained using both NSF and R13 solutions for and are sketched in Figures 2(a)–2(d). We see that the plots are almost similar and both solutions predict a main vortex provided by the walls motion. The secondary vortex observed at the low rarefied case in the cavity centre disappears at the rarefied one . To clarify the velocity variation more, especially, the longitudinal component, the profiles of are plotted as a function of coordinate for different values of , Figures 3(a) and 3(b). In the classical NSF solution, without velocity-slip and temperature-jump BCs, is a linear function of -coordinate between . Using the first order of slip and jump (BCs) we prove the validity NSF solution for velocity profiles up , slip regime, but above this value the solution is not sensitive to the rarefaction effects contrary to the new approach of R13. On the other hand, both approaches predict that the length of so-called Knudsen layer which corresponds to the nonequilibrium area near the walls is about .

(a)

(b)

(c)

(d)

Figure 2: (a, b) Streamlines obtained by NSF (a) and R13 (b) solutions for . (c, d) Streamlines obtained by NSF (c) and R13 (d) solutions for .

(a)

(b)

Figure 3: (a, b) profiles obtained by NSF and R13 solutions for (square), (diamond), (triangle), and (circle).

The thermal behaviour of a gas is studied through the temperature and heat flux changes. Figures 4(a) and 4(b) show the temperature profiles as a function of coordinate. As in the classical study of Couette flow problem [29, 30], the parabolic nature, , with positive coefficient , of temperature profile is shown by both solutions for . Note that the bearing observed in the temperature profiles curvature is due to the low speed or signal problem which is also the major difficulty encountered in the kinetic method like DSMC. With the increase of rarefaction effects, R13 solution predicts a strong inversion in the sharp of profiles especially for ; however, the NSF one is not well sensitive to the Knudsen value. This inversion can be explained by the prevalence of particle-wall collisions against to the particle-particle ones which take place at low rarefied cases.

(a)

(b)

Figure 4: (a, b) Temperature profiles obtained by NSF (a) and R13 (b) solutions for (square), (diamond), and (triangle).

The heat flux transfer without the temperature gradient is another important result provided by the present paper. Figures 5(a) and 5(b) give -component of heat flux vector as a function of -coordinate. At the continuum limit, , both NSF and R13 solutions give almost a constant value; . But, with the increasing of rarefaction degree the R13 approach, unlike the NSF one, leads to horizontal heat flux profiles compatible with the opposite motion of the top/bottom walls. This 1D plots of can be interpreted as an extension of the parallel heat flux component expression given, for gas flow between two infinite plates in opposite motion, by [29]

where is a constant relative to the particles-interaction model. Note that this expression is obtained from the Burnett equations, second order approximation of CE expansion [31, 32]. At the continuum limit, velocity component and temperature profiles show that and . This agrees with the linear variation of predicted by (14). Regarding values given by both solutions (Figures 5(a) and 5(b)), we conclude also that NSF approach predicts almost constant parallel heat flux opposite to R13 one which gives a strong variation. For more details, a two-dimensional plot (Figures 6(a)–6(d)) is used to describe both heat flux transfer and viscous stress behaviours in the whole cavity domain. In Figures 6(a), 6(b), 6(c), and 6(d), the heatlines are sketched on the shear stress contours. These figures confirm the breakdown of NSF solution and the accuracy of R13 one for gas flow with no neglected rarefaction effects.

(a)

(b)

Figure 5: (a, b) profiles obtained by NSF (a) and R13 (b) solutions for (square), (diamond), and (triangle).

(a)

(b)

(c)

(d)

Figure 6: (a, b) Heatlines superimposed on the shear stress contours for . (c, d) Heatlines superimposed on the shear stress contours for .

6. Conclusion

In this paper, we have analysed, numerically, the steady planar Couette flow of a dilute and monatomic gas beyond the scope of the classical Navier-Stokes description using velocity-slip and temperature-jump boundary conditions as a function of rarefaction degree. The behaviour of the gas in rectangular cavity is described by means of extended macroscopic approach R13. The results confirm that with slip and jump boundary conditions, the NS description can be improved. The difference between the results is due to the rarefaction effect which becomes important and the local equilibrium is broken down. The NSF equations cannot explain the variation of heat flux transfer and -component of shear stress tensor. Therefore, the R13 approach is more convenient to capture rarefaction effects heat transfer problem of gas flow in the slip and early transition regimes. Thus, this method can be considered as a good alternative to the kinetic approaches, like DSMC, which are expensive in the time cost.

Competing Interests

The authors declare that they have no competing interests.

References

I. A. Graur and A. P. Polikarpov, “Comparison of different kinetic models for the heat transfer problem,” Heat and Mass Transfer, vol. 46, no. 2, pp. 237–244, 2009.

H. A. Stone, A. D. Stroock, and A. Ajdari, “Engineering flows in small devices: microfluidics toward a lab-on-a-chip,” Annual Review of Fluid Mechanics, vol. 36, pp. 381–411, 2004.

H.-A. Yang, M. Wu, and W. Fang, “Localized induction heating solder bonding for wafer level MEMS packaging,” Journal of Micromechanics and Microengineering, vol. 15, no. 2, pp. 394–399, 2005.

M. Gad-el-Hak, The MEMS Handbook, vol. 1-2, CRC Press, Boca Raton, Fla, USA, 2nd edition, 2005.

M. Hssikou, J. Baliti, Y. Bouzineb, and M. Alaoui, “DSMC method for a two-dimensional flow with a gravity field in a square cavity,” Monte Carlo Methods and Applications, vol. 21, no. 1, pp. 59–67, 2015.

G. E. Karniadakis and A. Beskok, Micro Flows: Fundamentals and Simulation, Springer, New York, NY, USA, 2002.

S. K. Loyalka, “Temperature jump and thermal creep slip: rigid sphere gas,” Physics of Fluids A, vol. 1, no. 2, pp. 403–408, 1989.

S. K. Loyalka and J. W. Cipolla Jr., “Thermal creep slip with arbitrary accommodation at the surface,” Physics of Fluids, vol. 14, no. 8, pp. 1656–1661, 1971.

Y. Sone, Kinetic Theory and Fuid Dynamics, Birkhauser, 2002.

H. Struchtrup and M. Torrilhon, “Higher-order effects in rarefied channel flows,” Physical Review E, vol. 78, no. 4, Article ID 046301, 2008.

S. K. Loyalka, “Velocity profile in the Knudsen layer for the Kramer's problem,” Physics of Fluids, vol. 18, no. 12, pp. 1666–1669, 1975.

T. Ohwada, Y. Sone, and K. Aoki, “Numerical analysis of the Poiseuille and thermal transpiration flows between two parallel plates on the basis of the Boltzmann equation for hard-sphere molecules,” Physics of Fluids A, vol. 1, article 2042, 1989.

Y. Sone, Theoretical and Numerical Analysis of the Boltzmann Equation-Theory and Analysis of Rarefied Gas Flows-Part I, Lecture Notes Department of Aeronautics and Astronautics, Kyoto University, 1998.

Y. Sone, Y. Waniguchi, and K. Aoki, “One-way flow of a rarefied gas induced in a channel with a periodic temperature distribution,” Physics of Fluids, vol. 8, no. 8, pp. 2227–2235, 1996.

Y. Sone, Kinetic Theory and Fluid Dynamics, Birkhäauser, Boston, Mass, USA, 2002.

G. A. Bird, Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Clarendon, Oxford, UK, 1994.

N. G. Hadjiconstantinou, A. L. Garcia, M. Z. Bazant, and G. He, “Statistical error in particle simulations of hydrodynamic phenomena,” Journal of Computational Physics, vol. 187, no. 1, pp. 274–297, 2003.

R. Anirudh, Numerical simulation of rarefied gas flow in micro and vacuum devices [Ph.D. dissertation], 2014.

H. Struchtrup, Macroscopic Transport Equations for Rarefied Gas Flows: Approximation Methods in Kinetic Theory, Springer, Berlin, Germany, 2005.

H. Grad, “Principles of the kinetic theory of gases,” in Handbuch der Physik XII: Thermodynamik der Gase, S. Flugge, Ed., pp. 205–294, Springer, Berlin, Germany, 1958.

H. Struchtrup, “Stable transport equations for rarefied gases at high orders in the Knudsen number,” Physics of Fluids, vol. 16, no. 11, pp. 3921–3934, 2004.

H. Struchtrup, “Scaling and expansion of moment equations in kinetic theory,” Journal of Statistical Physics, vol. 125, no. 3, pp. 569–591, 2006.

H. Struchtrup and M. Torrilhon, “Regularization of Grad's 13 moment equations: derivation and linear analysis,” Physics of Fluids, vol. 15, no. 9, pp. 2668–2680, 2003.

A. Rana, M. Torrilhon, and H. Struchtrup, “A robust numerical method for the R13 equations of rarefied gas dynamics: application to lid driven cavity,” Journal of Computational Physics, vol. 236, no. 1, pp. 169–186, 2013.

H. Struchtrup and M. Torrilhon, “Higher-order effects in rarefied channel flows,” Physical Review E, vol. 78, no. 4, Article ID 046301, 11 pages, 2008.

A. Rana, M. Torrilhon, and H. Struchtrup, “Heat transfer in micro devices packaged in partial vacuum,” Journal of Physics: Conference Series, vol. 362, no. 1, Article ID 012034, 2012.

M. Torrilhon and H. Struchtrup, “Boundary conditions for regularized 13-moment-equations for micro-channel-flows,” Journal of Computational Physics, vol. 227, no. 3, pp. 1982–2011, 2008.

X. J. Gu and D. R. Emerson, “A computational strategy for the regularized 13 moment equations with enhanced wall-boundary conditions,” Journal of Computational Physics, vol. 225, no. 1, pp. 263–283, 2007.

W. Marques Jr., G. M. Kremer, and F. M. Sharipov, “Couette flow with slip and jump boundary conditions,” Continuum Mechanics and Thermodynamics, vol. 12, no. 6, pp. 379–386, 2000.

F. Sharipov, “Data on the velocity slip and temperature jump on a gas-solid interface,” Journal of Physical and Chemical Reference Data, vol. 40, no. 2, Article ID 023101, pp. 1–28, 2011.

C. S. Wang-Chang and G. E. Uhlenbeck, “On the transport phenomena in monatomic gases,” in Studies in Statistical Mechanics, J. de Boer and G. E. Uhlenbeck, Eds., vol. 5, pp. 17–26, North-Holland, Amsterdam, The Netherlands, 1970.

S. Reinecke and G. M. Kremer, “Burnett's equations from a -field theory,” Continuum Mechanics and Thermodynamics, vol. 8, no. 2, pp. 121–130, 1996.