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 109 records · Page 6

Double cylinder target design for study of hydrodynamic instabilities in multi-shell ICF

Cylindrical implosions are used to study hydrodynamic instability growth for inertial confinement fusion (ICF) applications, as the cylindrical geometry allows for easier diagnostic access while retaining convergence effects. In this work, we use the established cylindrical implosion platform [Palaniyappan et al., Phys. Plasmas 27, 042708 (2020)] to inform the double shell ICF campaign [Montgomery et al., Phys. Plasmas 25, 092706 (2018)]. We present a design for a double cylindrical target as an analogue to the double shell ICF capsule in order to study hydrodynamic instability growth on the high-Z inner shell. Our design work is done with two-dimensional (2D) Eulerian radiation-hydrodynamics simulations, considering the axial uniformity of the implosion and feasibility of measuring the instability growth of pre-seeded single mode sinusoidal perturbations. We discuss in depth the design for a target to be directly driven at the OMEGA laser facility [Boehly et al., Opt. Commun. 133, 495 (1997)]. We evaluate the design for axial implosion symmetry and visibility of instability growth using synthetic radiographs constructed from the simulations, as the instability growth on the inner cylinder is experimentally measured using x-ray radiography of the implosion. We find that the seeded perturbation growth on the inner cylinder should be visible in an experiment, even with axial implosion asymmetry and preheat. We compare our 2D simulations with linear theory predictions for perturbation growth and show that a cylinder with lower azimuthal mode number (mode-20) perturbations compares more favorably with linear theory, while a cylinder with higher azimuthal mode number (mode-40) perturbations at the same starting amplitude saturates and is over-predicted by linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrodynamic expansion and near-infrared absorption of x-ray heated aluminum plasmas

We use x-ray pulses from dense argon plasmas at the Z Machine (Sandia National Laboratories) to generate hypersonic aluminum plasmas akin to material ejecta during proposed planetary defense missions, fusion reactor wall excursions, and other high-energy density processes. Near-infrared absorption is used to diagnose the controlled expansion of the plasmas through cylindrical cavities following their generation from x-ray heating of solid aluminum 7075 alloy. The data are compared to multidimensional radiation hydrodynamics simulations utilizing the ALEGRA multiphysics code, accounting for the dynamics of radiation scattering, material phase change, plasma expansion, thermal re-irradiation, and interactions with the cavity and with the infrared beams. To allow for accurate simulation, density functional theory is used to apply the Hagen–Rubens relation for the far-infrared and is adjoined with a detailed configuration accounting model using the Propaceos code, producing opacities spanning 10 −1 –10 4 eV photon energy for aluminum 7075 alloy, and in comparison with pure aluminum. The model is found to agree with experimental data in the higher-fluence regime when the Hagen–Rubens relation is applied. The ejected material, which is observed to travel up to 55 km/s, is comprised of a strongly ionized, non-LTE plasma front at ∼10 eV temperature followed by a weakly ionized LTE gas at higher density. The present findings lend some confidence to the broad-range equation of state and infrared opacity models for weakly ionized aluminum plasmas while demonstrating an approach to their future refinement, with potential application to astrophysical plasmas and other extreme processes.

Adiabatic process↗

Radiation drive designed to extend the pressure ranges measured in Gbar equation of state experiments at the National Ignition Facility

We present the design and demonstration of a Shock-Strengthening hohlraum radiation temperature drive in the Gbar experimental platform at the National Ignition Facility intended to increase the pressure range measured in a single experiment. Previously published experiments by Döppner et al. measured the equation of state in polystyrene from 25 to 60 Mbar. Recent experimental data of the Shock-Strengthening drive initially demonstrated a much larger pressure range from 15 to 110 Mbar using the same peak radiation temperature and experimental platform. The Shock-Strengthening drive starts with a low temperature foot that launches a weak shock into the sample and is followed by a continuous increase in radiation temperature to strengthen the leading shock. The additional strengthening increases the pressure within the sample beyond what is achievable by convergence alone. Design features of the Shock-Strengthening drive and accompanying radiation hydrodynamics simulations are used to illustrate the method by which the pressure range is increased from previous experiments. This method of modifying the radiation temperature drive can be used on the Gbar platform to significantly increase the range for equation of state data collected in a single experiment for many materials.

Physics - Plasma physics↗

Experiments and simulations of isochorically heated warm dense carbon foam at the Texas Petawatt Laser

An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.

36 MATERIALS SCIENCE↗

Transition in ICF Capsule Implosions

Longstanding design and reproducibility challenges in inertial confinement fusion (ICF) capsule implosion experiments involve recognizing the need for appropriately characterized and modeled three-dimensional initial conditions and high-fidelity simulation capabilities to predict transitional flow approaching turbulence, material mixing characteristics, and late-time quantities of interest—e.g., fusion yield. We build on previous coarse graining simulations of the indirect-drive national ignition facility (NIF) cryogenic capsule N170601 experiment-a precursor of N221205 which resulted in net energy gain. We apply effectively combined initialization aspects and multiphysics coupling in conjunction with newly available hydrodynamics simulation methods, including directional unsplit algorithms and low Mach-number correction-key advances enabling high fidelity coarse grained simulations of radiation-hydrodynamics driven transition. Our presentation includes discussion of the capsule initialization and implosion dynamics, analysis of the vorticity production budget, transition signatures, quantities of interest—late-time ion temperature and fusion-neutron yield, numerical uncertainty quantification, and comparisons with NIF data.

42 ENGINEERING↗

Modeling radiative-shocks created by laser–cluster interactions

Radiative-shocks induced by laser–cluster interactions are modeled using radiation-hydrodynamic simulations. A good agreement—in both shock velocity and density profiles—is obtained between experiment and simulations, indicating that non-local thermodynamic equilibrium (NLTE) radiative effects are important in the experimental regime examined, particularly at early times (≤30 ns) due to the elevated temperatures (≥35 eV). The enhanced NLTE radiative emission causes the shock to be reduced in amplitude, increased in width, and reduced in propagation velocity, while the amplitude of the radiative precursor is increased. As the density and temperature conditions are relatively modest, this potentially has important implications for the scalings that are used in laboratory–astrophysics to transform between laboratory and astrophysical scales, which do not hold for non-LTE systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced laser-energy coupling with small-spot distributed phase plates (SG5-650) in OMEGA DT cryogenic target implosions

Cryogenic deuterium–tritium ice target implosions on OMEGA with new small-spot (SG5-650) distributed phase plates (DPPs) achieved an (11 ± 4)% increase in energy coupling compared to implosions with larger-spot SG5-850 DPPs by decreasing the ratio of the laser spot diameter to the target diameter from 0.93 to 0.75. The SG5-650 DPPs provide a focus spot size of 674 μm, which is defined as the diameter that encircles 95% of the measured beam energy compared to 834 μm for the SG5-850, which are the standard DPPs in cryogenic target implosions on OMEGA. The hydrodynamic efficiency, defined as the ratio of the kinetic energy in the imploding shell to the laser energy, increased from 4.5% to 5.0% based on radiation-hydrodynamic calculations benchmarked to shell trajectory and bang-time measurements. Furthermore, the higher coupling came with a trade-off of an increased hot-electron production as well as increased hydrodynamic instabilities seeded by a larger mode-10 amplitude from the beam port geometry, both of which may have reduced the fusion neutron production and areal density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhanced Laser–Energy Coupling with Small-Spot Distributed Phase Plates (SG5-650) in OMEGA DT Cryogenic Target Implosions

Cryogenic deuterium–tritium ice target implosions on OMEGA with new small-spot ("SG5-650") distributed phase plates (DPP's) achieved an (11 +/- 4)% increase in energy coupling compared to implosions with standard-spot DPP's by decreasing the ratio of the laser spot diameter to the target diameter from 0.93 to 0.75. The SG5-650 DPP's provide a focus spot size of 674 um, which is defined as the diameter that encircles 95% of the measured beam energy compared to 834 um for the SG5-850. The hydrodynamic effciency, defined as the ratio of the kinetic energy in the imploding shell to the laser energy, increased from 4.5% to 5.0% based on radiation-hydrodynamic calculations benchmarked to shell trajectory and bang-time measurements. The higher coupling came with a trade-off of an increased hot-electron production as well as increased hydrodynamic instabilities seeded by a larger mode-10 amplitude from the beam port geometry, both of which may have affected the fusion neutron production and areal density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct Experimental Proof of the Principal Role of Reduced High-Mode Hydrodynamic Mix in Recent Ignition Success on NIF

This Letter details the role of reduced hydrodynamic mix in achieving the first fusion plasma to satisfy the Lawson criterion for ignition. Using novel 3D measurements of the temperature and density of the burning deuterium-tritium plasma across a series of experiments, we demonstrate that hydrodynamic mixing from small-scale capsule defects was a key degradation mechanism inhibiting ignition. A series of five layered deuterium-tritium experiments using 1.9 MJ of laser energy demonstrates a 2 × reduction of the thermonuclear yield when a lower quality capsule, typical of preignition experiments, is used. Radiation hydrodynamic calculations consistent with the observed mix profile show that the mixed fuel’s reduced compression and increased radiative loss explain the yield reduction. Furthermore, a hydroscaling argument suggests that a ∼ 3 MJ laser would be needed to reach ignition with such a capsule, well beyond the current capabilities of high energy laser systems. As such, this Letter gives critical insight into the physics of mix in self-heated fusion plasmas and establishes a foundation to evaluate capsule quality needs of future inertial fusion designs that aim at ignition and high gain.

High-energy-density plasmas↗

X-ray heating and electron temperature of laboratory photoionized plasmas

We discuss experimental and modeling results for the x-ray heating and temperature of laboratory photoionized plasmas. A method is used to extract electron temperature based on the analysis of transmission spectroscopy data that is independent of atomic kinetics modeling. Boltzmann kinetic and radiation-hydrodynamics simulations emphasized the critical interplay between atomic physics and heating in photoionized plasmas. They also demonstrated the dramatic impact of photoexcitation on excited state populations, line emissivity, and radiation cooling. Temperature calculations performed with astrophysical codes significantly overestimated the measurement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

National Opacity Program: Analysis of Opacity-relevant X-ray Emission Spectra (Milestone Report ID 7121)

We have produced high energy density iron plasmas at temperatures above 1 keV and electron densities exceeding 1023 cm -3 (~ 1 g/cm 3 ) using the Orion laser at the Atomic Weapons Establishment. These plasmas were created by irradiating 50 µm diameter layered targets with frequency doubled (λ = 527 nm), 1 ps laser pulses focused to a 100 µm diameter producing an irradiance of ~ 2 x 10 18 W/cm 2 . The buried layer targets consist of 160 nm iron sulfide (FeS), 60 nm potassium chloride (KCl), and 15 nm carbon. The combined layers are tamped on both sides with 3 µm of parylene-N. The x-ray emission from the plasma was measured using two time-resolved, and four time-integrated Bragg crystal spectrometers, as well as one time-integrated imaging system. One time-resolved x-ray spectrometer measured emission from L-shell transitions in highly charged iron, the other from K-shell transitions in helium-like S 14+ and hydrogen-like S 15+ . The time-integrated spectrometers are intensity-calibrated and measured emission from K-shell transitions in sulfur, potassium, chlorine, and both K-shell and L-shell transitions in iron. The density and temperature of the plasma were determined by modeling the x-ray spectra using different spectral and hydrodynamic modeling packages. A brief overview of the uncertainties associated with the measurements and models are presented. We also give an overview of our 1-D HYDRA-DCA radiation-hydrodynamics model and improvements for future work. Our results aid in assessing experimental uncertainties associated with plasma uniformity and with x-ray emission employed as diagnostics in opacity experiments at temperatures and densities not achievable elsewhere and represent a significant step in creating and diagnosing plasmas near LTE. These results are summarized as part of the completion requirements for milestone 7121

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutron backscatter edges as a diagnostic of burn propagation

High gain in hotspot-ignition inertial confinement fusion (ICF) implosions requires the propagation of thermonuclear burn from a central hotspot to the surrounding cold dense fuel. As ICF experiments enter the burning plasma regime, diagnostic signatures of burn propagation must be identified. In previous work [A. J. Crilly et al., Phys. Plasmas 27(1), 012701 (2020)], it has been shown that the spectral shape of the neutron backscatter edges is sensitive to the dense fuel hydrodynamic conditions. The backscatter edges are prominent features in the ICF neutron spectrum produced by the 180° scattering of primary deuterium–tritium fusion neutrons from ions. In this work, synthetic neutron spectra from radiation-hydrodynamics simulations of burning ICF implosions are used to assess the backscatter edge analysis in a propagating burn regime. Significant changes to the edge's spectral shape are observed as the degree of burn increases, and a simplified analysis is developed to infer scatter-averaged fluid velocity and temperature. The backscatter analysis offers direct measurement of the increased dense fuel temperatures that result from burn propagation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of Self-shielding Minihalos on the Lyα Forest at High Redshift

Abstract Dense gas in minihalos with masses of 10 6 −10 8 M ⊙ can shield themselves from reionization for ∼100 Myr after being exposed to the UV background. These self-shielded systems, often unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Ly α forest. We first derive the H i column density profile of photoevaporating minihalos by conducting 1D radiation–hydrodynamics simulations. We utilize these results to estimate the Ly α opacity from minihalos in a large-scale simulation that cannot resolve self-shielding. When the ionization rate of the background radiation is 0.03 × 10 −12 s −1 , as expected near the end of reionization at z ∼ 5.5, we find that the incidence rate of damped Ly α absorbers increases by a factor of ∼2−4 compared to at z = 4.5. The Ly α flux is, on average, suppressed by ∼3% of its mean due to minihalos. The absorption features enhance the 1D power spectrum up to ∼5% at k ∼ 0.1 h Mpc −1 (or 10 −3 km −1 s), which is comparable to the enhancement caused by inhomogeneous reionization. The flux is particularly suppressed in the vicinity of large halos along the line-of-sight direction at separations of up to 10 h −1 Mpc at r ⊥ ≲ 2 h −1 Mpc. However, these effects become much smaller for higher ionizing rates (≳0.3 × 10 −12 s −1 ) expected in the post-reionization Universe. Our findings highlight the need to consider minihalo absorption when interpreting the Ly α forest at z ≳ 5.5. Moreover, the sensitivity of these quantities to the ionizing background intensity can be exploited to constrain the intensity itself.

79 ASTRONOMY AND ASTROPHYSICS↗

Beryllium–tungsten graded density inner shells in double shell capsules for improved hydrodynamic stability

The outer surface of the high-Z inner shell in the double shell configuration of inertial confinement fusion experiments experiences Rayleigh–Taylor instability growth during the implosion process due to inverted density and pressure gradients between a highly compressed foam interstitial layer and the accelerating dense inner shell. Graded density layers have long been known to reduce instability growth rates. In this study, we employ high-fidelity radiation hydrodynamic simulations to demonstrate this improved stability when grading beryllium into tungsten. We first characterize the response to L-band preheat of these layers using a newly calibrated radiation drive. While graded layer capsules suffer reduced performance (here, measured as DD neutron yield from a CD foam fuel) in 1D simulations due to reduced kinetic energy coupling and reduced fuel compression, they suffer less of a performance drop when 2D instabilities are accounted for. With the improved stability of graded layers, we explore the performance of capsules with larger fuel radii and thinner shells as a preliminary study to find new designs in which graded layers produce the highest yields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Go with the Radflow: Thermal radiation experiments in the high energy-density regime [Slides]

Los Alamos National Laboratory is a nuclear weapons laboratory supporting our nation's defense. In support of the mission is a high energy-density physics program in which we design and execute experiments to study radiation-hydrodynamics phenomena and improve the predictive capability of our large-scale multi-physics software codes on our supercomputers. The Radflow Project is maturing a unique spectroscopic measurement that yields a spatially dependent supersonic radiation wavefront profile. The spectroscopic measurement is much more constraining than the typically used radiography, which measures the density changes due to the shock as the radiation wave goes subsonic. Recently, the spectroscopic measurement has been applied to advanced targets with a detailed heterogeneity representing a single realization of a stochastic medium. The Radflow Project also conducts opacity experiments to explain the current conundrum of the experimental iron opacity differing from theory, which is important in modeling the sun. More than that, though, we need to be able to answer these science questions and be able to validate our theory and codes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A temperature profile diagnostic for radiation waves on OMEGA-60

Predicting and matching radiation wave propagation with computational models has proven difficult. Information provided by experiments studying radiation flow has been limited when only radiation breakout is measured. We have developed the COAX (co-axial) diagnostic platform to provide spatial temperature profiles of a radiation wave through low density foams as a more detailed constraint for simulations. COAX uses a standard, laser-driven OMEGA-60 halfraum to drive radiation down a titanium-laden silicon oxide foam. Point-projection X-ray absorption spectroscopy perpendicular to the radiation flow measures the spatial profile of titanium ionization. The spectroscopic measurement utilizes a broadband capsule backlighter. Imaging and streak spectroscopy are used to characterize the size and spectrum of this source. Radiography provides an additional constraint by capturing the developing shock as the radiation flow becomes subsonic. The DANTE diagnostic is used to measure the halfraum temperature. Here, we provide a spectroscopic analysis of COAX data to determine temperature, and we describe experimental sources of uncertainty. The temperature is obtained by comparison to multi-temperature synthetic spectra post-processed from radiation-hydrodynamics simulations. Quantitative comparison between data and synthetic spectra generated from temperature profiles at relevant simulation times enable determination of a peak temperature of 114 ± 8 eV at 265 ± 22.4 μm from the halfraum. This represents an improvement over the temperature uncertainties of previous radiation flow experiments. Further refinements to the spectroscopic analysis could achieve ±4 eV. The combination between space-resolved spectroscopy and radiography enables us to determine the distance from the halfraum of both the radiation front and the shock front at the time of measurement. For the example shown in this paper the radiation front position is 600–630 μm at 3.43 ± 0.16 ns and the shock front position is 633 μm at 3.3 ± 0.24 ns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Moment-based adaptive time integration for thermal radiation transport

Here, in this paper we develop a framework for moment-based adaptive time integration of deterministic multifrequency thermal radiation transpot (TRT). We generalize our recent semi-implicit-explicit (IMEX) integration framework for gray TRT to multifrequency TRT, and also introduce a semi-implicit variation that facilitates higher-order integration of TRT, where each stage is implicit in all components except opacities. To appeal to the broad literature on adaptivity with Runge–Kutta methods, we derive new embedded methods for four asymptotic preserving IMEX Runge–Kutta schemes we have found to be robust in our previous work on TRT and radiation hydrodynamics. We then use a moment-based high-order-low-order representation of the transport equations. Due to the high dimensionality, memory is always a concern in simulating TRT. We form error estimates and adaptivity in time purely based on temperature and radiation energy, for a trivial overhead in computational cost and memory usage compared with the base second order integrators. We then test the adaptivity in time on the tophat and Larsen problem, demonstrating the ability of the adaptive algorithm to naturally vary the timestep across 4–5 orders of magnitude, ranging from the dynamical timescales of the streaming regime to the thick diffusion limit.

97 MATHEMATICS AND COMPUTING↗