skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Abstract

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations.more » Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.« less

Authors:
ORCiD logo [1];  [2];  [2]
  1. Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
  2. Colorado School of Mines, Golden, CO (United States)
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1623433
Alternate Identifier(s):
OSTI ID: 1776276
Report Number(s):
LA-UR-20-22401
Journal ID: ISSN 0045-7930
Grant/Contract Number:  
89233218CNA000001
Resource Type:
Accepted Manuscript
Journal Name:
Computers and Fluids
Additional Journal Information:
Journal Volume: 205; Journal Issue: C; Journal ID: ISSN 0045-7930
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
42 ENGINEERING; detonation; reactive burn; shock capturing; Ghost Fluid Method

Citation Formats

Lozano Sanchez, Jose Eduardo, Jackson, Gregory S, and Petr, Vilem. An Eulerian multimaterial framework for simulating high-explosive aquarium tests. United States: N. p., 2020. Web. https://doi.org/10.1016/j.compfluid.2020.104524.
Lozano Sanchez, Jose Eduardo, Jackson, Gregory S, & Petr, Vilem. An Eulerian multimaterial framework for simulating high-explosive aquarium tests. United States. https://doi.org/10.1016/j.compfluid.2020.104524
Lozano Sanchez, Jose Eduardo, Jackson, Gregory S, and Petr, Vilem. Tue . "An Eulerian multimaterial framework for simulating high-explosive aquarium tests". United States. https://doi.org/10.1016/j.compfluid.2020.104524. https://www.osti.gov/servlets/purl/1623433.
@article{osti_1623433,
title = {An Eulerian multimaterial framework for simulating high-explosive aquarium tests},
author = {Lozano Sanchez, Jose Eduardo and Jackson, Gregory S and Petr, Vilem},
abstractNote = {Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.},
doi = {10.1016/j.compfluid.2020.104524},
journal = {Computers and Fluids},
number = C,
volume = 205,
place = {United States},
year = {2020},
month = {5}
}

Works referenced in this record:

The dependence of Ammonium-Nitrate Fuel-Oil (ANFO) detonation on confinement
journal, January 2017


High Explosive Detonation–Confiner Interactions
journal, January 2018


Steady Non-ideal Detonations in Cylindrical Sticks of Explosives
journal, September 2005


Performance Properties of Commercial Explosives
journal, February 1983

  • Johnson, J. N.; Mader, C. L.; Goldstein, S.
  • Propellants, Explosives, Pyrotechnics, Vol. 8, Issue 1
  • DOI: 10.1002/prep.19830080103

A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method)
journal, July 1999

  • Fedkiw, Ronald P.; Aslam, Tariq; Merriman, Barry
  • Journal of Computational Physics, Vol. 152, Issue 2
  • DOI: 10.1006/jcph.1999.6236

A Level-Set Algorithm for Tracking Discontinuities in Hyperbolic Conservation Laws
journal, March 2001


Ghost fluid method for strong shock impacting on material interface
journal, September 2003


A Real Ghost Fluid Method for the Simulation of Multimedium Compressible Flow
journal, January 2006

  • Wang, C. W.; Liu, T. G.; Khoo, B. C.
  • SIAM Journal on Scientific Computing, Vol. 28, Issue 1
  • DOI: 10.1137/030601363

Equation of State of Detonation Products Based on Exponential‐6 Potential Model and Analytical Representation of the Excess Helmholtz Free Energy
journal, March 2019

  • Suceska, Muhamed; Braithwaite, Martin; Klapötke, Thomas M.
  • Propellants, Explosives, Pyrotechnics, Vol. 44, Issue 5
  • DOI: 10.1002/prep.201800339

Stability of Chapman–Jouguet detonations for a stiffened-gas model of condensed-phase explosives
journal, March 2006

  • Short, Mark; Bdzil, John B.; Anguelova, Iana I.
  • Journal of Fluid Mechanics, Vol. 552, Issue -1
  • DOI: 10.1017/S0022112005008347

JWL++: A Simple Reactive Flow Code Package for Detonation
journal, April 2000


Numerical modelling of underwater detonation of non-ideal condensed-phase explosives
journal, January 2015

  • Schoch, Stefan; Nikiforakis, Nikolaos
  • Physics of Fluids, Vol. 27, Issue 1
  • DOI: 10.1063/1.4905337

The Compressibility of Media under Extreme Pressures
journal, September 1944

  • Murnaghan, F. D.
  • Proceedings of the National Academy of Sciences, Vol. 30, Issue 9
  • DOI: 10.1073/pnas.30.9.244

Spherical Solutions of an Underwater Explosion Bubble
journal, January 1998

  • Wardlaw, Andrew B.; Mair, Hans U.
  • Shock and Vibration, Vol. 5, Issue 2
  • DOI: 10.1155/1998/690105

Underwater Explosions
journal, October 1948

  • Cole, Robert H.; Weller, Royal
  • Physics Today, Vol. 1, Issue 6
  • DOI: 10.1063/1.3066176

Shock waves in water
journal, July 1988


High resolution numerical simulation of ideal and non-ideal compressible reacting flows with embedded internal boundaries
journal, January 1997

  • Xu, S.; Aslam, T.; Stewart, D. S.
  • Combustion Theory and Modelling, Vol. 1, Issue 1
  • DOI: 10.1080/713665233

Convex ENO High Order Multi-dimensional Schemes without Field by Field Decomposition or Staggered Grids
journal, May 1998

  • Liu, Xu-Dong; Osher, Stanley
  • Journal of Computational Physics, Vol. 142, Issue 2
  • DOI: 10.1006/jcph.1998.5937

An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws
journal, March 2008

  • Borges, Rafael; Carmona, Monique; Costa, Bruno
  • Journal of Computational Physics, Vol. 227, Issue 6
  • DOI: 10.1016/j.jcp.2007.11.038

An approximate linearised riemann solver for the Euler equations for real gases
journal, February 1988


Inviscid flux-splitting algorithms for real gases with non-equilibrium chemistry
journal, October 1990


Positivity-preserving method for high-order conservative schemes solving compressible Euler equations
journal, June 2013

  • Hu, Xiangyu Y.; Adams, Nikolaus A.; Shu, Chi-Wang
  • Journal of Computational Physics, Vol. 242
  • DOI: 10.1016/j.jcp.2013.01.024

Reaction-zone structure of a steady-state detonation wave in a cylindrical charge
journal, April 1997


Multiple interacting liquids
journal, July 2006

  • Losasso, Frank; Shinar, Tamar; Selle, Andrew
  • ACM Transactions on Graphics, Vol. 25, Issue 3
  • DOI: 10.1145/1141911.1141960

Motion of Multiple Junctions: A Level Set Approach
journal, June 1994

  • Merriman, Barry; Bence, James K.; Osher, Stanley J.
  • Journal of Computational Physics, Vol. 112, Issue 2
  • DOI: 10.1006/jcph.1994.1105

A multiple level set method for modeling grain boundary evolution of polycrystalline materials
journal, June 2008

  • Zhang, Xinwei; Chen, Jiun-Shyan; Osher, Stanley
  • Interaction and multiscale mechanics, Vol. 1, Issue 2
  • DOI: 10.12989/imm.2008.1.2.191

A fast marching level set method for monotonically advancing fronts.
journal, February 1996

  • Sethian, J. A.
  • Proceedings of the National Academy of Sciences, Vol. 93, Issue 4
  • DOI: 10.1073/pnas.93.4.1591

A partial differential equation approach to multidimensional extrapolation
journal, January 2004


A static PDE Approach for MultiDimensional Extrapolation Using Fast Sweeping Methods
journal, January 2014

  • Aslam, Tariq; Luo, Songting; Zhao, Hongkai
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 6
  • DOI: 10.1137/140956919

A fast sweeping method for Eikonal equations
journal, May 2004


Convergence analysis of the fast sweeping method for static convex Hamilton–Jacobi equations
journal, December 2016


An Isobaric Fix for the Overheating Problem in Multimaterial Compressible Flows
journal, January 1999

  • Fedkiw, Ronald P.; Marquina, Antonio; Merriman, Barry
  • Journal of Computational Physics, Vol. 148, Issue 2
  • DOI: 10.1006/jcph.1998.6129

A ghost fluid method for compressible reacting flows with phase change
journal, February 2013


Efficient Characteristic Projection in Upwind Difference Schemes for Hyperbolic Systems
journal, March 1998

  • Fedkiw, Ronald P.; Merriman, Barry; Osher, Stanley
  • Journal of Computational Physics, Vol. 141, Issue 1
  • DOI: 10.1006/jcph.1998.5913

Multi-phase simulation of ammonium nitrate emulsion detonations
journal, September 2013


Numerical resolution of pulsating detonation waves
journal, September 2000


Characteristic path analysis of confinement influence on steady two-dimensional detonation propagation
journal, January 2019


Analysis of algorithms for velocity estimation from discrete position versus time data
journal, January 1992

  • Brown, R. H.; Schneider, S. C.; Mulligan, M. G.
  • IEEE Transactions on Industrial Electronics, Vol. 39, Issue 1
  • DOI: 10.1109/41.121906