×

非定常流高阶模拟中时间误差的影响。 (英语) Zbl 1453.76071号

摘要:本文重点研究时间误差对用高空间阶有限元流线迎风/Petrov-Galerkin解算器建模的非定常解的影响。讨论了六种L稳定的一阶到四阶时间积分方法,并在对流等熵涡和横流中的二维圆柱体这两个典型测试案例中进行了实践。所有结果均显示在由P2(三阶)元素组成的网格上。等熵涡案例用于验证精度的时间顺序,评估每种方法在不同时间步长下耗散误差和色散误差之间的相对分裂,并量化每种方法的相对计算成本。将二维圆柱在不同雷诺数下的横流模拟结果与实验结果和公认的涡脱落频率(斯特鲁哈尔数)、平均阻力系数和均方根升力和阻力系数的数值结果进行了比较。显示了这些参数相对于时间步长的收敛性。这些结果突出了高阶(大于二阶)时间积分方法的优点,以及对非定常模拟进行时间步长收敛研究的必要性。

MSC公司:

76M10个 有限元方法在流体力学问题中的应用
76号06 可压缩Navier-Stokes方程

软件:

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

参考文献:

[1] 王,Z。;Fidkowski,K。;Abgrall,R。;Bassi,F。;Caraeni,D。;A.卡里。;Deconick,H。;哈特曼,R。;Hillewaert,K。;Huynh,H。;北卡罗尔。;May,G。;佩尔松,P.-O。;van Leer,B。;Visbal,M.,《高阶CFD方法:现状和前景》,国际期刊数值。《液体方法》,72,8,811-845(2013)·Zbl 1455.76007号
[2] Cockburn,B。;Hou,S。;Shu,C.-W.,守恒定律的Runge-Kutta局部投影间断Galerkin有限元方法。四、 多维案例,数学。计算。,54190545(1990年)·Zbl 0695.65066号
[3] Bassi,F。;Rebay,S.,二维欧拉方程的高精度间断有限元解,J.Compute。物理。,138, 2, 251-285 (1997) ·Zbl 0902.76056号
[4] Bassi,F。;Rebay,S.,数值求解可压缩Navier-Stokes方程的高精度间断有限元方法,J.Compute。物理。,131, 2, 267-279 (1997) ·Zbl 0871.76040号
[5] Diosady,L.T.(迪奥萨迪,L.T.)。;Darmofal,D.L.,Navier-Stokes方程间断Galerkin解的预处理方法,J.Compute。物理。,228, 11, 3917-3935 (2009) ·Zbl 1185.76812号
[6] Mavrilis,D。;纳斯塔斯,C。;王,L。;沙巴齐,K。;Burgess,N.,《航空航天应用高阶非连续Galerkin方法的进展》,(第47届AIAA航空航天科学会议,包括新视野论坛和航空航天博览会(2009年),美国航空航天研究所:美国航空航天研究所,弗吉尼亚州雷斯顿)
[7] 里德·W·H。;Hill,T.R.,《中子输运方程的三角网格法》(1973),洛斯阿拉莫斯报告LA-UR-73-479(836)
[8] Bassi,F。;克里韦里尼,A。;Rebay,S。;Savini,M.,雷诺平均Navier-Stokes和k-ω湍流模型方程的间断Galerkin解,计算。液体,34,4-5规范发行。,507-540 (2005) ·Zbl 1138.76043号
[9] Remacle,J.-F。;Flaherty,J.E。;Shephard,M.S.,《应用正交基处理可压缩流动问题的自适应间断Galerkin技术》,SIAM Rev.,45,1,53-72(2003)·Zbl 1127.65323号
[10] 北切沃荣。;Xin,J。;胡,P。;李,X。;克莱尔,D。;Flaherty,J.E。;Shephard,M.S.,《应用于冲击和爆炸问题的间断Galerkin方法》,科学杂志。计算。,22-23、227-243(2005年6月)·Zbl 1065.76134号
[11] 布鲁克斯,A.N。;Hughes,T.J.,对流主导流的Streamline迎风/Petrov-Galerkin公式,特别强调不可压缩Navier-Stokes方程,计算。方法应用。机械。工程师,32,1-3,199-259(1982)·Zbl 0497.76041号
[12] Kröner,D。;Ohlberger,M。;Rohde,C.,《保护法理论和数值最新发展导论》(1999年),(alk.论文)·Zbl 0904.00045号
[13] Bonhaus,D.L.,粘性可压缩流动的高精度有限元方法(1998),弗吉尼亚理工学院和州立大学,博士论文
[14] Franca,L.P。;Frey,S.L。;Hughes,T.J.,稳定有限元方法:I.对流扩散模型的应用,计算。方法应用。机械。工程师,95,2,253-276(1992)·Zbl 0759.76040号
[15] Hughes,T.J.,SUPG方法发展和理解的最新进展,特别是可压缩Euler和Navier-Stokes方程,国际期刊Numer。方法流体,7,11,1261-1275(1987)·Zbl 0638.76080号
[16] 柯克,B.S。;Carey,G.F.,使用改进的无粘通量离散化为可压缩Navier-Stokes方程开发和验证SUPG有限元格式,国际J数值。《液体方法》,57,3,265-293(2008)·兹比尔1144.76032
[17] 沙基布,F。;休斯·T·J。;Johan,Z.,计算流体动力学的新有限元公式:X.可压缩Euler和Navier-Stokes方程,计算。方法应用。机械。工程,89,1-3,141-219(1991)
[18] Venkatakrishnan,V。;奥尔马拉斯,S。;卡梅内茨基,D。;Johnson,F.,可压缩Navier-Stokes方程的高阶格式,(第16届AIAA计算流体动力学会议(2003年),美国航空航天研究所:美国航空航天学会,弗吉尼亚州雷斯顿)
[19] 纽曼,J.C。;Anderson,W.K.,《非定常流动和运动区域的非结构化高阶方法研究》,(第22届AIAA计算流体动力学会议(2015年),美国航空航天研究所:美国航空航天学会,弗吉尼亚州雷斯顿)
[20] 尼格罗,A。;Ghidoni,A。;Rebay,S。;Bassi,F.,非定常可压缩流动间断Galerkin解的修正扩展BDF格式,国际J数值。《液体方法》,76,9,549-574(2014)
[21] 霍尔斯特,K.R。;格拉斯比,R.S。;Erwin,J.T。;Stefanski,D.L。;Bond,R.B.,使用HPCMP CREATE(TM)-AV Kestrel组件COFFE的高阶时间精确模拟,(2018年美国航空航天学会航空航天科学会议(2018年1月),美国航空航天学会:美国航空航天学会,弗吉尼亚州莱斯顿)
[22] 西川,H。;Liu,Y.,非定常问题的基于三阶边的方案,(2018年流体动力学会议(2018年6月),美国航空航天学会:美国航空航天学会,弗吉尼亚州莱斯顿)
[23] Bijl,H。;卡彭特,M。;Vatsa,V.,非定常Navier-Stokes方程的时间积分方案,(第15届AIAA计算流体动力学会议(2001年6月),美国航空航天研究所:美国航空航天学会,弗吉尼亚州雷斯顿)
[24] Bijl,H。;Carpenter,M.H。;Vatsa,V.N。;Kennedy,C.A.,《非定常可压缩Navier-Stokes方程的隐式时间积分格式:层流》,J.Compute。物理。,179, 1, 313-329 (2002) ·Zbl 1060.76079号
[25] Carpenter,M.H。;肯尼迪,C.A。;Bijl,H。;Viken,S.A。;Vatsa,V.N.,流体力学应用的四阶Runge-Kutta格式,科学杂志。计算。,25, 1, 157-194 (2005) ·Zbl 1203.76112号
[26] Bijl,H.,使用隐式Runge-Kutta积分方案进行非定常流计算的迭代方法,(第44届AIAA航空航天科学会议和展览(2006年1月),美国航空航天研究所:美国航空航天学会,弗吉尼亚州雷斯顿)
[27] 布隆,D.S。;伯肯,P。;Bijl,H。;Kessels,F。;梅斯特,A。;van Zuijlen,A.H.,《Rosenbrock和ESDIRK方法与非定常可压缩流动迭代求解器的比较》,高级计算。数学。,42, 6, 1401-1426 (2016) ·Zbl 1388.76166号
[28] 佩尔松,P.-O。;Peraire,J.,用于时间相关问题的高效低内存隐式DG算法,(第44届AIAA航空航天科学会议和展览(2006年1月)),1-11
[29] Persson,P.-O.,使用隐式显式Runge-Kutta方案的高阶LES模拟,(第49届AIAA航空航天科学会议,包括新视野论坛和航空航天博览会(2011年),美国航空航天学会:美国航空航天研究所,弗吉尼亚州雷斯顿)
[30] Pazner,W。;Persson,P.O.,《用于不连续伽辽金流体模拟的阶段并行全隐式龙格-库塔解算器》,J.Comput。物理。,335, 700-717 (2017) ·Zbl 1375.76164号
[31] Pazner,W。;Persson,P.-O.,使用隐式张量积非连续Galerkin方法的高阶DNS和LES模拟,(第23届AIAA计算流体动力学会议(2017年6月),美国航空航天学会:美国航空航天研究所,弗吉尼亚州雷斯顿),1-17
[32] Nigro,A。;De Bartolo,C。;伦达,S.M。;Bassi,F.,应用于间断Galerkin离散Euler方程的高阶时间积分方法的比较,能量程序。,45, 518-527 (2014)
[33] Cash,J.R.,ODE和DAE中刚性初值问题数值解的修正后向微分公式,J.Compute。申请。数学。,125、1-2、117-130(2000年)·兹伯利0971.65063
[34] Psihoyios,G.,一组隐式高级步长(IAS)方法稳定性函数的通用公式,数学。计算。型号。,46, 1-2, 214-224 (2007) ·Zbl 1132.65074号
[35] 萨阿德,Y。;Schultz,M.H.,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7, 3, 856-869 (1986) ·Zbl 0599.65018号
[36] 格拉斯比,R.S。;Erwin,J.T.,《COFFE简介:下一代HPCMP CREATE(TM)-AV CFD求解器》,(第54届AIAA航空航天科学会议(2016年),美国航空航天学会:美国航空航天研究所,弗吉尼亚州雷斯顿)
[37] 哈利西,B.P。;格拉斯比,R.S。;Erwin,J.T.,《HPCMP CREATE(TM)-AV COFFE和Kestrel在二维和三维湍流情况下的比较》,(第54届AIAA航空航天科学会议(2016),美国航空航天研究所:美国航空航天学会,弗吉尼亚州雷斯顿)
[38] Erwin,J.T。;格拉斯比,R.S。;Stefanski,D.L。;Nichols,R.H.,《HPCMP CREATE(TM)-AV COFFE在三维湍流情况下的应用》,(第54届AIAA航空航天科学会议(2016年),美国航空航天研究所:美国航空航天学会,弗吉尼亚州雷斯顿)
[39] Dahlquist,G.G.,线性多步方法的一个特殊稳定性问题,BIT-Numer。数学。,3, 1, 27-43 (1963) ·Zbl 0123.11703号
[40] 苏丽,E。;Mayers,D.F.,《数值分析导论》(2003),剑桥大学出版社:剑桥大学出版社·Zbl 1033.65001号
[41] 库兰特,R。;弗里德里希(Friedrichs,K.)。;Lewy,H.,《数学差异》。《年鉴》,100,1,32-74(1928)·JFM 54.0486.01型
[42] 海尔,E。;Wanner,G.,《求解常微分方程II》,《计算数学中的Springer级数》,第14卷(1996),Springer Berlin Heidelberg:Springer Barlin Heitelberg Berlin,Heidelbrg·Zbl 0859.65067号
[43] 肯尼迪,C.A。;Carpenter,M.H.,《常微分方程的对角隐式Runge-Kutta方法》。《综述》(2016年3月),技术代表,NASA TM-2016-219173
[44] Alexander,R.,《刚性O.D.E.的对角隐式Runge-Kutta方法》,SIAM J.Numer。分析。,1006-1021年6月14日(1977年)·Zbl 0374.65038号
[45] Alexander,R.K.,刚性常微分方程Runge-Kutta方法的稳定性,SIAM J.Numer。分析。,31, 4, 1147-1168 (1994) ·Zbl 0811.65061号
[46] Alexander,R.,《刚性系统DIRK积分器的设计与实现》,应用。数字。数学。,46, 1, 1-17 (2003) ·Zbl 1023.65072号
[47] A.Kvrnö,More,and to be hope,better DIRK methods for solution of stiff ODE,未出版注释。
[48] Kværnø,A.,具有显式第一阶段的单对角隐式Runge-Kutta方法,BIT-Number。数学。,44, 3, 489-502 (2004) ·兹比尔1066.65077
[49] Verwer,J.G.,关于BN稳定性概念的实用价值,BIT-Numer。数学。,21, 3, 355-361 (1981) ·Zbl 0476.65055号
[50] 柯蒂斯,C.F。;Hirschfelder,J.O.,刚性方程的积分,Proc。国家。阿卡德。科学。,38, 3, 235-243 (1952) ·Zbl 0046.13602号
[51] Gear,C.W.,常微分方程中的数值初值问题,自动计算中的Prentice-Hall级数(1971),Prentice-Hall·Zbl 1145.65316号
[52] 瓦萨,V。;卡彭特,M。;Lockard,D.,非定常流应用优化二阶后向差分(BDF2OPT)方案的重新评估,(第48届AIAA航空航天科学会议,包括新视野论坛和航空航天博览会(2010年),美国航空航天学会:美国航空航天研究所,弗吉尼亚州雷斯顿)
[53] Chen,D.J.L.,刚性微分方程的单隐式Runge-Kutta方法的效率,数值。算法,65,3,533-554(2014)·Zbl 1291.65199号
[54] 海尔,E。;Wanner,G.,《求解常微分方程I》,《计算数学中的Springer级数》,第8卷(1993),Springer Berlin Heidelberg:Springer Barlin Heitelberg Berlin,Heidelbrg·Zbl 0789.65048号
[55] 梅辛格,F。;Arakawa,A.,《大气模型中使用的数值方法》,第一卷,GARP Publ。序列号。,第17卷(1976年),气象组织-国际科学理事会
[56] Kultler,P。;Warming,R。;Lomax,H.,使用非中心有限差分格式计算航天飞机流场,AIAA J.,11,2,196-204(1973)·Zbl 0268.76050号
[57] Bouche,D。;邦诺,G。;Ramos,D.,求解对流方程的数值格式比较,应用。数学。莱特。,16, 2, 147-154 (2003) ·Zbl 1047.65064号
[58] Wicker,L.J。;Skamarock,W.C.,使用正向时间方案的弹性模型的时间分割方法,周一。《天气评论》,130、8、2088-2097(2002)
[59] Takacs,L.L.,耗散和弥散误差最小的平流方程的两步格式,Mon。《天气评论》,113、6、1050-1065(1985)
[60] Yee,H.C。;北卡罗来纳州桑德姆。;Djomehri,M.J.,使用基于特征的滤波器的低密度高阶冲击捕获方法,J.Compute。物理。,150, 1, 199-238 (1999) ·Zbl 0936.76060号
[61] 王,L。;Mavrilis,D.J.,高精度间断Galerkin离散的非定常Euler方程的隐式解,J.Compute。物理。,225, 2, 1994-2015 (2007) ·Zbl 1343.76022号
[62] 明镜,南卡罗来纳州。;Huynh,H。;DeBonis,J.R.,《使用高阶方法对等熵Euler涡问题的调查》,(第22届AIAA计算流体动力学会议(2015年),美国航空航天研究所:美国航空航天学会,弗吉尼亚州雷斯顿)
[63] Zdravkovich,M.M.,《圆柱绕流:流动现象、实验、应用、数学模型和计算机模拟的综合指南》(1997),牛津大学出版社:牛津大学出版社
[64] 米塔尔·R。;Balachandar,S.,《三维性对名义上二维圆柱体升力和阻力的影响》,Phys。流体,7,8,1841-1865(1995)·Zbl 1032.76530号
[65] Henderson,R.D.,湍流尾流过渡中的非线性动力学和模式形成,J.流体力学。,352, 65-112 (1997) ·Zbl 0903.76070号
[66] 拉贾尼,B.N。;Kandasamy,A。;Majumdar,S.,通过圆柱的层流数值模拟,应用。数学。型号。,331228-1247(2009年)·Zbl 1168.76305号
[67] Williamson,C.H.K.,《尾迹过渡过程中斑点状“涡位错”的自然和强制形成》,J.流体力学。,243, 393-441 (1992)
[68] 威廉姆森,C.H.K.,《圆柱尾迹中的旋涡动力学》,年。流体力学版次。,28, 1, 477-539 (1996)
[69] Mittalt,R.等人。;Balachandar,S.,《关于在二维钝体尾流模拟中纳入三维效应》,(1997年ASME流体工程部夏季会议FEDSM’97(1997)),1-9
[70] Norberg,C.,《圆柱绕流:脉动升力方面》,J.流体结构。,15,3-4,459-469(2001年)
[71] Wieselsberger,C.,流体阻力定律的新数据(1922),技术代表,NACA-TN-84
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。