刘国灿



一类刚性问题的序列分裂方法求解

刘国灿

(湖南文理学院 数学与计算科学学院, 湖南 常德, 415000)

将序列分裂方法应用于一类刚性常微分方程初值问题的求解, 通过理论分析和数例说明分裂方法处理此类刚性问题是收敛的.

刚性问题;序列分裂方法; 相容性

1 问题类描述

考虑如下类型的刚性常微分方程的初值问题:

说明 a. 问题存在刚性, 这意味着实际过程包含多个子过程, 这些子过程相互影响但变化速度相差很大. 文献[1]指出, 人们对刚性问题的认识经历了从线性到非线性的逐步深化与完善的过程. 虽然目前仍无法给出其一个严格的数学定义, 但对其本质形成了如下共识:

(a) 刚性问题欲求的解是慢变的, 但存在快速衰减的干扰[2], 这种干扰使得数值计算复杂化, 在数值稳定性、计算精度及算法实现等方面带来一系列实质困难.

(b) 存在由快速衰减干扰引起的局部坏条件[3].

(c) 刚性问题具有巨大的 Lipschitz 常数.

b. 刚性问题的解包含暂态解和稳态解. 人们感兴趣的主要是稳态解, 因此暂态解至少要衰减到可以忽略不计的时刻. 若试图在整个计算过程中使用显式方法, 整个计算的积分步数将多到在量级上和刚性比相同, 尤其是强刚性问题, 计算量大得惊人, 甚至无法实现. 由此可见刚性问题一般不宜用显式方法求解, 而必须采用稳定性好的隐式方法, 以避免或减缓数值稳定性对步长的苛刻要求.

c. 问题(1)的右端函数既包含有刚性部……