×

球面三角坐标系下浅水方程的间断Galerkin方法。 (英语) Zbl 1218.76028号

摘要:提出了一个由浅水方程控制的全球大气模型,并在非结构三角形网格上用Runge-Kutta间断伽辽金方法离散。球体上的浅水方程,一个在\(\mathbb R^3\)中的二维表面,用球面三角坐标局部表示,即三角形上的适当局部坐标映射。在每个三角形网格元素上,这导致切向动量的二维表示,因此只有两个离散动量方程。
间断Galerkin方法由一个积分公式组成,该公式需要面积(元素)积分和线(元素面)积分。这里,我们使用Rusanov数值通量来解决元素面上的不连续通量。时间离散采用了一种强稳定的三阶Runge-Kutta方法。网格中每个曲线三角形上的多项式空间(k)阶由拉格朗日基表示,需要高阶求积规则来积分元素和元素面。对于所提出的方法,除了在预处理步骤中,不需要质量矩阵反演。
考虑到Williamson等人的标准试验,对大气模型进行了验证[D.L.Williamson、J.B.Drake、J.J.Hack、R.Jakob、P.N.Swarztrauber,J.计算。物理。102,第1期,211-224(1992年;Zbl 0756.76060号)]非线性浅水方程的非定常解析解和射流初始扰动引起的正压不稳定性。在模型实验中观察到了收敛速度为O(Deltax^{k+1})。此外,还进行了数值实验,三阶时间积分方法限制了模型误差。因此,时间步长(Delta t)受到CFL条件和精度要求的限制。质量守恒与机器精度有关,能量守恒在增加网格分辨率和增加多项式阶数时收敛。

MSC公司:

76M10个 有限元方法在流体力学问题中的应用
65平方米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
76B07型 不可压缩无粘流体的自由表面势流

软件:

阿马托斯冠军
PDF格式BibTeX公司 XML格式引用
全文: 内政部 链接

参考文献:

[1] 鲍姆加德纳,J.R。;Frederickson,P.O.,两个球体的二十面体离散化,SIAM J.Numer。分析。,22, 1107-1115 (1985) ·Zbl 0601.65084号
[2] 贝伦斯,J。;拉科斯基,N。;希勒,W。;Handorf,D。;Läuter,M。;Päpke,J。;Dethloff,K.,amatos:用于大气和海洋模拟的并行自适应网格生成器,海洋模型。,10, 171-183 (2005)
[3] Ciarlet,P.G。;Raviart,P.-A,曲线单元上的插值理论,及其在有限元方法中的应用,计算。方法应用。机械。工程,1217-249(1972)·Zbl 0261.65079号
[4] Cockburn,B。;Shu,C.-W.,标量守恒律的Runge-Kutta局部投影(p^1)-间断Galerkin有限元法,数学。国防部。数字分析。,25337-361(1991年)·Zbl 0732.65094号
[5] Cockburn,B。;Shu,C.-W.,对流占优问题的Runge-Kutta间断Galerkin方法,科学杂志。计算。,16, 173-261 (2001) ·Zbl 1065.76135号
[6] Cools,R.,《Stroud以来的单项式容积规则:汇编-第2部分,J.Comput》。申请。数学。,112, 21-27 (1999) ·Zbl 0954.65021号
[7] 冷却,R。;Rabinowitz,P.,《Stroud以来的单项式容积规则:汇编》,J.Compute。申请。数学。,48, 309-326 (1993) ·Zbl 0799.65027号
[8] Cóté,J.,球面上半拉格朗日模型公制项的Lagrange乘数法,Q.J.R.Meteorol。《社会学杂志》,1141347-1352(1988)
[9] 丹尼斯·J。;Fournier,A。;斯伯茨,W.F。;St-Cyr,A。;Taylor,医学硕士。;托马斯,S.J。;Tufo,H.,光谱元素大气动力核心的高分辨率网格收敛特性和并行效率,国际高性能杂志。计算。申请。,19, 225-235 (2005)
[10] Dolejsi,V。;Feistauer,M.,无粘可压缩流动数值解的半隐式间断galerkin有限元方法,J.Compute。物理。,198, 727-746 (2004) ·Zbl 1116.76386号
[11] Fournier,A。;Taylor,医学硕士。;Tribbia,J.J.,光谱元素大气模型(SEAM):区域动力学的高分辨率并行计算和局部分辨率,蒙大拿州。我们。修订版,132726-748(2004)
[12] Galewsky,J。;斯科特·R·K。;Polvani,L.M.,测试全球浅水方程数值模型的初值问题,Tellus A,56,429-440(2004)
[13] Giraldo,F.X.,球面测地线网格上的Lagrange-Galerkin方法,J.Compute。物理。,136, 197-213 (1997) ·Zbl 0909.65066号
[14] Giraldo,F.X.,可扩展光谱元素大气模型的半隐式时间积分器,Q.J.R.Meteorol。《社会学杂志》,1312431-2454(2005)
[15] Giraldo,F.X.,旋转球体上双曲方程的基于高阶三角角的间断Galerkin方法,J.Compute。物理。,214, 447-465 (2006) ·Zbl 1089.65096号
[16] 吉拉尔多·F·X。;赫塞文,J.S。;Warburton,T.,球面浅水方程的Nodal高阶间断Galerkin方法,J.Compute。物理。,181, 499-525 (2002) ·Zbl 1178.76268号
[17] 吉拉尔多·F·X。;Warburton,T.,球上浅水方程的基于节点三角形的谱元方法,J.Compute。物理。,207, 129-150 (2005) ·Zbl 1177.86002号
[18] 吉拉尔多·F·X。;Warburton,T.,《高阶三角形间断Galerkin海洋浅水模型》,国际期刊数值。方法。流体,56899-925(2008)·Zbl 1290.86002号
[19] 哥特利布,S。;舒,C.-W。;Tadmor,E.,强稳定性保持高阶时间离散化方法,SIAM Rev.,43,89-112(2001)·Zbl 0967.65098号
[20] Heikes,R。;Randall,D.A.,扭曲二十面体网格上浅水方程的数值积分。第一部分:基本设计和试验结果,周一。我们。第123版,1862-1880(1995)
[21] 海因策,T。;Hense,A.,球面上的浅水方程及其Lagrange-Galerkin解,Meteorol。大气。物理。,81, 129-137 (2002)
[22] Hesthaven,J.S.,《从静电学到单纯形多项式插值的几乎最优节点集》,SIAM J.Numer。分析。,35, 655-676 (1998) ·Zbl 0933.41004号
[23] Krivodonova,L。;Xin,J。;Remacle,J.-F。;Chevaugeon,N。;Flaherty,J.E.,双曲守恒律的冲击检测和间断Galerkin方法限制,应用。数字。数学。,48, 323-338 (2004) ·Zbl 1038.65096号
[24] Läuter,M。;Handorf,D。;Dethloff,K.,球面浅水方程的非定常分析解,J.Compute。物理。,210, 535-553 (2005) ·Zbl 1078.86001号
[25] Läuter,M。;Handorf,D。;拉科斯基,N。;贝伦斯,J。;弗里肯豪斯,S。;贝斯特,M。;德特洛夫,K。;Hiller,W.,《大气的并行自适应正压模式》,J.Compute。物理。,223, 609-628 (2007) ·Zbl 1162.35437号
[26] LeSaint,P。;Raviart,P.A.,《关于求解中子输运方程的有限元方法》,(de Boor,C.,《偏微分方程中有限元的数学方面》(1974),学术出版社),89-145·Zbl 0341.65076号
[27] 林,S.-J。;Rood,R.B.,球体上的显式通量形式半拉格朗日浅水模型,Q.J.R.Meteorol。Soc.,123,2477-2498(1997)
[28] Lyness,J。;Cools,R.,《三角形上数值容积的调查》,Proc。交响乐团。申请。数学。,48, 127-150 (1994) ·Zbl 0820.41026号
[29] 穆勒,德特列夫,旋转球体上浅水的角动量预算,物理学。版本A,45,8,5545-5555(1992)
[30] 奈尔·R·D。;托马斯,S.J。;Loft,R.D.,《间断Galerkin全球浅水模型》,周一。我们。修订版,133876-888(2005)
[31] Rancic,M。;帕斯尔,R.J。;Mesinger,F.,《使用扩展球面立方体的全球浅水模型:gnomonic与共形坐标》,Q.J.R.Meteorol。Soc.,122,959-982(1996)
[32] M.Restelli,F.X.Giraldo,非静力中尺度模拟中Navier-Stokes方程的保守间断Galerkin sem-implicit公式,SIAM J.Sci。计算。,正在审查中。;M.Restelli,F.X.Giraldo,非静力中尺度模拟中Navier-Stokes方程的保守间断Galerkin sem-implicit公式,SIAM J.Sci。计算。,审查中·兹比尔1405.65127
[33] Ronchi,C。;拉科诺(Iacono,R.)。;Paolucci,P.S.,《立方球面:球面几何中偏微分方程的一种新的求解方法》,J.Compute。物理。,124, 93-114 (1996) ·Zbl 0849.76049号
[34] Rossmanith,J.A.,球面上双曲线系统的波传播方法,J.Compute。物理。,213, 629-658 (2006) ·Zbl 1089.65088号
[35] Rossmanith,J.A。;贝尔,D.S。;LeVeque,R.J.,曲面流形上双曲系统的波传播算法,J.Compute。物理。,199, 631-662 (2004) ·Zbl 1126.76350号
[36] Ruuth,S.J。;Spiteri,R.J.,关于强稳定保持时间离散化方法的两个障碍,科学杂志。计算。,17, 211-220 (2002) ·Zbl 1003.65107号
[37] Sadourny,R.,拟均匀球面网格上原始方程的保守有限差分近似,Mon。我们。修订版,100136-144(1972)
[38] 萨多尔尼,R。;荒川,A。;Mintz,Y.,非发散正压涡度方程与球面二十面体六角网格的积分,Mon。我们。修订版,6351(1968)
[39] 舒,C.-W。;Osher,S.,《本质上非振荡冲击捕获方案的有效实现》,J.Compute。物理。,77, 439-471 (1988) ·Zbl 0653.65072号
[40] A.St-Cyr,C.Jablonowski,J.M.Dennis,H.M.Tufo,S.J.Thomas,《两种浅水模型与非一致自适应网格的比较》,Mon。我们。版次,印刷中。;A.St-Cyr,C.Jablonowski,J.M.Dennis,H.M.Tufo,S.J.Thomas,《两种浅水模型与非一致自适应网格的比较》,Mon。我们。版次,正在印刷中。
[41] Stroud,A.H.,《多重积分的近似计算》(1971),Prentice-Hall:Prentice-Hall London·Zbl 0379.65013号
[42] Stuhne,G.R。;Peltier,W.R.,球面上浅水方程的新二十面体网格点离散,J.Compute。物理。,148, 23-58 (1999) ·Zbl 0930.76067号
[43] Swarztrauber,P.N.,球面上向量函数及其导数的近似,SIAM J.Numer。分析。,18, 191-210 (1981) ·Zbl 0492.65010号
[44] Swarztrauber,P.N。;威廉姆森,D.L。;Drake,J.B.,求解球面几何中偏微分方程的笛卡尔方法,Dyn。大气。海洋,27679-706(1997)
[45] 泰勒,M。;Tribbia,J。;Iskandarani,M.,球面上浅水方程的谱元法,J.Compute。物理。,130, 92-108 (1997) ·兹伯利0868.76072
[46] Taylor,医学硕士。;Wingate,B.A。;Vicent,R.E.,计算三角形中Fekete点的算法,SIAM J.Numer。分析。,38, 5, 1707-1720 (2000) ·Zbl 0986.65017号
[47] Thuburn,J.,基于Pv的六边形椭球网格浅水模型,蒙大拿州。我们。修订版,125,2328-2347(1997)
[48] Tomita,H。;Tsugawa,M。;佐藤,M。;Goto,K.,利用弹簧动力学在修改的二十面体测地线网格上建立浅水模型,J.Compute。物理。,174, 579-613 (2001) ·Zbl 1056.76058号
[49] 威廉姆森,D.L.,球面测地线网格上正压流的积分,布尔。Am.流星。《社会学杂志》,48,589(1967)
[50] 威廉姆森,D.L。;Drake,J.B。;哈克·J·J。;雅各布,R。;Swarztrauber,P.N.,球面几何中浅水方程数值近似的标准测试集,J.Compute。物理。,102, 211-224 (1992) ·Zbl 0756.76060号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。