×

FLEXI:双曲抛物守恒律的高阶间断Galerkin框架。 (英语) Zbl 1461.76347号

摘要:高阶(HO)格式是流体力学及相关学科中多尺度问题数值解的理想候选格式。在HO离散化变体中,非连续Galerkin格式提供了一系列有利的特征,这些特征在过去十年中引起了人们对其和相关公式的极大兴趣。这些方法已经足够成熟,可用于一系列问题,例如湍流的直接数值模拟和大涡模拟。然而,为了充分利用这些方法的潜在优势,必须在设计和执行仿真链中的所有步骤时考虑HO。特别是在这一领域,许多商用的封闭源代码解决方案都不满足要求。因此,在这项工作中,我们提出FLEXI公司框架,一个HO一致的开源仿真工具链,用于在CPU集群上求解可压缩的Navier-Stokes方程。我们描述了数值算法和实现细节,并概述了框架所有部分的特性和功能。除了这些技术细节,我们还讨论了代码稳定性、再现性和用户友好性等重要但经常被忽视的问题。本文讨论了开发开放源码框架所获得的好处,特别关注开放源码社区的可用性。我们以示例应用程序结束,这些应用程序演示了广泛的用例和FLEXI公司以及当前和未来发展的概述。

MSC公司:

76个M12 有限体积法在流体力学问题中的应用
76M20码 有限差分方法在流体力学问题中的应用
76号06 可压缩Navier-Stokes方程
76升05 流体力学中的冲击波和爆炸波
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Silberman,I.,《大气中的行星波》,J.Meteorol。,11, 1, 27-34 (1954)
[2] Orszag,S.A.,湍流模拟的数值方法,物理学。流体,12,12,II-250(1969)·Zbl 0217.25803号
[3] 达尔伯格,R。;Picone,J.,可压缩磁流体动力学湍流的伪谱模拟,计算。方法应用。机械。工程,80,1-3,409-416(1990)·Zbl 0722.76056号
[4] 横川,M。;Itakura,K。;Uno,A。;石原慎太郎。;Kaneda,Y.,16.4-T在地球模拟器上用傅里叶谱方法直接数值模拟湍流,(SC'02:2002 ACM/IEEE超级计算会议论文集(2002),IEEE),50
[5] 杜姆布瑟。;Käser,M.,线性双曲方程组非结构网格上的任意高阶非振荡有限体积格式,J.Compute。物理。,221, 2, 693-723 (2007) ·Zbl 1110.65077号
[6] Shu,C.-W.,计算流体力学的高阶有限差分和有限体积WENO格式以及不连续伽辽金方法,国际计算杂志。流体动力学。,17, 2, 107-118 (2003) ·Zbl 1034.76044号
[7] 斯瓦德,M。;Carpenter,M.H。;Nordström,J.,可压缩Navier-Stokes方程的稳定高阶有限差分格式,远场边界条件,J.Compute。物理。,225, 1, 1020-1038 (2007) ·Zbl 1118.76047号
[8] Maday,Y。;Patera,A.T.,不可压缩Navier-Stokes方程的谱元方法,(计算力学最新研究(1989)),71-143
[9] 巴西,F。;Rebay,S.,数值求解可压缩Navier-Stokes方程的高精度间断有限元方法,J.Compute。物理。,131, 2, 267-279 (1997) ·Zbl 0871.76040号
[10] Cockburn,B。;Karniadakis,G.E。;Shu,C.-W.,《间断Galerkin方法:理论、计算和应用》,第11卷(2012年),施普林格科学与商业媒体
[11] 坎特韦尔,C.D。;莫西·D·。;Comerford,A。;Bolis,A。;罗科·G。;Mengaldo,G。;德格拉齐亚,D。;雅科夫列夫,S。;伦巴第,J.-E。;Ekelschot,D。;Jordi,B。;Xu,H。;穆罕默德,Y。;埃斯基尔森,C。;B.纳尔逊。;Vos,P。;比奥托,C。;Kirby,R.M。;Sherwin,S.J.,Nektar plus plus:开源光谱/hp元素框架,计算。物理学。Comm.,192,205-219(2015)·Zbl 1380.65465号
[12] H.Klimach,K.Jain,S.Roller,《使用APES进行端到端并行模拟》,载于:PARCO,2013年,第703-711页。
[13] F.D.威瑟顿。;法林顿,A.M。;Vincent,P.E.,PyFR:一个开源框架,用于使用流量重建方法解决流架构上的对流-扩散类型问题,Compute。物理学。通信,185,11,3028-3040(2014)·Zbl 1348.65005号
[14] Hillewaert,K.,《工业几何学中高分辨率、大规模CFD和声学的间断Galerkin方法的发展》(2013),卢浮大学出版社
[15] Gassner,G。;Kopriva,D.A.,《高斯和高斯-洛巴托间断Galerkin谱元方法的色散和耗散误差比较》,SIAM J.Sci。计算。,33, 5, 2560-2579 (2011) ·Zbl 1255.65089号
[16] 莫拉,R.C。;Mengaldo,G。;佩罗,J。;Sherwin,S.J.,《关于欧拉湍流隐式LES/欠分辨率DNS的高阶间断Galerkin方法的涡流分辨率》,J.Compute。物理。,330, 615-623 (2017) ·Zbl 1378.76036号
[17] Drikakis,D。;哈恩,M。;莫塞代尔,A。;Thornber,B.,《使用高分辨率和高阶方法进行大涡模拟》,Phil.Trans。R.Soc.A,36718992985-2997(2009)·Zbl 1185.76719号
[18] 巴西,F。;Rebay,S.,二维欧拉方程的高精度间断有限元解,J.Compute。物理。,138, 2, 251-285 (1997) ·Zbl 0902.76056号
[19] M.J.Brazell,D.J.Mavrilis,NACA 0012翼型的高阶间断Galerkin网格解析湍流模拟,载于:第53届AIAA航空航天科学会议,2015年,第1529页。
[20] A.Uranga,P.-O.Persson,M.Drela,J.Peraire,翼型和机翼过渡流动的隐式大涡模拟,载于:第19届AIAA计算流体动力学,2009年,第4131页·Zbl 1242.76085号
[21] L.T.Diosady,S.M.Murman,使用非连续Galerkin谱元法对周期性山丘上的水流进行DNS,载于:2014年第44届AIAA流体动力学会议,第2784页。
[22] 伦巴第,J.-E.W。;莫西·D·。;Sherwin,S.J。;霍斯勒,J.F。;丹达帕尼,S。;Taylor,M.J.,翼尖涡的隐式大涡模拟,AIAA J.,54,2,506-518(2016)
[23] 博姆,M。;温特斯,A.R。;Gassner,G.J。;Derigs,D。;辛登朗,F。;Saur,J.,电阻MHD方程的熵稳定节点间断Galerkin方法。第一部分:理论和数值验证,J.Compute。物理学。(2018)
[24] Hindenlang,F.,非结构网格高阶并行模拟的网格弯曲技术(2014),斯图加特大学(博士论文)
[25] 辛登朗,F。;Bolemann,T。;Munz,C.-D.,高阶非连续Galerkin模拟的网格弯曲技术,(IDIHOM:高阶方法的工业化——自顶向下方法(2015),Springer),133-152
[26] Geuzaine,C。;Remacle,J.-F.c.,Gmsh:一个内置预处理和后处理设施的三维有限元网格生成器,国际。J.数字。方法工程,79,11,1309-1331(2009)·Zbl 1176.74181号
[27] Touloge,T。;Geuzaine,C。;Remacle,J.-F。;Lambrechts,J.,《曲线网格的稳健解缠》,J.Compute。物理。,254,8-26(2013)·Zbl 1349.65670号
[28] 福图纳托,M。;Persson,P.-O.,使用Winslow方程生成高阶非结构化曲面网格,J.Compute。物理。,307, 1-14 (2016) ·Zbl 1352.65607号
[29] Marcon,J。;科普里瓦,D.A。;Sherwin,S.J。;Peiró,J.,四边形网格生成的高分辨率PDE方法,J.Compute。物理学。(2019) ·Zbl 1453.65045号
[30] HDF Group,J.,分层数据格式,第5版(1997)
[31] Farin,G.,《CAGD的曲线和曲面:实用指南》(2002),摩根考夫曼出版社:摩根考夫文出版社,美国加利福尼亚州旧金山
[32] 法林,G。;Hansford,D.,离散Coons补丁,计算。辅助Geom。设计,16,7,691-700(1999)·兹比尔0997.65033
[33] de Boer,A。;van Zuijlen,A。;Bijl,H.,界面插值和网格变形的径向基函数,(科学与工程高级计算方法(2009),Springer),143-178·Zbl 1187.68668号
[34] Zeifang,J。;Kaiser,K。;贝克,A。;舒茨,J。;Munz,C.-D.,低马赫数流的高效高阶间断Galerkin计算,Commun。申请。数学。计算。科学。,13, 2, 243-270 (2018) ·Zbl 1408.65091号
[35] Vangelatos,S.,关于非定常可压缩Navier-Stokes方程隐式间断Galerkin谱元方法的效率(2019),斯图加特大学,(博士论文)
[36] 阿塔克,M。;拉尔森,J。;Munz,C.-D.,《多核挑战:利用DGSEM对高雷诺数超音速湍流边界层进行空间发展的Petascale DNS》,(Resch,M.M.;Bez,W.;Focht,E.;Kobayashi,H.;Qi,J.;Roller,S.,《2015年持续模拟性能》(2015),Springer International Publishing:Springer国际出版公司Cham), 171-183
[37] Schulz-Rinne,C.W.,《二维气体动力学黎曼问题的分类》,SIAM J.Math。分析。,24, 1, 76-88 (1993) ·Zbl 0811.35082号
[38] Gassner,G.J。;温特斯,A.R。;Kopriva,D.A.,可压缩Euler方程具有逐部分求和性质的分裂形式节点间断Galerkin格式,J.Compute。物理。,327, 39-66 (2016) ·Zbl 1422.65280号
[39] 贝克,A.D。;弗拉德·D·G。;Tonhäuser,C。;加斯纳,G。;Munz,C.-D.,《多项式去混叠对亚网格模型的影响》,Flow Turbul。库布斯特。,97, 2, 475-511 (2016)
[40] 弗拉德·D。;Gassner,G.,《关于在大涡模拟中使用动能保持DG-方案》,J.Compute。物理。,350, 782-795 (2017) ·Zbl 1380.76019号
[41] 贝克,A.D。;Bolemann,T。;弗拉德·D。;弗兰克·H。;Gassner,G.J。;辛登朗,F。;Munz,C.-D.,过渡流和湍流模拟的高阶间断Galerkin谱元方法,国际。J.数字。方法流体,76,8,522-548(2014)
[42] Toro,E.F.,Riemann Solvers and Numerical Methods for Fluid Dynamics:a Practical Introduction(2013),施普林格科学与商业媒体
[43] 巴西,F。;Rebay,S。;Mariotti,G。;Pedinotti,S。;Savini,M.,《无粘和粘性涡轮机械流动的高精度间断有限元法》(Decuypere,R.;Dibelius,G.,第二届欧洲涡轮机械流体动力学和热力学会议(1997)),99-108
[44] Smagorinsky,J.,《原始方程的一般循环实验:I.基础实验》,Mon。《天气评论》,91,3,99-164(1963)
[45] Nicoud,F。;托达·H·B。;Cabrit,O。;Bose,S。;Lee,J.,《使用奇异值构建用于大涡模拟的亚脊尺度模型》,Phys。流体,23,8,第085106条pp.(2011)
[46] Vreman,A.,湍流剪切流的涡粘性亚网格尺度模型:代数理论和应用,Phys。流体,16,10,3670-3681(2004)·Zbl 1187.76543号
[47] 弗拉德·D。;贝克,A。;Munz,C.-D.,使用高阶间断Galerkin谱元方法通过自适应滤波模拟欠分辨率湍流,J.Compute。物理。,313, 1-12 (2016) ·Zbl 1349.76118号
[48] 拉尔森,J。;卡瓦伊,S。;Bodart,J。;Bermejo-Moreno,I.,《模拟墙应力的大涡模拟:最新进展和未来方向》,Mech。工程版次,3,1,15-00418(2016)
[49] P.Spalart,S.Allmaras,《气动流动的单方程湍流模型》,载于:第30届航空航天科学会议和展览,1992年,第439页。
[50] Collis,S.,湍流模拟的间断Galerkin方法,(2002年夏季计划论文集,湍流研究中心),155-167
[51] De Laage De Meux,B。;Audebert,B。;Manceau,R。;Perrin,R.,大涡模拟和RANS/LES混合建模中合成湍流生成的各向异性线性强迫,Phys。流体,27,3,035115(2015)
[52] 普鲁特,C。;Gatski,T。;格罗施,C.E。;Thacker,W.,《时间滤波的Navier-Stokes方程:残余应力的特性》,Phys。流体,15,8,2127-2140(2003)·Zbl 1186.76434号
[53] Ahrens,J。;Geveci,B。;Law,C.,《Paraview:大型数据可视化的最终用户工具》(The visualization Handbook,Vol.717(2005),Elsevier München)
[54] 贝克,A。;Krais,N.,Blender管道(2018年)
[55] Blender Online Community,A.,《Blender-3D建模和渲染包》(2018年),Blender Foundation,Stichting Blender Foundation:Blender基金会,Sticht Blender Foundation Amsterdam
[56] Kopriva,D.A.,曲线网格上的度量恒等式和间断谱元方法,J.Sci。计算。,26301(2006年)·Zbl 1178.76269号
[57] 科普里瓦,D.A。;Gassner,G.,关于配点型间断Galerkin谱元方法中的求积和弱形式选择,J.Sci。计算。,44, 2, 136-155 (2010) ·Zbl 1203.65199号
[58] Kopriva,D.A.,可压缩流的保守交错网格Chebyshev多域方法。二、。半结构化方法,J.Compute。物理。,128, 2, 475-488 (1996) ·Zbl 0866.76064号
[59] 费希尔,T.C。;Carpenter,M.H.,非线性守恒律的高阶熵稳定有限差分格式:有限域,J.Compute。物理。,252, 518-557 (2013) ·Zbl 1349.65293号
[60] Tadmor,E.,守恒定律系统熵稳定格式的数值粘性。I、 数学。公司。,49, 179, 91-103 (1987) ·Zbl 0641.65068号
[61] Jameson,A.,《气体动力学动能守恒格式的制定以及使用熵和动能守恒方案对激波管中一维粘性可压缩流动进行直接数值模拟》,J.Sci。计算。,34, 2, 188-208 (2008) ·Zbl 1133.76031号
[62] 布莱斯德尔,G。;Spyropoulos,E。;秦,J.,谱方法中非线性项公式对混叠误差的影响,应用。数字。数学。,21, 3, 207-219 (1996) ·Zbl 0858.76060号
[63] 肯尼迪,C.A。;Gruber,A.,可压缩流体Navier-Stokes方程中对流项的简化混叠公式,J.Compute。物理。,227, 3, 1676-1700 (2008) ·Zbl 1290.76135号
[64] 冯·诺依曼,J。;Richtmyer,R.D.,流体动力学冲击的数值计算方法,J.Appl。物理。,21, 3, 232-237 (1950) ·Zbl 0037.12002号
[65] P.-O.Persson,J.Peraire,《不连续Galerkin方法的亚细胞冲击捕获》,载于:第44届美国航空航天局航空科学会议和展览,2006年,第112页。
[66] 邱,J。;Shu,C.-W.,Hermite WENO格式及其作为Runge-Kutta间断Galerkin方法限制器的应用:一维情况,J.Compute。物理。,193, 1, 115-135 (2004) ·Zbl 1039.65068号
[67] A.韦尔塔。;卡索尼,E。;Peraire,J.,《高阶间断Galerkin方法的简单冲击捕获技术》,国际。J.数字。《液体方法》,69,10,1614-1632(2012)·Zbl 1253.76058号
[68] 桑塔格,M。;Munz,C.-D.,《使用有限体积子单元的非连续Galerkin方法激波捕获的高效并行化》,J.Sci。计算。,70, 3, 1262-1289 (2017) ·兹比尔1366.65089
[69] 桑塔格,M。;Munz,C.-D.,使用有限体积子单元的非连续Galerkin方法的激波捕获,(复杂应用的有限体积VII-椭圆。复杂应用的无限体积VII-椭圆型、抛物型和双曲型问题(2014),Springer),945-953·Zbl 1426.76429号
[70] Harten,A.,双曲守恒律的高分辨率格式,J.Compute。物理。,49, 3, 357-393 (1983) ·Zbl 0565.65050号
[71] Van Leer,B.,走向最终的保守差分方案。V.戈杜诺夫方法的二阶续集,J.Compute。物理。,32, 1, 101-136 (1979) ·Zbl 1364.65223号
[72] Sonntag,M.,不连续Galerkin方法中Navier-Stokes方程的形状导数和激波捕获(2017),斯图加特大学(博士论文)
[73] 杜姆布瑟。;O.扎诺蒂。;卢布埃,R。;Diot,S.,《双曲守恒律间断Galerkin有限元法的后验子单元限制》,J.Compute。物理。,278,47-75(2014)·Zbl 1349.65448号
[74] A.Jameson,W.Schmidt,E.Turkel,使用龙格-库塔时间步进格式通过有限体积方法对欧拉方程进行数值求解,载于:第14届流体和等离子体动力学会议,1981年,第1259页。
[75] Ducros,F。;费兰德,V。;Nicoud,F。;韦伯,C。;Darracq,D。;Gacherieu,C。;Poinsot,T.,冲击/湍流相互作用的大涡模拟,J.Compute。物理。,152, 2, 517-549 (1999) ·Zbl 0955.76045号
[76] Williamson,J.,低存储Runge-Kutta方案,J.Compute。物理。,35, 1, 48-56 (1980) ·Zbl 0425.65038号
[77] 库兰特,R。;弗里德里希斯(Friedrichs,K.)。;Lewy,H.,《数学差异》。《年鉴》,100,1,32-74(1928)·JFM 54.0486.01型
[78] Cockburn,B。;Shu,C.-W.,对流占优问题的Runge-Kutta间断Galerkin方法,科学杂志。计算。,16, 3, 173-261 (2001) ·Zbl 1065.76135号
[79] 沃伯顿,T。;Hagstrom,T.,《结构网格上间断Galerkin方法的CFL数驯服》,SIAM J.Numer。分析。,46, 6, 3151-3180 (2008) ·兹比尔1181.35010
[80] C.Rumsey,B.Wedan,T.Hauser,M.Poinot,计算流体动力学通用表示系统(CGNS)的最新更新,载于:第50届美国航空航天协会航空航天科学会议,包括新视野论坛和航空航天博览会,2012年,第1264页。
[81] Fehn,N。;Wall,W.A。;Kronbichler,M.,《无矩阵高阶间断Galerkin可压缩Navier-Stokes解算器:湍流不可压缩流动的可压缩和不可压缩公式的性能比较》,国际。J.数字。方法流体,89,3,71-102(2019)
[82] Bolemann,T.,《湍流高阶间断Galerkin方法的工业化》(2018),斯图加特大学(博士论文)
[83] Resch,M。;Kaminski,A.,《技术在计算机模拟和机器学习中的认知重要性》,Minds Mach。,29, 1, 9-17 (2019)
[84] Peng,R.D.,《计算科学中的可再现研究》,《科学》,334,6060,1226-1227(2011)
[85] Chacon,S。;Straub,B.,Pro Git(2014年),Apress
[86] Gassner,G.J。;Beck,A.D.,关于欠分辨率湍流模拟的高阶离散精度,Theor。计算。流体动力学。,27, 3-4, 221-237 (2013)
[87] 辛登朗,F。;Gassner,G.J。;Altmann,C。;贝克,A。;Staudenmaier,M。;Munz,C.-D.,非定常问题的显式间断Galerkin方法,计算与流体,61,86-93(2012)·Zbl 1365.76117号
[88] D.Flad,A.D.Beck,G.Gassner,C.-D.Munz,直接数值模拟气动声学的非连续Galerkin谱元方法,载于:第20届AIAA/CEAS气动声学会议,2014年,第2740页。
[89] 阿塔克,M。;贝克,A。;Bolemann,T。;弗拉德·D。;弗兰克·H。;Munz,C.-D.,使用高阶不连续伽辽金谱元方法的高保真度尺度解析计算流体动力学,(科学与工程中的高性能计算́15(2016),施普林格),511-530
[90] 弗兰克·H·M。;Munz,C.-D.,侧视镜上声反馈现象的直接气动声学模拟,J.Sound Vib。,371, 132-149 (2016)
[91] Pirozzoli,S.,《高速流动的数值方法》,年。流体力学版次。,43, 163-194 (2011) ·Zbl 1299.76103号
[92] Drela,M.,XFOIL:低雷诺数翼型的分析和设计系统,(Mueller,T.J.,《低雷诺数空气动力学》(1989),施普林格-柏林-海德堡:施普林格–柏林-海德堡-柏林,海德堡),1-12
[93] Shu,C.-W.,双曲守恒律的本质非振荡和加权本质非振荡格式,(非线性双曲方程的高级数值逼近(1998),Springer),325-432·Zbl 0927.65111号
[94] Chandrashekar,P.,可压缩Euler和Navier-Stokes方程的动能守恒和熵稳定有限体积格式,Commun。计算。物理。,14, 5, 1252-1286 (2013) ·Zbl 1373.76121号
[95] Lund,T.S。;吴,X。;Squires,K.D.,《为空间发展边界层模拟生成湍流流入数据》,J.Compute。物理。,140,2233-258(1998年)·Zbl 0936.76026号
[96] 卡瓦伊,S。;Larsson,J.,《大涡模拟中的壁面模型:长度尺度、网格分辨率和精度》,Phys。流体,24,1,第015105条pp.(2012)
[97] L.Erbig,N.Hu,S.Lardeau,被动间隙噪声的实验和数值研究,载于:2018年AIAA/CEAS气动声学会议,2018年,第3595页。
[98] T.Kuhn,使用高阶非连续Galerkin方法对主动空腔噪声进行分区大涡模拟,载于:第25届AIAA/CEAS气动声学会议,2019年,第2465页。
[99] 费希特,S。;蒙兹,C.-D。;罗德,C。;Zeiler,C.,《具有相变和表面张力的可压缩液-气流动的锐界面法》,J.Compute。物理。,336, 347-374 (2017) ·Zbl 1375.76194号
[100] 亨佩特,F。;Boblest,S。;厄特尔,T。;萨德洛,F。;奥芬哈·user,P。;玻璃,C。;霍夫曼,M。;贝克,A。;蒙兹,C.-D。;Iben,U.,使用不连续Galerkin谱元法的状态方程表模拟超音速甲烷射流中的真实气体效应,计算与流体,145167-179(2017)·Zbl 1390.76273号
[101] 贝克,A。;奥特温,P。;Kopper,P。;北卡罗来纳州克莱斯。;Kempf,D。;Koch,C.,《朝向高保真侵蚀预测:叶轮机械中的时间精确颗粒追踪》,《国际热流杂志》,79,108457(2019)
[102] Schnücke,G。;北克雷斯。;Bolemann,T。;Gassner,G.J.,双曲守恒律移动网格上熵稳定的间断Galerkin格式,J.Sci。计算。(2020) ·Zbl 1437.65142号
[103] 科普莱斯通,S.M。;Pfeiffer,M。;Fasoulas,S。;Munz,C.-D.,激光与等离子体相互作用的高阶粒子-细胞模拟,《欧洲物理》。J.规格顶部。,227, 14, 1603-1614 (2019)
[104] 奥特温,P。;铜石,S。;蒙兹,C.-D。;粘合剂T。;米尔扎。;Nizenkov,P。;Pfeiffer,M。;Reschke,W。;Fasoulas,S.,Piclas:用于模拟反应等离子体流的高度灵活的粒子代码,(2017 IEEE国际等离子体科学会议,2017 IEEE等离子体科学国际会议,ICOPS(2017),IEEE),1
[105] 库恩,T。;Dürrwächter,J。;梅耶,F。;贝克,A。;罗德,C。;Munz,C.-D.,空腔流直接气动声学模拟的不确定性量化,J.Compute。阿库斯特。(2018)
[106] H.J.本加兹。;林德纳,F。;Gatzhammer,B。;梅尔,M。;Scheufele,K。;Shukaev,A。;Uekermann,B.,PreCICE-多物理表面耦合的完全并行库,流体-结构相互作用进展。流体-结构相互作用、计算与研究进展流体,141250-258(2016)·Zbl 1390.76004号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。