Brian M. Haines

First Liquid Layer Inertial Confinement Fusion Implosions at the National Ignition Facility.

The first cryogenic deuterium and deuterium-tritium liquid layer implosions at the National Ignition Facility (NIF) demonstrate D_{2} and DT layer inertial confinement fusion (ICF) implosions that can access a low-to-moderate hot-spot convergence ratio (1230) DT ice layer implosions. Although high CR is desirable in an idealized 1D sense, it amplifies the deleterious effects of asymmetries. To date, these asymmetries prevented the achievement of ignition at the NIF and are the major cause of simulation-experiment disagreement. In the initial liquid layer experiments, high neutron yields were achieved with CRs of 12-17, and the hot-spot formation is well understood, demonstrated by a good agreement between the experimental data and the radiation hydrodynamic simulations. These initial experiments open a new NIF experimental capability that provides an opportunity to explore the relationship between hot-spot convergence ratio and the robustness of hot-spot formation during ICF implosions.

The effects of plasma diffusion and viscosity on turbulent instability growth

We perform two-dimensional simulations of strongly–driven compressible Rayleigh–Taylor and Kelvin–Helmholtz instabilities with and without plasma transport phenomena, modeling plasma species diffusion, and plasma viscosity in order to determine their effects on the growth of the hydrodynamic instabilities. Simulations are performed in hydrodynamically similar boxes of varying sizes, ranging from 1 μm to 1 cm in order to determine the scale at which plasma effects become important. Our results suggest that these plasma effects become noticeable when the box size is approximately 100 μm, they become significant in the 10 μm box, and dominate when the box size is 1 μm. Results suggest that plasma transport may be important at scales and conditions relevant to inertial confinement fusion, and that a plasma fluid model is capable of representing some of the kinetic transport effects.

Reynolds-averaged Navier–Stokes initialization and benchmarking in shock-driven turbulent mixing

We investigate a strategy for benchmarking Reynolds-averaged Navier–Stokes (RANS) models by comparing moments extracted from averaged large eddy simulation (LES) data and those predicted directly by RANS. We consider the Besnard–Harlow–Rauenzahn (BHR) RANS approach designed for variable-density compressible flows, which has been previously applied to a wide variety of turbulence problems of interest. We focus on the models ability to predict moments relevant to shock-driven material mixing. A prototypical inverse chevron shock tube configuration is considered, for which laboratory and previous LES studies are available for comparison and validation. We show that when appropriately initialized, BHR is capable of accurately capturing various characteristic integral measures of the flow; strategies for initialization are demonstrated while addressing sensitivity of BHR predictions to closure and initialization specifics, initial material interface conditions, and grid resolution. The reference simulations are performed using implicit LES based on the Los Alamos National Laboratory RAGE hydrodynamics code.

Detailed high-resolution three-dimensional simulations of OMEGA separated reactants inertial confinement fusion experiments

We present results from the comparison of high-resolution three-dimensional (3D) simulations with data from the implosions of inertial confinement fusion capsules with separated reactants performed on the OMEGA laser facility. Each capsule, referred to as a “CD Mixcap,” is filled with tritium and has a polystyrene (CH) shell with a deuterated polystyrene (CD) layer whose burial depth is varied. In these implosions, fusion reactions between deuterium and tritium ions can occur only in the presence of atomic mix between the gas fill and shell material. The simulations feature accurate models for all known experimental asymmetries and do not employ any adjustable parameters to improve agreement with experimental data. Simulations are performed with the RAGE radiation-hydrodynamics code using an Implicit Large Eddy Simulation (ILES) strategy for the hydrodynamics. We obtain good agreement with the experimental data, including the DT/TT neutron yield ratios used to diagnose mix, for all burial depths of the deuterated shell layer. Additionally, simulations demonstrate good agreement with converged simulations employing explicit models for plasma diffusion and viscosity, suggesting that the implicit sub-grid model used in ILES is sufficient to model these processes in these experiments. In our simulations, mixing is driven by short-wavelength asymmetries and longer-wavelength features are responsible for developing flows that transport mixed material towards the center of the hot spot. Mix material transported by this process is responsible for most of the mix (DT) yield even for the capsule with a CD layer adjacent to the tritium fuel. Consistent with our previous results, mix does not play a significant role in TT neutron yield degradation; instead, this is dominated by the displacement of fuel from the center of the implosion due to the development of turbulent instabilities seeded by long-wavelength asymmetries. Through these processes, the long-wavelength asymmetries degrade TT yield more than the DT yield and thus bring DT/TT neutron yield ratios into agreement with experiment. Finally, we present a detailed comparison of the flows in 2D and 3D simulations.

Simulations of material mixing in laser-driven reshock experiments

We perform simulations of a laser-driven reshock experiment [Welser-Sherrill et al., High Energy Density Phys. (unpublished)] in the strong-shock high energy-density regime to better understand material mixing driven by the Richtmyer–Meshkov instability. Validation of the simulations is based on direct comparison of simulation and radiographic data. Simulations are also compared with published direct numerical simulation and the theory of homogeneous isotropic turbulence. Despite the fact that the flow is neither homogeneous, isotropic nor fully turbulent, there are local regions in which the flow demonstrates characteristics of homogeneous isotropic turbulence. We identify and isolate these regions by the presence of high levels of turbulent kinetic energy (TKE) and vorticity. After reshock, our analysis shows characteristics consistent with those of incompressible isotropic turbulence. Self-similarity and effective Reynolds number assessments suggest that the results are reasonably converged at the fine...

Simulation ensemble for a laser–driven shear experiment

We perform an ensemble of simulations of a laser-driven shear experiment [L. Welser-Sherrill et al., “Two laser-driven mix experiments to study reshock and shear,” High Energy Density Phys. J. 9(3), 496–499 (2013)] in the strong-shock high energy-density regime to better understand material mixing driven by the Kelvin–Helmholtz instability. Each simulation uses a different realization of random initial interface perturbations based on data from targets used in experiments. Validation of the simulations is based on direct comparison of simulation and radiographic data. Simulations are also compared with published direct numerical simulation and the theory of homogeneous isotropic turbulence. Despite the fact that the flow is neither homogeneous, isotropic, nor fully turbulent, there are local regions in which the flow demonstrates characteristics of homogeneous isotropic turbulence. Our analysis shows characteristics consistent with those of incompressible isotropic turbulence. Our results show that turbul...

The effects of convergence ratio on the implosion behavior of DT layered inertial confinement fusion capsules

The wetted foam capsule design for inertial confinement fusion capsules, which includes a foam layer wetted with deuterium-tritium liquid, enables layered capsule implosions with a wide range of hot-spot convergence ratios (CR) on the National Ignition Facility. We present a full-scale wetted foam capsule design that demonstrates high gain in one-dimensional simulations. In these simulations, increasing the convergence ratio leads to an improved capsule yield due to higher hot-spot temperatures and increased fuel areal density. High-resolution two-dimensional simulations of this design are presented with detailed and well resolved models for the capsule fill tube, support tent, surface roughness, and predicted asymmetries in the x-ray drive. Our modeling of these asymmetries is validated by comparisons with available experimental data. In 2D simulations of the full-scale wetted foam capsule design, jetting caused by the fill tube is prevented by the expansion of the tungsten-doped shell layer due to prehe...

High-resolution modeling of indirectly driven high-convergence layered inertial confinement fusion capsule implosions

In this paper, we present the results of high-resolution simulations of the implosion of high-convergence layered indirect-drive inertial confinement fusion capsules of the type fielded on the National Ignition Facility using the xRAGE radiation-hydrodynamics code. In order to evaluate the suitability of xRAGE to model such experiments, we benchmark simulation results against available experimental data, including shock-timing, shock-velocity, and shell trajectory data, as well as hydrodynamic instability growth rates. We discuss the code improvements that were necessary in order to achieve favorable comparisons with these data. Due to its use of adaptive mesh refinement and Eulerian hydrodynamics, xRAGE is particularly well suited for high-resolution study of multi-scale engineering features such as the capsule support tent and fill tube, which are known to impact the performance of high-convergence capsule implosions. High-resolution two-dimensional (2D) simulations including accurate and well-resolved ...

Exponential yield sensitivity to long-wavelength asymmetries in three-dimensional simulations of inertial confinement fusion capsule implosions

In this paper, we perform a series of high-resolution 3D simulations of an OMEGA-type inertial confinement fusion (ICF) capsule implosion with varying levels of initial long-wavelength asymmetries in order to establish the physical energy loss mechanism for observed yield degradation due to long-wavelength asymmetries in symcap (gas-filled capsule) implosions. These simulations demonstrate that, as the magnitude of the initial asymmetries is increased, shell kinetic energy is increasingly retained in the shell instead of being converted to fuel internal energy. This is caused by the displacement of fuel mass away from and shell material into the center of the implosion due to complex vortical flows seeded by the long-wavelength asymmetries. These flows are not fully turbulent, but demonstrate mode coupling through non-linear instability development during shell stagnation and late-time shock interactions with the shell interface. We quantify this effect by defining a separation lengthscale between the fue...

Analysis of the effects of energy deposition on shock-driven turbulent mixing

We perform simulations of laser-driven turbulence experiments with energy deposition, modeling situations where energy is deposited in a mixing layer before or after it is reshocked. Such situations are experienced in, e.g., inertial confinement fusion capsules. We show that the timing of the energy deposition has a significant impact on the development of turbulent flow and corresponding turbulent material mixing. In particular, if the energy is deposited before the shock wave begins interacting with the mixing layer, the development of turbulence and turbulent mixing are inhibited. If, however, the energy is deposited after the shock wave has interacted with a portion of the mixing layer, turbulence generation and turbulent mixing are enhanced.