A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations
- Clemson Univ., SC (United States)
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Univ. of New Mexico, Albuquerque, NM (United States)
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- 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
Web of Science
Similar Records
Matrix-free subcell residual distribution for Bernstein finite element discretizations of linear advection equations
Sequential limiting in continuous and discontinuous Galerkin methods for the Euler equations