Positivity-Preserving Analysis of Numerical Schemes for Ideal Magnetohydrodynamics

DOI10.1137/18M1168017zbMATH Open1391.76369arXiv1802.02278WikidataQ129501635 ScholiaQ129501635MaRDI QIDQ3174822FDOQ3174822

Kailiang Wu

Publication date: 18 July 2018

Published in: SIAM Journal on Numerical Analysis (Search for Journal in Brave)

Abstract: Numerical schemes provably preserving the positivity of density and pressure are highly desirable for MHD, but the rigorous positivity-preserving (PP) analysis remains challenging. The difficulties mainly arise from the intrinsic complexity of the MHD equations as well as the indeterminate relation between the PP property and the divergence-free condition on magnetic field. We present the first rigorous PP analysis of conservative schemes with Lax-Friedrichs (LF) flux for ideal MHD. The significant innovation is the discovery of theoretical connection between PP property and a discrete divergence-free (DDF) condition. This connection is established through the generalized LF splitting properties, which are alternatives of the usually-expected LF splitting property that does not hold for ideal MHD. The generalized LF splitting properties involve a number of admissible states strongly coupled by DDF condition, making their derivation very difficult. We derive these properties via a novel equivalent form of the admissible state set and an important inequality skillfully constructed by technical estimates. Rigorous PP analysis is presented for finite volume and discontinuous Galerkin schemes with LF flux on uniform Cartesian meshes. In 1D case, PP property is proved for the first-order scheme with proper numerical viscosity, and also for arbitrarily high-order schemes under conditions accessible by a PP limiter. In 2D case, we show that the DDF condition is necessary and crucial for achieving PP property. It is observed that even slightly violating the proposed DDF condition may cause failure to preserve the positivity of pressure. We prove that the 2D LF type scheme with proper numerical viscosity preserves the positivity and DDF condition. Sufficient conditions are derived for 2D PP high-order schemes, and extension to 3D is discussed. Numerical examples confirm the theoretical findings.

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

zbMATH Keywords

hyperbolic conservation lawspositivity-preservingcompressible magnetohydrodynamicsdiscrete divergence-free conditiongeneralized Lax-Friedrichs splittingadmissible states

Mathematics Subject Classification ID

PDEs in connection with fluid mechanics (35Q35) Hyperbolic conservation laws (35L65) Finite volume methods for initial value and initial-boundary value problems involving PDEs (65M08) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Finite element methods applied to problems in fluid mechanics (76M10) Magnetohydrodynamics and electrohydrodynamics (76W05) Finite volume methods applied to problems in fluid mechanics (76M12)

Cites Work

• Maximum-principle-satisfying and positivity-preserving high-order schemes for conservation laws: survey and new developments
• Positivity-preserving method for high-order conservative schemes solving compressible Euler equations
• Strong stability-preserving high-order time discretization methods
• On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes
• On maximum-principle-satisfying high order schemes for scalar conservation laws
• High order strong stability preserving time discretizations
• A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations
• A positive conservative method for magnetohydrodynamics based on HLL and Roe methods
• The \(\nabla \cdot B=0\) constraint in shock-capturing magnetohydrodynamics codes
• Parametrized maximum principle preserving flux limiters for high order schemes solving multi-dimensional scalar hyperbolic conservation laws
• Parametrized positivity preserving flux limiters for the high order finite difference WENO scheme solving compressible Euler equations
• Parametrized maximum principle preserving flux limiters for high order schemes solving hyperbolic conservation laws: one-dimensional scalar problem
• A multiwave approximate Riemann solver for ideal MHD based on relaxation. II: Numerical implementation with 3 and 5 waves
• Hyperbolic divergence cleaning for the MHD equations
• High-order accurate physical-constraints-preserving finite difference WENO schemes for special relativistic hydrodynamics
• Arbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
• Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field
• The effect of nonzero \(\bigtriangledown\cdot B\) on the numerical solution of the magnetohydrodynamic equations
• An Unstaggered, High‐Resolution Constrained Transport Method for Magnetohydrodynamic Flows
• Locally divergence-free discontinuous Galerkin methods for MHD equations
• Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics
• Maintaining pressure positivity in magnetohydrodynamic simulations
• On the divergence-free condition in Godunov-type schemes for ideal magnetohydrodynamics: the upwind constrained transport method.
• Positivity-preserving DG and central DG methods for ideal MHD equations
• A multiwave approximate Riemann solver for ideal MHD based on relaxation. I: Theoretical framework
• Finite difference weighted essentially non-oscillatory schemes with constrained transport for ideal magnetohydrodynamics
• Locally Divergence-preserving Upwind Finite Volume Schemes for Magnetohydrodynamic Equations
• On high order strong stability preserving Runge-Kutta and multi step time discretizations
• Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms
• High order central scheme on overlapping cells for magneto-hydrodynamic flows with and without constrained transport method
• A positive MUSCL-Hancock scheme for ideal magnetohydrodynamics
• On positivity-preserving high order discontinuous Galerkin schemes for compressible Navier-Stokes equations
• Constraint-Preserving Upwind Methods for Multidimensional Advection Equations
• Entropy stable finite volume scheme for ideal compressible MHD on 2-D Cartesian meshes
• Bound-Preserving High-Order Schemes
• Positivity-Preserving Finite Difference Weighted ENO Schemes with Constrained Transport for Ideal Magnetohydrodynamic Equations
• Admissible states and physical-constraints-preserving schemes for relativistic magnetohydrodynamic equations
• Fast estimation from above of the maximum wave speed in the Riemann problem for the Euler equations
• A Provably Positive Discontinuous Galerkin Method for Multidimensional Ideal Magnetohydrodynamics
• Increasing the accuracy in locally divergence-preserving finite volume schemes for MHD
• Positivity-preserving discontinuous Galerkin methods with Lax-Wendroff time discretizations
• A high-order positivity-preserving single-stage single-step method for the ideal magnetohydrodynamic equations

Cited In (45)

• Minimum Principle on Specific Entropy and High-Order Accurate Invariant-Region-Preserving Numerical Methods for Relativistic Hydrodynamics
• Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes
• Implementation of the HLL-GRP solver for multidimensional ideal MHD simulations based on finite volume method
• The finite volume scheme preserving maximum principle for diffusion equations with discontinuous coefficient
• A new locally divergence-free path-conservative central-upwind scheme for ideal and shallow water magnetohydrodynamics
• Provably physical-constraint-preserving discontinuous Galerkin methods for multidimensional relativistic MHD equations
• On high order positivity-preserving well-balanced finite volume methods for the Euler equations with gravitation
• Positivity-Preserving Finite Difference Weighted ENO Schemes with Constrained Transport for Ideal Magnetohydrodynamic Equations
• Uniformly High-Order Structure-Preserving Discontinuous Galerkin Methods for Euler Equations with Gravitation: Positivity and Well-Balancedness
• High order asymptotic preserving finite difference WENO schemes with constrained transport for MHD equations in all sonic Mach numbers
• GQL-based bound-preserving and locally divergence-free central discontinuous Galerkin schemes for relativistic magnetohydrodynamics
• Unconditional optimal first-order error estimates of a full pressure segregation scheme for the magnetohydrodynamics equations
• On physical-constraints-preserving schemes for special relativistic magnetohydrodynamics with a general equation of state
• A New Discretely Divergence-Free Positivity-Preserving High-Order Finite Volume Method for Ideal MHD Equations
• Provably Positive Central Discontinuous Galerkin Schemes via Geometric Quasilinearization for Ideal MHD Equations
• High-Order Accurate Entropy Stable Finite Difference Schemes for One- and Two-Dimensional Special Relativistic Hydrodynamics
• Small collaboration: Advanced numerical methods for nonlinear hyperbolic balance laws and their applications. Abstracts from the small collaboration held August 29 -- September 4, 2021 (hybrid meeting)
• Positivity-preserving well-balanced central discontinuous Galerkin schemes for the Euler equations under gravitational fields
• A Provably Positive Discontinuous Galerkin Method for Multidimensional Ideal Magnetohydrodynamics
• A physical-constraint-preserving finite volume WENO method for special relativistic hydrodynamics on unstructured meshes
• Provably convergent Newton-Raphson methods for recovering primitive variables with applications to physical-constraint-preserving Hermite WENO schemes for relativistic hydrodynamics
• Geometric Quasilinearization Framework for Analysis and Design of Bound-Preserving Schemes
• Asymptotic-preserving discretization of three-dimensional plasma fluid models
• Bound-preserving framework for central-upwind schemes for general hyperbolic conservation laws
• Physical-constraints-preserving Lagrangian finite volume schemes for one- and two-dimensional special relativistic hydrodynamics
• Is the classic convex decomposition optimal for bound-preserving schemes in multiple dimensions?
• On Optimal Cell Average Decomposition for High-Order Bound-Preserving Schemes of Hyperbolic Conservation Laws
• Positivity-Preserving Numerical Schemes for Lubrication-Type Equations
• High-order accurate positivity-preserving and well-balanced discontinuous Galerkin schemes for ten-moment Gaussian closure equations with source terms
• Error analysis of vector penalty-projection method with second order accuracy for incompressible magnetohydrodynamic system
• A positivity-preserving Lagrangian discontinuous Galerkin method for ideal magnetohydrodynamics equations in one-dimension
• A class of positive semi-discrete Lagrangian-Eulerian schemes for multidimensional systems of hyperbolic conservation laws
• A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations
• Bound-preserving OEDG schemes for Aw-Rascle-Zhang traffic models on networks
• A class of bound-preserving high order schemes: The main ideas and recent developments
• Entropy-stable Gauss collocation methods for ideal magneto-hydrodynamics
• Uniformly high-order bound-preserving OEDG schemes for two-phase flows
• Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations
• A fully decoupled linearized finite element method with second-order temporal accuracy and unconditional energy stability for incompressible MHD equations
• A flux-differencing formula for split-form summation by parts discretizations of non-conservative systems. Applications to subcell limiting for magneto-hydrodynamics
• A provable positivity-preserving local discontinuous Galerkin method for the viscous and resistive MHD equations
• A WENO-Based Stochastic Galerkin Scheme for Ideal MHD Equations with Random Inputs
• A new locally divergence-free WLS-ENO scheme based on the positivity-preserving finite volume method for ideal MHD equations
• Constraint preserving discontinuous Galerkin method for ideal compressible MHD on 2-D Cartesian grids
• Limiting and divergence cleaning for continuous finite element discretizations of the MHD equations

This page was built for publication: Positivity-Preserving Analysis of Numerical Schemes for Ideal Magnetohydrodynamics

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