Burnett simulation of flow and heat transfer in micro Couette flow using second-order slip conditions

The conventional Burnett equations with second-order velocity slip and temperature jump conditions were applied to the steady-state micro Couette flow of a Maxwellian monatomic gas. An analytical approach as well as a relaxation method was used to determine the velocity slip and temperature jump at the wall. Convergent solutions to the Burnett equations were obtained on arbitrary fine numerical grids for all Knudsen numbers (Kn) up to the limit of the equations’ validity. The Burnett equations with second-order slip conditions indicate a much better agreement with DSMC data over the first-order slip conditions at high Kn. The convergent Burnett solutions were obtained in orders of magnitude quicker than that with the corresponding DSMC simulation. The augmented Burnett equations were also introduced to model the flow but no obvious improvement in the results was found.     

微尺度泊肃叶流的高阶连续模型分析

分别从分子运动论及连续流理论出发,对体积力驱动的微尺度平面泊肃叶(Poiseuille)流的横向分布特征进行了分析. 分子水平模拟采用直接模拟蒙特卡罗(direct simulation Monte Carlo, DSMC)方法;连续流理论则主要考察了伯内特(Burnett)及超伯内特(Super-Burnett)等高阶连续模型,在平行流假设下,获得一组高阶非线性常微分方程,补充完整的边界条件,并应用龙格-库塔(Runge-Kutta)方法求解. 结果表明,即使对于过渡领域流动,高阶连续模型可以给出与DSMC 结果完全相符的压力分布,而速度分布当努森(Knudsen)数约为0.2时即在壁面开始出现偏差;对于温度的横向分布,伯内特模型回复到纳维-斯托克斯(Navier-Stokes)水平,不能得到与DSMC一致的双峰结构,而超伯内特模型在滑移流动领域与DSMC定性相符,在过渡领域却仅能正确预测主流区温度分布,壁面附近差异明显;横向热流与纳维-斯托克斯模型预测接近,但机理上存在本质区别. 本文结果提示选用连续模型时,不仅要根据流动参数来判断,还可以根据所关注的物理量来进行调整,适度扩大连续模型的适用范围. 但即使采用高阶本构关系,连续模型仍然不能完全描述壁面附近区域的非平衡效应(如努森层效应),这是试图扩大连续模型适用范围时必然会遇到的困难.     

Navier-Stokes方程二阶速度滑移边界条件的检验

对微尺度气体流动,Navier-Stokes方程和一阶速度滑移边界条件的结果与实验数据相比,在滑移区相互符合,在过渡领域则显著偏离.为改善Navier-Stokes方程在过渡领域的表现,有些研究者尝试引入二阶速度滑移边界条件,如Cercignani模型,Deissler模型和Beskok-Karniadakis模型.以微槽道气体流动为例,将Navier-Stokes方程在不同的二阶速度滑移模型下的结果与动理论的直接模拟Monte Carlo(DSMC)方法和信息保存(IP)方法以及实验数据进行比较.在所考察的3种具有代表性的二阶速度滑移模型中,Cercignani模型表现最好,其所给出的质量流率在Knudsen数为0.4时仍与DSMC和IP结果相符;然而,细致比较表明,Cercignani模型给出的物面滑移速度及其附近的速度分布在滑流区和过渡领域的分界处(Kn=0.1)已明显偏离DSMC和IP的结果.     

Asymptotic theory of the Boltzmann system,for a steady flow of a slightly rarefied gas with a finite Mach number: General theory

A steady rarefied gas flow with Mach number of the order of unity around a body or bodies is considered. The general behaviour of the gas for small Knudsen numbers is studied by asymptotic analysis of the boundary-value problem of the Boltzmann equation for a general domain. The effect of gas rarefaction (or Knudsen number) is expressed as a power series of the square root of the Knudsen number of the system. A series of fluid-dynamic type equations and their associated boundary conditions that determine the component functions of the expansion of the density, flow velocity, and temperature of the gas is obtained by the analysis. The equations up to the order of the square root of the Knudsen number do not contain non-Navier–Stokes stress and heat flow, which differs from the claim by Darrozes (in Rarefied Gas Dynamics, Academic Press, New York, 1969). The contributions up to this order, except in the Knudsen layer, are included in the system of the Navier–Stokes equations and the slip boundary conditions consisting of tangential velocity slip due to the shear of flow and temperature jump due to the temperature gradient normal to the boundary.     

Mass flow rate measurements in gas micro flows

The main objective of this experimental investigation on the gas flow slip regime is to measure the mass flow rate in isothermal steady flows through cylindrical micro tubes. Two technical procedures devoted to mass flow rate measurements are compared, and the measured values are also compared with the results yielded by different approximated analytical solutions of the gas dynamics continuum equations. Satisfactory results are obtained and the way is clearly opened to measuring mass flow rates for higher Knudsen numbers, over all the micro flow transitional regime.     

Kramers’ problem and the Knudsen minimum: a theoretical analysis using a linearized 26-moment approach

A set of linearized 26 moment equations, along with their wall boundary conditions, are derived and used to study low-speed gas flows dominated by Knudsen layers. Analytical solutions are obtained for Kramers’ defect velocity and the velocity-slip coefficient. These results are compared to the numerical solution of the BGK kinetic equation. From the analysis, a new effective viscosity model for the Navier–Stokes equations is proposed. In addition, an analytical expression for the velocity field in planar pressure-driven Poiseuille flow is derived. The mass flow rate obtained from integrating the velocity profile shows good agreement with the results from the numerical solution of the linearized Boltzmann equation. These results are good for Knudsen numbers up to 3 and for a wide range of accommodation coefficients. The Knudsen minimum phenomenon is also well captured by the present linearized 26-moment equations.     

二阶滑移边界对微型气浮轴承稳态性能的影响

考虑微型气浮轴承的尺寸特征,内部气流不再满足连续流的假设,根据Knudsen数可确定内 部气流为滑移流. 分别利用一阶速度滑移模型和二阶速度滑移模型对连续流的状态方程进行 修改,得到一阶滑移流和二阶滑移流机制下修正的雷诺方程. 利用有限差分法对连续流、一 阶滑移流和二阶滑移流的雷诺方程分别求解,得出相应的承载力和偏位角. 经过对比分析, 发现采用滑移流模型得到的轴承的稳态力学性能与连续流机制的结果存在较大差异,一阶滑 移流与二阶滑移流间的差异随偏心率增加而增加. 说明在MEMS环境下必须考虑滑移流效应 对微型气浮轴承稳态力学性能的影响. 在大偏心率工作状态下,二阶滑移流模型能够得到最 好的结果.     

Failures of the Burnett and super-Burnett equations in steady state processes

Linearized Burnett and super-Burnett equations are considered for steady state Couette flow. It is shown that the linear super-Burnett equations lead to periodic velocity and temperature curves, i.e. unphysical solutions. The problem is discussed as well for the so-called augmented Burnett equations by Zhong et al. (AIAA Journal 31, 1036-1043 (1993)), and for the recently introduced regularized 13 moment equations (R13) of Struchtrup and Torrilhon (Phys. Fluids 15(9), 2668-2680 (2003) ). It is shown that both theories exhibit proper Knudsen boundary layers for velocity and temperature. However, the heat flux parallel to the wall has different signs for the Burnett and the R13 equations, and a comparison with DSMC results shows that only the R13 equations predict the proper sign.Received: 30 March 2004, Accepted: 22 April 2004, Published online: 22 February 2005PACS: 51.10.-y, 47.45.-n     

Knudsen’s Permeability Correction for Tight Porous Media

Various flow regimes including Knudsen, transition, slip and viscous flows (Darcy’s law), as applied to flow of natural gas through porous conventional rocks, tight formations and shale systems, are investigated. Data from the Mesaverde formation in the United States are used to demonstrate that the permeability correction factors range generally between 1 and 10. However, there are instances where the corrections can be between 10 and 100 for gas flow with high Knudsen number in the transition flow regime, and especially in the Knudsen’s flow regime. The results are of practical interest as gas permeability in porous media can be more complex than that of liquid. The gas permeability is influenced by slippage of gas, which is a pressure-dependent parameter, commonly referred to as Klinkenberg’s effect. This phenomenon plays a substantial role in gas flow through porous media, especially in unconventional reservoirs with low permeability, such as tight sands, coal seams, and shale formations. A higher-order permeability correlation for gas flow called Knudsen’s permeability is studied. As opposed to Klinkenberg’s correlation, which is a first-order equation, Knudsen’s correlation is a second-order approximation. Even higher-order equations can be derived based on the concept used in developing this model. A plot of permeability correction factor versus Knudsen number gives a typecurve. This typecurve can be used to generalize the permeability correction in tight porous media. We conclude that Knudsen’s permeability correlation is more accurate than Klinkenberg’s model especially for extremely tight porous media with transition and free molecular flow regimes. The results from this study indicate that Klinkenberg’s model and various extensions developed throughout the past years underestimate the permeability correction especially for the case of fluid flow with the high Knudsen number.     

纳米尺度下气体薄膜润滑理论研究

提出了依赖于逆Knudsen数的纳米尺度影响因子Np，修正了目前应用较为普遍的FK-Boltzmann模型，并给出了其数值计算结果。结果表明：当膜厚小于单个空气分子平均自由程时，由于纳米尺度效应的影响，流量因子有所减小，承载能力略微增加。     

Gasdynamic equations of higher approximations of the chapman-enskog method. Review

Some features of the derivation of the gasdynamic equations of higher (Burnett and super-Burnett) approximations of the Chapman-Enskog method are discussed for non-dense monatomic gases in the absence of external forces. The equations are given, and their properties and the results of applying them to various problems of gas dynamics are described. The problems of deriving, justifying, and testing macroscopic models are among the fundamental problems of not only continuum mechanics but physical kinetics in general. In the latter case it is a question of reducing more general molecular-kinetic models to simpler macroscopic ones. Historically, the classical Chapman-Enskog method of solution (as the Knudsen number Kn → 0) of the kinetic Boltzmann-Maxwell equation is of primary importance. Its application yields a chain of Euler and Navier-Stokes equations and equations of higher (Burnett, super-Burnett, etc.) approximations of the method. The publication of this review is a result of the increasing interest in the equations of higher approximations of the Chapman-Enskog method, in particular, in the Burnett equations. Primarily, this is a consequence of advances in the application of the latter equations to the problem of shock wave structure and the possible (by virtue thereof) expansion of the region of applicability of the macroscopic rarefied gas flow models. In addition, using the Burnett approximation, a series of new interesting effects has been established. In particular, it has been shown that in the classical problem of heat transfer between bodies heated to different temperatures it is necessary to take into account not only the temperature creep but also the Burnett temperature stresses. At the same time, a number of fundamental problems concerning the status of these equations has not yet been solved. An extensive literature is devoted to this problem. In the present paper, we shall mainly consider studies carried out since the end of the sixties and concerned with the best developed case of neutral monatomic (structureless) non-dense gas flow in the absence of external forces. Purely mathematical problems are not considered. Preference is given to final or review studies appearing in the most accessible publications. In Section 1 we consider the problem of taking out-of-order terms in the Chapman-Enskog method into account and make a brief comparison with other perturbation methods. Section 2 is devoted to the expressions for the transport properties in the Burnett and super-Burnett approximations. The relations between the Burnett approximation and axiomatic continuum mechanics theory, the thermodynamics of irreversible processes, and nonequilibrium statistical mechanics are discussed in Section 3. New effects due to specific terms of the formulas for the Burnett transport properties as Kn → 0 are considered in Section 4. The basic features of formulations of the boundary value problems for the systems of gasdynamic equations considered are outlined in Section 5. Finally, we describe the results of applying these equations to particular problems: some exact results of the kinetic theory (Section 6), the propagation of sound waves and cylindrical Couette flow (Section 7), shock wave structure and hypersonic hydrodynamics (Section 8). Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 3–28, July–August, 1998. The work was carried out with financial support from the Russian Foundation for Basic Research (project No. 96-01-01244).     

Weakly rarefied gas flows between nonuniformly heated parallel plates

The aim of this research is to establish the validity of the predictions of the theory of slow nonisothermal flows, to study the limits of applicability (with respect to the Knudsen number) of the conclusions reached and to determine the effect of the Knudsen layers on these flows on the basis of a numerical investigation of slow nonisothermal weakly rarefied gas flow in a plane infinite channel with weakly nonequilibrium heating of the walls and a finite wall temperature difference. The gas flow is described by a relaxation transport equation. The results obtained show how quickly, as the Knudsen number decreases, the solutions of the transport equation outside the Knudsen layers tend to the solution of the equations of gas dynamics of slow nonisothermal flows (and not to the solution of the Navier-Stokes equations).Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 115–121, January–February, 1988.     

Effects of Knudsen number and geometry on gaseous flow and heat transfer in a constricted microchannel

A flow and heat transfer numerical simulation is performed for a 2D laminar incompressible gas flow through a constricted microchannel in the slip regime with constant wall temperature. The effects of rarefaction, creeping flow, first order slip boundary conditions and hydrodynamically/thermally developing flow are assumed. The effects of Knudsen number and geometry on thermal and hydrodynamic characteristics of flow in a constricted microchannel are explored. SIMPLE algorithm in curvilinear coordinate is used to solve the governing equations including continuity, energy and momentum with the temperature jump and velocity slip conditions at the solid walls in discretized form. The resulting velocity and temperature profiles are then utilized to obtain the microchannel C _f Re and Nusselt number as a function of Knudsen number and geometry. The results show that Knudsen number has declining effect on the C _f Re and Nusselt number in the constricted microchannel. In addition, the temperature jump on wall and slip velocity increase with increasing Knudsen number. Moreover, by decreasing the throttle area, the fluid flow characteristics experience more intense variations in the constricted region. To verify the code a comparison is carried out with available results and good agreement is achieved.     

复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

Empirical slip and viscosity model performance for microscale gas flow

For the simple geometries of Couette and Poiseuille flows, the velocity profile maintains a similar shape from continuum to free molecular flow. Therefore, modifications to the fluid viscosity and slip boundary conditions can improve the continuum based Navier–Stokes solution in the non‐continuum non‐equilibrium regime. In this investigation, the optimal modifications are found by a linear least‐squares fit of the Navier–Stokes solution to the non‐equilibrium solution obtained using the direct simulation Monte Carlo (DSMC) method. Models are then constructed for the Knudsen number dependence of the viscosity correction and the slip model from a database of DSMC solutions for Couette and Poiseuille flows of argon and nitrogen gas, with Knudsen numbers ranging from 0.01 to 10. Finally, the accuracy of the models is measured for non‐equilibrium cases both in and outside the DSMC database. Flows outside the database include: combined Couette and Poiseuille flow, partial wall accommodation, helium gas, and non‐zero convective acceleration. The models reproduce the velocity profiles in the DSMC database within an L₂ error norm of 3% for Couette flows and 7% for Poiseuille flows. However, the errors in the model predictions outside the database are up to five times larger. Copyright © 2005 John Wiley & Sons, Ltd.     

Flow of a viscous fluid in an annular gap with intensive blowing from the walls

Self-similar solutions are obtained in [1, 2] to the Navier-Stokes equations in gaps with completely porous boundaries and with Reynolds number tending to infinity. Approximate asymptotic solutions are also known for the Navier-Stokes equations for plane and annular gaps in the neighborhood of the line of spreading of the flow [3, 4]. A number of authors [5–8] have discovered and studied the effect of increase in the stability of a laminar flow regime in channels of the type considered and a significant increase in the Reynolds number of the transition from the laminar regime to the turbulent in comparison with the flow in a pipe with impermeable walls. In the present study a numerical solution is given to the system of Navier-Stokes equations for plane and annular gaps with a single porous boundary in the neighborhood of the line of spreading of the flow on a section in which the values of the local Reynolds number definitely do not exceed the critical values [5–8]. Generalized dependences are obtained for the coefficients of friction and heat transfer on the impermeable boundary. A comparison is made between the solutions so obtained and the exact solutions to the boundary layer equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 21–24, January–February, 1987.     

Mass flow rate prediction of pressure–temperature-driven gas flows through micro/nanoscale channels

In this paper, we study mass flow rate of rarefied gas flow through micro/nanoscale channels under simultaneous thermal and pressure gradients using the direct simulation Monte Carlo (DSMC) method. We first compare our DSMC solutions for mass flow rate of pure temperature-driven flow with those of Boltzmann-Krook-Walender equation and Bhatnagar-Gross-Krook solutions. Then, we focus on pressure–temperature-driven flows. The effects of different parameters such as flow rarefaction, channel pressure ratio, wall temperature gradient and flow bulk temperature on the thermal mass flow rate of the pressure–temperature-driven flow are examined. Based on our analysis, we propose a correlated relation that expresses normalized mass flow rate increment due to thermal creep as a function of flow rarefaction, normalized wall temperature gradient and pressure ratio over a wide range of Knudsen number. We examine our predictive relation by simulation of pressure-driven flows under uniform wall heat flux (UWH) boundary condition. Walls under UWH condition have non-uniform temperature distribution, that is, thermal creep effects exist. Our investigation shows that developed analytical relation could predict mass flow rate of rarefied pressure-driven gas flows under UWH condition at early transition regime, that is, up to Knudsen numbers of 0.5.     

Bulk equations and Knudsen layers for the regularized 13 moment equations

The order of magnitude method offers an alternative to the Chapman-Enskog and Grad methods to derive macroscopic transport equations for rarefied gas flows. This method yields the regularized 13 moment equations (R13) and a generalization of Grad’s 13 moment equations for non-Maxwellian molecules. Both sets of equations are presented and discussed. Solutions of these systems of equations are considered for steady state Couette flow. The order of magnitude method is used to further reduce the generalized Grad equations to the non-linear bulk equations, which are of second order in the Knudsen number. Knudsen layers result from the linearized R13 equations, which are of the third order. Superpositions of bulk solutions and Knudsen layers show good agreement with DSMC calculations for Knudsen numbers up to 0.5.      

Some self-similar solutions of Grad's equations

It is shown that the self-similar solutions of the Navier-Stokes and Burnett equations found earlier by the authors [1–9] can be extended to the case of two-dimensional flows of a weakly rarefied gas described by Grad's equations. Examples are given of numerical realization of self-similar solutions for flow in an expanding planar channel. It is found that there are appreciable differences between the behavior of the self-similar solutions of the Navier-Stokes, Burnett, and Grad equations in the neighborhood of a channel wall.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 88–94, May–June, 1982.