×

隐式显式多速率无穷小阶段重新启动方法。 (英语) 兹比尔1522.65113

摘要:隐式-显式(IMEX)方法是一种灵活的数值时间积分方法,用于解决被划分为刚性和非刚性过程的初值问题(IVP),与纯隐式或显式方法相比,其目标是降低计算成本。灵活IVP求解器的一种补充形式是将问题分解为快变化和慢变化动力学的多速率无穷小方法,该方法通过使用任何适当精确的算法进化“快”IVP序列来求解多速率IVP。本文介绍了一类新的高阶隐式-显式多速率方法,这些方法是为多速率IVP设计的,其中慢变化动态以IMEX方式进一步拆分。这个新类,我们称之为隐式显式多速率无穷小阶段重启(IMEX-MRI-SR),通过允许更容易地创建新的嵌入方法,改进了以前的隐式显示多速率无限小广义结构加法Runge-Kutta(IMEX-MRI-GARK)方法,并扩展了多速率指数龙格-库塔(MERK)方法,允许快速变化的动力学是非线性的,且方法是隐式的。我们利用GARK理论推导出精度高达四阶的条件,并提供了二阶和三阶精确的示例方法,这是已知的第一种具有IMEX结构的嵌入式MRI方法。然后,我们进行了数值模拟,以证明在固定步长和自适应步长设置下的收敛速度和计算性能。

MSC公司:

65升05 常微分方程初值问题的数值方法
65升06 常微分方程的多步、Runge-Kutta和外推方法
65L20英寸 常微分方程数值方法的稳定性和收敛性

软件:

罗德斯;犰狳
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 阿舍尔,U.M。;Ruuth,S.J。;Spiteri,R.J.,含时偏微分方程的隐式显式Runge-Kutta方法,应用。数字。数学。,25, 2-3, 151-167 (1997) ·Zbl 0896.65061号
[2] 库珀,G.J。;Sayfy,A.,常微分方程数值解的加法方法,数学。公司。,35, 152, 1159-1172 (1980) ·Zbl 0454.65053号
[3] 库珀,G.J。;Sayfy,A.,刚性常微分方程的加性Runge-Kutta方法,数学。公司。,40, 161, 207-218 (1983) ·Zbl 0525.65053号
[4] 肯尼迪,C.A。;Carpenter,M.H.,对流-扩散-反应方程的可加Runge-Kutta格式,应用。数字。数学。,44, 1-2, 139-181 (2003) ·Zbl 1013.65103号
[5] 肯尼迪,C.A。;Carpenter,M.H.,常微分方程的高阶可加Runge-Kutta格式,应用。数字。数学。,136, 183-205 (2019) ·Zbl 1408.65045号
[6] 三都。;Günther,M.,加法Runge-Kutta方法的广义结构方法,SIAM J.Numer。分析。,53, 1, 17-42 (2015) ·Zbl 1327.65132号
[7] Hochbruck,M。;Ostermann,A.,半线性抛物问题的显式指数Runge-Kutta方法,SIAM J.Numer。分析。,43, 3, 1069-1090 (2005) ·Zbl 1093.65052号
[8] 卢安,V.T。;奇诺莫纳,R。;Reynolds,D.R.,多速率微分方程的一类新的高阶方法,SIAM J.Sci。计算。,42、2、A1245-A1268(2020)·Zbl 1441.65061号
[9] Luan,V.T。;奇诺莫纳,R。;Reynolds,D.R.,多速率指数Rosenbrock方法(2021),arXiv:2106.05385·Zbl 1507.65119号
[10] Luan,V.T。;Ostermann,A.,抛物问题的高阶显式指数Runge-Kutta方法,J.Compute。申请。数学。,256, 168-179 (2014) ·Zbl 1314.65103号
[11] Luan,V.T。;Ostermann,A.,平行指数Rosenbrock方法,计算。数学。申请。,71, 5, 1137-1150 (2016) ·Zbl 1443.65093号
[12] Günther,M。;Sandu,A.,多速率广义可加龙格库塔方法,数值。数学。,133, 3, 497-524 (2016) ·Zbl 1344.65064号
[13] Sandu,A.,一类多速率无穷小GARK方法,SIAM J.Numer。分析。,57, 5, 2300-2327 (2019) ·Zbl 1435.65109号
[14] Sarshar,A。;罗伯茨,S。;Sandu,A.,高阶解耦多速率GARK方案的设计,SIAM J.Sci。计算。,41、2、A816-A847(2019)·Zbl 1411.65091号
[15] Wensch,J。;诺思,O。;Galant,A.,大气流动模拟的多速率无穷小步长方法,BIT Numer。数学。,49, 2, 449-473 (2009) ·Zbl 1419.76510号
[16] 奇诺莫纳,R。;Reynolds,D.R.,隐式显式多速率无穷小GARK方法,SIAM J.Sci。计算。,43、5、A3082-A3113(2021)·Zbl 1530.65076号
[17] 埃斯特普博士。;金廷,V。;罗普,D。;沙迪德,J.N。;Tarede,S.,《多尺度算子分裂的后验先验分析》,SIAM J.Numer。分析。,46, 3, 1116-1146 (2008) ·Zbl 1166.65367号
[18] Ropp,D.L。;Shadid,J.N.,含不定算子系统的算子分裂方法的稳定性:反应扩散系统,J.Compute。物理。,203, 2, 449-466 (2005) ·Zbl 1073.65088号
[19] Marchuk,G.I.,分裂方法在数学物理问题解决中的一些应用,Apl。材料,13,2,103-132(1968)·Zbl 0159.44702号
[20] Strang,G.,《关于差分格式的构造和比较》,SIAM J.Numer。分析。,5, 3, 506-517 (1968) ·Zbl 0184.38503号
[21] 斯皮特里,R.J。;Wei,S.,分步Runge-Kutta方法:表示和线性稳定性分析,J.Compute。物理。,476,第111900条pp.(2023)·Zbl 07652787号
[22] Goldman,D。;Kaper,T.J.,N阶算子分裂方案和不可逆系统,SIAM J.Numer。分析。,33, 1, 349-367 (1996) ·兹比尔0849.65070
[23] 卡尔沃,M。;德弗鲁托斯,J。;Novo,J.,对流-作用-扩散方程的线性隐式Runge-Kutta方法,应用。数字。数学。,37, 4, 535-549 (2001) ·Zbl 0983.65106号
[24] 肯尼迪,C.A。;Carpenter,M.H.,刚性常微分方程的对角隐式Runge-Kutta方法,应用。数字。数学。,146, 221-244 (2019) ·Zbl 1437.65061号
[25] 罗伯茨,S。;Loffeld,J。;Sarshar,A。;伍德沃德,C.S。;Sandu,A.,隐式多速率GARK方法,科学杂志。计算。,87, 1, 1-32 (2021) ·Zbl 1471.65084号
[26] 罗伯茨,S。;Sarshar,A。;Sandu,A.,多时间尺度刚性系统的耦合多速率无穷小GARK格式,SIAM J.Sci。计算。,42、3、A1609-A1638(2020)·Zbl 1441.65055号
[27] A.N.兰维尔。;Stewart,W.J.,《Kronecker积与随机自动机网络》,J.Compute。申请。数学。,167, 2, 429-447 (2004) ·Zbl 1104.68062号
[28] W.R.Inc.,Mathematica,13.2版,伊利诺伊州香槟市,2022年,URL。
[29] 普林斯,P。;Dormand,J.,高阶嵌入Runge-Kutta公式,J.计算。申请。数学。,7, 1, 67-75 (1981) ·Zbl 0449.65048号
[30] 桑德森,C。;Curtin,R.,Armadillo:基于模板的线性代数C++库,J.开源软件。,1, 2, 26 (2016)
[31] 桑德森,C。;Curtin,R.,《C++中一个用户友好的混合稀疏矩阵类》,(Davenport,J.H.;Kauers,M.;Labahn,G.;Urban,J.,《数学软件-ICMS 2018(2018)》,Springer International Publishing:Springer国际出版公司Cham),422-430·Zbl 1397.65343号
[32] Bogacki,P。;Shampine,L.,《Runge-Kutta公式3(2)对》,应用。数学。莱特。,2,4321-325(1989年)·Zbl 0705.65055号
[33] Hairer,大肠杆菌。;Wanner,G.,《求解常微分方程II》。刚性和微分代数问题,第14卷(1996),Springer·Zbl 0859.65067号
[34] 鱼类,交流。;Reynolds,D.R.,多速率无穷小方法的自适应时间步长控制,SIAM J.Sci。计算。,45、2、A958-A984(2023)·Zbl 1515.65210号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。