Solution of the temperature jump (Knudsen layer flow) and linear heat-transfer problems for two parallel plates in a rarefied gas

The Monte Carlo method [1, 2] is used to solve the linearized Boltzmann equation for the problem of heat transfer between parallel plates with a wall temperature jump (Knudsen layer flow). The linear Couette problem can be separated into two problems: the problem of pure shear and the problem of heat transfer between two parallel plates. The Knudsen layer problem is also linear [3] and, like the Couette problem, can be separated into the velocity slip and temperature jump problems. The problems of pure shear and velocity slip have been examined in [2].The temperature jump problem was examined in [4] for a model Boltzmann equation. For the linearized Boltzmann equation the problems noted above have been solved either by expanding the distribution function in orthogonal polynomials [5–7], which yields satisfactory results for small Knudsen numbers, or by the method of moments, with an approximation for the distribution function selected from physical considerations in the form of polynomials [8–10]. The solution presented below does not require any assumptions on the form of the distribution function.The concrete calculations were made for a molecular model that we call the Maxwell sphere model. It is assumed that the molecules collide like hard elastic spheres whose sections are inversely proportional to the relative velocity of the colliding molecules. A gas of these molecules is close to Maxwellian or to a gas consisting of pseudo-Maxwell molecules [3].     

Flow stability in a channel with porous walls

We study the stability of the flow which forms in a plane channel with influx of an incompressible viscous fluid through its porous parallel walls. Under certain assumptions the study of the stability reduces to the solution of modified Orr-Sommerfeld equation accounting for the transverse component of the main-flow velocity. As a result of numerical integration of this equation we find the dependence of the local critical Reynolds number on the blowing Reynolds number R₀, which may be defined by two factors: the variation of the longitudinal velocity profile with R₀ and the presence of the transverse velocity component. A qualitative comparison is made of the computational results with experimental data on transition from laminar to turbulent flow regimes in channels with porous walls, which confirms that it is necessary to take into account the effect of the transverse component of the main-flow velocity on the main-flow stability in the problem in question.Flows in channels with porous walls are of interest for hydrodynamic stability theory in view of the fact that they can be described by the exact solutions of the Navier-Stokes equations by analogy with the known Poiseuille and Couette flows. However, in contrast with the latter, the flows in channels with porous walls (studies in [1], for example) will be nonparallel.The theory of hydrodynamic stability of parallel flows has frequently been applied to nonparallel flows (in the boundary layer, for example). In so doing the nonparallel nature of the flow has been taken into account only by varying the longitudinal velocity component profiles. A study was made in [2, 3] of the effect of the transverse component of the main flow on its stability. In the case of the boundary layer in a compressible gas, a considerable influence of the transverse velocity component on the critical Reynolds number was found in [2] and confirmed experimentally. A strong influence of the transverse velocity component on the instability region was also found in [3] in a study of the flow stability in a boundary layer with suction for an incompressible fluid.     

Diffusion and thermal slip of a binary gas mixture

Let us note that the phenomenon of diffusion slip at a constant gas-mixture temperature has been considered in [1], for example, and thermal slip for a single-component gas in [2]. The slip velocity of a binary gas mixture has been calculated in a field of the temperature gradient and of the partial pressure gradients. The kinetic equation is solved by an approximate method based on physical considerations. A formula has been obtained analytically for the slip velocity for arbitrary accommodation coefficients as well as for arbitrary gas concentrations and arbitrary molecule masses. The results agree to 1% accuracy with the numerical computations of other authors.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 51–55, July–August, 1970.     

Influence of hall effect and ion slip on velocity and temperature fields in an MHD channel

The combined influence of viscosity, Hall effect and ion slip on hydrodynamic fields and on heat transfer is investigated. The exact solutions for velocity, induced magnetic field and temperature are derived for the laminar MHD flow in a flat channel assuming a small magnetic Reynolds number, finely segmented electrodes, fully developed flow and uniform heat flux at channel walls. The internal generation of heat is not considered. The Kantorowitsch method of variational calculus is employed to approximate the complicated velocity distribution.     

Effects of slip velocity and surface acoustic wave on the laminar flow stability

The stability of slip flows when a surface acoustic wave (SAW) propagating along the walls of a microchannel in the laminar flow regime is investigated. The governing equation which was derived by considering the weakly nonlinear coupling between the deformable wall and streaming slip flow is linearized and then the eigenvalue problem is solved by a numerical code together with the associated interface and slip velocity boundary conditions. The value of the critical Reynolds number was found to be near 1,441 for a Knudsen number being 0.001 (associated with a physical parameter K ₀ characterizing the SAW effect) which is much smaller than the static-wall case for conventional pressure-driven flows.     

Kinetic theory of recondensation in a binary gas mixture with arbitrary Knudsen numbers

A rigorous solution of the problem of recondensation between two surfaces with arbitrary Knudsen numbers is possible only on the basis of a consecutive kinetic consideration. For the single-component case, this problem was solved in [1] using the BGK model of the collision integral in the kinetic equation. In [2], for the same purpose, the method of moments for Maxwellian molecules was used. The case of a binary mixture, in which one of the components is a noncondensing gas was discussed in [3, 4]. Under these circumstances, in [3], a single-relaxation lumped model was used for each component; the model did not reflect many of the properties of the exact collision integral. A more rigorous model (the collision integral in the Hamel form) was applied in [4]. Here there was written a system of integral equations for the hydrodynamic quantities, and its numerical solution is examined in several specific partial cases. In the present article, the problem of recondensation in a binary mixture is treated by the method of moments for Maxwellian molecules. For the case of small relative difference in the temperature of the surfaces, analytical expressions are obtained for the rates of mass transfer and the heat fluxes, making it possible to shed light on the principal special characteristics of the process of recondensation in a binary mixture.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 150–155, July–August, 1975.The authors thank R. Ya. Kucherov for his useful observations.     

A lattice Boltzmann model for diffusion of binary gas mixtures that includes diffusion slip

This paper describes the development of a lattice Boltzmann (LB) model for a binary gas mixture, and applications to channel flow driven by a density gradient with diffusion slip occurring at the wall. LB methods for single component gases typically use a non‐physical equation of state in which the relationship between pressure and density varies according to the scaling used. This is fundamentally unsuitable for extension to multi‐component systems containing gases of differing molecular masses. Substantial variations in the species densities and pressures may exist even at low Mach numbers; hence, the usual linearized equation of state for small fluctuations is unsuitable. Also, existing methods for implementing boundary conditions do not extend easily to novel boundary conditions, such as diffusion slip. The new model developed for multi‐component gases avoids the pitfalls of some other LB models. A single computational grid is shared by all the species, and the diffusivity is independent of the viscosity. The Navier–Stokes equation for the mixture and the Stefan–Maxwell diffusion equation are both recovered by the model. Diffusion slip, the non‐zero velocity of a gas mixture at a wall parallel to a concentration gradient, is successfully modelled and validated against a simple one‐dimensional model for channel flow. To increase the accuracy of the scheme, a second‐order numerical implementation is needed. This may be achieved using a variable transformation method that does not increase the computational time. Simulations were carried out on hydrogen and water diffusion through a narrow channel for varying total pressure and concentration gradients. Copyright © 2011 John Wiley & Sons, Ltd.     

Rarefied gas flow into a vacuum from a plane long channel closed at one end

The time-dependent problem of rarefied gas flow into a vacuum from a plane long channel closed at one end is solved on the basis of the kinetic S-model. The effect of diffuse molecular reflection from the channel walls on the flow velocity and the process of channel cavity vacuumization is studied as a function of the channel length and the extent of gas rarefaction under the condition that the wall temperature is maintained to be constant. The kinetic equation is solved numerically using a conservative finite-difference method of the second order of accuracy in spatial coordinates. The possibility of simplification of the problem for long times by means of reduction to the diffusion process is considered.     

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.     

Numerical investigation of the structure of nonequilibrium three-dimensional hypersonic flow past blunt bodies

Three-dimensional dissociating air flow past blunt bodies is investigated within the framework of the parabolized Navier-Stokes equations in the thin layer approximation. Multicomponent diffusion, barodiffusion and homogeneous chemical reactions, including dissociation-recombination and exchange reactions, are taken into account. The boundary conditions are assigned in the free stream and at the surface of the body with allowance for heterogeneous catalytic reactions and slip effects. The problem of flow at zero angle of attack past blunt bodies possessing two planes of symmetry is investigated numerically for flow patterns varying from smeared layer structure to almost ideal flow (Re=50-10⁵). The flow conditions corresponded to the motion of a body with lift along a re-entry trajectory [1]. The contribution of the chemical reactions in the shock wave as compared to the diffusion flux at the edge of the shock wave was estimated. The edge of the shock wave is assumed to correspond to the point at which the density profile has the greatest slope. The influence of slip effects and barodiffusion on the flow characteristics is demonstrated. The results of the calculations are compared with calculations based on the thin viscous shock layer model [2].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 143–150, September–October, 1987.The author wishes to express his thanks to G. A. Tirskii and V. V. Lunev for useful discussions and valuable advice.     

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.     

Influence of a hydrodynamic field on laminar flame stability

The hydrodynamic stability of a plane flame front has been studied by L. D. Landau in [1], where it was shown that it is absolutely unstable.The aim of this note is to clarify the influence of the hydrodynamic field curvature on flame stability. Flame stability is analyzed within the framework of Landau's theory, in which the flame front is represented by a surface on which the velocity, density, and temperature values experience discontinuities. Viscosity, diffusion, and heat conduction are neglected. The front moves at a given constant velocity relative to the gas. The gas is assumed to be incompressible in front of and behind the front.It is shown that fields exist which will both stabilize and destabilize the gas. A cylindrical flame formed by a concentrated source of given intensity is examined (two-dimensional problem). Flame stability is studied for the case of a perturbed flame front. It is shown that in this case, the hydrodynamic field has the effect of stabilizing the flame. For the first perturbation harmonics, the destabilizing effect of gas expansion appears to be relatively weak compared to the stabilizing effect of the velocity field. The first perturbation harmonics attenuate. The destabilizing effect of the velocity field is demonstrated by an example in which the radial flow of the fresh mixture is applied externally, and there exists a concentrated sink flow for the combustion products.The author wishes to express his gratitude to G. I. Barenblatt, A. G. Istratova, and V. B. Librovich for their attention to the paper and valuable advice.     

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.     

Viscous incompressible flow in a narrow channel with one free wall and fluid supply through a porous insert

The problem investigated relates the plane unsteady flow of a viscous incompressible fluid in a narrow channel one of whose walls is free and acted upon by a given load, while the other is rigidly fixed. The fluid enters the channel through a porous insert in the stationary wall. A model of the flow of a thin film of viscous incompressible fluid and Darcy's law for flow in a porous medium are used to find the distribution of fluid pressure and velocity in the channel and the porous insert in the two-dimensional formulation for fairly general boundary conditions in the case where the length of the porous insert exceeds the length of the free wall. In the particular case where the length of the porous insert is equal to the length of the free wall an exact stationary solution of the problem is obtained for a given value of the channel height. The stability of the equilibrium position of the free wall supported on a hydrodynamic fluid film is examined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 16–24, January–February, 1986.     

Rotational inviscid flow with circulation zones in a channel of variable cross section

The development of the theory of rotational motion of inviscid fluids for the purposes of describing channel flow encounters certain difficulties in connection with the appearance of viscosity effects near the walls. In the potential-rotational model [1], in which the vorticity is nonzero only in a closed circulation zone surrounded by potential flow, it is assumed that the separation and attachment points are known in advance. For example, for flow around a cavity these points coincide with the extreme corner points of the contour. The problem of determining the vorticity in a closed zone for the potential-rotational model has been investigated in a number of studies [2, 3], etc. In the case of an incompressible fluid the vorticity in the circulation zone is constant for two-dimensional flow and proportional to the distance from the axis for axisymmetric flow. The value of the constant is found from the steady-state condition for the adjoining viscous layers. If the channel walls have a smooth profile without corner points, then for determining the boundaries of the circulation zones additional conditions must be used. This study employs another scheme, in which the vorticity is formed outside the region of flow and in a particular problem is specified in the form of a boundary condition. An analytic solution describing the rotational flow of an inviscid fluid in a channel with a slightly varying cross section is obtained. Three types of entrance flow nonuniformity are considered: 1) uniform shear flow, 2) wake-type flow, and 3) potential flow with a narrow wall boundary layer. Streamline patterns with circulation zones are constructed for flows in diffuser channels with the above-mentioned types of entrance nonuniformity. A model of flow separation in a channel with a turbulent boundary layer on the walls is discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 31–37, March–April, 1985.In conclusion the author wishes to thank E. Yu. Shal'man, A. N. Kraiko, and A. B. Vatazhin for useful discussions and advice.     

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.     

Flow of Particulate-Fluid Suspension in a Channel with Porous Walls

The problem of the dispersed particulate-fluid two-phase flow in a channel with permeable walls under the effect of the Beavers and Joseph slip boundary condition is concerned in this paper. The analytical solution has been derived for the longitude pressure difference, stream functions, and the velocity distribution with the perturbation method based on a small width to length ratio of the channel. The graphical results for pressure, velocity, and stream function are presented and the effects of geometrical coefficients, the slip parameter and the volume fraction density on the pressure variation, the streamline structure and the velocity distribution are evaluated numerically and discussed. It is shown that the sinusoidal channel, accompanied by a higher friction factor, has higher pressure drop than that of the parallel-plate channel under fully developed flow conditions due to the wall-induced curvature effect. The increment of the channel’s width to the length ratio will remarkably increase the flow rate because of the enlargement of the flow area in the channel. At low Reynolds number ranging from 0 to 65, the fluids move forward smoothly following the shape of the channel. Moreover, the slip boundary condition will notably increase the fluid velocity and the decrease of the slip parameter leads to the increment of the velocity magnitude across the channel. The fluid-phase axial velocity decreases with the increment of the volume fraction density.     

A moving mesh BGK scheme for multi‐material flows

In this paper, a moving mesh BGK scheme (MMBGK) for multi‐material flow computations is proposed. The basic idea of constructing the MMBGK is to couple the Lagrangian method, which tracks material interfaces and keeps the interfaces sharp, with a remapping‐free ALE‐type kinetic method within each single material region, where the kinetic method is based on the BGK (Bhatnagar–Gross–Krook) model. Within each single material region, a numerical flux formulation is developed on moving meshes from motion of microscope particles, and the mesh velocity is determined by requiring both mesh adaptation for accuracy and robustness, such that the grids are moving towards to the regions with high flow gradients in a way of diffusive mechanism (velocity) to adjust the distances between neighboring cells, thus increasing the numerical accuracy. To keep the sharpness of material interfaces, the Lagrangian velocity and flux are constructed at the interfaces only. Consequently, a BGK‐scheme‐based ALE‐type method (i.e., the MMBGK scheme) for multi‐material flows is constructed. Numerical examples in one and two dimensions are presented to illustrate the accuracy and robustness of the MMBGK scheme. Copyright © 2011 John Wiley & Sons, Ltd.