Oblique shock reflection from an axis of symmetry: shock dynamics and relation to the Guderley singularity

Oblique shock reflection from an axis of symmetry is studied using Whitham's theory of geometrical shock dynamics, and the results are compared with previous numerical simulations of the phenomenon by Hornung (2000). The shock shapes (for strong and weak shocks), and the location of the shock-shock (for strong shocks), are in good agreement with the numerical results, though the detail of the shock reflection structure is, of course, not resolved by shock dynamics. A guess at a mathematical form of the shock shape based on an analogy with the Guderley singularity in cylindrical shock implosion, in the form of a generalized hyperbola, fits the shock shape very well. The smooth variation of the exponent in this equation with initial shock angle from the Guderley value at zero to 0.5 at 90° supports the analogy. Finally, steady-flow shock reflection from a symmetry axis is related to the self-similar flow

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

Moving Overlapping Grids with Adaptive Mesh Refinement for High-Speed Reactive and Non-reactive Flow

We consider the solution of the reactive and non-reactive Euler equations on two-dimensional domains that evolve in time. The domains are discretized using moving overlapping grids. In a typical grid construction, boundary-fitted grids are used to represent moving boundaries, and these grids overlap with stationary background Cartesian grids. Block-structured adaptive mesh refinement (AMR) is used to resolve fine-scale features in the flow such as shocks and detonations. Refinement grids are added to base-level grids according to an estimate of the error, and these refinement grids move with their corresponding base-level grids. The numerical approximation of the governing equations takes place in the parameter space of each component grid which is defined by a mapping from (fixed) parameter space to (moving) physical space. The mapped equations are solved numerically using a second-order extension of Godunov's method. The stiff source term in the reactive case is handled using a Runge-Kutta error-control scheme. We consider cases when the boundaries move according to a prescribed function of time and when the boundaries of embedded bodies move according to the surface stress exerted by the fluid. In the latter case, the Newton-Euler equations describe the motion of the center of mass of the each body and the rotation about it, and these equations are integrated numerically using a second-order predictor-corrector scheme. Numerical boundary conditions at slip walls are described, and numerical results are presented for both reactive and non-reactive flows in order to demonstrate the use and accuracy of the numerical approach

Parallel Computation of Three-Dimensional Flows using Overlapping Grids with Adaptive Mesh Refinement

This paper describes an approach for the numerical solution of time-dependent partial differential equations in complex three-dimensional domains. The domains are represented by overlapping structured grids, and block-structured adaptive mesh refinement (AMR) is employed to locally increase the grid resolution. In addition, the numerical method is implemented on parallel distributed-memory computers using a domain-decomposition approach. The implementation is flexible so that each base grid within the overlapping grid structure and its associated refinement grids can be independently partitioned over a chosen set of processors. A modified bin-packing algorithm is used to specify the partition for each grid so that the computational work is evenly distributed amongst the processors. All components of the AMR algorithm such as error estimation, regridding, and interpolation are performed in parallel. The parallel time-stepping algorithm is illustrated for initial-boundary-value problems involving a linear advection-diffusion equation and the (nonlinear) reactive Euler equations. Numerical results are presented for both equations to demonstrate the accuracy and correctness of the parallel approach. Exact solutions of the advection-diffusion equation are constructed, and these are used to check the corresponding numerical solutions for a variety of tests involving different overlapping grids, different numbers of refinement levels and refinement ratios, and different numbers of processors. The problem of planar shock diffraction by a sphere is considered as an illustration of the numerical approach for the Euler equations, and a problem involving the initiation of a detonation from a hot spot in a T-shaped pipe is considered to demonstrate the numerical approach for the reactive case. For both problems, the solutions are shown to be well resolved on the finest grid. The parallel performance of the approach is examined in detail for the shock diffraction problem

A High-Resolution Godunov Method for Compressible Multi-Material Flow on Overlapping Grids

A numerical method is described for inviscid, compressible, multi-material flow in two space dimensions. The flow is governed by the multi-material Euler equations with a general mixture equation of state. Composite overlapping grids are used to handle complex flow geometry and block-structured adaptive mesh refinement (AMR) is used to locally increase grid resolution near shocks and material interfaces. The discretization of the governing equations is based on a high-resolution Godunov method, but includes an energy correction designed to suppress numerical errors that develop near a material interface for standard, conservative shock-capturing schemes. The energy correction is constructed based on a uniform pressure-velocity flow and is significant only near the captured interface. A variety of two-material flows are presented to verify the accuracy of the numerical approach and to illustrate its use. These flows assume an equation of state for the mixture based on Jones-Wilkins-Lee (JWL) forms for the components. This equation of state includes a mixture of ideal gases as a special case. Flow problems considered include unsteady one-dimensional shock-interface collision, steady interaction of an planar interface and an oblique shock, planar shock interaction with a collection of gas-filled cylindrical inhomogeneities, and the impulsive motion of the two-component mixture in a rigid cylindrical vessel

A Study of Detonation Propagation and Diffraction with Compliant Confinement

A previous computational study of diffracting detonations with the ignition-and-growth model demonstrated that contrary to experimental observations, the computed solution did not exhibit dead zones. For a rigidly confined explosive it was found that while diffraction past a sharp corner did lead to a temporary separation of the lead shock from the reaction zone, the detonation re-established itself in due course and no pockets of unreacted material were left behind. The present investigation continues to focus on the potential for detonation failure within the ignition-and-growth (IG) model, but now for a compliant confinement of the explosive. The aim of the present paper is two fold. First, in order to compute solutions of the governing equations for multi-material reactive flow, a numerical method of solution is developed and discussed. The method is a Godunov-type, fractional-step scheme which incorporates an energy correction to suppress numerical oscillations that would occur near the material interface separating the reactive material and the inert confiner for standard conservative schemes. The numerical method uses adaptive mesh refinement (AMR) on overlapping grids, and the accuracy of solutions is well tested using a two-dimensional rate-stick problem for both strong and weak inert confinements. The second aim of the paper is to extend the previous computational study of the IG model by considering two related problems. In the first problem, the corner-turning configuration is re-examined, and it is shown that in the matter of detonation failure, the absence of rigid confinement does not affect the outcome in a material way; sustained dead zones continue to elude the model. In the second problem, detonations propagating down a compliantly confined pencil-shaped configuration are computed for a variety of cone angles of the tapered section. It is found, in accord with experimental observation, that if the cone angle is small enough, the detonation fails prior to reaching the cone tip. For both the corner-turning and the pencil-shaped configurations, mechanisms underlying the behavior of the computed solutions are identified. It is concluded that disagreement between computation and experiment in the corner-turning case lies in the absence, in the model, of a mechanism that allows the explosive to undergo desensitization when subjected to a weak shock