Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Journal Article · · Journal of Applied Physics
DOI:https://doi.org/10.1063/5.0285135· OSTI ID:3003171
 [1];  [2]
  1. Los Alamos National Laboratory (LANL), Los Alamos, NM (United States); Univ. of Colorado, Boulder, CO (United States)
  2. Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

Research Organization:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
Grant/Contract Number:
NA0003962; 89233218CNA000001
OSTI ID:
3003171
Report Number(s):
LA-UR--25-25358; 10.1063/5.0285135; DOPSR-25-T-2136
Journal Information:
Journal of Applied Physics, Journal Name: Journal of Applied Physics Journal Issue: 15 Vol. 138; ISSN 0021-8979; ISSN 1089-7550
Publisher:
American Institute of Physics (AIP)Copyright Statement
Country of Publication:
United States
Language:
English

References (32)

Calibration and uncertainty quantification for Davis Equation of State models for the High Explosive PBX 9501 products journal January 2024
pyOpt: a Python-based object-oriented framework for nonlinear constrained optimization journal May 2011
Thermochemical functions for gas-phase, 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), its condensed phases, and its larger reaction products journal August 2002
Complete equation of state for unreacted solid explosive journal February 2000
Array programming with NumPy journal September 2020
SciPy 1.0: fundamental algorithms for scientific computing in Python journal February 2020
Equation of state and reaction rate for condensed-phase explosives journal September 2005
Shock temperature dependent rate law for plastic bonded explosives journal April 2018
Magpie: A new thermochemical code
  • Ticknor, C.; Andrews, S. A.; Leiding, J. A.
  • SHOCK COMPRESSION OF CONDENSED MATTER - 2019: Proceedings of the Conference of the American Physical Society Topical Group on Shock Compression of Condensed Matter, AIP Conference Proceedings https://doi.org/10.1063/12.0000785
conference January 2020
AWSD calibration for the HMX based explosive PBX 9501
  • Aslam, Tariq D.; Price, Matthew A.; Ticknor, Christopher
  • SHOCK COMPRESSION OF CONDENSED MATTER - 2019: Proceedings of the Conference of the American Physical Society Topical Group on Shock Compression of Condensed Matter, AIP Conference Proceedings https://doi.org/10.1063/12.0000891
conference January 2020
PBX 9012 Davis reactant and product EOS calibration with uncertainty quantification conference January 2023
Uncertainty quantification for a multi-phase magnesium equation of state
  • Lindquist, Beth A.; Jadrich, Ryan B.
  • SHOCK COMPRESSION OF CONDENSED MATTER - 2022: Proceedings of the Conference of the American Physical Society Topical Group on Shock Compression of Condensed Matter https://doi.org/10.1063/12.0020372
conference January 2023
Uncertainty quantified reactant and product equation of state for Composition B conference January 2023
Uncertainty quantification for equation of state using heterogeneous data sources conference January 2024
A reactive flow model for the 3,3′-diamino-4,4′-azoxyfurazan based plastic bonded explosive (PBX 9701) journal December 2021
Uncertainty quantification for a multi-phase carbon equation of state model journal April 2022
Uncertainty quantification for high explosive reactant and product equations of state journal August 2023
Weighted Bayesian uncertainty quantification for the high explosive reactants using limited data journal February 2025
AWSD reactive flow model for PBX 9404 journal April 2025
Sensitivity of an integrated experiment to uncertainty in the high explosive equations of state journal March 2025
Detonation waves in PBX 9501 journal December 2006
emcee : The MCMC Hammer
  • Foreman-Mackey, Daniel; Hogg, David W.; Lang, Dustin
  • Publications of the Astronomical Society of the Pacific, Vol. 125, Issue 925 https://doi.org/10.1086/670067
journal March 2013
Learning about physical parameters: the importance of model discrepancy journal October 2014
Matplotlib: A 2D Graphics Environment journal January 2007
Bayesian calibration of computer models journal August 2001
Estimating Physics Models and Quantifying Their Uncertainty Using Optimization With a Bayesian Objective Function journal March 2019
Posterior Covariance Matrix Approximations journal March 2024
Exploring the Uncertainty in the Equation of State for a High Explosive Fit to Heterogeneous Data conference November 2019
Variational Bayesian Calibration of a PTW Material Strength Model for OFHC Copper conference May 2023
SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization journal January 2005
Algorithm 733; TOMP---Fortran modules for optimal control calculations journal September 1994
pyOptSparse: A Python framework for large-scale constrained nonlinear optimization of sparse systems journal October 2020

Similar Records

Posterior Covariance Matrix Approximations
Journal Article · Mon May 13 00:00:00 EDT 2024 · Journal of Verification, Validation and Uncertainty Quantification · OSTI ID:2406578

Calibration and uncertainty quantification for Davis Equation of State models for the High Explosive PBX 9501 products
Journal Article · Tue Jan 30 23:00:00 EST 2024 · Propellants, Explosives, Pyrotechnics · OSTI ID:2283470