Stability Analysis of Inverse Lax-Wendroff Procedure for a High order Compact Finite Difference Schemes

Abstract

This paper considers the finite difference (FD) approximations of diffusion operators and the boundary treatments for different boundary conditions. The proposed schemes have the compact form and could achieve arbitrary even order of accuracy. The main idea is to make use of the lower order compact schemes recursively, so as to obtain the high order compact schemes formally. Moreover, the schemes can be implemented efficiently by solving a series of tridiagonal systems recursively or the fast Fourier transform (FFT). With mathematical induction, the eigenvalues of the proposed differencing operators are shown to be bounded away from zero, which indicates the positive definiteness of the operators. To obtain numerical boundary conditions for the high order schemes, the simplified inverse Lax-Wendroff (SILW) procedure is adopted and the stability analysis is performed by the Godunov-Ryabenkii method and the eigenvalue spectrum visualization method. Various numerical experiments are provided to demonstrate the effectiveness and robustness of our algorithms.

Keywords

Compact scheme / Diffusion operators / Inverse Lax-Wendroff (ILW) / Fourier analysis / Eigenvalue analysis

References

[1.] Ashcroft G, Zhang X. Optimized prefactored compact schemes. J. Comput. Phys., 2003, 190: 459-477,

[2.] 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: 1241-1252,

[3.] Chu PC, Fan C. A three-point combined compact difference scheme. J. Comput. Phys., 1998, 140: 370-399,

[4.] Chu PC, Fan C. A three-point sixth-order nonuniform combined compact difference scheme. J. Comput. Phys., 1999, 148: 663-674,

[5.] Ciment M, Leventhal SH. Higher order compact implicit schemes for the wave equation. Math. Comput., 1975, 29: 985-994,

[6.] Cockburn B, Shu C-W. Nonlinearly stable compact schemes for shock calculations. SIAM J. Numer. Anal., 1994, 31: 607-627,

[7.] Deng X, Zhang H. Developing high-order weighted compact nonlinear schemes. J. Comput. Phys., 2000, 165: 22-44,

[8.] Düring B, Heuer C. High-order compact schemes for parabolic problems with mixed derivatives in multiple space dimensions. SIAM J. Numer. Anal., 2015, 53: 2113-2134,

[9.] Gustafsson, B.: High Order Difference Methods for Time Dependent PDE. Springer Series in Comput. Math. 38. Springer-Verlag, Berlin (2008)

[10.] Gustafsson, B., Kreiss, H.-O., Oliger, J.: Time-Dependent Problems and Difference Methods. Wiley-Interscience (1972)

[11.] Hixon R. Prefactored small-stencil compact schemes. J. Comput. Phys., 2000, 165: 522-541,

[12.] Ju L, Zhang J, Zhu L, Du Q. Fast explicit integration factor methods for semilinear parabolic equations. J. Sci. Comput., 2015, 62: 431-455,

[13.] Kopal, Z.: Numerical Analysis. 2nd edition. Wiley, New York (1961)

[14.] Lele SK. Compact finite difference schemes with spectral-like resolution. J. Comput. Phys., 1992, 103: 16-42,

[15.] 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,

[16.] 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,

[17.] 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,

[18.] Liao H, Sun Z, Shi H. Error estimate of fourth-order compact scheme for linear Schrödinger equations. SIAM J. Numer. Anal., 2010, 47: 4381-4401,

[19.] Liu X, Zhang S, Zhang H, Shu C-W. A new class of central compact schemes with spectral-like resolution I: Linear schemes. J. Comput. Phys., 2013, 248: 235-256,

[20.] Liu X, Zhang S, Zhang H, Shu C-W. A new class of central compact schemes with spectral-like resolution II: Hybrid weighted nonlinear schemes. J. Comput. Phys., 2015, 284: 133-154,

[21.] 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,

[22.] Mahesh K. A family of high order finite difference schemes with good spectral resolution. J. Comput. Phys., 1998, 145: 332-358,

[23.] Orszag SA, Israeli M. Numerical simulation of viscous incompressible flows. Annu. Rev. Fluid Mech., 1974, 6: 281-318,

[24.] Settle SO, Douglas CC, Kim I, Sheen D. On the derivation of highest-order compact finite difference schemes for the one- and to-dimensional Poisson equation with Dirichlet boundary conditions. SIAM J. Numer. Anal., 2013, 51: 2470-2490,

[25.] Shang JS. High-order compact-difference schemes for time-dependent Maxwell equations. J. Comput. Phys., 1999, 153: 312-333,

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

[27.] Shukla RK, Zhong X. Derivation of high-order compact finite difference schemes for non-uniform grid using polynomial interpolation. J. Comput. Phys., 2005, 204: 404-429,

[28.] Sousa E. Stability analysis of difference methods for parabolic initial value problems. J. Sci. Comput., 2006, 26: 45-66,

[29.] Strikwerda JC. Initial boundary value problems for the method of lines. J. Comput. Phys., 1980, 34: 94-107,

[30.] Sutmann G. Compact finite difference schemes of sixth order for the Helmholtz equation. J. Comput. Appl. Math., 2007, 203: 15-31,

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

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

[33.] Tan S, Wang C, Shu C-W, Ning J. Efficient implementation of high order inverse Lax-Wendroff boundary treatment for conservation laws. J. Comput. Phys., 2012, 231: 2510-2527,

[34.] 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., 2015, 49: 39-67,

[35.] Wang YM. Fourth-order compact finite difference methods and monotone iterative algorithms for quasi-linear elliptic boundary value problems. SIAM J. Numer. Anal., 2015, 53: 1032-1057,

[36.] Wang Y, Zhang J. Sixth order compact scheme combined with multigrid method and extrapolation technique for 2D Poisson equation. J. Comput. Phys., 2009, 228: 137-146,

[37.] Wang T, Guo B, Xu Q. Fourth-order compact and energy conservative difference schemes for the nonlinear Schrödinger equation in two dimensions. J. Comput. Phys., 2013, 243: 382-399,

[38.] Wang Z, Li J, Wang B, Xu Y, Chen X. A new central compact finite difference scheme with high spectral resolution for acoustic wave equation. J. Comput. Phys., 2018, 366: 191-206,

[39.] Zhang S, Jiang S, Shu C-W. Development of nonlinear weighted compact schemes with increasingly higher order accuracy. J. Comput. Phys., 2008, 227: 7294-7321,

Funding

National Natural Science Foundation of China(11901213); National Natural Science Foundation of China(11901213); National Natural Science Foundation of China(12171177); Young Elite Scientists Sponsorship Program by Henan Association for Science and Technology(2022HYTP0009); Program for Young Key Teacher of Henan Province(2021GGJS067); National Natural Science Foundation of China(118011140)