×

粘性流体的修正浅水模型和保正数值逼近。 (英语) Zbl 1481.76091号

摘要:对于水平长度尺度远大于垂直长度尺度的地球物理流,浅水方程被广泛应用于模拟。受熔岩流示例的启发,我们在这里考虑了一个修正模型,其中添加了标量(例如温度)的输运方程,并从深度平均的Navier-Stokes方程推导了浅水方程。一些模型变量的恒定垂直剖面假设被放宽,允许存在垂直剖面,因此通量项的非线性导致引入适当的形状系数。
所得到的双曲型偏微分方程组的空间离散化是通过引入的有限体积中心迎风格式的修正版本获得的A.库尔加诺夫G.彼得罗娃【公共数学科学5,第1期,133–160(2007;Zbl 1226.76008号)]. 时间离散化基于隐式显式Runge-Kutta方法,该方法适当耦合双曲部分和刚性源项,避免了使用非常小的时间步长;复数运算的使用提高了刚性项隐式处理的准确性。证明了整个方案保持了流动厚度的正性和稳态。
通过一些数值实验验证了所提方法,并显示了模型中引入的形状系数数值解的关联性。

MSC公司:

76D99型 不可压缩粘性流体
76个M12 有限体积法在流体力学问题中的应用

软件:

冠军
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 谢里丹,M.F。;Stinton,A.J。;Patra,A。;皮特曼,E.B。;Bauer,A。;Nichita,C.C.,使用1963年小塔霍马峰雪崩评估Titan2D质量流模型,华盛顿雷尼尔山,J.火山。地热。第139号、第1-2号、第89-102号决议(2005年)
[2] 科斯塔,A。;Macedonio,G.,基于深度平均方程的熔岩流数值模拟,Geophys-Res-Lett,32,5(2005)·Zbl 1082.76039号
[3] Kelfoun,K。;Vargas,S.V.,VolcFlow模拟熔岩流的能力和潜在发展,地质学会,伦敦,专门出版物,426,1,337-343(2016)
[4] Procter,J.N。;克罗宁,S.J。;富勒,I.C。;谢里丹,M。;尼尔,V.E。;Keys,H.,《利用Titan2D对地貌快速变化的冲积扇进行拉哈尔灾害评估:黄盖湖河,鲁阿佩胡山,地貌学,116,1-2,162-174》(2010)
[5] 威廉姆森,D.L。;Drake,J.B。;哈克·J·J。;雅各布,R。;Swarztrauber,P.N.,球面几何中浅水方程数值近似的标准测试集,《计算物理杂志》,102,1,211-224(1992)·Zbl 0756.76060号
[6] Behrens,J.,浅水方程的自适应有限元求解器的大气和海洋建模,Appl。数字。数学。,26, 1-2, 217-226 (1998) ·Zbl 0897.76046号
[7] Lyard,F。;Lefevre,F。;Letellier,T。;Francis,O.,《全球海潮建模:来自FES2004的现代见解》,ocean Dyn,56,5-6,394-415(2006)
[8] Mader,C.L.,《水波数值模拟》(2004),CRC出版社·Zbl 1081.76001号
[9] 费尔南德斯·尼托,E.D。;Bouchut,F。;布雷施博士。;迪亚兹,M.J.C。;Mangeney,A.,海底雪崩和产生海啸的新型Savage-Hutter模型,《计算物理杂志》,227,16,7720-7754(2008)·Zbl 1156.76015号
[10] Boussinesq,J.,Essai sur la théorie des eaux courates(1877年),Impr。民族的
[11] 霍格,A.J。;Pritchard,D.,水力阻力对溃坝和其他浅层惯性流的影响,《流体力学杂志》,501179-212(2004)·Zbl 1074.76010号
[12] 布拉德福德,S.F。;Sanders,B.F.,任意地形浅水淹没的有限体积模型,J.Hydraul。工程师,128,3289-298(2002)
[13] 库加诺夫,A。;Levy,D.,圣维南系统的中心迎风方案,ESAIM:数学建模和数值分析,36,3,397-425(2002)·Zbl 1137.65398号
[14] 毕,S。;周,J。;刘,Y。;Song,L.,在自适应笛卡尔网格上模拟干湿锋浅层流动的有限体积方法,工程中的数学问题,2014(2014)·Zbl 1407.76082号
[15] Begnudelli,L。;Sanders,B.F.,《四边形网格有限体积模型的保守湿润和干燥方法》,J.Hydraul。工程师,133,312-322(2007)
[16] 梁,Q。;Borthwick,A.G.L.,复杂地形上具有干湿锋的浅层水流自适应四叉树模拟,计算机与流体,38,2,221-234(2009)·Zbl 1237.76017号
[17] 库加诺夫,A。;Petrova,G.,《圣维南体系的二阶平衡正守恒中心迎风方案》,Commun。数学。科学。,5, 1, 133-160 (2007) ·Zbl 1226.76008号
[18] 刘,X。;奥尔布赖特,J。;Epshteynb,Y。;Kurganov,A.,圣维南系统三角网格上新型湿/干重建的保正中心迎风格式,计算物理杂志,374,213-236(2018)·Zbl 1416.65295号
[19] Liu,X.,具有湿/干界面的浅水控制溶质运移的鲁棒数值模型,计算。方法应用。机械。工程,351,85-108(2019)·Zbl 1441.76027号
[20] Bouchut,F.,双曲守恒律有限体积方法的非线性稳定性:源的平衡方案(2004),Springer Science&Business Media·Zbl 1086.65091号
[21] 曼宁,R。;格里菲斯,J.P。;Pigot,T.F。;Vernon Harcourt,L.F.,《明渠和管道中的水流》(1890年)
[22] Russo,G.,《平衡定律的中心方案》(Freistȭhler H.,Warnecke G.,双曲型问题:理论,数值,应用。双曲型:理论,数字,应用,ISNM国际数值数学系列(2001))
[23] Russo,G.,《应用于浅水方程的守恒定律的中心方案》(Rionero S.,R.G.e.,《数学在力学中的趋势和应用》(2005),Springer),225-246
[24] W.Squire。;Trapp,G.,《使用复变量估计实函数导数》。,SIAM版本,40,1,110-112(1998)·Zbl 0913.65014号
[25] Castra Día,M.J。;费尔南德斯·尼托,E.D。;Ferreiro,A.,《虚水方程中的泥沙输移模型和高阶有限体积法的数值方法》,《计算机与流体》,37,3,299-316(2008)·兹比尔1237.76082
[26] 李,S。;Duffy,C.J.,《非结构化网格上浅水流动和污染物传输的完全耦合建模》,《Procedia Environ Sci》,第13、3、2098-2121页(2012年)
[27] de’Michieli Vitturi,M。;Esposti Ongaro,T。;Lari,G。;Aravena,A.,IMEX_SfloW2D 1.0。火山碎屑崩塌的深度平均数值流模型,Geosci。模型开发讨论。(2018)
[28] Wirth,A.,《物理海洋学导览》,Phys。Oceanogr公司。(2015)
[29] 克里格,I.M。;Dougherty,T.J.,刚性球体悬浮液中非牛顿流动的机制,《流变学会学报》,第3期,第1期,第137-152页(1959年)·Zbl 0100.21502号
[30] Dragoni,M。;Borsari,I。;Tallarico,A.,熔岩流锋形状模型,《地球物理研究杂志:固体地球》,110,B9(2005)
[31] 杰博,J。;Perthame,B.,层流浅水粘性圣维南系统的推导:数值验证,离散和连续动力系统-系列B,1,1,89-102(2001)·Zbl 0997.76023号
[32] LeVeque,R.J.,《双曲问题的有限体积方法》(2004),剑桥应用数学教材
[33] Toro,E.F.,《浅水流动中的干床问题》,航空学院,第9007号报告(1990年)
[34] 帕雷希,L。;Russo,G.,刚性微分方程组的隐式显式Runge-Kutta格式,数值分析的最新趋势,269-288(2000),新星科学·兹比尔1018.65093
[35] 哥特利布,S。;舒,C.W。;Tadmor,E.,强稳定性保持高阶时间离散化方法,SIAM Rev.,43,89-112(2001)·Zbl 0967.65098号
[36] 奥尔特加,J.M。;Rheinboldt,W.C.,多变量非线性方程的迭代解(1987),剑桥大学·Zbl 0241.65046号
[37] Lyness,J.N。;Moler,C.B.,解析函数的数值微分。,SIAM J数字分析,4202-210(1967)·Zbl 0155.48003号
[38] Lyness,J.N.,基于复变量理论的数值算法,1967年第22届全国会议论文集,125-133(1967),美国计算机学会
[39] 安德森,W.K。;纽曼,J.C。;Whitfield,D.L。;Nielsen,E.J.,使用复变量对非结构网格上Navier-Stokes方程的敏感性分析。,AIAA杂志,39,1,56-63(2001)
[40] 马丁斯,J.R.R.A。;克鲁,I.M。;Alonso,J.J.,使用复杂变量进行灵敏度分析的自动化方法。,AIAA论文2000-0689,第38届航空航天科学会议,2000-0688(2000)
[41] 马丁斯,J.R.R.A。;Sturdza,P。;Alonso,J.J.,《复阶导数近似》。,ACM事务处理。数学。软件,29,3,245-262(2003)·Zbl 1072.65027号
[42] M.de’Michieli Vitturi,G.Lari,IMEX_SfloW2D v.1.0.0(版本1.0.0),Zenodo,。
[43] Ziv,A.,“相对距离——舍入误差分析中的误差度量”,《数学计算》,39,160,563-569(1982)·兹比尔0533.65021
[44] 新泽西州巴尔姆福思。;Craster,R.V.公司。;佩罗纳,P。;交流生锈。;Sassi,R.,《粘塑性溃坝和bostwick稠度计》,《非牛顿流体力学杂志》,142,63-78(2007)·Zbl 1143.76319号
[45] 萨拉米托,P。;斯穆特克,C。;Cordonnier,B.,《浅层非牛顿流的数值模拟:第一部分:重访一维水平溃坝问题》,国际数值分析与建模杂志,B辑,4283-298(2013)·Zbl 1302.35324号
[46] 科丹尼尔,B。;列夫,E。;Garel,F.,基准熔岩流模型,伦敦地质学会,专门出版物,426425-445(2015)
[47] 李斯特,J.R.,《粘性从点源和线源沿斜面向下流动》,《流体力学杂志》,242631-653(1992)·兹比尔0825.76146
[48] 勒库纳,R。;德尔加多,F。;奥尔蒂斯,A。;卡斯特罗,P.B。;费尔南德斯,I。;Renedo,C.J.,电力变压器替代液体的热流体特性:数值方法,IEEE Trans。电介质。选举人。胰岛素。,22(5), 2522-2529 (2015)
[49] Wang,Y。;Hutter,K。;Pudasaini,S.P.,《Savage-Hutter理论:雪崩、泥石流和泥流的偏微分方程系统》,ZAMM,84(8),507-527(2004)·Zbl 1149.74323号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。