The behaviour of the local error in splitting methods applied to stiff problems

JOURNAL OF COMPUTATIONAL PHYSICS

卷:195

Article

Kozlov, R ; Kværno, A ; Owren, B

关键词: time integration;    geometric integration;    splitting methods;    numerical analysis;    order reduction;    singular perturbation;   

DOI  :  10.1016/j.jcp.2003.10.011

来源: Elsevier

【 摘 要 】

Splitting methods are frequently used in solving stiff differential equations and it is common to split the system of equations into a stiff and a nonstiff part. The classical theory for the local order of consistency is valid only for stepsizes which are smaller than what one would typically prefer to use in the integration. Error control and stepsize selection devices based on classical local order theory may lead to unstable error behaviour and inefficient stepsize sequences. Here, the behaviour of the local error in the Strang and Godunov splitting methods is explained by using two different tools, Lie series and singular perturbation theory. The two approaches provide an understanding of the phenomena from different points of view, but both are consistent with what is observed in numerical experiments. (C) 2003 Elsevier Inc. All rights reserved.

【 授权许可】

Free