Inertial Confinement Fusion Implosions Research Articles

Investigation of effects of initial interface conditions on the two-dimensional single-mode compressible Rayleigh–Taylor instability: Based on the discrete Boltzmann method

The Rayleigh–Taylor (RT) instability in inertial confinement fusion implosions evolves at the unstable interface of two fluids when the light fluid is pushing the heavy one. The effects of the initial amplitude and transition layer on the compressible RT instability are investigated numerically by using the discrete Boltzmann method. On the one hand, during the RT evolution, higher initial amplitudes initially increase the global density gradient and non-equilibrium area, with a subsequent reversal. The increasing initial amplitude leads to an initial rise followed by a decline in the system’s maximum Mach number. On the other hand, the impact of the transition layer is generally opposite to the one of initial amplitude in the RT process. These findings offer significant insights into controlling and understanding RT instability in fusion implosion scenarios, emphasizing novel aspects relative to existing literature.

Effects of drive pulse shape on graded metal pushered single shell capsule implosions on the National Ignition Facility

Graded metal pushered single shells (PSS) are a viable alternative to low-Z capsules (Z is the atomic number) for indirect drive inertial confinement fusion implosions due to enhanced core tamping and radiation trapping, but they can be compromised by the pusher mixing with the fuel. We compare 2-shock and 3-shock laser pulses for Be/Cr PSS capsules filled with deuterium–tritium gas fuel at 6 mg/cc density. 1D radiation-hydrodynamic simulations predict higher core compression and, hence, ∼2× higher fusion yield for the 3-shock drive than for 2-shock. Nevertheless, we observe similar core ion temperatures and fusion yields for both drives. The implosion burn duration is 25% shorter and the core volume is ∼2.5× smaller for the 3-shock drive than for 2-shock, consistent with a higher compression. 1D LASNEX mix simulations using a buoyancy-drag model matching the measured yields also agree with the observed core sizes and burn durations and suggest ∼40% and ∼70% yield degradations for 2-shock and 3-shock drives due to hydrodynamic instabilities and atomic mix at the pusher–fuel interface. At the same time, 2D HYDRA simulations show that mid-mode (2–250) instability degradations are negligible for the 2-shock implosion (9%) and significant (45%) for 3-shock. Subtracting these from the 1D mix simulations, we infer similar degradations from high-mode instabilities and atomic mix for both drives. Due to its robustness to mid-mode instabilities, future pusher–gas mix studies will use the 2-shock drive.

Validation of implosion modeling through direct-drive shock timing experiments at the National Ignition Facility.

Precise modeling of shocks in inertial confinement fusion implosions is critical for obtaining the desired compression in experiments. Shock velocities and postshock conditions are determined by laser-energy deposition, heat conduction, and equationsof state. This paper describes experiments at the National Ignition Facility (NIF) [E. M. Campbell and W. J. Hogan, Plasma Phys. Control. Fusion 41, B39 (1999)10.1088/0741-3335/41/12B/303] where multiple shocks are launched into a cone-in-shell target made of polystyrene, using laser-pulse shapes with two or three pickets and varying on-target intensities. Shocks are diagnosed using the velocity interferometric system for any reflector (VISAR) diagnostic [P.M. Celliers et al., Rev. Sci. Instrum. 75, 4916 (2004)0034-674810.1063/1.1807008]. Simulated and inferred shock velocities agree well for the range of intensities studied in this work. These directly-driven shock-timing experiments on the NIF provide a good measure of early-time laser-energy coupling. The validated models add to the credibility of direct-drive-ignition designs at the megajoule scale.

Modeling stimulated Brillouin backscatter from outer-cone quads across multiple inertial confinement fusion hohlraum designs

Stimulated Brillouin scattering (SBS) is a potential risk for laser damage in the experiments carried out at the National Ignition Facility (NIF), and by altering the energy deposition pattern in hohlraums, it affects the symmetry of indirect-drive inertial confinement fusion implosions. We have surveyed backward SBS on outer-cone quads across NIF integrated hohlraums of various platforms numerically, using three-dimensional (3D) simulations with the backscatter code pF3D [Berger et al. Phys. Plasmas 5, 4337 (1998)] and ray-based gain calculations. Simulated reflected powers and energies, as well as the spectrum of reflected light all compare favorably with measurements. Ray-based calculations of exponential SBS amplification (“gain”), which assume a strongly damped plasma wave and steady-state response, are performed using a novel method that includes the 3D speckled field of the laser that drives SBS. This approach is useful for understanding qualitative differences between hohlraum designs and identifying regions susceptible to SBS within hohlraums. Quantitatively, gains are not found to correlate with SBS reflectivities in 3D, necessitating fully wave-based calculations that naturally include diffraction and various temporal dependencies.

Enhanced sensitivity to target offset when using cross-beam energy transfer mitigation techniques in direct-drive inertial confinement fusion implosions

In direct-drive inertial confinement fusion, target offset from the target chamber center (or center of beam convergence) may lead to significant implosion asymmetry and fusion yield degradation. In addition, cross-beam energy transfer (CBET) has been shown to be a significant source of laser energy scattering and leads to a reduction in implosion velocity and yield. To improve energy coupling and implosion performance, several techniques for CBET mitigation have been proposed. Recent simulations, however, have shown that CBET also substantially mitigates the effect of target offset on implosion asymmetry and yield [Anderson et al., Phys. Plasmas 27, 112713 (2020)]. Furthermore, the inclusion of CBET models in radiation-hydrodynamics codes was shown to greatly improve agreement between simulations and experiments involving substantial target offset distances. This paper explores the intensity dependence of this CBET–offset effect. In addition, it is shown that enhanced sensitivity to target offset can be expected when CBET-mitigation techniques are used in direct-drive implosions. This is shown through simulations of two such CBET-mitigation techniques on the OMEGA laser: (1) decreased beam-to-target radius, and (2) beam-to-beam frequency detuning. For the typical target offset distances (<15 μm) observed in experiments on OMEGA, however, overall yield is still anticipated to be substantially higher when CBET-mitigation techniques are employed.

Coupling dynamics of capsule interior defects and its impact on hydrodynamic instabilities at ablation fronts for inertial confinement fusion implosions

The micrometer-scale internal defect in the capsule is one of the most important factors that limit implosion performance in inertial confinement fusion (ICF) experiments, which creates instability seeds as shocks propagate through the capsule shell. Here, we report the generation mechanism of vortex pairs resulting from the interaction of shock waves with multiple bubbles, as well as the origin of more intricate perturbation waves than those observed in the case of single defects. Based on the subsequent evolution of hydrodynamic instability, it is evident that the vortex pairs induce the emergence of low-density (light-bubble case) or high-density (referred to as heavy-bubble case) jets on the ablative front. The presence of multiple side-by-side defects can rapidly amplify the dimensions of the jet. These jets could be responsible for the “meteor shower” observed in implosion experiments. Converging disturbed waves between vertically aligned defects lead to a more complex nonlinear flow field evolution compared to the scenario with a single defect. A systematic study of localized perturbation growth as a function of defect placement is presented. We investigate the dependence of circulation in the flow field on the locations of the defects. The scanning results of defect scenes with different sizes revealed the reason why the depth of fluid penetration is affected by the position and size, and found that the effects of the position and size on the perturbation expansion width can be equivalent to a certain extent. The extension of the perturbation width when the defect is off-axis limits the degree of penetration of the perturbation depth. The results contribute to a more comprehensive understanding of physical processes, such as the seeding mechanism, shell integrity, and mass injection into the central region, which may be applied to inform the development of more effective strategies to mitigate implosion degradation in ICF implosion experiments.

Demonstration of hot-spot fuel gain exceeding unity in direct-drive inertial confinement fusion implosions

Demonstration of hot-spot fuel gain exceeding unity in direct-drive inertial confinement fusion implosions

Krypton He[formula omitted] line spectrum including n = 2, 3 Li-like satellites with detailed Stark broadened line shapes

We discuss the krypton Heβ line spectrum including Li-like satellites with a spectator electron in n = 2 and n = 3 and detailed line shapes computed using standard Stark broadening theory for hot dense plasma conditions relevant to X-ray tracer spectroscopy of inertial confinement fusion implosion cores. The results show that the interference term in the electron broadening does not produce a significant effect for these satellite transitions. However, the effect of the electric field mixing of the energy levels driven by the ion’s microfield distribution does produce a significant change in the line shape. Level populations calculated with a collisional radiative atomic kinetics model were employed to obtain the photon energy resolved emissivity and opacity using the Stark line shapes, and the emergent intensity distribution was calculated by integrating the radiation transport equation along chords assuming a uniform spherical plasma source. The line spectrum has electron temperature (Te) and density (ne) sensitivity due to the temperature and density dependence of level populations and the density dependence of the Stark line shapes. Hence, this spectrum is suitable for a simultaneous temperature and density plasma diagnostic of implosion cores.

Effect of hot-electron preheating on the multimode bubble-front growth of the ablative Rayleigh–Taylor instability

Hot electrons (HEs) generated via parametric instabilities at high laser intensities are a critical concern of laser-driven inertial confinement fusion (ICF), which can significantly impact the ICF performance by preheating the target. In this paper, the effects of HE preheating with moderate HE energy on the evolution of two-dimensional multimode ablative Rayleigh–Taylor instability (ARTI) up to the self-similar growth stage are studied through numerical simulations with a multigroup diffusion model. It is found that HE preheating stabilizes the linear growth of multimode ARTI and delays the onset of the self-similar growth regime. This time delay is more significant for the short-wavelength mode ARTI and higher energy HE cases. It is also shown that the variation of self-similar growth coefficients under HE preheating is not very significant. The delay to the onset of the nonlinear stage of multimode ARTI by HE preheating with moderate energy may be beneficial to ICF implosions.

Ultrahigh-energy density state in nanowire arrays irradiated with picosecond kilojoule-class ultra-intense laser

Ultrahigh-energy density (UHED) states greater than 1 Gbar pressure are typically observed under extreme conditions, such as in the core of an inertial confinement fusion implosion. A novel alternative approach for generating volumetric UHED states is to use nanowire arrays irradiated with a femtosecond ultrahigh-intensity laser. In this paper, we present an experimental investigation on laser absorption and energy transport in nanowire arrays irradiated with a picosecond kilojoule petawatt laser. The laser–target interactions were studied by measuring the x-ray emission and escaping hot electrons from a bare Cu foil and a foil with a nanowire array grown on its surface. The measured Cu-Kα and He-α emissions from the nanowire array were higher than those from the flat foil. In addition, hot electrons observed from the front surface of the nanowire array were enhanced. On the other hand, despite the stronger Kα emission from the nanowire array and the enhancement of hot electrons escaping from the front surface of the nanowire array, the number of hot electrons observed from the rear side of the flat foil target was slightly lower than that of the flat foil. A comparison of the experimental results with the results of a two-dimensional particle-in-cell simulation code suggested that the magnetic fields generated around the periodic nanowire array trap hot electrons, improving the electron-to-target energy coupling and efficiently producing UHED states.

X-ray continuum spectroscopy of inertial confinement fusion implosions at the National Ignition Facility

A methodology for measuring x-ray continuum spectra of inertial confinement fusion (ICF) implosions is described. The method relies on the use of ConSpec, a high-throughput spectrometer using a highly annealed pyrolytic graphite crystal [MacDonald et al., J. Instrum. 14, P12009 (2019)], which measures the spectra in the ≃20–30 keV range. Due to its conical shape, the crystal is sagittally focusing a Bragg-reflected x-ray spectrum into a line, which enhances the recorded x-ray emission signal above the high neutron-induced background accompanying ICF implosions at the National Ignition Facility. To improve the overall measurement accuracy, the sensitivity of the spectrometer measured in an off-line x-ray laboratory setting was revised. The error analysis was expanded to include the accuracy of the off-line measurements, the effect of the neutron-induced background, as well as the influence of possible errors in alignment of the instrument to the ICF target. We demonstrate how the improved methodology is applied in the analysis of ConSpec data with examples of a relatively low-neutron-yield implosion using a tritium–hydrogen–deuterium mix as a fuel and a high-yield deuterium–tritium (DT) implosion producing high level of the background. In both cases, the shape of the measured spectrum agrees with the exponentially decaying spectral shape of bremsstrahlung emission to within ±10%. In the case of the high-yield DT experiment, non-monotonic deviations slightly exceeding the measurement uncertainties are observed and discussed.

Measuring higher-order moments of neutron-time-of-flight data for cryogenic inertial confinement fusion implosions on OMEGA

Ion temperatures serve as an important diagnostic for inertial confinement fusion (ICF) implosions. In direct-drive ICF experiments on OMEGA, neutron-time-of-flight (nTOF) data are used to infer the ion temperature of the fusing plasma produced in the implosion experiment. The analysis of the nTOF data requires an assumption about the shape of the underlying source signal. Since the source nTOF signal is a near-replica of the neutron energy spectrum, an ideal Gaussian shape, corresponding to the neutron energy spectrum of a uniform temperature plasma, is routinely employed. However, spatial and temporal variations of the ion temperature in the plasma give rise to higher-order moments, which were first described by Munro [Nucl. Fusion 56, 036001 (2016)]. In this work, we show a simpler alternative analysis to derive moments of the neutron energy spectrum for a plasma with variations in ion temperature. We also present a revised analysis of measured nTOF signals that uses a model with an additional degree of freedom to take into account the effect of ion temperature variations on the shape of the spectrum. Compared to presently used nTOF analysis, the revised analysis yields on average ≈2× more accurate fits to the data and up to 15% higher ion temperatures for cryogenic experiments. Furthermore, we quantify the ion temperature inflation caused by radially symmetric fluid flows, which are present even in a symmetric implosion, and which serve as a lower bound on the ion temperature inflation in real implosions.

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications.

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. We utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code lilac to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. Results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

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. Both effects explain the longer precursor layer in the helium shock.

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

Correlations between asymmetric compression, burn amplification, and hot-spot velocities in inertial confinement fusion implosions

This manuscript examines the correlations between the hot-spot velocity (an observable signature of residual kinetic energy), low-mode implosion asymmetries, and burn amplification in inertial confinement fusion implosions on the National Ignition Facility (NIF). Using a combination of two-dimensional axis-symmetric and three-dimensional radiation-hydrodynamic simulations coupled to neutronics, we find that for typical NIF implosions, the stagnation asymmetry multiplies the observed hot-spot velocity anywhere from 80% to 120%, while burn amplification always increases it. Additionally, we find stagnation asymmetry typically deflects the observed hot-spot flow. The two mechanisms (low-mode implosion asymmetries and burn amplification) can be decoupled, and application of a simple model to a database of cryogenic implosions on the NIF infers the total hot-spot velocity amplification. This finding modifies the interpretation of data collected from inertial confinement fusion experiments and impacts the magnitude and origin of low-mode asymmetries.

Large-incidence-angle multiple-beam two-plasmon decay instability in inertial confinement fusion

A multi-dimensional code FLAME-MD solving fluid-like laser-plasma instabilities' (LPIs) equations has been developed and is used to study multiple-beam two-plasmon decay (TPD) instability relevant to a laser-entrance-hole window burn-off scenario in inertial confinement fusion (ICF) experiments. It is found that TPD can be collectively driven by multiple beams incident at large incidence angles with respect to the electron density gradient at a very low threshold. The polarization configuration of the laser beams is a key factor determining the way of sharing common daughter electron plasma waves (EPWs). The p-polarized beams arranged on the same incidence cone can collectively drive common EPWs along the cone axis. The common-wave sharing mechanisms among the s-polarized beams are largely dependent on the geometry of the beams and are less robust. The simulation results also show that the p-polarized beams are dominating the multiple-beam TPD processes. The common EPWs along the cone axis can accelerate energetic electrons toward the capsule inside the hohlraum and, therefore, pose a fuel-preheat risk to ICF implosions.

Temperature relaxation and generalized Coulomb logarithm in two-temperature dense plasmas relevant to inertial confinement fusion implosions

Detailed knowledge of energy exchange between electrons and ions is of fundamental importance for the description of temperature relaxation and also other nonequilibrium physics in Inertial Confinement Fusion (ICF). We present a theoretical model for the temperature relaxation rate and the related generalized Coulomb logarithm based on the Quantum Lenard–Balescu (QLB) kinetic equation, where no special cutoffs are needed to be introduced. To describe the collective modes characterizing the ionic acoustic waves, a single-pole approximation is introduced for the ionic dielectric response. The proposed model for the generalized Coulomb logarithm is examined over a wide range of plasma conditions for electron temperatures between 102 and 105 eV and electron densities between 1022 and . The values of the generalized Coulomb logarithm are demonstrated to be in excellent agreement with the ones evaluated using the original QLB kinetic approach but with much less computational cost. Comparisons with molecular dynamics simulations and other theoretical approaches are presented. For further applications of our model, we present results for the recently measured experimental Coulomb logarithm. Compared to other widely applied models such as the Landau–Spitzer Coulomb logarithm, our model provides more consistent descriptions for the results of molecular dynamics simulations and also for the experimental outcomes. Our model for the generalized Coulomb logarithm is easy to calculate and can benefit efficient and reliable simulations for the ICF implosions.

Reduced mixing in inertial confinement fusion with early-time interface acceleration.

In inertial confinement fusion (ICF) implosions, the interface between the cryogenic DT fuel and the ablator is unstable to shock acceleration (the Richtmyer-Meshkov instability, RM) and constant acceleration (Rayleigh-Taylor instability, RT). Instability growth at this interface can reduce the final compression, limiting fusion burnup. If the constant acceleration is in the direction of the lighter material (negative Atwood number), the RT instability produces oscillatory motion that can stabilize against RM growth. Theory and simulations suggest this scenario occurred at early times in some ICF experiments on the National Ignition Facility, possibly explaining their favorable performance compared to one-dimensional simulations. This characteristic is being included in newer, lower adiabat designs, seeking to improve compression while minimizing ablator mixing into the fuel.

Dynamics and Power Balance of Near Unity Target Gain Inertial Confinement Fusion Implosions.

The change in the power balance, temporal dynamics, emission weighted size, temperature, mass, and areal density of inertially confined fusion plasmas have been quantified for experiments that reach target gains up to 0.72. It is observed that as the target gain rises, increased rates of self-heating initially overcome expansion power losses. This leads to reacting plasmas that reach peak fusion production at later times with increased size, temperature, mass and with lower emission weighted areal densities. Analytic models are consistent with the observations and inferences for how these quantities evolve as the rate of fusion self-heating, fusion yield, and target gain increase. At peak fusion production, it is found that as temperatures and target gains rise, the expansion power loss increases to a near constant ratio of the fusion self-heating power. This is consistent with models that indicate that the expansion losses dominate the dynamics in this regime.