4.6 Article

Globally Divergence-Free Discontinuous Galerkin Methods for Ideal Magnetohydrodynamic Equations

Journal

JOURNAL OF SCIENTIFIC COMPUTING
Volume 77, Issue 3, Pages 1621-1659

Publisher

SPRINGER/PLENUM PUBLISHERS
DOI: 10.1007/s10915-018-0750-6

Keywords

MHD equations; Divergence-free magnetic field; Discontinuous Galerkin methods; H(div)-conforming finite element spaces; Fourier analysis

Funding

  1. NSF [DMS-1318409, DMS-1719942]
  2. NSFC [11722112, 91630207]

Ask authors/readers for more resources

Ideal magnetohydrodynamic (MHD) equations are widely used in many areas in physics and engineering, and these equations have a divergence-free constraint on the magnetic field. In this paper, we propose high order globally divergence-free numerical methods to solve the ideal MHD equations. The algorithms are based on discontinuous Galerkin methods in space. The induction equation is discretized separately to approximate the normal components of the magnetic field on elements interfaces, and to extract additional information about the magnetic field when higher order accuracy is desired. This is then followed by an element by element reconstruction to obtain the globally divergence-free magnetic field. In time, strong-stability-preserving Runge-Kutta methods are applied. In consideration of accuracy and stability of the methods, a careful investigation is carried out, both numerically and analytically, to study the choices of the numerical fluxes associated with the electric field at element interfaces and vertices. The resulting methods are local and the approximated magnetic fields are globally divergence-free. Numerical examples are presented to demonstrate the accuracy and robustness of the methods.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.6
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Mathematics, Applied

Asymptotically compatible approximations of linear nonlocal conservation laws with variable horizon

Xiaoxuan Yu, Yan Xu, Qiang Du

Summary: The study focuses on a quadrature-based finite difference discretization of one-dimensional scalar linear nonlocal conservation laws with spatially varying range of nonlocal interactions. It addresses the convergence of the discrete approximation in both nonlocal and local settings, presenting the first complete proof of numerical discretization convergence to feasible kernels in the nonlocal regime and local limit, establishing the asymptotic compatibility of the scheme. Numerical results are also provided to illustrate the impact of variable horizon on wave propagation in the nonlocal model.

NUMERICAL METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS (2022)

Article Physics, Mathematical

High-Order Positivity-Preserving Well-Balanced Discontinuous Galerkin Methods for Euler Equations with Gravitation on Unstructured Meshes

Weijie Zhang, Yulong Xing, Yinhua Xia, Yan Xu

Summary: This paper proposes a high-order accurate DG method for the compressible Euler equations on unstructured meshes under gravitational fields, which preserves a general hydrostatic equilibrium state and guarantees the positivity of density and pressure. Through a special way to recover the equilibrium state and the design of novel interface variables, the scheme achieves well-balanced and positivity-preserving properties.

COMMUNICATIONS IN COMPUTATIONAL PHYSICS (2022)

Article Physics, Mathematical

High Order Finite Difference WENO Methods with Unequal-Sized Sub-Stencils for the Degasperis-Procesi Type Equations

Jianfang Lin, Yan Xu, Huiwen Xue, Xinghui Zhong

Summary: In this paper, two finite difference weighted essentially non-oscillatory (WENO) schemes with unequal-sized sub-stencils are developed for solving the Degasperis-Procesi (DP) and mu-Degasperis-Procesi (mu DP) equations. The proposed schemes are simple in stencil choice and allow for arbitrary positive linear weights. Numerical tests demonstrate the high order accuracy and non-oscillatory properties of the schemes.

COMMUNICATIONS IN COMPUTATIONAL PHYSICS (2022)

Article Mathematics, Applied

Numerical Simulation of Singularity Propagation Modeled by Linear Convection Equations with Spatially Heterogeneous Nonlocal Interactions

Xiaoxuan Yu, Yan Xu, Qiang Du

Summary: We study the propagation of singularities in solutions of linear convection equations with spatially heterogeneous nonlocal interactions. We are interested in understanding the impact on singularity propagation due to the heterogeneities of the nonlocal horizon and the local and nonlocal transition. We derive equations to characterize the propagation of different types of singularities for various forms of nonlocal horizon parameters in the nonlocal regime, and use numerical simulations to illustrate the propagation patterns in different scenarios.

JOURNAL OF SCIENTIFIC COMPUTING (2022)

Article Mathematics, Applied

A Parallel Eigensolver for Photonic Crystals Discretized by Edge Finite Elements

Zhongjie Lu, Yan Xu

Summary: This paper introduces a method to compute the discrete Maxwell eigenproblem by discretizing and constructing a penalty term to deal with its huge kernel. The algorithm is efficient and capable as demonstrated by numerical examples.

JOURNAL OF SCIENTIFIC COMPUTING (2022)

Article Mathematics, Applied

Accuracy-Enhancement of Discontinuous Galerkin Methods for PDEs Containing High Order Spatial Derivatives

Qi Tao, Liangyue Ji, Jennifer K. Ryan, Yan Xu

Summary: In this paper, the authors discuss the enhancement of accuracy in discontinuous Galerkin methods for solving PDEs with high order spatial derivatives. They introduce the use of a Smoothness-Increasing Accuracy-Conserving (SIAC) filter to create a superconvergence filtered solution. The authors provide theoretical proofs for the accuracy improvement and demonstrate the applicability of the SIAC filter for both linear and nonlinear higher order PDEs.

JOURNAL OF SCIENTIFIC COMPUTING (2022)

Article Computer Science, Interdisciplinary Applications

Structure-preserving finite volume arbitrary Lagrangian-Eulerian WENO schemes for the shallow water equations

Jiahui Zhang, Yinhua Xia, Yan Xu

Summary: This paper develops ALE-WENO schemes for the shallow water equations under the ALE framework, adopting the WENO hybrid reconstruction on moving meshes. The schemes demonstrate high order accuracy and positivity-preserving property, as well as the preservation of equilibrium and capturing of small perturbations without numerical oscillations. The ALE-WENO hybrid schemes have an advantage over static mesh simulations due to higher resolution interface tracking.

JOURNAL OF COMPUTATIONAL PHYSICS (2023)

Article Mathematics, Applied

Error estimates of the local discontinuous Galerkin methods for two-dimensional (μ)-Camassa-Holm equations

Jinyang Lu, Yan Xu, Chao Zhang

Summary: This paper presents a uniform framework of local discontinuous Galerkin (LDG) methods for two-dimensional Camassa-Holm equations and two-dimensional mu-Camassa-Holm equations. The energy stability and semi-discrete error estimates are derived based on this framework. Numerical experiments demonstrate the accuracy and stability of these schemes.

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (2023)

Article Mathematics, Applied

Moving Water Equilibria Preserving Discontinuous Galerkin Method for the Shallow Water Equations

Jiahui Zhang, Yinhua Xia, Yan Xu

Summary: In this paper, a new well-balanced discontinuous Galerkin (DG) method is proposed to preserve moving-water equilibria in the shallow water equations. Instead of approximating the conservative variables, the scheme approximates the equilibrium variables in the DG piecewise polynomial space. The numerical fluxes are modified based on the generalized hydrostatic reconstruction to achieve moving water equilibrium preservation.

JOURNAL OF SCIENTIFIC COMPUTING (2023)

Article Mathematics, Applied

A discontinuous Galerkin method for the Camassa-Holm-Kadomtsev-Petviashvili type equations

Qian Zhang, Yan Xu, Yue Liu

Summary: This paper develops a high-order discontinuous Galerkin (DG) method for the Camassa-Holm-Kadomtsev-Petviashvili (CH-KP) type equations on Cartesian meshes. The proposed method efficiently deals with the integration operator partial differential delta(-1) element by element and is applicable to most solutions. The DG scheme is proven to be an energy stable numerical scheme using the instinctive energy of the original PDE as a guiding principle, and the semi-discrete error estimates are derived for the nonlinear case without any priori assumption. Several numerical experiments demonstrate the capability of the proposed schemes for various types of solutions.

NUMERICAL METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS (2023)

Article Mathematics, Applied

High-Order Decoupled and Bound Preserving Local Discontinuous Galerkin Methods for a Class of Chemotaxis Models

Wei Zheng, Yan Xu

Summary: In this paper, bound preserving and high-order accurate LDG schemes are used to solve chemotaxis models. First-order temporal accurate schemes are constructed using the gradient flow structure of the models, and can be extended to high accuracy schemes. Bound preserving is achieved by using Lagrange multipliers and an efficient active set semi-smooth Newton method.

COMMUNICATIONS ON APPLIED MATHEMATICS AND COMPUTATION (2023)

Article Mathematics, Applied

High order structure-preserving arbitrary Lagrangian-Eulerian discontinuous Galerkin methods for the Euler equations under gravitational fields

Weijie Zhang, Yulong Xing, Yinhua Xia, Yan Xu

Summary: In this study, high-order arbitrary Lagrangian-Eulerian discontinuous Galerkin (ALE-DG) methods for the Euler equations with gravitational fields on moving meshes are proposed. The goal is to show that ALE-DG methods can achieve structure-preserving properties of DG methods for the Euler equations on arbitrary moving meshes. Two well-balanced and positivity-preserving ALE-DG schemes are presented, which can preserve the explicitly given equilibrium state on any moving grids. The schemes are established in both one and two dimensions on unstructured triangular meshes.

COMPUTERS & MATHEMATICS WITH APPLICATIONS (2023)

Article Mathematics, Applied

SEMI-IMPLICIT SPECTRAL DEFERRED CORRECTION METHODS BASED ON SECOND-ORDER TIME INTEGRATION SCHEMES FOR NONLINEAR PDES

Ruihan Guo, Yan Xu

Summary: In this paper, we propose a class of semi-implicit spectral deferred correction (SDC) methods based on second-order time integration methods, where the order of accuracy is increased by two for each additional iteration. The spatial discretization is done using the local discontinuous Galerkin (LDG) method to obtain fully-discrete schemes that are high-order accurate in both space and time. Numerical experiments are conducted to demonstrate the accuracy, efficiency, and robustness of the proposed semi-implicit SDC methods in solving complex nonlinear PDEs.

JOURNAL OF COMPUTATIONAL MATHEMATICS (2023)

Article Engineering, Mechanical

Adaptive local discontinuous Galerkin methods with semi-implicit time discretizations for the Navier-Stokes equations

Xiangyi Meng, Yan Xu

Summary: In this paper, a mesh adaptation algorithm for the unsteady compressible Navier-Stokes equations is presented. The algorithm combines local discontinuous Galerkin methods with implicit-explicit Runge-Kutta or spectral deferred correction time discretization methods. The algorithm shows high accuracy, efficiency, and capability in capturing flow structures through numerical experiments.

ADVANCES IN AERODYNAMICS (2022)

Article Mathematics, Applied

Negative Norm Estimates for Arbitrary Lagrangian-Eulerian Discontinuous Galerkin Method for Nonlinear Hyperbolic Equations

Qi Tao, Yan Xu, Xiaozhou Li

Summary: In this paper, the important of negative norm estimates for the arbitrary Lagrangian-Eulerian discontinuous Galerkin (ALE-DG) method solving nonlinear hyperbolic equations with smooth solutions is presented. The smoothness-increasing accuracy-conserving (SIAC) filter is used as a post-processing technique to enhance the accuracy of the discontinuous Galerkin (DG) solutions. This work is a crucial step in extending the SIAC filter to moving mesh for nonlinear problems.

COMMUNICATIONS ON APPLIED MATHEMATICS AND COMPUTATION (2022)

No Data Available