×

大型刚性非自治线性问题的最小残差多步方法。 (英语) Zbl 1458.65074号

摘要:本文的目的是介绍一种新的思路,即在求解刚性常微分方程组时如何避免大矩阵的因式分解。从显式线性多步方法的一般形式出发,我们建议在每个积分步骤上自适应地选择其系数,以最小化隐式BDF公式的残差范数。因此,我们将每个步骤上的未知数从(n)减少到(O(1),其中(n)是ODE系统的维数。我们称这种方法为最小剩余多步(MRMS)方法。在线性非自治问题的情况下,除了对ODE右侧的求值外,所得到的数值方案还需要一个线性最小二乘问题的解,每个步骤有一个薄矩阵。我们表明,该方法的阶数及其零稳定性质与使用的基本BDF公式的阶数一致。
对于隐式Euler方法的最简单模拟,研究了线性稳定性的性质。虽然经典的绝对稳定性分析与MRMS方法并不完全相关,但表明这种一步方法适用于刚性情况。在数值实验部分,我们考虑使用MRMS方法及其经典BDF对应方法对二维非自治热方程进行固定步长积分。起始值取自预设的缓慢变化的精确解。比较表明,这两种方法给出了相似的数值解,但在大系统的情况下,MRMS方法更快,并且其优势随着维数的增加而显著增加。带有实验代码的Python代码可以从GitHub存储库下载https://github.com/bfaleichik/mrms网站.

MSC公司:

65升04 刚性方程的数值方法
65升05 常微分方程初值问题的数值方法
65升06 常微分方程的多步、Runge-Kutta和外推方法
65L20英寸 常微分方程数值方法的稳定性和收敛性
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] J.Van Der Houwen,P。;P.Sommeuer,B.,一类特殊的具有扩展实稳定区间的多步Runge-Kutta方法,IMA J.Numer。分析。,2183-209年(1982年)·Zbl 0481.65038号
[2] Lebedev,V.,《求解刚性方程组的变时间步长显式差分格式》,(Vulkov,L.;Wasniewski,J.;Yalamov,P.,《数值分析及其应用》,计算机科学讲义,第1196卷(1997),施普林格-柏林-海德堡),274-283
[3] Faleichik,B。;邦达尔,I。;Byl,V.,《广义picard迭代:刚性问题的一类迭代Runge-Kutta方法》,J.Compute。申请。数学。,262、37-50(2014),NUMDIFF-13论文集·Zbl 1301.65059号
[4] Brown,P.N。;Hindmarsh,A.C.,《ODE刚性系统的无矩阵方法》,SIAM J.Numer。分析。,23, 3, 610-638 (1986) ·Zbl 0615.65078号
[5] Rosenbrock,H.H.,微分方程数值解的一些一般隐式过程,计算。J.,5,4,329-330(1963年)·兹比尔0112.07805
[6] 柯蒂斯,C。;Hirschfelder,J.O.,刚性方程的积分,Proc。国家。阿卡德。科学。,38, 3, 235-243 (1952) ·Zbl 0046.13602号
[7] 奥尔特加,J.M。;Rheinboldt,W.C.,多变量非线性方程的迭代解,(计算机科学和应用数学(1970),学术出版社),I-XV,1-572·Zbl 0949.65053号
[8] Saad,Y.,《稀疏线性系统的迭代方法》(2003),工业和应用数学学会:美国宾夕法尼亚州费城·Zbl 1002.65042号
[9] 海尔,E。;Wanner,G.,《求解常微分方程II:刚性和微分代数问题》(1996),Springer-Verlag:Springer-Verlag纽约·Zbl 0859.65067号
[10] Demmel,J.W.,《应用数值线性代数》(1997),SIAM·Zbl 0879.65017号
[11] Faleichik,B.,MRMS方法实验代码,GitHub知识库(2019)
[12] Li,X.S。;Shao,M.,部分枢轴不完全LU因式分解的超节点方法,ACM-Trans。数学。软件,37,4,43:1-43:20(2011)·兹比尔1365.65078
[13] Dahlquist,G.,常微分方程数值积分的收敛性和稳定性,数学。扫描。,4, 1, 33-53 (1956) ·Zbl 0071.11803号
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。