×

耦合可压缩Navier-Stokes方程的质量守恒隐式显式方法。 (英语) Zbl 1506.76149号

概述:地球系统模型由耦合组件组成,这些组件分别对全球大气、海洋和陆地表面等系统进行建模。虽然这些组件开发得很好,但将它们耦合到单个系统中可能是一项重大挑战。计算效率、准确性和稳定性是主要关注点。在本研究中,我们关注这些问题。特别地,研究了用于处理不同时间尺度的隐式显式(IMEX)紧耦合和松耦合策略。对于气-气相互作用问题的简化模型,我们考虑带界面条件的耦合可压缩Navier-Stokes方程。在刚性线假设下,水平动量和热流通过界面进行交换。通过数值实验验证了耦合方案的稳定性。我们在数值和理论上都表明,对于具有刚性盖条件的耦合可压缩Navier-Stokes系统,我们的IMEX耦合方法是质量守恒的。

MSC公司:

76号06 可压缩Navier-Stokes方程
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Golaz,J.-C。;考德威尔,P.M。;Van Roekel,L.P。;彼得森,M.R。;唐奇。;Wolfe,J.D。;Abeshu,G。;阿南塔拉吉,V。;Asay-Davis,X.S。;Bader,D.C.,《DOE E3SM耦合模型版本1:标准分辨率下的概述和评估》,J.Adv.建模地球系统。,11, 7, 2089-2129 (2019)
[2] 黑泽勒格,W。;Severijns,C。;塞姆勒,T。;cStefuanescu,S。;Yang,S。;王,X。;Wyser,K。;杜特拉,E。;巴尔达萨诺,J.M。;Bintanja,R.,《EC-Earth:无缝地球系统预测方法的实际应用》,公牛。Am.Meteorol公司。Soc.,91,10,1357-1364(2010年)
[3] 庆大霉素,P。;格林,J.K。;盖林,M。;汉弗莱,V。;Seneviratne,S.I。;Zhang,Y。;周,S.,《陆地碳和水循环之间的耦合综述》,环境。Res.Lett.公司。,第14、8条,第083003页(2019年)
[4] Burrows,S。;马尔特鲁德,M。;Yang,X.先生。;朱,Q。;Jeffery,N。;施,X。;Ricciuto,D。;王,S。;比什特,G。;Tang,J.,DOE E3SM v1。1生物地球化学结构:描述和模拟生态系统气候对历史强迫变化的响应,J.Adv.建模地球系统。,第12、9条,第2019MS001766页(2020年)
[5] 柯林斯,W.D。;比茨,C.M。;布莱克蒙,M.L。;Bonan,G.B。;Bretherton,C.S。;Carton,J.A。;Chang,P。;南卡罗来纳州多尼。;哈克·J·J。;Henderson,T.B.,社区气候系统模型第3版(CCSM3),J.Clim。,19, 11, 2122-2143 (2006)
[6] 安德森·T·R。;霍金斯,E。;Jones,P.D.,《二氧化碳、温室效应和全球变暖:从Arrhenius和Callendar的开创性工作到今天的地球系统模型》,Endeavour,40,3,178-187(2016)
[7] 霍夫曼,M.J。;Asay-Davis,X。;价格,S.F。;Fyke,J。;Perego,M.,海底融化变化对Thwaites冰川海平面上升贡献的影响,南极洲,地球物理杂志。研究:地球表面(2019)
[8] 克雷格,A.P。;Vertenstein,M。;Jacob,R.,为CCSM4和CESM1开发的用于地球系统建模的新型挠性耦合器,国际高性能计算杂志。申请。,26, 1, 31-42 (2012)
[9] 雅各布,R。;Larson,J。;Ong,E.,《使用模型耦合工具包的社区气候系统模型第3版中的M(时间N)通信和并行插值》,国际期刊《高性能计算》。申请。,19, 3, 293-307 (2005)
[10] Valcke,S.,《OASIS3耦合器:欧洲气候建模社区软件》,Geosci。模型开发,6,2,373(2013)
[11] R.Hallberg,冰/海洋耦合系统的数值不稳定性,载于:CLIVAR WGOMD High研讨会,2014年,第38页。
[12] Lemarié,F。;布莱奥,E。;Debreu,L.,《海洋-大气耦合算法分析:一致性和稳定性》,Procedia Compute。科学。,51, 2066-2075 (2015)
[13] Beljaars,A。;杜特拉,E。;Balsamo,G。;Lemarié,F.,关于天气和气候模型中表面-大气耦合的数值稳定性,Geosci。模型开发,10,2,977-989(2017)
[14] 张,H。;刘,Z。;Constantinescu,E。;Jacob,R.,海洋-大气耦合界面条件的稳定性分析,科学杂志。计算。,84, 3, 1-25 (2020) ·兹比尔1452.65183
[15] 彼得森,K。;Bochev,P。;Kuberry,P.,基于拉格朗日乘子的界面问题显式同步分区算法,计算。数学。申请。,78, 2, 459-482 (2019) ·Zbl 1442.65148号
[16] Sockwell,K.C。;彼得森,K。;Kuberry,P。;Bochev,P。;Trask,N.,海洋-大气系统界面通量恢复耦合方法,结果应用。数学。,第100110条pp.(2020)·Zbl 1443.76239号
[17] Lemarié,F。;Marchesello,P。;德布鲁,L。;Blayo,E.,海洋-大气耦合模型对耦合方法的敏感性:热带气旋Erica的例子(2014)
[18] Bresch,D。;Koko,J.,两种耦合Navier-Stokes流体数值模拟的算子分裂和拉格朗日乘子区域分解方法,国际期刊应用。数学。计算。科学。,16419-429(2006年)·Zbl 1119.35346号
[19] Formaggia,L。;杰布·J·F。;Nobile,F。;Quarteroni,A.,关于柔顺容器中流动问题的3D和1D Navier-Stokes方程耦合,计算。方法应用。机械。工程,191,6-7,561-582(2001)·Zbl 1007.74035号
[20] 康纳斯,J.M。;Dolan,R.D.,两种保守的高阶流体-流体耦合方法的稳定性,高级应用。数学。机械。,11, 6, 1287-1338 (2019) ·Zbl 1488.65415号
[21] Carpenter,M.H。;费希尔,T.C。;尼尔森,E.J。;Frankel,S.H.,Navier-Stokes方程的熵稳定谱配置格式:不连续界面,SIAM J.Sci。计算。,36、5、B835-B867(2014)·Zbl 1457.65140号
[22] 阿舍尔,U.M。;Ruuth,S.J。;Spiteri,R.J.,含时偏微分方程的隐式显式Runge-Kutta方法,应用。数字。数学。,25, 2, 151-167 (1997) ·Zbl 0896.65061号
[23] Pareschi,L。;Russo,G.,隐式显式Runge-Kutta格式及其在带松弛双曲方程组中的应用,J.Sci。计算。,25, 1-2, 129-155 (2005) ·Zbl 1203.65111号
[24] Constantinescu,E。;Sandu,A.,外推IMplicit-EXplicit时间步进,SIAM J.Sci。计算。,31, 6, 4452-4477 (2010) ·Zbl 1209.65069号
[25] 肯尼迪,C.A。;Carpenter,M.H.,对流-扩散-反应方程的可加Runge-Kutta格式,应用。数字。数学。,44, 1-2, 139-181 (2003) ·Zbl 1013.65103号
[26] Roldan,T。;Higueras,I.,低存储要求的高效隐式显式Runge-Kutta方法,SciCADE 2013,45,1,174(2013)
[27] 博斯卡里诺,S。;Pareschi,L。;Russo,G.,《具有多尺度松弛的双曲型系统的统一IMEX Runge-Kutta方法》,SIAM J.Numer。分析。,2008年4月55日至109日(2017年)·Zbl 1372.65261号
[28] Kanevsky,A。;Carpenter,M.H。;Gottlieb,D。;Hesthaven,J.S.,隐式显式高阶Runge-Kutta方法在间断Galerkin格式中的应用,J.Compute。物理。,225, 2, 1753-1781 (2007) ·Zbl 1123.65097号
[29] Kang,S。;吉拉尔多·F·X。;Bui-Thanh,T.,IMEX HDG-DG:浅水系统的耦合隐式混合间断Galerkin和显式间断Galergin方法,J.Compute。物理。,第109010条pp.(2019)·Zbl 1453.65328号
[30] Restelli,M。;Giraldo,F.X.,非静力中尺度模拟中Navier-Stokes方程的保守间断Galerkin半隐式公式,SIAM J.Sci。计算。,31, 3, 2231-2257 (2009) ·Zbl 1405.65127号
[31] Giraldo,F。;Restelli,M.,三角间断Galerkin海洋浅水模型的高阶半隐式时间积分器,国际。J.数字。《液体方法》,63,9,1077-1102(2010)·Zbl 1267.76010号
[32] 吉拉尔多·F·X。;Kelly,J.F。;Constantinescu,E.,《三维非静力统一大气模型(NUMA)的隐式显式公式》,SIAM J.Sci。计算。,35、5、B1162-B1194(2013)·兹比尔1280.86008
[33] 加德纳·D·J。;Guerra,J.E。;哈蒙,F.P。;雷诺兹,D.R。;Ullrich,P.A。;Woodward,C.S.,非静力大气模型的隐式显式(IMEX)Runge-Kutta方法,《地球科学》。模型开发,11,4,1497-1515(2018)
[34] Vogl,C.J。;Steyer,A。;雷诺兹,D.R。;Ullrich,P.A。;Woodward,C.S.,HOMME-NH动力堆芯隐显加性Runge-Kutta积分器的评估,J.高级建模地球系统。,1124228-4244(2019)
[35] Abdi,D.S。;吉拉尔多·F·X。;Constantinescu,E.M。;卡尔·L·E。;Wilcox,L.C.公司。;Warburton,T.C.,在多核处理器上加速IMplicit EXplicit大气非流体静力统一模型,Int.J.高性能计算。申请。,33, 2, 242-267 (2019)
[36] Froehle,B。;Persson,P.-O.,高效隐式-显式时间步长的流体-结构相互作用的高阶间断Galerkin方法,J.Compute。物理。,272, 455-470 (2014) ·Zbl 1349.76203号
[37] 黄,D.Z。;佩尔松,P.-O。;Zahr,M.J.,基于隐式显式Runge-Kutta格式的一般多物理问题的高阶线性稳定分区解算器,计算。方法应用。机械。工程,346674-706(2019)·Zbl 1440.65072号
[38] 华盛顿,W.M。;帕金森,C.,《三维气候建模导论》(2005),大学科学书籍
[39] Müller,P.,《海洋运动方程》(2006),剑桥大学出版社·Zbl 1134.76001号
[40] 康纳斯,J.M。;豪厄尔,J.S。;Layton,W.J.,流体-流体相互作用的解耦时间步进方法,SIAM J.Numer。分析。,50, 3, 1297-1319 (2012) ·Zbl 1457.65102号
[41] Liu,W.T。;Katsaros,K.B。;Businger,J.A.,包括界面分子约束在内的热和水蒸气的气-气交换的整体参数化,J.Atmos。科学。,36, 9, 1722-1735 (1979)
[42] Smith,S.D.,《海面风应力、热通量和风廓线系数与风速和温度的关系》,J.Geophys。研究:海洋,93,C12,15467-15472(1988)
[43] Fairall,C.W。;布拉德利,E.F。;罗杰斯,D.P。;埃德森,J.B。;Young,G.S.,热带海洋-全球大气耦合海洋-大气响应实验的大气-海洋通量的整体参数化,J.Geophys。研究:海洋,101,C2,3747-3764(1996)
[44] Bao,J。;Wilczak,J。;Choi,J。;Kantha,L.,使用耦合模型对大风条件下的气-气相互作用进行数值模拟:飓风发展研究,孟买。《天气评论》,128,7,2190-2210(2000)
[45] Panosfsky,H。;Dutton,J.,《大气湍流:工程应用的模型和方法》(1984),John Wiley&Sons
[46] 维氏,D。;马哈特,L。;Andreas,E.L.,《用平均风和整体稳定性表示海面摩擦速度的公式》,J.Appl。美托洛尔。气候。,54, 3, 691-703 (2015)
[47] 雅各布斯,G.B。;科普里瓦,D.A。;Mashayek,F.,可压缩Navier-Stokes方程的守恒等温壁面边界条件,J.Sci。计算。,30, 2, 177-192 (2007) ·Zbl 1109.76046号
[48] Roe,P.,欧拉方程基于特征的方案,年。Rev.流体机械。,18, 1, 337-365 (1986) ·Zbl 0624.76093号
[49] Toro,E.F.,Riemann Solvers and Numerical Methods for Fluid Dynamics:A Practical Introduction(2013),施普林格科学与商业媒体
[50] H.Nishikawa,将扩散方案扩展到Navier-Stokes方案的两种方法:梯度公式或迎风通量,载于:第20届AIAA计算流体动力学会议,2011年,第3044页。
[51] Syrakos,A。;Varchanis,S。;Dimakopoulos,Y。;Goulas,A。;Tsamopoulos,J.,有限体积法中梯度算子离散化的一些流行方法的批判性分析,Phys。流体,29,12,第127103条pp.(2017)
[52] 克里斯蒂尼,P。;Komatitsch,D.,《在海洋声学中使用光谱元素方法的一些示例》,J.Acoust。《美国社会》,131,3,EL229-EL235(2012)
[53] Constantinescu,E.M。;Sandu,A.,双曲守恒律的多速率时间步长方法,科学杂志。计算。,33, 3, 239-278 (2007) ·Zbl 1127.76033号
[54] 施莱格尔,M。;诺思,O。;阿诺德,M。;Wolke,R.,对流方程的多速率Runge-Kutta格式,J.Compute。申请。数学。,226, 2, 345-357 (2009) ·Zbl 1167.65428号
[55] Constantinescu,E.M。;Sandu,A.,《多时间尺度微分方程的外推多速率方法》,《科学杂志》。计算。,56, 1, 28-44 (2013) ·Zbl 1269.65069号
[56] 塞尼,B。;兰布雷奇斯,J。;康布伦,R。;Legat,V。;Remacle,J.-F.,用于加速显式不连续Galerkin计算并应用于地球物理流的多速率时间步进,国际。J.数字。《液体方法》,71,1,41-64(2013)·Zbl 1431.65178号
[57] Sandu,A.,一类多速率无穷小GARK方法,SIAM J.Numer。分析。,57, 5, 2300-2327 (2019) ·Zbl 1435.65109号
[58] 韦勒,H。;锁,S.-J。;Wood,N.,波动方程水平显式/垂直隐式(HEVI)解的Runge-kutta IMEX格式,J.Compute。物理。,252365-381(2013年)·Zbl 1349.86015号
[59] Thuburn,J.,《数值预报和气候模型动力核心的一些守恒问题》,J.Compute。物理。,227, 7, 3715-3730 (2008) ·Zbl 1132.86314号
[60] Ghosh,D。;Constantinescu,E.M.,《基于特征通量分割的大气流半隐式时间积分》,SIAM J.Sci。计算。,38、3、A1848-A1875(2016)·Zbl 1419.76478号
[61] 泰勒,G.I。;Green,A.E.,大漩涡产生小漩涡的机理,Proc。R.Soc.A,158,895,499-521(1937)·JFM 63.1358.03标准
[62] Young,W.R.,《动态焓、保守温度和海水Boussinesq近似》,J.Phys。海洋学家。,40394-400(2010年)
[63] 文莱,C.-H。;Saad,M.,《重访二维眼睑驱动腔问题》,计算与流体,35,3,326-348(2006)·Zbl 1099.76043号
[64] Drazin,P.G。;Reid,W.H.,《流体动力稳定性》(2004),剑桥大学出版社·Zbl 1055.76001号
[65] Springel,V.,E pur si muove:运动网格上的伽利略非变宇宙流体动力学模拟,Mon。不是。R.阿斯顿。Soc.,401,2791-851(2010年)
[66] Lecoanet,D。;McCourt,M。;Quataert,E。;Burns,K.J。;瓦西尔,G.M。;Oishi,J.S。;Brown,B.P。;斯通,J.M。;O'Leary,R.M.,《数值流体力学的经验证非线性Kelvin-Helmholtz基准》,蒙大拿州。不是。R.阿斯顿。Soc.,455,4,4274-4288(2016)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。