×

使用地形跟随坐标,针对二维Euler和Navier-Stokes方程的水平显式、垂直隐式(HEVI)间断Galerkin格式。 (英语) Zbl 07516454号

摘要:提出了一种基于间断伽辽金(DG)方法求解具有任意扩散的欧拉方程的方法。为了适用于从全球尺度到中尺度的大气模拟问题,将水平显式、垂直隐式(HEVI)方法应用于DG离散化,以避免薄网格单元的微小时间步长。为了以正确的近似顺序考虑地形,使用了地形跟随坐标,但为了保持DG的局部守恒特性,方程以强守恒形式表示。IMEX-Runge-Kutta时间积分可以实现时间上的至少三阶近似。以类似的方式,Bassi,Rebay将扩散视为物理过程(湍流参数化)和物理驱动的稳定机制[F.巴斯西S.返利,J.计算。物理学。131,第2期,267–279页(1997年;Zbl 0871.76040号)]以地形坐标和保护方式接近。几个与大气相关的测试案例证明了该方法的有效性。尽管所提出的大多数理论都适用于二维和三维(2D和3D),但这些测试用例都是纯二维的。他们特别表明,拟议方案似乎能够承受非常陡峭的地形。

MSC公司:

6500万 偏微分方程、初值和含时初边值问题的数值方法
7.6亿 流体力学基本方法
86轴 地球物理学

软件:

偶像
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 西卡罗来纳州斯卡马洛克。;Klemp,J.B。;杜迪亚,J。;吉尔,D.O。;刘,Z。;伯纳,J。;Wang,W。;鲍尔斯,J.G。;杜达,M.G。;巴克,D.M。;Huang,X.-Y.,《WRF第4版高级研究说明》(2019年),NCAR:NCAR Boulder,技术代表。
[2] 西卡罗来纳州斯卡莫罗克。;Klemp,J.B.,《用于天气研究和预测应用的时间分割非静力大气模型》,J.Compute。物理。,227, 3465-3485 (2008) ·Zbl 1132.86312号
[3] 西卡罗来纳州斯卡莫罗克。;Klemp,J.B。;杜达,M.G。;福勒,L。;公园,S.-H。;Ringler,T.D.,使用质心Voronoi细分和C网格交错的多尺度非静力大气模型,Mon。《天气评论》,240,3090-3105(2012)
[4] Doms,G。;Baldauf,M.,《非静力区域COSMO模型的描述》,第一部分:动力学和数值(2018年2月),德国Wetterdienst:Deutscher Wetterdeest Offenbach,德国,v5.05
[5] Zängl,G。;Reinert,D。;Ripodas,P。;Baldauf,M.,DWD和MPI-M的ICON(ICOShedral非静水力学)建模框架:非静水动力学核心的描述,Q.J.R.Meteorol。Soc.,141563-579(2015)
[6] Cockburn,B。;Shu,C.-W.,TVB Runge-Kutta局部投影非连续Galerkin守恒定律有限元方法II:一般框架,数学。计算。,52, 186, 411-435 (1989) ·Zbl 0662.65083号
[7] Cockburn,B。;Lin,S.-Y。;Shu,C.-W.,TVB Runge-Kutta局部投影非连续Galerkin有限元方法,守恒定律III:一维系统,J.Compute。物理。,84, 90-113 (1989) ·Zbl 0677.65093号
[8] Cockburn,B。;Hou,S。;Shu,C.-W.,守恒定律的Runge-Kutta局部投影间断Galerkin有限元方法IV:多维情况,数学。计算。,54, 190, 545-581 (1990) ·兹伯利0695.65066
[9] 赫塞文,J.S。;Warburton,T.,Nodal间断Galerkin方法(2008),Springer·Zbl 1134.65068号
[10] Marras,S。;凯利,J.F。;莫拉格,M。;米勒,A。;Kopera,医学硕士。;巴斯克斯,M。;吉拉尔多·F·X。;Houzeaux,G。;Jorba,O.,数值天气预测的基于单元的伽辽金方法综述:有限元、谱元和不连续伽辽金,Arch。计算。方法工程,23,673-722(2016)·Zbl 1360.86004号
[11] Bao,L。;Kloefkorn,R。;Nair,R.D.,非连续Galerkin非静力模型的水平显式和垂直隐式(HEVI)时间离散格式,Mon。《天气评论》,143972-990(2015)
[12] 奈尔·R·D。;Bao,L。;玩具,医学博士。;Kloefkorn,R.,高阶多尺度全球大气模式,AIAA Aviat。,1-13 (2016)
[13] Pareschi,L。;Russo,G.,隐式显式Runge-Kutta格式及其在带松弛双曲方程组中的应用,J.Sci。计算。,25, 1/2, 129-155 (2005) ·Zbl 1203.65111号
[14] 锁,S.-J。;木材,N。;Weller,H.,大气模型水平显式和垂直隐式解的Runge-Kutta隐式显式格式的数值分析,Q.J.R.Meteorol。Soc.,140,1654-1669(2014)
[15] 吉拉尔多·F·X。;Kelly,J。;Constantinescu,E.M.,《三维非静力统一大气模型(NUMA)的隐式显式公式》,SIAM J.Sci。计算。,35、5、B1162-B1194(2013)·Zbl 1280.86008号
[16] 布莱斯,S。;兰布雷希茨,J。;Deleersnijder,E.,《应用于非静力大气模拟的三维间断Galerkin离散化的稳定性》,国际期刊Numer。《液体方法》,81,9,558-585(2016)
[17] Orgis,T。;Läuter,M。;Handorf,D。;Dethloff,K.,使用低阶间断Galerkin离散的β平面斜压波,J.Compute。物理。,339, 461-481 (2017)
[18] Abdi,D.S。;吉拉尔多·F·X。;Constantinescu,E.M。;卡尔,L.E。;Wilcox,L.C.公司。;Warburton,T.C.,在多核处理器上加速隐式-显式非静力统一大气模型(NUMA),国际高性能计算杂志。申请。,33, 2, 242-267 (2019)
[19] 克莱因,R。;贝茨,K.R。;Nikiforakis,N.,《大气流动的井平衡可压缩切割细胞模拟》,Phil.Trans。R.Soc.A,367,4559-4575(2009)·Zbl 1192.86013号
[20] 很好,B。;Gadian,A。;锁,S.-J。;Ross,A.,在理想模拟中表示地形的切分法的性能,Atmos。科学。莱特。,15, 44-49 (2014)
[21] 肖,J。;Weller,H.,使用非流体静力模型比较地形跟随网格和切割网格,Mon。天气评论,1442085-2099(2016)
[22] 巴尔达夫,M。;Brdar,S.,《地形坐标下的三维扩散:水平显式、垂直隐式离散化的测试和稳定性》,Q.J.R.Meteorol。Soc.,1422087-2101(2016)
[23] Stephani,H.,《广义相对论——引力场理论导论》(1990),剑桥大学出版社:剑桥大学出版社·Zbl 0712.53031号
[24] 查伦,M。;Zadra,A。;Girard,C.,《外部引力场中经典流体的四维张量方程》,Q.J.R.Meteorol。Soc.,140908-916(2014)
[25] Jorgensen,B.H.,一般坐标下不可压缩流强守恒形式模型方程的张量公式(2003年12月),里索国家实验室:里索国家实验,丹麦罗斯基勒,技术代表。
[26] Baldauf,M.,球面和椭球体上协变形式浅水方程的间断Galerkin解算器,J.Compute。物理。,410, 109384, 1-26 (2020) ·Zbl 1436.65136号
[27] 吉拉尔多·F·X。;Restelli,M。;Läuter,M.,《Navier-Stokes方程的半隐式公式:在非静力大气建模中的应用》,SIAM J.Sci。计算。,32, 6, 3394-3425 (2010) ·Zbl 1237.76153号
[28] Landau,L.D。;Lifshitz,E.M.,《流体力学,理论物理课程》,第6卷(1987),佩加蒙出版社,iSBN 978-0-08-033933-7·Zbl 0655.76001号
[29] Smagorinsky,J.,《原始方程的一般循环实验》,孟买。《天气评论》,91,3,99-164(1963)
[30] Pope,S.B.,《湍流》(2001),剑桥大学出版社:剑桥大学出版社
[31] Clark,T.L.,使用地形跟随坐标变换的小规模动态模型,J.Compute。物理。,24, 186-215 (1977) ·Zbl 0369.76025号
[32] 北卡罗来纳州韦迪。;Smolarkiewicz,P.K.,在时间相关曲线边界上进行坐标变换后扩展Gal-Chen和Somerville地形,J.Comput。物理。,193, 1-20 (2003) ·Zbl 1117.76307号
[33] 舒斯特,D。;布尔达尔,S。;巴尔达夫,M。;Dedner,A。;Klöfkorn,R。;Kröner,D.,《关于有和无水分的中尺度大气流动的间断Galerkin方法》,Meteorol。Z.,23,4,449-464(2014)
[34] Bassi,Rebay,数值求解可压缩Navier-Stokes方程的高精度间断有限元方法,J.Compute。物理。,131, 267-279 (1997) ·兹比尔0871.76040
[35] 吉拉尔多·F·X。;赫塞文,J.S。;Warburton,T.,球面浅水方程的Nodal高阶间断Galerkin方法,J.Compute。物理。,181, 499-525 (2002) ·Zbl 1178.76268号
[36] Läuter,M。;吉拉尔多·F·X。;Handorf,D。;Dethloff,K.,球面三角坐标系中浅水方程的间断Galerkin方法,J.Compute。物理。,227, 10226-10242 (2008) ·Zbl 1218.76028号
[37] Kirby,R.M。;Karniadakis,G.E.,《非均匀网格上的去混叠:算法和应用》,J.Compute。物理。,191, 249-264 (2003) ·Zbl 1161.76534号
[38] Gassner,G.J。;Beck,A.D.,关于欠分辨率湍流模拟的高阶离散精度,Theor。计算。流体动力学。,27, 221-237 (2013)
[39] 韦勒,H。;锁,S.-J。;Wood,N.,波动方程水平显式/垂直隐式(HEVI)解的Runge-Kutta IMEX格式,J.Compute。物理。,252, 365-381 (2013) ·Zbl 1349.86015号
[40] Cockburn,B。;Shu,C.-W.,含时对流扩散系统的局部间断Galerkin方法,SIAM J.Numer。分析。,35, 6, 2440-2463 (1998) ·Zbl 0927.65118号
[41] 斯特拉卡,J.M。;威廉姆森,R.B。;威克,L.J。;安德森,J.R。;Droegemeier,K.K.,《非线性密度电流的数值解:基准解和比较》,《国际数值杂志》。液体方法,17,1-22(1993)
[42] Davies,H.C.,多级预测模型的横向边界公式,Q.J.R.Meteorol。Soc.,102,432,405-418(1976)
[43] Klemp,J.B。;杜迪亚,J。;Hassiotis,A.D.,NWP应用的上部重力波吸收层,周一。《天气评论》,1363987-4004(2008)
[44] Baldauf,M.,Runge-Kutta时间积分平流方程线性离散化的稳定性分析,J.计算。物理。,227, 6638-6659 (2008) ·Zbl 1149.65064号
[45] 西卡罗来纳州斯卡莫罗克。;Klemp,J.B.,《Klemp-Wilhelmson时间分裂方案的效率和准确性》,周一。《天气评论》,1222623-2630(1994)
[46] 巴尔达夫,M。;Brdar,S.,《河道中线性重力波的解析解,作为使用非静水压可压缩Euler方程的数值模型的测试》,Q.J.R.Meteorol。Soc.,1391977-1989(2013)
[47] Schär,C。;Leuenberger,D。;福勒,O。;吕蒂,D。;Girard,C.,大气预测模型的新地形垂直坐标公式,Mon。《天气评论》,130,10,2459-2480(2002)
[48] Baldauf,M.,山区流量的线性解及其与COSMO模型的比较,COSMO-Newsl。,9, 19-24 (2008)
[49] Smith,R.B.,通过孤立山体的分层静水流线性理论,Tellus,32348-364(1980)
[50] Zängl,G.,《在陡坡上扩展地形跟随坐标模型的数值稳定性极限》,Mon。《天气评论》,140、11、3722-3733(2012)
[51] Tachikawa,T。;卡库,M。;川崎,A。;Gesch,D.B。;Oimoen,M.J。;张,Z。;丹尼尔森·J·J。;Krieger,T。;柯蒂斯,B。;Haase,J。;艾布拉姆斯,M。;Carabajal,C.,ASTER全球数字高程模型第2版-验证结果摘要(2011年),NASA,地球资源观测与科学(EROS)中心,联邦政府丛书,技术代表。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。