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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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 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 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↗
DOE OSTI2022

3D Simulations Capture the Persistent Low-Mode Asymmetries Evident in Laser-Direct-Drive Implosions on OMEGA

Spherical implosions in Inertial Confinement Fusion (ICF) are inherently sensitive to perturbations that may arise from experimental constraints and errors. Control and mitigation of low-mode (long wavelengths) perturbations is a key milestone to improving implosion performances. Here, we present the first 3-D radiation-hydrodynamic simulations of directly driven ICF implosions with an inline package for polarized Crossed-Beam Energy Transfer (CBET). Simulations match bang times, yields (separately accounting for laser-induced high modes and fuel age), hot spot flow velocities and direction, for which polarized CBET contributes to the systematic flow orientation evident in the OMEGA implosion database. Current levels of beam mispointing, imbalance, target offset and asymmetry from polarized CBET degrade yields by more than 40%. The effectiveness of two mitigation strategies for low-modes is explored.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Inverse Bremsstrahlung Absorption

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Reduced Model of Ionization Lag in Intense Laser-Produced Plasmas

A physics-based empirical formula is derived to predict the ionization lag in underdense plasmas generated by an intense laser. Time-dependent nonlocal thermodynamic equilibrium calculations demonstrate significantly delayed ionization, due to rapid changes in plasma conditions, which critically impacts plasma properties such as opacity, emissivity, and heat transport. The reduced model, based on these calculations, enables the estimation of ionization lag without requiring in-depth knowledge of nonlocal thermodynamic equilibrium modeling. Furthermore, modeling reveals that the two-step ionization process—collisional excitation followed by photoionization—plays a crucial role in this ionization delay, with collisional excitation setting the timescale for ionization. Simulations across a range of elements, from beryllium to germanium, demonstrate that ionization lag is a widespread phenomenon, underscoring the importance of incorporating such processes into ionization modeling in radiation hydrodynamic simulations for various laser-plasma experiments.

direct drive↗
DOE OSTI2024

Imaging the dynamics of initial laser-driven shocks and blowoff plasmas in polystyrene under laser-direct-drive fusion conditions

The dynamics of laser-driven shock propagation in solid ablator material, along with the concomitant processes of shinethrough and blowoff plasma plume expansion, at conditions relevant to laser-direct-drive fusion are measured using a tabletop system capable of supporting rapid dataset development with significantly higher resolution compared to existing methods. The experimental results are directly compared with current radiation-hydrodynamic simulations using the two-dimensional code . Discrepancies between experiments and simulations are evidenced and attributed to limitation of simulation to account for important processes involved in the solid-to-plasma transition. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

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

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

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

Exploring composition mixing in kilonova ejecta with ray-by-ray simulations

Binary neutron star merger (BNSM) ejecta are considered a primary repository of r-process nucleosynthesis and a source of the observed heavy-element abundances. We implement composition mixing into ray-by-ray radiation-hydrodynamic simulations of BNSM ejecta, coupled with an online nuclear network (NN). We model mixing via a gradient-based mixing approximation that evolves simultaneously with the hydrodynamics. Here, we find that mixing occurs in regions where the electron fraction changes rapidly. While mixing smooths composition gradients in transition regions, it has a negligible impact on the heavy-element yields. This is because the primary r-process site (the equatorial ejecta) is initially homogeneous in free neutrons, leaving no strong gradients for mixing to act upon. In each angular ray, the abundances of the most produced elements are robust under mixing, while the less abundant ones are more affected. The total global abundances change only slightly from mixing, since each angular ray contributes its most abundant elements. Furthermore, the predicted kilonova light curves show only minor reddening, with differences below the detectability of state-of-the-art telescopes. In general, we do not observe significant effects from mixing in the time span of the r-process. Consequently, mixing only leads to minor variations in abundances and light curves in ray-by-ray simulations.

Explosive burning↗