×

非定常不可压Navier-Stokes方程的高效高阶半隐式时间步长格式。 (英语) Zbl 1410.76186号

总结:我们重点研究了求解不可压缩Navier-Stokes方程的高阶半隐式时间步长方法的发展。泰勒-胡德混合有限元法{P} _2-\mathbb{P} 1个\))用于速度和压力的空间近似。半隐式时间离散可以为非平稳流提供非常精确的近似,同时显式地处理非线性项,避免了对复杂非线性迭代求解器的需要。我们首先介绍了三种高阶半隐式时间步长格式:最容易实现的多步半隐式后向差分公式(SBDF);一步缺陷校正(DC),在给定阶次的速度和压力(时间)的数值误差方面,它比SBDF方法产生更好的近似值;由提出的J.-L.格蒙德P.Minev公司(GM)[SIAM J.Sci.Compute.37,No.6,A2656–A2681(2015;Zbl 1382.65315号)]一种基于缺陷修正和人工压缩方法的解耦方法。SBDF和DC方法生成鞍点线性系统,而GM生成两个具有对称正定矩阵的线性系统。我们提出了对GM表示为GM-SRM方法的修改,该方法实现了Navier-Stokes方程的顺序正则化方法。与GM方法相比,GM-SRM方法在启动期间对压力的收敛速度更快,同时需要更小的稳定参数。我们通过数值测试案例展示了这些方法的数值准确性、稳定性和效率:两个制造的解决方案;流经圆柱体的流量;和眼睑驱动的腔流,都在2D设置中。所有方法均与理论收敛速度一致。SBDF方法是CPU效率最高的方法,其次是DC、GM和GM-SRM方法。我们观察到,GM和GM-SRM中的(text{grad}-\text{div})项的存在提高了数值稳定性,从而产生了更高的CFL界。此外,相关参数,如斯特鲁哈尔数、升力系数和阻力系数与公布的值非常接近。

MSC公司:

76M10个 有限元方法在流体力学问题中的应用
65平方米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
76D05型 不可压缩粘性流体的Navier-Stokes方程
35季度30 Navier-Stokes方程
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Akrivis,G。;克鲁泽克斯,M。;Makridakis,C.,非线性抛物问题的隐式显式多步有限元方法,数学计算,67,222,457-477,(1998)·Zbl 0896.65066号
[2] Auteri,F。;北卡罗来纳州帕罗里尼。;Quartapelle,L.,奇异驱动空腔流稳定性的数值研究,《计算物理杂志》,183,1-25,(2002)·兹比尔1021.76040
[3] 阿舍尔,U.M。;Ruuth,S.J。;Wetton,B.T.R.,含时偏微分方程的隐式显式方法,SIAM J Numer Anal,32,3,797-823,(1995)·Zbl 0841.65081号
[4] Blevins,R.D.,《流致振动》,377,(1990),Van Nostrand Reinhold Co.,Inc.,美国纽约州纽约市
[5] Boffi,D。;布雷齐,F。;Fortin,M.,《混合有限元方法和应用》,691,(2013),柏林施普林格出版社,海德堡·Zbl 1277.65092号
[6] Baek,H。;Karniadakis,G.E.,Sub-iteration导致Navier-Stokes方程半隐式格式的精度和稳定性增强,计算物理杂志,230,4384-4402,(2011)·Zbl 1416.76093号
[7] 贝克,G.A。;Dougalis,V.A。;Karakashian,O.A.,关于Navier-Stokes方程的高阶精确全离散Galerkin近似,数学计算,39,160,339-375,(1982)·Zbl 0503.76038号
[8] 文莱,C.H。;Saad,M.,《重新审视二维眼睑驱动腔问题》,《计算流体》,35,326-348,(2005)·Zbl 1099.76043号
[9] Chorin,A.J.,《解决不可压缩粘性流动问题的数值方法》,《计算物理杂志》,2,1,12-26,(1967)·Zbl 0149.44802号
[10] 戴,X。;Sun,J。;Cheng,X.L.,Navier-Stokes方程算子分裂方法的误差估计:二阶格式,计算应用数学杂志,231,696-704,(2009)·兹比尔1167.76024
[11] Ethier,M。;Bourgault,Y.,bidomain模型的半隐式时间离散化方案,SIAM J Numer Ana,46,5,2443-2468,(2008)·Zbl 1182.92009年
[12] Elman,H。;西尔维斯特,D。;Wathen,A.,有限元和快速迭代求解器在不可压缩流体动力学中的应用,400,(2008),牛津大学出版社,纽约
[13] Fortin A.、Jardak M.、Gervais J.J.、Pierre R.,1997年。流体流动问题中Hopf分岔的局部化。24.第1185-1210页。;Fortin A.、Jardak M.、Gervais J.J.、Pierre R.,1997年。流体流动问题中Hopf分岔的局部化。24.第1185-1210页·Zbl 0886.76042号
[14] Fortin,M。;Glowinski,R.,增广拉格朗日方法:边界值问题数值解的应用,340,(2000),Elsevier
[15] 戈尔茨坦S.,1938年。流体动力学的现代发展。35.第421-447页。;戈尔茨坦S.,1938年。流体动力学的现代发展。35.第421-447页。
[16] Girault,V。;Raviart,P.A.,《Navier-Stokes方程的有限元方法:理论和算法》,376,(2011),施普林格-弗拉格-海德堡
[17] 格洛温斯基,R。;Pironneau,O.,Navier-Stokes方程的有限元方法,流体力学年鉴,24167-204,(1992)·Zbl 0743.76051号
[18] Ghia,美国。;Ghia,K.N。;Shin,C.T.,使用Navier-Stokes方程和多重网格方法求解不可压缩流的高分辨率解,《计算物理杂志》,48,387-411,(1982)·Zbl 0511.76031号
[19] Gresho,P.M。;格里菲斯,D.F。;Silvester,D.J.,《不可压缩流的自适应时间步进第i部分:标量平流扩散》,SIAM科学计算杂志,30,42018-2054,(2008)·Zbl 1167.76026号
[20] Guermond,J.L。;Minev,P.D.,不可压缩Navier-Stokes方程的高阶时间步进,SIAM科学计算杂志,37,6,A2656-A2681,(2015)·Zbl 1382.65315号
[21] Guermond,J.L。;Minev,P。;Jie,S.,《不可压缩流投影方法概述》,计算方法应用机械工程,195,6011-6045,(2006)·Zbl 1122.76072号
[22] Hecht F..FreeFEM++:第三版,3.32版。2014http://www.freefem.org/ff++/ftp/freefem++doc.pdf;Hecht F..FreeFEM++:第三版,3.32版。2014http://www.freefem.org/ff++/ftp/freefem++文档.pdf
[23] John,V.,圆柱体周围二维时间相关流的阻力和升力参考值,《国际数值方法流体》,44,777-788,(2004)·Zbl 1085.76510号
[24] Kupferman,R.,不可压缩粘性流纯流函数公式的中心差分格式,SIAM科学计算杂志,23,1,1-18,(2001)·Zbl 0995.76060号
[25] 克劳切克,P。;Rys,F.S.,分数阶跃的稳定性θ-非平稳Navier-Stokes方程的格式,SIAM J Numer Ana,31,5,1312-1335,(1994)·Zbl 0810.76036号
[26] Kress,W。;Gustafsson,B.,初始边值问题的延迟校正方法,科学计算杂志,17,1-4241-251,(2001)·Zbl 1003.65106号
[27] Kress,W。;Lötstedt,P.,使用半隐式方法求解不可压缩Navier-Stokes方程的时间步长限制,计算方法应用机械工程,195,4433-4447,(2006)·Zbl 1124.76040号
[28] Kim,J。;Moin,P.,分步法在不可压缩Navier-Stokes方程中的应用,计算物理杂志,59,308-323,(1985)·Zbl 0582.76038号
[29] Labovschii,A.,含时Navier-Stokes方程的缺陷修正方法,数值方法部分微分Equ,25,1,1-25,(2009)·Zbl 1394.76087号
[30] Lin,P.,含时不可压缩Navier-Stokes方程的序贯正则化方法,SIAM J Numer Ana,34,3,1051-1071,(1997)·Zbl 0882.76018号
[31] 拉博夫斯基,A。;莱顿,W.J。;马尼卡,哥伦比亚特区。;内达,M。;Rebholz,L.G.,Navier-Stokes方程的稳定外推梯形有限元方法,计算方法应用机械工程,198,958-974,(2009)·Zbl 1229.76051号
[32] 卢,X.L。;林,P。;Liu,J.G.,非定常Navier-Stokes方程的序列正则化方法分析,数学计算,77,263,1467-1494,(2008)·Zbl 1285.76018号
[33] 米塔尔,S。;Raghuvanshi,A.,《低雷诺数流动中圆柱后涡旋脱落的控制》,《内部J数值方法流体》,35,421-447,(2001)·兹比尔1013.76049
[34] Mittal S.,Tezduyar T.E.,1992年。关于非定常不可压流的稳定时空有限元公式的注记。73.第93-112页。;Mittal S.,Tezduyar T.E.,1992年。关于非定常不可压流的稳定时空有限元公式的注记。73.第93-112页。
[35] Olshanskii,硕士。;润滑油,G。;Heister,T。;Löwe,J.,不可压缩Navier-Stokes方程的梯度-径向稳定和亚网格压力模型,计算方法应用-机械工程,198,3975-3988,(2009)·Zbl 1231.76161号
[36] Olshanskii,硕士。;Reusken,A.,《斯托克斯方程的梯度-径向稳定化》,数学计算,73,248,1699-1718,(2003)·Zbl 1051.65103号
[37] 潘·T·W。;Glowinski,R.,用投影/波动方程方法对Navier-Stokes方程模拟的不可压缩粘性流体流动进行数值模拟,计算流体动力学J,9,28-42,(2000)
[38] Quarteroni,A。;萨列里,F。;Veneziani,A.,Navier-Stokes方程数值逼近的因式分解方法,计算方法应用机械工程,188,505-526,(1998)·兹比尔0976.76044
[39] 拉贾尼,B.N。;Kandasamy,A。;Majumdar,S.,通过圆柱体的层流数值模拟,应用数学模型,331228-1247,(2009)·Zbl 1168.76305号
[40] Rangan,A.,低指数微分代数方程的延迟校正方法,BIT数值数学,43,1-18,(2003),Kluwer学术出版社
[41] Ruuth,S.J.,模式形成中反应扩散问题的隐式-显式方法,《数学生物学杂志》,34,148-176,(1995)·Zbl 0835.92006号
[42] Schäfer,M。;Turek,S。;杜斯特,F。;克劳斯,E。;Rannacher,R.,圆柱周围层流的基准计算,547-566,(1996),Vieweg+Teubner Verlag·Zbl 0874.76070号
[43] Shankar,P.N。;Deshpande,M.D.,驱动腔中的流体力学,《流体力学年鉴》,32,93-136,(2000)·Zbl 0988.76006号
[44] Simo,J.C.(西莫,J.C.)。;Armero,F.,不可压缩Navier-Stokes和Euler方程瞬态算法的无条件稳定性和长期行为,计算方法应用机械工程,111,111-154,(1993)·Zbl 0846.76075号
[45] Stetter,H.J.,《缺陷补偿原理和离散化方法》,《数值数学》,29,425-443,(1978),施普林格出版社·Zbl 0362.65052号
[46] Temam,R.,Navier-Stokes方程:理论和数值分析,407,(2001),AMS Chelsea Publishing Providence,罗德岛·Zbl 0981.35001号
[47] 蒂辛加,G。;伍布斯,F.B。;Veldman,A.E.P.,用Newton-Picard方法对驱动腔内不可压缩流动进行分叉分析,《计算应用数学杂志》,140,751-772,(2001)·Zbl 1063.76040号
[48] Turek,S.,《不可压缩Navier-Stokes方程时间步长技术的比较研究:从完全隐式非线性方案到半隐式投影方法》,《国际数值方法流体》,22,987-1011,(1996)·兹伯利0864.76052
[49] Turek,S.,《不可压缩流问题的高效求解器:算法和计算方法》,《计算科学与工程讲义》,358,(1999),柏林施普林格出版社·Zbl 0930.76002号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。