两类两步runge-two kinds of two-step runge.docxVIP

第一章引言刚性问题广泛出现于控制系统、生物学、电子网络、物理及化学等重要的科学和工程领域。二十年多来，刚性问题算法理论的研究迅速发展,现已形成一套相对成熟的理论—B-理论（详见文献[11，15，18，22，23，26～28]）。但该理论只适用于单边Lipschitz常数具有适度大小的刚性微分方程。单参数奇异摄动初值问题?x′=?f(x,y,t),?0000?εy′=g(x,y,t),t∈[t0,T],(1.1)这里常数ε满足0ε口?x(t)=x,y(t)=y,1，Jacobian矩阵gy的对数范数μ(gy)满足μ(gy)≤?c0。它是一类特殊的刚性初值问题，广泛出现在流体力学、弹性力学、量子力学、声学、光学、化学反应和最优控制等重要领域。其特征是，在讨论的微分方程中含有摄动参数ε，这种参数可以是反映一定物理特性而自然出现的，也可以是人为引进的。由于奇异摄动问题的经典Lipschitz常数和单边Lipschitz常数都具有O(ε?1)的量级[9]，而已有的B-理论主要适合于单边Lipschitz常数具有适度大小的刚性微分方程，因此奇异摄动问题不能被B-理论覆盖。鉴于这种特殊性，有必要对其数值方法的稳定性、收敛性作专门的研究。早在1988年,Hairer、Lubich和Roche[6]利用ε-渐进展开方式间接地获得了Runge-Kutta方法关于问题(1.1)的收敛性结论；1989年，他们[7]讨论了Rosenbrock方法收敛性结果；1991年Lubich[5]采用不同于ε-渐进展开方式的方法得到线性多步方法求解奇异摄动问题的整体误差估计；1993年，Schneider[30]将文献[6]中的研究推广到一般线性方法。此外，关于奇异摄动问题数值方法的收敛性研究还有单支方法、分离线性隐式Runge-Kutta方法求解奇异摄动问题的定量收敛结果。以上系列结论都是讨论单刚性奇异摄动问题，即初值问题(1.1)中的f(x,y,t)关于x的经典Lipschitz常数具有适度大小。在某些领域我们还会遇到多刚性奇异摄动问题，即f(x,y,t)关于x的经典Lipschitz常数非常巨大。肖爱国等率先研究了多刚性奇异摄动初值问题的数值方法的收敛性，直接获得了单支方法、Runge?Kutta方法及线性多步方法、分离线性隐式Runge?Kutta方法定量收敛结果（详见文献[35～37]）。这一系列的理论结果都是针对单参数的奇异摄动问题，即含单个摄动参数的奇异摄动问题。文献[40]考虑了两参数奇异摄动边值问题的渐近解，并作了简单的余项估计；肖爱国等人[38]将文献[5]中相关结论拓展到了求解两参数奇异摄动问题，并得出A(α)-稳定的线性多步法的定量收敛结果。1995年，Jackiewicz和Tracogna[44]引入了两步Runge?Kutta方法（简称TSRK方法）?mjj?yn+1=θyn?1+(1?θ)yn+h∑(vjf(Yn?1)+wjf(Yn)),?j=1?m(1.2)nin?1in∑ijn?1ijn?Y=uy+(1?u)y+h(af(Y)+bf(Y))。ijj??j=1它不同于文献[22]中提到的s级r步Runge-Kutta方法?sri∑ijnjj∑ijn?1+j?Y(n)=hc11f(t+μh,Y(n))+c12y,i=1,2,L,s,?j=1j=1?sr?y=hγf(t+μh,Y(n))+αy。??n+r∑j=1jnjj∑j=1jn?1+jn它不仅依赖于当前步的阶段值Yj(j=1,2,L,s)，还涉及前一步的阶段值jYn?1(j=1,2,L,s)。正因为这一特点，当级数s相等时两步Runge-Kutta方法较经典Runge-Kutta方法具有更高的阶且具有更好的稳定性[31]。例如文献[43]构造了具有较好绝对稳定性的半隐式的二步二级4阶Runge-Kutta方法，而对一般Runge-Kutta方法来说，唯一的4阶二级方法不是半隐式的。Albrecht[25],Butcher[16]和Hairer[10]等人分别用Albrecht′s逼近和B-级数讨论了两步Runge-Kutta方法的阶条件；随后，人们开始讨论两步Runge-Kutta方法关于常微分方程的局部和整体误差，用不同方法获得了在定步长和变步长两种情况下该方法的误差估计（详见文献[31,32,45]）；后来，Chollam和Jackiewicz[17]构造了两类具有较大绝对稳定域的两步Runge-Kutta方法；此外，文献[43,44]还讨论了该方法的A-稳定性和L-稳定性。Jackiewicz[42,45]等人将两步Runge-Kutta方法拓展到延迟微分方程。但遗憾的是，在对TSRK方法的稳定性分析上,人们着重讨论的是它的线性稳定性。相比之下，对它的非线性稳定性讨论的却很少