Time-discretization of a plasma-neutral MHD model with a semi-implicit leapfrog algorithm

DOI10.1016/J.CPC.2022.108288arXiv2110.02328OpenAlexW3203631010MaRDI QIDQ6159630FDOQ6159630

Authors: Sina Taheri, Jacob R. King, U. Shumlak 

Publication date: 20 June 2023

Published in: Computer Physics Communications (Search for Journal in Brave)

Abstract: The semi-implicit leapfrog time-discretization is a workhorse algorithm for initial-value MHD codes to bridge between vastly separated time scales. Inclusion of atomic interactions with neutrals breaks the functional structure of the MHD equations that exploited by the leapfrog. We address how to best integrate atomic physics into the semi-implicit leapfrog. Following the Crank-Nicolson method, one approach is to time-center the atomic interactions in the linear solver and use a Newton method to include the nonlinear contributions. Alternatively, another family of methods are based on operator-splitting the terms associated with the atomic interactions using a Strang-splitting technique. These methods naturally break equations into constituent ODE and PDE parts and preserve the structure exploited by the semi-implicit leapfrog. We study the accuracy and efficiency of these methods through a battery of 0D and 1D cases and show that a second-order-in-time Douglas-Rachford inspired coupling between the ODE and PDE advances is effective in reducing the time-discretization error to be comparable to that of Crank-Nicolson with Newton iteration of the nonlinear terms. Splitting ODE and PDE parts results in independent matrix solves for each field which reduces the computational cost considerably and provides parallelization over species relative to Crank-Nicolson.

Full work available at URL: https://arxiv.org/abs/2110.02328

Mathematics Subject Classification ID

Numerical analysis (65-XX) Partial differential equations (35-XX)

Cites Work

• On the Numerical Solution of Heat Conduction Problems in Two and Three Space Variables
• Finite Volume Methods for Hyperbolic Problems
• Solving Nonlinear Equations with Newton's Method
• Numerical analysis of time integration errors for nonequilibrium radiation diffusion
• On the Construction and Comparison of Difference Schemes
• Jacobian-free Newton-Krylov methods: a survey of approaches and applications.
• Nonlinear magnetohydrodynamics simulation using high-order finite elements.
• Lectures in magnetohydrodynamics. With an appendix on extended MHD
• Spectral element spatial discretization error in solving highly anisotropic heat conduction equation
• Analysis of a mixed semi-implicit/implicit algorithm for low-frequency two-fluid plasma modeling
• Semi-implicit magnetohydrodynamic calculations
• A class of basis functions for non-ideal magnetohydrodynamic computations

Cited In (3)

• An unconditionally stable modified leapfrog method for Maxwell's equation in Kerr-type nonlinear media
• An analysis and successful benchmarking of the Chapman-Enskog-like (CEL) continuum kinetic closure approach algorithm in NIMROD
• An indirect ALE discretization of single fluid plasma without a fast magnetosonic time step restriction

This page was built for publication: Time-discretization of a plasma-neutral MHD model with a semi-implicit leapfrog algorithm

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6159630)