Computation of rarefied diatomic gas flows through a plane microchannel

A numerical method based on a model kinetic equation was developed for computing diatomic rarefied gas flows in two dimensions. Nitrogen flows through a plane microchannel were computed, and the gas flow rate was constructed as a function of the Knudsen number for various channel lengths.     

Application of moment equations to the mathematical simulation of gas microflows

An approach to the simulation of moderately rarefied gas flows in a transition zone is developed. The applicability of the regularized Grad 13-moment (R13) equations to the numerical simulation of a transition flow between the continual and free-molecular gas flow regimes is explored. For the R13 equations, a numerical method is proposed that is a higher order accurate version of Godunov’s explicit method. A numerical procedure for implementing solid-wall boundary conditions is developed. One- and two-dimensional test problems are solved, including the shock wave structure and the Poiseuille flow in a plane channel. The possibility of applying the R13 equations to the simulation of plane channel and jet flows in a transition regime is explored. To this end, the flow in a square cavity generated by the motion of one of the walls is studied and the operation of the Knudsen pump is analyzed.     

Zu einigen mathematischen problemen in der theorie der transonischen strömung

The work deals with a definition of a weak solution of steady plane transonic flows past a thin profile, with the properties of the solution across a shock wave, and with a derivation of a conservative difference scheme suitable for numerical solution of the above mentioned problem by a finite difference method. The work presents several examples of numerical solution of transonic flows past a profile, through a plane cascade and some three-dimensional results. The numerical results presented are compared with experimental results or with numerical results by other authors.     

Kinetic model of the Boltzmann equation for a diatomic gas with rotational degrees of freedom

A system of model kinetic equations is proposed to describe flows of a diatomic rarefied gas (nitrogen). A conservative numerical method is developed for its solution. A shock wave structure in nitrogen is computed, and the results are compared with experimental data in a wide range of Mach numbers. The system of model kinetic equations is intended to compute complex-geometry three-dimensional flows of a diatomic gas with rotational degrees of freedom.     

Supersonic flow past a flat lattice of cylindrical rods

Two-dimensional supersonic laminar ideal gas flows past a regular flat lattice of identical circular cylinders lying in a plane perpendicular to the free-stream velocity are numerically simulated. The flows are computed by applying a multiblock numerical technique with local boundary-fitted curvilinear grids that have finite regions overlapping the global rectangular grid covering the entire computational domain. Viscous boundary layers are resolved on the local grids by applying the Navier–Stokes equations, while the aerodynamic interference of shock wave structures occurring between the lattice elements is described by the Euler equations. In the overlapping grid regions, the functions are interpolated to the grid interfaces. The regimes of supersonic lattice flow are classified. The parameter ranges in which the steady flow around the lattice is not unique are detected, and the mechanisms of hysteresis phenomena are examined.     

带激波的一维非定常两相流动的数值模拟

本文将处理带激波的单相气体非定常流动问题的两种高分辨数值方法(随机取样法和二阶GRP有限差分法)推广应用于气固悬浮体(亦称含灰气体)两相情况,计算了含灰气体激波管中两相激波特性、波后流场结构及气固两相流动参数随时间的变化.数值结果表明:这两种方法均能给出带有尖锐间断阵面的两相激波松弛结构.二阶GRP方法在计算精度和机时耗用等方面优于随机取样法.     

Nonlinear nonequilibrium kinetic model of the boltzmann equation for monatomic gases

A model kinetic equation approximating the Boltzmann equation in a wide range of nonequilibrium gas states was constructed to describe rarefied gas flows. The kinetic model was based on a distribution function depending on the absolute velocity of the gas particles. Highly efficient in numerical computations, the model kinetic equation was used to compute a shock wave structure. The numerical results were compared with experimental data for argon.     

Numerical Simulation of Detonation Excitation at Shock Wave Interaction with Dust Cloud in Channel

The study of detonation ability of reactive particle gas mixtures is necessary to prevent industrial explosions in industries where dispersed powders are used. The present paper focuses on numerical simulation of the shock wave interaction with semiinfinite aluminum dust cloud, which is situated inside a plane channel. The cloud fills entirely or partly the channel cross‐section and has initially a rectangular shape. The scenarios of detonation initiation in the cloud are determined depending on the incident shock wave amplitude values. The processes of transformation and spreading of finite width clouds under weak incident shock wave action (when the particles do not ignite) are investigated. The types of an oblique shock wave reflection from the plane of symmetry in the cloud are analyzed. The processes of particle ignition and detonation structure formation at strong incident shock wave action are investigated. Nonstationary periodic fuctuations take place in the detonation flow due to transversal wave effect. Nevertheless the detonation structure established propagates in quasistationary regime. If the incident shock wave is attenuated with a rarefaction wave then the detonation formation fails at clouds of insufficient width.     

Dusty shock wave interaction – formation of fully dispersed waves and particle focusing effect

The problem of regular (symmetric and asymmetric) interaction of plane shock waves in a steady-state dusty gas flow is considered. For near-sonic flows with a fairly high particle mass loading, the possibility of the formation of wave structures is revealed, in which either all or only some of the incident or reflected waves degenerate into so-called fully dispersed waves, i.e. zones in which no discontinuities appear in the parameters of each phase. For stronger shock waves and low particle mass concentration, the effect of aerodynamic particle focusing and the formation of a narrow high-concentration beam of particles behind the point of the interaction of the waves are detected on the basis of parametric numerical calculations. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Kinetic model of the Boltzmann equation for a power-law intermolecular interaction potential

A model kinetic equation approximating the Boltzmann equation with a linearized collision integral is constructed to describe rarefied gas flows at moderate and low Knudsen numbers. The kinetic model describes gas flows with a power-law intermolecular interaction potential and involves five relaxation parameters. The structure of a shock wave is computed, and the results are compared with an experiment for argon.     

Calculation of shock compression of porous media in conical solid-state targets with an outlet hole

The Godunov-type method that was earlier developed for the hydrodynamic equations of compressible media on regular grids is extended to the case of shock wave flows in elastoplastic porous media. The results of numerical simulations of a hypothetical experimental plant for checking the effect of an increased intensity of converging shock waves in porous media with increasing degree of porosity are presented. This effect was earlier discovered by numerical simulations.     

Solution of the Boltzmann equation for unsteady flows with shock waves in narrow channels

Unsteady rarefied gas flows in narrow channels accompanied by shock wave formation and propagation were studied by solving the Boltzmann kinetic equation. The formation of a shock wave from an initial discontinuity of gas parameters, its propagation, damping, and reflection from the channel end face were analyzed. The Boltzmann equation was solved using finite differences. The collision integral was calculated on a fixed velocity grid by a conservative projection method. A detector of shock wave position was developed to keep track of the wave front. Parallel computations were implemented on a cluster of computers with the use of the MPI technology. Plots of shock wave damping and detailed flow fields are presented.     

Multiscale approach to computation of three-dimensional gas mixture flows in engineering microchannels

A multiscale approach to computing real gas flows in engineering microchannels on high-performance computer systems in a wide range of Knudsen numbers is described. The numerical implementation of the approach combines the solution of quasigasdynamic equations and the molecular dynamics method. Following the approach, the parameters of the real gas equation of state are found at the molecular level, the kinetic gas properties are calculated, and the form of boundary conditions on the microchannel walls are determined. The technique is verified by computing several test problems. The results agree well with available theoretical and experimental data.     

Nonself-similar flow with a shock wave reflected from the center of symmetry and new self-similar solutions with two reflected shocks

In some problems concerning cylindrically and spherically symmetric unsteady ideal (inviscid and nonheat-conducting) gas flows at the axis and center of symmetry (hereafter, at the center of symmetry), the gas density vanishes and the speed of sound becomes infinite starting at some time. This situation occurs in the problem of a shock wave reflecting from the center of symmetry. For an ideal gas with constant heat capacities and their ratio γ (adiabatic exponent), the solution of this problem near the reflection point is self-similar with a self-similarity exponent determined in the course of the solution construction. Assuming that γ on the reflected shock wave decreases, if this decrease exceeds a threshold value, the flow changes substantially. Assuming that the type of the solution remains unchanged for such γ, self-similarity is preserved if a piston starts expanding from the center of symmetry at the reflection time preceded by a finite-intensity reflected shock wave propagating at the speed of sound. To answer some questions arising in this formulation, specifically, to find the solution in the absence of the piston, the evolution of a close-to-self-similar solution calculated by the method of characteristics is traced. The required modification of the method of characteristics and the results obtained with it are described. The numerical results reveal a number of unexpected features. As a result, new self-similar solutions are constructed in which two (rather than one) shock waves reflect from the center of symmetry in the absence of the piston.     

Unsteady rarefied gas flow in a microchannel driven by a pressure difference

The kinetic S-model is used to study the unsteady rarefied gas flow through a plane channel between two parallel infinite plates. Initially, the gas is at rest and is separated by the plane x = 0 with different pressure values on opposite sides. The gas deceleration effect of the channel walls is studied depending on the degree of gas rarefaction and the initial pressure drop, assuming that the molecules are diffusely reflected from the boundary. The decay of the shock wave and the disappearance of the uniform flow region behind the shock wave are monitored. Special attention is given to the gas mass flux through the cross section at x = 0, which is computed as a function of time. The asymptotic behavior of the solution at unboundedly increasing time is analyzed. The kinetic equation is solved numerically by applying a conservative finite-difference method of second-order accuracy in space.     

Spherical shock wave generated by a moving piston in mixture of a non-ideal gas and small solid particles under a gravitational field

Self-similar solutions are obtained for one-dimensional isothermal and adiabatic unsteady flows behind a strong spherical shock wave propagating in a dusty gas. The shock is assumed to be driven out by a moving piston and the dusty gas to be a mixture of a non-ideal (or perfect) gas and small solid particles, in which solid particles are continuously distributed. It is assumed that the equilibrium flow-conditions are maintained and variable energy input is continuously supplied by the piston. The medium is under the influence of the gravitational field due to a heavy nucleus at the origin (Roche model). The effects of an increase in the mass concentration of solid particles, the ratio of the density of the solid particles to the initial density of the gas, the gravitational parameter and the parameter of non-idealness of the gas in the mixture, are investigated. It is shown that due to presence of gravitational field the compressibility of the medium at any point in the flow-field behind the shock decreases and all other flow-variables and the shock strength increase. A comparison has also been made between the isothermal and adiabatic flows. It is investigated that the singularity in the density and compressibility distributions near the piston in the case of adiabatic flow are removed when the flow is isothermal.     

On the acceleration of shock waves and concentration of energy

By a series of simple examples related to exact solutions of problems in gas dynamics and magnetohydrodynamics, possible mechanisms of acceleration of shock waves and concentration of energy are elucidated. The acceleration of a shock wave is investigated in the problem of motion of a plane piston at a constant velocity in the case when the initial density of the medium drops in the presence of constant counterpressure. It is shown that in this situation a “blow-up” regime is induced by a shock wave going to infinity in finite time even for limited work of the piston. A simple spherically symmetric solution with a converging shock wave is constructed and shown to lead to the concentration of energy. A general method for solving one-dimensional non-self-similar problems related to matching the equilibrium state to a motion with homogeneous deformation on a shock wave is discussed; this method leads to a solution in quadratures.     

Numerical verification of delta shock waves for pressureless gas dynamics

The subject of this paper is theoretical analysis and numerical verification of delta shock wave existence for pressureless gas dynamic system. The existence of overcompressive delta shock wave solution in the framework of Colombeau generalized functions is proved. This result is verified numerically by specially designed procedure that is based on wave propagation method implemented in CLAWPACK. The method is coupled with dynamic refinement mesh. We also consider a strictly hyperbolic system obtained from the original one by perturbation and change of variables. The same numerical procedure is applied to the perturbed problem. The obtained numerical results in both cases confirm theoretical expectations.     

On behavior of two-dimensional non-plane parallel turbulent flows of inviscid gas

The time behavior of two-dimensional flows of inviscid gas in which the velocity component normal to the plane of independent variables and the vorticity components parallel to this plane are different from zero, is investigated. Equations of such flows form two different subsystems. The first subsystem describes a plane parallel (“primary”) flow without the third velocity component, and is independent of the second subsystem consisting of a single equation for the third velocity component and determining the “secondary” flow. The flows are analyzed with sufficient detail without using numerical integration which carries with it unavoidable errors, and without linearization, both of which are employed to a lesser or greater degree in the study of the evolution of vortical structures (see /1–6/).     

Stability of a supersonic flow past a wedge with adjoint weak neutrally stable shock wave

We study the classical problem of a supersonic stationary flow of a nonviscous nonheat-conducting gas in local thermodynamic equilibrium past an infinite plane wedge. Under the Lopatinski? condition on the shock wave (neutral stability), we prove the well-posedness of the linearized mixed problem (the main solution is a weak shock wave), obtain a representation of the classical solution, where, in this case (in contrast to the case of the uniform Lopatinski? condition—an absolutely stable shock wave), plane waves additionally appear in the representation. If the initial data have compact support, the solution reaches the given regime in infinite time.