Abstract
We investigate the stability of a temporal discretization of interface terms in grid overlapping methods. A matrix stability analysis is performed on a model problem of the one-dimensional diffusion equation on overlapping grids. The scheme stability is first analyzed theoretically, and a proof of the unconditional stability of the first-order interface extrapolation scheme with the firstand second-order time integration for any overlap size is presented. For the higher-order schemes, we obtain explicit estimates of the spectral radius of the corresponding discrete matrix operator and document the values of the stability threshold depending on the number of grid points and the size of overlap. The influence of iterations on stability properties is also investigated. Numerical experiments are then presented relating the obtained stability bounds to the observed numerical values. Semidiscrete analysis confirms the derived scaling for the stability bounds.
Original language | English (US) |
---|---|
Pages (from-to) | 3375-3401 |
Number of pages | 27 |
Journal | SIAM Journal on Numerical Analysis |
Volume | 50 |
Issue number | 6 |
DOIs | |
State | Published - 2012 |
Externally published | Yes |
Keywords
- Backward-differentiation scheme
- Explicit interface extrapolation
- Grid overlapping methods
- Matrix analysis
- Temporal stability
ASJC Scopus subject areas
- Numerical Analysis
- Computational Mathematics
- Applied Mathematics