Modeling of vibration-dissociation oxygen kinetics at temperatures of 4,000-11,000 K

The time profiles of vibrational molecular oxygen temperature T _v measured earlier in experiments behind a strong shock wave were used for testing the theoretical and empirical models of thermal nonequilibrium dissociation of molecules. To do this, dissociating gas flows behind the strong shock wave front were calculated with account for these models. If the initial gas temperature behind the wave front T ₀ < 6.5×10³ K, the models well describe changing the temperature with time. However, for T ₀ > 7×10³ K neither of the models tested describes the measured temperature profiles satisfactorily. Using the empirical model proposed in the present study made it possible to satisfactorily describe the vibrational temperature evolution observed in experiments at temperatures up to 11×10³ K.     

Effect of precursor radiation on the flow structure and ionization behind the shock front in inert gases

The problem of the propagation of strong, intensely radiating shock waves in inert gases is considered. It is shown that the heating of the shock tube walls by the precursor radiation, accompanied by an increase in the temperature of the adjacent gas, leads to the transverse stratification of the medium and to the disturbance of the one-dimensionality of the flow of shock-heated gas behind the wave front. Ionization kinetics calculations which take this into account indicate an acceleration of ionization near the tube walls, which is consistent with experiment. On the basis of the gas heating values obtained it is possible to establish critical values of the gas pressures ahead of the front and the shock wave Mach numbers, transition through which is accompanied by a radical restructuring of the flow with the formation of a configuration.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 124–131, May–June, 1991.     

Effect of leading radiation on the stability of strong shock waves in gases with an arbitrary equation of state

The effect of leading radiation on the stability of a strong shock wave in an ideal gas with an arbitrary equation of state is investigated. The ionization ahead of and behind the shock front and the radiation are assumed to be in equilibrium. The investigation is carried out in the linear approximation with respect to amplitude for disturbances with a wavelength much greater than the width of the relaxation zones ahead of and behind the shock. The conditions under which the leading radiation has a destabilizing effect on the shock wave are established. It is shown, in particular, that neutrally stable shock waves become unstable. The conditions under which the onset of instability is of the threshold type with respect to the radiation intensity are determined. It is found that the radiation also has a destabilizing effect on stable shock waves, including shock waves in a perfect gas. However, in this case instability can develop only when the disturbances have a wavelength comparable with the width of the relaxation zone. A simple physical mechanism of the onset of instability under the influence of leading radiation is proposed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 125–133, May–June, 1990.The authors are grateful to A. G. Kulikovskii and A. A. Barmin for their constant interest and useful discussions.     

Specific features of incident and reflected blast waves

On the basis of numerical modeling specific features of shock wave reflections were analyzed. It was found, that after diaphragm rupture self-modeling pressure and velocity distributions nearby the shock front establish. But in some special cases the temperature behind the shock front can rise. This peculiarity should be taken into account when performing experiments with high reactive gaseous mixtures. The temperature on the shock front and the velocity gradient behind it are uniform in the case of strong blast wave reflections. This effect is observed in the zone with an elevated temperature profile behind the incident blast wave. The reflected triangular waves conserve a quasi-self-modeling character by pressure. Typical experiments were carried out to verify the theoretical predictions. The effects of reflected wave acceleration in the case of triangular waves and the self-similar character of the pressure profiles were observed.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

The quasi-stationary structure of radiating shock waves. I. The one-temperature fluid

We calculate the quasi-stationary structure of a radiating shock wave propagating through a spherically symmetric shell of cold gas by solving the time-dependent equations of radiation hydrodynamics on an implicit adaptive grid. We show that this code successfully resolves the shock wave in both the subcritical and supercritical cases and, for the first time, we have reproduced all the expected features – including the optically thin temperature spike at a supercritical shock front – without invoking analytic jump conditions at the discontinuity. We solve the full moment equations for the radiation flux and energy density, but the shock wave structure can also be reproduced if the radiation flux is assumed to be proportional to the gradient of the energy density (the diffusion approximation), as long as the radiation energy density is determined by the appropriate radiative transfer moment equation. We find that Zel'dovich and Raizer's (1967) analytic solution for the shock wave structure accurately describes a subcritical shock but it underestimates the gas temperature, pressure, and the radiation flux in the gas ahead of a supercritical shock. We argue that this discrepancy is a consequence of neglecting terms which are second order in the minimum inverse shock compression ratio [, where is the adiabatic index] and the inaccurate treatment of radiative transfer near the discontinuity. In addition, we verify that the maximum temperature of the gas immediately behind the shock is given by , where is the gas temperature far behind the shock. Received 21 September 1998/ Accepted 2 February 1999     

Distribution of radiant thermal flux over the surface of a sphere for hypersonic flow of a nonviscous radiating gas past the sphere

In the present study using the Newtonian approximation [1] we obtain an analytical solution to the problem of flow of a steady, uniform, hypersonic, nonviscous, radiating gas past a sphere. The three-dimensional radiative-loss approximation is used. A distribution is found for the gasdynamic parameters in the shock layer, the withdrawal of the shock wave and the radiant thermal flux to the surface of the sphere. The Newtonian approximation was used earlier in [2, 3] to analyze a gas flow with radiation near the critical line. In [2] the radiation field was considered in the differential approximation, with the optical absorption coefficient being assumed constant. In [3] the integrodifferential energy equation with account of radiation was solved numerically for a gray gas. In [4–7] the problem of the flow of a nonviscous, nonheat-conducting gas behind a shock wave with account of radiation was solved numerically. To calculate the radiation field in [4, 7] the three-dimensional radiative-loss approximation was used; in [5, 6] the self-absorption of the gas was taken into account. A comparison of the equations obtained in the present study for radiant flow from radiating air to a sphere with the numerical calculations [4–7] shows them to have satisfactory accuracy.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 44–49, November–December, 1972.In conclusion the author thanks G. A. Tirskii and É. A. Gershbein for discussion and valuable remarks.     

Initial stage of the collision of blast waves

The collision of two blast waves is analyzed for the case of variable parameters of the gas behind the wave front and wave reflection at a plane, a cylindrical, and a spherical obstacle. The reflection of a blast wave from a nonmoving obstacle is investigated in detail. The problem of the collision of two shock waves with constant parameters behind the front is solved both in the symmetrical case (reflection from a nonmoving wall) and in the case of waves of different amplitudes by a system of algebraic relations for the compression shocks. The reflection of a strong point-source spherical shock wave from a wall has been treated in [1, 2]. The present article examines the initial stage of wave collision for an arbitrary distribution of the parameters behind the front.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 41–48, September–October, 1971.The authors are grateful to V. P. Korobeinikov for a discussion of the results and to V. P. Kolgan for furnishing the numerical solutions.     

Shock Compression of a Plate on a Wedged–Shaped Target

Problems of compression of a plate on a wedge–shaped target by a strong shock wave and plate acceleration are studied using the equations of dissipationless hydrodynamics of compressible media. The state of an aluminum plate accelerated or compressed by an aluminum impactor with a velocity of 5—15 km/sec is studied numerically. For a compression regime in which a shaped–charge jet forms, critical values of the wedge angle are obtained beginning with which the shaped–charge jet is in the liquid or solid state and does not contain the boiling liquid. For the jetless regime of shock–wave compression, an approximate solution with an attached shock wave is constructed that takes into account the phase composition of the plate material in the rarefaction wave. The constructed solution is compared with the solution of the original problem. The temperature behind the front of the attached shock wave was found to be considerably (severalfold) higher than the temperature behind the front of the compression wave. The fundamental possibility of initiating a thermonuclear reaction is shown for jetless compression of a plate of deuterium ice by a strong shock wave.     

Radiative cooling of plasma behind a reflected shock front

Calculation of gas flow in a shock tube on the basis of ideal theory [1] leads to results that differ from the real picture. In particular, the calculated velocity of the reflected shock wave exceeds the experimentally measured velocity [2] by about 20%. The calculated parameters of shock-heated gas agree well with the experimental results only directly behind the shock front [3]. The present paper reports a theoretical and experimental investigation of the variation of the plasma parameters behind the front of a reflected shock wave in argon. A picture of the gas-dynamic processes taking place after reflection of the incident shock wave by the end of the shock tube is determined. A method is developed for approximate analytic calculation, this making it possible to determine not only the parameters of the gas directly behind the front of the reflected shock wave for different positions of the wave relative to the end of the shock tube but also the variation of these parameters in other regions behind the reflected shock wave. The calculation takes into account the influence of the boundary layer and radiative cooling in the approximation of a low degree of ionization of the plasma and persistence of equilibrium conditions in the entire region behind the reflected shock wave. The experimental and theoretical profiles of the radiation behind the reflected shock wave are compared.     

Unsteady viscous shock layer in the case of supersonic motion through an inhomogeneity

The model of a perfect gas is used in a numerical simulation of unsteady effects in a viscous shock layer near the stagnation streamline near the front part of a rotating blunt body in an inhomogeneous external flow with pressure difference P _t8, temperature difference T _t8, and vorticity difference . The evolution of nonlinear disturbances due to the passage of a heated region and a change of the injection regime is followed. A divergence finite-difference scheme of second order of approximation across the shock layer, realized by vector sweeps with allowance for the boundary conditions on the surface of the body and behind the separated bow shock wave, is used.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 138–145, January–February, 1992.     

Supersonic flow of combustible gas mixture past sphere with account for ignition delay time

Assume an axisymmetric blunt body or a symmetric profile is located in a uniform supersonic combustible gas mixture stream with the parameters M₁, p₁, and T₁. A detached shock is formed ahead of the body and the mixture passing through the, shock is subjected to compression and heating. Various flow regimes behind the shock wave may be realized, depending on the freestream conditions. For low velocities, temperatures, or pressures in the free stream, the mixture heating may not be sufficient for its ignition, and the usual adiabatic flow about the body will take place. In the other limiting case the temperature behind the adiabatic shock and the degree of gas compression in the shock are so great that the mixture ignites instantaneously and burns directly behind the shock wave in an infinitesimally thin zone, i. e., a detonation wave is formed. The intermediate case corresponds to the regime in which the width of the reaction zone is comparable with the characteristic linear dimension of the problem, for example, the radius of curvature of the body at the stagnation point.The problem of supersonic flow of a combustible mixture past a body with the formation of a detonation front has been solved in [1, 2]. The initial mixture and the combustion products were considered perfect gases with various values of the adiabatic exponent .These studies investigated the effect of the magnitude of the reaction thermal effect and flow velocity on the flow pattern and the distribution of the gasdynamic functions behind the detonation wave.In particular, the calculations showed that the strong detonation wave which is formed ahead of the sphere gradually transforms into a Chapman-Jouguet wave at a finite distance from the axis of symmetry. For planar flow in the case of flow about a circular cylinder it is shown that the Chapman-Jouguet regime is established only asymptotically, i. e., at infinity.This result corresponds to the conclusions of [3, 4], in which a theoretical analysis is given of the asymptotic behavior of unsteady flows with planar, spherical, and cylindrical detonation waves.Available experimental data show that in many cases the detonation wave does not degenerate into a Chapman-Jouguet wave as it decays, bur rather at some distance from the body it splits into an adiabatic shock wave and a slow combustion front.The position of the bifurcation point cannot be determined within the framework of the zero thickness detonation front theory [1], and for the determination of the location of this point we must consider the structure of the combustion zone in the detonation wave. Such a study was made with very simple assumptions in [5].The present paper presents a numerical solution of the problem of combustible mixture flow about a sphere with a very simple model for the structure of the combustion zone, in which the entire flow behind the bow shock wave consists of two regions of adiabatic flow-an induction region and a region of equilibrium flow of products of combustion separated by the combustion front in which the mixture burns instantaneously. The solution is presented only for subsonic and transonic flow regions.     

Experimental investigation of the emission in the unsteady flow region behind strong shock fronts in mixtures of co and CO2 with nitrogen and argon

Mixtures of CO (or CO₂) gases and N₂ behind strong shock fronts at temperatures 4000–10 000 ° K have been investigated with a view to elucidating the mechanism of the physicochemical processes in the unsteady region of the gas flow behind a shock front leading to the behavior of strongly radiating CN and C₂ molecules and C atoms and also determining the quantitative characteristics of the chemical reactions. A shock tube was used in the investigations, which made it possible to obtain the intensity distribution of the radiation of several components — CN, C₂, and C — behind shock fronts.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 120–129, March–April, 1981.We thank S. A. Losev and O. P. Shatalov for assisting in the work and for valuable discussions.     

Shock wave structure in a mixture of gas,vapour and droplets

The structure of the relaxation zone behind a shock wave of moderate strength in a mixture of gas, vapour and droplets is analysed. A model is presented for shock induced evaporation, which is based on wet-bulb equilibrium and on the absence of relative motion between droplets and gas. Experimental and numerical data on heterogeneous condensation induced by an unsteady rarefaction wave and on re-evaporation due to shock wave passage are reported for a mixture of water vapour, nitrogen gas and condensation nuclei. Pressure, temperature, saturation ratio and droplet size are experimentally obtained and are very well predicted by a numerical simulation based on the non-linear quasisteady wet-bulb model for phase transition, as well for the expansion wave as for the shock wave. During expansion, droplet number density decays much faster than predicted, which is not yet satisfactorily explained. Shock induced droplet evaporation is studied for post-shock saturation ratios ranging from 5×10^–3 to 0.2, corresponding to shock Mach numbers of 1.2 to 1.9. The evaporation times are well predicted by the theoretical model. No evidence is found for droplet break-up for Weber numbers up to 13, and droplet radii of the order of 1m.On leave at Institute of Fluid Science, Shock Wave Research Center, Tohoku University, Sendai 980, JapanThis article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Nonequilibrium molecular radiation behind the front of a strong shock wave in the CO<Subscript>2</Subscript>-N<Subscript>2</Subscript>-O<Subscript>2</Subscript> mixture

Models of population of some radiating electron-vibrational states of CO, CN, and C₂ molecules are developed. The characteristics of radiation in a chemically nonequilibrium flow behind the front of a strong shock wave in a mixture of gases constituting the Martian atmosphere are calculated. The numerical data are compared with experimental results.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 2, pp. 13–22, March–April, 2005     

Influence of the temperature of the surface of an axisymmetrical body and forward radiation on the distribution of a radiant flux on its surface in hypersonic flow-past

A solution is obtained of the flow-past problem for an axisymmetrical body with steady-state hypersonic nonviscous, space-radiating gas flow in a hypersonic approximation. It is shown as illustrated by the example of flow-past of a sphere by an air flow, that the relative distribution of the radiant flux weakly depends on a calculation of surface re-radiation, while the size of the radiant flux substantially depends on body temperature T_W at a critical point. The distributions of radiant flux for sphere flow-past by a CO₂-N₂ gas mixture (at T_W = 0) are calculated using a previously developed method. It is shown that different CO₂ contents in the initial mixture of the incident gas flow weakly affect this distribution. The dependence of the distribution of the radiant flux and departure of the shock wave on the boundary condition for gas enthalpy in the pressure shock, taking into account forward radiation, is investigated. Asymptotic expressions are obtained for sphere flow-past for the case of a strongly radiating gas. Distributions of the radiant flux for different assumptions for the boundary conditions in shocks are calculated.     

Calculation of the motion of a gas behind the front of a luminous detonation wave taking account of the lateral expansion of the plasma column

Nonsteady-state gasdynamic processes are considered in a plasma column (hot channel), formed behind the front of a shock wave moving toward a laser beam. A quasi-one-dimensional approximation is used-the parameters in the channel are assumed to be compensated with respect to cross section, but depend on the time and distance along the axis. Motion in the cold dense casing surrounding the channel is assumed to be one-dimensional and cylindrically symmetrical. The solutions of the corresponding systems of equations in partial derivatives permit the parameters to be determined approximately both in the case when the mean free path of the radiation is small in comparison with the radius of the beam (luminous detonation) and also in the case when the mean free path is comparable with the radius. Examples are given of the corresponding numerical calculations. It is shown that in the latter case, a cycle of incomplete absorption can be achieved when behind the shock wave front, moving with constant velocity up to the Jouguet plane, only a part of the radiation energy incident on the front is released.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 18–28, May–June, 1976.     

Regular reflection of a shock wave from a solid wall in a stream of fuel mixture

The two-dimensional stationary problem of regular reflection of a shock wave from a plane solid wall in a fuel gas mixture is examined in the case when the mixture is ignited at the intersection of the incident wave with the wall and a flame front is formed behind the reflected shock wave. The shock waves and the flame front are considered plane surfaces of discontinuity. The fuel mixture and the reaction products are considered perfect, inviscid, and non-heat-conducting gases.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 161–163, July–August, 1978.     

Unsteady-state flow of a gas in a channel,resulting from the motion of a piston,with magnetohydrodynamic offtake of energy and luminescence of the gas

The present paper discusses the one-dimensional unsteady-state flow of a gas resulting from the motion of a piston in the presence of weak perturbing factors, with which the investigation of the perturbed (with respect to the usual self-similar conditions) motion reduces to the solution of ordinary differential equations, is indicated. The distributions of the parameters of the gas between the piston and the shock wave are found. The conditions under which there is acceleration or slowing down of the shock front are clarified. As an example, this paper considers the unsteady-state motion of a conducting gas in a channel with solid electrodes under conditions where electrical energy is generated, and the flow of a gas taking radiation into account, under the assumption of optical transparency of the medium. The theory developed is used to solve the problem of the motion of a thin wedge with a high supersonic velocity in an external axial magnetic field, taking account of the luminescence of the layer of heated gas between the wedge and the shock wave.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 17–25, September–October, 1970.     

Propagation of shock waves in free space

The problem of the exit of a shock wave from an axisymmetric channel and its propagation in a free space occupied by an ideal gas is examined. This problem has been studied earlier in [1], in which the shock wave front was considered planar, as well as in [2], in which the wave front was regarded as a surface of an ellipsoid of revolution. The solutions obtained in these studies assumed the presence of two regions in the wave-front surface: the region of the original shock wave and a region stemming from the decomposition of an infinitesimally thin annular discontinuity of the gas parameters, with the wave intensity over the front surface in each region being considered constant, i.e., the wave character of the process over the front was not considered. In this study a solution will be achieved by the method of characteristics [3–5] of the equations of motion of the shock-wave front, as obtained in [6, 7]. Flow fields are determined for the region immediately adjacent to the shock-wave front for a wide range of shock-wave Mach numbers M _a =1.6–20.0 for = 1.4. On the basis of the data obtained, by introduction of variables connected with the length of the undisturbed zone, as calculated from the channel cross-section along the x axis, together with the pressure transition at the wave front, relationships are proposed which approach self-similarity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 163–166, September–October, 1971.In conclusion, the author thanks S. S. Semenov for his valuable advice on this study.     

Parameters of shock waves with emission in gases

Parameters of emitting shock waves in gases are investigated in the limiting case when there is no screening of emission from the shock front by the precursory layer. The one-dimensional quasi-steady-state formulation of the problem with deceleration of high-speed gas flow against a plane fixed obstacle under conditions of strong emission is given. The case of the shock waves of large optical thickness is analytically considered over a wide range of variation of the obstacle reflectivity. The parameters of emitting shock waves generated in experiments in shock tubes in the inert argon gas are estimated using the methods developed and compared with the measurement results. The shock “adiabats” of optically thick shock waves are considered with allowance for the radiation energy losses. The calculations are carried out for aluminium plasma.