Stability of the method of lines |
| |
Authors: | Satish C Reddy Lloyd N Trefethen |
| |
Institution: | (1) Department of Mathematics, Massachusetts Institute of Technology, 02139 Cambridge, MA, USA;(2) Present address: Courant Institute of Mathematical Sciences, 251 Mercer Street, 10012 NY, New York, USA;(3) Present address: Department of Computer Science, Cornell University, 14853 Ithaca, NY, USA |
| |
Abstract: | Summary It is well known that a necessary condition for the Lax-stability of the method of lines is that the eigenvalues of the spatial discretization operator, scaled by the time stepk, lie within a distanceO(k) of the stability region of the time integration formula ask 0. In this paper we show that a necessary and sufficient condition for stability, except for an algebraic factor, is that the -pseudo-eigenvalues of the same operator lie within a distanceO( )+O(k) of the stability region ask, ![epsi](/content/j6q3hr77v3627870/xxlarge949.gif) 0. Our results generalize those of an earlier paper by considering: (a) Runge-Kutta and other one-step formulas, (b) implicit as well as explicit linear multistep formulas, (c) weighted norms, (d) algebraic stability, (e) finite and infinite time intervals, and (f) stability regions with cusps.In summary, the theory presented in this paper amounts to a transplantation of the Kreiss matrix theorem from the unit disk (for simple power iterations) to an arbitrary stability region (for method of lines calculations).Work supported by an NSF Presidential Young Investigator Award to L.N. Trefethen |
| |
Keywords: | 65M20 15A18 |
本文献已被 SpringerLink 等数据库收录! |
|