UC Berkeley

Plasmas are made up of charged particles whose short-range and long-range interactions give rise to complex behavior that can be difficult to fully characterize experimentally. One of the most complete theoretical descriptions of a plasma is that of kinetic theory, which treats each particle species as a probability distribution function in a six-dimensional position-velocity phase space. Drawing on statistical mechanics, these distribution functions mathematically represent a system of interacting particles without tracking individual ions and electrons. The evolution of the distribution function(s) is governed by the Boltzmann equation coupled to Maxwell's equations, which together describe the dynamics of the plasma and the associated electromagnetic fields. When collisions can be neglected, the Boltzmann equation is reduced to the Vlasov equation. High-fidelity simulation of the rich physics in even a subset of the full six-dimensional phase space calls for low-noise high-accuracy numerical methods. To that end, this dissertation investigates a fourth-order finite-volume discretization of the Vlasov-Maxwell equation system, and addresses some of the fundamental challenges associated with applying these types of computationally intensive enhanced-accuracy numerical methods to phase space simulations.

The governing equations of kinetic theory are described in detail, and their conservation-law weak form is derived for Cartesian and cylindrical phase space coordinates. This formulation is well known when it comes to Cartesian geometries, as it is used in finite-volume and finite-element discretizations to guarantee local conservation for numerical solutions. By contrast, the conservation-law weak form of the Vlasov equation in cylindrical phase space coordinates is largely unexplored, and to the author's knowledge has never previously been solved numerically. Thereby the methods described in this dissertation for simulating plasmas in cylindrical phase space coordinates present a new development in the field of computational plasma physics.

A fourth-order finite-volume method for solving the Vlasov-Maxwell equation system is presented first for Cartesian and then for cylindrical phase space coordinates. Special attention is given to the treatment of the discrete primary variables and to the quadrature rule for evaluating the surface and line integrals that appear in the governing equations. The finite-volume treatment of conducting wall and axis boundaries is particularly nuanced when it comes to phase space coordinates, and is described in detail. In addition to the mechanics of each part of the finite-volume discretization in the two different coordinate systems, the complete algorithm is also presented.

The Cartesian coordinate discretization is applied to several well-known test problems. Since even linear analysis of kinetic theory governing equations is complicated on account of velocity being an independent coordinate, few analytic or semi-analytic predictions exist. Benchmarks are particularly scarce for configurations that have magnetic fields and involve more than two phase space dimensions. Ensuring that simulations are true to the physics thus presents a difficulty in the development of robust numerical methods. The research described in this dissertation addresses this challenge through the development of more complete physics-based benchmarks based on the Dory-Guest-Harris instability. The instability is a special case of perpendicularly-propagating kinetic electrostatic waves in a warm uniformly magnetized plasma. A complete derivation of the closed-form linear theory dispersion relation for the instability is presented. The electric field growth rates and oscillation frequencies specified by the dispersion relation provide concrete measures against which simulation results can be quantitatively compared. Furthermore, a specialized form of perturbation is shown to strongly excite the fastest growing mode. The fourth-order finite-volume algorithm is benchmarked against the instability, and is demonstrated to have good convergence properties and close agreement with theoretical growth rate and oscillation frequency predictions.

The Dory-Guest-Harris instability benchmark extends the scope of standard test problems by providing a substantive means of validating continuum kinetic simulations of warm magnetized plasmas in higher-dimensional 3D (x,v_x,v_y) phase space. The linear theory analysis, initial conditions, algorithm description, and comparisons between theoretical predictions and simulation results are presented.

The cylindrical coordinate finite-volume discretization is applied to model axisymmetric systems. Since mitigating the prohibitive computational cost of simulating six dimensions is another challenge in phase space simulations, the development of a robust means of exploiting symmetry is a major advance when it comes to numerically solving the Vlasov-Maxwell equation system. The discretization is applied to a uniform distribution function to assess the nature of the singularity at the axis, and is demonstrated to converge at fourth-order accuracy. The numerical method is then applied to simulate electrostatic ion confinement in an axisymmetric Z-pinch configuration. To the author's knowledge this presents the first instance of a conservative finite-volume discretization of the cylindrical coordinate Vlasov equation. The computational framework for the Vlasov-Maxwell solver is described, and an outlook for future research is presented.