Stability at Nonconforming Grid Interfaces for a High Order Discretization of the Schrödinger Equation

Anna Nissen, Gunilla Kreiss, and Margot Gerritsen

August 2011

Abstract:

In this paper we extend the results from our earlier work on stable boundary closures for the Schrödinger equation using the summation-by-parts-simultaneous approximation term (SBP-SAT) method to include stability and accuracy at non-conforming grid interfaces. Stability at the grid interface is shown by the energy method, and the estimates are generalized to multiple dimensions. The accuracy of the grid interface coupling is investigated using normal mode analysis for operators of 2nd and 4th order formal interior accuracy. We show that full accuracy is retained for the 2nd and 4th order operators. The accuracy results are extended to 6th and 8th order operators by numerical simulations, in which case two orders of accuracy is gained with respect to the lower order approximation close to the interface.

Note: Typos corrected October 7, 2011.

Available as PDF (432 kB, no cover)

Download BibTeX entry.