Generalization to a wider class of entropy split methods for compressible ideal MHD. (English) Zbl 07832203
Summary: The high order entropy split methods of the authors [J. Sci. Comput. 81, No. 3, 1359–1385 (2019; Zbl 1448.76110); Commun. Appl. Math. Comput. 5, No. 2, 653–678 (2023; Zbl 1524.76248)], by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [A. Harten, J. Comput. Phys. 49, 151–164 (1983; Zbl 0503.76088); M. Gerritsen and P. Olsson, J. Comput. Phys. 129, No. 2, 245–262 (1996; Zbl 0899.76281); the second author et al., J. Comput. Phys. 162, No. 1, 33–81 (2000; Zbl 0987.65086)]. Their formulation have been proven entropy conserving and stable by taking advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, dispersion relation-preserving (DRP) [C. Bogey and C. Bailly, “A family of low dispersive and low dissipative explicit schemes for computing aerodynamic noise”, in: Proceedings of the 8th AIAA/CEAS aeroacoustics conference & exhibit, 17–19 June 2002, Breckenridge, Colorado. Reston, VA: ARC (2002; doi:10.2514/6.2002-2510); E. J. Brambley, J. Comput. Phys. 324, 258–274 (2016; Zbl 1371.76101); S. Johansson, High order summation by parts operator based on a DRP scheme applied to 2D aeroacoustics 2004, Techn. Rep., Uppsala University (2004)] or Padé (compact) spatial discretizations (Hirsh, 1975) with summation-by-parts (SBP) operators (Strand, 1994). These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. The present work employs a new approach to obtain a wider class of high order entropy split methods that do not have the homogeneous property. The new method consists of a two-point numerical flux portion and a non-conservative portion in such a way that entropy conservation holds without requiring the homogeneity property of the compressible inviscid flux. For high order classical spatial central, DRP or Padé spatial discretizations, this new approach can be proven to be entropy conservative while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. More importantly, we extend this new approach to derive a new entropy split method for the equations of ideal MHD. Representative test cases are illustrated with comparison among existing methods of the same high order of accuracy.
MSC:
| 76-XX | Fluid mechanics |
Keywords:
conservative entropy split methods; high order physical preserving methods; entropy conserving methods; compressible gas dynamics; MHDCitations:
Zbl 1448.76110; Zbl 1524.76248; Zbl 0503.76088; Zbl 0899.76281; Zbl 0987.65086; Zbl 1371.76101References:
| [1] | Harten, A., On the symmetric form of systems of conservation laws with entropy. J Comput Phys, 151-164 (1983) · Zbl 0503.76088 |
| [2] | Yee, H.; Vinokur, M.; Djomehri, M., Entropy splitting and numerical dissipation. J Comput Phys, 1, 33-81 (2000) · Zbl 0987.65086 |
| [3] | Sjögreen, B.; Yee, H., High order entropy conservative central schemes for wide ranges of compressible gas dynamics and MHD flows. J Comput Phys, 153-185 (2018) · Zbl 1392.76045 |
| [4] | Sjögreen, B.; Yee, H., Entropy stable method for the Euler equations revisited: Central differencing via entropy splitting and SBP. J Sci Comput, 1359-1385 (2019) · Zbl 1448.76110 |
| [5] | Sjögreen, B.; Yee, H.; Kotov, D.; Kristsuk, A., Skew-symmetric splitting for multiscale gas dynamics and MHD turbulence flows. J Sci Comput, 43 (2020) · Zbl 1444.76073 |
| [6] | Sjögreen, B.; Yee, H., Construction of conservative numerical fluxes for the entropy split method. Comm Appl Math Comput (2021) |
| [7] | Gerritsen, M.; Olsson, P., Designing an efficient solution strategy for fluid flows. 1. A stable high order finite difference scheme and sharp shock resolution for the Euler equations. J Comput Phys, 245-262 (1996) · Zbl 0899.76281 |
| [8] | Vinokur, M.; Yee, H., Extension of efficient low dissipation high-order schemes for 3D curvilinear moving grids. Front Comput Fluid Dyn, 129D164 (2002), Also, Proceedings of the Robert MacCormack 60th Birthday Conference, June 26-28, 2000, Half Moon Bay, CA, NASA/TM-2000-209598 · Zbl 1047.76559 |
| [9] | Sjögreen, B.; Yee, H.; Vinokur, M., On high order finite-difference metric discretizations satifying GCL on moving and deforming grids. J Comput Phys, 211-220 (2014) · Zbl 1349.65333 |
| [10] | Sandham, N.; Li, Q.; Yee, H., Entropy splitting for high-order numerical simulation of compressible turbulence. J Comput Phys, 307-322 (2002) · Zbl 1139.76332 |
| [11] | Sjögreen, B.; Yee, H., On skew-symmetric splitting and entropy conservation schemes for the Euler equations, 817-827 · Zbl 1431.35117 |
| [12] | Kotov, D.; Yee, H.; Wray, A.; Hadjadj, A.; Sjögreen, B., High order numerical methods for the dynamic SGS model of turbulent flows with shocks. Commun Comput Phys, 273-300 (2016) · Zbl 1373.76061 |
| [13] | Kotov, D.; Yee, H.; Wray, A.; Sjögreen, B.; Kritsuk, A., Numerical dissipation control in high order shock-capturing schemes for LES of low speed flows. J Comput Phys, 189-202 (2016) · Zbl 1351.76037 |
| [14] | Yee, H.; Sjögreen, B., Recent developments in accuracy and stability improvement of nonlinear filter methods for DNS and LES of compressible flows. Comput Fluids, 331-348 (2018) · Zbl 1410.76128 |
| [15] | Yee, H.; Sjögreen, B., On entropy conservation and kinetic energy preservation methods. J Phys: Conf Ser (2020) |
| [16] | Sjögreen, B.; Yee, H.; Kotov, D.; Kristsuk, A., Skew-symmetric splitting for multiscale gas dynamics and MHD turbulence flows. J Sci Comput, 43 (2020) · Zbl 1444.76073 |
| [17] | Ducros, F.; Laporte, F.; Soulères, T.; Guinot, V.; Moinat, P.; Caruelle, B., High-order fluxes for conservative skew-symmetric-like schemes in structured meshes: Application to compressible flows. J Comput Phys, 114-139 (2000) · Zbl 0972.76066 |
| [18] | Kennedy, C. A.; Gruber, A., Reduced aliasing formulations of the convective terms within the Navier-Stokes equations for a compressible fluid. J Comput Phys, 3, 1676-1700 (2008) · Zbl 1290.76135 |
| [19] | Coppola, G.; Capuano, F.; Pirozzoli, S.; de Luca, L., Numerically stable formulations of convective terms for turbulent compressible flows. J Comput Phys, 86-104 (2019) · Zbl 1451.76081 |
| [20] | Johansson, S., High order summation by parts operator based on a DRP scheme applied to 2D aeroacousticsTech. rep. 2004-050 (2004), Dept. of Information Technology, Uppsala University: Dept. of Information Technology, Uppsala University Sweden, URL http://www.it.uu.se/research/publications/reports/2004-050 |
| [21] | Bogey, C.; Bailly, C., A family of low dispersive and low dissipative explicit schemes for computing aerodynamic noise, AIAA-paper 2002-2509 |
| [22] | Brambley, E., Optimized finite-difference (DRP) schemes perform poorly for decaying or growing oscillations. J Comput Phys, 258-274 (2016) · Zbl 1371.76101 |
| [23] | Tam, C. |
| [24] | Sjögreen, B.; Yee, H., Accuracy consideration by DRP schemes for DNS and LES of compressible flow computations. Comput Fluids, 123-136 (2017) · Zbl 1390.76220 |
| [25] | Tadmor, E., The numerical viscosity of entropy stable schemes for systems of conservation laws. I. Math Comp, 91-103 (1987) · Zbl 0641.65068 |
| [26] | Sjögreen, B.; Yee, H., High order compact central spatial discretization under the framework of entropy split methods, 439-454 |
| [27] | Godunov, S., Symmetric form of the equations of magnetohydrodynamics. Numer Methods Mech Continuum Medium, 1, 26-34 (1972) |
| [28] | Yee, H.; Sjögreen, B., Recent advancement of entropy split methods for compressible gas dynamics and MHD. J Appl Math Comput (2023) · Zbl 07689952 |
| [29] | FLASH user’s guide (2016), Flash Center for Computational Science, University of Chicago, URL https://usermanual.wiki/Document/FLASHmanual.1352128691/help |
| [30] | Castro, M.; Gallardo, J.; Marquina, A., Jacobian-free incomplete Riemann solvers, 295-307 · Zbl 1407.65142 |
| [31] | Colella, P.; Woodward, P. R., The Piecewise Parabolic Method (PPM) for gas-dynamical simulations. J Comput Phys, 174-201 (1984) · Zbl 0531.76082 |
| [32] | Li, S., An HLLC Riemann solver for magneto-hydrodynamics. J Comput Phys, 1, 344-357 (2005) · Zbl 1299.76302 |
| [33] | Gurski, K., An HLLC-type approximate Riemann solver for ideal magnetohydrodynamics. SIAM J Sci Comput, 2165-2187 (2004) · Zbl 1133.76358 |
| [34] | Yee, H.; Sjögreen, B., Efficient low dissipative high order schemes for multiscale MHD flows, II: Minimization of \(\nabla \cdot \mathbf{B}\) numerical error. J Sci Comput, 115-164 (2006) · Zbl 1149.76648 |
| [35] | Sjögreen, B.; Yee, H., Development of low dissipative high order filter schemes for multiscale Navier-Stokes/MHD systems. J Comput Phys, 910-934 (2007) · Zbl 1343.76053 |
| [36] | Shu, C.-W.; Osher, S., Efficient implementation of essentially non-oscillatory shock capturing schemes. J Comput Phys, 32-78 (1989) · Zbl 0674.65061 |
| [37] | Balsara, D. S.; Shu, C.-W., Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy. J Comput Phys, 405-452 (2000) · Zbl 0961.65078 |
| [38] | Sjögreen, B.; Yee, H. C., Multiresolution wavelet based adaptive numerical dissipation control for high order methods. J Sci Comput, 2, 211-255 (2004) · Zbl 1106.76411 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
