High Order Structure Preserving Numerical Methods for Euler Equations with Gravitation
Dr. Yulong Xing
, Ohio State University
Hydrodynamical evolution in a gravitational field arises in many astrophysical and atmospheric problems. In this presentation, we will talk about arbitrary order structure preserving discontinuous Galerkin finite element methods for the Euler equations under gravitational fields, which can exactly capture the non-trivial steady state solutions, and at the same time maintain the non-negativity of some physical quantities. In addition, we consider the Euler–Poisson equations in spherical symmetry with an equilibrium state governed by the Lane–Emden equation, and design well-balanced and total-energy-conserving discontinuous Galerkin methods. High order semi-implicit, well-balanced asymptotic preserving finite difference scheme, for all Mach Euler equations with gravitation, may also be discussed. Extensive numerical examples, including a toy model of stellar core-collapse with a phenomenological equation of state that results in core-bounce and shock formation, are provided to verify the well-balanced property, positivity-preserving property, high-order accuracy, total energy conservation, and good resolution for both smooth and discontinuous solutions.