×

波动方程水平显式/垂直隐式(HEVI)解的Runge-Kutta IMEX格式。 (英语) Zbl 1349.86015号

摘要:许多业务天气预报中心使用半隐式时间步长方案,因为其效率很高。然而,随着计算机变得越来越并行,由于隐式方法的额外处理器间通信,大气运动方程的水平显式解可能成为一种有吸引力的替代方案。隐式和显式(IMEX)时间步长方案长期以来一直使用半隐式、分裂-显式或HEVI分裂在大气模型中进行组合。然而,对IMEX格式的精度和稳定性的研究大多局限于对流扩散方程的抛物型情形。我们演示了如何使用一些Runge-Kutta IMEX格式来半隐式或HEVI求解双曲波方程。提出了一种新的HEVI分裂形式UfPreb,它大大提高了分层流中重力波模拟的准确性和稳定性。因此,我们发现,与半隐式格式相比,HEVI格式的精确度不会降低。
找到了梯形隐式和一些Runge-Kutta-IMEX格式的许多变体的稳定性极限,并使用可压缩Boussinesq方程在两个垂直切片情况下对这些格式进行了测试,该方程分为隐式和显式项的各种组合。一些Runge-Kutta格式被发现在梯形上是有益的,特别是因为它们可以在不降低一阶精度的情况下抑制高频。对于刚性系统,我们测试的方案在形式上并不精确,但在刚性极限(几乎不可压缩)下,发现它们可以很好地执行。方案ARK2(2,3,2)在测试中表现最佳。

MSC公司:

86-08 地球物理问题的计算方法
76M20码 有限差分方法在流体力学问题中的应用
6500万06 含偏微分方程初值和初边值问题的有限差分方法
86A10美元 气象学和大气物理学

软件:

开放式泡沫
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 阿舍尔,美国。;Ruuth,S。;Spiteri,R.,含时偏微分方程的隐式显式Runge-Kutta方法,创新时间积分器研讨会,Ctr。数学与公司。科学。。创新时间积分器研讨会,Ctr。数学与公司。科学。,荷兰阿姆斯特丹,1996年10月30日至11月1日。创新时间积分器研讨会,Ctr。数学与公司。科学。。创新时间积分器研讨会,Ctr。数学与公司。科学。,荷兰阿姆斯特丹,1996年10月30日至11月1日,申请。数字。数学。,25, 2-3, 151-167 (1997) ·Zbl 0896.65061号
[2] Asselin,R.,用于时间积分的频率滤波器,周一。《天气评论》,100487-490(1972)
[3] Cullen,M.,全可压缩非静水模型的半隐式积分技术测试,Quart。J.罗伊。流星。Soc.,116,495,1253-1258(1990)
[4] Dahlquist,G.,线性多步方法的一个特殊稳定性问题,BIT-Numer。数学。,3, 1, 27-43 (1963) ·Zbl 0123.11703号
[5] Davies,T。;M.卡伦。;马尔科姆,A。;Mawson,M。;Staniforth,A。;怀特,A。;Wood,N.,气象局全球和区域大气建模的新动力核心,夸特。J.罗伊。流星。《社会学杂志》,131,608,1759-1782(2005)
[6] Durran,D.,地球物理流体动力学中波动方程的数值方法(1999),Springer
[7] Durran,D.,流体动力学数值方法及其在地球物理中的应用,文本应用。数学。,第32卷(2010),施普林格出版社·Zbl 1214.76001号
[8] Durran,D.,稳定快波,(全球大气模型的数值技术。全球大气模型的数值技术,Lect.Notes Comput.Sci.Eng.,vol.80(2011),Springer),第6章·Zbl 1215.86001号
[9] 杜兰,D。;Blossey,P.,快波-慢波问题的隐式-显式多步方法,Mon。《天气评论》,140,4,1307-1325(2012)
[10] Gassmann,A.,全球六边形C网格非静水动力堆芯(ICON-IAP),设计用于能量一致性,夸脱。J.罗伊。流星。Soc.139670152-175(2013年1月)
[11] Giraldo,F.,可扩展光谱元素大气模型的半隐式时间积分器,Quart。J.罗伊。流星。《社会学杂志》,1316102431-2454(2005)
[12] Giraldo,F。;Kelly,J。;Constantinescu,E.,《三维非静力统一大气模型(NUMA)的隐式显式公式》,SIAM J.Sci。计算。(2013),出版中·兹比尔1280.86008
[13] Higdon,R.,《分层海洋环流模型的两级时间步长方法:进一步开发和测试》,J.Compute。物理。,206, 2, 463-504 (2005) ·Zbl 1121.86002号
[14] 肯尼迪,C。;Carpenter,M.,对流-扩散-反应方程的可加Runge-Kutta格式,应用。数字。数学。,44, 139-181 (2003) ·Zbl 1013.65103号
[15] 库普卡,F。;Happenhofer,N。;Higueral,I。;Kock,O.,《天体物理学中双扩散对流模拟的全变分-减小隐式-显式Runge-Kutta方法》,J.Compute。物理。,231, 3561-3586 (2012) ·Zbl 1241.85010号
[17] Mesinger,F.,前向反向方案及其在有限区域模型中的应用,Contrib.Atmos。物理。,50, 200-210 (1977)
[19] Pareschi,L。;Russo,G.,隐式显式Runge-Kutta格式及其在松弛双曲方程组中的应用,J.Sci。计算。,25, 1/2, 129-155 (2005) ·Zbl 1203.65111号
[20] 斯卡马洛克,W。;Klemp,J.,静力学和非静力学弹性方程的时间分割数值方法的稳定性,Mon。《天气评论》,1202109-2127(1992)
[21] Szmelter,J。;Smolarkiewicz,P.,地球物理框架中基于边缘的非结构化网格离散化,J.Compute。物理。,229, 4980-4995 (2010) ·Zbl 1346.76097号
[22] 塔普,M。;White,P.,非静力中尺度模式,夸特。J.罗伊。流星。Soc.,102,432,277-296(1976)
[23] Ullrich,P。;Jablonowski,C.,非静力大气模型的Operator-split Runge-Kutta-Rosenbrock方法,蒙大拿州。《天气评论》,140,4,1257-1284(2012)
[24] Yeh,K.-S。;科特,J。;砾石,S。;Méthot,A。;帕托因,A。;Roch,M。;Staniforth,A.,CMC-MRB全球环境多尺度(GEM)模型。第三部分:非静力公式,周一。《天气评论》,130,2,339-356(2002)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。