High order well-balanced and total-energy-conserving local discontinuous Galerkin methods for compressible self-gravitating Euler equations

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Pan, Liang, Chen, Wei, Qiu, Jianxian, Xiong, Tao
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909825550516224
author Pan, Liang
Chen, Wei
Qiu, Jianxian
Xiong, Tao
author_facet Pan, Liang
Chen, Wei
Qiu, Jianxian
Xiong, Tao
contents In this paper, we develop a high order structure-preserving local discontinuous Galerkin (DG) scheme for the compressible self-gravitating Euler equations, which pose great challenges due to the presence of time-dependent gravitational potential. The designed scheme is well-balanced for general polytropic equilibrium state and total energy conserving for multiple spatial dimensions without an assumption of spherical symmetry. The well-balanced property is achieved by decomposing the gravitational potential into equilibrium and perturbation parts, employing a modified Harten-Lax-van Leer-contact flux and a modification of the discretization for the source term. Conservation of total energy is particularly challenging in the presence of self-gravity, especially when aiming for high order accuracy. To address this, we rewrite the energy equation into a conservative form, and carefully design an energy flux with the aid of weak formulation from the DG method to maintain conservation as well as high order accuracy. The resulting scheme can be extended to high order in time discretizations. Numerical examples for two and three dimensional problems are provided to verify the desired properties of our proposed scheme, including shock-capturing, high order accuracy, well balance, and total energy conservation.
format Preprint
id arxiv_https___arxiv_org_abs_2510_04112
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle High order well-balanced and total-energy-conserving local discontinuous Galerkin methods for compressible self-gravitating Euler equations
Pan, Liang
Chen, Wei
Qiu, Jianxian
Xiong, Tao
Numerical Analysis
35L65, 65M60, 76L05
In this paper, we develop a high order structure-preserving local discontinuous Galerkin (DG) scheme for the compressible self-gravitating Euler equations, which pose great challenges due to the presence of time-dependent gravitational potential. The designed scheme is well-balanced for general polytropic equilibrium state and total energy conserving for multiple spatial dimensions without an assumption of spherical symmetry. The well-balanced property is achieved by decomposing the gravitational potential into equilibrium and perturbation parts, employing a modified Harten-Lax-van Leer-contact flux and a modification of the discretization for the source term. Conservation of total energy is particularly challenging in the presence of self-gravity, especially when aiming for high order accuracy. To address this, we rewrite the energy equation into a conservative form, and carefully design an energy flux with the aid of weak formulation from the DG method to maintain conservation as well as high order accuracy. The resulting scheme can be extended to high order in time discretizations. Numerical examples for two and three dimensional problems are provided to verify the desired properties of our proposed scheme, including shock-capturing, high order accuracy, well balance, and total energy conservation.
title High order well-balanced and total-energy-conserving local discontinuous Galerkin methods for compressible self-gravitating Euler equations
topic Numerical Analysis
35L65, 65M60, 76L05
url https://arxiv.org/abs/2510.04112