Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Direct-drive double-shell implosion: A platform for burning-plasma physics studies

Double-shell ignition designs have been investigated with the indirect-drive inertial confinement fusion (ICF) scheme in both simulations and experiments in which the inner-shell kinetic energy was limited to ~10 to 15 kJ, even driven by megajoule-class lasers such as the National Ignition Facility. Since direct-drive ICF can couple more energy to the imploding shells, we have performed a detailed study on direct-drive double-shell (D 3 S) implosions with state-of-the-art physics models implemented in radiation-hydrodynamic codes (LILAC and DRACO), including nonlocal thermal transport, cross-beam energy transfer (CBET), and first-principles-based material properties. To mitigate classical unstable interfaces, we have proposed the use of a tungsten/beryllium–mixed inner shell with gradient-density layers that can be made by magnetron sputtering. In our D 3 S designs, a 70-μm-thick beryllium outer shell is driven symmetrically by a high-adiabat (α ≥ 10), 1.9-MJ laser pulse to a peak velocity of ~240 km/s. Upon spherical impact, the outer shell transfers ~30 to 40 kJ of kinetic energy to the inner shell filled with deuterium–tritium gas or liquid, giving neutron-yield energies of ~6 MJ in 1-D simulations. Two-dimensional high-mode DRACO simulations indicated that such high-adiabat D 3 S implosions are not susceptible to laser imprint, but the long-wavelength perturbations from the laser port configuration along with CBET can be detrimental to the target performance. Yet, neutron yields of ~0.3- to 1.0-MJ energies can still be obtained from our high-mode DRACO simulations. The robust α-particle bootstrap is readily reached, which could provide a viable platform for burning-plasma physics studies. Once CBET mitigation and/or more laser energy becomes available, we believe that breakeven or moderate energy gain might be feasible with the proposed D 3 S scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Diagnostic signature of the compressibility of the inertial-confinement-fusion pusher

In this study, carbon shell areal density measurements from many types of inertial confinement fusion implosions at the National Ignition Facility (NIF) demonstrate that the final state of the outside portion of the shell is set primarily by capsule coast time, the coasting period between main laser shut off and peak fusion output. However, the fuel areal density does not correlate with the increasing carbon compression. While two-dimensional (2D) radiation-hydrodynamic simulations successfully capture the carbon compression, energy must be added to the simulated fuel-ice layer to reproduce fuel areal density measurements. The data presented demonstrates that the degradation mechanisms that reduce the compressibility of the fuel do not reduce the compressibility of the ablator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Observations of anomalous x-ray emission at early stages of hot-spot formation in deuterium-tritium cryogenic implosions

In DT cryogenic implosions, hot-spot x-ray self-emission is observed to begin at a larger shell radius than is predicted by a 1-D radiation-hydrodynamic implosion model. Laser-imprint is shown to explain the observation for a low-adiabat implosion. Furthermore, for more-stable implosions the data are not described by the imprint model and suggest there are additional sources of decompression of the dense fuel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/ cm 3 and 0 to 10 9 K) equation-of-state (EOS) table for a CH 1.72 O 0.37 N 0.086 quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. Here, the simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental evidence of early-time saturation of the ion-Weibel instability in counterstreaming plasmas of CH, Al, and Cu

The collisionless ion-Weibel instability is a leading candidate mechanism for the formation of collisionless shocks in many astrophysical systems, where the typical distance between particle collisions is much larger than the system size. Multiple laboratory experiments aimed at studying this process utilize laser-driven (I≳10 15 W/cm 2 ), counterstreaming plasma flows (V≲ 2000 km/s) to create conditions unstable to Weibel-filamentation and growth. This technique intrinsically produces temporally varying plasma conditions at the midplane of the interaction where Weibel-driven B fields are generated and studied. Experiments discussed herein demonstrate robust formation of Weibel-driven B fields under multiple plasma conditions using CH, Al, and Cu plasmas. Linear theory based on benchmarked radiation-hydrodynamic FLASH calculations is compared with Fourier analyses of proton images taken ~5–6 linear growth times into the evolution. Our new analyses presented here indicate that the low-density, high-velocity plasma-conditions present during the first linear-growth time (~300–500 ps) sets the spectral characteristics of Weibel filaments during the entire evolution. It is shown that the dominant wavelength (~300 μm) at saturation persists well into the nonlinear phase, consistent with theory under these experimental conditions. However, estimates of B-field strength, while difficult to determine accurately due to the path-integrated nature of proton imaging, are shown to be in the ~10–30 T range, an order of magnitude above the expected saturation limit in homogenous plamas but consistent with enhanced B fields in the midplane due to temporally varying plasma conditions in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Structures of strong shocks in low-density helium and neon gases

Strong shocks are essential components in many high-energy-density environments such as inertial confinement fusion implosions. However, the experimental measurements of the spatial structures of such shocks are sparse. In this paper, the soft x-ray emission of a shock front in a helium gas mixture (90% helium, 10% neon) and a pure neon gas was spatially resolved using an imaging spectrometer. We observe that the shock width in the helium mixture gas is about twice as large as in the pure neon gas. Moreover, they exhibit different precursor layers, where electron temperature greatly exceeds ion temperature, extending for more than ~350 µm with the helium gas mixture but less than 30 µm in the pure neon. At the shock front, calculations show that the electrons are strongly collisional with mean-free path two orders of magnitude shorter than the characteristic length of the shock. However, the helium ions can reach a kinetic regime as a consequence of their mean-free path being comparable to the scale of the shock. A radiation-hydrodynamic simulation demonstrates the impact of thermal conduction on the formation of the precursors with charge state, Z, playing a major role in heat flow and the precursor formation in both the helium mixture and the pure neon gases. Particle-in-cell simulations are also performed to study the ion kinetic effects on the formation of the observed precursors. A group of fast-streaming ions is observed leading the shock only in the helium gas mixture. Additionally, both effects explain the longer precursor layer in the helium shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Thermal conductivity of a laser plasma

Here, we present a model of the electron thermal conductivity of a laser-produced plasma. The model, supported by Vlasov–Fokker–Planck simulations, predicts that laser absorption reduces conductivity by forcing electrons out of a Maxwell–Boltzmann equilibrium, which results in depletion of both low-velocity bulk electrons and high-velocity tail electrons. We show that both the bulk and tail electrons approximately follow super-Gaussian distributions, but with distinct exponents that each depend on the laser intensity and wavelength through the parameter a = Zv$^{2}_{E}$/v$^{2}_{T}$. For a value of a = 0.5, tail depletion reduces the thermal conductivity to half its zero-intensity value. We present our results as simple analytic fits that can be readily implemented in any radiation-hydrodynamics code or used to correct the local limit of non-local conduction models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inverse bremsstrahlung absorption rate for super-Gaussian electron distribution functions including plasma screening

Here we provide analytic expressions for the effective Coulomb logarithm for inverse bremsstrahlung absorption which predict significant corrections to the Langdon effect and overall absorption rate compared to previous estimates. The calculation of the collisional absorption rate of laser energy in a plasma by the inverse bremsstrahlung mechanism usually makes the approximation of a constant Coulomb logarithm. We dispense with this approximation and instead take into account the velocity dependence of the Coulomb logarithm, leading to a more accurate expression for the absorption rate valid in both classical and quantum conditions. In contrast to previous work, the laser intensity enters into the Coulomb logarithm. In most laser-plasma interactions the electron distribution function is super-Gaussian [Langdon, Phys. Rev. Lett. 44, 575 (1980)], and we find the absorption rate under these conditions is increased by as much as ≈ 30% compared to previous estimates at low density. In many cases of interest the correction to Langdon's predicted reduction in absorption is large; for example at Z = 6 and Te = 400 eV the Langdon prediction for the absorption is in error by a factor of ≈ 2. However, we also account for the additional effect of plasma screening, which predicts a reduction in absorption by a similar amount (up to ≈ 30%). These two effects compete to determine the overall absorption, which may be increased or decreased, depending on the conditions. The corrections can be incorporated into radiation-hydrodynamics simulation codes by replacing the familiar Coulomb logarithm with an analytic expression which depends on the super-Gaussian order “M” and the screening length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivity of the dynamic-shell target to laser drive nonuniformities

The dynamic-shell concept for inertial confinement fusion (ICF) uses an initially homogeneous target and a carefully shaped laser pulse to form a shell and implode it. The laser pulse consists of a series of pickets that drive shocks into the target. The first few shocks converge inwards and rebound from the center of the target, creating an expanding, low-density plasma. Subsequent shocks are launched into the expanding plasma and eventually coalesce to form a shell, which is then imploded with a traditional ICF laser pulse. This study describes radiation-hydrodynamic simulations that investigate the sensitivity of dynamic-shell targets to imperfections in the laser drive. A one-dimensional (1D) study looks at mistiming and power variations in the pickets and a two-dimensional (2D) study examines irradiation perturbations imposed by the laser-beam geometry. Simulations show that less than ~2% power imbalance or 200 ps timing variation in the pickets is sufficient to keep the yield above 90% of the maximum. Additionally, the 2D simulations show that 72 or more beams are required to keep irradiation nonuniformities low enough to obtain fusion yields close to that of 1D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Direct-drive measurements of laser-imprint-induced shock velocity nonuniformities [First direct-drive measurements of laser-imprint-induced shock velocity nonuniformities]

Perturbations in the velocity profile of a laser-ablation-driven shock wave seeded by speckle in the spatial beam intensity (i.e., laser imprint) have been measured. Direct measurements of these velocity perturbations were recorded using a two-dimensional high-resolution velocimeter probing plastic material shocked by a 100-ps picket laser pulse from the OMEGA laser system. The measured results for experiments with one, two, and five overlapping beams incident on the target clearly demonstrate a reduction in long-wavelength (>25– μ m) perturbations with an increasing number of overlapping laser beams, consistent with theoretical expectations. Furthermore, these experimental measurements are crucial to validate radiation-hydrodynamics simulations of laser imprint for laser direct drive inertial confinement fusion research since they highlight the significant (factor of 3) underestimation of the level of seeded perturbation when the microphysics processes for initial plasma formation, such as multiphoton ionization are neglected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Rarefaction Flows and Mitigation of Imprint in Direct-Drive Implosions

Using highly resolved 3-D radiation-hydrodynamic simulations, we identify a novel mechanism by which the deleterious impact of laser imprinting is mitigated in direct-drive inertial confinement fusion. Unsupported shocks and associated rarefaction flows, commonly produced with short laser bursts, are found to reduce imprint modulations prior to target acceleration. Optimization through the choice of laser pulse with picket(s) and target dimensions may improve the stability of lower adiabat designs, thus providing the necessary margin for ignition-relevant implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Anomalous Absorption by the Two-Plasmon Decay Instability

Radiation-hydrodynamic simulations of directly driven fusion experiments at the Omega Laser Facility predict absorption accurately when targets are driven at low overlapped laser intensity. Discrepancies appear at increased intensity, however, with higher-than-expected laser absorption on target. Strong correlations with signatures of the two-plasmon decay (TPD) instability—including half-harmonic and hard-x-ray emission—indicate that TPD is responsible for this anomalous absorption. Scattered light data suggest that up to ≈ 30 % of the laser power reaching quarter-critical density can be absorbed locally when the TPD threshold is exceeded. A scaling of absorption versus TPD threshold parameter was empirically determined and validated using the laser–plasma simulation environment code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Shock Ignition Laser-Plasma Interactions in Ignition-Scale Plasmas

Here, we use a subignition scale laser, the 30 kJ Omega, and a novel shallow-cone target to study laser-plasma interactions at the ablation-plasma density scale lengths and laser intensities anticipated for direct drive shock-ignition implosions at National Ignition Facility scale. Our results show that, under these conditions, the dominant instability is convective stimulated Raman scatter with experimental evidence of two plasmon decay (TPD) only when the density scale length is reduced. Particle-in-cell simulations indicate this is due to TPD being shifted to lower densities, removing the experimental back-scatter signature and reducing the hot-electron temperature. The experimental laser energy-coupling to hot electrons was found to be 1%–2.5%, with electron temperatures between 35 and 45 keV. Radiation-hydrodynamics simulations employing these hot-electron characteristics indicate that they should not preheat the fuel in MJ-scale shock ignition experiments.

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↗