skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

Journal Article · · Journal of Computational Physics
ORCiD logo [1];  [2];  [3];  [3];  [4]
  1. Clemson Univ., SC (United States)
  2. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Univ. of New Mexico, Albuquerque, NM (United States)
  3. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
  4. TU Dortmund Univ. (Germany)

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

Research Organization:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); German Research Foundation (DFG)
Grant/Contract Number:
AC04-94AL85000; NA0003525; KU 1530/15-2
OSTI ID:
1618098
Alternate ID(s):
OSTI ID: 1605331
Report Number(s):
SAND-2019-3661J; 674276; TRN: US2106794
Journal Information:
Journal of Computational Physics, Vol. 410, Issue C; ISSN 0021-9991
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 10 works
Citation information provided by
Web of Science

References (34)

Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics journal December 2001
High-resolution FEM–FCT schemes for multidimensional conservation laws journal November 2004
Locally Divergence-Free Discontinuous Galerkin Methods for MHD Equations journal June 2005
Local bounds preserving stabilization for continuous Galerkin discretization of hyperbolic systems journal May 2018
Flux Correction Tools for Finite Elements journal January 2002
Notes on the Eigensystem of Magnetohydrodynamics journal February 1996
A Staggered Mesh Algorithm Using High Order Godunov Fluxes to Ensure Solenoidal Magnetic Fields in Magnetohydrodynamic Simulations journal March 1999
Numerical magetohydrodynamics in astrophysics: Algorithm and tests for one-dimensional flow` journal March 1995
Positivity-Preserving Analysis of Numerical Schemes for Ideal Magnetohydrodynamics journal January 2018
Monolithic convex limiting for continuous finite element discretizations of hyperbolic conservation laws journal April 2020
Gradient-based nodal limiters for artificial diffusion operators in finite element schemes for transport equations: Gradient-based nodal limiters for artificial diffusion operators in finite element schemes for transport equations journal February 2017
Hyperbolic Divergence Cleaning for the MHD Equations journal January 2002
The group finite element formulation journal April 1983
Limiting and divergence cleaning for continuous finite element discretizations of the MHD equations journal April 2020
Small-scale structure of two-dimensional magnetohydrodynamic turbulence journal January 1979
Sequential limiting in continuous and discontinuous Galerkin methods for the Euler equations journal March 2018
Convergence of locally divergence-free discontinuous-Galerkin methods for the induction equations of the 2D-MHD system journal November 2005
Fully multidimensional flux-corrected transport algorithms for fluids journal June 1979
A Provably Positive Discontinuous Galerkin Method for Multidimensional Ideal Magnetohydrodynamics journal January 2018
Second-Order Invariant Domain Preserving Approximation of the Euler Equations Using Convex Limiting journal January 2018
A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure journal July 2016
Efficient MHD Riemann solvers for simulations on unstructured triangular grids journal January 2002
A second-order unsplit Godunov scheme for cell-centered MHD: The CTU-GLM scheme journal March 2010
Differentiable monotonicity-preserving schemes for discontinuous Galerkin methods on arbitrary meshes journal June 2017
High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics journal November 2017
Monotonicity-preserving finite element schemes based on differentiable nonlinear stabilization journal January 2017
A Higher-Order Godunov Method for Multidimensional Ideal Magnetohydrodynamics journal March 1994
Comparison of Some Flux Corrected Transport and Total Variation Diminishing Numerical Schemes for Hydrodynamic and Magnetohydrodynamic Problems journal October 1996
Provably positive high-order schemes for ideal magnetohydrodynamics: analysis on general meshes journal May 2019
An upwind differencing scheme for the equations of ideal magnetohydrodynamics journal April 1988
An FCT finite element scheme for ideal MHD equations in 1D and 2D journal June 2017
Linearity-preserving monotone local projection stabilization schemes for continuous finite elements journal August 2017
The ∇·B=0 Constraint in Shock-Capturing Magnetohydrodynamics Codes journal July 2000
Edge-based nonlinear diffusion for finite element approximations of convection–diffusion equations and its relation to algebraic flux-correction schemes journal May 2016

Cited By (1)

On differentiable local bounds preserving stabilization for Euler equations text January 2019