×

外推稳定的显式Runge-Kutta方法。 (英语) Zbl 1422.65233号

摘要:提出了求解多维非线性偏微分方程的外推稳定显式Runge-Kutta方法(ESERK)。在这种方法中,有必要计算每一步的函数(n_t)次,但稳定区域是(O(n_t^2))。因此,计算成本比传统显式算法低(O(n_t)倍。这样,可以通过使用简单的显式评估来集成刚性问题,在这种情况下,通常必须使用隐式方法。因此,它们特别适合于抛物线非线性多维偏微分方程的线(MOL)离散方法。在这项工作中,获得了沿负实轴具有扩展稳定性的一阶一阶方法。与其他传统的一阶稳定显式Runge-Kutta算法(也称为Runge-Kutta-Chebyshev码)相比,它们的稳定区域略短。随后,它们被用于通过Richardson外推导出(n_t)级二阶和四阶格式。这些四阶码的稳定域包括区间([-0.01n_t^2,0])((n_t)是总函数求值的数量),例如,它比ROCK4方法的稳定域短。然而,新算法既不受错误传播的影响(如其他Runge-Kutta-Chebyshev代码ROCK4或DUMKA),也不存在内部不稳定性。此外,许多其他类型的高阶(和低阶)方法也可以以类似的方式轻松获得。这些方法还允许在不增加额外成本的情况下调整长度步长。因此,稳定域以最优的方式精确地适应当前积分时间的问题谱,即用最少的额外阶段。我们将新技术与其他著名算法进行了比较,这些算法在非常僵化的扩散或反应扩散多维非线性方程中具有良好的结果。

MSC公司:

65平方米 含偏微分方程初值和初边值问题离散方程的数值解
65升06 常微分方程的多步、Runge-Kutta和外推方法
65L20英寸 常微分方程数值方法的稳定性和收敛性
65升04 刚性方程的数值方法
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Brown,P.N。;拜恩,G.D。;Hindmarsh,A.C.,VODE:可变效率ODE解算器,SIAM J.Sci。统计计算。,10, 5, 1038-1051 (1989) ·Zbl 0677.65075号
[2] Hindmarsh,A.C.,LSODE和LSODI,两个新的初值常微分方程解算器,SIGNUM Newsl。,15, 4, 10-11 (1980)
[3] Petzold,L.,《DASSL的描述:微分代数系统求解器》,IMACS Trans。科学。计算。,65-68 (1993)
[4] Cash,J.,使用改进的扩展后向微分公式在常微分方程中集成刚性初值问题,计算。数学。申请。,9, 5, 645-657 (1983) ·兹伯利0526.65052
[5] 海尔,E。;Wanner,G.,《求解常微分方程》。二: 刚性和微分代数问题(1996),Springer:Springer Berlin·Zbl 0859.65067号
[7] Alexander,R.,刚性常微分方程的对角隐式Runge-Kutta方法,SIAM J.Numer。分析。,1006-1021年(1977年)·Zbl 0374.65038号
[8] Asante-Asamani,E.O。;Wade,B.A.,反应扩散系统ETD方案的维数分裂,Commun。计算。物理。,19, 1343-1356 (2016) ·Zbl 1388.65045号
[9] 巴特,H。;Khaliq,A.,多维反应扩散系统的局部外推指数时间差分LOD方案,J.Compute。申请。数学。,285, 256-278 (2015) ·兹比尔1329.65205
[10] Hochbruck先生。;卢比奇,C。;Selhofer,H.,大型微分方程组的指数积分器,SIAM J.Sci。计算。,19, 5, 1552-1574 (1998) ·Zbl 0912.65058号
[11] 伊克萨鲁,L。;Berghe,G.V。;Meyer,H.D.,ODE指数拟合多步算法中的频率评估,J.Compute。申请。数学。,140, 1-2, 423-434 (2002) ·Zbl 0996.65075号
[12] Khaliq,A。;Martín-Vaquero,J。;韦德,B。;Yousuf,M.,《具有非光滑数据的反应扩散系统的平滑方案》,J.Compute。申请。数学。,223, 1, 374-386 (2009) ·Zbl 1155.65062号
[13] Martín-Vaquero,J。;Vigo-Aguiar,J.,指数拟合BDF算法:显式和隐式0-稳定方法,J.计算。申请。数学。,192, 1, 100-113 (2006) ·兹比尔1091.65074
[14] 阿舍尔,U.M。;Ruuth,S.J。;Spiteri,R.J.,含时偏微分方程的隐式显式Runge-Kutta方法,应用。数字。数学。,25, 2-3, 151-167 (1997) ·Zbl 0896.65061号
[15] 诺尔斯,L。;斯泰尼尔,L。;Ponthot,J.-P.,《板料成形数值模拟的隐式/显式组合时间积分算法》,J.Compute。申请。数学。,168, 1-2, 331-339 (2004) ·Zbl 1107.74344号
[16] Abdulle,A.,具有递推关系的四阶Chebyshev方法,SIAM J.Sci。计算。,23, 6, 2041-2054 (2001) ·Zbl 1009.65048号
[17] 阿卜杜勒。;Medovikov,A.A.,基于正交多项式的二阶切比雪夫方法,数值。数学。,90, 1, 1-18 (2001) ·Zbl 0997.65094号
[18] Hundsdorfer,W.H。;Verwer,J.G.,《与时间相关的对流扩散反应方程的数值解》(2007),Springer:Springer Berlin,Heidelberg,New York
[19] Ketcheson,D.I。;Ahmadia,A.J.,初值问题数值积分的最优稳定性多项式,Commun。申请。数学。计算。科学。,7, 2, 247-271 (2012) ·Zbl 1259.65114号
[20] Lebedev,V.I.,一种新方法,用于确定具有权重且受附加条件约束的线段上最小偏差多项式的根。第一部分,Russ.J.Numer。分析。数学。型号。,8, 3, 195-222 (1993) ·Zbl 0818.65035号
[21] Lebedev,V.I.,一种确定具有权重和附加条件的线段上最小偏差多项式根的新方法。第二部分,Russ.J.Numer。分析。数学。型号。,8, 5, 397-426 (1993) ·兹伯利0818.65036
[22] Medovikov,A.A.,抛物方程的高阶显式方法,BIT Numer。数学。,38, 2, 372-390 (1998) ·兹比尔0909.65060
[23] Sommeijer,B。;Shampine,L。;Verwer,J.,RKC:抛物线偏微分方程的显式解算器,J.Compute。申请。数学。,88, 2, 315-326 (1998) ·Zbl 0910.65067号
[24] van Der Houwen,P.J。;Sommeijer,B.P.,《关于大m值显式m阶段Runge-Kutta方法的内部稳定性》,ZAMM-J.Appl。数学。机械。,60, 10, 479-485 (1980) ·Zbl 0455.65052号
[25] Skvortsov,L.,显式稳定Runge-Kutta方法,计算。数学。数学。物理。,51, 7, 1153-1166 (2011) ·Zbl 1249.65156号
[26] 托里伦,M。;Jeltsch,R.,《基本最优显式Runge-Kutta方法及其在双曲抛物方程中的应用》,数值。数学。,106, 2, 303-334 (2007) ·Zbl 1113.65074号
[27] 列别捷夫,V。;Finogenov,S.,切比雪夫迭代法中参数排序问题的解决,苏联计算。数学。数学。物理。,13, 1, 21-41 (1974) ·Zbl 0285.65039号
[28] 海尔,E。;诺塞特,S.P。;Wanner,G.,《求解常微分方程I:非刚性问题》(1993),Springer-Verlag纽约公司:Springer-Verlag纽约有限公司,美国纽约州纽约市(第二修订版)·Zbl 0789.65048号
[29] Abdulle,A.,基于正交多项式的切比雪夫方法(2001),日内瓦大学博士论文·Zbl 0997.65094号
[30] O'Sullivan,S.,一类高阶Runge-Kutta-Chebyshev稳定性多项式,J.Compute。物理。,300, 665-678 (2015) ·Zbl 1349.65225号
[31] Alexiades,V。;阿米兹·G。;Gremaud,P.-A,抛物线问题显式格式的超时间步进加速,Commun。数字。方法工程,12,1,31-42(1996)·Zbl 0848.65061号
[32] 奥沙利文,S。;Downes,T.P.,《模拟弱电离等离子体的三维数值方法》,周一。不是。R.阿斯顿。《社会学杂志》,3761648-1658(2007)
[33] 奥沙利文,S。;O'Sullivan,C.,《关于期权定价显式有限差分方法的加速》,Quant。《金融》,11,8,1177-1191(2011)·Zbl 1266.91108号
[34] 杜阿尔特,M。;博纳文图拉,Z。;马索特,M。;Bourdon,A。;Descombes,S。;Dumont,T.,一种用于多尺度流光放电模拟的具有时空自适应和误差控制的新数值策略,J.Comput。物理。,231, 3, 1002-1019 (2012) ·Zbl 1406.76063号
[35] 杜蒙特,T。;杜阿尔特,M。;Descombes,S。;医学硕士Dronne。;Massot,M。;Louvet,V.,《真实三维几何中人类缺血性中风的模拟》,Commun。非线性科学。数字。模拟。,18, 6, 1539-1557 (2013) ·Zbl 1322.92019年
[36] Verwer,J.G。;Sommeijer,B.P.,扩散反应方程的隐式显式Runge-Kutta-Chebyshev格式,SIAM J.Sci。计算。,25, 5, 1824-1835 (2004) ·Zbl 1061.65090号
[37] Verwer,J.G。;Sommeijer,B.P。;Hundsdorfer,W.,《对流-扩散-反应问题的RKC时间步长》,J.Compute。物理。,201, 1, 61-79 (2004) ·Zbl 1059.65085号
[38] 郑,Z。;Petzold,L.,Runge-Kutta-Chebyshev投影法,J.计算。物理。,219, 2, 976-991 (2006) ·Zbl 1103.76048号
[39] 阿卜杜勒。;Cirilli,S.,S-ROCK:刚性随机微分方程的Chebyshev方法,SIAM J.Sci。计算。,30, 2, 997-1014 (2008) ·Zbl 1159.60329号
[40] 阿卜杜勒。;Li,T.,S-ROCK方法用于刚性Ito SDE,Commun。数学。科学。,6, 4, 845-868 (2008) ·Zbl 1162.60330号
[41] 杜阿尔特,M。;Descombes,S。;特诺,C。;坎德尔,S。;Massot,M.,燃烧前锋模拟的时空自适应数值方法,Combust。《火焰》,160,6,1083-1101(2013)
[42] Kleefeld,B。;Martín-Vaquero,J.,Serk2v2:刚性问题的一种新的二阶稳定显式Runge-Kutta方法,Numer。方法部分差异。Equ.、。,29, 1, 170-185 (2013) ·Zbl 1258.65066号
[43] Martín-Vaquero,J。;Janssen,B.,刚性问题的二阶稳定显式Runge-Kutta方法,计算。物理学。社区。,180, 10, 1802-1810 (2009) ·Zbl 1197.65006号
[44] Martín-Vaquero,J。;Khaliq,A。;Kleefeld,B.,多资产美式期权的稳定显式Runge-Kutta方法,计算。数学。申请。,67,61293-1308(2014)·Zbl 1386.91165号
[45] 斯托尔,J。;Buliersch,R.,用外推法对常微分方程的数值处理,Numer。数学。,8, 1-13 (1966) ·Zbl 0135.37901号
[46] Kleefeld,B。;Martín-Vaquero,J.,Serk2v3:求解轻度刚性非线性偏微分方程,J.Comput。申请。数学。,299, 194-206 (2016) ·Zbl 1333.65104号
[47] Verwer,J.G.,抛物型偏微分方程的显式Runge-Kutta方法,应用。数字。数学。,22, 1-3, 359-379 (1996) ·兹伯利0868.65064
[48] 穆索,V。;Knoll,D。;Rider,W.,基于物理的预处理和非平衡辐射扩散的Newton-Krylov方法,J.Compute。物理。,160, 2, 743-765 (2000) ·Zbl 0949.65092号
[49] Branco,J。;费雷拉,J。;de Oliveira,P.,广义Fisher-Kolmogorov-Petrovskii-Piskunov方程的数值方法,应用。数字。数学。,57, 1, 89-102 (2007) ·Zbl 1105.65123号
[50] 马伟(Ma,W.)。;Fuchssteiner,B.,Kolmogorov-Petrovskii-Piskunov方程的显式精确解,国际非线性力学杂志。,31, 3, 329-338 (1996) ·Zbl 0863.35106号
[51] 沃尔珀特,V。;Petrovskii,S.,《生物学中的反应扩散波》,《物理学》。生活评论,6,4,267-310(2009)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。