Inverse Lax-Wendroff Boundary Treatment for Solving Conservation Laws with Finite Volume Methods

Abstract

In this paper, we concentrate on high-order boundary treatments for finite volume methods solving hyperbolic conservation laws. The complex geometric physical domain is covered by a Cartesian mesh, resulting in the boundary intersecting the grids in various fashions. We propose two approaches to evaluate the cell averages on the ghost cells near the boundary. Both of them start from the inverse Lax-Wendroff (ILW) procedure, in which the normal spatial derivatives at inflow boundaries can be obtained by repeatedly using the governing equations and boundary conditions. After that, we can get an accurate evaluation of the ghost cell average by a Taylor expansion joined with high-order extrapolation, or by a Hermite extrapolation coupling with the cell averages on some “artificial” inner cells. The stability analysis is provided for both schemes, indicating that they can avoid the so-called “small-cell” problem. Moreover, the second method is more efficient under the premise of accuracy and stability. We perform numerical experiments on a collection of examples with the physical boundary not aligned with the grids and with various boundary conditions, indicating the high-order accuracy and efficiency of the proposed schemes.

Keywords

Inverse Lax-Wendroff (ILW) method / Numerical boundary treatment / Finite volume method / High-order accuracy / Fixed Cartesian mesh / Hyperbolic conservation laws

References

[1] BassiF, RebayS. High-order accurate discontinuous finite element solution of the 2D Euler equations. J. Comput. Phys., 1997, 138: 251-285.

[2] BergerMJ, HelzelC, LeVequeRJ. H\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$H$$\end{document}-box methods for the approximation of hyperbolic conservation laws on irregular grids. SIAM J. Numer. Anal., 2003, 41: 893-918.

[3] BoularasA , ClainS, BaudoinF. A sixth-order finite volume method for diffusion problem with curved boundaries. Appl. Math. Modell., 2017, 42: 401-422.

[4] CarpenterMH, GottliebD, AbarbanelS, DonW-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: 1241-1252.

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

[6] CostaR, ClainSL, LoubèreR, MachadoGJ. Very high-order accurate finite volume scheme on curved boundaries for the two-dimensional steady-state convection-diffusion equation with Dirichlet condition. Appl. Math. Modell., 2018, 54: 752-767.

[7] CostaR, ClainSL, MachadoGJ, NóbregaJM. Very high-order accurate finite volume scheme for the steady-state incompressible Navier-Stokes equations with polygonal meshes on arbitrary curved boundaries. Comput Methods Appl. Mech. Eng., 2022, 396. 115064

[8] CostaR, NóbregaJM, ClainSL, MachadoGJ. Efficient very high-order accurate polyhedral mesh finite volume scheme for 3D conjugate heat transfer problems in curved domains. J. Comput. Phys., 2021, 445. 110604

[9] DingS, ShuC-W, ZhangM. 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

[10] FilbetF, YangC. An inverse Lax-Wendroff method for boundary conditions applied to Boltzmann type models. J. Comput. Phys., 2013, 245: 43-61.

[11] GustafssonB, KreissH-O, SundströmA. Stability theory of difference approximations for mixed initial boundary value problem. II. Math. Comput., 1972, 26: 649-686.

[12] HelzelC, BergerMJ, LeVequeRJ. A high-resolution rotated grid method for conservation laws with embedded geometries. SIAM J. Sci. Comput., 2005, 26: 785-809.

[13] KrivodonovaL, BergerM. High-order accurate implementation of solid wall boundary conditions in curved geometries. J. Comput. Phys., 2006, 211: 492-512.

[14] LiT, LuJ, ShuC-W. Stability analysis of inverse Lax-Wendroff boundary treatment of high order compact difference schemes for parabolic equations. J. Comput. Appl. Math., 2021, 400. 113711

[15] LiT, ShuC-W, ZhangM. 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.

[16] LiT, ShuC-W, ZhangM. 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.

[17] LiuS, ChengZ, JiangY, ZhangM, ZhangS. Numerical simulation of a complex moving rigid body under the impingement of a shock wave in 3D. Adv. Aerodyn., 2022, 4: 8.

[18] LiuS, JiangY, ShuC-W, ZhangM, ZhangS. A high order moving boundary treatment for convection-diffusion equations. J. Comput. Phys., 2022, 473. 111752

[19] Liu, S., Li, T., Cheng, Z., Jiang, Y., Shu, C.-W., Zhang, M.: A new type of simplified inverse Lax-Wendroff boundary treatment I: hyperbolic conservation laws. arXiv:2402.10152 (2024)

[20] LiuX-D, OsherS, ChanT. Weighted essentially nonoscillatory schemes. J. Comput. Phys., 1994, 115: 200-212.

[21] LuJ, FangJ, TanS, ShuC-W, ZhangM. Inverse Lax-Wendroff procedure for numerical boundary conditions of convection-diffusion equations. J. Comput. Phys., 2016, 317: 276-300.

[22] LuJ, ShuC-W, TanS, ZhangM. An inverse Lax-Wendroff procedure for hyperbolic conservation laws with changing wind direction on the boundary. J. Comput. Phys., 2021, 426. 109940

[23] Mavriplis, D.: Results from the 3rd drag prediction workshop using the NSU3D unstructured mesh solver. In: 45th AIAA Aerospace Sciences Meeting and Exhibit. AIAA 2007–256 (2007)

[24] PeskinCS. Numerical analysis of blood flow in the heart. J. Comput. Phys., 1977, 25: 220-252.

[25] TanS, ShuC-W. Inverse Lax-Wendroff procedure for numerical boundary conditions of conservation laws. J. Comput. Phys., 2010, 229: 8144-8166.

[26] TanS, ShuC-W. A high order moving boundary treatment for compressible inviscid flows. J. Comput. Phys., 2011, 230: 6023-6036.

[27] TanS, WangC, ShuC-W, NingJ. Efficient implementation of high order inverse Lax-Wendroff boundary treatment for conservation laws. J. Comput. Phys., 2012, 231: 2510-2527.

[28] Vilar, F.: Development and stability analysis of the inverse Lax-Wendroff boundary treatment for central compact schemes. ESAIM Math. Modell. Numer. Anal. 49, 39–67 (2015)

[29] ZhuJ, ShuC-W. A new type of multi-resolution WENO schemes with increasingly higher order of accuracy. J. Comput. Phys., 2018, 375: 659-683.

Funding

National Natural Science Foundation of China(12271499)

R &D project of Pazhou Lab (Huangpu)(2023K0609)

RIGHTS & PERMISSIONS

Shanghai University