Abstract
Temporal, or "strict," stability of approximation to PDEs is much more difficult to achieve than the "classical" Lax stability. In this paper, we present a class of finite-difference schemes for hyperbolic initial boundary value problems in one and two space dimensions that possess the property of strict stability. The approximations are constructed so that all eigenvalues of corresponding differentiation matrix have a nonpositive real part. Boundary conditions are imposed by using penalty-like terms. Fourth- and sixth-order compact implicit finite-difference schemes are constructed and analyzed. Computational efficacy of the approach is corroborated by a series of numerical tests in 1-D and 2-D scalar problems.
Original language | English |
---|---|
Pages (from-to) | 42-66 |
Number of pages | 25 |
Journal | Journal of Computational Physics |
Volume | 160 |
Issue number | 1 |
DOIs | |
State | Published - 1 May 2000 |
Externally published | Yes |
Keywords
- Accuracy
- Boundary conditions
- Error bounds
- Hyperbolic PDEs
- Stability