A Discontinuous Galerkin Method for Ideal Two-Fluid Plasma Equations

A discontinuous Galerkin method for the ideal 5 moment two-fluid plasma system is presented. The method uses a second or third order discontinuous Galerkin spatial discretization and a third order TVD Runge-Kutta time stepping scheme. The method is benchmarked against an analytic solution of a dispersive electron acoustic square pulse as well as the two-fluid electromagnetic shock and existing numerical solutions to the GEM challenge magnetic reconnection problem. The algorithm can be generalized to arbitrary geometries and three dimensions. An approach to maintaining small gauge errors based on error propagation is suggested.Comment: 40 pages, 18 figures

Phase space eigenfunctions with applications to continuum kinetic simulations

Continuum kinetic simulations are increasingly capable of resolving high-dimensional phase space with advances in computing. These capabilities can be more fully explored by using linear kinetic theory to initialize the self-consistent field and phase space perturbations of kinetic instabilities. The phase space perturbation of a kinetic eigenfunction in unmagnetized plasma has a simple analytic form, and in magnetized plasma may be well approximated by truncation of a cyclotron-harmonic expansion. We catalogue the most common use cases with a historical discussion of kinetic eigenfunctions and by conducting nonlinear Vlasov-Poisson and Vlasov-Maxwell simulations of single- and multi-mode two-stream, loss-cone, and Weibel instabilities in unmagnetized and magnetized plasmas with one- and two-dimensional geometries. Applications to quasilinear kinetic theory are discussed and applied to the bump-on-tail instability. In order to compute eigenvalues we present novel representations of the dielectric function for ring distributions in magnetized plasmas with power series, hypergeometric, and trigonometric integral forms. Eigenfunction phase space fluctuations are visualized for prototypical cases such as the Bernstein modes to build intuition. In addition, phase portraits are presented for the magnetic well associated with nonlinear saturation of the Weibel instability, distinguishing current-density-generating trapping structures from charge-density-generating ones.Comment: 51 pages, 26 figures, 4 appendice

Hybridizable discontinuous Galerkin methods for solving the two-fluid plasma model

The two-fluid plasma model has a wide range of timescales which must all be numerically resolved regardless of the timescale on which plasma dynamics occurs. The answer to solving numerically stiff systems is generally to utilize unconditionally stable implicit time advance methods. Hybridizable discontinuous Galerkin (HDG) methods have emerged as a powerful tool for solving stiff partial differential equations. The HDG framework combines the advantages of the discontinuous Galerkin (DG) method, such as high-order accuracy and flexibility in handling mixed hyperbolic/parabolic PDEs with the advantage of classical continuous finite element methods for constructing small numerically stable global systems which can be solved implicitly. In this research we quantify the numerical stability conditions for the two-fluid equations and demonstrate how HDG can be used to avoid the strict stability requirements while maintaining high order accurate results

Numerical approximation of the Euler-Poisson-Boltzmann model in the quasineutral limit

This paper analyzes various schemes for the Euler-Poisson-Boltzmann (EPB) model of plasma physics. This model consists of the pressureless gas dynamics equations coupled with the Poisson equation and where the Boltzmann relation relates the potential to the electron density. If the quasi-neutral assumption is made, the Poisson equation is replaced by the constraint of zero local charge and the model reduces to the Isothermal Compressible Euler (ICE) model. We compare a numerical strategy based on the EPB model to a strategy using a reformulation (called REPB formulation). The REPB scheme captures the quasi-neutral limit more accurately

Robust and conservative dynamical low-rank methods for the Vlasov equation via a novel macro-micro decomposition

Dynamical low-rank (DLR) approximation has gained interest in recent years as a viable solution to the curse of dimensionality in the numerical solution of kinetic equations including the Boltzmann and Vlasov equations. These methods include the projector-splitting and Basis-update & Galerkin (BUG) DLR integrators, and have shown promise at greatly improving the computational efficiency of kinetic solutions. However, this often comes at the cost of conservation of charge, current and energy. In this work we show how a novel macro-micro decomposition may be used to separate the distribution function into two components, one of which carries the conserved quantities, and the other of which is orthogonal to them. We apply DLR approximation to the latter, and thereby achieve a clean and extensible approach to a conservative DLR scheme which retains the computational advantages of the base scheme. Moreover, our approach requires no change to the mechanics of the DLR approximation, so it is compatible with both the BUG family of integrators and the projector-splitting integrator which we use here. We describe a first-order integrator which can exactly conserve charge and either current or energy, as well as an integrator which exactly conserves charge and energy and exhibits second-order accuracy on our test problems. To highlight the flexibility of the proposed macro-micro decomposition, we implement a pair of velocity space discretizations, and verify the claimed accuracy and conservation properties on a suite of plasma benchmark problems.Comment: 33 pages, 6 figure

2-D Magnetohydrodynamic Simulations of Induced Plasma Dynamics in the Near-Core Region of a Galaxy Cluster

We present results from numerical simulations of the cooling-core cluster A2199 produced by the two-dimensional (2-D) resistive magnetohydrodynamics (MHD) code MACH2. In our simulations we explore the effect of anisotropic thermal conduction on the energy balance of the system. The results from idealized cases in 2-D axisymmetric geometry underscore the importance of the initial plasma density in ICM simulations, especially the near-core values since the radiation cooling rate is proportional to ${n_e}^2$. Heat conduction is found to be non-effective in preventing catastrophic cooling in this cluster. In addition we performed 2-D planar MHD simulations starting from initial conditions deliberately violating both thermal balance and hydrostatic equilibrium in the ICM, to assess contributions of the convective terms in the energy balance of the system against anisotropic thermal conduction. We find that in this case work done by the pressure on the plasma can dominate the early evolution of the internal energy over anisotropic thermal conduction in the presence of subsonic flows, thereby reducing the impact of the magnetic field. Deviations from hydrostatic equilibrium near the cluster core may be associated with transient activity of a central active galactic nucleus and/or remnant dynamical activity in the ICM and warrant further study in three dimensions.Comment: 16 pages, 13 figures, accepted for publication in MNRA

Sheared Flow As A Stabilizing Mechanism In Astrophysical Jets

It has been hypothesized that the sustained narrowness observed in the asymptotic cylindrical region of bipolar outflows from Young Stellar Objects (YSO) indicates that these jets are magnetically collimated. The j cross B force observed in z-pinch plasmas is a possible explanation for these observations. However, z-pinch plasmas are subject to current driven instabilities (CDI). The interest in using z-pinches for controlled nuclear fusion has lead to an extensive theory of the stability of magnetically confined plasmas. Analytical, numerical, and experimental evidence from this field suggest that sheared flow in magnetized plasmas can reduce the growth rates of the sausage and kink instabilities. Here we propose the hypothesis that sheared helical flow can exert a similar stabilizing influence on CDI in YSO jets.Comment: 13 pages, 2 figure

Whole Device Modeling of the FuZE Sheared-Flow-Stabilized Z Pinch

The FuZE sheared-flow-stabilized Z pinch at Zap Energy is simulated using whole-device modeling employing an axisymmetric resistive magnetohydrodynamic formulation implemented within the discontinuous Galerkin WARPXM framework. Simulations show formation of Z pinches with densities of approximately 10^22 m^-3 and total DD fusion neutron rate of 10^7 per {\mu}s for approximately 2 {\mu}s. Simulation-derived synthetic diagnostics show peak currents and voltages within 10% and total yield within approximately 30% of experiment for similar plasma mass. The simulations provide insight into the plasma dynamics in the experiment and enable a predictive capability for exploring design changes on devices built at Zap Energy.Comment: 8 pages, 9 figures, IAEA FEC 202

Numerical approximation of the Euler-Maxwell model in the quasineutral limit

International audienceWe derive and analyze an Asymptotic-Preserving scheme for the Euler-Maxwell system in the quasi-neutral limit. We prove that the linear stability condition on the time-step is independent of the scaled Debye length $\lambda$ when $\lambda \to 0$. Numerical validation performed on Riemann initial data and for a model Plasma Opening Switch device show that the AP-scheme is convergent to the Euler-Maxwell solution when $\Delta x/ \lambda \to 0$ where $\Delta x$ is the spatial discretization. But, when $\lambda /\Delta x \to 0$, the AP-scheme is consistent with the quasi-neutral Euler-Maxwell system. The scheme is also perfectly consistent with the Gauss equation. The possibility of using large time and space steps leads to several orders of magnitude reductions in computer time and storage

The Kadomtsev pinch revisited for sheared-flow-stabilized Z-pinch modeling

The Kadomtsev pinch, namely the Z-pinch profile marginally stable to interchange modes, is revisited in light of observations from axisymmetric MHD modeling of the FuZE sheared-flow-stabilized Z-pinch experiment. We show that Kadomtsev's stability criterion, cleanly derived by the minimum energy principle but of opaque physical significance, has an intuitive interpretation in the specific entropy analogous to the Schwarzschild-Ledoux criterion for convective stability of adiabatic pressure distributions in the fields of astrophysics, meteorology, and oceanography. By analogy, the Kadomtsev profile may be described as magnetoadiabatic in the sense that plasma pressure is polytropically related to area-averaged current density from the ideal MHD stability condition on the specific entropy. Further, the non-ideal stability condition of the entropy modes is shown to relate the specific entropy gradient to the ideal interchange stability function. Hence, the combined activity of the ideal interchange and non-ideal entropy modes drives both the specific entropy and specific magnetic flux gradients to zero in the marginally stable state. The physical properties of Kadomtsev's pinch are reviewed in detail and following from this the localization of pinch confinement, i.e., pinch size and inductance, is quantified by the ratio of extensive magnetic and thermal energies. In addition, results and analysis of axisymmetric MHD modeling of the FuZE Z-pinch experiment are presented where pinch structure is found to consist of a near-marginal flowing core surrounded by a super-magnetoadiabatic low-beta sheared flow.Comment: Author version of accepted article for IEEE Transactions on Plasma Scienc