Stability of Finite-Difference Models Containing Two Boundaries or Interfaces

Lloyd N. Trefethen
1985 Mathematics of Computation  
It is known that the stability of finite-difference models of hyperbolic initial-boundary value problems is connected with the propagation and reflection of parasitic waves. Here the waves point of view is applied to models containing two boundaries or interfaces, where repeated reflection of trapped wave packets is a potential new source of instability. Our analysis accounts for various known instability phenomena in a unified way and leads to several new results, three of which are as
more » ... hich are as follows. (1) Dissipativity does not ensure stability when three or more formulas are concatenated at a boundary or internal interface. (2) Algebraic "GK.S instabilities" can be converted by a second boundary to exponential instabilities only when an infinite numerical reflection coefficient is present. (3) "GKS-stability" and "P-stability" can be established in certain problems by showing that the numerical reflection coefficient matrices have norm less than one.
doi:10.2307/2008125 fatcat:7wzbknhlqvamrod3btm56szfci