Inverse Lax-Wendroff Boundary Treatment of Discontinuous Galerkin Method for 1D Conservation Laws

Abstract

In this paper, we propose a new class of discontinuous Galerkin (DG) methods for solving 1D conservation laws on unfitted meshes. The standard DG method is used in the interior cells. For the small cut elements around the boundaries, we directly design approximation polynomials based on inverse Lax-Wendroff (ILW) principles for the inflow boundary conditions and introduce the post-processing to preserve the local conservation properties of the DG method. The theoretical analysis shows that our proposed methods have the same stability and numerical accuracy as the standard DG method in the inner region. An additional nonlinear limiter is designed to prevent spurious oscillations if a shock is near the boundary. Numerical results indicate that our methods achieve optimal numerical accuracy for smooth problems and do not introduce additional oscillations in discontinuous problems.

Keywords

Discontinuous Galerkin (DG) method / Hyperbolic conservation laws / Numerical boundary conditions / Inverse Lax-Wendroff (ILW) method / High-order accuracy / Stability analysis / 65M60 / 65M12

References

[1] Bastian P, Engwer C, Fahlke J, Ippisch O. An unfitted discontinuous Galerkin method for pore-scale simulations of solute transport. Math. Comput. Simul.. 2011, 81(10): 2051-2061

[2] Burman Erik, Hansbo P. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method. Appl. Numer. Math.. 2012, 62(4): 328-341

[3] Carpenter MH, Gottlieb D, Abarbanel S, Don W-S. The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: a study of the boundary error. SIAM J. Sci. Comput.. 1995, 16(6): 1241-1252

[4] Cheng, Z., Liu, S., Jiang, Y., Lu, J., Zhang, M., Zhang, S.: A high order boundary scheme to simulate complex moving rigid body under impingement of shock wave. Appl. Math. Mech. 42(6), 841–854 (2021)

[5] 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(190): 545-581

[6] 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(1): 90-113

[7] 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(186): 411-435

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

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

[10] Ding S, Shu C-W, Zhang M. On the conservation of finite difference WENO schemes in non-rectangular domains using the inverse Lax-Wendroff boundary treatments. J. Comput. Phys.. 2020, 415 109516

[11] Fu P, Frachon T, Kreiss G, Zahedi S. High order discontinuous cut finite element methods for linear hyperbolic conservation laws with an interface. J. Sci. Comput.. 2022, 90(3): 84

[12] Fu P, Kreiss G. High order cut discontinuous Galerkin methods for hyperbolic conservation laws in one space dimension. SIAM J. Sci. Comput.. 2021, 43(4): A2404-A2424

[13] Giuliani A. A two-dimensional stabilized discontinuous Galerkin method on curvilinear embedded boundary grids. SIAM J. Sci. Comput.. 2022, 44(1): A389-A415

[14] Gurkan C, Sticko S, Massing A. Stabilized cut discontinuous Galerkin methods for advection-reaction problems. SIAM J. Sci. Comput.. 2020, 42(5): A2620-A2654

[15] Hansbo P, Larson MG, Zahedi S. A cut finite element method for a Stokes interface problem. Appl. Numer. Math.. 2014, 85: 90-114

[16] Li T, Lu J, Shu C-W. Stability analysis of inverse Lax-Wendroff boundary treatment of high order compact difference schemes for parabolic equations. J. Comput. Appl. Math.. 2022, 400 113711

[17] Li T, Lu J, Wang P. Stability analysis of inverse Lax-Wendroff procedure for a high order compact finite difference schemes. Commun. Appl. Math. Comput.. 2024, 6(1): 142-189

[18] Li T, Shu C-W, Zhang M. Stability analysis of the inverse Lax-Wendroff boundary treatment for high order upwind-biased finite difference schemes. J. Comput. Appl. Math.. 2016, 299: 140-158

[19] Li T, Shu C-W, Zhang M. Stability analysis of the inverse Lax-Wendroff boundary treatment for high order central difference schemes for diffusion equations. J. Sci. Comput.. 2017, 70: 576-607

[20] Liu, S., Cheng, Z., Jiang, Y., Lu, J., Zhang, M., Zhang, S.: Numerical simulation of a complex moving rigid body under the impingement of a shock wave in 3D. Adv. Aerodyn. 4(1), 1–29 (2022)

[21] Liu S, Jiang Y, Shu C-W, Zhang M, Zhang S. A high order moving boundary treatment for convection-diffusion equations. J. Comput. Phys.. 2023, 473 111752

[22] Lu J, Fang J, Tan S, Shu C-W, Zhang M. Inverse Lax-Wendroff procedure for numerical boundary conditions of convection-diffusion equations. J. Comput. Phys.. 2016, 317: 276-300

[23] Lu J, Shu C-W, Tan S, Zhang M. An inverse Lax-Wendroff procedure for hyperbolic conservation laws with changing wind direction on the boundary. J. Comput. Phys.. 2021, 426 109940

[24] Müller B, Krämer-Eis S, Kummer F, Oberlack M. A high-order discontinuous Galerkin method for compressible flows with immersed boundaries. Int. J. Numer. Methods Eng.. 2017, 110(1): 3-30

[25] Qin R, Krivodonova L. A discontinuous Galerkin method for solutions of the Euler equations on Cartesian grids with embedded geometries. J. Comput. Sci.. 2013, 4(1/2): 24-35

[26] Reed, W.H., Hill, T.R.: Triangular mesh methods for the neutron transport equation. Technical report, Los Alamos Scientific Lab., N. Mex. (USA) (1973)

[27] Schoeder S, Sticko S, Kreiss G, Kronbichler M. High-order cut discontinuous Galerkin methods with local time stepping for acoustics. Int. J. Numer. Methods Eng.. 2020, 121(13): 2979-3003

[28] Shu C-W. Bertoluzza S, Falletta S, Russo G, Shu C-W. Discontinuous Galerkin methods: general approach and stability. Numerical Solutions of Partial Differential Equations. Advanced Courses in Mathematics-CRM Barcelona. 2009, Basel, Birkhäuser: 149201

[29] Song T, Main A, Scovazzi G, Ricchiuto M. The shifted boundary method for hyperbolic systems: embedded domain computations of linear waves and shallow water flows. J. Comput. Phys.. 2018, 369: 45-79

[30] Sticko S, Kreiss G. A stabilized Nitsche cut element method for the wave equation. Comput. Methods Appl. Mech. Eng.. 2016, 309: 364-387

[31] Sticko S, Kreiss G. Higher order cut finite elements for the wave equation. J. Sci. Comput.. 2019, 80: 1867-1887

[32] Tan S, Shu C-W. Inverse Lax-Wendroff procedure for numerical boundary conditions of conservation laws. J. Comput. Phys.. 2010, 229(21): 8144-8166

[33] Tan S, Shu C-W. A high order moving boundary treatment for compressible inviscid flows. J. Comput. Phys.. 2011, 230(15): 6023-6036

[34] Tan S, Shu C-W. Melnik R, Kotsireas I. Inverse Lax-Wendroff procedure for numerical boundary conditions of hyperbolic equations: survey and new developments. Advances in Applied Mathematics, Modeling and Computational Science. Fields Institute Communications. 2013, New York, Springer: 416366

[35] Vilar F, Shu C-W. Development and stability analysis of the inverse Lax-Wendroff boundary treatment for central compact schemes. ESAIM Math. Model. Numer. Anal.. 2014, 49(1): c115

[36] Zhao W, Huang J, Ruuth SJ. Boundary treatment of high order Runge-Kutta methods for hyperbolic conservation laws. J. Comput. Phys.. 2020, 421 109697

Funding

National Natural Science Foundation of China(12271499)

National Science Foundation(DMS-2309249)

Cyrus Tang Foundation(E2013-T2-456-67)

RIGHTS & PERMISSIONS

Shanghai University