Accurate treatment of energetic fusion byproducts in laboratory plasmas often requires a kinetic description, owing to their large birth kinetic energy and long mean-free-paths compared with the characteristic system scale lengths. For example, alpha particles produced by deuterium–tritium fusion reactions are born at high energies (3.5MeV) and predominantly slow down through interactions with electrons traveling at comparable speeds. As an alpha particle slows, its distribution collapses near the background ion-thermal speed, forming a sharp structure in velocity space. Such sharp features pose numerical challenges in grid-based Eulerian methods: capturing the full alpha-particle energies demands a large velocity domain, while resolving the near-thermal region requires a sufficiently fine mesh. Inspired by the work of Peigney et al.[J. Comput. Phys. 278 (2014)], we present a two-grid approach that splits the alpha-particle distribution into energetic (suprathermal) and ash (thermal) components. A Gaussian-based sink term transfers particles from the energetic population to the ash population as they slow to the thermal regime, and a conservative projection scheme ensures that mass, momentum, and energy of the alpha and ash interactions are preserved. Unlike the formulation of Peigney, our method does not require a strict asymptotic separation of velocity scales, which can, in principle, be arbitrary. We demonstrate the robustness of this approach on challenging multiscale problems, including a surrogate for an igniting inertial confinement fusion capsule.

On December 5, 2022, an indirect drive fusion implosion on the National Ignition Facility (NIF) achieved a target gain G_{target} of 1.5. This is the first laboratory demonstration of exceeding "scientific breakeven" (or G_{target}>1) where 2.05 MJ of 351 nm laser light produced 3.1 MJ of total fusion yield, a result which significantly exceeds the Lawson criterion for fusion ignition as reported in a previous NIF implosion [H. Abu-Shawareb et al. (Indirect Drive ICF Collaboration), Phys. Rev. Lett. 129, 075001 (2022)PRLTAO0031-900710.1103/PhysRevLett.129.075001]. This achievement is the culmination of more than five decades of research and gives proof that laboratory fusion, based on fundamental physics principles, is possible. This Letter reports on the target, laser, design, and experimental advancements that led to this result.

We report on simulations of counter-propagating laser-produced plasmas in an inertial confinement fusion (ICF) hohlraum surrogate, aiming to replicate observations reported by Le Pape et. al in recent work. The conditions of the colliding plasmas are relevant to ICF hohlraums used for indirect-drive ignition, and are obtained both with and without low-density He-gas fill. We compare experimental diagnostics to outputs from simulations using the 1D-2V Vlasov-Fokker-Planck kinetic code iFP and the xRAGE radiation-hydrodynamics code. These include the inferred radial lineouts of inferred ion number fraction and ion and electron temperatures, as well as the reported experimental Thomson-scattering (TS) spectra (compared via synthetic TS diagnostics). We observe that 1D kinetic simulations capture the plasma states reported in the experimental diagnostics quite well. Counter-intuitively, the kinetic simulations capture the gas-fill experiment (expected to be more `hydro-like') better than the vacuum experiment, while the reverse is observed for hydrodynamic simulations. This is attributed to the presence of non-trivial multi-dimensional hydrodynamic effects which are more dominant in the vacuum experiment. These effects are somewhat inhibited in the gas-fill experiment, permitting quasi-1D kinetic plasma transport to play more of a role in producing plasma interpenetration. Differences between the effects of Maxwellian vs. non-Maxwellian (`full f') synthetic TS diagnostics are investigated for the kinetic simulations. We find non-Maxwellian TS spectra differ non-trivially from Maxwellian spectra, which suggests caution may be warranted when applying Maxwellian TS models to infer plasma conditions via backward modeling when kinetic effects may be present.

The classical work of Braginskii [Zh. Eksp. Teor. Fiz. 33, 459 (1957)] published almost 65 years ago was the first to provide a complete, closed fluid description of a weakly coupled, fully ionized, collisional plasma immersed in a magnetic field. While this fact is not widely known or appreciated, the Braginskii expressions for the electron transport coefficients can under- or overestimate the said coefficients by up to a factor of two for the electron Hall parameter of order unity (with the Hall parameter being proportional to the electron gyro-frequency over the electron–ion collision frequency) and can provide incorrect Hall-parameter scalings for its large values. Starting with the work of Epperlein and Haines [Phys. Fluids 29, 1029 (1986)], several papers attempted to correct the Braginskii electron results with varying degrees of success. Herein, we present our own effort with a hope to finally put this problem to rest.

We develop a conservative phase-space grid-adaptivity strategy for the Vlasov-Fokker-Planck equation in a planar geometry. The velocity-space grid is normalized to the thermal speed and shifted by the bulk-fluid velocity. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The grid adaptation ensures discrete conservation of the collisional invariants (mass, momentum, and energy). The conservative grid-adaptivity strategy provides an efficient scheme which resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity.

excited. Therefore, simulating these interactions requires tools that include wave particle interactions and that include wave nonlinearities. One methodology for studying such interactions is particle-in-cell (PIC) simulations. While PIC codes include most of the relevant physics they are also the most computer intensive. However, with the development of sophisticated software and the use of massively parallel computers, PIC codes can now be used to accurately study a wide range of problems in HEDS. The research in this project involved building, maintaining, and using the UCLA parallel computing infrastructure. This infrastructure includes the codes OSIRIS and UPIC which have been improved or developed during this grant period. Specifically, we used this PIC infrastructure to study laser-plasma interactions relevant to future NIF experiments and high-intensity laser and beam plasma interactions relevant to fast ignition fusion. The research has led to fundamental knowledge in how to write parallel PIC codes and use parallel PIC simulations, as well as increased the fundamental knowledge of HEDS. This fundamental knowledge will not only impact Inertial Confinement Fusion but other fields such as plasma-based acceleration and astrophysics.

We report on the first steady-state simulations of strong plasma shocks with fully kinetic ions and electrons, independently confirmed by two fully kinetic codes (an Eulerian continuum and a Lagrangian particle-in-cell). While kinetic electrons do not fundamentally change the shock structure as compared with fluid electrons, we find an appreciable rearrangement of the preheat layer, associated with nonlocal electron heat transport effects. The electron heat flux profile qualitatively agrees between kinetic and fluid electron models, suggesting a certain level of "stiffness", though substantial nonlocality is observed in the kinetic heat flux. We also find good agreement with nonlocal electron heat-flux closures proposed in the literature. Finally, in contrast to the classical hydrodynamic picture, we find a significant collapse in the "precursor" electric-field shock at the preheat layer edge, which correlates with the electron-temperature gradient relaxation.

hydrodynamic ion model, designed for nanosecond time-scale laser-plasma interactions. Heat-flux effects in Ohm’s law under non-local conditions was investigated; physics that is not well captured by standard numerical models but is nevertheless important in fusion-related scenarios. Under such conditions there are numerous interesting physical effects, such as collisional magnetic instabilities, amplification of magnetic fields, re-emergence of non-locality through magnetic convection, and reconnection of magnetic field lines and redistribution of thermal energy. In this project highlights included the first full-scale kinetic simulations of a magnetized hohlraum and the discovery of a new magnetic reconnection mechanism, as well as a completed PhD thesis and the production of a new code for Inertial Fusion research.

We present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

We consider the solution of the fully kinetic (including electrons) Vlasov-Ampère system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian algorithm. The problem of velocity-space meshing for disparate thermal and bulk velocities is dealt with by an adaptive coordinate transformation of the Vlasov equation for each species, which is then discretized, including the resulting inertial terms. Mass, momentum, and energy are conserved, and Gauss's law is enforced to within the nonlinear convergence tolerance of the iterative solver through a set of nonlinear constraint functions while permitting significant flexibility in choosing discretizations in time, configuration, and velocity space. We mitigate the temporal stiffness introduced by, e.g., the plasma frequency through the use of high-order/low-order (HOLO) acceleration of the iterative implicit solver. We present several numerical results for canonical problems of varying degrees of complexity, including the multiscale ion-acoustic shock wave problem, which demonstrate the efficacy, accuracy, and efficiency of the scheme.

Revolver and Double Shell Inertial Confinement Fusion capsule designs hope to achieve a robust volumetric thermonuclear burn via the use of a high-Z pusher shell filled with a cryogenic D–T fuel. Unfortunately, mix of the pusher material into the fuel (gas) may adversely impact the burn performance. Hydrodynamic instability of the metal/gas interface as the mix source is an obvious concern, but 1D effects may also be detrimental. Such effects include plasma diffusion at material interfaces, which has been the subject of numerous theoretical, computational, and experimental investigations. However, other 1D mix mechanisms may exist, which have yet to be thoroughly explored. In particular, plasma kinetic effects may drive the mix when a shock breaks out of the metal/gas interface. Using the state-of-the-art, hybrid (kinetic-ion/fluid electron), multi-ion Vlasov–Fokker–Planck code, iFP, we show that shock-driven kinetic effects can reconfigure the interface and the interfacial width subsequently grows diffusively. Finally, we consider any implications for high-Z pusher designs.

A multi-species plasma ion transport model has been added to the adaptive mesh refinement radiation hydrodynamics code, xRage, to include kinetic transport effects when the particle distributions are near Maxwellian, with deviations proportional to a Knudsen number smaller than one. The model is first verified against self-similar solutions reported previously for the pressure equilibrium case, and next shown to be relatively insensitive to the choice of equation of state for the ions. Simulations are then used to examine Inertial Confinement Fusion dynamics in a 1D spherical geometry characteristic of an Omega implosion with a plastic (CH) shell containing a deuterium-tritium (DT) fuel, and with an added heavy ion impurity, argon. Even in this simplified 1D geometry, several interesting results are apparent. Ion stratification occurs similarly to that reported previously in purely kinetic simulations. The hydrogen in the plastic shell is transported radially inward, carried with the main drive shock, and thus migrates away from the C ions. The fuel D and T ions show the expected stratification with an increase in the lighter species concentration during the shock implosion and a reversal, with heavier species concentrations enhanced after shock expansion from the center. This stratification during burn yields different burn weighted ion temperatures, Ti, for the reactions, Ti[DD] < Ti[DT] < Ti[TT], consistent in their ordering with experiments. The mix widths per ion, measured where concentrations fall to 10% of their interfacial value, are evaluated as a function of time, and these are seen to be significant (of order 10 μm) even at early times, well before the main shock converges and before the shell deceleration. The 1D geometry may be a reasonable approximation for this early time mix and implies that this transport may play a role in reducing or modifying the instabilities driven by initial perturbations, ablation, and Rayleigh-Taylor unstable deceleration. An apparent depletion of the heavier ions seen at the incoming ion shock front warrants further investigation.

Fuel-ion species dynamics in hydrodynamiclike shock-driven DT^{3}He-filled inertial confinement fusion implosion is quantitatively assessed for the first time using simultaneously measured D^{3}He and DT reaction histories. These reaction histories are measured with the particle x-ray temporal diagnostic, which captures the relative timing between different nuclear burns with unprecedented precision (∼10  ps). The observed 50±10  ps earlier D^{3}He reaction history timing (relative to DT) cannot be explained by average-ion hydrodynamic simulations and is attributed to fuel-ion species separation between the D, T, and ^{3}He ions during shock convergence and rebound. At the onset of the shock burn, inferred ^{3}He/T fuel ratio in the burn region using the measured reaction histories is much higher as compared to the initial gas-filled ratio. As T and ^{3}He have the same mass but different charge, these results indicate that the charge-to-mass ratio plays an important role in driving fuel-ion species separation during strong shock propagation even for these hydrodynamiclike plasmas.

Submitted for the DPP19 Meeting of The American Physical Society The Effects of Kinetically Enhanced Interfacial Mix on Reactivity Variation in Omega Capsules1 WILLIAM TAITANO, ANDREI SIMAKOV, LUIS CHACON, BRETT KEENAN, STEVEN ANDERSON, Los Alamos National Laboratory — Analyzing the effects of the hydrodynamic mix has been an essential part of designing ICF targets. In contrast, atomic mixing processes due to kinetic effects have largely been overlooked until recent years. Various studies have been performed both theoretical and computational to assess the role of these effects on the yield. A particular experiment (amongst many) which has challenged our understanding of both kinetic and hydrodynamic theory is the so-called Rygg experiment [1]. The targets were comprised of a direct drive CH capsule with a D-3He fill that was varied in concentration while ensuring hydro-equivalence (the initial total mass and pressure are kept fixed). In the experiments, anomalous (non-hydrodynamic) yield variations were observed that had eluded explanations thus far. A recent fuel (D-3He) only Vlasov-Fokker-Planck (VFP) simulation [2] of a similar experiment has shown how kinetically enhanced fuel-species separation can partially predict the observational yield trend. In this work, we present our results on similar simulations with the CH pusher partially modeled to assess the interplay of interface dynamics on the yield. [1] J. R. Rygg et al., Phys. Plasmas, 13, 052702 (2006) [2] W. T. Taitano et al., Phys. Plasmas, 25, 056310 (2018) 1This work was sponsored by the Thermonuclear Burn Initiative of ASC, and the LANL Institutional Computing. This work was performed under the NNSA of the USDOE at LANL under contract 89233218CNA000001. William Taitano Los Alamos National Laboratory Date submitted: 02 Jul 2019 Electronic form version 1.4

Presentation Conclusions: We have derived manifold-preserving algorithms for kinetic plasma simulation; Collisionless (PIC): we have solved the 40-year-old algorithmic challenge of developing accurate implicit PIC algorithms; Collisional (VFP): we have demonstrated a truly multiscale algorithm that has enabled routine simulation of ICF spherical capsule implosions with a few hundred cores for a couple of days.; In both cases: Strict conservation properties have been shown to be critical for long-term accuracy.; Significant algorithmic acceleration has been achieved by using nested asymptotic models.; We have seen similar benefits in other applications