Abstract:This paper proposes and analyzes a novel efficient high-order finite volume method for the ideal magnetohydrodynamics (MHD). As a distinctive feature, the method simultaneously preserves a discretely divergence-free (DDF) constraint on the magnetic field and the positivity-preserving (PP) property, which ensures the positivity of density, pressure, and internal energy. To enforce the DDF condition, we design a new discrete projection approach that projects the reconstructed point values at the cell interface into a DDF space, without using any approximation polynomials. This projection method is highly efficient, easy to implement, and particularly suitable for standard high-order finite volume WENO methods, which typically return only the point values in the reconstruction. Moreover, we also develop a new finite volume framework for constructing provably PP schemes for the ideal MHD system. The framework comprises the discrete projection technique, a suitable approximation to the Godunov--Powell source terms, and a simple PP limiter. We provide rigorous analysis of the PP property of the proposed finite volume method, demonstrating that the DDF condition and the proper approximation to the source terms eliminate the impact of magnetic divergence terms on the PP property. The analysis is challenging due to the internal energy function's nonlinearity and the intricate relationship between the DDF and PP properties. To address these challenges, the recently developed geometric quasilinearization approach is adopted, which transforms a nonlinear constraint into a family of linear constraints. Finally, we validate the effectiveness of the proposed method through several benchmark and demanding numerical examples. The results demonstrate that the proposed method is robust, accurate, and highly effective, confirming the significance of the proposed DDF projection and PP techniques.

Fundamental Higher Order Weaker HIE-PML Algorithm for Magnetized Plasma

One-step Leapfrog ADI-FDTD Method for Anisotropic Magnetized Plasma

Perfectly Matched Layer for Accurate FDTD for Anisotropic Magnetized Plasma

Combined Piecewise Linear Recursive Convolution-Bilinear Transform Implementation of the CFS-PML for Unmagnetized Plasma

Improved HIE-FDTD Method for Simulating Thin Layer Structure in Linear Dispersive Media

Efficient Implementation of the CPML in 3-D Hybrid Implicit-Explicit FDTD Method

Complex Envelope Weaker HIE with Higher Order PML Algorithm for FDTD Simulation

High-order Implicit Particle-in-cell Method for Plasma Simulations at Solid Densities.

A mimetic finite difference based quasi-static magnetohydrodynamic solver for force-free plasmas in tokamak disruptions

Three-Dimensional Direct-Implicit Particle-in-Cell Model Using Trilinear Anisotropic Immersed-Finite-Element for Plasma Propulsion

Tridiagonalized WLP-FDTD Implementation for Wave Propagation in Magnetized Plasma

Accurate and Robust Surface MLFMA Solution for Scattering from Large Plasma Object with Large Negative Permittivity

A new discretely divergence-free positivity-preserving high-order finite volume method for ideal MHD equations

Hybrid fluid–particle modeling of shock-driven hydrodynamic instabilities in a plasma

Multi-Physics Mathematical Model for Weakly-Ionized Plasma Flows

Geometric Particle-In-Cell discretizations of a plasma hybrid model with kinetic ions and mass-less fluid electrons

Nonlinear magnetohydrodynamics simulation using high-order finite elements

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

An Extended Hybrid Magnetohydrodynamics Gyrokinetic Model For Numerical Simulation Of Shear Alfven Waves In Burning Plasmas

Nonlinear extended magnetohydrodynamics simulation using high-order finite elements

Magnetohydrodynamic with Adaptively Embedded Particle-in-Cell model: MHD-AEPIC