×

用高阶龙格-库塔积分器构造的积分延迟校正方法。 (英语) Zbl 1209.65073号

本文利用高阶Runge-Kutta(RK)积分器构造了一类积分延迟校正(IDC)方法来求解常微分方程初值问题。本质上,作者证明,在温和的正则性假设下,使用RK阶来求解校正回路中的误差方程,IDC的精度提高了(r)阶。然而,与谱延迟校正方法相比,在本文中,作者需要假设正交节点均匀分布,而正交节点非均匀分布的扩展将在下一篇论文中详细讨论。在本文中,只讨论了关于误差向量的离散Sobolev范数测量的不光滑性的几个关键细节。
本文的核心主要由两个定理组成。在第一篇文章中,本文提供了使用前向欧拉积分器构造的IDC方法的局部误差估计。它指出,如果一个常见初值问题的解(y(t))在连续意义上至少具有(M+2)个平滑度,那么由均匀分布的节点构造的IDC方法的局部误差为(O(h^{k+2})。在第二个主要定理中,本文提供了IDC方法的局部误差估计,这些方法在预测和校正步骤中包含了高阶RK。它指出,如果一个常见初值问题的解(y(t))在连续意义上至少具有(M+2)个平滑度,那么由均匀分布的节点构造的IDC方法的局部误差,在预测步骤中使用RK阶方法(r_0)和RK阶法(r_1,dots,r_k)在校正循环中,是\(O(h^{s+1})\),其中\(s=\sum_{j=0}^{k} rj(参考号)\leq M+1)。
这两个定理的证明遵循相似的路线(预测和校正步骤的命题和引理),但第二个更具技术性。作者还声称,使用类似的步骤证明第二个主要定理,可以在预测和校正中使用RK积分器证明延迟校正方法的类似定理。
本文最后用一些数值例子检验了一些IDC方法的精度阶数,并与理论结果相一致。为此,作者考虑了初值问题\(y'(t)=-2\pi\sin(2\pit)-2(y-\cos(2\pyt)),\四y(0)=1\)。使用问题\(y'(t)=\lambda y(t),\quad y(0)=1\)也可以考虑相同数量的函数求值的稳定性和准确性。

MSC公司:

65升06 常微分方程的多步、Runge-Kutta和外推方法
34A34飞机 非线性常微分方程和系统
65L20英寸 常微分方程数值方法的稳定性和收敛性
65升05 常微分方程初值问题的数值方法
65升70 常微分方程数值方法的误差界
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] W.Auzinger、H.Hofstätter、W.Kreuzer和E.Weirmüller,ODEs的改进缺陷校正算法。I.一般理论,数值。《算法36》(2004),第2期,135–155·Zbl 1058.65068号 ·doi:10.1023/B:NUMA.000033129.73715.7f
[2] Anne Bourlioux、Anita T.Layton和Michael L.Minion,反应流问题的高阶多隐式光谱延迟校正方法,J.Compute。物理学。189(2003),第2期,651-675·Zbl 1061.76053号 ·doi:10.1016/S0021-9991(03)00251-1
[3] Andrew Christlieb、Benjamin Ong和Jing Mei Qiu,关于嵌入积分延迟校正方法中的高阶积分器的评论,正在准备中·Zbl 1167.65389号
[4] Alok Dutt、Leslie Greengard和Vladimir Rokhlin,《常微分方程的谱延迟校正方法》,BIT 40(2000),第2期,241–266·兹比尔0959.65084 ·doi:10.1023/A:1022338906936
[5] Thomas Hagstrom和Ruhai Zhou,关于初值问题分裂方法的谱延迟校正,Commun。申请。数学。计算。科学。1 (2006), 169 – 205. ·Zbl 1105.65076号 ·doi:10.2140/camcos.2006.1.169
[6] A.Hansen和J.Strain,光谱延迟校正的收敛理论,预印本,加州大学伯克利分校,2005年2月。
[7] A.C.Hansen和J.Strain,关于延迟修正的顺序,http://www.damtp.cam。ac.uk/user/na/people/Anders/Deferred.pdf·Zbl 1217.65132号
[8] 黄敬芳,贾军,迈克尔·米尼恩,加速光谱延迟校正方法的收敛,J.Compute。物理学。214(2006),第2期,633–656·Zbl 1094.65066号 ·doi:10.1016/j.jcp.2005.10.004
[9] 黄敬芳,贾军,迈克尔·米尼恩,微分代数方程的任意阶Krylov延迟校正方法,J.Compute。物理学。221(2007),第2期,739–760·Zbl 1110.65076号 ·doi:10.1016/j.jcp.2006.06.040
[10] Anita T.Layton和Michael L.Minion,反应气体动力学的保守多隐式光谱延迟校正方法,J.Compute。物理学。194(2004),第2期,697–715·Zbl 1100.76048号 ·doi:10.1016/j.jcp.2003.09.010
[11] Anita T.Layton和Michael L.Minion,常微分方程Picard积分延迟校正方法正交节点选择的含义,BIT 45(2005),第2期,341-373·Zbl 1078.65552号 ·doi:10.1007/s10543-005-0016-1
[12] Anita T.Layton,关于半隐式Picard延迟校正方法的校正器选择,应用。数字。数学。58(2008),第6845-858号·Zbl 1143.65057号 ·doi:10.1016/j.apnum.2007.03.003
[13] Anita T.Layton和Michael L.Minion,半隐式Picard积分延迟校正方法预测因子选择的含义,Commun。申请。数学。计算。科学。2 (2007), 1 – 34. ·Zbl 1131.65059号 ·doi:10.2140/camcos.2007.2.1
[14] 刘元,舒志旺,张孟平,延迟修正时间离散化的强稳定性保持性,计算机学报。数学。26(2008),第5期,633–656·Zbl 1174.65036号
[15] Michael L.Minion,常微分方程的半隐式谱延迟校正方法,Commun。数学。科学。1(2003),第3期,471-500·Zbl 1088.65556号
[16] Michael L.Minion,基于光谱延迟修正的不可压缩流半隐式投影方法,应用。数字。数学。48(2004),第3-4、369–387号。PDE创新时间积分器研讨会·Zbl 1035.76040号 ·doi:10.1016/j.apnum.2003.11.005
[17] Robert D.Skeel,《证明延迟修正精度结果的理论框架》,SIAM J.Numer。分析。19(1982),第1期,171-196·Zbl 0489.65051号 ·doi:10.1137/0719009
[18] 夏银华,徐彦,舒志旺,局部间断Galerkin方法的有效时间离散,离散Contin。动态。系统。序列号。B 8(2007),第3期,677–693·Zbl 1141.65076号 ·doi:10.3934/dcdsb.2007.8.677
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。