Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiation Hydrodynamics”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission

Abstract In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor or other means of deflection proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's X-ray energy deposition into the asteroid and the resulting material ablation. The relevant physics in these simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed X-ray energy deposition model developed using the radiation-hydrodynamics code Kull that can be used to initiate a nuclear mitigation mission calculation. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials at a range of porosities, two different source spectra, and a broad range of radiation fluences, source durations, and angles of incidence. Using blowoff momentum as the primary metric, the model-initiated simulation results match the full radiation-hydrodynamics results to within 10%.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Verification and Validation of Two Hydrodynamic Methods for Simulations of High Energy Density Physics Problems

A 3D verification and validation suite of test problems is presented and used to evaluate hydrodynamic methods within a radiation hydrodynamics code, xRAGE. These test problems exercise different levels of complexity, building towards ICF problems which in addition to hydrodynamics also include three temperature plasma physics, thermal conduction, and radiation diffusion. Among the problems in the test suite are the Kidder ball problem, the Verney shell problem, and a 5-material compression problem, which exercise different purely hydrodynamic methods implemented within xRAGE. There is excellent agreement between 2D and 3D XRAGE simulation results and between the xRAGE results and the benchmark solutions. Two 3D ICF test problems are also presented, based on an OMEGA direct drive capsule experiment and on a NIF indirect drive capsule experiment. It is demonstrated that the newer unsplit hydrodynamic method in xRAGE produces more vorticity relative to the older default method. For the indirect drive capsule, the 3D simulations are in reasonable agreement with the experimental values of ion temperature and neutron production.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analytic insights into nonlocal energy transport: Steady state Fokker Planck theory in arbitrary Z plasmas

The generation of energetic electrons in laser fusion in an important issue. The electrons may either arise from a laser plasma instability, or from the uncoupled high temperature tail of a Maxwellian distribution. To study these in a laser fusion context, it is important to find a method accurate enough to be useful, and simple enough to be incorporated into a radiation hydrodynamics numerical simulation, the main workhorse for studying the laser fusion target. That is why analytic insights become important, they allow one to simplify the Fokker Planck theory so that a solution of it can be incorporated into a radiation hydrodynamic simulation. This work develops and analyzes a steady state Fokker Planck theory for plasmas of arbitrary Z. It develops a method of solving the simplified Fokker Planck method with a technique called sparse eigenfunction analysis. As a result, this method appears to work reasonably well when compared to the experimental results from the Rochester/NIF on plastic spherical targets with and without a silicon layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Species Separation and Hydrogen Streaming upon Shock Release from Polystyrene under Inertial Confinement Fusion Conditions

Shock release from inertial confinement fusion (ICF) shells poses a great challenge to single-fluid hydrodynamic equations, especially for describing materials composed of different ion species. This has been evidenced by a recent experiment [Haberberger et al., Phys. Rev. Lett. 123, 235001 (2019)], in which low-density plasmas ( 10 19 to 10 20 cm – 3 ) are measured to move far ahead of what radiation-hydrodynamic simulations predict. To understand such experimental observations, in this work we have performed large-scale nonequilibrium molecular-dynamics simulations of shock release in polystyrene (CH) at experimental conditions. These simulations revealed that upon shock releasing from the back surface of a CH foil, hydrogen can stream out of the bulk of the foil due to its mass being lighter than carbon. This released hydrogen, exhibiting a much broader velocity distribution than carbon, forms low-density plasmas moving in nearly constant velocities ahead of the in-flight shell, which is in quantitative agreement with the experimental measurements. Such kinetic effect of species separation is currently missing in single-fluid radiation-hydrodynamics codes for ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing

The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

NED_Asteroid_Energy_Deposition

In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's x-ray energy deposition into the asteroid and the resulting material ablation. These simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed x-ray energy deposition model developed using the radiation-hydrodynamics code Kull which can be used to initiate a nuclear mitigation mission hydrocode. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials (Silicon Dioxide, Forsterite, Iron, and Ice), two different source spectra (1 and 2 keV blackbodies), and then a broad range of radiation fluences (0.0001 to 1 kt per square meter), source durations (10 to 100 ns), and asteroid porosities (0 to 80 percent). Using blowoff momentum as the primary metric, the modelinitiated simulation results match the full radiation-hydrodynamics results to within 10 percent. Please reference the journal article: Burkey et al., X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission, Planetary Science Journal, (2023) for more information.

Managan, RobertA↗
DOE OSTI2022

First measurement of the 10B(α,n)13N reaction in an inertial confinement fusion implosion at the National Ignition Facility: Initial steps toward the development of a radiochemistry mix diagnostic

We report the first measurement of the 10B(α,n)13N reaction in a polar-direct-drive exploding pusher (PDXP) at the National Ignition Facility (NIF). This work is motivated by the need to develop alternative mix diagnostics, radiochemistry being the focus here. The target is composed of a 65/35 at. % deuterium–tritium (DT) fill surrounded by a roughly 30 μm thick beryllium ablator. The inner portion of the beryllium ablator is doped with 10 at. % of 10B. Radiation-hydrodynamics calculations were performed in 1D to optimize both the remaining boron rho-R and the DT neutron yield. A charged-particle transport post-processor has been developed to study α-induced reactions on the ablator material. Results indicate a large 13N production from α-induced reactions on 10B, measurable by the radiochemical analysis of a gaseous samples system at the NIF. The PDXP target N201115-001 was successfully fielded on the NIF, and nitrogen from the 10B(α,n)13N reaction was measured. The 13N production yield, as well as the DT neutron yield, was, however, lower than expected. Some of the reduced yields can be explained by the oblate shape, but the ratios of the various radiochemical signals are not commensurate with expectations based on a simple reduction of the 1D results. Preliminary 2D radiation-hydrodynamics computations are consistent with the experimental measurements, and work is ongoing to extend the radiochemistry analysis into higher dimensions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An accelerated approach to inline non-LTE modeling

A goal across many physics and engineering disciplines is to use radiation-hydrodynamics codes to accurately simulate systems that are expected to deviate from local thermodynamic equilibrium (LTE), such as laser-driven hohlraums, stellar coronae, and supernova ejecta. The computational expense of traditional non-LTE (NLTE) opacity calculations can dominate the cost of simulating these systems and render NLTE physics untenable in most cases. In this paper, we present a new parallel implementation for NLTE opacity calculations that employs advanced load balancing, vectorization, and memory management to significantly reduce the computational costs of these calculations. Using the radiation-hydrodynamics code xRAGE , we demonstrate this new implementation on three test problems and present order-of-magnitude speed improvements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Design of Cylindrical Implosion Experiments to Demonstrate Scale-Invariant Rayleigh-Taylor Instability Growth

Radiation-hydrodynamics simulations are used to design laser-driven cylindrical implosion experiments to directly measure hydrodynamic instability growth in convergent geometry. Designs for two different size targets, varying in radial dimension by a factor of three, are presented. A set of beam pointings and powers are identified for each scale design that result in a nearly axially uniform implosion of an embedded marker layer. The implosion trajectories are shown to be scale-invariant between designs, with nearly identical scaled acceleration profiles. Finally, linear theory and radiation-hydrodynamics simulations predict that Rayleigh-Taylor instability growth of an azimuthal perturbation, machined on the inner surface of the embedded marker, is also scale-invariant between designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

X-ray phase-contrast imaging of strong shocks on OMEGA EP

The ongoing improvement in laser technology and target fabrication is opening new possibilities for diagnostic development. An example is x-ray phase-contrast imaging (XPCI), which serves as an advanced x-ray imaging diagnostic in laser-driven experiments. In this work, we present the results of the XPCI platform that was developed at the OMEGA EP Laser-Facility to study multi-Mbar single and double shocks produced using a kilojoule laser driver. Two-dimensional radiation-hydrodynamic simulations agree well with the shock progression and the spherical curvature of the shock fronts. It is demonstrated that XPCI is an excellent method to determine with high accuracy the front position of a trailing shock wave propagating through an expanding CH plasma that was heated by a precursor Mbar shock wave. Furthermore, the interaction between the rarefaction wave and the shock wave results in a clear signature in the radiograph that is well reproduced by radiation-hydrodynamic simulations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Characterization of thermal transport and evolution of Au plasma in ICF experiments by Thomson scattering

Here, this paper demonstrates the capability of optical Thomson scattering (OTS) to measure thermal transport, and provides support to radiation hydrodynamic and kinetic simulations of electron thermal transport and plasma evolution. OTS theory and plasma simulations are applied to the interpretation of experimental measurements of laser-produced plasma from spherical gold targets on the OMEGA laser facility. The dynamical form factor, S($\vec{k}$, ω), of electron density fluctuations that is used in the fitting of Thomson scattering spectra includes ion–ion collisions and effects of non-Maxwellian distribution functions. OTS measurements and their interpretation are consistent with the nonlocal transport model in radiation hydrodynamic simulations as well as with kinetic simulations in the second half of the probe pulse duration. In particular, the reversal of heat transport during cooling is observed to be consistent with simulations, while some discrepancies are noted during the initial heating of the Au targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Three-dimensional reconstruction of inertial confinement fusion hot-spot plasma from x-ray and nuclear diagnostics on OMEGA

Multidimensional effects degrade the neutron yield and the compressed areal density of laser-direct-drive inertial confinement fusion implosions of layered deuterium–tritium cryogenic targets on the OMEGA Laser System with respect to 1D radiation-hydrodynamic simulation predictions. A comprehensive physics-informed 3D reconstruction effort is under way to infer hot-spot and shell conditions at stagnation from four x-ray and seven neutron detectors distributed around the OMEGA target chamber. Neutron diagnostics, providing measurements of the neutron yield, hot-spot flow velocity, and apparent ion-temperature distribution, are used to infer the mode-1 perturbation at stagnation. The x-ray imagers record the shape of the hot-spot plasma to diagnose mode-1 and mode-2 perturbations. A deep-learning convolutional neural network trained on an extensive set of 3D radiation-hydrodynamic simulations is used to interpret the x-ray and nuclear measurements to infer the 3D profiles of the hot-spot plasma conditions and the amount of laser energy coupled to the hot-spot plasma. A 3D simulation database shows that larger mode-1 asymmetries are correlated with higher hot-spot flow velocities and reduced laser-energy coupling and neutron yield. Three-dimensional hot-spot reconstructions from x-ray measurements indicate that higher amounts of residual kinetic energy are correlated with higher measured hot-spot flow velocities, consistent with 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Charged particle transport coefficient challenges in high energy density plasmas

High energy density physics (HEDP) and inertial confinement fusion (ICF) research typically relies on computational modeling using radiation-hydrodynamics codes in order to design experiments and understand their results. These tools, in turn, rely on numerous charged particle transport and relaxation coefficients to account for laser energy absorption, viscous dissipation, mass transport, thermal conduction, electrical conduction, non-local ion (including charged fusion product) transport, non-local electron transport, magnetohydrodynamics, multi-ion-species thermalization, and electron-ion equilibration. In many situations, these coefficients couple to other physics, such as imposed or self-generated magnetic fields. Furthermore, how these coefficients combine are sensitive to plasma conditions as well as how materials are distributed within a computational cell. Uncertainties in these coefficients and how they couple to other physics could explain many of the discrepancies between simulation predictions and experimental results that persist in even the most detailed calculations. This paper reviews the challenges faced by radiation-hydrodynamics in predicting the results of HEDP and ICF experiments with regard to these and other physics models typically included in simulation codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Long-lived neutron-star remnants from asymmetric binary neutron star mergers: element formation, kilonova signals and gravitational waves

We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ∼10 −3 –10 −2 M ⊙ followed by a persistent spiral-wave/neutrino-driven wind of ≳ 10 −2 M ⊙ with material spanning a wide range of electron fractions, ∼0.1–0.55. Dynamical ejecta (winds) have fast velocity tails up to ∼0.8 (∼0.4) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete r-process yields and identify the production of 56 Ni and the subsequent decay chain to 56 Co and 56 Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from r-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak r-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.

abundances↗