×

浅水方程的局部指数时间差分方法:算法和数值研究。 (英语) Zbl 1473.65141号

小结:在本文中,我们研究了指数时间差分(ETD)方法应用于旋转浅水方程的性能。与在时间上具有相同阶精度的显式时间步长相比,ETD算法在许多情况下可以减少计算时间,因为它允许使用较大的时间步长,同时仍然保持数值稳定性。为了加速ETD仿真,我们提出了一种综合ETD方法和重叠区域分解的局部方法。通过将原问题分解为许多较小的子域问题并局部求解,该方法可以加快矩阵指数向量积的计算。考虑了单层或多层浅水方程的几个标准测试案例。结果表明,局部化ETD方法在高性能计算中具有很大的潜力,因为每个子域问题都可以在每个时间步长自然地并行求解。

MSC公司:

6500万08 含偏微分方程初值和初边值问题的有限体积法
65M55型 多重网格方法;涉及偏微分方程初值和初边值问题的区域分解
76U05型 旋转流体的一般理论
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] M.Hochbruck,C.Lubich,H.Selhofer,大型微分方程组的指数积分器,SIAM科学计算杂志19(5)(1998)1552-1574·Zbl 0912.65058号
[2] S.Cox,P.Matthews,刚性系统的指数时间差分,计算物理杂志176(2)(2002)430-455·Zbl 1005.65069号
[3] 杜春华,朱文伟,指数时间差分格式及其轮廓积分修正的分析与应用,BIT数值数学45(2)(2005)307-328·Zbl 1080.65074号
[4] S.Krogstad,刚性偏微分方程的广义积分因子方法,计算物理杂志203(1)(2005)72-88·Zbl 1063.65097号
[5] M.Tokman,使用指数传播迭代(EPI)方法有效集成大型刚性ODE系统,计算物理杂志213(2)(2006)748-776·Zbl 1089.65063号
[6] M.Hochbruck,A.Ostermann,J.Schweitzer,指数Rosenbrock型方法,SIAM数值分析杂志47(1)(2009)786-803·Zbl 1193.65119号
[7] L.Ju,J.Zhang,L.Zhu,Q.Du,半线性抛物方程的快速显式积分因子方法,科学计算杂志62(2015)431-455·Zbl 1317.65172号
[8] L.Ju,J.Zhang,Q.Du,《模拟Cahn-Hilliard方程粗化动力学的快速准确算法》,计算材料科学108(2015)272-282。
[9] M.Hochbruck,A.Ostermann,指数积分器,《数值学报》19(2010)209-286·Zbl 1242.65109号
[10] C.Moler,C.Van Loan,《计算矩阵指数的十九种可疑方法》,SIAM Review 20(4)(1978)801-836·Zbl 0395.65012号
[11] C.Moler,C.Van Loan,《计算矩阵指数的十九种可疑方法》,二十五年后,《SIAM评论》45(1)(2003)3-49·Zbl 1030.65029号
[12] R.B.Sidje,Expokit:计算矩阵指数的软件包,ACM数学软件交易24(1)(1998)130-156·Zbl 0917.65063号
[13] J.Niesen,W.M.Wright,算法919:用于评估指数积分器中出现的φ-函数的Krylov子空间算法,ACM数学软件事务(TOMS)38(3)(2012)22·Zbl 1365.65185号
[14] S.Gaudreault,G.Rainwater,M.Tokman,KIOPS:指数积分器的快速自适应Krylov子空间解算器,计算物理杂志372(2018)236-255·Zbl 1418.65074号
[15] V.T.Luan,J.A.Pudykiewicz,D.R.Reynolds,《气象模型高效准确时间积分方案的进一步发展》,《计算物理杂志》376(2019)817-837·Zbl 1416.86011号
[16] E.Gallopoulos,Y.Saad,用Krylov近似方法有效求解抛物方程,SIAM科学与统计计算杂志13(5)(1992)1236-1264·Zbl 0757.65101号
[17] Y.Saad,矩阵指数算子的一些Krylov子空间近似的分析,SIAM数值分析杂志29(1)(1992)209-228·Zbl 0749.65030号
[18] M.Hochbruck,C.Lubich,关于矩阵指数算子的Krylov子空间逼近,SIAM数值分析杂志34(5)(1997)1911-1925·Zbl 0888.65032号
[19] S.Gaudreault,J.A.Pudykiewicz,一种用于球体上浅水方程数值解的有效指数时间积分方法,计算物理杂志322(2016)827-848·Zbl 1351.76106号
[20] J.Loffeld,M.Tokman,刚性问题的并行自适应Krylov指数解算器的实现,SIAM科学计算杂志36(5)(2014)C591-C616·Zbl 1307.65096号
[21] J.Zhang、C.Zhou、Y.Wang、L.Ju、Q.Du、X.Chi、D.Xu、D.Chen、Y.Liu、Z.Liu,Sunway Tahulight超级计算机上粗化动力学的极端尺度相场模拟,载《高性能计算、网络、存储和分析国际会议论文集》(SC'16),2016年,第4:1-4:12页。
[22] T.-T.-P.Hoang,L.Ju,Z.Wang,扩散问题的重叠局部指数时间差分方法,《数学科学中的通信》16(6)(2018)1531-1555·Zbl 1410.65361号
[23] X.Li,L.Ju,T.-T.-P.Hoang,基于重叠区域分解的半线性抛物方程指数时间差分方法,BIT:数值数学,待审·Zbl 1460.65115号
[24] T.-T.-P.Hoang,L.Ju,Z.Wang,扩散问题的非重叠局部指数时间差分方法,科学计算杂志82#37(2020)。https://doi.org/10.1007/s10915-020-01136-w。 ·Zbl 1431.65063号 ·doi:10.1007/s10915-020-01136-w
[25] T.Ringler,M.Petersen,R.L.Higdon,D.Jacobsen,P.W.Jones,M.Maltrud,《全球海洋建模的多分辨率方法》,《海洋建模》69(2013)211-232。
[26] J.Thuburn,T.Ringler,W.Skamarock,J.Klemp,任意结构C网格上地转模式的数值表示,计算物理杂志228(22)(2009)8321-8335·Zbl 1173.86304号
[27] T.D.Ringler,J.Thuburn,J.B.Klemp,W.C.Skamarock,任意结构C型网格能量守恒和位涡动力学的统一方法,《复合物理杂志》229(9)(2010)3065-3090·Zbl 1307.76054号
[28] Q.Du,V.Faber,M.Gunzburger,《形心voronoi细分:应用和算法》,SIAM Review 41(4)(1999)637-676·Zbl 0983.65021号
[29] Q.Du,M.D.Gunzburger,L.Ju,曲面的约束形心voronoi细分,SIAM科学计算杂志24(5)(2003)1488-1506·Zbl 1036.65101号
[30] M.Ben-Artzi,P.G.LeFloch,流形上几何相容双曲守恒律的井势理论,Annales De Linstitut Henri Poincare Non-Linear Analys-is 24(6)(2007)989-1008·Zbl 1138.35055号
[31] M.Ben-Artzi,J.Falcovitz,P.G.LeFloch,球面双曲守恒定律。几何相容有限体积格式,《计算物理杂志》228(16)(2009)5650-5668·Zbl 1168.65372号
[32] A.Beljadid,A.Mohammadian,H.M.Qiblawey,《使用四阶Adams格式的大规模浅层流动的非结构化有限体积法》,《计算机与流体》88(2013)579-589·兹比尔1390.86004
[33] T.D.Ringler、D.Jacobsen、M.Gunzburger、L.Ju、M.Duda、W.Skamarock,《探索浅水方程中的多分辨率建模方法》,《月度天气回顾》139(11)(2011)3348-3368。
[34] C.Clancy,J.A.Pudykiewicz,《关于在大气模型中使用指数时间积分方法》,Tellus A:Dynamic Meteorology and Oceanography 65(1)(2013)20898·Zbl 1349.65205号
[35] M.R.Petersen,D.W.Jacobsen,T.D.Ringler,M.W.Hecht,M.E.Maltrud,MPAS-海洋模型中任意拉格朗日-欧拉垂直坐标法的评估,海洋建模86(2015)93-113。
[36] Q.Chen,T.Ringler,P.R.Gent,《扩展潜在涡度输送涡度闭合以包括空间变异系数》,《计算机与数学应用》71(11)(2016)2206-2217·Zbl 1443.86018号
[37] K.Pieper,K.C.Sockwell,M.Gunzburger,模拟多层海洋模型的指数时间差分,计算物理杂志398(2019)108900·Zbl 1453.86019号
[38] S.Sari,T.Rowan,M.Seaid,F.Benkhaldoun,使用二维多层浅水方程模拟三维自由表面流动,通信假设物理学27(2020)1413-1442·Zbl 1473.65143号
[39] J.Thuburn,C.J.Cotter,任意多边形网格上旋转浅水方程的模拟离散化框架,SIAM科学计算杂志34(3)(2012)B203-B225·兹比尔1246.65155
[40] L.Ju,T.Ringler,M.Gunzburger,《Voronoi细分及其在气候和全球建模中的应用》,摘自:《全球大气模型的数值技术》,Springer,2010年,第313-342页。
[41] D.L.Williamson,J.B.Drake,J.J.Hack,R.Jakob,P.N.Swarztrauber,球面几何中浅水方程数值近似的标准测试集,计算物理杂志102(1)(1992)211-224·Zbl 0756.76060号
[42] M.Hochbruck,A.Ostermann,半线性抛物问题的显式指数Runge-Kutta方法,SIAM数值分析杂志43(3)(2005)1069-1090·Zbl 1093.65052号
[43] Y.Saad,Arnoldi计算大型非对称矩阵特征元方法的变化,线性代数及其应用34(1980)269-295·Zbl 0456.65017号
[44] A.Koskela,用不完全正交化方法近似平流扩散算子的矩阵指数,载于:数值数学与高级应用——ENUMATH 2013,Springer,2015,第345-353页·Zbl 1328.65107号
[45] G.Karypis,V.Kumar,用于划分非结构化图、划分网格和计算稀疏矩阵的填充-约简顺序的软件包,明尼苏达大学计算机科学与工程系,陆军HPC研究中心,明尼阿波利斯。
[46] P.J.Wolfram,T.D.Ringler,M.E.Maltrud,D.W.Jacobsen,M.R.Petersen,《通过拉格朗日(lagrangian)原位全球高性能粒子追踪(light)诊断涡流理想中纬度洋盆中的等密度扩散率》,《物理海洋学杂志》45(8)(2015)2114-2133。
[47] J.Thuburn,Y.Li,Rossby-Haurwitz波的数值模拟,Tellus A 52(2)(2000)181-189。
[48] R.K.Smith,D.G.Dritschel,《带等高线平流的Rossby-Haurwitz波试验案例再探讨》,《计算物理杂志》217(2)(2006)473-484·Zbl 1160.86305号
[49] J.-G.Li,球形多芯网格上的浅水方程,《皇家气象学会季刊》144(2018)1-12。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。