×

非结构网格上不可压缩流动的一类新的半隐式有限体积/虚拟元方法。 (英语) Zbl 07725685号

摘要:我们引入了一类新的高精度半隐式格式,用于求解非结构化多边形网格上非线性含时偏微分方程组(PDE)。时间离散化是基于显式项和隐式项之间的分裂,这些分裂可能来自控制方程的多尺度性质,包括慢尺度和快尺度,或者在投影方法的背景下,其中数值解被投影到物理意义的解流形上。我们建议对显式项使用高阶有限体积(FV)格式,从而确保激波的守恒性和鲁棒性,而使用虚拟元方法(VEM)处理隐式项的离散化,这通常需要解决椭圆问题。然后通过适当的\(\mathrm传递数值解{五十} _2\)从FV到VEM解决方案空间和vice-versa的投影运算符。然后使用半隐式IMEX Runge-Kutta格式实现了高阶时间精度,并且证明了新方案是渐近保持(AP)和良好平衡(WB)的。作为代表性模型,我们选择了浅水方程(SWE),从而处理了具有不同弗劳德数的多个时间尺度,以及不可压缩的Navier-Stokes方程(INS),这些方程是借助投影方法求解的,以满足速度场的螺线约束。此外,还设计了基于VEM技术的INS模型的粘性项隐式离散化。因此,关于最大容许时间步长的CFL型稳定性条件仅基于流体速度,而不基于速度或粘性特征值。一组大型测试案例证明了新方案系列的准确性和能力,可以解决不可压缩流体领域的相关基准问题。

MSC公司:

76个M12 有限体积法在流体力学问题中的应用
6500万08 含偏微分方程初值和初边值问题的有限体积法

软件:

HE-E1GODF公司
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 阿舍尔,U.M。;Ruuth,S.J。;Spiteri,R.J.,含时偏微分方程的隐式显式Runge-Kutta方法,应用。数字。数学。,25, 151-167 (1982) ·Zbl 0896.65061号
[2] 博斯卡里诺,S。;Pareschi,L.,关于双曲平衡律IMEX Runge-Kutta格式的渐近性质,J.Compute。申请。数学。,316, 60-73 (2017) ·Zbl 1375.65119号
[3] 帕雷希,L。;Russo,G.,隐式显式Runge-Kutta格式及其在带松弛双曲方程组中的应用,J.Sci。计算。,25, 129-155 (2005) ·Zbl 1203.65111号
[4] 博斯卡里诺,S。;费尔贝特,F。;Russo,G.,含时偏微分方程的高阶半隐式格式,J.Sci。计算。,68, 975-1001 (2016) ·Zbl 1353.65075号
[5] Casulli,V.,《二维浅水方程的半隐式有限差分方法》,J.Compute。物理。,86, 56-74 (1990) ·Zbl 0681.76022号
[6] 卡苏利,V。;Cheng,R.,三维浅水流动的半隐式有限差分方法,国际。J.数字。液体方法,15629-648(1992)·Zbl 0762.76068号
[7] Dumbser,M。;Casulli,V.,浅水方程的交错半隐式谱间断Galerkin格式,Appl。数学。计算。,219, 8057-8077 (2013) ·Zbl 1366.76050号
[8] Park,J。;Munz,C.-D.,《所有马赫数下流体流动的多压力变量方法》,国际。J.数字。《液体方法》,49,905-931(2005)·Zbl 1170.76342号
[9] Boscheri,W。;Pareschi,L.,所有马赫数下3D Navier-Stokes方程的基于高阶压力的半隐式IMEX格式,J.Compute。物理。,434,第110206条pp.(2021)·Zbl 07508522号
[10] Boscheri,W。;Tavelli,M.,三维粘性可压缩流动的高阶半隐式格式,应用。数学。计算。,434,第127457条pp.(2022)·Zbl 1510.76036号
[11] Fambri,F.,粘性和电阻磁流体力学的一种新的保持结构的半隐式有限体积方法,国际。J.数字。液体方法,93,3447-3489(2021)
[12] Ioriati,M。;Dumbser,M.,弹性管中轴对称粘性可压缩流动的半隐式交错间断Galerkin格式,计算与流体,167166-179(2018)·Zbl 1390.76322号
[13] Ioriati,M。;Dumbser,M.,浅水方程交错半隐式间断Galerkin格式的后验子单元有限体积限制,应用。数字。数学。,135, 443-480 (2019) ·Zbl 1403.76046号
[14] Casulli,V.,非静水自由表面流的半隐式有限差分方法,Int.J.Numer。液体方法,30,425-440(1999)·Zbl 0944.76050号
[15] Tumolo,G。;Bonaventura,L。;Restelli,M.,《浅水方程的半隐式、半拉格朗日、p-自适应间断Galerkin方法》,J.Compute。物理。,232, 46-67 (2013) ·Zbl 1291.65305号
[16] Tumolo,G。;Bonaventura,L.,自适应数值天气预测的半隐式、半拉格朗日间断伽辽金框架,Q.J.R.Meteorol。Soc.,141,692,2582-2601(2015)
[17] 奥兰多,G。;Benacchio,T。;Bonaventura,L.,《自适应网格细化大气动力学模拟的IMEX-DG解算器》,J.Compute。申请。数学。,427,第115124条pp.(2023)·Zbl 07698169号
[18] Fambri,F。;Dubser,M.,自适应交错笛卡尔网格上不可压缩Navier-Stokes方程的半隐式间断Galerkin方法,计算。方法应用。机械。工程,324170-203(2017)·Zbl 1439.76066号
[19] Boscheri,W。;Dumbser,M。;Righetti,M.,非结构化Voronoi网格上三维自由表面流动的高阶速度重建半隐式格式,国际。J.数字。《液体方法》,72,607-631(2013)·Zbl 1455.76104号
[20] Boscheri,W。;Pisaturo,G。;Righetti,M.,非结构化Voronoi网格上自由表面流动的高阶无发散速度重建,国际。J.数字。《液体方法》,90,296-321(2019)
[21] Boscheri,W.,用于自由表面流动的交错Voronoi网格上的时空半拉格朗日平流格式,计算与流体,202,第104503条pp.(2020)·Zbl 1519.65033号
[22] 塔维利,M。;Dumbser,M.,交错非结构网格上二维浅水方程的高阶半隐式间断Galerkin方法,应用。数学。计算。,234, 623-644 (2014) ·Zbl 1298.76120号
[23] 塔维利,M。;Dumbser,M.,非结构化四面体网格上三维不可压缩Navier-Stokes方程的交错时空间断Galerkin方法,J.Compute。物理。,319, 294-323 (2016) ·Zbl 1349.76271号
[24] 奥兰多,G。;巴班特,P.F。;Bonaventura,L.,非理想气体可压缩Navier-Stokes方程的有效IMEX-DG求解器,J.Comput。物理。,471,第111653条pp.(2022)·Zbl 07605612号
[25] 巴斯托,S。;费林,J。;托罗,E。;Vázquez-Cendón,M.,《不可压缩湍流的投影混合高阶有限体积/有限元法》,J.Compute。物理。,353, 169-192 (2018) ·Zbl 1380.76052号
[26] 贝穆德斯,A。;巴斯托,S。;Dumbser,M。;费林,J。;Saavedra,L。;Vázquez-Cendón,M.,弱可压缩流动的交错半隐式FV/FE混合投影方法,J.Compute。物理。,421,第109743条pp.(2020)·兹伯利07508366
[27] 巴斯托,S。;Río,L.D。;Vázquez-Cendón,M。;Dumbser,M.,交错非结构网格上所有马赫数流动的半隐式混合有限体积/有限元格式,应用。数学。计算。,402,第126117条pp.(2021)·Zbl 1510.76082号
[28] 巴斯托,S。;Dumbser,M.,适用于所有Froude数的浅水方程的交错半隐式混合有限体积/有限元格式,应用。数字。数学。,175, 108-132 (2022) ·Zbl 1485.76056号
[29] Boscheri,W。;塔维利,M。;Castro,C.,非结构化Voronoi网格上浅水方程的全Froude高阶IMEX格式,应用。数字。数学。,185, 311-335 (2023) ·Zbl 07699012号
[30] 路易斯安那州贝朗·达维加。;布雷齐,F。;Cangiani,A。;马里尼,L。;Manzini,G。;Russo,A.,《虚拟元素方法的基本原理》,数学。模型方法应用。科学。,23, 1, 199-214 (2013) ·兹伯利1416.65433
[31] 艾哈迈德,B。;Alsadei,A。;布雷齐,F。;马里尼,L。;Russo,A.,虚拟元素方法的等效投影仪,计算。数学。申请。,66, 3, 376-391 (2013) ·兹比尔1347.65172
[32] 路易斯安那州贝朗·达维加。;Manzini,G.,《具有任意正则性的虚拟单元法》,IMA J.Numer。分析。,34, 2, 759-781 (2013) ·Zbl 1293.65146号
[33] 路易斯安那州贝朗·达维加。;布雷齐,F。;Marini,D。;Russo,A.,多边形网格上一般二阶椭圆问题的虚拟单元方法,数学。模型方法应用。科学。,26, 729-750 (2016) ·Zbl 1332.65162号
[34] Lipnikov,K。;Manzini,G。;Shashkov,M.,《模拟有限差分法》,J.Compute。物理。,257,B部分,1163-1227(2014)·Zbl 1352.65420号
[35] 路易斯安那州贝朗·达维加。;布雷齐,F。;Marini,L.,线性弹性问题的虚拟元,SIAM J.Numer。分析。,51, 2, 794-812 (2013) ·Zbl 1268.74010号
[36] 增益A。;Talischi,C。;Paulino,G.,关于任意多面体网格上三维线性弹性问题的虚拟元方法,计算。方法应用。机械。工程,282132-160(2014)·Zbl 1423.74095号
[37] Artioli,E。;da Veiga,L.B。;Dassi,F.,《二维实体力学应用的曲线虚拟元素》,计算。方法应用。机械。工程,359,第112667条pp.(2020)·Zbl 1441.74229号
[38] 达西,F。;罗瓦迪纳,C。;Visinoni,M.,线性弹性问题的三维hellinger-Reissner虚拟单元法,计算。方法应用。机械。工程,364,第112910条pp.(2020)·Zbl 1442.65358号
[39] Chi,H。;da Veiga,L.B。;Paulino,G.,有限变形虚拟单元法VEM的一些基本公式,计算。方法应用。机械。工程,318142-190(2017)·Zbl 1439.74397号
[40] Antonietti,P。;Manzini,G。;Mazzieri,I。;穆拉德,H。;Verani,M.,线性弹性动力学模型的任意阶虚拟元方法:收敛性、稳定性和色散损耗分析,国际。J.数字。方法工程师,122934-971(2021)
[41] Benvenuti,E。;Chiozzi,A。;Manzini,G。;Sukumar,N.,具有奇点和不连续性的拉普拉斯问题的扩展虚拟元方法,计算。方法应用。机械。工程,356571-597(2019)·Zbl 1441.74230号
[42] Benvenuti,E。;Chiozzi,A。;Manzini,G。;Sukumar,N.,二维线弹性断裂的扩展虚拟元法,计算。方法应用。机械。工程,390,第114352条pp.(2022)·Zbl 1507.74452号
[43] Nguyen-Thanh,V。;庄,X。;Ngyyen-Xuan,H。;Rabczuk,T。;Wriggers,P.,二维线弹性断裂分析的虚拟单元法,计算。方法应用。机械。工程,340,366-395(2018)·Zbl 1440.74430号
[44] 侯赛因,A。;Aldakheel,F。;Hudobivnik,B。;Wriggers,P。;Guidault,P。;Allix,O.,基于有效虚拟单元法的脆性裂纹扩展计算框架,Comput。方法应用。机械。工程,159,15-32(2019)
[45] A.鲍里奥。;哈蒙,F。;北卡罗来纳州卡斯特莱托。;怀特,J。;Settgast,R.,耦合孔隙力学的混合模拟有限差分和虚拟元公式,计算。方法应用。机械。工程师,383,第113917条,第(2021)页·Zbl 1506.74390号
[46] 贝内代托,M。;A.鲍里奥。;Kyburg,F。;莫利卡,J。;Scialó,S.,耦合多维流动问题的任意阶混合虚拟元公式,计算。方法应用。机械。工程,391,第114204条pp.(2022)·兹比尔1507.76095
[47] 王,G。;Mu,L。;Wang,Y。;He,Y.,斯托克斯问题的压力-破裂虚拟元方法,计算。方法应用。机械。工程,382,第113879条pp.(2021)·Zbl 1506.76102号
[48] da Veiga,L.B。;罗瓦迪纳,C。;Vacca,G.,多边形网格上Navier-Stokes问题的虚拟元素,SIAM J.Numer。分析。,56, 3, 1210-1242 (2018) ·Zbl 1397.65302号
[49] Chernov,A。;马卡蒂,C。;Mascotto,L.,P-和hp-斯托克斯问题的虚拟元素,高级计算。数学。,47, 24 (2021) ·Zbl 1475.65187号
[50] Wang,Y。;王,G。;Wang,F.,不可压缩流动的自适应虚拟元方法,计算。数学。申请。,101, 63-73 (2021) ·Zbl 1491.76057号
[51] 路易斯安那州贝朗·达维加。;达西,F。;Di Pietro,D.A。;Droniou,J.,多面体网格上Stokes问题的任意阶压力鲁棒DDR和VEM方法,计算。方法应用。机械。工程,397,第115061条pp.(2022)·Zbl 1507.76093号
[52] Antonietti,P。;瓦卡,G。;Verani,M.,Navier-Stokes方程与热方程耦合的虚拟元方法,IMA J.Numer。分析。(2022)
[53] Boscheri,W。;Dimarco,G.等人。;Tavelli,M.,可压缩Navier-Stokes方程的高效二阶全马赫有限体积解算器,计算。方法应用。机械。工程,374,第113602条pp.(2021)·Zbl 1506.76109号
[54] Boscheri,W。;塔维利,M。;Paoluzzi,N.,具有隐式粘性的不可压缩Navier-Stokes方程的高阶有限差分/间断Galerkin格式,Commun。申请。Ind.数学。,13, 21-38 (2022) ·Zbl 1492.76077号
[55] 托罗,E。;卡斯特罗,C。;Vanzo,D。;Siviglia,A.,将浅水方程的通量矢量分裂格式扩展到非结构网格上的高阶,Internat。J.数字。液体方法,1,1-27(2022)
[56] Vater,S。;Klein,R.,低弗劳德数下浅水流动的半隐式多尺度格式,Commun。申请。数学。计算。科学。,13, 303-336 (2018) ·Zbl 1402.65094号
[57] 托罗,E。;Vázquez-Cendón,M.,欧拉方程的通量分裂方案,计算与流体,70,1-12(2012)·Zbl 1365.76243号
[58] Dumbser,M。;Casulli,V.,带一般状态方程的可压缩Navier-Stokes方程的保守、弱非线性半隐式有限体积格式,Appl。数学。计算。,272, 479-497 (2016) ·Zbl 1410.76220号
[59] 奥舍,S。;Solomon,F.,《一种适用于大型刚性Ode系统的部分隐式方法,只有几个方程引入了小时间常数,SIAM J.Numer。分析。,13, 645-663 (1976) ·Zbl 0399.65043号
[60] 吉拉德,H。;Murrone,A.,关于逆风方案在低马赫数限制下的行为:II。Godunov类型方案、计算和流体,33,655-675(2004)·Zbl 1049.76040号
[61] 吉拉德,H。;Viozat,C.,关于低马赫极限下迎风格式的行为,计算与流体,28,63-86(1999)·兹伯利0963.76062
[62] Dellacherie,S.,低马赫数下可压缩欧拉系统的Godunov型格式分析,J.Compute。物理。,229, 978-1016 (2010) ·Zbl 1329.76228号
[63] Stroud,A.,《多重积分的近似计算》(1971年),普伦蒂斯·霍尔公司:普伦蒂斯霍尔公司,新泽西州恩格尔伍德克利夫斯·Zbl 0379.65013号
[64] Sommariva,A。;Vianello,M.,《任意几何体上的高斯-格林体积和力矩计算》,J.Compute。申请。数学。,231, 886-896 (2009) ·Zbl 1170.65017号
[65] 马丁内斯。;Piazzon,F。;Sommariva,A。;Vianello,M.,基于四次多项式优化,Optim。莱特。,14, 1027-1036 (2020) ·Zbl 1445.90076号
[66] 塔维利,M。;Dumbser,M.,二维不可压缩Navier-Stokes方程的交错半隐式间断Galerkin方法,应用。数学。计算。,248, 70-92 (2014) ·兹比尔1338.76068
[67] 萨阿德,Y。;Schultz,M.,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7, 856-869 (1986) ·Zbl 0599.65018号
[68] LeVeque,R.,《高分辨率Godunov方法中平衡源项和通量梯度:准静态波传播算法》,J.Compute。物理。,146, 346-365 (1998) ·Zbl 0931.76059号
[69] Schlichting,H。;Gersten,K.、Grenzschichtheorie(2005)、Springer Verlag
[70] 利维,D。;Puppo,G。;Russo,G.,双曲守恒律系统的中央WENO格式,M2AN数学。模型。数字。分析。,33, 3, 547-571 (1999) ·Zbl 0938.65110号
[71] 利维,D。;Puppo,G。;Russo,G.,《多维守恒定律的紧凑中央WENO方案》,SIAM J.Sci。计算。,22, 2, 656-672 (2000) ·Zbl 0967.65089号
[72] 博斯切里,W。;Dimarco,G.,《一般多边形网格上BGK模型的高阶中心WENO-Implicit-Explicit Runge-Kutta格式》,J.Compute。物理。,422,第109766条pp.(2020)·Zbl 07508389号
[73] Boscheri,W。;Dimarco,G.,非结构化网格上Boltzmann模型的带IMEX时间步进的高阶有限体积格式,计算。方法应用。机械。工程,387,第114180条pp.(2021)·兹比尔1507.65147
[74] 路易斯安那州贝朗·达维加。;达西,F。;Russo,A.,多面体网格上的高阶虚拟元方法,计算。方法应用。机械。工程,741110-1122。(2017) ·Zbl 1448.65215号
[75] 马斯科托,L.,《虚元法中的病态:稳定性和基》,数值。方法偏微分方程,34,4,1258-1281(2018)·兹比尔1407.65294
[76] Canestrelli,A。;Dumbser,M。;Siviglia,A。;Toro,E.,《固定床和移动床浅水方程非结构网格上的井平衡高阶中心格式》,Adv.water Resour。,33, 3, 291-303 (2010)
[77] Toro,E.,Riemann解算器和流体动力学数值方法:实用简介(2009),Springer·Zbl 1227.76006号
[78] Toro,E.,《自由表面浅流的冲击捕获方法》(2001),John Wiley&Sons·Zbl 0996.76003号
[79] 伯内蒂,R。;Titarev,V。;Toro,E.,具有间断底部几何的浅水方程Riemann问题的精确解,J.Compute。物理。,227, 3212-3243 (2008) ·Zbl 1132.76027号
[80] Bassi,F。;Rebay,S.,二维欧拉方程的高精度间断有限元解,J.Compute。物理。,138, 251-285 (1997) ·Zbl 0902.76056号
[81] Womersley,J.R.,《已知压力梯度时动脉中流速、流速和粘性阻力的计算方法》,J.Physiol。,127, 3, 553 (1955)
[82] Loudon,C。;Tordesillas,A.,《使用无量纲Womersley数表征内部流动的非定常性》,J.Theoret。生物学,191,63-78(1998)
[83] 贝尔,J.B。;科莱拉,P。;Glaz,H.M.,《不可压缩Navier-Stokes方程的二阶投影法》,J.Compute。物理。,85, 2, 257-283 (1989) ·Zbl 0681.76030号
[84] Dumbser,M。;佩什科夫,I。;Romenski,E。;Zanotti,O.,连续介质力学统一一阶双曲线公式的高阶ADER格式:粘性导热流体和弹性固体,J.Compute。物理。,314, 824-862 (2016) ·Zbl 1349.76324号
[85] Boscheri,W。;Dumbser,M。;Ioriati,M。;佩什科夫,I。;Romenski,E.,《连续体力学的结构保护交错半隐式有限体积格式》,J.Compute。物理。,424,第109866条pp.(2021)·兹伯利07508471
[86] Ghia,美国。;Ghia,K.N。;Shin,C.T.,《使用Navier-Stokes方程和多重网格法求解不可压缩流的高雷诺数解》,J.Compute。物理。,48, 387-411 (1982) ·Zbl 0511.76031号
[87] 塔维利,M。;Dumbser,M.,二维三角网格上不可压Navier-Stokes方程的交错时空间断Galerkin方法,计算与流体,119,235-249(2015)·Zbl 1390.76360号
[88] Qu,L。;诺伯格,C。;戴维森,L。;彭,S.-H。;Wang,F.,雷诺数在50和200之间时圆柱绕流的定量数值分析,《流体结构杂志》。,39, 347-370 (2013)
[89] Pareschi,L。;Russo,G.,具有刚性松弛的双曲型系统的高阶渐近强稳定保持方法,(Hou,T.;Tadmor,E.,双曲型问题:理论,数值,应用(2003),Springer:Springer-Belin),241-251·Zbl 1064.65105号
[90] 博斯卡里诺,S。;邱,J。;Russo,G。;Xiong,T.,全马赫全欧拉气体动力学系统的高阶半隐式WENO格式,SIAM J.Sci。计算。,44,B368-B394(2022)·兹伯利07511035
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。