Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
4967669 | Journal of Computational Physics | 2017 | 21 Pages |
Abstract
We describe a unique averaging procedure to design an entropy stable dissipation operator for the ideal magnetohydrodynamic (MHD) and compressible Euler equations. Often in the derivation of an entropy conservative numerical flux function much care is taken in the design and averaging of the entropy conservative numerical flux. We demonstrate in this work that if the discrete dissipation operator is not carefully chosen as well it can have deleterious effects on the numerical approximation. This is particularly true for very strong shocks or high Mach number flows present, for example, in astrophysical simulations. We present the underlying technique of how to construct a unique averaging technique for the discrete dissipation operator. We also demonstrate numerically the increased robustness of the approximation.
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Andrew R. Winters, Dominik Derigs, Gregor J. Gassner, Stefanie Walch,