×

\不可压缩湍流欠分辨率模拟的(p)-多重网格无矩阵间断Galerkin解策略。 (英语) Zbl 1519.76144号

摘要:近年来,几项研究工作集中于发展用于湍流尺度分解模拟的高阶间断Galerkin(dG)方法。然而,在不可压缩流计算的背景下,求解以不定雅可比矩阵为特征的大规模方程组的计算费用往往阻碍了工业相关计算的模拟。在这项工作中,我们试图通过设计健壮、可扩展和内存级的解决策略来提高Rosenbrock型线性隐式Runge-Kutta方法的效率。特别地,我们引入了内存节省的多重网格预条件,耦合了无矩阵和基于矩阵的Krylov迭代平滑器。多重网格预处理程序依赖于精细空间上廉价的元素-线路块-对角平滑器,以降低组装成本和内存分配,并使用加性Schwarz平滑器确保多重网格迭代的最粗空间有足够的分辨率,以获得满意的收敛速度和最佳的并行效率。此外,还探讨了如何使用特制的重定标继承粗算子来克服精细空间算子的标准继承所提供的过度稳定性。进行了广泛的数值验证。Rosenbrock公式被应用于越来越复杂的测试用例:二维圆柱体周围的层流非定常流动(Re=200)和球体周围的层压非定常流(Re=300),ERCOFTAC T3L测试用例套件中具有不同自由流紊流水平的过渡流问题。作为概念的证明,本文报道了波音公司在(R e=10^6)条件下基本起落架试验案例的数值解。记录了解决方案与实验数据的良好一致性,而据报告,相对于最先进的解决方案策略,内存占用减少了约92%,执行时间增加了高达3.5。

MSC公司:

76M10个 有限元方法在流体力学问题中的应用
65M60毫米 涉及偏微分方程初值和初边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65M55型 多重网格方法;涉及偏微分方程初值和初边值问题的区域分解
76层65 湍流的直接数值模拟和大涡模拟

软件:

硝基苯砜
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Bassi,F。;博蒂,L。;科伦坡,A。;Crivellini,A。;Ghidoni,A。;Massa,F.,《关于湍流DNS和隐式LES隐式高阶间断Galerkin方法的发展》,《Eur J Mech-B/Fuids》,55,367-379(2016)·Zbl 1408.76360号
[2] Chapelier,J.-B。;De La Llave Plata,M。;雷纳克,F。;Lamballais,E.,湍流DNS的高阶非连续Galerkin方法评估,计算流体,95,210-226(2014)·Zbl 1391.76209号
[3] 维亚特,C.C。;Hillewaert,K。;Bricteux,L。;Winckelmans,G.,基于间断Galerkin/对称内罚方法的自由和边界湍流的隐式大涡模拟,国际J数值方法流体,78,6,335-354(2015)
[4] Bassi,F。;Crivellini,A。;Di Pietro,D。;Rebay,S.,不可压缩定常和非定常流动的隐式高阶不连续伽辽金方法,Comput Fluids,36,1011529-1546(2007)·Zbl 1194.76102号
[5] Uranga,A.公司。;佩尔松,P.-O。;Drela,M。;Peraire,J.,使用间断Galerkin方法对低雷诺数下湍流过渡的隐式大涡模拟,国际J数值方法工程,87,1-5,232-261(2011)·Zbl 1242.76085号
[6] Bassi,F。;博蒂,L。;科伦坡,A。;Ghidoni,A。;Massa,F.,应用于可压缩和不可压缩非定常流的间断Galerkin解的线性隐式Rosenbrock型Runge-Kutta格式,计算流体,118,305-320(2015)·Zbl 1390.76833号
[7] Crivellini,A。;Bassi,F.,粘性和湍流空气动力学模拟的隐式无矩阵间断Galerkin解算器,计算流体,50,1,81-93(2011)·Zbl 1271.76164号
[8] Sarshar,A。;特兰基利,P。;Pickering,B。;McCall,A。;罗伊·C·J。;Sandu,A.,大型CFD模拟中无矩阵隐式时间步长方法的数值研究,计算流体,159,53-63(2017)·Zbl 1390.76688号
[9] Franciolini,M。;克里韦里尼,A。;Nigro,A.,关于不可压缩湍流高阶间断Galerkin解的无矩阵线性隐式时间积分策略的效率,计算流体,159276-294(2017)·兹比尔1390.76312
[10] Crivellini,A。;Franciolini,M。;Nigro,A.,湍流模拟的无矩阵不可压缩DG算法,湍流进展VII,153-159(2017),Springer
[11] 迪奥萨迪,L.T。;Murman,S.M.,高阶时空间断Galerkin方法的张量积预条件,计算物理杂志,330,296-318(2017)·Zbl 1378.76046号
[12] Pazner,W。;Persson,P.-O.,非常高阶间断Galerkin方法的近似张量积预条件,计算物理杂志,354,344-369(2018)·Zbl 1380.65067号
[13] 迪奥萨迪,L.T。;Murman,S.M.,《高阶有限元方法的可伸缩张量积预条件:标量方程》,《计算物理杂志》(2019)·Zbl 1452.65050号
[14] 海伦布鲁克,B。;Mavrilis,D。;Atkins,H.,连续和非连续有限元离散化的p-多重网格分析,第16届AIAA计算流体动力学会议,3989(2003)
[15] 海伦布鲁克,B.T。;Atkins,H.L.,p-多重网格在泊松方程间断Galerkin公式中的应用,AIAA J,44,3,566-575(2006)
[16] Bassi,F。;Rebay,S.,用多阶间断有限元法数值求解欧拉方程,计算流体动力学2002,199-204(2003),Springer·Zbl 1140.76360号
[17] Fidkowski,K.J。;奥利弗,T.A。;卢,J。;Darmofal,D.L.,可压缩Navier-Stokes方程高阶间断Galerkin离散的p-多重网格解,计算物理杂志,207,1,92-113(2005)·Zbl 1177.76194号
[18] 普里尔,F。;卢卡科夫-梅德维奥娃,M。;Hartmann,R.,不连续Galerkin方法的平滑聚合多重网格,SIAM科学计算杂志,31,5,3503-3528(2009)·兹比尔1200.35509
[19] Wallraff,M。;哈特曼,R。;Leicht,T.,DG方法的多重网格解算器算法及其在空气动力流中的应用,IDIHOM:高阶方法的工业化——自顶向下方法,153-178(2015),Springer
[20] Antonietti,P.F。;萨蒂,M。;Verani,M.,椭圆问题hp-连续Galerkin离散化的多重网格算法,SIAM J Numer Ana,53,1,598-618(2015)·Zbl 1312.65181号
[21] 沙巴齐,K。;Mavrilis,D.J。;Burgess,N.K.,可压缩Navier-Stokes方程高阶间断Galerkin离散的多重网格算法,计算物理杂志,228,21,7917-7940(2009)·Zbl 1391.65181号
[22] 迪奥萨迪,L.T。;Darmofal,D.L.,Navier-Stokes方程间断Galerkin解的预处理方法,计算物理杂志,228,11,3917-3935(2009)·Zbl 1185.76812号
[23] 博蒂,L。;科伦坡,A。;Bassi,F.,H-不可压缩流动问题间断Galerkin离散的基于多重网格聚集的求解策略,《计算物理杂志》,347382-415(2017)·Zbl 1380.65251号
[24] Crivellini,A。;Franciolini,M。;Nigro,A.,一种用于湍流模拟的具有减少记忆足迹的隐式不连续伽辽金方法,直接和大涡模拟XI,69-74(2019),施普林格
[25] Bassi,F。;Rebay,S.,可压缩Navier-Stokes方程数值解的高精度间断有限元方法,计算物理杂志,131267-279(1997)·Zbl 0871.76040号
[26] 布雷齐,F。;Manzini,G。;Marini,D。;Pietra,P。;Russo,A.,椭圆问题的间断Galerkin逼近,数值方法部分微分Equ,16,365-378(2000)·Zbl 0957.65099号
[27] Arnold,D.N。;布雷齐,F。;Cockburn,B。;Marini,L.D.,椭圆问题间断Galerkin方法的统一分析,SIAM J Numer Ana,39,5,1749-1779(2002)·Zbl 1008.65080号
[28] Bassi,F。;Crivellini,A。;Di Pietro,D。;Rebay,S.,不可压缩Navier-Stokes方程间断Galerkin解的人工压缩性通量,计算物理杂志,218,2794-815(2006)·Zbl 1158.76313号
[29] Crivellini,A。;D’Alessandro,V。;Bassi,F.,《不可压缩三维Navier-Stokes方程的高阶间断Galerkin方法评估:通过re=500球体的流量基准结果》,计算流体,86,0,442-458(2013)·兹比尔1290.76022
[30] Bassi,F。;Ghidoni,A。;Perbellini,A。;Rebay,S。;Crivellini,A。;Franchina,N。;Savini,M.,《不可压缩RANS和(k-\omega)湍流模型方程的高阶间断Galerkin解算器》,计算流体,98,54-68(2014)·兹比尔1391.76301
[31] Rang,J。;Angermann,L.,《指数2 PDAE的3阶新Rosenbrock方法》,《Equadiff-11 2005年学报》,385-394(2007),夸美纽斯大学出版社(布拉迪斯拉发)
[32] Knoll,D。;Keyes,D.,无雅可比的Newton-Krylov方法:方法和应用综述,J Comput Phys,193,2357-397(2004)·Zbl 1036.65045号
[33] 佩妮斯,M。;Walker,H.F.,NITSOL:非线性系统的牛顿迭代求解器,SIAM科学计算杂志,19,1,302-318(1998)·Zbl 0916.65049号
[34] Crivellini,A。;Bassi,F.,粘性和湍流空气动力学模拟的隐式无矩阵间断Galerkin解算器,计算流体,50,1,81-93(2011)·Zbl 1271.76164号
[35] http://www.mcs.anl.gov/petsc
[36] 博蒂,L。;科伦坡,A。;Crivellini,A。;弗朗西奥利尼,M。,{h-p-hp}-多级椭圆算子的间断Galerkin解策略,国际计算流体力学杂志,1-9(2019)·Zbl 07474500号
[37] Antonietti,P.F。;阿尤索·德·迪奥斯(Ayuso de Dios),B。;Brenner,S。;Sung,L.,Schwarz方法用于椭圆问题的预处理WOPSIP方法,计算方法应用数学,12,3,241-272(2012)·Zbl 1284.65162号
[38] Saad,Y.,一种灵活的内外预处理GMRES算法,SIAM科学计算杂志,14,2,461-469(1993)·Zbl 0780.65022号
[39] 梅内吉尼,J.R。;萨尔塔拉,F。;西奎拉,C.L.R。;Ferrari Jr,J.A.,串联和并排布置的两个圆柱体之间流动干扰的数值模拟,《流体结构杂志》,15,2,327-350(2001)
[40] Botti,L.,参考-物理框架映射对不连续分段多项式空间逼近性质的影响,科学计算杂志,52,3,675-703(2012)·Zbl 1255.65222号
[41] 约翰逊,T.A。;Patel,V.C.,《雷诺数达300的球体绕流》,《流体力学杂志》,378,19-70(1999)
[42] Tomboulides,A.G。;Orszag,S.A.,通过球体的过渡和弱湍流的数值研究,《流体力学杂志》,41645-73(2000)·Zbl 1156.76419号
[43] Ploumhans,P。;温克尔曼,G。;Salmon,J.K。;伦纳德(A.Leonard)。;Warren,M.S.,《三维钝体流动直接数值模拟的旋涡方法:在re=300、500和1000时球体的应用》,《计算物理杂志》,178、2、427-463(2002)·兹比尔1045.76030
[44] Crivellini,A。;Franciolini,M。;科伦坡,A。;Bassi,F.,非连续Galerkin解算器的OpenMP并行化策略,国际J并行程序,178,2,427-463(2018)
[45] Langari,M。;Yang,Z.,高架自由流湍流下分离边界层转捩初级不稳定性的数值研究,《物理流体》,25,7,074106(2013)
[46] 杨,Z。;Voke,P.R.,表面曲率变化时边界层分离和过渡的大涡模拟,《流体力学杂志》,439305333(2001)·Zbl 1020.76023号
[47] 希利尔,R。;Cherry,N.,《水流紊流对分离气泡的影响》,《风力工程与空气动力学杂志》,第8期,第149-58页(1981年)
[48] Y.中村。;Ozono,S.,《湍流对分离流和再附着流的影响》,《流体力学杂志》,178477-490(1987)
[49] 第三届国际微尺度传热与流体流动会议
[50] 兰巴拉斯,E。;西尔维斯特里尼,J。;Laizet,S.,《圆形前缘后流动分离的直接数值模拟:曲率效应研究》,《国际热流学杂志》,31,3,295-306(2010)
[51] Spalart,P.R.,《湍流边界层的直接模拟》,《流体力学杂志》,187,61-98(1988)·Zbl 0641.76050号
[52] Cutrone,L。;德帕尔马,P。;帕斯卡齐奥,G。;Napolitano,M.,《预测二维和三维分离流动中的转变》,《国际热流体流动杂志》,29,2,504-526(2008)
[53] 斯帕拉特,P。;舒尔,M。;Strelets,M。;Travin,A.,基本起落架的初始rans和ddes,混合rans-LES建模进展,101-110(2010),Springer
[54] 文卡塔克里希南。;Karthikeyan,N。;Mejia,K.,基本四轮起落架的实验研究,AIAA J,50,11,2435-2447(2012)
[55] S.-H.Peng,M.Strelets,测试案例ST12,波音基本起落架。,单位:ATAAC。
[56] 董庆林。;Xu,H.-Y。;Ye,Z.-Y.,使用DDES、LES和URANS对基本起落架非定常流动的数值研究,工程应用计算流体力学,12,1,689-710(2018)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。