# Continuum Analysis of Rarefaction Effects on a Thermally Induced Gas Flow

Authors: Mohamed Hssikou; Jamal Baliti; Mohammed Alaoui
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

A Maxwell gas confined within a micro cavity with nonisothermal walls is investigated in the slip and early transition regimes using the classical and extended continuum theories. The vertical sides of the cavity are kept at the uniform and environmental temperatureT0, while the upper and bottom ones are linearly heated in opposite directions from the cold value T0 to the hot one TH. The gas flow is, therefore, induced only by the temperature gradient created along the longitudinal walls. The problem is treated from a macroscopic point of view by solving numerically the so-called regularized 13-moment equations (R13) recently developed as an extension of Grad 13-moment theory to the third order of the Knudsen number powers in the Chapman-Enskog expansion. The gas macroscopic properties obtained by this method are compared with the classical continuum theory results (NSF) using the first and second order of velocity slip and temperature jump boundary conditions. The gas flow behavior is studied as a function of the Knudsen number (Kn), nonlinear effects, for different heating rates T0/TH. The micro cavity aspect ratio effect is also evaluated on the flow fields in this study.


---

## Body

## 1. Introduction

The frequent use of small-size devices in several applications needs types of equipment that dissipate maximum amount of heat per unit area. With the development of micro-electro-mechanical systems (MEMS), the heat transfer mechanism has gained recently a great interest. The micro pumps, micro ducts, micro nozzles, micro turbines, and micro valves are some typical examples of these small devices involving liquid and gas flows [1–3]. Regarding the dimensions of these devices, the continuum assumptions of the gas flow break down and some deviations are observed in respect of the macro scale case. The use of velocity slip and temperature jump boundary conditions proves that the rarefaction degree has a significant effect on the gas micro flows or low-pressure gas problems (Karniadakis and Beskok 2005, [4]). To estimate the degree of rarefaction, the so-called Knudsen number defined as Kn~λ/L is used to classify the flow regimes, where λ and L represent, respectively, the mean free path and a characteristic length of the system. For Kn<10-3 the flow can be considered as a continuum medium; hence the classical Navier-Stokes and Fourier equations are sufficient to describe the gas flow behavior. In the range of 10-3≤Kn<0.1 (slip regime), the gas flow fields can be described using these equations if the typical boundary conditions of velocity slip and temperature jump are applied at the walls [5]. After the slip limit, 0.1≤Kn<10, i.e., transition regime, the nonequilibrium effects become more important even in the bulk of flow and, therefore, the continuum theory breaks down [6, 7]. In this case, the gas flow must be described from a kinetic point of view by following the governing equation of Boltzmann (3). However, the high-dimensionality of this equation makes its direct solving extremely expensive in the discrete phase space. So far, the direct simulation Monte Carlo (DSMC) method is the main practical method used to characterize the strong nonequilibrium flows in NEMS/MEMS [8, 9]. However, the DSMC method is very expensive both in computational time and memory requirements, especially for low-speed flow in microelectromechanical devices (MEMS) and nanodevices [10]. In thermal cavities, when Kn number increases, the intermolecular- and molecular-surface interactions decrease [11]. Therefore, a huge number of particles are needed to capture the flow features. In the macroscopic point of view, the problems involving Couette and Fourier flows have been treated enough with isothermal walls in the previous study [12–15]. However, minor attempts have been observed for treating the problems with nonisothermal walls boundary conditions [16–18]. In this case, a thermal slip gas flow is induced under the effects oftranspirationorthermal creep flows; i.e., the gas is forced into motion at the boundaries. The thermal creep flow of a rarefied gas within an enclosure is treated deterministically by solving, on the one hand, the Boltzmann-Shakhov model equation and using DSMC as a stochastic approach in the range of 10-2<Kn<10 [17]. On the other hand, in such problem, a competition between thermal creep flow and thermal stress flow induced by a gradient of heat flux in the bulk is investigated as a function of rarefaction degree (Kn) [19]. At the micro scale and rarefied conditions, the temperature inhomogeneity of a gas system can lead to a variety of flow phenomena such as theghost effects observed at the continuum limit [20]. The fact that many of these micro devices involve slip and early transition regime gas flows is the first motivation of this continuum study. The second one is to evaluate thethermal slip contribution in the velocity slip and temperature jump phenomena near the heated walls, which is proportional to the tangential temperature gradient [21]. This effect is the main correction made by the second order of velocity slip and temperature jump in the NSF approach. Investigating thermal behaviors under rarefied condition due to high power densities requires the application of efficient techniques to predict allowable performance limits of these systems [22]. By comparing the DSMC method and R13 solution, it is shown that the R13 equations can be solved in short computation time as compared to those for DSMC method [23]. Thus, the main goal of this paper is to investigate the rarefaction effects on the macroscopic proprieties of a dilute gas flow induced only by nonisothermal walls effects, thermal creep flow, using the classical and extended continuum-based model. The regularized 13-moment equations (R13), developed by [24], are solved numerically to capture the strong rarefaction or nonlinear effects that occur in the early transition regime. In this study, we discuss also the validity of the first (NSF1) and second (NSF2) orders of velocity slip and temperature jump boundary conditions in the Navier-Stokes and Fourier solution, within the rarefaction range of 0.05≤Kn≤0.3. In this study, we assume that the gas is not subjected to any external body force. In the next section, we present a brief description of the problem, after we recall the set of basic equations and the boundary conditions used in both solutions of NSF and R13.


## 2. Statement of Problem

A Maxwell and monatomic gas is confined within a two-dimensional heated microcavity with orthogonal cross section ofH×L being shown in Figure 1. The left and right walls of the cavity are kept at a uniform and environmental temperature T0=273K while the bottom and upper sides are linearly heated from the cold value T0 to hot one TH in opposite directions; i.e., Txy=0,H=T0,H±TH-T0x/L. The macroscopic properties of the gas are evaluated for different values of Knudsen number in the slip and early transition regimes for the hot temperatures 1.25T0,1.5T0,1.75T0, and 2T0. Two values of the cavity aspect ratio, AR=H/L, are considered in this study (AR=2 and 1). Figure 1 shows the configuration of the cavity in a square geometry (AR=1). The flow fields are evaluated as a function of Knudsen number Kn defined as(1)Kn=μ0v0p01Lwhere R, μ0, and p0 denote, respectively, the gas specific constant, the reference viscosity at the temperature T0, and the hydrostatic pressure given by p0=ρ0RT0. Also v0=2RT0 is the most probable molecular velocity which is taken as the characteristic velocity.Figure 1
Micro cavity geometry.


## 3. Continuum Model Based Description

The evolution equations for density, velocity, and internal energyε are given by the conservation laws of mass, momentum, and energy that can be written in the general form as follows:(2a)∂ρ∂t+∂ρvi∂xi=0,(2b)∂ρvi∂t+∂Pik+ρvivk∂xk=0,(2c)∂ρε+1/2v2∂t+∂ρε+1/2v2vk+Pikvi+qk∂xk=0,where ρ,vk,Pik, and qk denote, respectively, mass density, macroscopic velocity, pressure tensor, and heat flux vector components. Note that the pressure tensor is related to the hydrostatic pressure and stress tensor by the relation: Pij=pδij+σij where δij represents the Kronecker delta. In the frame of ideal gas approximation, the hydrostatic pressure is given by p=ρRT=ρθ where θ corresponds to the temperature in energy units. To close this set of equations, one needs to define the constitutive equations including the relationship between the stress tensor σ, the heat flux vector q, and the macroscopic variables derivatives. As an extension of the classical NSF theory, Grad introduces the balance equations corresponding to higher order moments of the particle velocity distribution function fr,c,t [26] in which the evolution is governed by the Boltzmann equation given by (3)∂∂t+c.∇rf=Jf,f′,where J[f,f′] and c are, respectively, the collisions operator and the molecular velocity.Due to the collision’s random nature, the exact analytic solution of (3) remains difficult. For this reason, some approximations are used to simplify this term for the small deviations in respect of the equilibrium state. In this context, it has been shown recently that the Grad’s approach accuracy can be further improved up to O(Kn3) of the Knudsen number powers in the Chapman-Enskog expansion by regularization of the moment equations of Grad [27, 28]. Note that the NSF solution with velocity slip and temperature jump boundary conditions leads to the first order in this expansion.In the13-moment theory of Grad, the gas state is described by the set of 13 moments of functions φA=m{1,ci,1/2C2,CiCj,1/2C2Ci}. The Multiplication of (3), at the steady state, by φA and subsequently by integration over velocity space yields the governing equations of moments: mass density, hydrodynamic (bulk) velocity, hydrostatic pressure, stress tensor, temperature, and heat flux vector related to the distribution function by(4)ρx,y=m∫fr,cdc,(5)vx,y=mρx,y∫cfr,cdc,(6)px,y=m3∫C2fr,cdc,(7)Pijx,y=m∫CiCjfr,cdc,(8)Tx,y=m3nKB∫C2fr,cdc,(9)qx,y=m2∫C2Cfr,cdcwhere C=c-v is the peculiar velocity.By settingφA=m,mci and 1/2mC2, respectively, and using the invariance property of the collisions operator, the right side of (3), one can find the balance equations ((2a), (2b), (2c)). The two other values of φA lead to the nonlinear constitutive equations of stress tensor σij and heat flux vector qi, respectively, that are given in the so-called R13 approach by [27] (10)∂σij∂t+∂σijvk∂xk+45∂qi∂xj+2p∂vi∂xj+2σki∂vj∂xk+∂mijk∂xk=-pμσij,(11)∂qi∂t+∂qivk∂xk+52p∂θ∂xi+52σik∂θ∂xk+θ∂σik∂xk-σikθρ∂ρ∂xk-σikρ∂σkl∂xl+75qk∂vi∂xk+25qk∂vk∂xi+25qi∂vk∂xk+12∂Rik∂xk+16∂∆∂xi+mikl∂vk∂xl=-23pμqi,where the indices in the angular brackets denote the symmetric trace-free parts of tensors. The R13 equations closure contains, therefore, additional quantitiesmijk, Rik, and ∆, corresponding to the higher order moments, given by [24, 28, 29](12a)∆=-σklσklρ+6σklσklNSFρ+566qkqkNSFp-12μpθ∂qk∂xk-θqk∂ln⁡p∂xk,(12b)Rij=-47σkiσjkρ+247σkiσjkNSFρ+19275qiqjNSFp-245μpθ∂qi∂xj-θqi∂ln⁡p∂xj,(12c)mijk=815σijqkNSFp+45qiσjkNSFp-2μpθ∂σij∂xk-θσij∂ln⁡p∂xk.If one neglects these quantities, the R13 set of equations will be reduced to the well-known set of Grad’s 13-moment equations. At the small values ofKn (Kn<<1), the stress tensors and heat flux components can be written in the Chapman-Enskog expansion of Knudsen number powers as follows.(13)σij=σij0+Knσij1+Kn2σij2…(14)qi=qi0+Knqi1+Kn2qi2…The zeroth order of this expansion leads to the Euler equation while the first one corresponds to the classical NSF theory. In this case, the Newton and Fourier laws, respectively, give the stress and heat flux vector for Maxwell gas:(15a)σijNFS=-2μ∂vi∂xj,(15b)qiNFS=-154μ∂θ∂xi.As in any problem involving partial differential equations, the solving of R13-moment equations requires a set of boundary conditions that one must specify. In this context, Torrilhon and Struchtrup have shown kinetically thatvn,στn,qn,Rτn,mnnn, and mττn are the only components that can be prescribed [30]. The subscripts τ and n denote the tangential and normal components of the tensors, respectively, where the wall normal is pointing toward the gas. Using the Maxwell accommodation model for boundary conditions results in a kinetic link between the moments in front of the wall and the tangential wall velocity vτw and the wall temperature θw; the R13-BCs are written as(16a)vn=0(16b)στn=-χ/2-χ2πθPVτ+15qτ+12mτnn,(16c)qn=-χ/2-χ2πθ2PT-12PVτ2+12θσnn+115∆+528Rnn,(16d)Rτn=χ/2-χ2πθ6PTVτ+PθVτ-PVτ3-115θqτ-12θmτnn,(16e)mnnn=χ/2-χ2πθ25PT-35PVτ2+75θσnn+175∆-114Rnn,(16f)mττn=-χ/2-χ2πθ15PT-45PVτ2+114Rττ+θσττ-1/5θσnn+1155∆.Here P=ρθ+1/2σττ-1/120∆/θ-1/28Rττ/θ.Note thatVτ=vτ-vτw and T=θ-θw represent, respectively, the velocity slip and temperature jump at the vicinity of the walls. The parameter χ is the Maxwell accommodation coefficient whose value is χ=1 for the full diffuse reflection and χ=0 for the specular one. The effect of such coefficient on thermal transpiration phenomena is evaluated using velocity dependent Maxwell (VDM) boundary condition and by means of DSMC method [31]. A heated microcavity with specular walls is also treated by means of NSF and moments theories [32, 33].


## 4. Numerical Scheme

To solve the above differential equations ((2a), (2b), (2c), (10), (11), (12a), (12b), and (12c)) in the steady state, it is more convenient to rewrite the equations in the condensed and matrix form:(17)Aω∂ω∂x+Bω∂ω∂y+1KnPωω=0.Here ω=[ρ,vx,vy,θ,qx,qy,σxx,σxy,σyy,Rxx,Rxy,Ryy, mxxx,mxxy,mxyy,myyy,∆] is the vector of field variables; Aω, Bω, and Pω are, respectively, the coefficient matrices in x,y- directions and the production matrix. The matrix equation (17) is solved by means of finite difference approach. The cavity domain is discretized on a network of nodes in both directions. The R13 solution is obtained using the boundary conditions prescribed by a set of equations ((16a)–(16f)). In the slip flow regime, the Navier-Stokes and Fourier equations solution must be obtained with slip and jump boundary conditions. By ignoring the higher order term smijk, Rik, ∆ and replacing the stress tensor and heat flux vector with their corresponding expressions of NSF, as indicated in the equations ((15a), (15b)), the first order of NSF-BCs can be written, therefore, as follows: (18a)vn=0,(18b)στnNSF1=-χ/2-χ2πθPVτ-34μ∂θ∂xτ,(18c)qnNSF1=-χ/2-χ2πθ2PT-12PVτ2-μθ∂vn∂xn,where P=ρθ+1/2σττ.Meanwhile, the second order of these BCs can be obtained using the second order of stress tensor and heat flux vector described by the Burnett equations [24, 34]:(19a)vn=0,(19b)στnNSF2=-χ/2-χ2πθPVτ+15qτNSF+mτnn22+15μ2ρ4516∂2vτ∂xk∂xk-134∂2vk∂xk∂xτ-118Rτn2θ,(19c)qnNSF2=-χ/2-χ2πθ2PT-12PVτ2+12θσnnNSF+∆215+13Rnn263-μ2ρ4516∂2vn∂xk∂xk-134∂2vk∂xk∂xn,where the terms with the superscript (2) refer to the second order corrections which are obtained from the R13 constitutive relations (12a), (12b), and (12c) by replacing stress tensor and heat flux vector with their NSF expressions; these terms are given as follows:(20a)Rij2=-245μρ∂qiNSF∂xj,(20b)mijk2=-2μρ∂σijNSF∂xk,(20c)∆2=-12μρ∂qkNSF∂xk.For more convenience, the results are shown with the dimensionless variables using the following normalizations:(21)x,y=x,yL,Vi=viv0,qi=qiρ0θ0v0,ρ-=ρρ0,θ-=θθ0≡TT0,σ-ij=σijρ0θ0,R-ij=Rijρ0θ02,∆-=∆ρ0θ02andm-ijk=mijkρ0θ0v0.


## 5. Results and Discussion

The behavior of gas flows in the slip and early transition regimes, often encountered in MEMS devices, is one of the challenging problems which need further investigations. In this regime, the mean free path of gas is at the same order of or even much larger than the system characteristic length, the continuum assumption will break down, and the traditional CFD techniques will, therefore, lead to large errors [35]. The rarefaction effects influence the flow, such as heat flux induced by gradients of stresses even if there is no temperature gradient [36]. For validation purpose, firstly, we compare the DSMC and R13 results by considering the test case treated by Vargas et al. [18]. In this case, the gas flow is induced, on the one hand, by the thermal creep force created along the nonisothermal walls and, on the other hand, by a temperature gradient between the other walls. To stay in the considerations of macroscopic approaches, a small temperature gradient and moderate Knudsen number values are considered in this test case, T0/TH=0.9, Kn=0.01 and 0.1. In this case, the gas flow is subjected to two thermal forces: atranspiration force due to the linear profile imposed on the walls-temperature andthermal stress induced by the temperature gradient in bulk of flow. These forces may be acting in the same or opposite direction of the gas flow. Figure 2 shows the streamlines overlaid on the temperature contours obtained by both approaches.Figure 2
Streamlines and temperature contours comparison forT0/TH=0.9 and Kn=0.1 between R13 (right) and DSMC method [18] (left) adapted to the Vargas’s problem (vertical, top-to-down, linearly heated cavity).R13 and DSMC predict two similar primary vortices in the flow field which push the flow from cold to hot region at the side walls. Kinetically, these patterns are interpreted by the fast motion of hot particles compared to particles coming from the colder region, which are slower. As a reaction, the wall is pushed toward the colder region, or, when the wall is at rest, the rarefied gas is driven from cold to hot, which is observed here. Both approaches give almost similar temperature distribution in the cavity.In contrast with the classical problems of viscous flows induced by the relative motion of walls, this study focuses on the pure thermal creep flow induced only, in opposite directions, by the effects of nonisothermal walls as shown in Figure1. To stay in the frame of small deviations from equilibrium state, we restrict our attention to the macroscopic study of gas flow fields for a moderate temperature TH=1.5T0 [37]. Then, we evaluate the effect of TH on the maximum of tangential thermal-velocity along the bottom walls. By analogy with the viscous Couette flow [38, 39], the present problem of “thermal Couette flow” is treated by means of classical approach of NSF, with first and second orders of slip and jump BCs, and extended continuum theory.The breakdown of the classical theory between the slip and transition regimes limit,Kn=0.1, is clarified by solving the regularized 13-moment equations (R13). At low rarefaction degree, Kn=0.05, the R13 approach predicts a primary vortex centered at x=y=0.5L while the NSF1 and NSF(2) solutions envisage, inside this main vortex, two symmetric secondary vortices corresponding to the flow stagnation. Near longitudinal sides, the tangential temperature gradient forces a creep-driven flow in a circulatory motion from the cold to hot region. In addition to these eddies induced by the thermal convection phenomenon, thermal creep flow, two secondary symmetric counter-circulating vortices at the hot cavity corners are observed; see Figures 3(a), 3(c), and 3(e). The inverted transpiration force induced at the hot corners is the origin of these secondary vortices.Figure 3
Streamlines overlaid on the temperature contours forKn=0.05 and 0.3 using NSF approaches (a-d) and R13 approach (e, f).
(a)
NSF1 for K n = 0.05 (b)
NSF1 for K n = 0.3 (c)
NSF2 for K n = 0.05 (d)
NSF2 for K n = 0.3 (e)
R13 for K n = 0.05 (f)
R13 for K n = 0.3ForKn=0.3, unlike first order solution which is not much sensitive to the rarefaction increase (Figure 3(b)), the second order solution (Figure 3(d)) predicts a little change in the center and size of vortices which correspond to the rarefaction effects contribution. In contrast to both NSF approaches, the R13 solution allows capturing the interplay between thermal stress and transpiration force, which are responsible for the streamlines inversion and disappearance observed at hot corners in the rarefied case. The third term in the R13-constitutive equation (10) expresses the thermal stress contribution. This is due to the predominance of nonequilibrium effects, with respect to thermal convection phenomenon, observed in the slightly rarefied case (Kn=0.05). Indeed, the central vortex vanishes and the hot-corners vortices agglomerate to give an inverted circling of gas flow. This is in good agreement with the simplestwindmill experiment realized by Sone [25] where the qualitative results are published through a brief communication [25] and recorded on a videotape (VHS). In this experiment, a rectangular glass plate of 70×200mm was set in space with its longer sides in the vertical direction and electrically heated from the lower end of the plate (in back). A windmill was placed in front of the plate to detect the vertical flow induced by the heating process. The system is placed within a glass cylindrical chamber in which the pressure change is well controlled; the experimental protocol is shown in Figure 4. At the high-pressure condition, low rarefaction degree, the thermal convection phenomenon provides the windmill motion. This motion becomes increasingly slow with the pressure decrease. At low pressure, ~40Pa, the windmill motion is reversed by rarefaction effects. Experimentally, it fined as the pressure decreases; downward flow, a flow in the direction of the temperature gradient along the plate, grows and finally overcomes the thermal creep flow. This experiment was performed later with the same qualitative results in various systems with different windmills and plates.Figure 4
Windmill experimental apparatus [25].The accuracy of theR13 solution, in the rarefied case, is also provided by a wide range of temperature variation, 1.04≤T/T0≤1.24, unlike the closed one predicted by NSF solutions, especially the second order. The gas flow behavior is, therefore, influenced by both rarefaction and nonisothermal wall effects. To more understand this critical change, observed in the flow streamlines at the rarefied case (Figure 3(f)), the velocity and heat flux fields are evaluated in the vertical center-line; i.e., x=L/2. Figures 5 and 6 show the velocity and heat flux components as a function of y-coordinate. The rarefaction effects strongly influence the velocity fields, Vx and Vy. This higher sensitivity of velocity components is due to the small velocity-magnitude induced only by the longitudinal temperature gradient along the walls. Even at the slightly rarefied case, the classical NSF solutions cannot capture the nonlinear phenomenon contribution, which is more highlighted by solving the regularized 13-moment equations.Figure 5
x -velocity and heat flux components profiles at the vertical center-line of the cavity for Kn=0.05 (a), 0.1 (b), and 0.3 (c) using NSF1, NSF2, and R13 approaches; the black and blue symbols denote x-velocity and heat flux components, respectively.
(a) (b) (c)Figure 6
y -velocity and heat flux components profiles at the vertical center-line for Kn=0.05 and 0.3 using NSF1, NSF2, and R13 approaches.
(a) (b) (c) (d)This velocity-rarefaction dependence is consistent with the experiment results for different pressure values. The combined effects of rarefaction and nonisothermal wall-driven is proved later by Sone through the improved version of his simple windmill experiment. The results show that the windmill rotation speed is strongly linked to the temperature gradient, i.e., the heating rateT0/TH, and the rarefaction degree (the pressure value in the chamber).TheR13 approach also provides a good clarification of the so-called Knudsen-layer which corresponds to the nonequilibrium area induced by the rarefaction effects near the heated walls, located at y=0,H. The rarefaction effects also influence the profiles of heat flux components. At the slightly rarefied case, Kn=0.05, the NSF approaches give rather similar profiles in the bulk region over the range of 0.25≤y/L≤0.75. However, the validity of first and second orders BCs of NSF solution decreases with the increase of Knudsen number value before its breakdown for Kn=0.3. Therefore, the enlargement of the Knudsen-layers in terms of Kn value, induced by the thermal creep flow, is more captured by the R13 solution, unlike the second order NSF2 which predicts only a small correction with respect to the first one NSF1.This nonequilibrium region, which appears in the profiles of longitudinal components as a function of the normal coordinate,y, does not take place in the vertical ones. This is illustrated by the Vy and qy profiles for Kn=0.05 and 0.3 (Figure 6). With the second order of velocity slip and temperature jump BCs, the NSF theory predicts an improved correction of y-velocity and heat flux components in the rarefied case; see Figures 6(c) and 6(d). Figure 6 shows that the three solutions are almost confounded at the slightly rarefied case Kn=0.05. However, for Kn=0.3, the classical NSF theory, especially the first order one, cannot describe correctly the thermally driven rarefied gas flow behavior. Using the first order of slip and jump BCs, the slip velocity phenomenon is only induced in the proximity of walls, i.e., the Knudsen-layer, by the viscous slip flow. The sufficiency of this order of correction has been proved recently by comparing the R13 results and those obtained by solving the linearized Boltzmann equation for viscous slip and transpiration flow problems [40, 41]. However, when the tangential temperature gradient is applied, i.e., ∂T/∂xτ≠0, contribution of the second order term (~Kn2) in the slip BCs must be considered. This explains the accuracy of NSF2 correction in the profile Vx=f(y), with respect to the first one, NSF1.To evaluate the combined effects of temperature gradient and rarefaction degree, we show the maximum value ofVx along the bottom wall as a function of Knudsen number for different temperatures, 1.25≤TH/T0≤2 (Figure 7). In conclusion, the maximum velocity along the wall depends on the gas flow regime defined by the Knudsen number value. In the slip regime, Kn≤0.1, Vx,max decreases with the increase of rarefaction degree. After this limit, transition regime, the gas particles undergo a strong thermal motion at the wall proportional to Kn and TH values.Figure 7
Maximumx-velocity along the bottom wall for different values of TH.To understand more about the heat transfer process induced by the temperature gradient, the heat-lines overlaid on the thermal shear stress contours obtained by NSF2 and R13 approaches are shown in Figure8. Comparing the figures, one can conclude that both approaches provide a rather similar distribution of heat flux isolines at the slightly rarefied case (Kn=0.05), while the longitudinal thermal stress magnitude is strongly influenced by the nonlinear mechanism of momentum transfer at the boundary. The nonlinear phenomena contribution, described with more details from the kinetic point of view, is well considered by the extended continuum theory of R13 by means of the terms Rij, mijk, and ∆, unlike the NSF2 solution which gives a small range of change [41]. However, when the Knudsen number becomes larger (Kn=0.3), the rarefaction effects influence the gas flow. The heat flux caused by the gradients of stresses and thermal creep is better described by R13 method [42, 43]. The rarefaction effect causes a change in the heat-lines curvature (Figure 9(d)), which is not predicted by NSF solution. This is due to the faster motion of particles coming from the hot corners than those coming from the cold ones. The stress quadratic flow, as well as nonlinear thermal stress effects, becomes larger in the flow bulk region [20].Figure 8
Heat-lines overlaid on the stressσxy contours for Kn=0.05 and 0.3 using NSF2 (a, b) and R13 (c, d) approaches.
(a)
NSF2 for K n = 0.05 (b)
NSF2 for K n = 0.3 (c)
R13 for K n = 0.05 (d)
R13 for K n = 0.3Figure 9
Streamlines overlaid on the temperature contours atKn=0.05 and 0.3 for AR=2 using NSF1 and R13 approaches.
(a)
NSF1 for K n = 0.05 (b)
R13 for K n = 0.05 (c)
NSF1 for K n = 0.3 (d)
R13 for K n = 0.3The mutual interaction between the longitudinal walls, subject to the temperature gradient, on the flow pattern is evaluated by considering higher cavity (H=2L), i.e., AR=2. Figure 9 shows the streamlines superimposed on the temperature contours for Kn=0.05 and 0.3 using both NSF1 and R13 approaches. Comparing these flow fields with the corresponding ones in the square case (AR=1) in Figure 5, one can conclude that both solutions are sensitive to the cavity aspect ratio and rarefaction effects. For Kn=0.05, the primary vortices observed between the cold corners in the square case, using NSF approach, are detached to form two secondary vortices on the side of lateral walls. However, at the rarefied case, the situation is reversed and the hot-corners vortices begin to agglomerate to form one primary vortex between the hot corners. However, the extended continuum method, R13, predicts rather similar fields but with small edges at the cold corners. For Kn=0.3, two secondary vortices appear, under the vertical flow elongation effect, within the inverted vortex obtained between the hot corners in the square cavity case. These vortices are made to the flow stagnation appearing at y=0.5L and 1.5L. On the other hand, the temperature contours show that the R13 approach predicts a larger range of temperature variation, T/T0, in agreement with the rarefaction degree: [1.04-1.24] in square case and [1.02-1.22] for AR=2. However, with the NSF solution, the temperature range varies from [1.08-1.18] to [1.065-1.13] in the higher cavity. The sensitivity of thermal behavior with the cavity aspect ratio is, therefore, better captured by the extended continuum theory of R13 than by the classical one (NSF).


## 6. Concluding Remarks

In the present study, the behavior of a rarefied gas flow, induced thermally by nonisothermal wall effects within a microcavity, is investigated in the slip and early transition regimes. The problem is treated in the frame of classical and extended macroscopic theory by solving numerically, on the one hand, the Navier-Stokes and Fourier equations using the first and second orders of velocity slip and temperature jump and, on the other hand, the R13 set of equations. The gas flow sensitivity to the rarefaction degree is clarified using the main macroscopic flow fields. The results prove the break downs of NSF solutions at slip-transition regimes limit. Thus, the extended continuum solution allows the capture of the inverted motion of rarefied gas flow created by the thermal creep flow at the walls. This critical change in the flow streamlines is in a good agreement with thewindmill experiment results realized by Sone. The nonlinear thermal stress and thermal creep co-contribution on the flow fields at the rarefied case have been illustrated by means of the R13 method. In view of different geometries that one can meet in several industrial applications, two different aspect ratios are evaluated in our study.



---

*Source: 5084098-2019-04-30.xml*