Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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.

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Search for Nova Presolar Grains: $γ$-Ray Spectroscopy of 34 Ar and its Relevance for the Astrophysical 33 Cl($p,γ$) Reaction

The discovery of presolar grains in primitive meteorites has launched a new era of research in the study of stellar nucleosynthesis. However, the accurate classification of presolar grains as being of specific stellar origins is particularly challenging. Recently, it has been suggested that sulfur isotopic abundances may hold the key to definitively identifying presolar grains with being of nova origins and, in this regard, the astrophysical Cl 33 ( p , γ ) Ar 34 reaction is expected to play a decisive role. As such, we have performed a detailed γ -ray spectroscopy study of Ar 34 . Excitation energies have been measured with high precision and spin-parity assignments for resonant states, located above the proton threshold in Ar 34 , have been made for the first time. Uncertainties in the Cl 33 ( p , γ ) reaction have been dramatically reduced and the results indicate that a newly identified ℓ = 0 resonance at E r = 396.9 ( 13 ) keV dominates the entire rate for T = 0.25 – 0.40 GK . Furthermore, nova hydrodynamic simulations based on the present work indicate an ejected S 32 / S 33 abundance ratio distinctive from type-II supernovae and potentially compatible with recent measurements of a presolar grain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Novel Hot-Spot Ignition Designs for Inertial Confinement Fusion with Liquid-Deuterium-Tritium Spheres

A new class of ignition designs is proposed for inertial confinement fusion (ICF) experiments. These designs are based on the hot-spot-ignition approach, but instead of conventional targets that comprise of spherical shells with thin frozen deuterium{tritium (DT) layers, homogeneous-density (solid or liquid) DT spheres are used where the lower-density central region and higher-density shell are created dynamically by appropriately shaping the laser pulse. These offer several advantages, including simplicity in target production and lower sensitivity to both laser imprint and physics uncertainty in shock interaction with the ice-vapor interface. The design evolution starts by launching an ~1-Mbar shock into a homogeneous DT sphere. After bouncing from the center, the reflected shock reaches the outer surface of the sphere and the shocked material starts to expand outward until its pressure drops below the ablation pressure. At this point, an adjustment shock is launched inward by supporting ablation pressure. This shock compresses the ablator and fuel, forming a shell. The shell is then accelerated and compressed by appropriately shaping the drive laser pulse, similar to the conventional thin-shell, hot-spot designs. Here, this paper demonstrates the feasibility of the new concept using hydrodynamic simulations and discusses the advantages and disadvantages of the concept compared with more-traditional ICF designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Constraining the P 30 ( p , γ ) S 31 Reaction Rate in ONe Novae via the Weak, Low-Energy, β -Delayed Proton Decay of Cl 31

The 30 P(p,γ) 31 S reaction plays an important role in understanding the nucleosynthesis of A ≥ 30 nuclides in oxygen-neon novae. The Gaseous Detector with Germanium Tagging was used to measure 31 Cl β-delayed proton decay through the key J π = 3/2 + , 260-keV resonance. The intensity $I$$^{260}_{βp}$ = $8.3$$^{+1.2}_{–0.9}$ × 10 –6 represents the weakest β-delayed, charged-particle emission ever measured below 400 keV, resulting in a proton branching ratio of Γ p /Γ = $2.5$$^{+0.4}_{–0.3}$ × 10 –4 . Here, by combining this measurement with shell-model calculations for Γγ and past work on other resonances, the total 30 P(p,γ) 31 S rate has been determined with reduced uncertainty. The new rate has been used in hydrodynamic simulations to model the composition of nova ejecta, leading to a concrete prediction of 30 Si: 28 Si excesses in presolar nova grains and the calibration of nuclear thermometers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Inverse Bremsstrahlung Absorption

Inverse bremsstrahlung absorption was measured based on transmission through a finite-length plasma that was thoroughly characterized using spatially resolved Thomson scattering. Expected absorption was then calculated using the diagnosed plasma conditions while varying the absorption model components. To match data, it is necessary to account for: (1) the Langdon effect; (2) laser-frequency (rather than plasma-frequency) dependence in the Coulomb logarithm, as is typical of bremsstrahlung theories but not transport theories; and (3) a correction due to ion screening. Radiation-hydrodynamic simulations of inertial confinement fusion implosions have to date used a Coulomb logarithm from the transport literature and no screening correction. Here, we anticipate that updating the model for collisional absorption will substantially revise our understanding of laser-target coupling for such implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hydrodynamic Mechanism for Clumping along the Equatorial Rings of SN1987A and Other Stars

An explanation for the origin and number of clumps along the equatorial ring of Supernova 1987A has eluded decades of research. Our linear analysis and hydrodynamic simulations of the expanding ring prior to the supernova reveal that it is subject to the Crow instability between vortex cores. Here, the dominant wave number is remarkably consistent with the number of clumps, suggesting that the Crow instability stimulates clump formation. Although the present analysis focuses on linear fluid flow, future nonlinear analysis and the incorporation of additional stellar physics may further elucidate the remnant structure and the evolution of the progenitor and other stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Direct Evidence of Multispecies Hydrodynamics in Ignition-Scale Hohlraums

A targeted experiment at the National Ignition Facility (NIF) confirms the presence of multispecies hydrodynamics in inertial confinement fusion hohlraums relevant to ignition. Here, the effects are identified by filling the gold hohlraum with a deuterium-tritium (DT) gas mixture instead of helium. As the hohlraum is heated by the NIF lasers, it implodes inward, compressing and heating the DT, which leads to fusion. The resulting DT-fusion neutrons are measured in space, time, yield, angle, and energy. A distinct, peaked, triangular shape in the radial neutron emission profile provides evidence of a “leaky piston” effect caused by the interpenetration of DT into the expanding gold. This process reduces the reversibility of DT compression and decreases neutron generation on the temporal trailing edge compared to the leading edge. These results are well described by multispecies hydrodynamics simulations, which reproduce the observed spatial and temporal features, as well as the total neutron yield and angularly resolved energy spectra. In contrast, conventional simulations that use only single-species physics fail to match the experimental data. They overpredict the neutron yield, produce a flatter-than-observed spatial profile, and show excessive emission at late times compared to the experimental measurements.

Higginson, Drew P. [Lawrence Livermore National La↗
DOE OSTI2024

Laser pulse-length dependent ablation and shock generation in silicon at 5 × 10 14 W/ cm 2 intensities

The effect of laser pulse duration on energy coupling into a planar silicon target is investigated in experiments at the OMEGA-EP facility by varying the laser pulse length τ —spanning 3 orders of magnitude from 100 ps to 10 ns—while maintaining a constant peak laser intensity, I 0 = 5 × 10 14 W/ cm 2 . In theoretical models, the ablation pressure primarily scales for a given material with laser intensity and wavelength, which are all fixed variables here, allowing us to explore the specific role of laser pulse duration. Two-dimensional radiation-hydrodynamics simulations benchmarked with optical probing of the expanding plasma show that the pulse duration is critical for the ablation pressure to reach a steady state. Moreover, the pulse duration impacts shock decay and multiple wave effects, which strongly dictate the evolving shock profile that propagates within the laser-shocked target as ultimately measured by rear-surface diagnostics. The shock velocities inferred from the theoretical model, after considering shock decay, impedance matching, and shock Hugoniot, are found to be in good agreement with velocimetry measurements. However, discrepancies are observed with simulations for the shorter (0.1 ns) and longer (10 ns) pulse durations, which are respectively attributed to unaccounted contributions of kinetic absorption mechanisms and instabilities in simulations. Published by the American Physical Society 2024

Bailly-Grandvaux, M. (ORCID:0000000175294013)↗
DOE OSTI2025

Dark Matter Velocity Distributions for Direct Detection: Astrophysical Uncertainties Are Smaller Than They Appear

The sensitivity of direct detection experiments depends on the phase-space distribution of dark matter near the Sun, which can be modeled theoretically using cosmological hydrodynamical simulations of Milky Way–like galaxies. However, capturing the halo-to-halo variation in the local dark matter speeds—a necessary step for quantifying the astrophysical uncertainties that feed into experimental results—requires a sufficiently large sample of simulated galaxies, which has been a challenge. In this Letter, we quantify this variation with nearly 100 Milky Way–like galaxies from the tng50 simulation, the largest sample to date at this resolution. Moreover, we introduce a novel phase-space scaling procedure that endows every system with a reference frame that accurately reproduces the local standard-of-rest speed of our Galaxy, providing a principled way of extrapolating the simulation results to real-world data. The ensemble of predicted speed distributions is well characterized by the standard halo model, a Maxwell-Boltzmann distribution truncated at the escape speed, though the individual distributions can deviate from it, especially at high speeds. The dark matter–nucleon cross section limits placed by these speed distributions vary by ∼ 60% about the median. This places the 1⁢𝜎 astrophysical uncertainty at or below the level of the systematic uncertainty of current ton-scale detectors, even down to the energy threshold. The predicted uncertainty remains unchanged when subselecting on those TNG 50 galaxies with merger histories similar to the Milky Way. Tabulated speed distributions, as well as Maxwell-Boltzmann fits, are provided for use in computing direct detection bounds or projecting sensitivities.

Milky Way↗
DOE OSTI2022

Logic operations with active topological defects

Logic operations performed by semiconductor-based transistors are the basis of modern computing. There is considerable interest in creating autonomous materials systems endowed with the capability to make decisions. In this work, we introduce the concept of using topological defects in active matter to perform logic operations. When an extensile active stress in a nematic liquid crystal is turned on, +1/2 defects can self-propel, in analogy to electron transport under a voltage gradient. By relying on hydrodynamic simulations of active nematics, we demonstrate that patterns of activity, when combined with surfaces imparting certain orientations, can be used to control the formation and transport of +1/2 defects. We further show that asymmetric high- and low-activity patterns can be used to create effective defect gates, tunnels, and amplifiers. The proposed active systems offer the potential to perform computations and transmit information in active soft materials, including actin-, tubulin-, and cell-based systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Predicting fusion ignition at the National Ignition Facility with physics-informed deep learning

Here, an inertial confinement fusion experiment, carried out at the National Ignition Facility, has achieved ignition by generating fusion energy exceeding the laser energy that drove the experiment. Prior to the experiment, a generative machine learning model that combines radiation hydrodynamics simulations, deep learning, experimental data, and Bayesian statistics was used to predict, with a probability greater than 70%, that ignition was the most likely outcome for this shot.

Spears, Brian K. [Lawrence Livermore National Labo↗
DOE OSTI2021

PeleLM

An adaptive mesh hydrodynamics simulation code for low Mach number reacting flows

Esclapez, Lucas↗
DOE OSTI2023

Feature Interpretability

The feature interpretability code is a python module that interprets and analyzes neural networks trained on hydrodynamic simulation output in the form of numpy arrays. The code takes trained neural networks and extracts internal model states in the form of images. Additionally, tools for covariance analysis of network weights and predictions are provided. This code is built on the TensorFlow and PyTorch python libraries, and includes trained networks and example input data for demonstration purposes.

Callis, Skylar↗
DOE OSTI2022

Physics-driven learning of Wasserstein GAN for density reconstruction in dynamic tomography

Object density reconstruction from projections containing scattered radiation and noise is of critical importance in many applications. Existing scatter correction and density reconstruction methods may not provide the high accuracy needed in many applications and can break down in the presence of unmodeled or anomalous scatter and other experimental artifacts. Incorporating machine-learning models could prove beneficial for accurate density reconstruction, particularly in dynamic imaging, where the time evolution of the density fields could be captured by partial differential equations or by learning from hydrodynamics simulations. In this work, we demonstrate the ability of learned deep neural networks to perform artifact removal in noisy density reconstructions, where the noise is imperfectly characterized. Here, we use a Wasserstein generative adversarial network (WGAN), where the generator serves as a denoiser that removes artifacts in densities obtained from traditional reconstruction algorithms. We train the networks from large density time-series datasets, with noise simulated according to parametric random distributions that may mimic noise in experiments. The WGAN is trained with noisy density frames as generator inputs, to match the generator outputs to the distribution of clean densities (time series) from simulations. A supervised loss is also included in the training, which leads to an improved density restoration performance. In addition, we employ physics-based constraints such as mass conservation during the network training and application to further enable highly accurate density reconstructions. Our preliminary numerical results show that the models trained in our frameworks can remove significant portions of unknown noise in density time-series data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

On Hydrodynamics Generated from Electrostatic Discharges: A Heat Source Perspective

A parameter space map of a heat source depositing energy into a cylindrical spark channel modeling the hydrodynamic effects of an electrostatic discharge event is explored. The four parameters of spatial weight, λ R , rise time, τ r , fall time, τ f , and energy per unit length deposited in to the system, E 0 /h, were varied to observe resulting hydrodynamic phenomena. The energy budget con-sequences of the associated hydrodynamic flow are discussed in terms of the remaining electrical energy available for the joule heating of sensitive energetic materials (or “victims”) are discussed,where hydrodynamic simulations are validated by experimental data. Future results and conclusions will be added to this report as they are produced.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

ALEGRA/Sceptre Code Coupling [Brief]

Researchers at Sandia have developed an advanced radiation hydrodynamic simulation capability by coupling the ALEGRA and Sceptre codes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Opacity Table Generation

Radiation opacity is a key property of material systems which describes the absorption and scattering of photons. Such processes are an important part of energy transport, diagnostics, and astrophysical phenomena. Recent experiments in iron opacity and continuum lowering have sparked renewed interest in and careful vetting of opacity models. Even within the local thermodynamic equilibrium approximation, opacity varies over a large parameter space: material, temperature, density, and energy. Opacity tables that span wide regimes of this parameter space are of prime importance to radiation-hydrodynamics codes. Accurate opacity calculations cannot be achieved inline during a radiation hydrodynamics simulation since they involve accounting for all possible electron transitions of all important atomic states of a thermal ensemble. Livermore’s current framework utilizes precalculated tables at many densities, temperatures, and energies. Last year, the opacity theory team completed an L2 milestone concerning the development of a new code, Opus, which served the dual purpose of creating a modern code infrastructure and version control as well as the basis for training new members of that team, a critical need as key members near or are passed retirement. For the current L2 milestone, we exercised and validated Opus on three elements: boron, carbon, and nitrogen, which have few enough electrons to allow a careful convergence study over many input parameters. Here, we report on the production of new tables, show their convergence with respect to important parameters, state some of their limitations, and document a mostly automated production process. In this section, we address the milestone completion criteria, found on the title page and the subsection headings. In section 2, we validate Opus against Tycho and compare with TabOp and Atomic. Section 3 details some of the physics implemented in Opus and relevant to this report. We then cover details of the density-temperature grid in section 4 and do a spectral comparison between different parameter settings and different codes in section 5. After summarizing these results in section 6, we give instructions how to run the opacity table scripts in appendix A and document the scripts and input/output files in appendix B.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On Two-Phase Pressure and Temperature Equilibration with Mie-Grüneisen Equations of State

In the course of reactive flow modeling, the phases of high explosive (HE) reactants and products introduce additional degrees of freedom to the underlying hydrodynamic equations of compressible reactive fluid mechanics. Hydrodynamic simulations track the bulk (mixture) conservation properties of mass, momentum and energy. When multiple phases are introduced, one must decide on appropriate closure rules to determine the relative phase volumes and energy distribution between the phases. One typical closure is to assume that the phases are in pressure and temperature equilibrium. Here, the requisite equations of p-T equilibration are explored for Mie-Grüneisen Equations of State (EOS). Examples with Davis reactants and Davis products EOSs will be given. The Davis forms are used in both the WSD and AWSD reactive flow models.

42 ENGINEERING↗