Similarity solution for a cylindrical shock wave in a rotational axisymmetric dusty gas with heat conduction and radiation heat flux

The propagation of shock waves in a rotational axisymmetric dusty gas with heat conduction and radiation heat flux, which has a variable azimuthally fluid velocity together with a variable axial fluid velocity, is investigated. The dusty gas is assumed to be a mixture of non-ideal (or perfect) gas and small solid particles, in which solid particles are continuously distributed. It is assumed that the equilibrium flow-condition is maintained and variable energy input is continuously supplied by the piston (or inner expanding surface). The fluid velocities in the ambient medium are assume to be vary and obey power laws. The density of the ambient medium is assumed to be constant, the heat conduction is express in terms of Fourier’s law and the radiation is considered to be of the diffusion type for an optically thick grey gas model. The thermal conductivity K and the absorption coefficient α_R are assumed to vary with temperature and density. In order to obtain the similarity solutions the angular velocity of the ambient medium is assume to be decreasing as the distance from the axis increases. The effects of the variation of the heat transfer parameter and non-idealness of the gas in the mixture are investigated. The effects of an increase in (i) the mass concentration of solid particles in the mixture and (ii) the ratio of the density of solid particles to the initial density of the gas on the flow variables are also investigated.     

Similarity solution for a spherical shock wave in a non‐ideal gas under the influence of gravitational field and monochromatic radiation with increasing energy

The propagation of a spherical shock wave in a non‐ideal gas with or without gravitational effects is investigated under the action of monochromatic radiation. Similarity solutions are obtained for adiabatic flow between the shock and the piston. The numerical solutions are obtained using the Runge‐Kutta method of the fourth order. The density of the gas is assumed to be constant. The total energy of the shock wave is non‐constant and varies with time. The effects of change in values of non‐idealness parameter, gravitational parameter, shock Mach number, radiation parameter, and adiabatic exponent of the gas on shock strength and flow variables are worked out in detail. It is investigated that the presence of gravitational field increases the compressibility of the medium, due to which it is compressed and, therefore, the distance between the inner contact surface and the shock surface is reduced. A comparison is also made between the solutions in the cases of the gravitating and the non‐gravitating media. It is manifested that the gravitational parameter and the radiation parameter have in general opposite behaviour on the flow variables and the shock strength.     

Wave front analysis of weak nonlinear waves in a two-phase flow of gas and dust particles

The non-linear behavior of waves including the characteristic front, the expansion wave front and the shock front, in a mixture of gas and dust particles has been studied. Such waves are conceived of as produced by a piston moving with a small velocity as compared with the speed of sound. The trajectories of these waves and the particle paths in the physical plane are determined. The effect of solid particles and the adiabatic heat exponent on the wave propagation is also investigated.     

Analytic solution of self-similar strong shocks in an exponential medium for zero temperature gradient

The self-similar one-dimensional propagation of a strong shock wave in a medium with an exponentially decreasing density is studied. The flow behind the shock is assumed to be spatially isothermal rather than adiabatic to simulate the conditions of large radiative transfer behind the shock. The solution in closed form is obtained. An analytic expression for the similarity exponent has also been obtained.     

Approximate analytical solutions for strong shocks with variable energy

Approximate analytical solutions are obtained for self-similar flows behind strong shocks with variable energy deposition or withdrawal at the wavefront in a perfect gas at rest with constant initial density. Numerical solutions are also obtained and the approximate solutions agree with these solutions. The effect of the adiabatic index on the solutions is investigated. The dependence of shock density ratio on the parameter characterizing the energy of the flow is studied. It is observed that the rate of deposition of energy at the wavefront decreases with increase of the parameter that specifies the total energy of the flow.     

Self-similar time-varying flows of an ideal gas with a change in the adiabatic exponent in a “reflected” shock wave

Self-similar one-dimensional time-varying problems are considered under the assumption that there is a change in the adiabatic exponent in a shock wave (SW) running (“reflected”) from a centre or axis of symmetry (later from a centre of symmetry, CS) or from a plane. The medium is an ideal (inviscid and non-heat-conducting) perfect gas with constant heat capacities. In problems with strong SW, the change in the adiabatic exponent in a gas approximately simulates physicochemical processes such as dissociation and ionization and, in the problem of the collapse of a spherical cavity in a liquid, the conversion of liquid into vapour. In both cases, the adiabatic exponent decreases on passing across a reflected SW. Problems of the collapse of a spherical cavity, the reflection of a strong SW from a centre of symmetry and a simpler problem with a self-similarity index of one are examined. When it is assumed that there is an increase in the adiabatic exponent, the self-similar solutions of the first two problems are rejected due to the decrease in entropy from the instant when the SW is reflected. When it is assumed that there is a decrease in the adiabatic exponent, the solutions of these problems only become unsuitable after a finite time has elapsed for the same reason. Up to this time when the decrease in the adiabatic exponent has not reached a certain threshold, the structure of the self-similar solution does not undergo qualitative changes. When the above-mentioned threshold is exceeded, a self-similar solution is possible if a cylindrical or spherical piston expands according to a special law from the instant of SW reflection from the CS. When there is no piston, the flow behind the reflected wave becomes non-self- similar. In the case of the deceleration of a plane flow, conditions are possible with the joining of SW from different sides to a centred rarefaction wave.     

Converging strong shock waves in magnetogasdynamics under isothermal condition

The similarity solutions to the problem of cylindrically symmetric strong shock waves in an ideal gas with a constant azimuthal magnetic field are presented. The flow behind the shock wave is assumed to spatially isothermal rather than adiabatic. We use the method of Lie group invariance to determine the possible class of self-similar solutions. Infinitesimal generators of Lie group transformations are determined by using the invariance surface conditions to the system and on the basis of arbitrary constants occurring in the expressions for the generators, four different possible cases of the solutions are reckoned and we observed that only two out of all possibilities hold self-similar solutions, one of which follows the power law and another follows the exponential law. To obtain the similarity exponents numerical calculations have been performed and comparison is made with the existing results in the literature. The flow patterns behind the shock are analyzed graphically.

     

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.     

Hypersonic flow past a sphere

The effects of dissociation or ionization of air on the analytical solution for hypersonic flow past a sphere are considered here, under certain assumptions. It has been assumed that the shock wave is in the shape of a sphere, that the density ratio across the shock is constant, that the flow behind the shock is at constant density and that dissociation or ionization only occurs behind the shock wave. Thus the effects of the compressibility of the air, variation of density ratio along the shock, and the department of the shock shape from being circular are not taken into account. Here the velocity, pressure, temperature, pressure coefficient and vorticity, etc., at any point between the shock and the surface of the sphere in the presence of dissociation or ionization are obtained. In addition, shock detachment distance, drag coefficient, stagnation point velocity gradient and sonic points on the shock and the surface have also been obtained. The results have been compared with the corresponding results obtained in the case when dissociation or ionization does not occur behind the shock.     

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

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

Self-similar convergence of a shock wave in a heat conducting gas

The problem of the convergence of a spherical shock wave (SW) to the centre, taking into account the thermal conductivity of the gas in front of the SW, is considered within the limits of a proposed approximate model of a heat conducting gas with an infinitely high thermal conductivity and a small temperature gradient, such that the heat flux is finite in a small region in front of the converging SW. In this model, there is a phase transition in the surface of the SW from a perfect gas to another gas with different constant specific heat and the heat outflow. The gas is polytropic and perfect behind the SW. Constraints are derived which are imposed on the self-similarity indices as a function of the adiabatic exponents on the two sides of the SW. In front of the SW, the temperature and density increase without limit. In the general case, a set of self-similar solutions with two self-similarity indices exists but, in the case of strong SW close to the limiting compression, there are two solutions, each of which is completely determined by the motion of the spherical piston causing the self-similar convergence of the SW.     

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)     

Approximate calculation method for one-dimensional gas flows induced by nonmonotonic motion of a piston : PMM vol. 40, n 6, 1976, pp. 1058–1064

The problem of one-dimensional piston which at the beginning moves with increasing velocity into a gas at rest, then is decelerated, and finally stops, is solved by means of special series. The gas flow field is constructed by a successive joining of three characteristic Cauchy problems in terms of their characteristic solutions. Generalized solution of the problem of instantaneous arrest of the piston is derived. Obtained equations are used for the approximate calculation of the motion of generated shock waves.Representation of solutions of certain boundary value problems for nonlinear equations of the hyperbolic kind in the form of special series was proposed in [1, 2], The problem of the piston moving into a gas at rest is solved there, and the obtained solution was used for an approximate determination of the generated shock wave. The piston velocity was assumed to be monotonically increasing. That problem is solved here with the use of similar series in the case when the piston velocity is nonmonotonous,Numerical methods make it possible at present to determine one-dimensional flows similar to that considered below, and multidimensional problems can be solved by the method proposed in [1, 2]. The use of the proposed scheme for solving the problem of the multidimensional piston, whose velocity is nonmonotonous, does not present theoretical difficulties, but except that the formulas are more cumbersome.     

Shock wave structure in an inviscid heat-conducting locally equilibrium medium of a perfect gas with radiation

The structure of a shock wave (SW) in an inviscid heat-conducting locally equilibrium medium of a perfect gas with radiation at temperatures up to hundreds of thousands and millions of degrees is analysed in the approximation of non-linear heat conduction. At such temperatures, heat conduction is far more important than viscosity, the medium consists of ions, electrons and radiation, and the electrons and radiation make a decisive contribution to the heat transfer and an appreciable contribution to the thermodynamic functions of the medium. Combinations with the dimensions of all the flow parameters are constructed from the gas density ahead of the SW and the constants appearing in the equations of state of such a medium. If they are taken as the corresponding scales, the equations of state for different media in dimensionless form only differ in the ratio of the specific heat capacities (the adiabatic exponent) of the gas. The SW structure is studied for adiabatic exponents from 1 to 3, dimensionless temperatures ahead of the SW from zero to infinity and any SW velocities exceeding the speed of sound in the unperturbed medium. The following cases are successively considered: 1) when the contribution due to radiation is neglected in the equations of state, 2) when radiation is taken into account and the SW moves through a gas at a zero temperature (through a cold background), and 3) the radiation contribution is allowed for in the case of a warm background. It is ascertained when the SW structure is continuous and when it contains a finite or infinite forerunner and an isothermal shock. The transition from a continuous structure to a structure with an isothermal shock and the intensity of this shock are independent of the form of the formulae for the thermal conductivity in Fourier's law. In the case of adiabatic exponents greater than unity, the structure is continuous for all SW velocities starting from a certain dimensionless background temperature.     

Self-similar strong shocks in non-uniform atmosphere

The propagation of strong shocks in an atmosphere of variable density at rest is studied. The energy gain of the flow enveloped by the shock is assumed to be time-dependent. Analytical and numerical solutions of the similarity flows behind such shocks are obtained.     

含灰气体激波沿平壁传播时诱导的边界层流动*

本文研究合灰气体激波沿平直壁面传播过程中在壁面附近形成的层流边界层流动。我们依照双连续介质双向耦合模型处理含灰气体激波的波后流动及其诱导的边界层问题,控制方程采用有限差分方法数值求解,给出了激波下游两相流场特性并考虑了含灰气体激波的松弛结构对边界层流动的影响。     

Rapid cylindrically and spherically symmetric strong compression of a perfect gas

The problem of the rapid cylindrically and spherically symmetric strong compression of a perfect (non-viscous and non-heat-conducting) gas is solved. The term “rapid” denotes that the compression time is much less than the run time of a sound wave across the initial cylindrical or spherical volume, while the term “strong” in this case means the simultaneous attainment of as large a density and temperature as desired. By definition, rapid compression must begin in a strong shock wave, which propagates to the axis or centre of symmetry. When the shock wave approaches the centre of symmetry this flow is described by the self-similar Guderley equation with an unbounded rise in temperature, pressure and velocity and a finite increase in the density at the centre of symmetry both behind the arriving and behind the reflected shock waves. To obtain as high an increase in the density as desired one must add on a centred compression wave with focus at the centre of symmetry to the overtaking shock wave at the instant it arrives at the centre of symmetry C⁻-characteristic. Outside a small neighbourhood of the focus one can calculate, by the method of characteristics, the centred wave and the trajectory of the piston which produces it. As for any centred wave, this calculation must be carried out from the centre of symmetry. Since some of the parameters at the focus (certainly the pressure, temperature and velocity of the gas) are unbounded, it is necessary to preface the calculation by the method of characteristics by constructing an analytic solution which holds in a small neighbourhood of the centre of symmetry. Below, after constructing the required solution, the centred waves corresponding to it and the trajectories of the piston producing them are calculated.     

Isothermal flow behind a shock wave propagating in the solar wind

A Parker-type blast wave, which is headed by a strong shock, driven out by a propelling contact surface, moving into an ambient solar wind having a strictly inverse square law radial decay in density, is studied. Assuming the self-similar flow behind the shock to be isothermal, approximate analytical and exact numerical solutions are obtained. There is a good agreement between the approximate analytical and exact numerical solutions. It is observed that the mathematical singularity in density at the contact surface is removed for the isothermal flow.     

Flow Structure behind a Shock Wave in a Channel with Periodically Arranged Obstacles

We study the propagation of a pressure wave in a rectangular channel with periodically arranged obstacles and show that a flow corresponding to a discontinuity structure may exist in such a channel. The discontinuity structure is a complex consisting of a leading shock wave and a zone in which pressure relaxation occurs. The pressure at the end of the relaxation zone can be much higher than the pressure immediately behind the gas-dynamic shock. We derive an approximate formula that relates the gas parameters behind the discontinuity structure to the average velocity of the structure. The calculations of the pressure, velocity, and density of the gas behind the structure that are based on the average velocity of the structure agree well with the results of gas-dynamic calculations. The approximate dependences obtained allow us to estimate the minimum pressure at which there exists a flow with a discontinuity structure. This estimate is confirmed by gas-dynamic calculations.