Hostname: page-component-78c5997874-g7gxr Total loading time: 0 Render date: 2024-11-17T19:39:41.119Z Has data issue: false hasContentIssue false
  • English
  • Français

A New Family of High Order Unstructured MOOD and ADER Finite Volume Schemes for Multidimensional Systems of Hyperbolic Conservation Laws

Published online by Cambridge University Press:  03 June 2015

Raphaël Loubère ,
Michael Dumbser  and
Steven Diot
Raphaël Loubère*
Affiliation:
CNRS and Institut de Mathématiques de Toulouse (IMT) Université Paul-Sabatier, Toulouse, France
Michael Dumbser*
Affiliation:
Department of Civil, Environmental and Mechanical Engineering, University of Trento, Via Mesiano, 77 - 38123 Trento, Italy
Steven Diot*
Affiliation:
Fluid Dynamics and Solid Mechanics (T-3), Los Alamos National Laboratory, NM 87545 USA
*
Corresponding author.Email:raphael.loubere@math.univ-toulouse.fr
Email:michael.dumbser@unitn.it
Email:diot@lanl.gov
Get access

Abstract

In this paper, we investigate the coupling of the Multi-dimensional Optimal Order Detection (MOOD) method and the Arbitrary high order DERivatives (ADER) approach in order to design a new high order accurate, robust and computationally efficient Finite Volume (FV) scheme dedicated to solve nonlinear systems of hyperbolic conservation laws on unstructured triangular and tetrahedral meshes in two and three space dimensions, respectively. The Multi-dimensional Optimal Order Detection (MOOD) method for 2D and 3D geometries has been introduced in a recent series of papers for mixed unstructured meshes. It is an arbitrary high-order accurate Finite Volume scheme in space, using polynomial reconstructions with a posteriori detection and polynomial degree decrementing processes to deal with shock waves and other discontinuities. In the following work, the time discretization is performed with an elegant and efficient one-step ADER procedure. Doing so, we retain the good properties of the MOOD scheme, that is to say the optimal high-order of accuracy is reached on smooth solutions, while spurious oscillations near singularities are prevented. The ADER technique permits not only to reduce the cost of the overall scheme as shown on a set of numerical tests in 2D and 3D, but it also increases the stability of the overall scheme. A systematic comparison between classical unstructured ADER-WENO schemes and the new ADER-MOOD approach has been carried out for high-order schemes in space and time in terms of cost, robustness, accuracy and efficiency. The main finding of this paper is that the combination of ADER with MOOD generally outperforms the one of ADER and WENO either because at given accuracy MOOD is less expensive (memory and/or CPU time), or because it is more accurate for a given grid resolution. A large suite of classical numerical test problems has been solved on unstructured meshes for three challenging multi-dimensional systems of conservation laws: the Euler equations of compressible gas dynamics, the classical equations of ideal magneto-Hydrodynamics (MHD) and finally the relativistic MHD equations (RMHD), which constitutes a particularly challenging nonlinear system of hyperbolic partial differential equation. All tests are run on genuinely unstructured grids composed of simplex elements.

Keywords

52B1065D1868U0568U0765Y0565Y2065Z05Finite Volumehigh-orderconservation lawpolynomial reconstructionADERMOODhyperbolic PDEunstructured meshesfinite volumeone-step time discretizationlocal continuous space-time Galerkin methodWENOEuler equationsMHD equationsrelativistic MHD equations
Type
Research Article
Information
Communications in Computational Physics , Volume 16 , Issue 3 , September 2014 , pp. 718 - 763
Copyright
Copyright © Global Science Press Limited 2014

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Abgrall, R.On essentially non-oscillatory schemes on unstructured meshes: analysis and implementation. Journal of Computational Physics, 144:45–58, 1994.Google Scholar
[2]Aboiyar, T., Georgoulis, E. H., and Iske, A.Adaptive ADER Methods Using KernelBased Polyharmonic Spline WENO Reconstruction. SIAM Journal on Scientific Computing, 32:3251–3277, 2010.Google Scholar
[3]Balsara, D.Total variation diminishing scheme for relativistic magneto-hydrodynamics. The Astrophysical Journal Supplement Series, 132:83–101, 2001.Google Scholar
[4]Balsara, D. and Shu, C. W.Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy. Journal of Computational Physics, 160:405–452, 2000.Google Scholar
[5]Balsara, D. and Spicer, D.Maintaining pressure positivity in magnetohydrodynamic simulations. Journal of Computational Physics, 148:133–148, 1999.Google Scholar
[6]Balsara, D. and Spicer, D.A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations. Journal of Com-putational Physics, 149:270–292, 1999.Google Scholar
[7]Balsara, D. S.Self-adjusting, positivity preserving high order schemes for hydrodynamics and magnetohydrodynamics. Journal of Computational Physics, 231:7504–7517, 2012.Google Scholar
[8]Balsara, D. S., Meyer, C., Dumbser, M., Du, H., and Xu, Z.Efficient implementation of ADER schemes for Euler and magnetohydrodynamical flows on structured meshes - speed comparisons with RungeKutta methods. Journal of Computational Physics, 235:934–969, 2013.Google Scholar
[9]Balsara, D. S., Rumpf, T., Dumbser, M., and Munz, C. D.Efficient, high accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics. Journal of Computational Physics, 228:2480–2516, 2009.Google Scholar
[10]Barth, T. J. and Frederickson, P. O.Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction. AIAA paper no. 90–0013, 28th Aerospace Sciences Meeting January 1990.Google Scholar
[11]Barth, T. J. and Jespersen, D. C.The design and application of upwind schemes on unstructured meshes. AIAA Paper 89–0366, pp. 1–12, 1989.Google Scholar
[12]Ben-Artzi, M. and Falcovitz, J.A second-order Godunov-type scheme for compressible fluid dynamics. Journal of Computational Physics, 55:1–32, 1984.CrossRefGoogle Scholar
[13]Benson, D. J.Momentum advection on a staggered mesh. Journal of Computational Physics, 100(1):143–162, 1992.CrossRefGoogle Scholar
[14]Boris, J. P. and Book, D. L.Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works. Journal of Computational Physics, 11(1):38–69, 1973.Google Scholar
[15]Boris, J. P. and Book, D. L.Flux-corrected transport. Journal of Computational Physics, 135(2):172–186, 1997.Google Scholar
[16]Casper, J. and Atkins, H. L.A finite-volume high-order ENO scheme for two-dimensional hyperbolic systems. Journal of Computational Physics, 106:62–76, 1993.CrossRefGoogle Scholar
[17]Castro, C. C. and Toro, E. F.Solvers for the high-order riemann problem for hyperbolic balance laws. Journal of Computational Physics, 227:2481–2513, 2008.Google Scholar
[18]Clain, S., Diot, S., and Loubere, R.A high-order finite volume method for systems of conservation lawsmulti-dimensional optimal order detection (mood). Journal of Computational Physics, 230(10):4028–4050, 2011.CrossRefGoogle Scholar
[19]Cockburn, B., Hou, S., and Shu, C. W.The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: the multidimensional case. Mathematics of Computation, 54:545–581, 1990.Google Scholar
[20]Cockburn, B., Karniadakis, G. E., and Shu, C. W.Discontinuous Galerkin Methods. Lecture Notes in Computational Science and Engineering. Springer, 2000.Google Scholar
[21]Cockburn, B., Lin, S. Y., and Shu, C. W.TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one dimensional systems. Journal of Computational Physics, 84:90–113, 1989.Google Scholar
[22]Cockburn, B. and Shu, C. W.TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II: general framework. Mathematics of Computation, 52:411–435, 1989.Google Scholar
[23]Cockburn, B. and Shu, C. W.The Runge-Kutta local projection P1-Discontinuous Galerkin finite element method for scalar conservation laws. Mathematical Modelling and Numerical Analysis, 25:337–361, 1991.Google Scholar
[24]Cockburn, B. and Shu, C. W.The Runge-Kutta discontinuous Galerkin method for conservation laws V: multidimensional systems. Journal of Computational Physics, 141:199–224, 1998.Google Scholar
[25]Colella, P. and Sekora, M. D.A limiter for PPM that preserves accuracy at smooth extrema. Journal of Computational Physics, 227:7069–7076, 2008.CrossRefGoogle Scholar
[26]Dahlburg, R. B. and Picone, J. M.Evolution of the orszagtang vortex system in a compressible medium. I. initial average subsonic flow. Phys. Fluids B, 1:2153–2171, 1989.Google Scholar
[27]Dedner, A., Kemm, F., Kröner, D., Munz, C.-D., Schnitzer, T., and Wesenberg, M.Hyperbolic divergence cleaning for the MHD equations. Journal of Computational Physics, 175:645–673, 2002.Google Scholar
[28]Desveaux, V. and Berthon, C.An entropic mood scheme for the Euler equations. Submitted, 2013.Google Scholar
[29]Diot, S., Clain, S., and Loubère, R.Improved detection criteria for the multi-dimensional optimal order detection (MOOD) on unstructured meshes with very high-order polynomials. Computers and Fluids, 64:43–63, 2012.CrossRefGoogle Scholar
[30]Diot, S., Loubère, R., and Clain, S.The MOOD method in the three-dimensional case: Very-high-order finite volume method for hyperbolic systems. International Journal of Numerical Methods in Fluids, 73:362–392, 2013.CrossRefGoogle Scholar
[31]Dubiner, M.Spectral methods on triangles and other domains. Journal of Scientific Computing, 6:345–390, 1991.CrossRefGoogle Scholar
[32]Dumbser, M.Arbitrary high order PNPM schemes on unstructured meshes for the compressible Navier-Stokes equations. Computers & Fluids, 39:60–76, 2010.Google Scholar
[33]Dumbser, M., Balsara, D., Toro, E. F., and Munz, C. D.A unified framework for the construction of one-step finite-volume and discontinuous Galerkin schemes. Journal of Computational Physics, 227:8209–8253, 2008.Google Scholar
[34]Dumbser, M., Balsara, D., Toro, E. F., and Munz, C. D.A unified framework for the construction of one-step finite volume and discontinuous galerkin schemes on unstructured meshes. Journal of Computational Physics, 227:8209–8253, 2008.CrossRefGoogle Scholar
[35]Dumbser, M., Castro, M., Parés, C., and Toro, E. F.ADER schemes on unstructured meshes for non-conservative hyperbolic systems: Applications to geophysical flows. Computers and Fluids, 38:1731–1748, 2009.Google Scholar
[36]Dumbser, M., Enaux, C., and Toro, E. F.Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws. Journal of Computational Physics, 227:3971–4001, 2008.Google Scholar
[37]Dumbser, M. and Käser, M.Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems. Journal of Computational Physics, 221:693–723, 2007.Google Scholar
[38]Dumbser, M., Käser, M., Titarev, V. A., and Toro, E. F.Quadrature-free non-oscillatory finite volume schemes on unstructured meshes for nonlinear hyperbolic systems. Journal of Com-putational Physics, 226:204–243, 2007.CrossRefGoogle Scholar
[39]Dumbser, M. and Munz, C. D.ADER discontinuous Galerkin schemes for aeroacoustics. Comptes Rendus Mécanique, 333:683–687, 2005.Google Scholar
[40]Dumbser, M. and Munz, C. D.Building blocks for arbitrary high order discontinuous Galerkin schemes. Journal of Scientific Computing, 27:215–230, 2006.CrossRefGoogle Scholar
[41]Dumbser, M. and Toro, E. F.On universal Oshertype schemes for general nonlinear hyperbolic conservation laws. Communications in Computational Physics, 10:635–671, 2011.Google Scholar
[42]Dumbser, M. and Toro, E. F.A simple extension of the Osher Riemann solver to non-conservative hyperbolic systems. Journal of Scientific Computing, 48:70–88, 2011.CrossRefGoogle Scholar
[43]Dumbser, M. and Zanotti, O.Very high order PNPM schemes on unstructured meshes for the resistive relativistic MHD equations. Journal of Computational Physics, 228:6991–7006, 2009.CrossRefGoogle Scholar
[44]Dyson, R. W.Tchnique for very high order nonlinear simulation and validation. Technical Report TM-2001-210985, NASA, 2001.Google Scholar
[45]Einfeldt, B., Munz, C. D., Roe, P. L., and Sjögreen, B.On Godunov-type methods near low densities. Journal of Computational Physics, 92:273–295, 1991.Google Scholar
[46]Gassner, G., Dumbser, M., Hindenlang, F., and Munz, C. D.Explicit one-step time discretizations for discontinuous Galerkin and finite volume schemes based on local predictors. Journal of Computational Physics, 230:4232–4247, 2011.Google Scholar
[47]Giacomazzo, B. and Rezzolla, L.The exact solution of the Riemann problem in relativistic magnetohydrodynamics. Journal of Fluid Mechanics, 562:223–259, 2006.Google Scholar
[48]Gottlieb, S. and Shu, C. W.Total variation diminishing Runge-Kutta schemes. Mathematics of Computation, 67:73–85, 1998.Google Scholar
[49]Harten, A., Engquist, B., Osher, S., and Chakravarthy, S.Uniformly high order essentially non-oscillatory schemes, III. Journal of Computational Physics, 71:231–303, 1987.CrossRefGoogle Scholar
[50]Harten, A., Engquist, B., Osher, S., and Chakravarthy, S. R.Uniformly high order accurate essentially non-oscillatory schemes III. Journal of Computational Physics, 71:231–303, 1987.CrossRefGoogle Scholar
[51]Harten, A., Lax, P. D., and Leer, B.van. On upstream differencing and Godunov-type schemes for hyperbolic conservation laws. SIAM Review, 25(1):35–61, 1983.Google Scholar
[52]Harten, A. and Osher, S.Uniformly highorder accurate nonoscillatory schemes I. SIAM J. Num. Anal., 24:279–309, 1987.CrossRefGoogle Scholar
[53]Hidalgo, A. and Dumbser, M.ADER schemes for nonlinear systems of stiff advection-diffusion-reaction equations. Journal of Scientific Computing, 48:173–189, 2011.Google Scholar
[54]Honkkila, V. and Janhunen, P.HLLC solver for ideal relativistic MHD. Journal of Computational Physics, 223:643–656, 2007.CrossRefGoogle Scholar
[55]Hu, C. and Shu, C. W.Weighted essentially non-oscillatory schemes on triangular meshes. Journal of Computational Physics, 150:97–127, 1999.Google Scholar
[56]Hu, X. Y., Adams, N. A., and Shu, C. W.Positivity-preserving method for high-order conservative schemes solving compressible Euler equations. Journal of Computational Physics, 242:169–180, 2013.CrossRefGoogle Scholar
[57]Jiang, G.-S. and Shu, C. W.Efficient implementation of weighted ENO schemes. Journal of Computational Physics, 126:202–228, 1996.Google Scholar
[58]Jiang, G. S. and Wu, C. C.A high-order WENO finite difference scheme for the equations of ideal magnetohydrodynamics. Journal of Computational Physics, 150:561–594, 1999.Google Scholar
[59]Käser, M. and Iske, A.ADER schemes on adaptive triangular meshes for scalar conservation laws. Journal of Computational Physics, 205:486–508, 2005.Google Scholar
[60]Kolgan, V. P.Application of the minimum-derivative principle in the construction of finite-difference schemes for numerical analysis of discontinuous solutions in gas dynamics. Transactions of the Central Aerohydrodynamics Institute, 3(6):68–77, 1972. in Russian.Google Scholar
[61]Loubère, R. and Shashkov, M. J.A subcell remapping method on staggered polygonal grids for arbitrary-Lagrangian-Eulerian methods. J. Comput. Phys., 209:105–138, 2005.Google Scholar
[62]Luo, H., Luo, L., Nourgaliev, R., Mousseau, V. A., and Dinh, N.A reconstructed discontinuous Galerkin method for the compressible NavierStokes equations on arbitrary grids. Journal of Computational Physics, 229:6961–6978, 2010.Google Scholar
[63]Luo, H., Xia, Y., Spiegel, S., Nourgaliev, R., and Jiang, Z.A reconstructed discontinuous Galerkin method based on a Hierarchical WENO reconstruction for compressible flows on tetrahedral grids. Journal of Computational Physics, 236:477–492, 2013.Google Scholar
[64]Millington, R. C., Toro, E. F., and Nejad, L. A. M.Arbitrary High Order Methods for Conservation Laws I: The One Dimensional Scalar Case. PhD thesis, Manchester Metropolitan University, Department of Computing and Mathematics, June 1999.Google Scholar
[65]Orszag, S. A. and Tang, C. M.Small-scale structure of two-dimensional magnetohydrodynamic turbulence. Journal of Fluid Mechanics, 90:129,1979.Google Scholar
[66]Palenzuela, C., Lehner, L., Reula, O., and Rezzolla, L.Beyond ideal MHD: towards a more realistic modeling of relativistic astrophysical plasmas. Mon. Not. R. Astron. Soc., 394:1727–1740, 2009.Google Scholar
[67]Picone, J. M. and Dahlburg, R. B.Evolution of the orszagtang vortex system in a compressible medium. II. supersonic flow. Phys. Fluids B, 3:29–44, 1991.Google Scholar
[68]Qiu, J. and Shu, C. W.Finite difference WENO schemes with Lax-Wendroff type time discretization. SIAM J. Sci. Comput., 24(6):2185–2198, 2003.Google Scholar
[69]Qiu, J. and Shu, C. W.Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method II: two dimensional case. Computers and Fluids, 34:642–663, 2005.CrossRefGoogle Scholar
[70]Qiu, J. and Shu, C. W.Runge-Kutta discontinuous Galerkin method using WENO limiters. SIAM Journal on Scientific Computing, 26:907–929, 2005.Google Scholar
[71]Reisner, J., Serencsa, J., and Shkoller, S.A spacetime smooth artificial viscosity method for nonlinear conservation laws. Journal of Computational Physics, 235(0):912–933, 2013.Google Scholar
[72]Rezzolla, L. and Zanotti, O.An improved exact riemann solver for relativistic hydrodynamics. Journal of Fluid Mechanics, 449:395–411, 2001.CrossRefGoogle Scholar
[73]Rusanov, V. V.Calculation of Interaction of Non-Steady Shock Waves with Obstacles. J. Comput. Math. Phys. USSR, 1:267–279, 1961.Google Scholar
[74]Schwartzkopff, T., Dumbser, M., and Munz, C. D.Fast high order ADER schemes for linear hyperbolic equations. Journal of Computational Physics, 197:532–539, 2004.Google Scholar
[75]Schwartzkopff, T., Munz, C. D., and Toro, E. F.ADER: A high order approach for linear hyperbolic systems in 2D. Journal of Scientific Computing, 17(1–4):231–240, 2002.Google Scholar
[76]Shi, J., Hu, C., and Shu, C. W.A technique of treating negative weights in WENO schemes. Journal of Computational Physics, 175:108–127, 2002.Google Scholar
[77]Shu, C. W.Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws. NASA/CR-97-206253 ICASE Report No.97-65, November 1997.Google Scholar
[78]Shu, C. W. and Osher, S.Efficient implementation of essentially non-oscillatory shock capturing schemes. Journal of Computational Physics, 77:439–471, 1988.Google Scholar
[79]Stroud, A. H.Approximate Calculation of Multiple Integrals. Prentice-Hall Inc., Englewood Cliffs, New Jersey, 1971.Google Scholar
[80]Suresh, A. and Huynh, H. T.Accurate monotonicity-preserving schemes with Runge-Kutta time stepping. Journal of Computational Physics, 136:83–99, 1997.Google Scholar
[81]Tang, H. and Liu, T.A note on the conservative schemes for the Euler equations. Journal of Computational Physics, 218(2):451–459, 2006.Google Scholar
[82]Taube, A., Dumbser, M., Balsara, D., and Munz, C. D.Arbitrary high order discontinuous Galerkin schemes for the magnetohydrodynamic equations. Journal of Scientific Computing, 30:441–464, 2007.Google Scholar
[83]Titarev, V. A. and Toro, E. F.ADER: Arbitrary high order Godunov approach. Journal of Scientific Computing, 17(1-4):609–618, December 2002.Google Scholar
[84]Titarev, V. A. and Toro, E. F.ADER schemes for three-dimensional nonlinear hyperbolic systems. Journal of Computational Physics, 204:715–736, 2005.Google Scholar
[85]Titarev, V. A., Tsoutsanis, P., and Drikakis, D.WENO schemes for mixed-element unstructured meshes. Communications in Computational Physics, 8:585–609, 2010.Google Scholar
[86]Toro, E. F. and Titarev, V. A.Derivative Riemann solvers for systems of conservation laws and ADER methods. Journal of Computational Physics, 212(1):150–165, 2006.Google Scholar
[87]Toro, E. F.Riemann Solvers and Numerical Methods for Fluid Dynamics. Springer, second edition, 1999.Google Scholar
[88]Toro, E. F. and Hidalgo, A.ADER finite volume schemes for nonlinear reaction-diffusion equations. Applied Numerical Mathematics, 59:73–100, 2009.Google Scholar
[89]Toro, E. F., Millington, R. C., and Nejad, L. A. M.Towards very high order Godunov schemes. In Toro, E.F., editor, Godunov Methods. Theory and Applications, pages 905–938. spluwer/Plenum Academic Publishers, 2001.Google Scholar
[90]Toro, E. F. and Titarev, V. A.Solution of the generalized Riemann problem for advection-reaction equations. Proc. Roy. Soc. London, pp. 271–281, 2002.Google Scholar
[91]Toro, E. F. and Titarev, V. A.ADER schemes for scalar hyperbolic conservation laws with source terms in three space dimensions. Journal of Computational Physics, 202:196–215, 2005.Google Scholar
[92]Tsoutsanis, P., Titarev, V. A., and Drikakis, D.WENO schemes on arbitrary mixed-element unstructured meshes in three space dimensions. Journal of Computational Physics, 230:1585–1601, 2011.Google Scholar
[93]Leer, B.van. Towards the ultimate conservative difference scheme II: Monotonicity and conservation combined in a second order scheme. Journal of Computational Physics, 14:361–370, 1974.Google Scholar
[94]Leer, B. van. Towards the ultimate conservative difference scheme V: A second order sequel to Godunov’s method. Journal of Computational Physics, 32:101–136, 1979.Google Scholar
[95]Woodward, P. and Colella, P.The numerical simulation of two-dimensional fluid flow with strong shocks. Journal of Computational Physics, 54:115–173, 1984.Google Scholar
[96]Zanna, L. Del, Bucciantini, N., and Londrillo, P.An efficient shock-capturing central-type scheme for multidimensional relativistic flows II. magnetohydrodynamics. Astronomy and Astrophysics, 400:397–413, 2003.Google Scholar
[97]Zanna, L. Del, Zanotti, O., Bucciantini, N., and Londrillo, P.ECHO: an Eulerian conservative high order scheme for general relativistic magnetohydrodynamics and magnetodynamics. Astronomy and Astrophysics, 473:11–30, 2007.Google Scholar
[98]Zhang, X. and Shu, C. W.Positivity-preserving high order finite difference WENO schemes for compressible Euler equations. Journal of Computational Physics, 231:2245–2258, 2012.Google Scholar