An Eulerian Vlasov-Fokker-Planck Algorithm for Spherical Implosion Simulations of Inertial Confinement Fusion Capsules

Abstract

We present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov-Fokker-Planck (VFP) simulation of an inertial confinement fusion (ICF) capsule implosion. 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.

62 pages, 34 figures

References in corpus (14)

• Thermo-diffusion in inertially confined plasmas
• Energy-conserving time propagation for a geometric particle-in-cell Vlasov--Maxwell solver
• Species separation and modification of neutron diagnostics in inertial-confinement fusion
• Deciphering the Kinetic Structure of Multi-Ion Plasma Shocks
• Ion Species Stratification Within Strong Shocks in Two-Ion Plasmas
• A Rescaling Velocity Method for Dissipative Kinetic Equations - Applications to Granular Media
• Simulation and assessment of ion kinetic effects in a direct-drive capsule implosion experiment
• An Adaptive, Implicit, Conservative 1D-2V Multi-Species Vlasov-Fokker-Planck Multiscale Solve in Planar Geometry
• Fokker Planck kinetic modeling of suprathermal alpha-particles in a fusion plasma
• Observation and modeling of interspecies ion separation in inertial confinement fusion implosions via imaging x-ray spectroscopy
• Kinetic simulations of fusion ignition with hot-spot ablator mix
• An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system
• Positivity-preserving scheme for two-dimensional advection-diffusion equations including mixed derivatives
• ADI type preconditioners for the steady state inhomogeneous Vlasov equation

Cited by in corpus (6)

• First principles simulations of dense hydrogen
• Fokker-planck Multi-species Equations In The Adiabatic Asymptotics
• Experimental Measurements of Ion Diffusion Coefficients and Heating in a Multi-Ion-Species Plasma Shock
• Optimal Transport for Parameter Identification of Chaotic Dynamics via Invariant Measures
• Hyperbolic Machine Learning Moment Closures for the BGK Equations
• Fully kinetic simulations of strong steady-state collisional planar plasma shocks