Vapor condensation behind a shock wave propagating through a large molecular-mass medium

Measurements of film condensation were made behind the incident shock wave propagating through a vapor-liquid two-phase medium. Major objective of the study is to identify condensing heat transfer rates of the vapor to the shock-tube side wall as well as to learn the condensing main flow field. Ethanol and E-10 (a heavy liquid named Afluid by the manufacturer) were extensively used as working fluid. Steady accumulation of the condensing vapor was confirmed on the wall surface, as similarly seen in the end-wall experiment conducted elsewhere. A most significant result is that "dual-step" shock pressurization was observed in E-10. The first pressure rise is a normal one created by an incident shock front, whereas the second pressure rise is taken place by some large disturbance in the main flow. The reason for this is not certain yet, but is speculated to be a long relaxation time or inefficient compressibility of the fluid. The visualized shock front and its vicinity of E-10 is completely different from those of normal gases. Received May 31, 1994 / Accepted April 20, 1995     

Computation and experimental validation of N_{2}-CH_{4}-Ar mixture flows behind normal shock wave

Different vibration-dissociation-vibration coupling models have been used to compute the nonequilibrium N-CH-Ar mixture flow behind a normal shock wave. A three-temperature model was used and the diffuse nature of vibrational relaxation at high temperatures was accounted for. The numerical results obtained with the Treanor and Marrone model (preferential or non preferential) and the Park model of vibration-dissociation-vibration coupling are compared. These results show that the temperatures and the concentrations are mainly affected by the value of the characteristic temperature U in the preferential model of Marrone and Treanor. An assessment of the more realistic model was realized by comparing numerical results with shock tube experiments. The influence of argon addition on the nonequilibrium emission of CN behind the shock wave was also numerically studied and compared to experimental measurements. Received 1 September 1995 / Accepted 10 December 1996     

Nonequilibrium ionizatlon of air behind a shock wave front at speeds of 5–10 km/sec

The solution of the problem of nonequilibrium ionization of dissociating air in a shock wave propagating with a speed of 5–10 km/sec has shown that the electron concentration distribution has a maximum behind the wave front for speeds below 9 km/sec.The formation of this maximum is caused by the high associative ionization rate in comparison with the nitrogen dissociation rate. In the nonequilibrium region behind the shock wave front there is formed a considerable concentration peak of the molecular ions which are formed as a result of associative ionization and charge exchange.The author wishes to thank V. P. Stulov for useful discussions.     

Vibrational temperature measurements in a shock layer using laser induced predissociation fluorescence

Single shot spatially and spectrally resolved laser induced predissociation fluorescence measurements in a shock layer around a cylinder in a pulsed supersonic free stream are presented. Fluoresence signals were produced using the tuned output of an argon fluoride excimer laser to excite a mixture of rovibrational transitions in molecular oxygen. The signals produced along a line inside the shock layer were focussed onto a two dimensional detector coupled to a spectrometer, thus allowing spectral and spatial resolution of the fluoresence. In this way, it was possible to detect two fluoresence signals from two different transitions simultaneously, allowing the determination of vibrational temperatures without the need for calibration. To minimize problems associated with low signal to noise ratios, background subtraction and spatial averaging was required.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 formation of a secondary shock wave behind a shock wave diffracting at a convex corner

This paper deals with the formation of a secondary shock wave behind the shock wave diffracting at a two-dimensional convex corner for incident shock Mach numbers ranging from 1.03 to 1.74 in air. Experiments were carried out using a 60 mm 150 mm shock tube equipped with holographic interferometry. The threshold incident shock wave Mach number () at which a secondary shock wave appeared was found to be = 1.32 at an 81° corner and = 1.33 at a 120° corner. These secondary shock waves are formed due to the existence of a locally supersonic flow behind the diffracting shock wave. Behind the diffracting shock wave, the subsonic flow is accelerated and eventually becomes locally supersonic. A simple unsteady flow analysis revealed that for gases with specific heats ratio the threshold shock wave Mach number was = 1.346. When the value of is less than this, the vortex is formed at the corner without any discontinuous waves accompanying above the slip line. The viscosity was found to be less effective on the threshold of the secondary shock wave, although it attenuated the pressure jump at the secondary shock wave. This is well understood by the consideration of the effect of the wall friction in one-dimensional duct flows. In order to interpret the experimental results a numerical simulation using a shock adaptive unstructured grid Eulerian solver was also carried out. Received 1 May 1996 / Accepted 12 September 1996     

Particle and drop dynamics in the flow behind a shock wave

The results of an investigation of the dynamics of hard particles and liquid drops in the flow behind a transmitted shock wave are presented. From the equation of motion of a particle in the shock wave, relations for the displacement, velocity and acceleration as functions of time and certain velocity-relaxation parameters taking into account the properties of the gas and the aerodynamic drag of the particles are obtained for unsteady flow around the particles at an acceleration of 10³–10⁴ m/s². It is shown that the velocity-relaxation parameters are universal. Approaches to finding the aerodynamic drag of freely-accelerating bodies from the dynamics of their acceleration after being suddenly exposed to the flow are considered. It is established that under these conditions the drop dynamics observed can be well described in terms of the same velocity-relaxation parameters with account for linear growth of the transverse drop size. All the kinematic functions obtained are confirmed experimentally.     

Formation of a jet and vortices behind a shock wave reflected from a concave body

A jet and vortices have been observed when a plane shock wave reflects from a concave body in a shock tube. If the cavity is deep enough then two reflected shocks appear near its edges. Air, carbon tetrafluoride (CF) and dichlorodifluoromethane (CClF) were chosen as test gases. The flow was visualized with the aid of a conventional shadow technique. Pressure measurements at the body surface were also obtained. Numerical studies have been conducted using a two-dimensional inviscid model. There is a good qualitative agreement between the experimental and numerical results. Received 8 February 1996 / Accepted 30 June 1997     

Experimental and numerical studies of rotational relaxation behind a strong shock wave in air

This paper describes the experimental and numerical investigations of unknown characteristics of the rotational nonequilibrium phenomena behind a strong shock wave in air. Experiments were carried out using a piston-driven shock tube with helium as driving gas and air as driven (test) gas, operated as a two-stage shock tube. In the experiments, emission spectra of NO were measured to evaluate the rotational temperature behind a strong shock wave. The numerical calculations use the computational code for the thermal and chemical nonequilibrium flow behind a strong shock wave developed by the present author's group, where 11 chemical species (N, O, NO, N, O, N, O, NO, N, O, e) and 48 chemical reactions of high-temperature air are considered. The thermal nonequilibrium is expressed by introducing an 8 temperature model composed of translational temperature, rotational and vibrational temperatures for N, O, NO, and electron temperature. The coupling of a rotation, vibration and dissociation (CRVD) model was incorporated to take sufficiently into account the rotational nonequilibrium. The calculations were conducted for the same conditions as the experimental ones. From the calculated flow properties, emission spectra were re-constructed using the code for computing spectra of high temperature air “SPRADIAN”. Furthermore, rotational and vibrational temperatures of NO (0,1) were determined from a curve fitting method and compared with the computed results. Received 12 September 2001 / Accepted 18 February 2002     

The flow structure of dilute gas-particle suspensions behind a shock wave moving along a flat surface

The asymptotic and numerical investigations of shock-induced boundary layers in gas-particle mixtures are presented.The Saffman lift force acting on a particle in a shear flow istaken into account.It is shown that particle migration across the boundary layer leads tointersections of particle trajectories.The corresponding modification of dusty gas model isproposed in this paper.The equations of two-phase sidewall boundary layer behind a shock wave moving at aconstant speed are obtained by using the method of matched asymptotic expansions.Themethod of the calculation of particle phase parameters in Lagrangian coordinates isdescribed in detail.Some numerical results for the case of small particle concentration aregiven.     

Investigation of the acceleration of aluminum particles behind a shock wave using instantaneous Laser Doppler Velocimetry

The acceleration of aluminum particles with a 5μm diameter in the flow field behind an incident shock wave was investigated experimentally in a 10-m long and 70 mm inner diameter shock tube. By means of instantaneous Laser Doppler Velocimetry (LDV) the velocity of the particles was observed directly. The light scattered by the moving particles is Doppler shifted and sent to the laser Doppler velocimeter. The velocimeter essentially consists of a phase-stabilized Michelson interferometer used as a sensitive spectrometer. An electro-optical circuit ensures the phase stabilization that results in a voltage signal independent of the scattered light intensity and proportional to the mean velocity of the particles at the measurement point. Because of the very short response time (1μs) of the LDV system used here, the latter gives a continuous real-time signal of the particle acceleration. To avoid particle oxidation the particles were accelerated by a high-speed nitrogen gas flow. From the measured velocity the dimensionless drag coefficient was calculated. The drag coefficient is related to the fluid dynamic force exerted by the gas on the particles. The experimental data were compared to theoretical models from the literature. A significant deviation between the model and the experimental data was observed. This deviation is supposed to be induced by the shock wave, which hits the particles and breaks them into pieces of a smaller diameter. Further experiments will be carried out in the future to check the size distribution of the particles after the shock has gone past them.      

The effect of an unsteady drag force on the structure of a non-equilibrium region behind a shock wave in a gas-particle mixture

For numerical analysis of shock wave propagation in gas-particle mixtures, drag coefficients of a sphere in steady flows are generally used. However, it is shown both experimentally and numerically that a shock loaded solid sphere experiences unsteady drag forces. The paper describes a model of unsteady drag force and its effect on the structure of the non-equilibrium region behind a shock front traveling in a dusty gas. The results are compared with those obtained by using a steady drag coefficient and are discussed. It is demonstrated that the large drag force at the early stage of the interaction between shock-wave induced flow and a solid particle affects the flow structure that is obtained with a steady drag force.      

Axisymmetric shock wave interaction with a cone: a benchmark test

Results of the benchmark test are presented of comparing numerical schemes solving shock wave of M_s = 2.38 in nitrogen and argon interacting with a 43∘ semi-apex angle cone and corresponding experiments. The benchmark test was announced in Shock Waves Vol. 12, No. 4, in which we tried to clarify the effects of viscosity and heat conductivity on shock reflection in conical flows. This paper summarizes results of ten numerical and two experimental applications. State of the art in studies regarding the shock/cone interaction is clarified. PACS 01.50.Kw, 47.15.Pn Communicated by K. Takayama     

Processes of hydrate formation and dissolution behind a shock wave in a liquid containing gas bubbles (mixture of nitrogen and carbon dioxide)

The processes of dissolution and hydrate formation behind a moderate-amplitude shock wave in water containing gas bubbles (mixture of nitrogen and carbon dioxide) are studied in experiments with different initial static pressures in the medium and concentrations of carbon dioxide in bubbles. An increase in static pressure in the gas-liquid medium is demonstrated to enhance the influence of the non-reacting gas (nitrogen) on the processes of dissolution and hydrate formation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 178–187, March–April, 2009.     

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     

Numerical simulation for nonstationary Mach reflection of a shock wave: A kinetic-model approach

A numerical simulation was performed for the process of formation of single Mach reflection on a wedge by solving a BGK type kinetic equation for the reduced distribution function with a finite difference scheme. The calculations were carried out for a shock Mach number 2.75 and wedge angle 25° in a monatomic gas, which corresponds to the conditions of single Mach reflection in the classical von Neumann theory. The calculations were performed for both diffuse and specular reflection of molecules at the wall surface. It is concluded that the diffuse reflection of molecules at the wall surface or the existence of the viscous or thermal layer is an essential factor for a nonstationary process at the initial stage of Mach reflection. Furthermore, the numerical results for diffuse reflection are found to simulate the experimental results very well, such as a transient process from regular reflection to Mach reflection along with shock propagation.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Dynamics of stress wave propagation in a chain of photoelastic discs impacted by a planar shock wave: Part II, numerical investigation

The propagation of stress waves through a chain of discs has been studied experimentally in Part I (Glam et al. [1]) and is completed here with numerical investigation using the standard package ABAQUS. A fair agreement is found between experimental findings and their simulations. Based on this agreement, parametric study of wave propagation through disc-chains was conducted. Specifically, effects associated with changes in the disc diameter, material density, stiffness/rigidity and the number of discs in the chain on the stressed chain have been studied. It was found that the propagation velocity of the evolved waves increases with improving contacts between the chain’s discs by exposing the chain to a static load before its dynamic loading. The wave- propagation velocity decreases with increase in the discs material density and it increases when its diameter increases. In case of a chain composed of small diameter discs and/or small material density, the transmitted stress wave is first strengthened and only at discs further down the chain it starts decaying. When checking the influence of the dynamic-loading duration it was found that long dynamic-load duration dissolves quickly into short pulses. It was also found that there is a ‘characteristic’ wave for a given chain. This wave propagates with minimal dispersion. Dynamic loads having shorter time duration than the ‘characteristic’ one experiences significant attenuation.