×

数值偏微分方程全隐式Runge-Kutta和间断Galerkin的快速时间解。一: 线性设置。 (英语) Zbl 1490.65178号

总结:全隐式Runge-Kutta(IRK)方法在精度和稳定性方面具有许多时间积分格式所需的特性,但由于求解阶段方程的困难,高阶IRK方法在实际中不常用于数值偏微分方程。本文介绍了一个理论和算法预处理框架,用于求解应用于线性数值偏微分方程(无代数约束)的IRK方法产生的方程组。该框架自然也适用于时间上的间断Galerkin离散。在得到稳定时间积分的空间离散化的相当一般的假设下,证明了预处理算子的条件数由一个小的一阶常数限定,与空间网格和时间步长无关,并且仅对级数/多项式阶数具有弱依赖性;例如,10阶高斯IRK的预处理算子的条件数小于2,与空间离散化和时间步长无关。新方法可以与任意现有的预条件器一起用于向后Euler型时间步长格式,并且当基础空间离散化是对称的时,可以使用三项递归Krylov方法。新方法被证明对线性抛物线和双曲线问题的各种高阶有限差分和有限元离散是有效的,证明了快速、可扩展的解决方案,精度高达10阶。新方法始终优于现有的块预处理方法,在某些情况下,使用高斯积分,新方法可以达到4阶精度,使用标准对角IRK方法所需的预处理应用程序和时钟时间约为一半。
第二部分见[B.S.Southworth公司等,SIAM J.Sci。计算。44,编号2,A636–A663(2022;Zbl 1490.65177号)].

MSC公司:

65平方米 含偏微分方程初值和初边值问题离散方程的数值解
65平方米 涉及偏微分方程的初值和初边值问题的线方法
65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65升06 常微分方程的多步、Runge-Kutta和外推方法
65F08个 迭代方法的前置条件
65层10 线性系统的迭代数值方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] G.Akrivis、C.Makridakis和R.H.Nochetto、Galerkin和Runge-Kutta方法:统一公式、后验误差估计和节点超收敛,数值。数学。,118(2011),第429-456页,https://doi.org/10.1007/s00211-011-0363-6。 ·Zbl 1228.65125号
[2] R.Anderson、J.Andrej、A.Barker、J.Bramwell、J.-S.Camier、J.Cerveny、V.Dobrev、Y.Dudouit、A.Fisher、T.Kolev、W.Pazner、M.Stowell、V.Tomov、J.Dahm、D.Medina和S.Zampini,MFEM:模块化有限元方法库,计算。数学。申请。,81(2021),第42-74页,https://doi.org/10.1016/j.camwa.2020.06.009。 ·Zbl 1524.65001号
[3] D.N.Arnold,具有间断单元的内部惩罚有限元法,SIAM J.Numer。分析。,19(1982),第742-760页,https://doi.org/10.1137/0719052。 ·Zbl 0482.65060号
[4] D.N.Arnold、F.Brezzi、B.Cockburn和L.D.Marini,椭圆问题间断Galerkin方法的统一分析,SIAM J.Numer。分析。,39(2002),第1749-1779页,https://doi.org/10.1137/S0036142901384162。 ·Zbl 1008.65080号
[5] S.Basting和E.Baánsch,任意阶非连续Galerkin时间步进方法的前提条件,ESAIM Math。模型。数字。分析。,51(2017),第1173-1195页,https://doi.org/10.1051/m2an/2016055。 ·Zbl 1375.65131号
[6] T.A.Bickart,隐式Runge-Kutta方法的有效求解过程,SIAM J.Numer。分析。,14(1977),第1022-1027页,https://doi.org/10.1137/0714069。 ·Zbl 0368.65037号
[7] W.C.Brown,交换环上的矩阵,Marcel Dekker,纽约,1993年·Zbl 0782.15001号
[8] K.Burrage,高效实现代数稳定的Runge-Kutta方法,SIAM J.Numer。分析。,19(1982),第245-258页,https://doi.org/10.1137/0719015。 ·Zbl 0483.65040号
[9] J.Butcher和D.Chen,一种新型的单隐式Runge-Kutta方法,应用。数字。数学。,34(2000),第179-188页,https://doi.org/10.1016/s0168-9274(99)00126-9. ·Zbl 0954.65058号
[10] J.C.Butcher,《关于隐式Runge-Kutta方法的实施》,BIT,16(1976),第237-240页,https://doi.org/10.1007/bf01932265。 ·Zbl 0336.65037号
[11] C.Canuto、P.Gervasio和A.Quarteroni,G-NI光谱方法的有限元预处理,SIAM J.Sci。计算。,31(2010),第4422-4451页,https://doi.org/10.1137/090746367。 ·Zbl 1223.65089号
[12] H.Chen,隐式Runge-Kutta和边值方法迭代解的分裂预条件,BIT,54(2014),第607-621页,https://doi.org/10.1007/s10543-014-0467-3。 ·Zbl 1302.65071号
[13] H.Chen,粘性波方程隐式Runge-Kutta离散的Kronecker积分裂预条件,应用。数学。型号。,40(2016),第4429-4440页,https://doi.org/10.1016/j.apm.2015.11.037。 ·Zbl 1459.65174号
[14] B.Cockburn和C.-W.Shu,对流占优问题的Runge-Kutta间断Galerkin方法,科学杂志。计算。,16(2001),第173-261页,https://doi.org/10.1023/a:1012873910884。 ·Zbl 1065.76135号
[15] R.D.Falgout和U.M.Yang,《Hypre:高性能预处理器库》,《欧洲并行处理会议记录》,2002年,第632-641页·Zbl 1056.65046号
[16] P.E.Farrell、R.C.Kirby和J.Marchena-Menendez,《Irksome:有限元方法的自动Runge-Kutta时间步进》,ACM Trans。数学。软件,47(2021),30·Zbl 1486.65074号
[17] M.J.Gander和M.Neumuler,抛物线问题的新时空并行多重网格算法分析,SIAM J.Sci。计算。,38(2016),第A2173-A2208页·Zbl 1342.65225号
[18] E.Hairer和G.Wanner,求解常微分方程II,刚性和微分代数问题,Springer Ser。计算。数学。,施普林格,柏林,1996年·Zbl 0859.65067号
[19] E.Hairer、G.Wanner和C.Lubich,《几何-数值积分:常微分方程的结构保持算法》,Springer,柏林,2002年,https://doi.org/10.1007/978-3-662-05018-7。 ·Zbl 0994.65135号
[20] W.Hoffmann和J.J.B.D.Swart,用三角矩阵逼近龙格-库塔矩阵,BIT,37(1997),第346-354页,https://doi.org/10.1007/bf02510217。 ·Zbl 0879.65053号
[21] P.J.v.d.Houwen和J.J.B.d.Swart,Runge-Kutta方法的并行线性系统求解器,高级计算。数学。,7(1997),第157-181页,https://doi.org/10.1023/a:1018990601750。 ·兹伯利0886.65078
[22] P.J.v.d.Houwen和J.J.B.d.Swart,ODE-IVP解算器的三角隐式迭代方法,SIAM J.Sci。计算。,18(1997),第41-55页,https://doi.org/10.1137/s1064827595287456。 ·Zbl 0872.65072号
[23] L.O.Jay,隐式Runge-Kutta方法的非精确简化牛顿迭代,SIAM J.Numer。分析。,38(2000),第1369-1388页,https://doi.org/10.1137/s0036142999360573。 ·Zbl 0989.65070号
[24] X.Jiao,X.Wang和Q.Chen,抛物线偏微分方程全隐式Runge-Kutta格式的最优和低记忆近最优预处理,SIAM J.Sci。计算。,43(2021),第A3527-A3551页·Zbl 07418132号
[25] C.Kennedy和M.H.Carpenter,常微分方程的对角隐式Runge-Kutta方法。《回顾》,技术报告,NASA兰利研究中心,2016年。
[26] M.Kronbichler和K.Ljungkvist,图形处理器上无矩阵高阶有限元计算的多重网格,ACM Trans。并行计算。,6(2019),第1-32页,https://doi.org/10.1145/3322813。
[27] P.Lasaint和P.Raviart,《关于求解中子输运方程的有限元方法》,《偏微分方程中有限元的数学方面》,学术出版社,纽约,1974年,第89-123页,https://doi.org/10.1016/b978-0-12-208350-1.50008-x。 ·Zbl 0341.65076号
[28] J.V.Lent和S.Vandewalle,隐式Runge-Kutta多重网格方法和抛物线偏微分方程的边值法离散,SIAM J.Sci。计算。,27(2005),第67-92页,https://doi.org/10.1137/030601144。 ·Zbl 1109.65081号
[29] R.J.LeVeque,《常微分方程和偏微分方程的有限差分方法:稳态和时间相关问题》,SIAM,费城,2007年·Zbl 1127.65080号
[30] C.Makridakis和R.H.Nochetto,进化问题高阶耗散方法的后验误差分析,Numer。数学。,104(2006),第489-514页,https://doi.org/10.1007/s00211-006-0013-6。 ·Zbl 1104.65091号
[31] T.A.Manteuffel、S.Muönzenmaier、J.Ruge和B.S.Southworth,基于非对称约简的代数多重网格,SIAM J.Sci。计算。,41(2019年),第S242-S268页,https://doi.org/10.1137/18M1193761。 ·Zbl 1436.65027号
[32] T.A.Manteuffel、J.Ruge和B.S.Southworth,基于局部近似理想约束的非对称代数多重网格,SIAM J.Sci。计算。,40(2018),第A4105-A4130页,https://doi.org/10.1137/17M1144350。 ·Zbl 1412.65130号
[33] K.A.Mardal、T.K.Nilssen和G.A.Staff,适用于抛物线偏微分方程的隐式Runge-Kutta格式的序最优预条件,SIAM J.Sci。计算。,29(2007),第361-375页,https://doi.org/10.1137/05064093x。 ·Zbl 1133.65071号
[34] J.Melenk、K.Gerdes和C.Schwab,《全离散有限元:快速求积,计算》。方法应用。机械。工程,190(2001),第4339-4364页,https://doi.org/10.1016/s0045-7825(00)00322-4. ·Zbl 0985.65141号
[35] T.K.Nilssen,G.A.Staff和K.Mardal,应用于双域方程的全隐式Runge-Kutta格式的最优预条件阶数,Numer。方法偏微分方程,27(2011),第1290-1312页,https://doi.org/10.1002/num.20582。 ·Zbl 1232.65128号
[36] S.P.Nørsett,仅具有多重实特征值的Runge Kutta方法,BIT,16(1976),第388-393页·Zbl 0345.65036号
[37] Y.Notay,《柔性共轭梯度》,SIAM J.Sci。计算。,22(2000),第1444-1460页,https://doi.org/10.1137/s1064827599362314。 ·Zbl 0980.65030号
[38] B.Orel,指数函数的实极点近似,BIT,31(1991),第144-159页,https://doi.org/10.1007/bf01952790。 ·Zbl 0737.41016号
[39] S.A.Orszag,复杂几何问题的谱方法,J.Compute。物理。,37(1980),第70-92页,https://doi.org/10.1016/0021-9991(80)90005-4. ·Zbl 0476.65078号
[40] W.Pazner,无矩阵高阶连续和间断Galerkin方法的高效低阶精细预处理,SIAM J.Sci。计算。,45(2020年),第A3055-A3083页·Zbl 1452.65051号
[41] W.Pazner和P.-O.Persson,不连续Galerkin流体模拟的分段并行全隐式Runge-Kutta解算器,J.Compute。物理。,335(2017),第700-717页,https://doi.org/10.1016/j.jcp.2017.01.050。 ·Zbl 1375.76164号
[42] M.M.Rana、V.E.Howle、K.Long、A.Meek和W.Milestone,抛物线PDE隐式Runge-Kutta方法的新块预条件,预印本,arXiv:2010.11377[math.NA],2020年·兹布尔07418116
[43] S.C.Reddy和L.N.Trefethen,直线法的稳定性,数值。数学。,62(1992),第235-267页,https://doi.org/10.1007/bf01396228。 ·Zbl 0734.65077号
[44] T.Richter、A.Springer和B.Vexler,抛物问题时间离散化的间断Galerkin方法的高效数值实现,Numer。数学。,124(2013),第151-182页,https://doi.org/10.1007/s00211-012-0511-7。 ·Zbl 1268.65138号
[45] R.R.Rosales、B.Seibold、D.Shirokoff和D.Zhou,初边值问题Runge-Kutta方法中降阶的空间表现,预印本,arXiv:1712.00897[math.NA],2017年。
[46] Y.Saad,一种灵活的内外预处理GMRES算法,SIAM J.Sci。计算。,14(1993),第461-469页·兹伯利0780.65022
[47] D.Schoátzau和C.Schwab,用间断Galerkin有限元方法的HP-版本对抛物问题进行时间离散,SIAM J.Numer。分析。,38(2000),第837-875页,https://doi.org/10.1137/s0036142999352394。 ·Zbl 0978.65091号
[48] K.Shahbazi,内部惩罚方法的惩罚参数的显式表达式,J.Compute。物理。,205(2005),第401-407页·Zbl 1072.65149号
[49] I.Smears,非连续Galerkin时间步进方法的稳健高效预处理,IMA J.Numer。分析。,37(2017),第1961-1985页,https://doi.org/10.1093/imanum/drw050。 ·Zbl 1433.65224号
[50] B.S.Southworth、O.A.Krzysik和W.Pazner,数值偏微分方程的完全隐式Runge-Kutta和间断Galerkin的快速时间解,第二部分:非线性和DAE,预印本,arXiv:2101.01776[math.NA],2021。
[51] G.A.Staff、K.-A.Mardal和T.K.Nilssen,抛物线偏微分方程全隐式Runge-Kutta格式的预处理,MIC J.,27(2006),第109-123页,https://doi.org/10.4173/mic.2006.2.3。
[52] L.N.Trefethen和M.Embree,《谱和伪谱:非正规矩阵和算子的行为》,普林斯顿大学出版社,新泽西州普林斯顿,2005年·Zbl 1085.15009号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。