Existence of solutions of Kraiko's problem

Three initial-boundary-value problems for the equations of gas dynamics are formulated. Successive solution of these problems yields a solution of Kraiko's problem of the isentropic transition of an ideal gas from a homogeneous state of rest to another state of rest with higher or lower density. Solutions are constructed for plane, cylindrical, and spherical layers of an ideal gas. The existence of locally analytic solutions is proved. Ural State Academy of Service, Ekaterinburg 620034. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 3, pp. 48–55, May–June, 2000.     

Explicit Rate–Time Solutions for Modeling Matrix–Fracture Flow of Single Phase Gas in Dual-Porosity Media

One of the widely used methods for modeling matrix–fracture fluid exchange in naturally fractured reservoirs is dual porosity approach. In this type of modeling, matrix blocks are regarded as sources/sinks in the fracture network medium. The rate of fluid transfer from matrix blocks into fracture medium may be modeled using shape factor concept (Warren and Root, SPEJ 3:245–255, 1963); or the rate–time solution is directly derived for the specific matrix geometry (de Swaan, SPEJ 16:117–122, 1976). Numerous works have been conducted to study matrix–fracture fluid exchange for slightly compressible fluids (e.g. oil). However, little attention has been taken to systems containing gas (compressible fluid). The objective of this work is to develop explicit rate–time solutions for matrix–fracture fluid transfer in systems containing single phase gas. For this purpose, the governing equation describing flow of gas from matrix block into fracture system is linearized using pseudopressure and pseudotime functions. Then, the governing equation is solved under specific boundary conditions to obtain an implicit relation between rate and time. Since rate calculations using such an implicit relation need iterations, which may be computationally inconvenient, an explicit rate–time relation is developed with the aid of material balance equation and several specific assumptions. Also, expressions are derived for average pseudopressure in matrix block. Furthermore, simplified solutions (originated from the complex general solutions) are introduced applicable in infinite and finite acting flow periods in matrix. Based on the derived solutions, expressions are developed for shape factor. An important observation is that the shape factor for gas systems is the same as that of oil bearing matrix blocks. Subsequently, a multiplier is introduced which relates rate to matrix pressure instead of matrix pseudopressure. Finally, the introduced equations are verified using a numerical simulator.     

“Simple” solutions of the equations of dynamics for a polytropic gas

The notion of a “simple” solution of a system of differential equations that admit a local Lie group G of transformations of the basic space is considered as an invariant H-solution of type (0, 0) with respect to the subgroup HυG. Such solutions are attractive since they are described by explicit formulas that provide a clear physical interpretation for them. For gas-dynamic equations with a polytropic gas law, all simple solutions that are not related to special forms of gas flow are listed. Examples of simple solutions are given and the collapse phenomenon, which has been previously studied for barochronic flows, is described. Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 2, pp. 5–12, March–April, 1999.     

Two-dimensional gas vortices and twisted gas jets

This paper studies an invariant solution of rank one of the equations of motion of a polytropic gas that describes two-dimensional gas vortices and twisted gas jets. Flow types are classified according to the governing parameter: vortices in the form of sources and sinks, unlimited expansion, and collapse. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 71–81, May–June, 2009.     

Direct numerical simulation of a liquid sheet in a compressible gas stream in axisymmetric and planar configurations

A thin liquid sheet present in the shear layer of a compressible gas jet is investigated using an Eulerian approach with mixed-fluid treatment for the governing equations describing the gas–liquid two-phase flow system, where the gas is treated as fully compressible and the liquid as incompressible. The effects of different topological configurations, surface tension, gas pressure and liquid sheet thickness on the flow development of the gas–liquid two-phase flow system have been examined by direct solution of the compressible Navier–Stokes equations using highly accurate numerical schemes. The interface dynamics are captured using volume of fluid and continuum surface force models. The simulations show that the dispersion of the liquid sheet is dominated by vortical structures formed at the jet shear layer due to the Kelvin–Helmholtz instability. The axisymmetric case is less vortical than its planar counterpart that exhibits formation of larger vortical structures and larger liquid dispersion. It has been identified that the vorticity development and the liquid dispersion in a planar configuration are increased at the absence of surface tension, which when present, tends to oppose the development of the Kelvin–Helmholtz instability. An opposite trend was observed for an axisymmetric configuration where surface tension tends to promote the development of vorticity. An increase in vorticity development and liquid dispersion was observed for increased liquid sheet thickness, while a decreasing trend was observed for higher gas pressure. Therefore surface tension, liquid sheet thickness and gas pressure factors all affect the flow vorticity which consequently affects the dispersion of the liquid.      

An asymptotically self-similar solution of the problem of three-dimensional gas flow through a porous medium

An exact solution of the Cauchy problem is constructed for the equation describing the three-dimensional molecular diffusion of a gas. The result obtained is a natural generalization of the solution of the analogous Boussinesq problem. Dnepropetrovsk. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 178–180, September–October, 1988.     

Negative Saturation Approach for Non-Isothermal Compositional Two-Phase Flow Simulations

This article deals with developing a solution approach, called the non-isothermal negative saturation (NegSat) solution approach. The NegSat solution approach solves efficiently any non-isothermal compositional flow problem that involves phase disappearance, phase appearance, and phase transition. The advantage of the solution approach is that it circumvents using different equations for single-phase and two-phase regions and the ensuing unstable procedure. This paper shows that the NegSat solution approach can also be used for non-isothermal systems. The NegSat solution approach can be implemented efficiently in numerical simulators to tackle modeling issues for mixed CO₂–water injection in geothermal reservoirs, thermal recovery processes, and for multicontact miscible and immiscible gas injection in oil reservoirs. We illustrate the approach by way of example to cold mixed CO₂–water injection in a 1D geothermal reservoir. The solution is compared with an analytical solution obtained with the wave-curve method (method of characteristics) and shows excellent agreement. A complete set of simulations is carried out, which identifies six bifurcations. The two main bifurcations are (1) when the most downstream compositional wave is replaced by a compositional shock and (2) when an extra Buckley–Leverett rarefaction appears. The plot of the useful energy (exergy) versus the CO₂ storage capacity shows a Z-shape. The top horizontal part represents a branch of high exergy recovery/relatively lower storage capacity, whereas the bottom part represents a branch of lower exergy recovery/higher storage capacity.     

Non-linear Evolution of P-waves in Viscous–Elastic Granular Saturated Media

The Frenkel–Biot P-wave of the first type is a seismic longitudinal wave observed in rocks fully saturated with oil, water or high-pressure gas. The P-wave of the second type is observed in unsaturated soils and other porous media saturated with gas of low pressure. Their models include properties of the skeleton, that is, its elastic modules and its own viscosity. If the non-linear terms are accounted for, the asymptotic analysis, usual for weak non-linear waves, might be applied to get the wave spectrum evolution. The wetness of grains contacts in soils and such components of oil as tars or bitumen, which attached to the skeleton, can be described by generalized viscous–elastic stress–strain connections. The latter are nominated in such a way that creates the narrow frequency interval of wave of negative dissipation where the non-linear terms begin to play the main role besides the neutral stability for waves of zero wave number. The corresponding case, relevant to single continuum model, was analyzed in the literature. Here it is shown that the interpenetrating continua with interaction of the Darcy type provide the dissipation sink in the wave evolution equation. This generalization, (Tribelsky, M.I.: Phys. Rev. Lett. (2007, submitted)), can stabilize the asymptotic solution of the evolution equation, where the dispersion terms are omitted. The asymptotic solution of the equation is invariant to initial conditions and it means a transformation of initial wave spectra to unique one while wave is spreading in the viscous–elastic medium under consideration. This explains the phenomenon, observed in wave tests at marine beach, when any dynamics action (impact, explosion, and ultrasound action) created at some distance a wave of a single frequency (~25 Hz).     

Diffusive dispersion of natural gas hydrates in ocean sediments

A mathematical model of the process of dissolving and dispersion of natural methane hydrates in ocean sediments is proposed. A condition which must be satisfied on the solution boundary of the hydrates in the region of their stable thermodynamic state is derived. The calculated solution rates are several centimeters per year. It is shown that the temperature perturbations initiated by the transition of a gas from the hydrate into the dissolved state are small and the process can be considered isothermal. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 170–173, January–February, 1999. The work was carried out with financial support from the Russian Foundation for Fundamental Research (project No. 96-01-00521).     

Heat transfer to blunt-edged swept wings at angle of attack in hypersonic flow

An approximate analytic solution is obtained for the relative heat fluxes on the lateral surface of swept wings of infinite span in a hypersonic viscous gas flow at angles of attack and yaw. The accuracy of the relations obtained is estimated on the basis of a comparison with numerical solutions. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 170–179, July–August, 1994.     

Navier�CStokes Equations with Nonhomogeneous Boundary Conditions in a Bounded Three-Dimensional Domain

Following the ideas developed in Girinon (Annales de l’Institut Poincaré. Analyse Non Linéaire 26:2025–2053, 2009), we prove the existence of a weak solution to Navier–Stokes equations describing the isentropic flow of a gas in a bounded region, W ì R³{\Omega\subset \mathbf{R}^{3}} , with nonhomogeneous Dirichlet boundary conditions on ∂Ω.     

Computational Modeling of the Structure of a High-Current Discharge in a Magnetogasdynamic Channel

A two-dimensional computational model is proposed to calculate radiative-convective heat transfer in gas flows with large gradients of physical properties. The model is based on the numerical solution of the unsteady dynamic equations for a compressible inviscid gas and the radiative transfer equations. Flow calculations for the magnetogasdynamic channel of a rail accelerator show that the dynamics of the process is substantially affected by the flow in the discharge region and hydrodynamic instability, resulting in the nonstationarity and nonuniformity of the flow and discharge structure. During the process, the discharge can exist both in the form of several current-carrying channels and in the form of a unified plasma formation. Results of the numerical calculations agree qualitative with experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 5–13, November–December, 2005.     

Optimum clearance of a gas hydrostatic thrust bearing with maximum load capacity

The optimum design of a gas hydrostatic thrust bearing clearance is obtained using the methods of calculus of variations. The variational problem of determining the clearance shape giving the maximum load capacity is solved for a given external pressurization and various journal speeds. The structure of the optimum solution is found on the basis of the gas lubrication approximation with and without constraints on the height of the bearing pad (pocket). The calculation results embrace all possible values of the parameters. A comparison with optimum liquid bearings is carried out. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 68–78, July–August, 2000. The work was carried out with support from the Russian Foundation for Basic Research (projects Nos. 99-01-01211 and 96-15-96158).     

On a gas source in a constant force field

The nonbarochronic regular partially invariant submodel of the equations of gas dynamics is studied. The submodel reduces to an implicit ordinary differential equation of the first order for an auxiliary function X = X(x). The physical quantities (velocity, density, and pressure) are expressed in terms of the function X. The properties of the solutions of the equation are investigated and interpreted physically in terms of gas motion. The existence of a shock-wave solution is proved. The properties of the shock adiabat are studied. It is shown that the results obtained are new and differ significantly from the results for the case of no constant force. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 3–16, November–December, 2006.     

Triple configurations of traveling shock waves in inviscid gas flows

Triple configurations of shock waves in supersonic inviscid flows of a perfect gas are considered. The basic parameters of triple configurations are determined, and the conditions of solution existence are analyzed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 3–10, May–June, 2008.     

Subsonic rarefied gas flow over a rack of flat transverse plates

Parameters of a rarefied gas flow through a rack of flat plates aligned across the flow are studied by means of the joint numerical solution of the Boltzmann and Navier-Stokes equations. A subsonic flow regime is considered. The changes in flow characteristics are calculated as functions of the free-stream velocity and plate temperature. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 1, pp. 59–67, January–February, 2008.     

Mathematical modeling of the motion of a portion of helium under pulsed injection over a fixed bed of cenospheres

A mathematical model is constructed and an analytical solution is obtained for the problem of a one-dimensional steady flow of a mixture of different gases with hollow permeable particles. The case of a one-dimensional unsteady flow of such a mixture is analyzed numerically. The numerical solutions are compared with experimental data on the motion of the peak concentration of helium in a fixed bed filled with cenospheres (solid hollow permeable spherical particles). The permeability of cenosphere walls and the drag coeficient of cenospheres in the gas flow are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 92–102, May–June, 2007.     

Near-resonance highly nonlinear gas oscillations in tubes

Near-resonance highly nonlinear ideal perfect gas oscillations in tubes are studied numerically for boundary conditions of various types. The oscillations are initiated by weak periodic perturbations at one end of the tube. As distinct from earlier studies [1–10], the oscillation amplitudes were not assumed to be small and the entropy increase at the shock waves formed was taken into account. Periodic flow regimes result as a limit of the solution of a Cauchy problem for one-dimensional time-dependent gasdynamic equations. The frequency responses of the oscillations under consideration are determined for boundary conditions of various types. It is shown that in specific cases the attainment of a periodic regime is accompanied by the appearance of long-wave modulations. The “repeated resonance” effect is revealed. This is due to the change in the tube's natural acoustic frequency, which takes place during the heating of the gas in the tube by the shock waves traveling in it. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 150–157, July–August, 1994.     

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).     

A phenomenological and extended continuum approach for modelling non-equilibrium flows

This paper presents a new technique that combines Grad’s 13-moment equations (G13) with a phenomenological approach to rarefied gas flows. This combination and the proposed solution technique capture some important non-equilibrium phenomena that appear in the early continuum-transition flow regime. In contrast to the fully coupled 13-moment equation set, a significant advantage of the present solution technique is that it does not require extra boundary conditions explicitly; Grad’s equations for viscous stress and heat flux are used as constitutive relations for the conservation equations instead of being solved as equations of transport. The relative computational cost of this novel technique is low in comparison to other methods, such as fully coupled solutions involving many moments or discrete methods. In this study, the proposed numerical procedure is tested on a planar Couette flow case, and the results are compared to predictions obtained from the direct simulation Monte Carlo method. This test case highlights the presence of normal viscous stresses and tangential heat fluxes that arise from non-equilibrium phenomena, which cannot be captured by the Navier–Stokes–Fourier constitutive equations or phenomenological modifications.