No abstract

No abstract

Geometrical Shock Dynamics (GSD) is a simplified model for nonlinear shock wave propagation. It is based on the decomposition of the shock front into elementary ray tubes with a simple relation linking its local curvature and velocity. This relation is obtained under the assumption of strong shock in order to neglect the effect of the post-shock flow on the front. More recently, a new simplified model, referenced as the Kinematic model, was proposed. This model is obtained by combining the threedimensional Euler equations and the Rankine-Hugoniot relations at the front, which leads to an equation for the normal variation of the shock Mach number at the wave front. In the same way as GSD, the Kinematic model is closed by neglecting the post-shock flow effects. Although each model's approach is different, we prove here their structural equivalence: the Kinematic model can be rewritten under the form of GSD with a specific A − M relation. Both models are thus compared through a wide variety of examples including experimental data or Eulerian simulations results when available. Attention is drawn to the simple cases of compression ramps and convex corners' diffraction. The analysis is completed by the more complex cases of the diffraction over a cylinder, a sphere, a mound and a trough.

We develop a new algorithm for the computation of the geometrical shock dynamics model (GSD). The method relies on the fast-marching paradigm and enables the discrete evaluation of the first arrival time of a shock wave and its local velocity on a cartesian grid. The proposed algorithm is based on a second order upwind finite difference scheme and reduces to a local nonlinear system of two equations solved by an iterative procedure. Reference solutions are built for a smooth radial configuration and for the 2D Riemann problem. The link between the GSD model and p-systems is given. Numerical experiments demonstrate the accuracy and the ability of the scheme to handle singularities.