Effects of nonmonotonicity and change in sign of the energy flux in the transient regime for the Couette problem with heat transfer

On the basis of the numerical solution of the Boltzmann equation the Couette flow is investigated. The effect of the upper plate velocity on the value and position of the extremum of the dependence of the energy flux to the lower plate on the plate temperature ratio in the free-molecular regime is considered. It is established that in the case of sufficiently small surface temperature ratio the dependence of the energy flux to the lower plate on the Knudsen number Kn has a local extremum. It is found that there is a plate temperature ratio range on which with increase in Kn the energy flux to the lower plate reverses direction. The Kn dependence of its absolute value then has a local maximum.     

Couette problem for the generalized Krook equation stress-peak effect

The Couette problem is the simplest problem of steady shear flow of rarefied gas in a region bounded by solid surfaces. This problem has been examined in the linear formulation by many authors, using either the linearized Krook equation or the moment methods (see [1]). It has recently been solved by the Monte Carlo method [2].The nonlinear problem of Couette flow with heat transfer for the Krook equation has been solved by reducing the problem to a system of integral equations [3] over a wide range of flat-plate velocities and temperature ratios and by the discrete-velocity method [4] for moderate plate velocities. In this article we solve the same problem for the generalized Krook equation [5] which approximates the Boltzmann equation for a pseudo-Maxwellian gas in accordance with the method suggested by the author [6, 7]. The generalized Krook equation was solved numerically by a modified discrete-velocity method which has been used by the author previously to solve the problem of shock wave structure [8].The primary case examined is that of pseudo-Maxwellian molecules, in which the viscosity is proportional to the temperature. The computations were made for Prandtl numbers of 1 and 2/3 over a wide range of Mach and Knudsen numbers as well as flat-plate temperature ratios. As we would expect, the Prandtl number effect is greatest for small Knudsen numbers. The flow velocity profiles are not very sensitive to variation of the Prandtl number (at least for pseudo-Maxwellian molecules).However, the most interesting result of the study is independent of the Prandtl number. Specifically, it was found that for any sufficiently high flat-plate velocities the friction stress, referred to the corresponding free molecular value, does not change monotonically with variation of the Knudsen number; instead, there is a peak. As far as the author is aware, this nonlinear effect has not been discussed previously in the literature (including articles [3, 4]).     

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.     

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

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

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.     

Flow instability of nanofuilds in jet

The flow instability of nanofluids in a jet is studied numerically under various shape factors of the velocity profile, Reynolds numbers, nanoparticle mass loadings,Knudsen numbers, and Stokes numbers. The numerical results are compared with the available theoretical results for validation. The results show that the presence of nanoparticles enhances the flow stability, and there exists a critical particle mass loading beyond which the flow is stable. As the shape factor of the velocity profile and the Reynolds number increase, the flow becomes more unstable. However, the flow becomes more stable with the increase of the particle mass loading. The wavenumber corresponding to the maximum of wave amplification becomes large with the increase of the shape factor of the velocity profile, and with the decrease of the particle mass loading and the Reynolds number. The variations of wave amplification with the Stokes number and the Knudsen number are not monotonic increasing or decreasing, and there exists a critical Stokes number and a Knudsen number with which the flow is relatively stable and most unstable,respectively, when other parameters remain unchanged. The perturbation with the first azimuthal mode makes the flow unstable more easily than that with the axisymmetric azimuthal mode. The wavenumbers corresponding to the maximum of wave amplification are more concentrated for the perturbation with the axisymmetric azimuthal mode.     

Aerodynamic coefficients of a spinning sphere in a rarefied-gas flow

A three-dimensional rarefied-gas flow past a spinning sphere in the transitional and near-continuum flow regimes is studied numerically. The rarefaction and compressibility effects on the lateral (Magnus) force and the aerodynamic torque exerted on the sphere are investigated for the first time. The coefficients of the drag force, the Magnus force, and the aerodynamic torque are found for Mach numbers ranging from 0.1 to 2 and Knudsen numbers ranging from 0.05 to 20. In the transitional regime, at a certain Knudsen number depending on the Mach number the Magnus force direction changes. This change is attributable to the increase in the role of normal stresses and the decrease in the contribution of the shear stresses to the Magnus force with decrease in the Knudsen number. A semi-empirical formula for the calculation of the Magnus force coefficient in the transitional flow regime is proposed.     

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.     

Direct Statistical Simulation of a Supersonic Flow of a Binary Mixture of Rarefied Gases around a Transversely Positioned Cylinder

A supersonic flow of a binary mixture of gases in a wide range of rarefaction (from a flow with a Knudsen number K _n = 0.1 to a free-molecular flow) around a cylinder is studied by means of direct statistical Monte Carlo simulations (DSMC method). The influence of a small fraction of heavy particles in a light gas flow on the region of significant nonequilibrium near the cylinder and on the heat flux is considered.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 5, pp. 53–59, September–October, 2005.     

Burnett simulations of gas flow in microchannels

The Burnett equations with slip boundary conditions are used to model the gas flow in microchannels in transition flow regime. As the Navier-Stokes equations are not appropriate to model the gas flow in this regime, the higher-order Burnett equations are adopted in the present study. In earlier studies, convergent solutions of the Burnett equations of microPoiseuille flow could only be obtained when Knudsen number is less than 0.2. By using a relaxation method on the boundary values, convergent solutions of the Burnett equations can be obtained even when Knudsen number reaches 0.4. The solutions of Burnett equations agree very well with experimental data and direct simulation Monte Carlo (DSMC) results. The pressure distributions and velocity profiles are then discussed in detail.     

Supersonic Rarefied Flow Past a Cascade of Transverse Flat Plates

The results of numerical calculations of the supersonic rarefied flow through an infinite, periodic cascade of flat plates set transverse to the incident flow are presented. The flow in the vicinity of the cascade is described on the basis of the kinetic Boltzmann equation. The Knudsen numbers, based on the plate span, the distances between the plates and the scalelength of the flow under consideration, range from 0.2 to 0.003. Both steady-state regimes with a shock attached to the cascade and time-dependent ones with an upstream-traveling shock are investigated.     

Numerical study of the strong evaporation of a binary mixture

The problem of the one-dimensional evaporation of a binary mixture is investigated by numerically solving a system of two coupled Boltzmann equations. The numerical method is based on the direct discretization of the Boltzmann equation and the Monte Carlo evaluation of the collision integrals. It is assumed that the fluid flows between an evaporating plate and a totally absorbing plate. The spatial profiles of macroscopic quantities as well as evaporation rates have been calculated for values of the Knudsen number between 1/2 and 1/20.     

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.     

多场作用下输流单层碳纳米管的动力学特性

以非局部弹性理论为基础,采用欧拉-伯努利梁模型,考虑管型区域内滑移边界条件以及碳纳米管的小尺度效应,应用哈密顿原理获得了温度场与轴向磁场共同作用下的输流单层固支碳纳米管（SWCNT）的振动控制方程以及边界条件,依靠微分变换法（DTM法）对此高阶偏微分方程进行求解,通过数值计算研究了多场中单层固支输流碳纳米管的振动与失稳问题。结果表明：温度场、轴向磁场强度、Knudsen数及小尺度参数都会对系统振动频率以及失稳临界流速产生影响。     

Direct Monte Carlo Simulation of Transverse Supersonic Rarefied Gas Flow around a Cylinder

Supersonic flow around a cylinder is investigated using the direct simulation Monte Carlo method over a wide rarefaction range: from the Knudsen number Kn = 0.1 to free-molecular flow. The effect of the cylinder temperature on the region of sharp nonequilibrium near the cylinder and the heat flux is studied.     

Numerical study of the generalized cylindrical Couette flow of rarefied gas

We study the cylindrical Couette flow of a rarefied gas between two cylinders in the generalized setup in which the inner of which not only rotates but also slides along its axis. The analysis is based on the numerical solution of the S-model kinetic equation. The influence of ratio of cylinder radiuses, velocities of the inner cylinder and Knudsen number on shear stresses, mass-flow rates as well as macroscopic parameters is investigated in the broad range of Knudsen numbers.     

Density distribution in rarefied gas flow behind a perforated plate

The relative density distribution in the rarefied CO₂ flow field behind a perforated plate has been experimentally investigated on the range of Knudsen numbers 10^–3 Kn 10^–1. The results of these experiments serve to demonstrate the validity of using the molecular velocity distribution function for determining the relative density on the Knudsen number interval in question. This distribution function was previously employed for calculating the parameters in molecular beams formed by capillary sieves in the free-molecular flow regime [1].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 183–184, May–June, 1989.     

Forced convection gaseous slip flow in circular porous micro-channels

Laminar forced convection of gaseous slip flow in a circular micro-channel filled with porous media under local thermal equilibrium condition is studied numerically using the finite difference technique. Hydrodynamically fully developed flow is considered and the Darcy–Brinkman–Forchheimer model is used to model the flow inside the porous domain. The present study reports the effect of several operating parameters (Knudsen number (Kn), Darcy number (Da), Forchhiemer number (Γ), and modified Reynolds number ) on the velocity slip and temperature jump at the wall. Results are given in terms of the velocity distribution, temperature distribution, skin friction , and the Nusselt number (Nu). It is found that the skin friction is increased by (1) decreasing Knudsen number, (2) increasing Darcy number, and (3) decreasing Forchheimer number. Heat transfer is found to (1) decrease as the Knudsen number, or Forchheimer number increase, (2) increase as the Peclet number or Darcy number increase.     

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.