Implicit numerical method for computing three-dimensional rarefied gas flows on unstructured meshes

The method based on the numerical solution of a model kinetic equation is proposed for analyzing three-dimensional rarefied gas flows. The basic idea behind the method is the use of a second-order accurate TVD scheme on hybrid unstructured meshes in physical space and a fast implicit time discretization method without iterations at the upper level. The performance of the method is illustrated by computing test examples of three-dimensional rarefied gas flows in variously shaped channels in a wide range of Knudsen numbers.     

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.     

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.     

Simulation of rarefied gas flows on the basis of the Boltzmann kinetic equation solved by applying a conservative projection method

Flows of a simple rarefied gas and gas mixtures are computed on the basis of the Boltzmann kinetic equation, which is solved by applying various versions of the conservative projection method, namely, a two-point method for a simple gas and gas mixtures with a small difference between the molecular masses and a multipoint method in the case of a large mass difference. Examples of steady and unsteady flows are computed in a wide range of Mach and Knudsen numbers.     

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

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

Mathematical modeling of shock-wave processes under gas solid boundary interaction

The results of numerical simulations are presented for planar air flows in a bounded volume of square cross section diminishing due to a uniform motion of the walls, for a flow of a propane-air mixture under sinusoidal variation of the size of the square domain, and for three-dimensional supersonic air and propane-air flows in channels of variable square cross section. Specific features of shock-wave processes that are associated with the piston effect and cumulation are established. The hypersonic analogy between planar and spatial flows is confirmed, which allows one to use two-dimensional solutions in estimating three-dimensional flows. The equations of a multicomponent ideal perfect gas and one-stage kinetics of chemical reactions are used to describe the flows. The method of numerical simulations is based on S.K. Godunov’s scheme and implemented within an original software package.     

Numerical study of unsteady rarefied diatomic gas flows in a plane microchannel

The numerical solution of a kinetic equation for a diatomic gas (nitrogen) is used to study two-dimensional unsteady gas flows in a plane microchannel caused by discontinuous in the initial distributions of macroscopic gas parameters. The plane discontinuity fronts are perpendicular to the walls of the channel. The arising flows are model ones for gas flows in a shock tube and a microchannel. The reflection of an incident shock wave from a flat end face is studied. It is found that the gas piles up at the cold wall, which slows down the shock wave detachment. The numerical results are in qualitative agreement with experimental data.     

Numerical method for computing two-dimensional unsteady rarefied gas flows in arbitrarily shaped domains

A high-order accurate method for analyzing two-dimensional rarefied gas flows is proposed on the basis of a nonstationary kinetic equation in arbitrarily shaped regions. The basic idea behind the method is the use of hybrid unstructured meshes in physical space. Special attention is given to the performance of the method in a wide range of Knudsen numbers and to accurate approximations of boundary conditions. Examples calculations are provided.     

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.     

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.     

Modeling of flow separation of assist gas as applied to laser cutting of thick sheet metal

The present paper describes the results of mathematical modeling of supersonic flows of a viscous compressible gas, obtained by numerically solving three-dimensional full Navier–Stokes equations, and also the results of experiments with visualization of gas jet flows in channels geometrically similar to the laser cut. Separation of the gas flow from the cut front is predicted numerically and then validated by experiments on a model setup. The gas flow structure arising in a narrow channel behind a sonic (conical) or supersonic nozzle is described. Specific features of originating in the flow separation on a smooth surface in a narrow channel are examined, and mechanisms controlling the separation are proposed. Flow separation directly affects the changes in the shape and structure of striations and is the one of main reason for the worse quality of the laser cut surface. It is shown that the changes in the structures of striations over the thickness of the sheet being cut are closely related to aerodynamic features of jet flows of the assisting gas in the cut channel.     

Unidirectional flows of viscoelastic fluids and their gasdynamic counterparts

We have shown that unidirectional flows of viscoelastic fluids are mathematically equivalent to plane potential flows of a fictitious gas. The classification of type and some exact solutions of the equations describing fictitious gas flows are presented. The Prandtl-Mayer flow and its viscoelastic counterpart is discussed in more details in connection with the Sternberg-Koiter paradox.     

The non-stationary invariant solution of the equations of gas dynamics describing the spreading of a gas into a vacuum

An invariant solution of the equations of gas dynamics, constructed on a one-dimensional subgroup (according to the classification in /1/) which is only allowed in the case of a polytropic gas with a special adiabatic index, is considered. The gas spreads out into a vacuum after a finite time. New solutions are constructed which describe one-dimensional flows from a source into a vacuum and the focussing of the gas within a sphere or a cylinder with shock waves. The spreading of a concentration of the gas with an arbitrary boundary when there is a contact discontinuity is also considered.

One-dimensional flows have been treated in detail in /2, 3/, mainly in the case of extended subgroups.     

Steady motions of a continuous medium, resonances and Lagrangian turbulence

A method which enables one to establish a non-regularity property of the motion of fluid particles (known as chaotic advection or Lagrangian turbulence) for typical steady flows is developed. The method is based on expanding solutions of the equations of motion of a continuous medium in powers of a small parameter and using the conditions for the destruction of invariant resonant tori when perturbations are added. It is shown that the velocity field, defined as the solution of the Burgers equations, generates a generally non-regular dynamical system. For an ideal barotropic fluid in an irrotational force field, the method proposed yields a well-known necessary condition for chaotization: the velocity field is collinear with its curl. Special attention is given to investigating the chaotization of typical steady flows of a heat-conducting perfect gas.     

Three-dimensional generalization for <Emphasis Type="Italic">W</Emphasis> modification of a Godunov method

A high-accuracy modification of Godunov’s method for three-dimensional unsteady ideal gas flows is proposed. For the linear advection equation, a fully three-dimensional second-order accurate monotone scheme is designed with corrections computed on a variable stencil whose orientation depends on the signs of the equation coefficients. For the linear scalar advection equation, the scheme is proved to possess the positive approximation property. The method is tested by computing the flow in a three-dimensional Ludwieg tube with a square cross section.     

动力学的流矢量分裂法模拟浅水波方程

Based on the analogy to gas dynamics, the kinetic flux vector splitting (KFVS) method is used to stimulate the shallow water wave equations. The flux vectors of the equations are split on the basis of the local equilibrium Maxwell-Boltzmann distribution. One dimensional examples including a dam breaking wave and flows over a ridge are calculated. The solutions exhibit second-order accuracy with no spurious oscillation.     

A generalized Riemann problem for quasi-one-dimensional gas flows

A generalization of the Riemann problem for gas dynamical flows influenced by curved geometry, such as flows in a variable-area duct, is solved. For this generalized Riemann problem the initial data consist of a pair of steady-state solutions separated by a jump discontinuity. The solution of the generalized Riemann problem is used as a basis for a random choice method in which steady-state solutions are used as an Ansatz to approximate the spatial variation of the solution between grid points. For nearly steady flow in a Laval nozzle, where this Ansatz is appropriate, this generalized random choice method gives greatly improved results.     

Implicit finite element schemes for stationary compressible particle-laden gas flows

The derivation of macroscopic models for particle-laden gas flows is reviewed. Semi-implicit and Newton-like finite element methods are developed for the stationary two-fluid model governing compressible particle-laden gas flows. The Galerkin discretization of the inviscid fluxes is potentially oscillatory and unstable. To suppress numerical oscillations, the spatial discretization is performed by a high-resolution finite element scheme based on algebraic flux correction. A multidimensional limiter of TVD type is employed. An important goal is the efficient computation of stationary solutions in a wide range of Mach numbers. This is a challenging task due to oscillatory correction factors associated with TVD-type flux limiters and the additional strong nonlinearity caused by interfacial coupling terms. A semi-implicit scheme is derived by a time-lagged linearization of the nonlinear residual, and a Newton-like method is obtained in the limit of infinite CFL numbers. The original Jacobian is replaced by a low-order approximation. Special emphasis is laid on the numerical treatment of weakly imposed boundary conditions. It is shown that the proposed approach offers unconditional stability and faster convergence rates for increasing CFL numbers. The strongly coupled solver is compared to operator splitting techniques, which are shown to be less robust.