SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
The essential idea of this algorithm is to compute the angular distribution of a vector quantity, then create an informative image. In our example, we will compute the angular distribution of linear momentum from an xRage simulation of an exploding shaped charge. We will then explore one possible method for adding information to the resulting image.
It has previously been shown that the rate parameters of the SURF reactive burn model can be adjusted to fit the gap-stick data for PBX 9501 [Johnson et al., 2018, fig 4]. A recent recalibration of the SURF model for PBX 9501 (lot 730-010 at ρ = 1.837 g/cc 3 ) used different reactants and products EOS and an updated fitting form for the burn rate [Menikoff, 2021]. The new calibration fit gap-stick data [Hill et al., 2021] together with data from SDT experiments and curvature effect data. Here the results of gap-stick simulations using the xRage code are shown. They provide insight into the reactive flow in the gap-stick experiment.
The primary goals for the May 2022 shot day at OMEGA are to: image an imploding inner cylinder as proof-of-principle for double cylinder experiments, measure the growth of pre-seeded perturbations on the inner cylinder, and measure the axial non-uniformity of the implosion with a sidelighter that provides a transverse view of the target. We will compare these experimental results to xRAGE (2D-only) and FLASH (2D and 3D) radiation-hydrodynamics calculations.
A detailed discussion about the PTW flow-stress model has been presented, including theoretical considerations, model calibrations to experimental data, and model implementations in LANL ASC hydro-codes. The main findings are: In order to reproduce the intended behavior of the PTW model at high strain rates, the following conditions on the model parameters need to be valid: y 1 ≥ s 0 and y 2 ≥ β ; When the conditions above are not satisfied, PTW and PTW_Mod1 have slightly different behavior in the phonon-drag regime, and should be treated as different models. PTW exhibits a discontinuity above the critical strain rate. PTW_Mod1 does not show such a discontinuity but implies a different strain-rate dependence at high strain rates. Examples of this behavior include older calibrations for stainless steel and depleted uranium; In the thermal-activation regime, some flexibility was introduced in the original PTW model to allow for different rate-dependent behaviors at intermediate strain rates. Example: copper. However, many material calibrations do not require this extra flexibility, and the complexity and ambiguity of the intended behavior can be partially blamed for the inconsistent model implementations across ASC codes; The limit for the material constant p going to zero is well defined and leads to a simplified flow-stress equation, discontinuity free; Modern calibrations of the PTW parameters are done using Impala. Constraints on the PTW parameters are enforced; Implementation of the PTW model in LANL ASC hydro-codes suffered of several shortcom ings. At the time of writing, FLAG and xRAGE contain consistent and correct implemen tations of the PTW and PTW_Mod1 models. Corresponding changes to Pagosa have been suggested.
The following slides summarize the effect that AMR has on results obtained using the inviscid version of xRAGE, i.e., when solving the Euler equations. The simulations are 2D and were done in the context of the ModCons experiments, where a shock travels from a heavy material (PAI of density 1.45 g/cm 3 ) to a light material (foam of density 0.1 g/cm 3 ). The shock is initialized using the 1D shock relations and is sustained by a freeze region during the entire simulation (0 ≤ t ≤ 6 ns). The initial perturbation has a wavelength of λ = 100 µm and an initial amplitude of a 0 = 10 µm. Five cases are considered: four with AMR with decreasing resolutions (from 1 µm to 1/8 µm at the interface) and one without AMR with 1/8 µm resolution in the entire domain.
This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.
This document summarizes the recent progress towards incorporating charged-particle radiochemical (RadChem) analysis in different LANL codes within the Advanced Scientific Computing (ASC) Program. The first section is focused on xRAGE and the implementation of the necessary physics packages to build an in-line RadChem capability. The second section provides a summary of the current RadChem postprocessing capabilities.
Purpose: Understand impact of heterogeneous dopants in high-Z metal coatings on instability seeding; Validate xRAGE modeling of fill tube physics and mitigation strategies
The objective of this report is threefold. First, it details a method for deriving a semi-analytic solution to the Noh Problem when using a “black-box” equation of state. Such capability allows us to perform verification on complicated, more realistic equations of state. Examples include Steinberg equations of state for materials and tabulated equations of state. The second objective is to apply the methodology to verify the singularity-eos equation of state library. We do so by solving the Rankine-Hugoinot jump conditions for the Noh Problem, ensuring singularity derives the correct solution and comparing the error to an exact implementation of the equation of state. The third objective is to perform verification of the xRAGE Eulerian hydrodynamics code when interfaced with singularity. We provide the theory, analysis, documentation for a python implementation of the proposed solver, and verification results.
Explore the source record for details and available documents.