On the existence of intermediate magnetohydrodynamic shock waves

We consider the questions related to the structure of shock waves for a system of magnetohydrodynamic equations. Using Conley's connection matrix, we recover and extend earlier results due to C. Conley and J. Smoller. In particular, we give a simpler proof of the existence of fast and slow shocks with structure. We also demonstrate that for some viscosity parameters intermediate shocks occur. Furthermore, under an assumption of transversality, we show that there exist multi-parameter families of these intermediate shocks.This research was done while both authors were visiting the Lefschetz Center for Dynamical Systems at Brown University.Supported in part by the NSF under Grant DMS-8507056.Supported in part by AFOSR 87-0347.     

Nonlinear wave propagation in binary mixtures of Euler fluids

In the present paper, we study the propagation of acceleration and shock waves in a binary mixture of ideal Euler fluids, assuming that the difference between the atomic masses of the constituents is negligible. We evaluate the characteristic speeds, proving that they can be separated into two groups: one is related to the case of a single Euler fluid, provided that an average ratio of specific heats is introduced; the other is new and related to the propagation speed due to diffusion. We evaluate the critical time for sound acceleration waves and compare its value to that of a single fluid. We then study shock waves, showing that three types of shock waves appear: sonic and contact shocks, which have counterparts in the single fluid case, and the diffusive shock, which is peculiar to the mixture. We discuss the admissibility of the shock waves using the Lax-Liu conditions and the entropy growth criterion. It is proved that the sonic and the characteristic shock obey the same properties as in the single fluid case, while for the diffusive shock there exists a locally exceptional case that is determined by a particular value of the concentration of the constituents, for which the genuine nonlinearity is lost and no shocks are admissible. For other values of the unperturbed concentration, the diffusive shock is stable in a bounded interval of admissibility.Received: 15 December 2002, Accepted: 28 June 2003 Correspondence to: T. RuggeriS. Simi: On leave from the Department of Mechanics, Faculty of Technical Sciences, University of Novi Sad, Serbia     

Wave interaction in a nonequilibrium gas flow

By employing the method of multiple time scales, we derive here the transport equations for the primary amplitudes of resonantly interacting high-frequency waves propagating into a non-equilibrium gas flow. Evolutionary behavior of non-resonant wave modes culminating into shocks or no shocks, together with their asymptotic decay behavior, is studied. Effects of non-linearity, which are noticeable over times of order O(ε^-1), are examined, and the model evolution equations for resonantly interacting multi-wave modes are derived.     

On the response of rubbers at high strain rates—III. Effect of hysteresis

In this series of papers, we examine the propagation of waves of finite deformation in rubbers through experiments and analysis; in the present paper, Part III, the effect of hysteretic material behavior on the free retraction of prestretched rubber is considered. A rubber strip stretched to many times its initial length is released at one end and the resulting unloading is examined. A high-speed video camera was used to monitor the motion and to determine the evolution of strain and particle velocity in rubber strips. Simple waves as well as shock waves are observed in these unloading experiments. The measurements are modeled using a power-law model of the material behavior. The hysteretic material response and the formation of shocks are characterized.     

Spallation of polycarbonate under laser-driven shocks

Laser driven shocks have been used to investigate dynamic failure (spallation) of polycarbonate under uniaxial tensile loading at very high strain rate, of the order of 10 s. First, uninstrumented recovery shots have been performed, post-test examination of the fracture damage has been carried out, and the influences of the experimental parameters (loading conditions and target thickness) have been analyzed. Then, an attempt to model the response of polycarbonate to plane shock loading has been made. On one hand, in-situ measurements have been performed in polycarbonate samples submitted to the plane detonation wave of a strong explosive, and the results have led to content with simple constitutive relations. On the other hand, piezoelectric measurements under laser shocks have provided a characterization of the loading pressure pulse, and comparisons of the measured and computed signals have confirmed the ability of the model to describe wave propagation in polycarbonate. Finally, the spallation experiments have been simulated. A spall strength has been estimated, on the basis of the experimental data, and the predictive capability of the model has been tested. Received: 18 February 1997 / Accepted 1 April 1997     

Numerical simulations of shock wave reflection phenomena in non-stationary flows using regularized smoothed particle hydrodynamics (RSPH)

In this paper we wish to demonstrate to what extent the numerical method regularized smoothed particle hydrodynamics (RSPH) is capable of modelling shocks and shock reflection patterns in a satisfactory manner. The use of SPH based methods to model shock wave problems has been relatively sparse, both due to historical reasons, as the method was originally developed for studies of astrophysical gas dynamics, but also due to the fact that boundary treatment in Lagrangian methods may be a difficult task. The boundary conditions have therefore been given special attention in this paper. Results presented for one quasi-stationary and three non-stationary flow tests reveal a high degree of similarity, when compared to published numerical and experimental data. The difference is found to be below 5, in the case where experimental data was found tabulated. The transition from regular reflection (RR) to Mach reflection (MR) and the opposite transition from MR to RR are studied. The results are found to be in close agreement with the results obtained from various empirical and semi-empirical formulas published in the literature. A convergence test shows a convergence rate slightly steeper than linear, comparable to what is found for other numerical methods when shocks are involved.     

Nonlinear waves in electromigration dispersion in a capillary

We construct exact solutions to an unusual nonlinear advection–diffusion equation arising in the study of Taylor–Aris (also known as shear) dispersion due to electroosmotic flow during electromigration in a capillary. An exact reduction to a Darboux equation is found under a traveling-wave ansatz. The equilibria of this ordinary differential equation are analyzed, showing that their stability is determined solely by the (dimensionless) wave speed without regard to any (dimensionless) physical parameters. Integral curves, connecting the appropriate equilibria of the Darboux equation that governs traveling waves, are constructed, which in turn are shown to be asymmetric kink solutions (i.e., non-Taylor shocks). Furthermore, it is shown that the governing Darboux equation exhibits bistability, which leads to two coexisting non-negative kink solutions for (dimensionless) wave speeds greater than unity. Finally, we give some remarks on other types of traveling-wave solutions and a discussion of some approximations of the governing partial differential equation of electromigration dispersion.     

Suppression of pull-in in a microstructure actuated by mechanical shocks and electrostatic forces

This work investigates the effect of a high-frequency voltage (HFV) on the pull-in instability in a microstructure actuated by mechanical shocks and electrostatic forces. The microstructure is modelled as a single-degree-of-freedom mass-spring-damper system. The method of direct partition of motion is used to split the fast and slow dynamics. Analysis of steady-state solutions of the slow dynamic allows the investigation of the influence of the HFV on the pull-in. The results show that adding HFV rigidifies the system, creates new stable equilibria and suppresses the pull-in instability for adequate high-frequency voltages. To illustrate the applicability of the result, a specific capacitive microelectromechanical system consisting of a clamped-clamped microbeam is considered.     

Numerical schemes with high order of accuracy for the computation of shock waves

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｏｍｐｕｔａｔｉｏｎｏｆｆｌｏｗｆｉｅｌｄｗｉｔｈｓｈｏｃｋｗａｖｅｓｈａｓｂｅｅｎｔｈｅｓｕｂｊｅｃｔｏｆｒｅｓｅａｒｃｈｆｏｒｍａｎｙｙｅａｒｓ.Ｔｈｅｒｅａｒｅｂａｓｉｃａｌｌｙｔｗｏｃａｔｅｇｏｒｉｅｓｏｆｍｅｔｈｏｄｓ,ｎａｍｅｌｙ,ｓｈｏｃｋｆｉｔｔｉｎｇｍｅｔｈｏｄａｎｄｓｈｏｃｋｃａｐｔｕｒｉｎｇｍｅｔｈｏｄ.Ｔｈｅｆｏｒｍｅｒｄｉｖｉｄｅｓｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｄｏｍａｉｎｉｎｔｏｓｕｂ＿ｄｏｍａｉｎｓｂｙｔｈｅｓｈｏｃｋｗａｖｅｓ.Ｉｎｅ…     

Constant rate production of geothermal fluid from a two-phase vertical column I: Theory

Previous theoretical results on geothermal two-phase flows in porous media are extended and applied to the case of withdrawal of fluid at a constant rate from a vertical column. Dimensional considerations show that pressure and saturation behaviour is controlled by a single parameter which is the ratio of the withdrawal speed to buoyancy speed. For large flows (large) fluid withdrawal is a mining process, and a vapour dominated zone spreads out from the production level. Production enthalpies tend towards steam values. However if is small then gravity dominates, and buoyancy forces can lead to the formation of a steam bubble which escapes from the production boundary and rises towards the surface. Production enthalpy may then remain at the liquid value over long periods. In addition certain saturation ranges at the sink may be forbidden as a consequence of the constant rate boundary condition. Then saturation shocks will form at the production boundary and travel out from the sink. Internally generated shocks may also occur. Pressure and saturation response to a steady withdrawal of fluid is more complicated than in a two-phase gravity-free situation. Since gravity is an essential component of even horizontal two-phase flow this suggests that two-phase studies which ignore the role of gravity may be too simplistic.     

On the response of rubbers at high strain rates—II. Shock waves

In this series of papers, we examine the propagation of waves of finite deformation in rubbers through experiments and analysis; in the present paper, Part II, attention is focused on the propagation of one-dimensional tensile shock waves in strips of latex and nitrile rubber. Tensile wave propagation experiments were conducted at high strain rates by holding one end fixed and displacing the other end at a constant velocity. A high-speed video camera was used to monitor the motion and to determine the evolution of strain and particle velocity in rubber strips. Shock waves have been generated under tensile impact in prestretched rubber strips; analysis of the response yields the tensile shock adiabat for rubbers. The propagation of shocks is analyzed by developing an analogy with the theory of detonation; it is shown that the condition for shock propagation can be determined using the Chapman-Jouguet shock condition.     

A simple physical theory of weak Mach reflection over plane surfaces

It is shown that simple physical principles coupled with the inviscid shock jump relations can be applied to the problem of weak Mach reflection to the extent that the triple point path can be predicted from the incident shock Mach number , gas specific heat ratio and the inclination angle of the reflecting surface to the shock normal. Comparison with the Euler code data and with experiments show close agreement for conditions both far and close to transition and that the general shape of the reflected and Mach stem shocks follow simple curves except in the neighbourhood of the triple point. The conflict at the triple point in matching the flow deflection angles and pressures across the contact discontinuity remains. It is shown however that the simple model presented here gives a close match to the cfd and experimental overall shock and contact surface shapes although it cannot predict these or the flow properties in any detail. Received 10 May 1999 / Accepted 17 December 1999     

Modeling plastic shocks in periodic laminates with gradient plasticity theories

Steady plastic shocks generated by planar impact on metal-polymer laminate composites are analyzed in the framework of gradient plasticity theories. The laminate material has a periodic structure with a unit cell composed of two layers of different materials. First- and second-order gradient plasticity theories are used to model the structure of steady plastic shocks. In both theories, the effect of the internal structure is accounted for at the macroscopic level by two material parameters that are dependent upon the layer's thickness and the properties of constituents. These two structure parameters are shown to be uniquely determined from experimental data. Theoretical predictions are compared with experiments for different cell sizes and for various shock intensities. In particular, the following experimental features are well-reproduced by the modeling:

•

the shock width is proportional to the cell size;

•

the magnitude of strain rate is inversely proportional to cell size and increases with the amplitude of applied stress following a power law.

While these results are equally described by both the plasticity theories, the first gradient plasticity approach seems to be favored when comparing the structure of the shock front to the experimental data.     

Collective Bow Shock Ahead of a Transverse System of Spheres in a Supersonic Flow Behind a Moving Shock Wave

The results of investigating the dynamics and physical conditions of formation of a collective bow shock ahead of a system of spheres with the line of centers normal to the supersonic flow behind a traveling shock wave are presented. Two types of shock-wave patterns that necessarily precede the formation of the collective shock wave and correspond to regular and Mach interaction of the bow shocks were detected experimentally. On the basis of a local gasdynamic-discontinuity interference theory, quantitative criteria of the existence of these regimes and of the formation of a common shock wave are determined. These criteria are confirmed in a series of experiments for the transitional regimes.     

关于气动声学数值计算的方法与进展

气动声学数值计算是近年才出现的研究领域。本文介绍了气动声学数值计算的方法和有关的问题、边界条件的处理以及计算非线性声波的数值方法和进展。讨论了计算气动声学(CAA)的特性及其与计算流体力学(CFD)的差异，指出气动声学数值方法的关键是建立能保持色散关系的差分方程和正确处理无反射边界条件。对于非线性声波传播的问题，为了得到正确的解，应注意提高差分格式对短波的分辨能力，同时发展能抑制“伪”振荡(短波)而对长波基本不起作用的数值方法。     

冲击波对含水炸药减敏作用的实验研究

对水胶炸药和乳化炸药在冲击波作用后的爆炸性能作了对比实验研究,结果表明水胶炸药的抗压性能明显优于乳化炸药,分析认为造成这种差别的原因是两种炸药不同的结构组分和敏化方式。用示波器捕捉两种炸药受冲击波作用前后的水中爆炸冲击波,利用冲击波参数计算炸药的减敏程度,将炸药的压力减敏作用数量化,克服了以往的研究方法只能用半爆率、拒爆率等概念计数表述的不足。     

On the distance to blow-up of acceleration waves in random media

We study the effects of material spatial randomness on the distance to form shocks from acceleration waves, , in random media. We introduce this randomness by taking the material coefficients and – that represent the dissipation and elastic nonlinearity, respectively, in the governing Bernoulli equation – as a stochastic vector process. The focus of our investigation is the resulting stochastic, rather than deterministic as in classical continuum mechanics studies, competition of dissipation and elastic nonlinearity. Quantitative results for are obtained by the method of moments in special simple cases, and otherwise by the method of maximum entropy. We find that the effect of even very weak random perturbation in and may be very significant on . In particular, the full negative cross-correlation between and $ results in the strongest scatter of , and hence, in the largest probability of shock formation in a given distance x. Received November 6, 2001 / Published online September 4, 2002 Dedicated to Professor Ingo Müller on the occasion of his 65th birthday Communicated by Kolumban Hutter, Darmstadt     

Ultrafast high repetition rate absorption spectroscopy of polymer shock compression

A dual-beam transient absorption spectrometer for high repetition rate (80 shocks per second) studies of shock compressed materials is described. The apparatus time response is 100 ps, so the time resolution of the shock compression process is generally limited by the shock transit time across the sample. In turn the sample thickness is limited by the sensitivity of the spectrometer. Using 400 nm thick samples of R640 dye aggregates in \textit{poly} methyl methacrylate (PMMA) and a 4.2 GPa laser-driven shock, transient absorption spectra show a shock induced absorption redshift occurring in 500 ps, considerably longer than the 200 ps shock front transit time (round trip) through the sample. This noninstantaneous shock compression is consistent with the $\sim 300$ ps viscoelastic response of PMMA at 4.2 GPa. Received 30 July 2001 / Accepted 13 March 2002 – Published online 17 June 2002     

Discrete element simulation of shock wave propagation in polycrystalline copper

The implementation of the characteristic of compressive plasticity into the Discrete Element Code, DM2, while maintaining its quasi-molecular scheme, is described. The code is used to simulate the shock compression of polycrystalline copper at 3.35 and 11.0 GPa. The model polycrystal has a normal distribution of grain sizes, with mean diameter 14 μm, and three distinct grain orientations are permitted with respect to the shock direction; 〈1 0 0〉, 〈1 1 0〉, and 〈1 1 1〉. Particle velocity dispersion (PVD) is present in the shock-induced flow, attaining its maximum magnitude at the plastic wave rise. PVD normalised to the average particle velocity of and are yielded for the 3.35 and 11.0 GPa shocks, respectively, and are of the same order as those seen in the experiment. Non-planar elastic and plastic wave fronts are present, the distribution in shock front position increasing with propagation distance. The rate of increase of the spread in shock front positions is found to be significantly smaller than that seen in probabilistic calculations on nickel polycrystals, and this difference is attributed, in the main, to grain interaction. Reflections at free surfaces yield a region of tension near to the target free surface. Due to the dispersive nature of the shock particle velocity and the non-planarity of the shock front, the tensile pressure is distributed. This may have implications for the spall strength, which are discussed. Simulations reveal a transient shear stress distribution behind the shock front. Such a distribution agrees with that put forward by Lipkin and Asay to explain the quasi-elastic reloading phenomenon. Simulation of reloading shocks show that the shear stress distribution can give rise to quasi-elastic reloading on the grain scale.