Home   >   Products   >    Rocpack: Particle Packing   >   RocSDT – Shock to Detonation Transition

Multiphysics Simulations Supersonic nozzle flow against a flat plate, showing density (grey) and vorticity (colors). Read More >

Rocstar Simulation Cutaway of a joint slot in the Space Shuttle solid motor showing inhibitor fluid-structure interaction. Read More >

RocSDT – Shock to Detonation Transition

Modeling Shock to Detonation Transition

Illinois Rocstar has developed a mesoscale simulation suite to model shock-to-detonation in energetic materials. Possible applications include the simulation of explosive fills during penetration events, simulation of new explosive formulations for sensitivity to insult for Insensitive Munitions (IM) modeling, and simulation of novel explosive formulations to determine performance. Our approach is unique in that we:

  1. Use our micro-scale simulations codes to generate pore collapse data for explosives of interest, resulting in physics-based hot-spot generation model
  2. Use our packing code, Rocpack, to generate a realistic mesoscale microstructure
  3. Use our hydrodynamic shock code, RocSDT3D, to perform simulations of shocks propagating through explosive packs


The RocSDT3D code is written  in an object-oriented C++ style. The code is fully three-dimensional and uses a Cartesian finite volume framework to solve the Euler equations. The physics models included in RocSDT3D are the following:



RocSDTX is an axisymmetric version of RocSDT3D used to simulate non-symmetric pore collapse in energetic materials. It is intended to bridge the gap between the fast-running 1-D spherical collapse code, RocSVC, and RocSDT3D. The code is used primarily to model the collapse of a single pore due to shock interaction. It contains the models outlined above in RocSDT3D as well as a unique, 3rd order, semi-implicit, Runge-Kutta time discretization model to capture the effects of numerically stiff viscous and heat conduction effects.


RocSVC was developed to examine the symmetric collapse of a spherical void. It includes the effects of mass and heat transfer in the pore, the surrounding condensate, and at the boundary between the two. The condensate is treated as incompressible and the pressure in the pore is treated as uniform. The governing equations are cast in spherical coordinates for a reference frame that follows an evolving boundary between the pore gas and condentate. The result is a fast running 1-D code that was used to generate hundreds of simulation results that form the basis for our physics-based analytic hot-spot model.