×

在多个图形处理单元上使用上下对称Gauss-Seidel预处理器求解超音速和高超音速三维流动问题的隐式方法。 (英语) Zbl 1453.76089号

摘要:本文描述了一种在非结构化网格上求解定常气体动力学三维空间方程的数值方法,该方法是为多图形处理单元(GPU)计算架构设计的。使用一阶和二阶TVD格式对控制方程进行离散化。采用带有简单伪时间步长同伦的牛顿法求解该问题。每个迭代步骤都涉及到源于气体动力学方程线性化的线性系统的求解。使用Krylov子空间迭代方法来求解线性系统。本文的主要目的是描述一种预处理的上下对称高斯-赛德尔(LU-SGS)方法及其在多个GPU计算系统上的适应性。研究表明,故意重新排序的矩阵具有重新排列的求解过程,可以产生上下三角线性系统,从而获得与原始单线程LU-SGS相近的代数性质。该方法以公布的结果为基准。对马赫数在1.2到25之间的不同流动进行了计算效率和加速度分析。

MSC公司:

76个M12 有限体积法在流体力学问题中的应用
65F08个 迭代方法的前置条件
65平方米 含偏微分方程初值和初边值问题离散方程的数值解
76K05美元 高超音速流
76J20型 超音速流动
2005年5月 并行数值计算
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] A.詹姆逊。;Yoon,S.,Euler方程的多网格Lower-upper隐式格式,AIAA J.,25,929-935(1987)
[2] Kim,J。;Kwon,O.,非结构网格Navier-Stokes计算中块Lu-sgs格式的改进,(第40届美国航空航天局航空科学会议和展览(2002))
[3] 沙洛夫,D。;Nakahashi,K.,《用于上下对称高斯-赛德尔计算的混合非结构网格的重新排序》,AIAA J.,36,484-486(1998)
[4] 沙洛夫,D。;罗,H。;约瑟夫·D·B。;Lohner,R.,《非结构化网格gmres+Lu-sgs方法在共享内存、基于缓存的并行计算机上的实现》(2012年),AIAA-2000-927-航空航天科学会议和展览
[5] 罗,H。;鲍姆·J·D。;Lohner,R.,《非结构网格上可压缩流动的快速无矩阵隐式方法》,J.Compute。物理。,146, 664-690 (1998) ·兹伯利0931.76045
[6] Elsen,E。;Legresley,P。;Darve,E.,使用gpu对高超声速飞行器上方的流量进行大型计算,J.Compute。物理。,227, 10148-10161 (2008) ·Zbl 1218.76035号
[7] 赖,J。;李,H。;田,Z。;Zhang,Y.,高超声速流动计算中的多gpu并行算法,数学。问题。工程,2019年,第2053156条pp.(2019)·Zbl 1435.76033号
[8] 埃梅利亚诺夫,V。;Karpenko,A。;Volkov,K.,在图形处理器单元上模拟具有平衡化学反应的高超音速流动,《宇航员学报》。,163,A部分,259-271(2019)
[9] Bocharov,A。;Evstigneev,N。;Ryabkov,O.,异质团簇结构上宽马赫数范围内的计算气体动力学,物理学杂志。Conf.序列号。,653,第012119条pp.(2015)
[10] 门肖夫,I。;Pavlukhin,P.,gpu加速集群上气体动力学隐式无矩阵解算器的高度可扩展实现,J.Supercomput。,73, 2, 631-638 (2017)
[11] 艾莎,M。;Verstraete,T。;Vuik,C.,关于结构网格上显式和隐式cfd模拟的gpu感知比较,计算。数学。申请。,74, 1, 201-217 (2017) ·Zbl 1394.65141号
[12] 嘉乐,Z。;M.志华。;C.洪泉。;Cheng,C.,可压缩流动的gpu加速隐式无网格方法,J.Compute。物理。,360, 39-56 (2018) ·Zbl 1391.76591号
[13] 沃特金斯,J。;罗梅罗,J。;Jameson,A.,Multi-gpu,非结构化网格上高阶方法的隐式时间步进,(第46届AIAA流体动力学会议(2016))
[14] Tuan,N。;卡斯通圭,P。;Laurendeau,E.,用于多块网格上三维空气动力学模拟的多重网格trans解算器的Gpu并行化,J.Supercomput。,75,52562-2583(2019)
[15] Bocharov,A。;比特尤林,V。;Evstigneev,N.,《在图形处理单元上使用完全隐式方法求解耦合气体热质传递问题》,J.Phys。Conf.序列号。,1147,第012053条pp.(2019)
[16] Sutherland,W.,《气体粘度和分子力》,Philos。Mag.,5,36,507-531(1893)
[17] 伯杰,M。;Aftosmis,M.,《不规则网格上坡度限制器的分析》,(第43届美国航空航天局航空科学会议和展览(2012年))
[18] Breuss,M.,《Lax-Friedrichs方法的正确使用》,ESAIM:数学。模式。编号。分析。,38, 3, 519-540 (2004) ·Zbl 1077.65089号
[19] Liou,M.S.,《ausm的续集,第二部分:所有速度的ausm+-up》,J.Compute。物理。,214, 1, 137-170 (2006) ·Zbl 1137.76344号
[20] 瓦达,Y。;Liou,M.-S.,《冲击和接触不连续性的精确而稳健的通量分裂方案》,SIAM J.Sci。计算。,18, 633-657 (1997) ·Zbl 0879.76064号
[21] Kim,K.H。;Kim,C。;Rho,O.-H.,高超声速气流精确计算方法:I.ausmpw+方案,J.Compute。物理。,174, 1, 38-80 (2001) ·Zbl 1106.76421号
[22] 鲁迪·D·H。;Strikwerda,J.,亚音速可压缩Navier-Stokes计算的边界条件,计算。流体,9327-338(1981)·Zbl 0452.76042号
[23] 黄,A.C。;奥尔马拉斯,S。;加尔布雷思,M。;Darmofal,D.,Navier-Stokes方程的井位亚音速流入/流出边界条件,(2018年AIAA航空航天科学会议(2018)),1-21
[24] 卡特希尔,E.H。;McKee,J.,《减少稀疏对称矩阵的带宽》,(《第24届美国国家会议论文集》,《第24次美国国家会议文献集》,纽约(1969年),157-172
[25] Gibbs,J。;吉布斯,E。;威廉·G·普尔(William G.Poole)。;Stockmeyer,P.K.,《减少稀疏矩阵带宽和轮廓的算法》,SIAM J.Numer。分析。,13, 2, 236-250 (1976) ·兹伯利0329.65024
[26] Karypis,G。;库马尔,V.,《不规则图的多级k-路划分方案》,J.并行分布计算。,48, 1, 96-129 (1998)
[27] Karypis,G。;Kumar,V.,《划分不规则图的快速高质量多级方案》,SIAM J.Sci。计算。,20, 1, 359-392 (1998) ·兹比尔0915.68129
[28] Tadjouddine,S.A.F.M。;秦,N.,《消除法》,适用于隐式可压缩cfd求解器的雅可比集合,国际J·数值。方法流体,4711315-1321(2005)·Zbl 1155.76384号
[29] Sleijpen,G.L。;Fokkema,D.R.,Bickstab(l),关于涉及具有复谱的非对称矩阵的线性方程组,电子。事务处理。数字。分析。,1, 11-32 (1993) ·兹比尔0820.65016
[30] 阿里亚斯,O。;O.Falcinelli。;菲科,N。;Elaskar,S.,使用Jameson、maccormack、shu和tvd esquemes对naca 0012上方流动的有限体积模拟,Mec。计算。,26, 3097-3116 (2007)
[31] M.,C.,5.2a,带凸点通道中的三维超声速流动(2008),[在线]
[32] Gnoffo,P.,《四面体网格上高超音速加热模拟的多维无粘通量重建》,(第47届美国航空航天局航空科学会议,包括新视野论坛和航空航天博览会(2009)),1-18
[33] 坎德勒,G。;Barnhardt,M。;Drayna,T。;Nompelis,I。;彼得森,D。;Subbareddy,P.,准确空气加热模拟的非结构化网格方法(2007年)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。