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

Title: (U) Physics Validation of the RMI-Based Ejecta Source Model Implementation in FLAG


This document describes LANL's source model for RMI-based ejecta.

  1. Los Alamos National Lab. (LANL), Los Alamos, NM (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:
Report Number(s):
DOE Contract Number:
Resource Type:
Technical Report
Country of Publication:
United States

Citation Formats

Tregillis, Ian Lee. (U) Physics Validation of the RMI-Based Ejecta Source Model Implementation in FLAG. United States: N. p., 2017. Web. doi:10.2172/1392822.
Tregillis, Ian Lee. (U) Physics Validation of the RMI-Based Ejecta Source Model Implementation in FLAG. United States. doi:10.2172/1392822.
Tregillis, Ian Lee. 2017. "(U) Physics Validation of the RMI-Based Ejecta Source Model Implementation in FLAG". United States. doi:10.2172/1392822.
title = {(U) Physics Validation of the RMI-Based Ejecta Source Model Implementation in FLAG},
author = {Tregillis, Ian Lee},
abstractNote = {This document describes LANL's source model for RMI-based ejecta.},
doi = {10.2172/1392822},
journal = {},
number = ,
volume = ,
place = {United States},
year = 2017,
month = 9

Technical Report:

Save / Share:
  • The Los Alamos Physics and Engineering Models (PEM) program has developed a model for Richtmyer-Meshkov instability (RMI) based ejecta production from shock-melted surfaces, along with a prescription for a self-similar velocity distribution (SSVD) of the resulting ejecta particles. We have undertaken an effort to validate this source model using data from explosively driven tin coupon experiments. The model’s current formulation lacks a crucial piece of physics: a method for determining the duration of the ejecta production interval. Without a mechanism for terminating ejecta production, the model is not predictive. Furthermore, when the production interval is hand-tuned to match time-integrated massmore » data, the predicted time-dependent mass accumulation on a downstream sensor rises too sharply at early times and too slowly at late times because the SSVD overestimates the amount of mass stored in the fastest particles and underestimates the mass stored in the slowest particles. The functional form of the resulting m(t) is inconsistent with the available time-dependent data; numerical simulations and analytic studies agree on this point. Simulated mass tallies are highly sensitive to radial expansion of the ejecta cloud. It is not clear if the same effect is present in the experimental data but if so, depending on the degree, this may challenge the model’s compatibility with tin coupon data. The current implementation of the model in FLAG is sensitive to the detailed interaction between kinematics (hydrodynamic methods) and thermodynamics (material models); this sensitivity prohibits certain physics modeling choices. The appendices contain an extensive analytic study of piezoelectric ejecta mass measurements, along with test problems, excerpted from a longer work (LA-UR-17-21218).« less
  • The BHR-2 turbulence model, developed at Los Alamos National Laboratory for variable density and compressible flows, is implemented in an Arbitrary Lagrangian-Eulerian hydrocode, FLAG. The BHR-2 formulation is discussed, with emphasis on its connection to multi-component flow formulations that underlie FLAG's treatment of multi-species flow. One-dimensional and two-dimensional validation tests are performed and compared to experiment and Eulerian simulations. Turbulence is an often studied and ubiquitous phenomenon in nature, and modeling its effects is essential in many practical applications. Specifically the behavior of turbulence in the presence of strong density gradients and compressibility is of fundamental importance in applications rangingmore » from Inertial Confinement Fusion (ICF) [1], supernovae [2], and atmospheric flows. The BHR closure approach [3] seeks to model the physical processes at work in variable density turbulence including Kelvin-Helmholtz (KH) [4], Rayleigh-Taylor (RT) [5], and Richtmyer-Meshkov (RM) [6], driven turbulence. The effectiveness of the BHR-2 implementation has been demonstrated for variable density mixing in the KH, RT, and RM cases in an Eulerian framework [7]. The primary motivation of the present work is to implement the BHR-2 turbulence model in the Arbitrary Lagrangian-Eulerian (ALE) hydrodynamics code FLAG. The goal is not only to demonstrate results in agreement with previous Eulerian calculations, but also document behavior that arises from the underlying differences in code philosophy.« less
  • This report describes the SURF reactive burn model and its use cases.
  • Properly validated and calibrated reactive burn models (RBM) can be useful engineering tools for assessing high explosive performance and safety. Experiments with high explosives are expensive. Inexpensive RBM calculations are increasingly relied on for predictive analysis for performance and safety. This report discusses the validation of Menikoff and Shaw’s SURF reactive burn model, which has recently been implemented in the FLAG code. The LANL Gapstick experiment is discussed as is its’ utility in reactive burn model validation. Data obtained from pRad for the LT-63 series is also presented along with FLAG simulations using SURF for both PBX 9501 and PBXmore » 9502. Calibration parameters for both explosives are presented.« less