×

低色散和耗散误差的龙格-库塔投影方法。 (英语) Zbl 1308.65115号

小结:本文提出了一种新的一步法,将龙格-库塔(RK)公式与步长后的适当投影相结合,用于初值问题的数值求解。这种投影的目的是在数值积分中保留一些第一积分。与标准正交投影相比,每个步骤的投影方向是从另一个合适的嵌入公式中获得的,因此整体方法是仿射不变量。对这些投影方法的局部误差进行了研究,结果表明,选择适当的嵌入公式可以提高谐振子的阶数。提供了Bogacki和Shampine(BS3)三阶方法的特定嵌入公式。给出了一般问题获得适当动力学方向的一些准则,以及根据这些准则保证BS3中嵌入RK方法存在的充分条件。最后,给出了一些数值实验来测试新投影方法的性能。

MSC公司:

65升05 常微分方程初值问题的数值方法
65升06 常微分方程的多步、Runge-Kutta和外推方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Bogacki,P.,Shampine,L.F.:三(2)对龙格库塔公式。申请。数学。莱特。2(4), 321-325 (1989) ·Zbl 0705.65055号 ·doi:10.1016/0893-9659(89)90079-7
[2] Brugnano,L.,Calvo,M.,Montijano,J.I.,Rández,L.:泊松系统的能量守恒方法。J.计算。申请。数学。236(16), 3890-3904 (2012) ·兹比尔1247.65092 ·doi:10.1016/j.cam.2012.02.033
[3] Brugnano,L.,Iavernaro,F.:保留保守问题所有不变量的线积分方法。J.计算。申请。数学。236(16), 3905-3919 (2012) ·Zbl 1246.65108号 ·doi:10.1016/j.cam.2012.03.026
[4] Brugnano,L.,Iavernaro,F.,Trigiante,D.:哈密顿BVM(HBVM):一类用于积分多项式哈密顿系统的“无漂移”方法。AIP确认程序。1168, 715-718 (2009) ·Zbl 1182.65188号 ·doi:10.1063/1.3241566
[5] Brugnano,L.,Iaverano,F.,Trigiante,D.:一个简单的框架,用于推导和分析ODE的有效一步方法类。申请。数学。计算。218(17), 8475-8485 (2012) ·Zbl 1245.65086号 ·doi:10.1016/j.amc.2012.01.074
[6] Calvo,M.,Franco,J.M.,Montijano,J.I.,Rández,L.:振动解初值问题的显式Runge-Kutta方法。J.公司。申请。数学。76(1-2), 195-212 (1996) ·Zbl 0871.65073号 ·doi:10.1016/S0377-0427(96)00103-3
[7] Calvo,M.,Hernández-Abreu,D.,Montijano,J.I.,Rándoz,L.:关于通过显式龙格库塔方法保存不变量。SIAM J.科学。计算。28(3), 868-885 (2006) ·Zbl 1118.65085号 ·数字对象标识码:10.1137/04061979X
[8] Calvo,M.,Laburta,M.P.,Montijano,J.I.,Rández,L.:通过显式runge-kutta方法近似保存二次一阶积分。高级计算。数学。32(3), 255-274 (2010) ·Zbl 1187.65077号 ·数字对象标识代码:10.1007/s10444-008-9105-4
[9] Calvo,M.,Laburta,M.P.,Montijano,J.I.,Rández,L.:保留lyapunov函数的投影方法。位数字。数学。50(2), 223-241 (2010) ·兹比尔1194.65090 ·doi:10.1007/s10543-010-0259-3
[10] Celledoni,E.、McLachlan,R.I.、McLaren,D.I.、Owren,B.、Quispel,G.R.W.、Wright,W.M.:节能龙格库塔方法。M2AN数学。模型。数字。分析。43(4), 645-649 (2009) ·Zbl 1169.65348号 ·doi:10.1051/m2安/2009020
[11] Cohen,D.,Hairer,E.:泊松系统的线性保能积分器。位数字。数学。51(1), 91-101 (2011) ·Zbl 1216.65175号 ·doi:10.1007/s10543-011-0310-z
[12] Cooper,G.J.:轨道问题的Runge-Kutta方法的稳定性。IMA J.数字。分析。7(1), 1-13 (1987) ·Zbl 0624.65057号 ·doi:10.1093/imanum/7.1.1
[13] Fassó,F.:刚体分裂算法的比较。J.计算。物理学。189(2), 527-538 (2003) ·Zbl 1041.70002号 ·doi:10.1016/S0021-9991(03)00232-8
[14] Grimm,V.,Quispel,G.R.W.:保持lyapunov函数的几何积分方法。钻头45(4),709-723(2005)·Zbl 1094.65070号 ·数字对象标识代码:10.1007/s10543-005-0034-z
[15] Hairer,E.:搭配方法的节能变体。J.数字。分析。Ind.申请。数学。(JNAIAM)5(1-2),73-84(2010)·Zbl 1432.65185号
[16] Hairer,E.,Lubich,C.,Wanner,G.:《几何数值积分:常微分方程的保结构算法》。柏林施普林格出版社(2002年)·Zbl 0994.65135号 ·doi:10.1007/978-3-662-05018-7
[17] van der Houwen,P.J.,Sommeijer,B.P.:计算振荡解的具有减少相位误差的显式runge kutta(Nyström)方法。SIAM J.数字。分析。24(3), 595-617 (1987) ·Zbl 0624.65058号 ·doi:10.1137/0724041
[18] van der Houwen,P.J.,Sommeijer,B.P.:隐式龙格库塔方法的相图分析。SIAM J.数字。分析。26(1), 214-229 (1989) ·Zbl 0669.65055号 ·doi:10.1137/0726012
[19] Iawernaro,F.,Pace,B.:多项式型哈密顿函数守恒的s级梯形方法。AIP确认程序。936, 603-606 (2007) ·兹比尔1152.65345 ·doi:10.1063/1.2790219
[20] Iawernaro,F.,Trigiante,D.:多项式哈密顿问题能量守恒的高阶对称格式。J.数字。分析。Ind.申请。数学。(JNAIAM)4(1-2),87-101(2009)·Zbl 1191.65169号
[21] McLachlan,R.I.,Quispel,G.R.W.:保守多项式ODE的几何积分。申请。数字。数学。45(4), 411-418 (2003) ·Zbl 1021.65064号 ·doi:10.1016/S0168-9274(03)00022-9
[22] Quispel,G.R.W.,Capel,H.W.:在数值上求解常微分方程,同时保留第一个积分。物理学。莱特。A 218(3-6),223-228(1996)·Zbl 0972.65507号 ·doi:10.1016/0375-9601(96)00403-3
[23] Quispel,G.R.W.,McLaren,D.I.:一类新的保能数值积分方法。《物理学杂志》。A: 数学。西奥。41(4), 045206,-7 (2008) ·Zbl 1132.65065号 ·doi:10.1088/1751-8113/41/4/045206
[24] Shampine,L.F.,Reichelt,M.W.:MATLAB ODE套件。SIAM J.科学。计算。18(1), 1-22 (1997) ·Zbl 0868.65040号 ·doi:10.1137/S1064827594276424
[25] Vilmart,G.:减少刚体动力学中的舍入误差。J.计算。物理学。227(15), 7083-7088 (2008) ·Zbl 1140.70002号 ·doi:10.1016/j.jcp.2008.04.013
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。