×

基于pyfr的混合非结构化网格上的异构计算。 (英语) Zbl 1390.76014号

概述:PyFR是一个用于非结构化网格的开源高精度计算流体动力学求解器。在本文中,我们详细介绍了如何将PyFR扩展到在混合元素网格和一系列硬件平台上运行,包括异构多节点系统。我们的实现的性能是使用圆柱周围空隙空间的纯六面体网格和混合棱柱四面体网格进行基准测试的。具体来说,在三个不同的硬件平台上,对每个网格的性能进行不同精度等级的评估;NVIDIA Tesla K40c GPU、Intel Xeon E5-2697 v2 CPU和AMD FirePro W9100 GPU。然后在由上述硬件组合而成的异构多节点系统上评估性能。结果表明,PyFR实现了跨各种硬件平台的性能可移植性。特别是,PyFR能够使用其“本地”语言针对各个平台,从而显著提高了仅使用OpenCL针对每个平台的性能。PyFR还被发现在异构多节点系统上具有良好的性能,实现了可用FLOP/s的很大一部分。最后,使用NVIDIA K20c GPU集群,使用每个网格对雷诺数为3900的圆柱体上方的非定常流动进行名义上五阶精度的长时模拟。对尾迹的长期动力学进行了详细研究,发现结果与先前的实验/数值数据非常吻合。所有结果都是通过PyFR v0.2.2获得的,该版本根据3条款新风格BSD许可证免费提供(请参见网址:http://www.pyfr.org).

MSC公司:

76-04 流体力学相关问题的软件、源代码等
65日元10 特定类别建筑的数值算法
76M10个 有限元方法在流体力学问题中的应用
76牛顿 可压缩流体和气体动力学
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] ReedWH,HillTR公司。中子输运方程的三角网格法。1973年,技术报告LA-UR-73-479,洛斯阿拉莫斯科学实验室。
[2] David A.Kopriva。;Kolias,John H.,可压缩流动的保守交错网格Chebyshev多域方法,计算物理杂志,125,1,244-261,(1996)·兹伯利0847.76069
[3] 孙玉芝;王志坚。;Liu,Yen,非结构化六面体网格上Navier-Stokes方程的高阶多域谱差分方法,Commun Comput Phys,2,2,310-333,(2007)·Zbl 1164.76360号
[4] Huynh,H.T.,高阶格式的通量重建方法,包括间断Galerkin方法,AIAA Pap,4079:2007,(2007)
[5] Hesthaven,Jan S。;Tim Warburton,Nodal间断Galerkin方法:算法、分析和应用,54,(2008),纽约:Springer Verlag·Zbl 1134.65068号
[6] 安德烈亚斯·科克纳;蒂姆·沃伯顿;杰夫·布里奇;Hesthaven,Jan S.,图形处理器上的Nodal非连续Galerkin方法,计算物理杂志,228,21,7863-7882,(2009)·Zbl 1175.65111号
[7] Patrice Castonguay;大卫·M·威廉姆斯(David M.Williams)。;Vincent,Peter E。;曼努埃尔·洛佩兹(Manuel Lopez);安东尼·詹姆逊(Antony Jameson),《关于为混合非结构网格开发高阶、多gpu、可压缩粘性流解算器》,AIAA Pap,3229:2011,(2011)
[8] F.D.威瑟顿。;Farrington,A.M。;Vincent,P.E.,Pyfr:一个开源框架,用于使用流量重建方法解决流架构上的对流-扩散型问题,Compute Phys Commun,185,11,3028-3040,(2014)·兹比尔1348.65005
[9] 文森特,体育。;卡斯通圭,P。;Jameson,A.,一类新的高阶能量稳定通量重建方案,J Sci-Comput,47,1,50-72,(2011)·Zbl 1433.76094号
[10] 卡斯通圭,P。;文森特,P。;Jameson,A.,三角形单元的一类新的高阶能量稳定通量重建方案,科学计算杂志,51,1,224-256,(2011)·Zbl 1457.65101号
[11] 威廉姆斯,D.M。;Jameson,A.,四面体平流扩散问题的能量稳定通量重建方案,科学计算杂志,59,3,721-759,(2013)·Zbl 1303.65076号
[12] 日本东京,Kazushige;Geijn,Robert A,《高性能矩阵乘法剖析》,ACM Trans Math Softw,34,3,12,(2008)·Zbl 1190.65064号
[13] 拜耳迈克尔。Mako:python模板。2013
[14] Klöckner,安德烈亚斯;尼古拉斯·平托(Nicolas Pinto);Lee,Yunsup;布赖恩·卡坦扎罗(Bryan Catanzaro);保罗·伊万诺夫(Paul Ivanov);Fasih、Ahmed、Pycuda和pyopencl:基于脚本的gpu运行时代码生成方法,并行计算,38,3,157-174,(2012)
[15] 佛梅尔,公元前。;Nadarajah,S.,《通过重建方案使用校正程序的ILES自适应IMEX时间步进》,AIAA Pap,2687:2013,(2013)
[16] 佛梅尔,公元前。;Nadarajah,S.,《高阶非结构化方法的自适应IMEX方案》,《计算物理杂志》,280,261-286,(2015)·Zbl 1349.76545号
[17] 罗斯科·阿纳托尔。涡街湍流尾迹的发展。1953.技术报告编号NACA TR 1191,加利福尼亚理工学院。
[18] Bloor,M.S.,《圆柱尾迹向湍流的过渡》,《流体力学杂志》,第19卷,第290-304页,(1964年)·Zbl 0119.41703号
[19] 威廉姆森,C.H.K,圆柱尾迹向三维过渡过程中两个阶段的存在,《物理流体》,313165-3168,(1988)
[20] 威廉姆森,C.H.K,《圆柱尾迹中的旋涡动力学》,《流体力学年鉴》,第28期,第477-539页,(1996年)
[21] Norberg,C.,圆柱近尾迹的Ldv测量,国际数值方法流体,28,9,1281-1302,(1998)·Zbl 0933.76041号
[22] 马,X。;卡拉马诺斯,G.S。;Karniadakis,G.E.,湍流近尾流的动力学和低维,流体力学杂志,310,29-65,(2000)·Zbl 0987.76041号
[23] Breuer,M.,《圆柱绕流亚临界流的大涡模拟》,《国际数值方法流体》,28,9,1281-1302,(1998)·兹伯利0933.76041
[24] Kravchenko,A.G。;Moin,P.,《(r e_d=3900)下圆柱绕流的数值研究》,《物理流体》,第12期,第403-417页,(2000)·Zbl 1149.76441号
[25] 菲利普·帕尔纳多;Johan Carlier;多米尼克·海茨;Lamballais,Eric,雷诺数为3900时圆柱绕流的实验和数值研究,Phys Fluids,20,8,085101-085101-14,(2008)·Zbl 1182.76591号
[26] O.莱姆库尔。;罗德里格斯,I。;博雷尔,R。;Oliva,A.,圆柱体涡旋形成区域的低频不稳定性,Phys Fluids,25,8,3165-3168,(2013)
[27] 克里斯托弗·肯尼迪(Christopher A.Kennedy)。;木匠,Mark H。;Lewis,R.Michael,可压缩Navier-Stokes方程的低存储显式Runge-Kutta格式,应用数值数学,35,3,177-219,(2000)·Zbl 0986.76060号
[28] 佛梅尔,公元前。;Cagnone,J.S。;Nadarajah,S.,ILES通过重建方案使用校正程序,AIAA Pap,1001:2013,(2013)
[29] 佛梅尔,公元前。;Nadarajah,S。;Tucker,P.G.,高阶非结构化隐式大涡模拟的标准试验案例,AIAA Pap,0935:2014,(2014)
[30] 乔治·卡里皮斯;Kumar,Vipin,一种快速高质量的不规则图分割多级方案,SIAM科学计算杂志,20,1,359-392,(1998)·Zbl 0915.68129号
[31] 弗朗索瓦佩莱格里尼;Roman,Jean,Scotch:《通过过程和架构图的双重递归双分区实现静态映射的软件包》,《高性能计算与网络》,493-498,(1996),Springer
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。