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 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↗
DOE OSTI2026

Utilizing the deuterium-tritium fusion resonance to diagnose thermal runaway in igniting plasmas

For high-efficiency inertial confinement fusion implosions, it is predicted that a burning hot spot will successfully encompass all surrounding fuel and then transition into a thermal runaway where the internal energy increase from fusion occurs on a timescale faster than the expansion of the fuel is able to quench the fusion chain reaction after ignition occurs. Observation of this dynamic phase transition would indicate distinct burn properties and indicate an implosion's robustness. A technique for diagnosing the presence of thermal runaway from measurements of nuclear reaction history is presented. The technique is based on taking the logarithmic derivative of the nuclear reaction history, called the 𝛼 curve, and allowing a mathematical decoupling of the mass, volume, and thermal reactivity in the fusion reaction rate equation. During thermal runaway, where the thermal temperature dominates the burn dynamics, a maximum in the 𝛼 curve is found where there is a maximum in the first derivative of the thermal fusion reactivity, an effect to the deuterium-tritium (DT) fusion cross-section resonance. This provides a distinct signature related to the fundamental nature of the DT fusion nuclear resonance and signifies the transition into the fusion thermal instability. Impacts of charged particle transport on the effect are also assessed and the analytical formulas are compared and found to be in agreement with radiation hydrodynamic codes.

high-energy-density plasmas↗
DOE OSTI2021

Revisiting Theoretical Limits for One Degree-of-Freedom Wave Energy Converters

This work revisits the theoretical limits of one degree-of-freedom wave energy converters (WECs). This work considers the floating sphere used in the OES Task 10 WEC modeling and verification effort for analysis. Analytical equations are derived to determine bounds on displacement amplitude, time-averaged power (TAP), and power-take-off (PTO) force. A unique result found shows that the TAP absorbed by a WEC can be defined solely by the inertial properties and radiation hydrodynamic coefficients. In addition, a unique expression for the PTO force was derived that provides lower and upper bounds when resistive control is used to maximize power generation. For complex conjugate control, this same expression only provides a lower bound, as there is theoretically no upper bound. These bounds assist in comparing the performance of the floating sphere if it were to extract energy using surge or heave motion. The analysis shows because of differences in hydrodynamic coefficients for each oscillating mode, there are different frequency ranges that provide better power capture efficiency. The influence of a motion constraint on TAP while utilizing a nonideal power take-off is examined and found to reduce the losses associated with bidirectional energy flow. The expression to calculate TAP with a nonideal PTO is modified by the electrical conversion efficiency and the ratio of the PTO spring and damping coefficients. Furthermore, the PTO spring and damping coefficients were separated in the expression, allowing for limits to be set on the PTO coefficients to ensure net power generation.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Revisiting Theoretical Limits for One-Degree-of-Freedom Wave Energy Converters

This work revisits the theoretical limits of one-degree-of-freedom wave energy converters. This paper considers the floating sphere used in the Ocean Energy Systems Task 10 Wave Energy Converter modeling and verification effort for analysis. Analytical equations are derived to determine bounds on the motion amplitude, time-averaged power, and power-take-off (PTO) force. A unique result was found that shows the time-averaged power absorbed by a wave energy converter can be defined solely by the inertial properties and radiation hydrodynamic coefficients. In addition, a unique expression for the PTO force amplitude was derived that has provided upper and lower bounds when resistive control is used to maximize power generation. For complex conjugate control, this same expression can only provide a lower bound, as there is theoretically no upper bound. These bounds are used to compare the performance of a floating sphere if it were to extract energy using surge or heave motion. The analysis shows that because of the differences in hydrodynamic coefficients of each oscillating mode, there will be different frequency ranges that provide better power capture efficiency. The influence of a motion constraint on power absorption while also utilizing a nonideal power take-off is examined and found to reduce the losses associated with bidirectional energy flow. The expression to calculate the time-averaged power with a nonideal PTO is modified by the mechanical-to-electrical efficiency and the ratio of the PTO spring and damping coefficients. The PTO spring and damping coefficients were separated in the expression, which allows for limits to be set on the possible values of PTO coefficients to ensure a net flow of power to the grid.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2025

Predicting fusion ignition at the National Ignition Facility with physics-informed deep learning

Here, an inertial confinement fusion experiment, carried out at the National Ignition Facility, has achieved ignition by generating fusion energy exceeding the laser energy that drove the experiment. Prior to the experiment, a generative machine learning model that combines radiation hydrodynamics simulations, deep learning, experimental data, and Bayesian statistics was used to predict, with a probability greater than 70%, that ignition was the most likely outcome for this shot.

Spears, Brian K. [Lawrence Livermore National Labo↗
DOE OSTI2024

artemis

Astrophysical multifluid radiation hydrodynamics code

Ryan, Ben↗
DOE OSTI2022

Multi-frame, ultrafast, x-ray microscope for imaging shockwave dynamics

Inertial confinement fusion (ICF) holds increasing promise as a potential source of abundant, clean energy, but has been impeded by defects such as micro-voids in the ablator layer of the fuel capsules. It is critical to understand how these micro-voids interact with the laser-driven shock waves that compress the fuel pellet. At the Matter in Extreme Conditions (MEC) instrument at the Linac Coherent Light Source (LCLS), we utilized an x-ray pulse train with ns separation, an x-ray microscope, and an ultrafast x-ray imaging (UXI) detector to image shock wave interactions with micro-voids. To minimize the high- and low-frequency variations of the captured images, we incorporated principal component analysis (PCA) and image alignment for flat-field correction. After applying these techniques we generated phase and attenuation maps from a 2D hydrodynamic radiation code (xRAGE), which were used to simulate XPCI images that we qualitatively compare with experimental images, providing a one-to-one comparison for benchmarking material performance. Moreover, we implement a transport-of-intensity (TIE) based method to obtain the average projected mass density (areal density) of our experimental images, yielding insight into how defect-bearing ablator materials alter microstructural feature evolution, material compression, and shock wave propagation on ICF-relevant time scales.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Understanding Transients: Needs from the Laboratory Astrophysics Community

A growing number of astrophysical transients are pushing astronomers to develop increasingly complex computational tools to model both radiation hydrodynamics and electron transport. Here we review the physics needs and simulation uncertainties associated with modeling these transients. Although this review will focus on theory and simulation aspects of this problem, we also review some experiments designed to study this physics.

79 ASTRONOMY AND ASTROPHYSICS↗