On the direct construction of the steady traveling solution to high explosive sandwich, cylinder and aquarium tests via a streamline finite volume approximation

Abstract The cylinder, sandwich and aquarium tests are common experiments used to investigate the performance of high explosives. This paper presents a solution technique which simulates these experiments as two dimensional problems in steady traveling coordinates. Using a coordinate system aligned with the streamlines of the flow, an efficient numerical scheme is developed which uses an approximate Riemann solver to model the change in state across streamlines. The technique can also accommodate arbitrary equations of state for the materials. Two verification exercises show that the algorithm converges at second order. Validation studies using the aquarium experiment agreed well with simulations, though the omission of a material strength model leads to some discrepancies between simulation and experiments for the cylinder test. This computational technique is able to predict the behavior of these common high explosive experiments on the order of minutes on a single processor with very fine spatial resolution, on the order of 10 μm.

[1]  J. Anderson,et al.  Modern Compressible Flow: With Historical Perspective , 1982 .

[2]  D. Scott Stewart,et al.  High-order numerical simulation and modelling of the interaction of energetic and inert materials , 2007 .

[3]  Tariq D. Aslam,et al.  iVCJ: A tool for Interactive Visualization of high explosives CJ states , 2016 .

[4]  DETONATION SHOCK DYNAMICS CALIBRATION OF PBX 9501 , 2007 .

[5]  Arash Kheradvar,et al.  Fundamental Fluid Mechanics , 2012 .

[6]  Tariq D. Aslam,et al.  High-order shock-fitted detonation propagation in high explosives , 2017, J. Comput. Phys..

[7]  Björn Engquist,et al.  Fast sweeping methods for hyperbolic systems of conservation laws at steady state , 2013, J. Comput. Phys..

[8]  William J. Rider,et al.  On sub-linear convergence for linearly degenerate waves in capturing schemes , 2008, J. Comput. Phys..

[9]  R. LeVeque Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[10]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

[11]  高山 和喜,et al.  Shock wave science and technology reference library , 2007 .

[12]  H. Glaz,et al.  A high-order Godunov scheme for steady supersonic gas dynamics , 1985 .

[13]  Jaw-Yen Yang,et al.  A high-order streamline Godunov scheme for steady supersonic/hypersonic equilibrium flows , 1998 .

[14]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[15]  P. Souers,et al.  Upgraded Analytical Model of the Cylinder Test , 2013 .

[16]  Daniel C. Reda Boundary-Layer Transition Experiments on Sharp, Slender Cones in Supersonic Free Flight , 1979 .

[17]  Geoffrey Ingram Taylor,et al.  The Air Pressure on a Cone Moving at High Speeds. II , 1933 .

[18]  F. Murnaghan The Compressibility of Media under Extreme Pressures. , 1944, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Kihong Kim,et al.  A particle level-set based Eulerian method for multi-material detonation simulation of high explosive and metal confinements , 2013 .

[20]  R. LeVeque Numerical methods for conservation laws , 1990 .

[21]  A LINEARIZED RIEMANN SOLVER FOR THE STEADY SUPERSONIC EULER EQUATIONS , 1993 .

[22]  S. Marsh Lasl Shock Hugoniot Data , 1980 .

[23]  R. Fedkiw,et al.  Coupling an Eulerian fluid calculation to a Lagrangian solid calculation with the ghost fluid method , 2002 .

[24]  M. Wilkins Calculation of Elastic-Plastic Flow , 1963 .

[25]  W. Fickett,et al.  Numerical Calculation of the Cylinder Test , 1975 .

[26]  Development of the Lanl Sandwich Test , 2002 .

[27]  Jaw-Yen Yang,et al.  A high-order streamline Godunov scheme for steady supersonic flow computation , 1995 .