Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems
Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems
Journal of Computational Physics
A General Framework for Deriving Integral Preserving Numerical Methods for PDEs
SIAM Journal on Scientific Computing
Hi-index | 7.29 |
We consider numerical integration of dissipative gradient systems. For such systems, a class of special, stable integrators that strictly maintain dissipation is known, but they generally yield expensive fully implicit schemes, and when the system is large, linearization is indispensable for practical efficiency. However, this can in turn destroy the originally expected stability, and so far no effective principle has been formulated for a stable linearization. In this note, we point out that the behavior of the linearized schemes can be understood from a dynamical systems theory viewpoint and propose a simple principle for a stable linearization.