A novel averaging technique for discrete entropy-stable dissipation operators for ideal MHD

详细信息    查看全文

• 作者：Dominik Derigs^a ; ^{derigs@ph1.uni-koeln.de} ; Andrew R. Winters^b ; Gregor J. Gassner^b ; Stefanie Walch^a
• 关键词：Magnetohydrodynamics ; Entropy stable ; Entropy Jacobian ; Kinetic energy preserving
• 刊名：Journal of Computational Physics
• 出版年：2017
• 出版时间：1 February 2017
• 年：2017
• 卷：330
• 期：Complete
• 页码：624-632
• 全文大小：392 K
• 卷排序：330

文摘

Entropy stable schemes can be constructed with a specific choice of the numerical flux function. First, an entropy conserving flux is constructed. Secondly, an entropy stable dissipation term is added to this flux to guarantee dissipation of the discrete entropy. Present works in the field of entropy stable numerical schemes are concerned with thorough derivations of entropy conservative fluxes for ideal MHD. However, as we show in this work, if the dissipation operator is not constructed in a very specific way, it cannot lead to a generally stable numerical scheme.The two main findings presented in this paper are that the entropy conserving flux of Ismail & Roe can easily break down for certain initial conditions commonly found in astrophysical simulations, and that special care must be taken in the derivation of a discrete dissipation matrix for an entropy stable numerical scheme to be robust.We present a convenient novel averaging procedure to evaluate the entropy Jacobians of the ideal MHD and the compressible Euler equations that yields a discretization with favorable robustness properties.