The moiré method for the study of explosions

Explosives have been used to produce weak shock waves in central cavities of Plexiglas spheres. Shock pressures engendered were less than the so-called stable shock threshold and they resulted in a decomposition of the shock wave into an elastic precursor preceding the so-called plastic shock wave. The differential moiré method has been used to measure simultaneously the precursor and shock-wave velocities, as well as the velocity of the cavity interface. Moreover, the displacement field given by the moiré pattern yields the components of strains and the particle velocity behind each shock.     

Unloading effects in the dynamic response of granular soil

Experimental studies are described in which two different types of air shock are used to induce stress waves in an elastically confined, cylindrical specimen of sand, the objective being to permit direct observation of the effects of unloading on the dynamic response of granular material. In one case, “pure” unloading waves are generated by suddenly subjecting the end of a prestressed specimen to pressures below atmospheric. In the second type of test, the specimens are impacted at one end by an air shock which has the form of a sharp front followed immediately by a rapid exponential decay of pressure. The decay of the shock pressure generates unloading waves which cause attenuation of the wave as it propagates in the granular material. As a result of the tests, it is suggested that seismic-wave velocity may be correlated to the modulus of unloading cycles in a quasi-static test and that, when strainrate effects are included, a relatively simple model may be used to predict the dynamic response of granular materials.     

Numerical simulation of shock wave intersections

A large number of papers, generalized and classified in [1, 2], have been devoted to unsteady gas flows arising in shock wave interaction. Experimental results [3–5] and theoretical analysis [6–9] indicate that the most interesting and least studied types of interaction arise in cases when there are several shock waves. At the same time, nonlinear effects, which depend largely on the nature of the shock wave intersections, become appreciable. Regions of existence of different types, of plane shock wave intersections have been analyzed in [10–13]. It has been shown that in a number of cases the simultaneous existence of different types of intersections is possible. The aim of the present paper is to study unsteady shock wave intersections in the framework of a numerical solution of the axisymmetric boundary-value problem that arises in the diffraction of a plane shock wave on a cone in a supersonic gas flow. Flow regimes that augment the experimental data of [3–5] and the theoretical analysis of [9] are considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 134–140, September–October, 1986.     

Determination of the minimum entropy configuration and velocity surfaces for elastic–plastic wave scattering at grain boundaries

A recent model based on full elastic anisotropy and crystal plasticity predicted the existence of multiple wave configurations during the interaction of stress waves with grain boundaries. Since the multiple wave configuration scenario cannot exist in nature, the principle of minimum entropy production is applied in the current work to find the most probable configuration. A large amplitude transmitted quasi-longitudinal wave is predicted for the given bicrystal orientation studied due to the wave propagating near a [001] direction and thus requiring large stress given the very low Schmid factor in this direction (for nickel aluminide (NiAl) as a model material). Anisotropic elastic–plastic velocity surfaces for quasi-longitudinal and quasi-shear waves in NiAl have also been constructed to gain an understanding of the general nature of plastic waves as a function of crystallographic direction.     

Dilatational-wave-induced pore-water pressure in soil

This paper examines peak and residual pore-water pressures in water-saturated soil induced by a dilatational stress wave. Our new laboratory testing device applies submillisecond, high pressure dilatational stress-wave loadings to water-saturated soil. The soil's initial effective stress, density, back pressure and saturation can be controlled with our device. Experimental results show that it is possible to induced residual excess pore-water pressure and liquefaction in water-saturated Monterey No. 0/30 sand. Liquefaction is induced with compressive strains exceeding 0.1 percent for loose samples consolidated at 172 kPa and 1 percent for dense samples consolidated at 690 kPa. Below a threshold compressive strain of about 0.005 percent, no significant residual excess pore-water pressures are developed.     

Establishment of circumfluence (steady-state flow) when a shock wave strikes a cylinder or a sphere

The problem of the diffraction of a shock wave at a stationary sphere or cylinder is considered. The finite-difference method proposed by S. K. Godunov [1, 2] is employed Numerical solutions are obtained for the stage of the diffraction of the shock wave and for the subsequent steady state of flow around the object (circumfluence). Cases of sub-, trans-, and supersonic flow behind the shock wave are considered. When strong shock waves undergo diffraction, zones of reverse flow appear in the neighborhood of the tail part of the obstacle.Moscow. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 97–103, September–October, 1972.     

Transmission of waves in a liquid through absorbing layers and subsequent interaction of the waves with an obstacle

The propagation of waves in porous media is investigated both experimentally [1, 2] and by numerical simulation [3–5]. The influence of the relaxation properties of porous media on the propagation of waves has been investigated theoretically and compared with experiments [3, 4]. The interaction of a wave in air that passes through a layer of porous medium before interacting with an obstacle has been investigated with allowance for the relaxation properties [5]. In the present paper, in which the relaxation properties are also taken into account, a similar investigation is made into the interaction with an obstacle of a wave in a liquid that passes through a layer of a porous medium before encountering the obstacle.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 53–53, March–April, 1983.     

Application of unsteady mixing equations to some aerodynamic problems

Tangential discontinuities [1] are introduced in solving several transient and steady-state problems of gas dynamics. These discontinuities are unstable [2] as a result of the effects of viscosity and thermal conductivity. Therefore it is advisable to replace the tangential discontinuity by a mixing region and account for its interaction with the inviscid flows, establishing on the boundaries of this region the conditions of vanishing friction stress and equality of the velocity and temperature components to the corresponding velocity and temperature components of the inviscid flows. This formulation improves the accuracy of the solution of such problems by posing them as problems with irregular reflection and intersection of shock waves [1].The consideration of the interaction of unsteady turbulent mixing regions with the inviscid flow also permits the formulation of several problems in which the effects of viscosity lead to complete rearrangement of the flow pattern (the lambda-configuration) with the interaction of the reflected shock wave with the boundary layer in the shock tube [3,4], the formation of zones of developed separation ahead of obstacles, etc.).In this connection, §1 presents an analysis of the self-similar solutions of the unsteady turbulent mixing equations (a corresponding analysis of the laminar mixing equations which coincide with the boundary layer equations is presented in [1]). It is shown that these self-similar solutions describe, along with the several problems noted above, the problems of the formation of steady jets and mixing zones in the base wake.As an example, §2 presents, within the framework of the proposed schematization, an approximate solution of the problem of the interaction of a shock wave reflected from a semi-infinite wall with the boundary layer on a horizontal plate behind the incident shock wave. The results obtained are applied to the analysis of reflection in a shock tube. Computational results are presented which are in qualitative agreement with experiment [3, 4].     

Gas content factor in the interaction of shock waves of different intensity in a two-phase gas-liquid medium

The transition from regular to Mach interaction is investigated in connection with the interaction of two plane weak or moderate shock waves of different intensity in a two-phase gas-liquid medium over the entire range of gas contents. A nonmonotonic dependence of the transition limit and the flow parameters on the gas content is detected. The investigation extends the results of [1] corresponding to the reflection of a shock wave from a wall. At intermediate gas contents in the case of opposing shock waves, analogous to the normal reflection of a shock wave from a solid wall, the results are in agreement with [2]. In the case of weak shock waves non-linear asymptotic expansions [3] are employed. In the extreme cases of single-phase media the results coincide with the findings of [3, 4]. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 172–174, November–December, 1986.     

The use of a slow sound speed fluorocarbon liquid for shock wave research

Fluorocarbon liquids have very low sound speeds in comparison with water. The measurement of shock waves in water is complicated by its relatively high sound speed. This paper presents a fluorocarbon liquid with a sound speed of 655 m/s for use in liquid shock experiments. Experimental and numerical results of shock wave reflection from various parabolas and wedges are given. Experiments were performed in a vertical liquid shock tube. The properties including an equation of state for the liquid are given. Numerical simulations using this equation of state are performed using a finite element program. It is shown that the investigation of non-linearities in water will require shock tubes that can withstand high pressures. Due to the high B/A parameter for this fluorocarbon liquid, it is demonstrated that non-linearities can be achieved and studied at much lower pressures. Received 1 July 1996 / Accepted 26 September 1996     

Half a century of continuous shock interaction investigations in the Joint Institute for High Temperatures of Russian Academy of Sciences

This article describes the history of the investigations of shock wave interactions at the Physical Gasdynamic Department, starting from the early 50s of the last century, when the first research related to missile reentry was made. The review focuses on a number of topics studied over more than 50 years and includes the study of strong shock waves, where it is necessary to take into account the physicochemical transformations in gases, shock wave reflection, diffraction, interaction with the boundary layer and with the nozzle, as well as detonation wave formation and interactions. The investigation of shock wave interactions is a current topic at the Joint Institute for High Temperatures of the Russian Academy of Sciences. Some new results are observed: the formation of impulse jets and the self-ignition of a cold hydrogen jet, diffraction of 3D shock waves, and the effect of an impulse jet and diffracted shock wave on an obstacle.     

Interaction of a wave with an obstacle in a multicomponent two-phase medium

The author's model [1] of a multicomponent liquid medium with nonlinear limiting compression diagrams and constant coefficient of viscosity is improved by the introduction of a coefficient of viscosity that varies during the deformation. The new model is used to obtain a numerical solution to the problem of the propagation of a plane wave produced by a shock load and the interaction of the wave with a fixed obstacle. Such a problem was solved earlier [2] in the case of a viscous medium for linear diagrams of static and dynamic compression and constant coefficient of viscosity. It is shown that the nonlinearity of the diagram of static compression leads with increasing pressure first to an increase in the reflection coefficient and then to a decrease of it. If the load has a sufficient duration, the initial section of the obstacle is subject to a succession of several waves, the number of which increases with increasing duration and amplitude of the load. The calculation was made for glycerine with air bubbles. It is assumed that at pressures up to 400·10⁵ N/m² glycerine is a linearly elastic medium In this case, the dynamic compression diagram of the two-component glycerine—gas-bubble medium is linear.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 183–187, May–June, 1981.I thank Yu. A. Sozonenko for discussion and valuable comments.     

Shock-Induced Borehole Waves and Fracture Effects

We perform wave experiments using a vertical shock tube setup. Shock waves are generated by the rupture of a thin membrane. In the test section the incident pressure waves generate borehole-guided waves along water-saturated samples. The tube is equipped with side wall gages and a mobile pressure probe, so that the attenuation and reflection of the wave can be measured. The computation for a single horizontal fracture intersecting a vertical borehole gives a quantitative prediction of reflection and transmission of borehole-guided waves. Three different fracture apertures are used for the calculation. Fracture aperture significantly affects both reflection and transmission coefficients. Large fractures increase reflectivity and decrease transmissivity. In the experiment, we found that both pressures above and below the fracture are influenced by the fracture aperture indeed, thus indicating the potential for fracture detection by borehole waves.     

Shock-Induced Borehole Waves and Fracture Effects

We perform wave experiments using a vertical shock tube setup. Shock waves are generated by the rupture of a thin membrane. In the test section, the incident pressure waves generate borehole-guided waves along water-saturated samples. The tube is equipped with side wall gages and a mobile pressure probe, so that the attenuation and reflection of the wave can be measured. The computation for a single horizontal fracture intersecting a vertical borehole gives a quantitative prediction of reflection and transmission of borehole-guided waves. Three different fracture apertures are used for the calculation. Fracture aperture significantly affects both reflection and transmission coefficients. Large fractures increase reflectivity and decrease transmissivity. In the experiment, we found that both pressures above and below the fracture are influenced by the fracture aperture indeed, thus indicating the potential for fracture detection by borehole waves.     

Interaction Between a Slow Magnetohydrodynamic Shock Wave and a Tangential Discontinuity

The self-similar problem of the oblique interaction between a slow MHD shock wave and a tangential discontinuity is solved within the framework of the ideal magnetohydrodynamic model. The constraints on the initial parameters necessary for the existence of a regular solution are found. Various feasible wave flow patterns are found in the steady-state coordinate system moving with the line of intersection of the discontinuities. As distinct from the problems of interaction between fast shock waves and other discontinuities, when the incident shock wave is slow the state ahead of it cannot be given and must to be determined in the process of solving the problem. As an example, a flow in which the slow shock wave incident on the tangential discontinuity is generated by an ideally conducting wedge located in the flow is considered. The basic features of the developing flows are determined.     

Study of unsteady flow past bodies behind shock waves

Several theoretical and experimental studies have been devoted to the problem of the nonstationary action of the stream behind a shock wave on bodies of varied shape. In particular, in [1], the pressure and density are calculated for flow about bodies of the more typical shapes in the initial stage of the process. The basic relations which accompany the interaction of shock waves are considered in [2, 3]. The analysis of the phenomena of diffraction of shock waves on the sphere, cylinder, and cone is presented in [4]. Problems of unsteady flow about a wing are examined in [5, 6]. A detailed review of the foreign studies on unsteady flow is given in [7]. Of great practical interest is the question of the time for flow formation and the magnitudes of the unsteady loads during this period. Experimental investigations have been made recently [8, 9] in which some criteria are presented for estimating the bow shock formation time for supersonic flow about the sphere and cylinder with flat blunting. However the question of the formation time of the stationary pressure on the body surface is not referred to in these studies and no relationship is shown between the transient position of the reflected wave and the corresponding unsteady pressure on the surface. Moreover, in [8] the dimensionless time criterion is determined very approximately, independently of the Mach number of the shock wave. The present study was undertaken with the object of determining the basic criteria which characterize unsteady flow about bodies behind a plane shock wave which has time-independent parameters, and clarification of the shock wave reflected from the body and the pressure on the surface of the body during the transient period. The most typical body shapes were studied: 1) a cylinder with flat face aligned with the stream; 2) a spherically-blunted cylinder; and 3) a cylinder transverse to the stream. The experiments were conducted in a conventional shock tube using the single-diaphragm scheme. The measurements of the pressure on the models and the velocity of the incident shock wave were made using the technique analogous to that of [10, 11]. A highspeed movie camera was used to record the pattern of the wave diffraction on the body. The Mach number of the incident shock wave varied in the range from M=1.5 to M≈6.0, which corresponded to a range of Mach numbers M_∞ of the stream behind the shock wave from 0.6 to 2.1. The calculations of the required gas dynamic parameters for high temperatures were made with account for equilibrium dissociation of the air on the basis of the data of [10, 12, 13]. The magnitude of the relative maximal shock wave standoff Δ at the stagnation point obtained in the present experiments was compared with the values of Δ from other studies. In the case of the flat-blunted cylinder it was in good agreement with the results of [8–14], and in the case of the spherically-blunted cylinder and the transverse cylinder it was in agreement with the results of [15].     

Shock wave propagation in a branched duct

The propagation of a planar shock wave in a 90° branched duct is studied experimentally and numerically. It is shown that the interaction of the transmitted shock wave with the branching segment results in a complex, two-dimensional unsteady flow. Multiple shock wave reflections from the duct's walls cause weakening of transmitted waves and, at late times, an approach to an equilibrium, one-dimensional flow. While at most places along the branched duct walls calculated pressures are lower than that existing behind the original incident shock wave, at the branching segment's right corner, where a head on-collision between the transmitted wave and the corner is experienced, pressures that are significantly higher than those existing behind the original incident shock wave are encountered. The numerically evaluated pressures can be accepted with confidence, due to the very good agreement found between experimental and numerical results with respect to the geometry of the complex wave pattern observed inside the branched duct. Received 15 July 1996 / Accepted 20 February 1997