Development and stability analysis of the inverse Lax−Wendroff boundary treatment for central compact schemes^∗

Division of Applied Mathematics, Brown University, Providence, RI 02912, USA.
francoisvilar@brown.edu; shu@dam.brown.edu

Received: 10 November 2013
Revised: 26 April 2014

Abstract

In this paper, we generalize the so-called inverse Lax−Wendroff boundary treatment [S. Tan and C.-W. Shu, J. Comput. Phys. 229 (2010) 8144–8166] for the inflow boundary of a linear hyperbolic problem discretized by the recently introduced central compact schemes [X. Liu, S. Zhang, H. Zhang and C.-W. Shu, J. Comput. Phys. 248 (2013) 235–256]. The outflow boundary is treated by the classical extrapolation and a stability analysis for the resulting scheme is provided. To ensure the stability of the considered schemes provided with the chosen boundaries, the G-K-S theory [B. Gustafsson, H.-O. Kreiss and A. Sundström, Math. Comput. 26 (1972) 649–686] is used, first in the semidiscrete case then in the fully discrete case with the third-order TVD Runge−Kutta time discretization. Afterwards, due to the high algebraic complexity of the G-K-S theory, the stability is analyzed by visualizing the eigenspectrum of the discretized operators. We show in this paper that the results obtained with these two different approaches are perfectly consistent. We also illustrate the high accuracy of the presented schemes on simple test problems.