APPLICATION OF DETONATION SHOCK DYNAMICS (DSD) TO YOUNGS-TYPE DISCONTINUOUS INTERFACE GEOMETRY

Detonation Shock Dynamics (DSD) describes the evolution of a two‐ or three‐dimensional detonation wave in a way that accounts for the finite detonation reaction zone thickness as long as it is sufficiently small compared to the radius of curvature of the detonation wave front. The current DSD solver obtains its input parameters by superimposing a rectangular Cartesian grid over the high explosive (HE) regions, determining the signed normal distance from each grid point to the nearest point on the HE boundary (negative on the inside of HE and positive on the outside) and assigning a material identification to each grid point. It has been shown previously to work well with a Lagrangian geometric description where the mesh entities, particularly cell faces, are contiguous and therefore distances to the HE boundary are precisely defined. In this paper a new scheme for the DSD driver code is presented that allows the HE boundary to be represented in a noncontiguous fashion, such as is obtained from a Youngs‐type material interface reconstruction, as often used in Eulerian hydrodynamics codes.