Numerical shock propagation using geometrical shock dynamics

A simple numerical scheme for the calculation of the motion of shock waves in gases based on Whitham's theory of geometrical shock dynamics is developed. This scheme is used to study the propagation of shock waves along walls and in channels and the self-focusing of initially curved shockfronts. The numerical results are compared with exact and numerical solutions of the geometrical-shock-dynamics equations and with recent experimental investigations

Non-resonant magnetohydrodynamics streaming instability near magnetized relativistic shocks

We present in this paper both a linear study and numerical relativistic MHD simulations of the non-resonant streaming instability occurring in the precursor of relativistic shocks. In the shock front restframe, we perform a linear analysis of this instability in a likely configuration for ultra-relativistic shock precursors. This considers magneto-acoustic waves having a wave vector perpendicular to the shock front and the large scale magnetic field. Our linear analysis is achieved without any assumption on the shock velocity and is thus valid for all velocity regimes. In order to check our calculation, we also perform relativistic MHD simulations describing the propagation of the aforementioned magneto-acoustic waves through the shock precursor. The numerical calculations confirm our linear analysis, which predicts that the growth rate of the instability is maximal for ultra-relativistic shocks and exhibits a wavenumber dependence $\propto k_x^{1/2}$. Our numerical simulations also depict the saturation regime of the instability where we show that the magnetic amplification is moderate but nevertheless significant ($\delta B/B\leq 10$). This latter fact may explain the presence of strong turbulence in the vicinity of relativistic magnetized shocks. Our numerical approach also introduces a convenient means to handle isothermal (ultra-)relativistic MHD conditions.Comment: 14 pages, 6 figures, MNRAS (in press

Shocks in non-loaded bead chains with impurities

We numerically investigate the problem of the propagation of a shock in an horizontal non-loaded granular chain with a bead interaction force exponent varying from unity to large values. When $\alpha$ is close to unity we observed a cross-over between a nonlinearity-dominated regime and a solitonic one, the latest being the final steady state of the propagating wave. In the case of large values of $\alpha$ the deformation field given by the numerical simulations is completely different from the one obtained by analytical calculation. In the following we studied the interaction of these shock waves with a mass impurity placed in the bead chain. Two different physical pictures emerge whether we consider a light or a heavy impurity mass. The scatter of the shock wave with a light impurity yields damped oscillations of the impurity which then behave as a solitary wave source. Differently an heavy impurity is just shifted by the shock and the transmitted wave loses its solitonic character being fragmented into waves of decreasing amplitudes.Comment: 9 pages, 18 figures, Accepted in European Physical Journal

Aeroacoustic Flow Phenomena Accurately Captured by New Computational Fluid Dynamics Method

One of the challenges in the computational fluid dynamics area is the accurate calculation of aeroacoustic phenomena, especially in the presence of shock waves. One such phenomenon is "transonic resonance," where an unsteady shock wave at the throat of a convergent-divergent nozzle results in the emission of acoustic tones. The space-time Conservation-Element and Solution-Element (CE/SE) method developed at the NASA Glenn Research Center can faithfully capture the shock waves, their unsteady motion, and the generated acoustic tones. The CE/SE method is a revolutionary new approach to the numerical modeling of physical phenomena where features with steep gradients (e.g., shock waves, phase transition, etc.) must coexist with those having weaker variations. The CE/SE method does not require the complex interpolation procedures (that allow for the possibility of a shock between grid cells) used by many other methods to transfer information between grid cells. These interpolation procedures can add too much numerical dissipation to the solution process. Thus, while shocks are resolved, weaker waves, such as acoustic waves, are washed out

Shock waves and drag in the numerical calculation of isentropic transonic flow

Properties of the shock relations for steady, irrotational, transonic flow are discussed and compared for the full and approximate governing potential in common use. Results from numerical experiments are presented to show that the use of proper finite difference schemes provide realistic solutions and do not introduce spurious shock waves. Analysis also shows that realistic drags can be computed from shock waves that occur in isentropic flow. In analogy to the Oswatitsch drag equation, which relates the drag to entropy production in shock waves, a formula is derived for isentropic flow that relates drag to the momentum gain through an isentropic shock. A more accurate formula for drag, based on entropy production, is also derived, and examples of wave drag evaluation based on these formulas are given and comparisons are made with experimental results

Non-linear effects and shock formation in the focusing of a spherical acoustic wave : Numerical simulations and experiments in liquid helium

The focusing of acoustic waves is used to study nucleation phenomena in liquids. At large amplitude, non-linear effects are important so that the magnitude of pressure or density oscillations is difficult to predict. We present a calculation of these oscillations in a spherical geometry. We show that the main source of non-linearities is the shape of the equation of state of the liquid, enhanced by the spherical geometry. We also show that the formation of shocks cannot be ignored beyond a certain oscillation amplitude. The shock length is estimated by an analytic calculation based on the characteristics method. In our numerical simulations, we have treated the shocks with a WENO scheme. We obtain a very good agreement with experimental measurements which were recently performed in liquid helium. The comparison between numerical and experimental results allows in particular to calibrate the vibration of the ceramics used to produce the wave, as a function of the applied voltage.Comment: 20 pages, 26 figures. Submitted to The European Physical Journal

An Equation of State for Anisotropic Solids under Shock Loading

An anisotropic equation of state is proposed for accurate extrapolation of high-pressure shock Hugoniot states to other thermodynamics states for shocked single crystals and polycrystalline alloys. The proposed equation of state represents mathematical and physical generalization of the Mie-Gr\"{u}neisen equation of state for isotropic material and reduces to this equation in the limit of isotropy. Using an anisotropic nonlinear continuum framework and generalized decomposition of a stress tensor [Int. J. Plasticity \textbf{24}, 140 (2008)], the shock waves propagation along arbitrary directions in anisotropic solids of any symmetry can be examined. The non-associated strength model includes the distortion effect of the yield surface which can be used to describe the anisotropic strength differential effect. A numerical calculation showed that the general pulse shape, Hugoniot Elastic Limits (HELs), and Hugoniot stress levels for aluminum alloy 7010-T6 agree with the experimental data. The results are presented and discussed, and future studies are outlined.Comment: 6 pages, 2 figure