×

使用FLASH并行Navier-Stokes模拟任意几何形状的高速可压缩流。 (英语) Zbl 1390.76676号

摘要:我们报告了对FLASH代码的扩展,以实现对任意二维和三维静止物体的高速可压缩粘性流模拟。通过实现适当的计算机图形算法,将人体形状嵌入到块结构的笛卡尔自适应网格细化网格中。精确的身体形状表示需要较高的网格细化级别,这会导致较大的网格尺寸,特别是在三维模拟中。在IBM Blue Gene/Q计算系统上并行进行了仿真,并在该系统上评估了纯MPI和混合MPI-OpenMP模式下的代码性能。我们还在FLASH中实现了适当的壁面边界条件,以模拟粘壁效应。文中给出了各种二维试验情况下的Navier-Stokes(NS)解,如激波边界层相互作用问题,以及钝锥-圆柱-火焰和双锥几何体的高超声速流。对用于高超声速流动控制的微涡发生器进行了三维NS模拟,计算结果与实验结果一致。

MSC公司:

76米25 其他数值方法(流体力学)(MSC2010)
2005年5月 并行数值计算
PDF格式BibTeX公司 XML格式引用
全文: DOI程序

参考文献:

[1] FLASH用户指南-版本4.0;2012
[2] A.杜比。;Antypas,K。;Ganapathy,M.K。;里德,L.B。;莱利,K.M。;Sheeler,D.,FLASH的可扩展组件架构,一种大规模并行多物理仿真代码,并行计算,35,512-522,(2009)
[3] Fryxell,B。;Olson,K。;Ricker,P。;蒂姆斯,F.X。;Zingale,M。;Lamb,D.Q.,FLASH:用于模拟天体物理热核闪光的自适应网格流体动力学代码,《天体物理学杂志》增刊,131,273-334,(2000)
[4] 伍德沃德,P。;Colella,P.,《强冲击下二维流体流动的数值模拟》,《计算物理杂志》,54,115-173,(1984)·Zbl 0573.76057号
[5] Gorobets,A。;Trias,F.X。;Oliva,A.,用于不可压缩流混合超级计算的并行MPI+openmp+opencl算法,计算流体,88,764-772,(2013)·Zbl 1391.76090号
[6] Ye,Y。;Li,K.,基于熵格子Boltzmann方法在GPU上使用CUDA进行高雷诺数流动模拟,计算流体,88,241-249,(2013)·Zbl 1391.76645号
[7] 李凯。;杨伟(Yang,W.)。;Li,K.,使用概率建模对GPU上spmv进行性能分析和优化,IEEE跨并行分布系统,99,(2014)
[8] 梁,J。;Li,K.,利用GPU和CUDA对大型液态金属原子系统的加速凝固过程模拟,计算物理杂志,257,521-535,(2014)·Zbl 1349.82093号
[9] 约翰·B。;艾默生·D·R。;Gu,X.J.,使用FLASH代码的并行可压缩粘性流模拟:任意三维几何体的实现,Procedia Eng,61,52-56,(2013)
[10] 伍德沃德,P。;Colella,P.,用于气体动力学模拟的分段抛物线方法(PPM),《计算物理杂志》,54,174-201,(1984)·兹比尔0531.76082
[11] 加德纳,T。;Stone,J.,通过三维约束传输实现理想MHD的非分裂Godunov方法,计算物理杂志,2274123-4141,(2008)·Zbl 1317.76057号
[12] Lee,D。;Deane,A.,《多维磁流体力学的非分裂交错网格格式》,《计算物理杂志》,228952-975,(2009)·Zbl 1330.76093号
[13] Lee,D.,三维磁流体力学解的精确、高效和稳定的非分裂交错网格方案,计算机物理,243269-292,(2013)·Zbl 1349.76494号
[14] Olson K.PARAMESH:一种并行自适应网格工具。在:并行计算流体动力学:理论和应用-并行计算流体动力学2005会议记录。爱思唯尔;2006年,第341页。
[15] 伯杰,M.J。;Colella,P.,《冲击流体动力学的局部自适应网格细化》,《计算物理杂志》,82,64-84,(1989)·Zbl 0665.76070号
[16] Carvalho P,Cavalcanti P。使用球形多边形进行多面体测试的点。Graphics Gems V.学术出版社;1995
[17] White,F.M.,《粘性流体流动》(2005),麦格劳-希尔
[18] Toro,E.F.,Riemann解算器和流体动力学数值方法,实用简介,(2009),Springer·Zbl 1227.76006号
[19] Zhang M,Yu J,Chang SC。用CESE方法求解Navier-Stokes方程。AIAA论文2004-0075;2004
[20] 江,C。;崔,S。;Feng,X.,用AMR-CESE方法求解Euler和Navier-Stokes方程,计算流体,54,105-117,(2012)·Zbl 1291.76244号
[21] Hakkinen RJ、Greber I、Trilling L、Abarbanel SS。斜激波与层流边界层的相互作用。NASA备忘录2-18-59W;1959
[22] 姚,Y。;Krishnan,L。;Sandham,N.D。;Roberts,G.,《马赫数对激波/边界层相互作用中不稳定扰动的影响》,《物理流体》,19054104,(2007)·Zbl 1146.76572号
[23] Gray J,Lindsay E.标准超高速弹道模型HB-1和HB-2在1.5至10马赫下的力测试,AEDC-TDR-64-246;1964
[24] Gilchrist AR,Williams MJ。M等于7.5时AGARD模型HB1和HB2上的压力分布和力。技术代表,空气动力学注释346,ARL;1974
[25] Tissera S.高超声速真实气流中高分辨率方法的评估。博士论文。克兰菲尔德大学;2010
[26] Birch T,Prince S,Ludlow DK,Qin N。抛物线化Navier-Stokes解算器在某些高超声速流动问题中的应用。AIAA论文2001-1753;2001
[27] Coblish JJ、Smith MS、Hand T、Candler GV、Nompelis I.AEDC第9号超高速风洞的双锥体实验和数值分析。AIAA论文2005-0902;2005
[28] Babinsky,H。;李毅。;Ford,C.W.P.,超音速斜激波/边界层相互作用的微斜面控制,AIAA J,47,668-675,(2009)
[29] Li Q,Liu C.LES,用于在M=2.5和Reθ=1440下使用MVG的超声速斜坡控制流。AIAA论文2010-0592;2010
[30] 霍尔顿,H。;Babinsky,H.,微体转运体发生器对分离的正常冲击/边界层相互作用的影响,J Aircr,44170-174,(2007)
[31] Anderson BH,Tinapple J,Surber L.使用微阵列驱动优化控制激波-湍流边界层相互作用。AIAA论文2006-3197;2006
[32] 莫德,R。;侯赛因,Z。;阿扎姆,C。;Konstantinos,K.,《用于高超音速流量控制的微型放大器》,《微型机器》,第3364-378页,(2012年)
[33] HDF5用户指南。版本1.8.7;2011
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。