Onsager's reciprocal relations in the linear problem of flow of rarefied gas past a heat conducting body

Symmetry of integral transport coefficients is established on the basis of a linear stationary Boltzmann equation for the problem of flow past a heat conducting body. An expression is obtained for the entropy production in the gas-body system, and this determines the thermodynamic fluxes and forces. Bakanov and Roldugin [2] have considered the problem of motion of a heat conducting sphere at small Knusden numbers Kn, using the symmetry of the Onsager coefficients to construct an asymptotic solution as Kn → 0. In the present paper, a general method is proposed for establishing the Onsager relations that does not require the actual construction of a solution to the problem and is applicable for bodies of arbitrary shape and all values of the Knudsen number.     

Forced convection of gaseous slip-flow in porous micro-channels under Local Thermal Non-Equilibrium conditions

Steady laminar forced convection gaseous slip-flow through parallel-plates micro-channel filled with porous medium under Local Thermal Non-Equilibrium (LTNE) condition is studied numerically. We consider incompressible Newtonian gas flow, which is hydrodynamically fully developed while thermally is developing. The Darcy–Brinkman–Forchheimer model embedded in the Navier–Stokes equations is used to model the flow within the porous domain. The present study reports the effect of several operating parameters on velocity slip and temperature jump at the wall. Mainly, the current study demonstrates the effects of: Knudsen number (Kn), Darcy number (Da), Forchheimer number (Γ), Peclet number (Pe), Biot number (Bi), and effective thermal conductivity ratio (K _R) on velocity slip and temperature jump at the wall. Results are given in terms of skin friction (C _f Re ^*) and Nusselt number (Nu). It is found that the skin friction: (1) increases as Darcy number increases; (2) decreases as Forchheimer number or Knudsen number increases. Heat transfer is found to (1) decreases as the Knudsen number, Forchheimer number, or K _R increases; (2) increases as the Peclet number, Darcy number, or Biot number increases.     

Slip coefficients and macroparameter jumps of a diatomic gas with rotational degrees of freedom on a weakly curved spherical surface

Using the half-space moment method, the problem of the slip of a diatomic gas along a rigid spherical surface is solved within the framework of a model kinetic equation previously proposed which takes into account the rotational degrees of freedom of the gas. Second-order slip coefficients (correctionsC _m ^′ , β _R ^′ , and β _R to the isothermal and thermal slip which are linear with respect to the Knudsen number Kn) are obtained. The gas macroparameter jump coefficientsC _v andC _q, which are of the second order in the Knudsen number and characterize the discontinuity of the normal mass and heat fluxes on the gas-rigid phase interface, are calculated. These coefficients are given as functions of the tangential momentum accommodation coefficient, the translational and rotational energy accommodation coefficients, and the Prandtl number. The coefficients are calculated for certain diatomic gases. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 163–173, January–February, 2000.     

Effect of nonuniform accommodation coefficient distribution on the couette flow

The Couette flow is considered for surfaces with nonuniformly distributed energy accommodation coefficients α. It is shown that at Knudsen numbers greater or of the order of unity heat fluxes and viscous stresses can be considerably optimized by varying the surface distribution of α at a fixed integral value. At the same time, for Kn ≪ 1 the flows with nonuniformly distributed α are similar with the flow with a constant accommodation coefficient equal to its mean value.     

Asymptotic Solutions of the Equations for a Viscous Heat-Conducting Compressible Medium

Small nonstationary perturbations in a viscous heat-conducting compressible medium are analyzed on the basis of the linearization of the complete system of hydrodynamic equations for small Knudsen numbers (Kn ≪ 1). It is shown that the density and temperature perturbations (elastic perturbations) satisfy the same wave equation which is an asymptotic limit of the hydrodynamic equations far from the inhomogeneity regions of the medium (rigid, elastic or fluid boundaries) as M _a = v/a → 0, where v is the perturbed velocity and a is the adiabatic speed of sound. The solutions of the new equation satisfy the first and second laws of thermodynamics and are valid up to the frequencies determined by the applicability limits of continuum models. Fundamental solutions of the equation are obtained and analyzed. The boundary conditions are formulated and the problem of the interaction of a spherical elastic harmonic wave with an infinite flat surface is solved. Important physical effects which cannot be described within the framework of the ideal fluid model are discussed.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, 2005, pp. 76–87.Original Russian Text Copyright © 2005 by Stolyarov.     

A kinetic theory of gas mixtures with regard to the laws of the thermodynamics of irreversible processes

The derivation of the hydrodynamic equations for a gaseous mixture from the system of kinetic Boltzmann equations is analyzed. The form of the hydrodynamic equations is a unique consequence of necessary and sufficient conditions for the solvability of systems of linear integral equations with symmetrical kernels, which define the terms in the expansion of the distribution functions in a series with respect to a parameter of spatial non-homogeneity (actually, the Knudsen number). The transport laws are presented in a form for which the Onsager symmetry relations hold. In deriving the Onsager relations use is made of symmetry properties of integral operators, which are a consequence of the invariance of the equations of mechanics with respect to a transformation involving changing the sign of the time and the impulses of the particles. The Onsager relations are also derived from expressions for the kinetic coefficients in terms of correlation functions.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 44–50, July–August, 1970.We wish to express our thanks to V. V. Struminskii for useful discussions of this work.     

Slow gas motions and the negative drag of a strongly heated spherical particle

In the slow flows of a strongly and nonuniformly heated gas, in the continuum regime (Kn → 0) thermal stresses may be present. The theory of slow nonisothermal continuum gas flows with account for thermal stresses was developed in 1969–1974. The action of the thermal stresses on the gas results in certain paradoxical effects, including the reversal of the direction of the force exerted on a spherical particle in Stokes flow. The propulsion force effect is manifested at large but finite temperature differences between the particle and the gas. This study is devoted to the thermal-stress effect on the drag of a strongly heated spherical particle traveling slowly in a gas for small Knudsen numbers (M ～ Kn → 0), small but finite Reynolds numbers (Re ≤ 1), a linear temperature dependence of the transport coefficients µ ∝ T, and large but finite temperature differences ((T _w ? T _∞ )/T _M8 ～ 1). Two different systems of equations are solved numerically: the simplified Navier-Stokes equations and the modified Navier-Stokes equations with account for the thermal stresses.     

Rarefied gas flow in a channel with specular‐diffuse boundary conditions for arbitrary knudsen numbers. Onsager relations

The problem of a rarefied gas flow in a channel for arbitrary Knudsen numbers has been solved analytically for the first time in the case where the scattering of gas molecules on the channel walls can be described by speculardiffuse boundary conditions. The mean free path of gas molecules is assumed to be constant, i.e., the collision frequency is proportional to molecular velocity. The gas moves under the action of a streamwise temperature gradient. Exact relations for heat and mass fluxes and for meanmass velocity are obtained. It is shown that the Onsager relations are valid within the entire range of Knudsen numbers in the problem of heat and mass transfer in a channel. The dependence of heat and mass fluxes on the Knudsen number (channel thickness) is analyzed. A comparison with available results is performed.     

Effect of the mode of gas-surface interaction on intensive monatomic gas evaporation

The direct Monte Carlo simulation method is used for investigating the effect of the thermal accommodation coefficient α _E on the relations on the Knudsen layer boundary in the presence of intensive evaporation. The model of mirror reflection of molecules from the surface is considered. It is shown that diffuse reflection with α _E = 0 leads to almost the same relations on the Knudsen layer boundary as mirror reflection. The accuracy of the moment method is estimated in application to the problems of intensive evaporation with diffuse and mirror reflection from the surface.     

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.     

Boundary kinetic effects in the problem of slow flow over an evaporating heat-conducting spherical particle

It is shown that for sufficiently large values of the thermal conductivity of the condensed phase' as compared with the thermal conductivity of the vapor (/' Kn) the effects associated with the presence of a Knudsen layer on the evaporating surface must be taken into account in order to obtain a solution of the problem of a spherical particle in a slow (Re, 1) continuum (Knudsen number Kn 1) flow of its own vapor. The drag is calculated for various types of boundary conditions on the particle surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 134–139, March–April, 1987.In conclusion the authors wish to thank V. S. Galkin and M. N. Kogan for useful discussions.     

Propagation of Elastic Waves in a Gas-Filled Poroelastic Medium: The Influence of the Boundary Conditions

We consider the propagation of elastic waves in gas-filled porous media at small but non-zero values of Knudsen numbers \( {\text{Kn}} \), where \( {\text{Kn}} = \lambda /l \), \( \lambda \) is the mean free path of gas molecules; \( l \) is the characteristic size of inclusion (the so-called slip regime). In this case, it is possible to apply the classic equations of hydrodynamics with modified boundary conditions at solid walls. We have assumed that the gas molecules distribution function is satisfied at the modified Maxwell boundary conditions (Struchtrup 2013; Mohammadzadeh and Struchtrup 2015). We have obtained the expressions for drag and added mass coefficients for the Biot equations of poroelasticity for a system of randomly oriented gas-filled cylindrical capillaries. Our calculations have shown that the drag and added mass coefficients depend considerably on the Knudsen number and the properties of the surface. The influence of the interfacial slip effect on the velocities of the compressional wave of the first kind and shear wave is small, but the velocity and attenuation of the compressional wave of the second kind are considerably influenced by this effect. The results obtained show the fundamental possibility of the determination of the accommodation coefficient by measuring the velocity of the compressional wave of the second kind for different values of the Knudsen number.     

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.     

Effect of the Nature of the Interaction between a Gas and a Surface on the Relations for the Knudsen Layer in the Case of Strong Evaporation

The effect of the temperature accommodation coefficient _T on the relations at the Knudsen layer edge is investigated for strong evaporation using the moment method. An explicit expression for the dimensionless density as a function of the temperature and the Mach number M is obtained for 0 < _T < 1. For _T = 0 the entire solution is obtained in explicit form. It is shown that for = 0 and a condensation coefficient << 1 the temperature outside the Knudsen layer changes sharply as M varies from 0 to a certain value much less than unity after which the temperature ceases to depend on . For the model of specular reflection of the molecules from the surface the density and the temperature outside the Knudsen layer are found in explicit form as functions of the Mach number.     

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.     

Stochastic Effects of 2D Random Arrays of Cylinders on Rarefied Gas Permeability Using the Lattice Boltzmann Method

A large set of 2D random arrays of circular cylinders is generated to perform a statistical study on rarefied gas flow through micro-porous media. The flow regimes in this work lie for Knudsen numbers (Kn) ranging from the continuum to the transition regimes. Arrays are built by randomly placing cylinders with constant diameter with a uniform distribution without overlapping, and are generated for three target porosities. Fluid flow is assumed to be incompressible and isothermal. A modified lattice Boltzmann model is adopted to account for discrete effects, with slip-velocity boundary conditions and a Kn-dependent multi-relaxation time collision operator. The apparent permeability is modeled with Darcy’s law with a Klinkenberg-type relationship and compared with existing correlations. Velocity fields highlight the increasing contribution of fluid flow through small pores with increasing Kn. Numerical results show that porous media randomness leads to an uncertainty on rarefied gas permeability calculation despite the same structural characteristics and may not strictly follow a specific correlation. The influence of a local collision operator based on a local Kn instead of a global one in the numerical model is also studied. Results show that the permeability in rarefied regimes undergoes significant deviation when applying the local collision operator compared to the global one. These differences could result from a more accurate capture of the pore-scale behavior with a local Kn. Thus, it emphasizes the sensitivity of the model and the apparent permeability calculation to the appropriate definition of Kn.

     

Study of stable and unstable jet flows on the basis of the boltzmann equation

Free supersonic underexpanded jets are studied using a direct method conservative splitting scheme for solving the Boltzmann equation. Numerical solutions for a jet flowing into a vacuum and into a fluid-filled space are presented for the following ranges of the parameters: Knudsen number 10⁻⁶<Kn<∞ and pressure ratio 10<n<∞. The solutions are compared with experimental data. Instabilities associated with free turbulence effects in the mixing layer are detected for low Kn numbers. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 153–157, March–April, 1998. The work was carried out with support from the Russian Foundation for Fundamental Research (project No. 96-01-00829).     

Vapor flows with evaporation and condensation in the continuum limit: effect of a trace of noncondensable gas

Steady flows of a vapor with evaporation and condensation on the boundary consisting of the condensed phase of the vapor are considered in the following situation: (i) the boundary is of arbitrary smooth shape; (ii) the Knudsen number Kn, the ratio of the typical mean free path of the vapor molecules to the characteristic length of the system, is small; (iii) a small amount of a noncondensable gas is contained in the system; more specifically, the amount is such that the average concentration of the noncondensable gas is of the order of Kn in the case of a closed domain (the case of an infinite domain is also discussed). The steady behavior of the vapor and the noncondensable gas, in particular, that in the continuum limit where Kn vanishes, is investigated by means of a systematic asymptotic analysis based on kinetic theory. In this situation, the average concentration of the noncondensable gas becomes infinitely small in the continuum limit in the case of a closed domain. However, it is shown that the noncondensable gas accumulates in the infinitely thin Knudsen layer on the boundary where condensation is taking place and has a significant effect on the global vapor flow in the continuum limit. An example demonstrating such an effect is also given.     

Boundary conditions for flows with chemical reactions occurring on a surface

With the use of a solution of a model Boltzmann equation for a binary mixture in the Knudsen layer, we obtain the boundary conditions for the equations of gas dynamics when the reactionl _iA_il _jA_j (l _i molecules of A_i change intol _j molecules of A_j, and vice versa) is occurring on a surface. The boundary condition that we obtain differs from those that are usually applicable by the presence of terms of the same order. This confirms the conclusion arrived at by the authors in [1], where it was shown that if the Knudsen layer is left out of account, which is precisely what is usually done, it is impossible to obtain correct boundary conditions.Moscow. Translated from Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 129–138, January–February, 1972.