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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Studies of localized mix and asymmetry caused by the stalk in inertial fusion experiments at OMEGA

In an inertially confined fusion experiment, a capsule of fusion fuel is compressed to high temperatures and densities using lasers, creating a hot plasma in which thermonuclear fusion reactions take place. To achieve efficient conversion of laser energy to fuel heating, the implosion must be as spherically symmetric as possible. One common source of asymmetry that degrades performance is the stalk and glue spot that hold the capsule in place before the experiment. The exact nature and magnitude of this perturbation had not been previously quantified. In this work, inertially confined fusion experiments were executed with an additional stalk on the capsule to measure and quantify the effect of the stalk and glue spot on the symmetry and performance. It is found that the stalk drives low-mode asymmetry by shadowing the laser drive on the part of the shell near the glue spot, enhancing the hot-spot velocity toward the stalk by 40–80 km/s and elongating the hot-spot along the stalk axis. It is also found that localized mix due to the additional stalk has a negligible effect on implosion performance. Two-dimensional radiation-hydrodynamic simulations using the xRAGE code reproduce the direction of the experimentally observed change to hot-spot velocity and implosion performance metrics, but not the observed elongation, suggesting that this feature results from an interaction between stalk perturbations and other preexisting asymmetry seeds, such as laser drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Wetted-foam target stability for laser-direct-drive inertial confinement fusion

Using highly resolved three-dimensional radiation-hydrodynamic simulations, the effects of wetted-foam (WF) ablators in direct-drive inertial confinement fusion cryogenic implosions on the OMEGA laser were investigated. The performance of these implosions was compared with that of implosions using targets with solid-plastic ablators. Simulations considering WF ablators as a homogeneous material show no mitigation of the effects of laser imprint, where the absence of this mitigation can be attributed to the use of picket pulses. Heterogeneous random-structure WF ablators introduce large short-wavelength perturbations, which are seeded by WF inhomogeneities. These perturbations exceed the perturbations seeded by the imprint and cause an additional reduction in implosion performance when the average pore size in the WF exceeds ~0.2 μm. Foam-seeded perturbations are expected to be less damaging in more stable high-adiabat (⁠ ≳ 5⁠) implosions. Here, the perturbations introduced by the WF ablators can be substantially reduced by using a foam matrix with a layered structure instead of a random one. Further study of the effects of WFs in ignition designs of the MJ laser energy scale is required in the perspective of future inertial fusion energy mass production.

Igumenshchev, I. V. [University of Rochester, NY (↗
DOE OSTI2026

Robustness of inertial fusion energy relevant implosions to low-mode asymmetries

Inertial confinement fusion implosions for power production need to be robust. Degradation sources need to be understood and considered for an economic power plant. In this work, the impact of low-mode asymmetries on yield is studied using 2D radiation hydrodynamics simulations for inertial fusion energy relevant implosions. First, a high-gain indirectly driven capsule design is developed. This design is subsequently degraded with low-mode drive asymmetries. The results show that yield is constant with increasing amounts of asymmetry until a cliff is reached where the yield quickly drops below ignition. Robustness arises because the degraded hotspot still has enough time to get to ignition before capsule decompression. This robustness occurs even though the asymmetry couples less energy to the hotspot and causes decompression to occur sooner. The robustness is linked to the generalized Lawson parameter, χ⁠. The value of χ steadily decreases with increasing asymmetry. Once χ is close to one, the hotspot is no longer able to ignite before the capsule decompresses. A simple piston model is used to predict where the robustness cliff is for any design. The trade-off between yield and robustness is also demonstrated.

Physics - Plasma physics↗
DOE OSTI2022

Combinatorial sputter deposition of ultrathick Au-Bi alloy films

We report gold-bismuth alloys are of interest as catalysts and catalytic sensing systems, electrochemical sensors, superconductors, and hohlraums for magnetically assisted inertial confinement fusion implosions. Radiation-hydrodynamics simulations with the Lasnex code of laser-driven hohlraums predict higher x-ray drive from Au-Bi alloys compared with cases of Au-Ta or pure Au and Bi hohlraums. Here, we use direct current magnetron sputtering in Ar gas, with co-sputtering from two elemental targets, to deposit Au-Bi alloys with Bi content of 9–77 at.% and thicknesses up to ~20 µm. Films are characterized by a combination of x-ray diffraction, Rutherford backscattering, scanning electron microscopy, substrate-curvature-based residual stress, and electronic transport measurements. Experiments are complemented by Monte Carlo simulations of ballistic sputtering and gas phase transport of depositing species and Ar gas atoms. Results show that all films are polycrystalline, with three distinct compositional regimes dominated by Au, Au 2 Bi, and Bi crystallographic phases. A metallic behavior of the temperature dependence of electrical resistivity is observed for all the films. Films with Bi content above ~30 at.% exhibit porosity, which is tolerable to hohlraum x-ray drive based on Lasnex simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulation studies of the interaction of laser radiation with additively manufactured foams

The interaction of laser radiation with foams of various porosities and low densities has been the subject of several numerical and experimental studies. In all cases, the modeling of low-Z under-dense foams as uniform gases of equivalent average density using standard radiation-hydrodynamics codes has resulted in heat-front velocities that are considerably faster than those observed experimentally. It has been theoretically conjectured that this difference may be attributed to the breakdown of the foam’s morphology, leading to a dynamics of filament expansion where the ion and electron energy partitions are significantly different from those calculated using the uniform gas model. We found that 3D computer simulations employing a disconnected representation of the foam’s microstructure which allowed for the dynamics of foam element heating, expansion, and stagnation largely supported the theoretical picture. Simulations using this model for laser experiments on under-dense 2 mg cc -1 SiO 2 aerogel foams reproduced the experimental data fairly well. We used the validated model in simulations of low-density structured foam-like materials (produced via additive manufacturing) with a variety of morphologies. We found that the log-pile configurations were consistent with the analytical propagation model of Gus’kov et al. Further validation of the model was obtained by simulating experiments performed at the Jupiter Laser Facility using the log-pile and octet-truss foam morphologies. Simulations of the foam-laser interaction using a wave propagation code showed that the microstructure was able to enhance stimulated Brillouin scattering (SBS) by concentrating the light energy into density holes. In turn, this promotes laser filamentation, reducing SBS and bringing the predicted values closer to the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

3D simulations of inertial confinement fusion implosions part 2: systematic flow anomalies and impact of low modes on performances in OMEGA experiments

We present the first 3D radiation-hydrodynamic simulations of directly driven inertial confinement fusion implosions with an inline package for polarized crossed-beam energy transfer, which were used to assess the impact of the current distributed polarization rotators (DPRs) on OMEGA as well as other known sources of asymmetry. Applied to OMEGA implosions, the simulations predict bang times with no need for ad hoc multipliers, as well as yields—if you separately account for the impacts of imprint and fuel age. The magnitude of the flow is well reproduced when the low mode sources are large, whereas the modeling of the stalk is thought to be required to match the flow magnitude in the remaining cases. For the cases explored in more detail, polarized cross-beam energy transfer (CBET)—the only known systematic drive asymmetry, brought the results closest to the measured flow vectors. The remaining discrepancies are shown to be stemming from the limited knowledge of the laser pointing modes. For typical current levels of beam mispointing, power imbalance, target offset, and asymmetry caused by polarized CBET, low modes degrade the yield by more than 40%. The current strategy of attempting to compensate the mode-1 asymmetry with a preimposed target offset recovers only about one-third of the losses caused by the low modes due to the dynamic nature of the multiple asymmetries and the presence of low modes other than l = 1. Furthermore, addressing the root causes of the drive asymmetries is apt to be more beneficial. To that end, one possible solution to the specific issue of polarized CBET (10 µm DPRs) is shown to work well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

3D simulations of inertial confinement fusion implosions part 1: inline modeling of polarized cross beam energy transfer and subsequent drive anomalies on OMEGA and NIF

Inertial confinement fusion experiments are sensitive to cross-beam energy transfer (CBET), a laser-plasma instability that redistributes laser energy in the coronal plasma through self-generated ion acoustic wave (IAW) gratings. The detailed CBET coupling depends on the polarization state of the crossing wavefields. CBET itself can also scramble the beam polarizations by inducing ellipticity through the IAW grating, and rotating the seed polarization toward that of the pump. We develop a ray-based model that describes the polarized CBET coupling and that is compatible with the framework of 3D inline radiative hydrodynamics simulations. The model is implemented in the ASTER/IFRIIT code and verified against an academic test case and an offline polarized CBET post-processor. It is then applied to the detailed configuration of the distributed polarization rotator system on OMEGA, where results highlight how polarized CBET induces significant low modes in the collisional absorption source term. Finally, the modeling is applied to a simple indirect-drive configuration, comparing CBET calculations with 96 unpolarized or polarized beams with 24 unpolarized quads. It is shown that these cases produce similar power amplification per cone of beams grouped with similar polar angles. However, the 96 beam geometry itself is found to reduce azimuthal variations in quad power after the interaction and favors beams with larger polar angles within the cones, an effect that is amplified by the polarized CBET. Application of the model to inline calculations of OMEGA implosions are presented in a companion paper.

3D inline modeling↗
DOE OSTI2021

Lyman-α transmission properties of the intergalactic medium in the CoDaII simulation

ABSTRACT The decline in abundance of Lyman-α (Lyα) emitting galaxies at z ≳ 6 is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of ∼94 comoving Mpc side length, to compute the Lyα transmission properties of the intergalactic medium (IGM) at z ∼ 5.8 to 7. Our results mainly confirm previous studies, i.e. we find a declining Lyα transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Lyα peaks at high redshift, we also analyse the IGM transmission on the blue side, which shows a rapid decline at z ≳ 6 of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\, \mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 1{{\ \rm per\ cent}}$ of sightlines towards M1600AB ∼ −21 galaxies at z ∼ 7 are transparent enough to allow a transmission of a blue Lyα peak. We discuss our results in the context of the interpretation of observations.

Gronke, Max (ORCID:000000032491060X)↗
DOE OSTI2021

The collapse and three-dimensional explosion of three-dimensional massive-star supernova progenitor models

ABSTRACT The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M⊙ model, evolved in three dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code fornax. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multidimensional initial models serves as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

r -process nucleosynthesis and kilonovae from hypermassive neutron star post-merger remnants

ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The force explosion condition is consistent with spherically symmetric CCSN explosions

ABSTRACT One of the major challenges in core-collapse supernova (CCSN) theory is to predict which stars explode and which collapse to black holes. The analytic force explosion condition (FEC) shows promise in predicting which stars explode in that the FEC is consistent with CCSN simulations that use the light-bulb approximation for neutrino heating and cooling. In this follow-up manuscript, we take the next step and show that the FEC is consistent with the explosion condition when using actual neutrino transport in gr1d simulations. Since most 1D simulations do not explode, to facilitate this test, we enhance the heating efficiency within the gain region. To compare the analytic FEC and radiation-hydrodynamic simulations, this manuscript also presents a practical translation of the physical parameters. For example: we replace the neutrino power deposited in the gain region, Lντg, with the net neutrino heating in the gain region; rather than assuming that $\dot{M}$ is the same everywhere, we calculate $\dot{M}$ within the gain region; and we use the neutrino opacity at the gain radius. With small, yet practical modifications, we show that the FEC predicts the explosion conditions in spherically symmetric CCSN simulations that use neutrino transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Cosmic variance and the inhomogeneous UV luminosity function of galaxies during reionization

When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from z ≃ 6–15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same z. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by z = 6 (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, ‘rolling-over’ at M UV ≳ -17. CoDa II’s UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at z = 6–8, including the faint end, except for the faintest data point at z = 6, based on one galaxy at M UV = -12.5. According to CoDa II, the probability of observing the latter is $\sim 5~{{\ \rm per\ cent}}$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. Here, we find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2–4, primarily by underestimating the variance in halo abundance when the sample volume is small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Neutrino signatures of 100 2D Axisymmetric Core-Collapse Supernova Simulations

ABSTRACT We present in this paper a public data release of an unprecedentedly large set of core-collapse supernova (CCSN) neutrino emission models, comprising 100 detailed 2D axisymmetric radiation-hydrodynamic simulations evolved out to as late as ∼5 s post-bounce and spanning an extensive range of massive-star progenitors. The motivation for this paper is to provide a physically and numerically uniform benchmark data set to the broader neutrino detection community to help it characterize and optimize subsurface facilities for what is likely to be a once-in-a-lifetime galactic supernova burst event. With this release, we hope to (1) help the international experiment and modelling communities more efficiently optimize the retrieval of physical information about the next galactic CCSN, (2) facilitate the better understanding of core-collapse theory and modelling among interested experimentalists, and (3) help further integrate the broader supernova neutrino community.

Astronomy & Astrophysics↗
DOE OSTI2020

From hydrogen to helium: the spectral evolution of white dwarfs as evidence for convective mixing

ABSTRACT We present a study of the hypothesis that white dwarfs undergo a spectral change from hydrogen- to helium-dominated atmospheres using a volume-limited photometric sample drawn from the Gaia-DR2 catalogue, the Sloan Digital Sky Survey (SDSS), and the Galaxy Evolution Explorer (GALEX). We exploit the strength of the Balmer jump in hydrogen-atmosphere DA white dwarfs to separate them from helium-dominated objects in SDSS colour space. Across the effective temperature range from 20 000 to 9000 K, we find that 22 per cent of white dwarfs will undergo a spectral change, with no spectral evolution being ruled out at 5σ. The most likely explanation is that the increase in He-rich objects is caused by the convective mixing of DA stars with thin hydrogen layers, in which helium is dredged up from deeper layers by a surface hydrogen convection zone. The rate of change in the fraction of He-rich objects as a function of temperature, coupled with a recent grid of 3D radiation-hydrodynamic simulations of convective DA white dwarfs – which include the full overshoot region – lead to a discussion on the distribution of total hydrogen mass in white dwarfs. We find that 60 per cent of white dwarfs must have a hydrogen mass larger than MH/MWD = 10−10, another 25 per cent have masses in the range MH/MWD = 10−14–10−10, and 15 per cent have less hydrogen than MH/MWD = 10−14. These results have implications for white dwarf asteroseismology, stellar evolution through the asymptotic giant branch and accretion of planetesimals on to white dwarfs.

Cunningham, Tim↗
DOE OSTI2024

Systematic Local Simulations of Fast Neutrino Flavor Conversions with Scattering Effects

We investigate the dynamics of fast neutrino flavor conversions (FFCs) in the one-dimensional (1D) and zero-dimensional (0D) models, in which spatial advection is considered and ignored, respectively. In this study, we employ snapshots obtained by our self-consistent, realistic Boltzmann-neutrino-radiation-hydrodynamics simulations. We show that the FFC growth rate is considerably larger in the 1D model than in the 0D model, as expected from the previous linear analysis results. We find that the momentum space dimension does not significantly influence the neutrino transition probability in 1D models. On the other hand, in the 0D model without collisions, the FFC depends on the momentum space, and the azimuthal angle dependence breaks the periodicity of the FFC. Our study demonstrates that collisional instability can lead to further flavor conversions on a long timescale in 1D models after the asymptotic state of FFC has been reached. Such an effect should be taken into consideration when the fast and collisional flavor instabilities coexist.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Direct drive with the argon fluoride laser as a path to high fusion gain with sub-megajoule laser energy

Argon fluoride (ArF) is currently the shortest wavelength laser that can credibly scale to the energy and power required for high gain inertial fusion. ArF's deep ultraviolet light and capability to provide much wider bandwidth than other contemporary inertial confinement fusion (ICF) laser drivers would drastically improve the laser target coupling efficiency and enable substantially higher pressures to drive an implosion. Our radiation hydrodynamics simulations indicate gains greater than 100 are feasible with a sub-megajoule ArF driver. Our laser kinetics simulations indicate that the electron beam-pumped ArF laser can have intrinsic efficiencies of more than 16%, versus about 12% for the next most efficient krypton fluoride excimer laser. We expect at least 10% ‘wall plug' efficiency for delivering ArF light to target should be achievable using solid-state pulsed power and efficient electron beam transport to the laser gas that was demonstrated with the U.S. Naval Research Laboratory's Electra facility. These advantages could enable the development of modest size and lower cost fusion power plant modules. This would drastically change the present view on inertial fusion energy as being too expensive and the power plant size too large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Gamma-ray signatures of 𝑟-process radioactivity from the collapse of magnetized white dwarfs

We predict the gamma-ray line emission from 𝑟-process nuclei synthesized in the ejecta of the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf. Using ejecta from a two-dimensional general-relativistic neutrino-magnetohydrodynamic simulation, further evolved with a radiation-hydrodynamics code coupled to an in situ nuclear reaction network, we construct angle-dependent gamma-ray spectra in the 0.01–10 MeV band via composition-dependent ray tracing through the ejecta. The emission between ∼1 and 10 d is dominated by 132 I (𝑡 1/2 = 2.3 h), continuously replenished by the decay of its parent 132 Te (𝑡 1/2 = 3.2 d), with additional contributions from 131 I, 133 Xe, and 132 Te . At 𝑡 ≳ 20 d, 56 Co (from 56 Ni decay) becomes the primary emitter. The simultaneous presence of 𝑟 process and iron-peak gamma-ray lines is distinctive of AIC ejecta and absent in binary neutron star mergers, where iron-peak nuclei are generally not synthesized. Comparing with the 3⁢𝜎 continuum sensitivities of planned MeV gamma-ray telescopes (COSI, AMEGO-X, e-ASTROGAM, GRAMS, GammaTPC), we find the brightest 𝑟-process lines detectable to ∼10 Mpc by GammaTPC and GRAMS, with the signal approaching their sensitivity threshold at 30 Mpc. As a result, the 𝑟-process spectral features survive time integration over ∼30 d exposures, demonstrating robustness against the long observation times required by gamma-ray detectors.

Nuclear reactions↗
DOE OSTI2021

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernovae (CCSN), since they substantially reduce the computational burden inherent in the multidimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zeroth and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. In this work, we find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\textit{κ}$). This issue originates from the insensitiveness of incoming neutrinos to $\textit{κ}$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of electron- neutrinos lepton number (ELN)-crossing searches with the two-moment neutrino transport. This analysis is complementary to our [L. Johns and H. Nagakura arXiv:2104.04106] that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\textit{κ}$ ≳ 0.5. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗