Evaporation from a surface and vapor flow through a plane channel into a vacuum

Two-dimensional steady rarefied-gas channel flow between two parallel walls, from an evaporating face to a perfectly absorbing plane end face, is studied. The vapor is considered to be a monatomic gas. The corresponding problem for the kinetic equation with collision integral in BGK form is formulated and solved numerically by two different finite-difference methods. Attention is focused on the calculation of the total gas flow rate through the channel cross-section. The structure of the gas channel flow as a function of the flow rarefaction, the channel length, and the ratio of the evaporation temperature to the wall temperature is studied.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 150–158, January–February, 1996.     

Solution to problem of strong evaporation of a monatomic gas by the Monte Carlo method

The Monte Carlo method has been used to obtain a numerical solution to the problem of strong evaporation of a monatomic gas in which the molecules are modeled by pseudo-Maxwellian and hard spheres. A comparison with the results of other authors is made. The results agree well with the solution of the problem obtained on the basis of the model Bhatnagar—Gross—Krook kinetic equation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 185–188, January–February, 1984.I should like to thank M. N. Kogan for discussing the results.     

A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

Strong recondensation in a one-component and a two-component rarefied gas for an arbitrary value of Knudsen number

As is known, surface phenomena such as evaporation, absorption, and reflection of molecules from the surface of a body depend strongly on its temperature [1–5]. This leads to the establishment of a flow of a substance between two surfaces maintained at different temperatures (recondensation). The phenomenon of recondensation was studied in kinetic theory comparatively long ago. However, up to the present, only the case of small mass flows in a onecomponent gas has been investigated completely [3,4]. Meanwhile it is clear that by the creation of appropriate conditions we can obtain considerable flows of the recondensing substance, so that the mass-transfer rate will be of the order of the molecular thermal velocity. Such a numerical solution of the problem with strong mass flows along the normal to the surface for small Knudsen numbers for a model Boltzmann kinetic equation was obtained in [7]. In this study we numerically solve the problem of strong recondensation between two infinite parallel plates over a wide range of Knudsen numbers for a one-component and a two-component gas, on the basis of the model Boltzmann kinetic equation [6] for a one-component gas and the model Boltzmann kinetic equation for a binary mixture in the form assumed by Hamel [8], for a ratio of the plate temperatures equal to ten. We also investigate the effect of the relative plate motion on the recondensation flow.Moscow. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 130–138, September–October, 1972.     

Calculation by the Krylov-Bogolyubov method of the velocity profile and fluxes in the nonisothermal flow of a rarefied gas in a cylindrical capillary

The Krylov-Bogolyubov numerical method is used to solve integral transfer equations obtained from the kinetic equation with the BHC (Bhatnagar—Cross—Krook) model of the collision operator. The velocity profiles and the thermal-creep flows and Poiseuille flow are calculated in different modes of flow under conditions of incomplete accommodation of the tangential momentum of molecules at the wall.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 143–150, November–December, 1978.     

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.     

The steady-state flow of a rarefied gas from a spherical source or sink

The method of characteristics is used to solve problems in the steady-state flows of a rarefied gas on the basis of approximating the kinetic equations. Numerical results are given for the solution of the problem of the flow from a spherical source or sink using the generalized Kruk equation [1] and approximating the Boltzmann equation by the method proposed by the author [2, 3], Various flow conditions are discussed: flow into a vacuum, flow into a flooded volume, flow from a sink.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 58–66, March–April, 1971.     

A linearization of Mieussens's discrete velocity model for kinetic equations

A linearization is developed for Mieussens's discrete velocity model (see, e.g., [L. Mieussens, Discrete-velocity models and numerical schemes for the Boltzmann-BGK equation in plane and axisymmetric geometries, J. Comput. Phys. 162 (2000) 429–466]) for kinetic equations. The basic idea is to use a linearized expression of the reference distribution function in the kinetic equation, instead of its exact expression, in the numerical scheme. This modified scheme is applied to various kinetic models, which include the BGK model, the ES-BGK model, the BGK model with velocity-dependent collision frequency, and the recently proposed ES-BGK model with velocity-dependent collision frequency. One-dimensional stationary shock waves and stationary planar Couette flow, which are two benchmark problems for rarefied gas flows, are chosen as test examples. Molecules are modeled as Maxwell molecules and hard sphere molecules. It is found that results from the modified scheme are very similar to results from the original Mieussens's numerical scheme for various kinetic equations in almost all tests we did, while, depending on the test case, 20–40 percent of computational time can be saved. The application of the method is not affected by the Knudsen number and molecular models, but is restricted to lower Mach numbers for the BGK (or the ES-BGK) model with velocity-dependent collision frequency.     

Numerical simulation of compressible turbulent flow via improved gas‐kinetic BGK scheme

In order to solve compressible turbulent flow problems, this study focuses on incorporating the Spalart–Allmaras turbulence model into gas‐kinetic BGK (Bhatnagar–Gross–Krook) scheme. The Spalart–Allmaras turbulence model is solved using finite difference discretization. The variables on the cell interface are interpolated via the van Leer limiter in the reconstruction stage. Simulation of subsonic and transonic flow over a NACA0012 airfoil has been implemented using two‐dimensional body‐fitted grids. The numerical results obtained appear in good agreement with the AGARD results, demonstrating the effectiveness and usefulness of the strategy of coupling the Spalart–Allmaras turbulence model with the BGK scheme for compressible turbulent flow simulation. Copyright © 2010 John Wiley & Sons, Ltd.     

The noneouilibrium dissociation of diatomic molecules in a boundary layer

A study is made of the influence of the process of diffusion of diatomic molecules on the rate of their thermal dissociation. On the basis of the ladder model of the dissociation of molecules — truncated harmonic oscillators — an expression is obtained for the macroscopic reaction rate which depends explicitly and exponentially on the square of the gradient of the gas temperature. The molecules constitute a small admixture in a monatomic inert gas. The effect of the diffusion on the rate of dissociation is illustrated by the results of numerical calculation of the flow in the boundary layer around the stagnation point of a blunt body.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Shidkosti i Gaza, No. 1, pp. 161–170, January–February, 1985.     

Numerical investigation from rarefied flow to continuum by solving the Boltzmann model equation

Based on the Bhatnagar–Gross–Krook (BGK) Boltzmann model equation, the unified simplified velocity distribution function equation adapted to various flow regimes can be presented. The reduced velocity distribution functions and the discrete velocity ordinate method are developed and applied to remove the velocity space dependency of the distribution function, and then the distribution function equations will be cast into hyperbolic conservation laws form with non‐linear source terms. Based on the unsteady time‐splitting technique and the non‐oscillatory, containing no free parameters, and dissipative (NND) finite‐difference method, the gas kinetic finite‐difference second‐order scheme is constructed for the computation of the discrete velocity distribution functions. The discrete velocity numerical quadrature methods are developed to evaluate the macroscopic flow parameters at each point in the physical space. As a result, a unified simplified gas kinetic algorithm for the gas dynamical problems from various flow regimes is developed. To test the reliability of the present numerical method, the one‐dimensional shock‐tube problems and the flows past two‐dimensional circular cylinder with various Knudsen numbers are simulated. The computations of the related flows indicate that both high resolution of the flow fields and good qualitative agreement with the theoretical, DSMC and experimental results can be obtained. Copyright © 2003 John Wiley & Sons, Ltd.     

Influence of real gas properties on the integral aerodynamic coefficients

The results of mathematically modeling axisymmetric hypersonic flow past an ellipsoid are presented. The calculation data are obtained on the basis of a numerical solution of the complete Navier—Stokes equations using a finite-difference method. The investigation is carried out for the problem of laminar flow over the windward side of the body. The effect of the elongation of the ellipsoid on the total heat flux and the viscous and pressure drag is considered. Results are obtained for three different gas (air) models: the perfect gas, chemical equilibrium and chemical nonequilibrium models. For the latter model various sets of catalyticity coefficients are considered. The effect of the real properties of air on the integral aerodynamic characteristics of ellipsoids with different elongations is analyzed.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 156–164, July–August, 1992.     

Blowing of a foreign gas in a hypersonic viscous shock layer

We consider the problem of a hypersonic viscous flow of a nonreactive mixture of ideal gases around smooth thick bodies in the framework of a two-layer model of a thin shock layer for moderately small Reynolds numbers. We investigate the effect of blowing of a foreign gas through a permeable surface in the bow region of a spherical blunt body. We introduce a transformation of variables that gives a number of important advantages in the numerical solution of the problem under consideration. The problem of mass blowing from the surface of a body into a boundary layer has an extensive literature. The effect of blowing for moderately small Reynolds numbers has been considerably less studied [1–5], and in the majority of papers on this question either the critical point of a blunt body or the blowing of a gas homogeneous with the gas in the incoming flow is investigated.Moscow. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 110–116, July–August, 1972.     

Rarefied gas flow past a sphere with injection

Hypersonic rarefied gas flow over the windward face of a sphere is considered in the presence of distributed injection from the surface of the body. A similar problem was previously solved in [1–3] within the framework of continuum mechanics and in [4] on the basis of model kinetic equations. In the present study the calculations were carried out using the Monte Carlo method of direct statistical modeling [5, 6]. The injected gas was the same as the free-stream gas. A simple monatomic gas model with a rigid sphere interaction potential was employed. The reflection of the molecules from the surface of the body was assumed to be diffuse with total energy accommodation. The calculation procedure using weighting factors is described in [7]. The influence of injection on the mechanical and thermal effect of the gas flow on the body is investigated for various degrees of rarefaction of the medium and injection rates.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 175–179, July–August, 1990.     

Approximate solution of the problem of the bottom zone in a rarefied gas

The problem of the steady-plane monatomic rarefied gas flow around a semiinfinite bar is considered (the plane stationary case of the problem about the bottom zone). The problem is solved numerically at the level of the Krook relaxation model [1, 2]. A depenence of the gas density, velocity, and temperature in the whole flow domain on the space coordinates is obtained for supersonic and subsonic gas streams flowing around a body by using calculations on an M-20 electronic calculator.Khar'kov. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 139–143, January–February, 1972.     

Boundary value problems for a model Boltzmann equation with frequency proportional to the molecule velocity

Exact solutions of a model Boltzmann equation with a collision frequency that depends on the molecule velocity and with a BGK (Bhatnagar-Gross-Krook) collision operator are constructed for the problems of weak evaporation and temperature jump in a rarefied vapor above a plane surface. The numerical calculations and a comparison with previous results are given.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 140–153, May–June, 1996.     

Eddy shocklet in coplanar molecular gas flow

The purpose of this study is to investigate the feasibility of applying a kinetic molecular model to the problem of turbulence-oriented computation of compressible flow. Consider a simple two-dimensional initial value problem of the Taylor-Green-type periodic flow in a square region. The Boltzmann equation with its collision term of the BGK model is used in its integral form along the characteristics. An approximation based on small time steps is utilized for actual computation. It takes several times more steps than the corresponding Euler computation. The results show in general that basically the same pattern for flow fields in the early stages is obtained independently by the Euler and the BGK models, with less violent motion with the BGK.     

Properties of the nonequilibrium high-velocity “Tails” of the molecular distribution function in plane couette flow of a compressible gas

The results of an analytic and numerical investigation of the properties of the high-velocity “tails” of the distribution function are given for the solution of the BGK model of the kinetic Boltzmann equation for plane Couette flow of a compressible gas. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 183–190, July–August, 1998. The work was carried out with financial support from the Russian Foundation for Basic Research (project No. 96-01-00573; grant in support of leading science schools No. 96-15-9603).     

Jumps in the density and temperature of one of the components of a binary gas mixture over the plane surface of a liquid

We consider a problem concerning the vaporization (or condensation) of one of the components of a binary gas mixture situated over the plane surface of a liquid. The kinetic equation in the model form of [1] is used to describe the system. As is well known, this model agrees well with experiment and is simpler than the Boltzmann equation so far as mathematical relations are concerned. This model fails to describe a number of effects because it is assumed that the collision time of the particles is independent of their velocity. This relates primarily to the phenomenon of thermal diffusion of the gases. Thus the formulas given below are applicable to gas mixtures having a small thermal-diffusion coefficient. The model equation is solved by an approximate method developed in [2]. In [3] the temperature jump of a single-component gas at a solid wall is calculated by this method, and the method is also employed in [4] in the calculation of the slipping rate of a binary gas mixture in the field of a temperature gradient. In both cases the results agree with numerical calculations by other authors within an accuracy of 1.5%.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 142–148, September–October, 1973.