DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A multiphase Lagrangian discontinuous Galerkin hydrodynamic method for high-explosive detonation physics

Abstract

We present a new multidimensional multiphase Lagrangian discontinuous Galerkin (DG) hydrodynamic method that supports programmed burn and reactive burn high-explosive (HE) detonation models. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. The HE reaction mass fraction is represented using a nodal quantity. A new limiting approach is presented along with new special limiting conditions for the reactive burn model. The HE models are applied to individual quadrature points in the DG elements to better resolve the physics with the modal method. In the programmed burn model, this means that reaction times are calculated for each point separately. The Scaled Unified Reactive Front (SURF) reactive burn model uses a lead shock pressure to determine the reaction rate at the quadrature points, with a new lead shock pressure detector that works with the Taylor series polynomials. The temporal evolution of the governing equations is achieved with the total variation diminishing Runge–Kutta (TVD RK) time integration method. The robustness of the new limiting approach is tested using 1D and 2D strong shock gamma-law gas problems. The accuracy of the point-wise programmed burn model is tested using amore » 1D gamma-law gas verification problem and a 2D radial problem. The implementation of the SURF reactive burn model within the DG method is tested using a 1D shock-to-detonation problem simulating the behavior of PBX-9501 and a 2D confined slab problem simulating the behavior of PBX-9502, both represented by the Davis equations of state with a multiphase modeling method. All 2D results use quadratic elements, which have faces that can bend.« less

Authors:
ORCiD logo; ; ;
Publication Date:
Research Org.:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
USDOE Laboratory Directed Research and Development (LDRD) Program; USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1681047
Alternate Identifier(s):
OSTI ID: 1739967
Report Number(s):
LA-UR-19-26917
Journal ID: ISSN 2666-4968; S2666496820300224; 100022; PII: S2666496820300224
Grant/Contract Number:  
89233218CNA000001
Resource Type:
Published Article
Journal Name:
Applications in Engineering Science
Additional Journal Information:
Journal Name: Applications in Engineering Science Journal Volume: 4 Journal Issue: C; Journal ID: ISSN 2666-4968
Publisher:
Elsevier
Country of Publication:
United Kingdom
Language:
English
Subject:
42 ENGINEERING; Lagrangian; Hydrodynamics; Discontinuous Galerkin; High-explosive; Programmed burn; Reactive burn

Citation Formats

Lieberman, Evan J., Liu, Xiaodong, Morgan, Nathaniel R., and Burton, Donald E. A multiphase Lagrangian discontinuous Galerkin hydrodynamic method for high-explosive detonation physics. United Kingdom: N. p., 2020. Web. doi:10.1016/j.apples.2020.100022.
Lieberman, Evan J., Liu, Xiaodong, Morgan, Nathaniel R., & Burton, Donald E. A multiphase Lagrangian discontinuous Galerkin hydrodynamic method for high-explosive detonation physics. United Kingdom. https://doi.org/10.1016/j.apples.2020.100022
Lieberman, Evan J., Liu, Xiaodong, Morgan, Nathaniel R., and Burton, Donald E. Tue . "A multiphase Lagrangian discontinuous Galerkin hydrodynamic method for high-explosive detonation physics". United Kingdom. https://doi.org/10.1016/j.apples.2020.100022.
@article{osti_1681047,
title = {A multiphase Lagrangian discontinuous Galerkin hydrodynamic method for high-explosive detonation physics},
author = {Lieberman, Evan J. and Liu, Xiaodong and Morgan, Nathaniel R. and Burton, Donald E.},
abstractNote = {We present a new multidimensional multiphase Lagrangian discontinuous Galerkin (DG) hydrodynamic method that supports programmed burn and reactive burn high-explosive (HE) detonation models. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. The HE reaction mass fraction is represented using a nodal quantity. A new limiting approach is presented along with new special limiting conditions for the reactive burn model. The HE models are applied to individual quadrature points in the DG elements to better resolve the physics with the modal method. In the programmed burn model, this means that reaction times are calculated for each point separately. The Scaled Unified Reactive Front (SURF) reactive burn model uses a lead shock pressure to determine the reaction rate at the quadrature points, with a new lead shock pressure detector that works with the Taylor series polynomials. The temporal evolution of the governing equations is achieved with the total variation diminishing Runge–Kutta (TVD RK) time integration method. The robustness of the new limiting approach is tested using 1D and 2D strong shock gamma-law gas problems. The accuracy of the point-wise programmed burn model is tested using a 1D gamma-law gas verification problem and a 2D radial problem. The implementation of the SURF reactive burn model within the DG method is tested using a 1D shock-to-detonation problem simulating the behavior of PBX-9501 and a 2D confined slab problem simulating the behavior of PBX-9502, both represented by the Davis equations of state with a multiphase modeling method. All 2D results use quadratic elements, which have faces that can bend.},
doi = {10.1016/j.apples.2020.100022},
journal = {Applications in Engineering Science},
number = C,
volume = 4,
place = {United Kingdom},
year = {Tue Dec 01 00:00:00 EST 2020},
month = {Tue Dec 01 00:00:00 EST 2020}
}

Works referenced in this record:

The Escape of High Explosive Products: An Exact-Solution Problem for Verification of Hydrodynamics Codes
journal, December 2016

  • Doebling, Scott W.
  • Journal of Verification, Validation and Uncertainty Quantification, Vol. 1, Issue 4
  • DOI: 10.1115/1.4035031

Robust high order discontinuous Galerkin schemes for two-dimensional gaseous detonations
journal, January 2012

  • Wang, Cheng; Zhang, Xiangxiong; Shu, Chi-Wang
  • Journal of Computational Physics, Vol. 231, Issue 2
  • DOI: 10.1016/j.jcp.2011.10.002

A Lagrangian staggered grid Godunov-like approach for hydrodynamics
journal, February 2014

  • Morgan, Nathaniel R.; Lipnikov, Konstantin N.; Burton, Donald E.
  • Journal of Computational Physics, Vol. 259
  • DOI: 10.1016/j.jcp.2013.12.013

Phenomenological model of shock initiation in heterogeneous explosives
journal, January 1980

  • Lee, E. L.; Tarver, C. M.
  • Physics of Fluids, Vol. 23, Issue 12
  • DOI: 10.1063/1.862940

A higher-order Lagrangian discontinuous Galerkin hydrodynamic method for elastic–plastic flows
journal, July 2019

  • Lieberman, Evan J.; Morgan, Nathaniel R.; Luscher, Darby J.
  • Computers & Mathematics with Applications, Vol. 78, Issue 2
  • DOI: 10.1016/j.camwa.2018.08.020

Equation of state and reaction rate for condensed-phase explosives
journal, September 2005

  • Wescott, B. L.; Stewart, D. Scott; Davis, W. C.
  • Journal of Applied Physics, Vol. 98, Issue 5
  • DOI: 10.1063/1.2035310

A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme
journal, June 2019

  • Liu, Xiaodong; Morgan, Nathaniel R.; Burton, Donald E.
  • Journal of Computational Physics, Vol. 386
  • DOI: 10.1016/j.jcp.2019.02.008

The SURF model and the curvature effect for PBX 9502
journal, December 2012


Level Set Methods Applied to Modeling Detonation Shock Dynamics
journal, July 1996

  • Aslam, Tariq D.; Bdzil, John B.; Stewart, D. Scott
  • Journal of Computational Physics, Vol. 126, Issue 2
  • DOI: 10.1006/jcph.1996.0145

Efficient implementation of essentially non-oscillatory shock-capturing schemes
journal, August 1988


The Sedov Blast Wave as a Radial Piston Verification Test
journal, June 2016

  • Pederson, Clark; Brown, Bart; Morgan, Nathaniel
  • Journal of Verification, Validation and Uncertainty Quantification, Vol. 1, Issue 3
  • DOI: 10.1115/1.4033652

A Cell-Centered Lagrangian Scheme for Two-Dimensional Compressible Flow Problems
journal, January 2007

  • Maire, Pierre-Henri; Abgrall, Rémi; Breil, Jérôme
  • SIAM Journal on Scientific Computing, Vol. 29, Issue 4
  • DOI: 10.1137/050633019

Scaling of detonation velocity in cylinder and slab geometries for ideal, insensitive and non-ideal explosives
journal, May 2015


Conservative high-order discontinuous Galerkin remap scheme on curvilinear polyhedral meshes
journal, November 2020


A cell-centered Lagrangian Godunov-like method for solid dynamics
journal, August 2013


Experimental validation of detonation shock dynamics in condensed explosives
journal, December 2005

  • Lambert, David E.; Stewart, D. Scott; Yoo, Sunhee
  • Journal of Fluid Mechanics, Vol. 546, Issue -1
  • DOI: 10.1017/S0022112005007160

Program Burn Algorithms Based on Detonation Shock Dynamics: Discrete Approximations of Detonation Flows with Discontinuous Front Models
journal, December 2001

  • Bdzil, J. B.; Stewart, D. S.; Jackson, T. L.
  • Journal of Computational Physics, Vol. 174, Issue 2
  • DOI: 10.1006/jcph.2001.6942

A Lagrangian discontinuous Galerkin hydrodynamic method
journal, February 2018


On the structure and accuracy of programmed burn
journal, April 2006

  • Kapila, A. K.; Bdzil, J. B.; Stewart, D. S.
  • Combustion Theory and Modelling, Vol. 10, Issue 2
  • DOI: 10.1080/13647830500436540

A discontinuous Galerkin discretization for solving the two-dimensional gas dynamics equations written under total Lagrangian formulation on general unstructured grids
journal, November 2014

  • Vilar, François; Maire, Pierre-Henri; Abgrall, Rémi
  • Journal of Computational Physics, Vol. 276
  • DOI: 10.1016/j.jcp.2014.07.030

An approach for treating contact surfaces in Lagrangian cell-centered hydrodynamics
journal, October 2013

  • Morgan, Nathaniel R.; Kenamond, Mark A.; Burton, Donald E.
  • Journal of Computational Physics, Vol. 250
  • DOI: 10.1016/j.jcp.2013.05.015

Monotonicity in high-order curvilinear finite element arbitrary Lagrangian-Eulerian remap: MONOTONICITY IN HIGH-ORDER CURVILINEAR FINITE ELEMENT ALE REMAP
journal, October 2014

  • Anderson, R. W.; Dobrev, V. A.; Kolev, Tz. V.
  • International Journal for Numerical Methods in Fluids, Vol. 77, Issue 5
  • DOI: 10.1002/fld.3965

A higher-order Lagrangian discontinuous Galerkin hydrodynamic method for solid dynamics
journal, August 2019

  • Lieberman, Evan J.; Liu, Xiaodong; Morgan, Nathaniel R.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 353
  • DOI: 10.1016/j.cma.2019.05.006

A new high-order discontinuous Galerkin spectral finite element method for Lagrangian gas dynamics in two-dimensions
journal, April 2011


A high-order conservative remap for discontinuous Galerkin schemes on curvilinear polygonal meshes
journal, December 2019


Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates
journal, November 2018

  • Liu, Xiaodong; Morgan, Nathaniel R.; Burton, Donald E.
  • Journal of Computational Physics, Vol. 373
  • DOI: 10.1016/j.jcp.2018.06.073

Slope limiting for discontinuous Galerkin approximations with a possibly non-orthogonal Taylor basis: SLOPE LIMITING FOR DG APPROXIMATIONS
journal, September 2012

  • Kuzmin, Dmitri
  • International Journal for Numerical Methods in Fluids, Vol. 71, Issue 9
  • DOI: 10.1002/fld.3707