Structure and evolution of shock waves in relativistic magnetohydrodynamics

We derive the existence conditions for relativistic shock waves propagating in a perfectly conducting fluid with a general equation of state that guarantees that the stationary wave has a continuous profile in the presence of weak viscosity. To this end we study the one-dimensional solutions of the magnetohydrodynamic equations with a relativistic viscosity tensor. We allow for anomalous regions of thermodynamic variables and do not use the well-known condition for the convexity of Poisson adiabats. The results lead to relationships among the velocities of magnetoacoustic, Alfvén, and shock waves in front of and behind the discontinuity that prove to be more stringent than the corollaries of the evolution conditions. In the nonrelativistic case and in parallel and perpendicular shock waves, any difference between the two conditions disappears. Zh. éksp. Teor. Fiz. 114, 881–891 (September 1998)     

Dynamics and Stability of a Weak Detonation Wave

One dimensional weak detonation waves of a basic reactive shock wave model are proved to be nonlinearly stable, i.e. initially perturbed waves tend asymptotically to translated weak detonation waves. This model system was derived as the low Mach number limit of the one component reactive Navier-Stokes equations by Majda and Roytburd [SIAM J. Sci. Stat. Comput. 43, 1086–1118 (1983)], and its weak detonation waves have been numerically observed as stable. The analysis shows in particular the key role of the new nonlinear dynamics of the position of the shock wave, The shock translation solves a nonlinear integral equation, obtained by Green's function techniques, and its solution is estimated by observing that the kernel can be split into a dominating convolution operator and a remainder. The inverse operator of the convolution and detailed properties of the traveling wave reduce, by monotonicity, the remainder to a small L ₁ perturbation. Received: 17 August 1998 / Accepted: 13 November 1998     

Numerical study of shock wave entry and propagation in a microchannel

The entry of a shock wave with the Mach number M_is = 2.03 into a microchannel and its further propagation is numerically studied with the use of kinetic and continuum approaches. Numerical simulations on the basis of the Navier ?? Stokes equations and the Direct Simulation Monte Carlo method are performed for different Knudsen numbers Kn = 8·10^?3 and 8·10^?2 based on the microchannel half-height. At the Knudsen number Kn = 8·10^?3, amplification of the shock wave after its entry into the microchannel is observed. Further downstream, the shock wave is attenuated, which is in qualitative agreement with experimental data. It is demonstrated that results predicted by a quasi-one-dimensional model (which ignores viscosity and heat conduction) of shock wave propagation over a channel with an abrupt change in the area agrees with results of numerical simulations on the basis of the Euler equations. In both cases, shock wave acceleration (amplification) after its entry into the microchannel is observed. At the Knudsen number Kn = 8·10^?2, the influence of the entrance shape on shock wave propagation over the microchannel is examined. Intense attenuation of the shock wave is observed in three cases: channel with sudden contraction, junction of two channels with an additional thin separating plate, and rounded junction in the form of a sector with an angle of 90° (quarter of a circumference). It is shown that the microchannel entrance shape can affect further propagation of the shock wave. The wave has the highest velocity in the case with a rounded entrance.     

Molecular Dynamic Simulation of High Speed Deformation Effect on Microstructure Change in Thermal Heated Metals

In the present paper computer simulation of high-speed deformation (shock wave propagation) by molecular dynamic method is performed in thin copper sample, having the form of rectangular parallelepiped (10 a ‐ 10 a ‐ 20 a , where a is the lattice constant) with 8000 atoms. On the surfaces Z 0 =0 and Z max =20 a the mirror boundary conditions with rigid walls and the periodic boundary conditions along X and Y directions corresponding to short sides of deformed crystal are used, which allows to investigate the reflection of shock wave from the surfaces in Z direction. The changes of microstructure have been investigated up to 12 ps. The numerical calculations of microstructure changes have been performed here taking into account the effect of thermal heating of crystal lattice before shock wave front. The numerical results show that comparing with the propagation of shock waves under room temperature in thermal heated structure additional displaced atoms (vacancies and interstitials) are produced. The obtained results show that the production of point defects during high-speed deformation is determined by the thermal softening of microstructure and generation rate of point defects very strong increases with an increasing of high speed deformation rate. The peculiarities of microstructure changes in deformed copper are analyzed here at the different initial temperatures and various high-speed deformations (average ion velocities behind shock wave).     

Problems of interpretation of holographic interferograms near shock wave fronts

A method of avoiding ambiguity in the interpretation of interferograms near a shock wave front is proposed. The method is based on combining the double-exposure schlieren method and holographic interferometry. Relations for calculating, on the basis of data obtained by analyzing double-exposure schlieren photographs, both the density at the shock wave front and the gradient of the density directly behind the front, which is necessary for calculating the shifts of the interference fringes near the shock wave front, are presented. Zh. Tekh. Fiz. 68, 88–91 (September 1998)     

激波统一来解实用公式

本文从麦克斯韦方程组和三个守恒定律出发,导出磁流体力学激波间断面两边所满足的跃变方程组,通过无量纲化处理,即可得到激波统一求解公式及其解析解。这些公式对理论研究和工程应用中的实际计算都很方便。最后,对垂直磁流体激波、平行磁流体激波以及气体动力学激波分别予以讨论。     

Numerical Study of Air Flow Induced by Shock Impact on an Array of Perforated Plates

This study is focused on the propagation behavior and attenuation characteristics of a planar incident shock wave when propagating through an array of perforated plates. Based on a density-based coupled explicit algorithm, combined with a third-order MUSCL scheme and the Roe averaged flux difference splitting method, the Navier–Stokes equations and the realizable k-ε turbulence model equations describing the air flow are numerically solved. The evolution of the dynamic wave and ring vortex systems is effectively captured and analyzed. The influence of incident shock Mach number, perforated-plate porosity, and plate number on the propagation and attenuation of the shock wave was studied by using pressure- and entropy-based attenuation rates. The results indicate that the reflection, diffraction, transmission, and interference behaviors of the leading shock wave and the superimposed effects due to the trailing secondary shock wave are the main reasons that cause the intensity of the leading shock wave to experience a complex process consisting of attenuation, local enhancement, attenuation, enhancement, and attenuation. The reflected shock interactions with transmitted shock induced ring vortices and jets lead to the deformation and local intensification of the shock wave. The formation of nearly steady jets following the array of perforated plates is attributed to the generation of an oscillation chamber for the inside dynamic wave system between two perforated plates. The vorticity diffusion, merging and splitting of vortex cores dissipate the wave energy. Furthermore, the leading transmitted shock wave attenuates more significantly whereas the reflected shock wave from the first plate of the array attenuates less significantly as the shock Mach number increases. The increase in the porosity weakens the suppression effects on the leading shock wave while increases the attenuation rate of the reflected shock wave. The first perforated plate in the array plays a major role in the attenuation of the shock wave.     

A boundary integral method to compute Green's functions for the radiative transport equation

The Green's function for the time-independent radiative transport equation in the whole space can be computed as an expansion in plane wave solutions. Plane wave solutions are a general class of solutions for the radiative transport equation. Because plane wave solutions are not known analytically in general, we calculate them numerically using the discrete ordinate method. We use the whole space Green's function to derive boundary integral equations. Through the solution of the boundary integral equations, we compute the Green's function for bounded domains. In particular we compute the Green's function for the half space, the slab, and the two-layered half space. The boundary conditions used here are in their most general form. Hence, this theory can be applied to boundaries with any kind of reflection and transmission law.     

Convective heat conduction and diffusion in one-dimensional hydrodynamics

The impurity concentration in localized structures is described on the basis of analytic solutions of model equations for convective diffusion in the one-dimensional hydrodynamic approximation without pressure. The simplicity of the derivation of the analytic results depends on the ratio of the kinetic coefficients of the liquid (the Prandtl numbers). For the same kinetic coefficients, any time-dependent problem can be reduced to problems for the conventional heat conduction equation. For integer Prandtl numbers the problem of time-dependent convective diffusion in the flow field of a uniformly moving shock wave likewise reduces to problems for the heat conduction equation. Relations are established between problems whose Prandtl numbers differ by an integer. Various representations of the Green’s functions for the equations of convective diffusion are analyzed. For integer Prandtl numbers they can be expressed in terms of error functions. The asymptotic character of the solutions depends strongly on the satisfaction of global conservation laws. For global conservation of the impurity mass, coalescence of shock waves corresponds to merging of impurity solitons, i.e., clustering. Zh. éksp. Teor. Fiz. 116, 1616–1629 (November 1999)     

Modelling detonation of heterogeneous explosives with embedded inert particles using detonation shock dynamics: Normal and divergent propagation in regular and simplified microstructure

This paper discusses the mathematical formulation of Detonation Shock Dynamics (DSD) regarding a detonation shock wave passing over a series of inert spherical particles embedded in a high-explosive material. DSD provides an efficient method for studying detonation front propagation in such materials without the necessity of simulating the combustion equations for the entire system. We derive a series of partial differential equations in a cylindrical coordinate system and a moving shock-attached coordinate system which describes the propagation of detonation about a single particle, where the detonation obeys a linear shock normal velocity-curvature (D_n–κ) DSD relation. We solve these equations numerically and observe the short-term and long-term behaviour of the detonation shock wave as it passes over the particles. We discuss the shape of the perturbed shock wave and demonstrate the periodic and convergent behaviour obtained when detonation passes over a regular, periodic array of inert spherical particles.     

The Self-Similarity Index of the Convergence of Strong Cylindrical Shock Waves in a Gas with a Uniform Density

A model of the convergence of cylindrical shock waves (SWs) in a gas with a uniform density has been considered. The partial differential equations of this model have been reduced to ordinary differential equations, from which the law of convergence of such shock waves and the dependence α = f(γ, γ_eff) of their self-similarity index α on the heat-capacity ratio in front of the shock wave (γ) and behind the shock wave front (γ_eff) of the gas have been found. This dependence for cylindrical shock waves has been shown to agree with the experimental data within the measurement error.     

Asymmetric interaction of blast waves and shock waves with a body flying at supersonic velocity

The problems of asymmetric interaction of a blunt wedge traveling at supersonic velocity with a cylindrical blast wave from a point explosion and with a plane shock wave are investigated by numerical simulation. The evolution of the interaction flow is analyzed, and data are obtained on how the structure of the shock layer changes. Zh. Tekh. Fiz. 69, 15–19 (May 1999)     

Superspherical cumulation. Converging shock waves with amplitudes increasing faster than in spherical cumulation

Experimental, numerical, and theoretical investigations are made of a gas flow generated by a pulsed high-current discharge in an axisymmetric cavity bounded by a spherical lens adjacent to a flat plate. It is shown that the shock wave forming in the discharge and converging toward the axis is accelerated and amplified as it converges. The amplitude of the shock wave increases faster than does that of a spherical converging shock wave. Zh. Tekh. Fiz. 69, 10–18 (March 1999)     

Dynamical coupling in the differential equations approach to atom-diatom exchange reactions

The atomic exchange reaction A + BC → AB + C is investigated quantum mechanically employing a coupled differential equations approach. The relative motion in reactant and product channels is described in the common coordinate R ₃ (the AC nuclear separation) and is developed in three-dimensional space. The total wave functions of the system are expressed as a superposition of valence bond electronic states of the initial (A, BC) and final (AB, C) configurations, with the coefficients describing the relative and internal (vibrational, rotational) nuclear motions. Choosing convenient trial functions with the appropriate boundary conditions and using the Kohn variational principle, a set of differential (rather than the usual integro-differential) equations is obtained for the relative motion wave functions in R ₃. The potential matrix elements turn out to be dynamical in that they depend on the initial k ₁ and final k ₂ wave vectors. Two-state coupled channel calculations of the differential and integral cross sections for the isotopic species D + H₂, H + H₂ and D + D₂ are presented for collision energies up to 0·8 eV.     

Going to 10 TPa: The Calculated Hugoniots for Cu,Ta, and Mo

The Hugoniot equations of state for shock compressed Cu, Ta, and Mo are calculated at pressures up to 4 TPa and then up to 10 TPa are obtained by extrapolation. The calculations are parameter-free in that the cold part of the Helmholtz free-energy is calculated using the first-principles full-potential linearized augmented plane wave method within the generalized gradient approximation, the thermal contribution to the Helmholtz free-energy due to the lattice oscillations is calculated using the recently developed classical mean-field potential approach, and that due to the thermal electrons is calculated using the one-dimensional numerical integration. The calculated results agree with the existing experimental values very well.     

Optical study of parameters of the passage of a shock wave from air to water

A study is made of the penetration of shock waves from air into water. The shock wave in air is generated as a result of dielectric breakdown induced by pulsed CO₂-laser radiation. A combination of the double-exposure shadow method and holographic interferometry is used to measure the shock-wave parameters. Density and pressure profiles behind the wave front are obtained at different times after onset of breakdown. It is shown experimentally that as the wave passes through the interface from the air to the water, there is a fourfold amplification of the pressure in the shock wave front. Estimates of the width of the shock wave front formed in the water are given in the context of studies of large-scale explosion processes. It is shown that simple empirical dependences, established in the course of studies of large-scale explosions, are also valid with certain corrections for microscopic laboratory experiments. Zh. Tekh. Fiz. 68, 39–43 (August 1998)     

平面定常超音速绕流问题的新的数值解法

用特征线法解平面定常超音速绕流问题虽然有效,但当激波很弱、几乎与特征线平行时则很难处理。用有限差分法计算此问题也比较复杂。本文把作者在文章中^[1]提出的新的数值解法,发展并应用到平面定常超音速绕流的问题。仍然采用了许为厚教授提出的新拉格朗日变量^[2],这使边界条件的提法大为简化。此新的数值解法按变量指标之和,一排排地往下计算,方法简单,可以处理各种形状物体的超音速绕流。本文对向上弯曲的抛物形固壁绕流向题的实例进行了具体计算。算出了激波的形状。当激波没有形成以前,相应的普朗特一迈耶气流是有准确分析解的,把数值解与准确解进行了l比较结果是满意的。     

State equation of condensed matter at high pressure: <Emphasis Type="Italic">D-U</Emphasis> diagram approach

For solids and liquids, an equation of state is suggested at high pressures up to a few megabars, for densities greater than that at normal conditions and for temperatures up to the melting point. Shock wave loading test data are analyzed for 40 basic chemical elements, and they prove the state equation suggested, within the limits of test error. The method is based on the analysis of D-U diagrams where D is the shock wave velocity and U is the material velocity behind the shock wave (both with respect to the material in front of the shock wave). Based on the state equation suggested the velocity of shock wave is shown to be a linear function of the material velocity behind the shock wave, the function being a specific characteristic of the material and its structure. Most significant anomaly belonging to carbon, iron, ice, and water is explained by the formation of new phases at high pressure, with two new phases of iron, and one phase in the case of water. For water, a simple nearly exact equation of state is suggested for pressures from 0.1 MPa to 150 GPa. For pressures from 0.1 to 300 MPa, it fits very well the extremely complicated state equation of the American standard obtained by static tests, and for pressures from 2 to 50 GPa it fits well the data of shock wave tests. In the pressure range from 45 to 1500 GPa liquid water becomes solid, which equation of state coincides with that of alkaline metal sodium. The model of ideal solid as contrary to ideal gas is introduced, with internal energy of ideal solid depending only on stresses or strains (and only on pressure or density, at high pressures). The equations of state for iron, diamond, pyrolithic graphite, and for several phases of ice are as well derived based on test data.     

Planarity and stability of shock driven directly by the ninth laser beam from “Shenguang-II" laser facility

Using the ninth laser beam (converted to 2ω) of “Shenguang-II” laser facility and the beam smoothing technology of lens-array [Appl. Opt. 25, 377 (1986); Phys. Plasmas. 9, 3201 (1995)], a shock wave with 700 μm (the root-mean-square of shock breakout time (RMS) RMS ≈ 6.32 ps) flat top was created. An Al-Al four-step target was designed to do research on shock wave stability in an Al target. And the shock stability experiment with the Al-Al four-step target indicated that the shock wave steadily propagated in the Al target of thickness of about 20–45 μm under the power density of ~ 1.0×10¹⁴ W/cm².