RKDG Methods with Multi-resolution WENO Limiters for Solving Steady-State Problems on Triangular Meshes

Abstract

In this paper, we design high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (multi-resolution WENO) limiters to compute compressible steady-state problems on triangular meshes. A troubled cell indicator extended from structured meshes to unstructured meshes is constructed to identify triangular cells in which the application of the limiting procedures is required. In such troubled cells, the multi-resolution WENO limiting methods are used to the hierarchical

L^2

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^2$$\end{document}

projection polynomial sequence of the DG solution. Through using the RKDG methods with multi-resolution WENO limiters, the optimal high-order accuracy can be gradually reduced to first-order in the triangular troubled cells, so that the shock wave oscillations can be well suppressed. In steady-state simulations on triangular meshes, the numerical residual converges to near machine zero. The proposed spatial reconstruction methods enhance the robustness of classical DG methods on triangular meshes. The good results of these RKDG methods with multi-resolution WENO limiters are verified by a series of two-dimensional steady-state problems.

Keywords

RKDG method / Steady-state problem / Multi-resolution WENO limiter / Triangular mesh / Machine zero

References

[1.] Balsara D, Altmann C, Munz C, Dumbser M. A sub-cell based indicator for troubled zones in RKDG schemes and a novel class of hybrid RKDG+HWENO schemes. J. Comput. Phys., 2007, 226: 586-620

[2.] Barth, T.J., Jespersen, D.: The design and application of upwind schemes on unstructured meshes. In: 27th Aerospace Sciences Meeting, 09 January 1989–12 January 1989, Reno, NV, USA (1989). https://doi.org/10.2514/6.1989-366

[3.] Biswas R, Devine KD, Flaherty J. Parallel, adaptive finite element methods for conservation laws. Appl. Numer. Math., 1994, 14: 255-283

[4.] Borges R, Carmona M, Costa B, Don WS. An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. J. Comput. Phys., 2008, 227: 3191-3211

[5.] Burbeau A, Sagaut P, Bruneau CH. A problem-independent limiter for high-order Runge-Kutta discontinuous Galerkin methods. J. Comput. Phys., 2001, 169: 111-150

[6.] Capdeville G. A central WENO scheme for solving hyperbolic conservation laws on non-uniform meshes. J. Comput. Phys., 2008, 227: 2977-3014

[7.] Castro M, Costa B, Don WS. High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws. J. Comput. Phys., 2011, 230: 1766-1792

[8.] Cockburn B, Hou S, Shu C-W. The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: the multidimensional case. Math. Comput., 1990, 54: 545-581

[9.] Cockburn B, Lin S-Y, Shu C-W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: one dimensional systems. J. Comput. Phys., 1989, 84: 90-113

[10.] Cockburn B, Shu C-W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II: general framework. Math. Comput., 1989, 52: 411-435

[11.] Cockburn B, Shu C-W. The Runge-Kutta local projection P1-discontinuous Galerkin finite element method for scalar conservation laws. RAIRO Model. Math. Anal. Numer., 1991, 25: 337-361

[12.] Cockburn B, Shu C-W. The Runge-Kutta discontinuous Galerkin method for conservation laws V: multidimensional systems. J. Comput. Phys., 1998, 141: 199-224

[13.] Cockburn B, Shu C-W. Runge-Kutta discontinuous Galerkin method for convection-dominated problems. J. Sci. Comput., 2001, 16: 173-261

[14.] Cummings, R.M., Mason, W.H., Morton, S.A., McDaniel, D.R.: Applied Computational Aerodynamics. A Modern Engineering Approach. Cambridge University Press, New York (2015)

[15.] Dumbser M. Arbitrary high order P N P M \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P_N}{P_M}$$\end{document} schemes on unstructured meshes for the compressible Navier-Stokes equations. Comput. Fluids, 2010, 39: 60-76

[16.] Dumbser M, Balsara DS, Toro EF, Munz CD. A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes. J. Comput. Phys., 2008, 227: 8209-8253

[17.] Dumbser M, Zanotti O, Loubere R, Diot S. A posteriori subcell limiting of the discontinuous Galerkin finite element method for hyperbolic conservation laws. J. Comput. Phys., 2014, 278: 47-75

[18.] Friedrichs O. Weighted essentially non-oscillatory schemes for the interpolation of mean values on unstructured grids. J. Comput. Phys., 1998, 144: 194-212

[19.] Harten A. High resolution schemes for hyperbolic conservation laws. J. Comput. Phys., 1983, 49: 357-393

[20.] Hu C, Shu C-W. Weighted essentially non-oscillatory schemes on triangular meshes. J. Comput. Phys., 1999, 150: 97-127

[21.] Jameson, A.: Steady-state solutions of the Euler equations for transonic flow by a multigrid method. In: Richard E.M. (ed) Transonic, Shock, and Multidimensional Flows: Advances in Scientific Computing, pp. 37–70. Academic Press, New York (1982). https://doi.org/10.2514/6.1993-3359

[22.] Jameson, A.: Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence in transonic and hypersonic flows. In: 11th Computational Fluid Dynamics Conference, 06 July 1993–09 July 1993, Orlando, FL, USA (1993). https://doi.org/10.2514/6.1993-3359

[23.] Jameson A. A perspective on computational algorithms for aerodynamic analysis and design. Prog. Aerosp. Sci., 2001, 37: 197-243

[24.] Jameson, A., Schmidt, W., Turkel, E.: Numerical solution of the Euler equations by finite volume methods using Runge-Kutta time-stepping schemes. In: 14th Fluid and Plasma Dynamics Conference, 23 June 1981–25 June 1981, Palo Alto, CA, USA (1981). https://doi.org/10.2514/6.1981-1259

[25.] Jiang G, Shu C-W. Efficient implementation of weighted ENO schemes. J. Comput. Phys., 1996, 126: 202-228

[26.] Levy D, Puppo G, Russo G. Central WENO schemes for hyperbolic systems of conservation laws. M2AN Math. Model. Numer. Anal., 1999, 33: 547-571

[27.] Levy D, Puppo G, Russo G. Compact central WENO schemes for multidimensional conservation laws. SIAM J. Sci. Comput., 2000, 22(2): 656-672

[28.] Liu X, Osher S, Chan T. Weighted essentially non-oscillatory schemes. J. Comput. Phys., 1994, 115: 200-212

[29.] Luo H, Baum JD, Lohner R. A Hermite WENO-based limiter for discontinuous Galerkin method on unstructured grids. J. Comput. Phys., 2007, 225: 686-713

[30.] Luo H, Baum JD, Lohner R. A discontinuous Galerkin method based on a Taylor basis for the compressible flows on arbitrary grids. J. Comput. Phys., 2008, 227: 8875-8893

[31.] Luo H, Luo L, Nourgaliev R, Mousseau VA, Dinh N. A reconstructed discontinuous Galerkin method for the compressible Navier-Stokes equations on arbitrary grids. J. Comput. Phys., 2010, 229: 6961-6978

[32.] Luo H, Xia Y, Li S, Nourgaliev R, Cai C. A Hermite WENO reconstruction-based discontinuous Galerkin method for the equations on tetrahedral grids. J. Comput. Phys., 2012, 231: 5489-5503

[33.] Luo H, Xia Y, Spiegel S, Nourgaliev R, Jiang Z. A reconstructed discontinuous Galerkin method based on a hierarchical WENO reconstruction for compressible flows on tetrahedral grids. J. Comput. Phys., 2013, 236: 477-492

[34.] Osher S, Chakravarthy C. High-resolution schemes and the entropy condition. SIAM J. Numer. Anal., 1984, 21: 955-984

[35.] Qiu J, Shu C-W. A comparison of troubled-cell indicators for Runge-Kutta discontinuous Galerkin methods using weighted essentially nonoscillatory limiters. SIAM J. Sci. Comput., 2005, 27: 995-1013

[36.] Reed, W.H., Hill, T.R.: Triangular mesh methods for the neutron transport equation. Tech. Report LA-UR-73-479, Los Alamos Scientific Laboratory, Los Alamos, NM (1973)

[37.] Serna S, Marquina A. Power ENO methods: a fifth-order accurate weighted power ENO method. J. Comput. Phys., 2004, 194: 632-658

[38.] Shida Y, Kuwahara K, Ono K, Takami H. Computation of dynamic stall of a NACA-0012 airfoil. AIAA J., 1987, 25: 408-413

[39.] Shu C-W. High order weighted essentially non-oscillatory schemes for convection dominated problems. SIAM Rev., 2009, 51: 82-126

[40.] Shu C-W, Osher S. Efficient implementation of essentially non-oscillatory shock capturing schemes. J. Comput. Phys., 1988, 77: 439-471

[41.] Venkatakrishnan V. Convergence to steady state solutions of the Euler equations on unstructured grids with limiters. J. Comput. Phys., 1995, 118: 120-130

[42.] Wu L, Zhang Y-T, Zhang S, Shu C-W. High order fixed-point sweeping WENO methods for steady state of hyperbolic conservation laws and its convergence study, Commun. Comput. Phys., 2016, 20: 835-869

[43.] Yee HC, Harten A. Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates. AIAA J., 1987, 25: 266-274

[44.] Yee HC, Warming RF, Harten A. Implicit total variation diminishing (TVD) schemes for steady-state calculations. J. Comput. Phys., 1985, 57: 327-360

[45.] Zhang L, Liu W, He L, Deng X, Zhang H. A class of hybrid DG/FV methods for conservation laws I: basic formulation and one-dimensional systems. J. Comput. Phys., 2012, 231: 1081-1103

[46.] Zhang L, Liu W, He L, Deng X, Zhang H. A class of hybrid DG/FV methods for conservation laws II: two-dimensional cases. J. Comput. Phys., 2012, 231: 1104-1120

[47.] Zhang S, Jiang S, Shu C-W. Improvement of convergence to steady state solutions of Euler equations with the WENO schemes. J. Sci. Comput., 2011, 47: 216-238

[48.] Zhang S, Shu C-W. A new smoothness indicator for WENO schemes and its effect on the convergence to steady state solutions. J. Sci. Comput., 2007, 31: 273-305

[49.] Zhu J, Shu C-W. A new type of multi-resolution WENO schemes with increasingly higher-order of accuracy on triangular meshes. J. Comput. Phys., 2019, 392: 19-33

[50.] Zhu J, Shu C-W, Qiu J. High-order Runge-Kutta discontinuous Galerkin methods with a new type of multi-resolution WENO limiters on triangular meshes, Applied. Numer. Math., 2020, 153: 519-539

[51.] Zhu J, Shu C-W, Qiu J. High-order Runge-Kutta discontinuous Galerkin methods with multi-resolution WENO limiters for solving steady-state problems, Applied. Numer. Math., 2021, 165: 482-499

[52.] Zhu J, Zhong X, Shu C-W, Qiu J. Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter on unstructured meshes. Commun. Comput. Phys., 2017, 21: 623-649

Funding

National Natural Science Foundation of China(11872210); http://dx.doi.org/10.13039/501100001809(12071392); AFOSR(FA9550-20-1-0055); NSF(DMS-2010107); NSFC(MCMS-I-0120G01)