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.
We present the first volumetric measurements of magnetic field, density, and temperature in laser-produced plasmas and shock waves. The results were obtained in one of the only rep-rated (1 Hz) HED experiment currently operating.
GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.
Traditional mesh generation approaches are labor intensive and have limited robustness when applied to parametric design exploration and optimization of complex geometries. While automatic mesh generation approaches exist, they tend to generate tetrahedral or mixed-hybrid meshes which are generally unsuitable for physics applications with strong shock waves, thin boundary layers, and strong gradients. In addition, simulations sizes in the billions of cells are becoming more common with traditional mesh generation methods quickly reaching scalability limits. CrossLink offers a topology-based mesh generation approach with unstructured block-filling methods and a scalable mesh generation engine. In addition, CrossLink incorporates a python based API for seamless workflow integration and robust repeatability of the geometry handling and mesh generation process. This makes it ideal for parametric design study and optimization of complex geometries. Finally, future versions of CrossLink will offer a parametric mesh capability that optimizes a high-order mesh and enables reconstruction of the final mesh in memory by the physics solver.
The purpose of this experimental series is to validate the reactive burn model in the PAGOSA hydrodynamic code. It will also help validate the EOS and Strength models for the aluminum and tantalum cylinders. C4 high explosive (HE) is used to drive an aluminum cylinder into an aluminum sandwiched PBX-9501. The follow-on experiment will drive the aluminum cylinder into a tantalum sandwiched PBX-9501. The experiment will provide time of arrival (TOA) data for the burn front in the HE through the metal cylinders using Fiber Bragg Gratings (FBG). Photon Doppler Velocimetry (PDV) will provide shock wave TOA and the velocity of the inward moving cylinder. The PAGOSA code will be used to predict TOA and cylinder wall velocity using PDV and tracers.
The essence of the project was to implement a program in python that demonstrated certain models put forth by fusion physicists. Across the time spent at Los Alamos National Laboratory, much was spent programming. Though fundamentally programmatical in nature, the task contained considerable difficulty hid within mathematically deriving and manipulating the models while maintaining data discretization. The result of solving these complications gave invaluable experience in applied mathematics. The final product, a piece of lengthy code, ended up having the full range of desired capabilities, with plans for more. Given the laser power history from a pulse shot on NIF, Omega, or any other Inertial Confinement Fusion facility, the code maps the terrain of radiation temperature inside the hohlraum throughout the period of the shot. Additionally, shock wave convergence prediction was implemented as additional capability. As the name of the project suggests, this tool acts as a compass of direction for researchers designing laser power profiles. Further development of this tool will undoubtedly assist in the research conducted on single and double-shell indirect ICF, consequently expanding the capabilities of operating stewardship over the national nuclear stockpile.
The purpose of this experimental series is to validate the reactive burn model in the PAGOSA hydrodynamic code. It will also help validate the EOS and Strength models for the aluminum and tantalum cylinders. C4 high explosive (HE) is used to drive an aluminum cylinder into a tantalum sandwiched PBX‐9501. The experiment will provide time of arrival (TOA) data for the burn front in the HE through the metal cylinders using Fiber Bragg Gratings (FBG). Photon Doppler Velocimetry (PDV) will provide shock wave TOA and the velocity of the inward moving cylinder. The PAGOSA code will be used to predict TOA and cylinder wall velocity using PDV and tracers.
Cerium ejecta particles form after a shock wave impacts a cerium plate with a rough surface or prescribed surface perturbation. When these particles are ejected into hydrogen or deuterium gas the material reacts exothermically and raises the temperature of the particle. Two models are developed to capture the heat and mass transfer that occurs during this process. One model assumes that for larger solid ejecta particles, the loss of mass is negligible while for smaller liquid ejecta particles, the loss of mass during reaction is accounted for. The loss of mass is assumed to occur from the hydride layer being extremely brittle and the fact that it will be relatively easy for the hydride material on the outside of the particle to go into tension, crack and flake off. The model is tested by performing simulations similar to the original experiments that motivated the development of these models and comparing the particle temperatures. While the model is able to capture some general features of the flowfield, inaccuracies clearly exist that must be addressed and the results underscore the need for more isolated experiments where fewer assumptions and physical processes are present simultaneously.
The Explosive Shock Image Analysis System (ESAIS) is designed to provide a semiautomated quantitative analysis of fast-frame imagery taken during the detonation of high explosive (HE) rate sticks. ESIAS does this by assisting the user in identifying the HE stick geometry on the image (i.e., sides of the undetonated HE stick, detonation front, an axial reference point to assist in the inter-comparison of images); scaling the images based on the known diameter of the undetonated HE stick; identifying the profiles of various surfaces resulting from the detonation process such as the shock wave and product gases or the outer pipe surface; calculating the detonation speed of the HE by tracking the distance it travels between two images with a known time lapse between the images; and determining the slit velocities resulting from the detonation. ESIAS will also import scope traces from oscilloscopes and assist in identifying the time of arrival of significant peaks.
Current and future high-power accelerators put severe requirements on materials used for target and beam windows and target facilities have been recognized as a critical challenge in development of future particle accelerators. In accelerators, window and target materials are exposed to extreme conditions, which include bombardment with very high energy protons (1- 100 GeV) and thermomechanical shock waves. Radiation can cause direct damage in the material, and it leads to production of transmutation products (especially helium), both phenomena having a potential adverse effect on the stability and durability of the target/window material. At high enough temperatures, He can aggregate to form gas bubbles, which in turn cause significant dimensional changes (swelling), enable easy crack propagation, and eventually cause failure by fracture. On the other hand, if the temperature is too low, radiation damage accumulates in the form of internal defects (e.g., dislocations), leading to hardening and a decreased ductility of the material. In this project, we will focus on beryllium since it is considered to be one of the candidate materials for beam windows and targets in the next-generation proton accelerators, e.g., the Long Baseline Neutrino Facility (LBNF). Radiation effects in Be have been studied in the context of nuclear fusion reactor applications. However, key differences exist between reactor and accelerator conditions, including neutron vs. proton irradiation, continuous vs. pulsed beam flux, much higher energies of bombarding particles in accelerators, and higher operating temperatures for typical reactors. For example, the impact of beam pulsing on the radiation damage and the He bubble kinetics is largely unknown. While results obtained on Be from fusion research might not be directly transferrable to understanding target materials, there is an opportunity to bring state-of-the-art tools from materials research in nuclear reactors to aid design of target and beam window materials in high-power accelerators. To this end, the overarching goal of this project are to develop an experimentally-validated computational framework capable of predicting radiation damage evolution in beryllium relevant to beam window and target conditions, focusing on He bubble formation and growth as a function of irradiation temperature. Our model will be based on the cluster dynamics formalism, where size distribution of defects and He bubbles is simulated as a function of time, temperature, and radiation dose. Parameters for the model will be taken from published experiments and from high-fidelity atomistic simulations proposed in this project. In addition, we will carry out a series of targeted ex-situ and in-situ dual-beam experiments using low-energy protons to provide critical data for validation of the model on the effects of radiation on He clustering, He bubble distribution, and dislocation loop density/size in proton irradiated Be.
This user guide supports a curated dataset of 71 meteor events recorded between 2006 and 2011 in Southwestern Ontario, Canada. Each event was simultaneously observed by ground-based optical cameras and an infrasound array, providing a rare opportunity to examine meteor trajectories and acoustic signals from the same atmospheric entry events. The dataset includes raw and processed optical data, meteor trajectories, photometric light curves, infrasound waveforms, and atmospheric specifications relevant for acoustic modeling. The archive is structured to support reproducible research in meteor physics, atmospheric acoustics, and shock wave analysis. It is organized following transparent file naming conventions and structured folders to facilitate scientific reuse, comparison, and integration across research domains. The dataset is freely available on Zenodo, doi: 10.5281/zenodo.15868512.
Radiation-hydrodynamics simulations are critical to the NNSA complex for understanding high energy density physics and for addressing problems in national security science. However, current rad-hydro codes typically do not account for viscosity, which generally leads to large uncertainties in problems involving mixing of materials under shock loading. This report details progress made towards experimental measurements of fluid viscosity at high pressures and temperatures achieved by shock wave compression techniques, as well as complementary material model development and hydrodynamic simulations.
The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.
Exploding Bridge Wire (EBW) is an important design for initiating detonators relevant to several LLNL systems. In such designs, a metal bridge wire is in direct contact with a porous secondary high explosive (HE), typically pressed to around 50% of its crystalline density. Initiation is believed to occur through shock waves generated when the bridge wire is made to explode via the rapid passing of a strong current pulse. Permeametry is an important characterization tool for the porous powdered HE (PETN, HMX, etc.) component in such devices. Such powders are well known to coarsen with age, with an accompanying deterioration in performance, which includes increased function time and sensitivity. Among various material characteristics, one that has been shown to correlate the most with such performance parameters is the flow-permeable surface area (FSSA), measured using air-permeametry. Commercially available permeametry apparatuses, such as that from Fisher Scientific or a more modern version from Micromeritics suffer from a few drawbacks, most notably as related to monitoring age-related changes in FSSA. This is because it requires re-compaction of the aged powder into the permeametry tube, which in most cases leads to an unwanted increase in FSSA. To address this and a few other shortcomings, we have been developing a modern version of the permeameter, Gee-Reinstein-Maiti Modern Air-Permeametry System (GRM-MAPS). In this report, we summarize progress made in FY25 toward perfecting and calibrating the instrument, which we aim at deploying within surveillance in FY26.
Chondrules are tiny particles that occur in stony meteorites and are considered as the building blocks of early asteroids and planets. It is believed that they were formed by the fast heating of the dust in the solar nebula. To date, there is no lab-scale experimental study of the formation of chondrules from the initial gas phase precursors following fast heating and crystallisation. The motivation of this work is a pre-trial study of the formation of chnodrule-like particles. The formation of meteorites in the space environment is associated with the aggregation of small particles or molecular clouds under the influence of shock waves or high-energy gas discharges in the solar nebula. In this work, the properties of product formation at the nanoscale-level were investigated using different feedstock materials which are the dominant elements in the meteorite. The structural and morphological properties of the synthesised Si-Fe nanomaterials were analysed by scanning/transmission electron microscopy (SEM/TEM), and chemical composition was analysed by X-ray energy-dispersive spectroscopy (EDS). The identification of crystalline phases was carried out by X-ray diffraction (XRD), whereas the presence of an Fe-Si system in the synthesised particles was demonstrated by Mössbauer spectroscopy. The obtained materials were exposed to the relatively high-energy pulsed plasma beam on the substrate with the aim to emulate the possible fast heating and melting of the formed nanoparticles. The formation steps of growing synthetic (engineered) chondro-like particles and nanostructures in laboratory conditions is discussed.
We present two-dimensional (2D) optical Thomson scattering measurements of electron density and temperature in laser-produced plasmas. The novel instrument directly measures ne(x,y) and Te(x,y) in two dimensions over large spatial regions (cm2) with sub-mm spatial resolution, by automatically translating the scattering volume while the plasma is produced repeatedly by irradiating a solid target with a high-repetition-rate laser beam (10 J, ∼1012 W/cm2, 1 Hz). In this paper, we describe the design and motorized auto-alignment of the instrument and the computerized algorithm that autonomously fits the spectral distribution function to the tens-of-thousands of measured scattering spectra, and captures the transition from the collective to the non-collective regime with distance from the target. As an example, we present the first 2D scattering measurements in laser-driven shock waves in ambient nitrogen gas at a pressure of 0.13 mbar.
We investigate the neutrino flavor change effects due to neutrino self-interaction and shock wave propagation, as well as the matter effects on the neutrino process in core-collapsing supernovae (CCSNe). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamics model for SN1987A. For the presupernova model, we take an updated model, adjusted to explain SN1987A, which employs recent developments in the (n, γ) reaction rates for nuclei near the stability line (A ~ 100). As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and nonequivalent luminosity models. The latter is taken from a synthetic analysis of CCSN simulation data, which quantitatively presented the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated using a shell model for light nuclei and a quasiparticle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ( 7 Li, 7 Be, 11 B, and 11 C) and the heavy nuclei ( 92 Nb, 98 Tc, 138 La, and 180 Ta) produced by the neutrino process. The light nuclei abundances turn out to be sensitive to the Mikheyev–Smirnov–Wolfenstein (MSW) region around O-Ne-Mg layer while the heavy nuclei are mainly produced prior to the MSW region. Through detailed analyses, we find that neutrino self-interaction becomes a key ingredient, in addition to the MSW effect, for understanding the neutrino process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. The main nuclear reactions for each nucleus are also investigated in detail.
It is argued that quantum states of geometry, like those of particles, should be coherent on light cones of any size. An exact classical solution, the gravitational shock wave of a relativistic point particle, is used to estimate gravitational drag from coherent energy flows, and the expected gravitational effect of virtual transverse vacuum energy fluctuations on surfaces of causal diamonds. It is proposed that the appropriately spacetime-averaged gravitational effect of the Standard Model vacuum state leads to the observed small nonzero value of the cosmological constant, dominated by gravitational drag of virtual gluonic strings at the strong interaction scale.
Methods for calibration of non-contact velocity measurements and systems for implementing the same are described. Generally, the method comprises inducing a shock wave into a sample at a stress intensity that varies across the sample's elastic limit, which corresponds to the elastic-plastic state transition of the sample. That transition state may be at the sample's Hugoniot elastic limit. The velocity of the sample is measured using a non-contact velocity measurement instrument such as a velocimeter. The measurement may be compared to a predicted velocity or a velocity measurement made by another system to determine calibration parameters.