基于GPU的非结构网格CFD求解器的设计与优化.pdf.pdfVIP

· 技术 ／TECHNOLOGY · 基于 GPU 的非结构 网格 CFD求解 器的设计与优化 宋慎义 ，王彦桐 ，刘冰 ，陆忠华 1．中国科学院计算机网络信息中心，北京 100190 2．中国科学院研究生院，北京 100190 摘 要： 本文介绍了在 GPU(图形处理器)上实现 CFD(计算流体力学)问题在非结构网格的有限体积法求 解，并针对非结构网格求解随机访问内存的瓶颈，结合 GPU的编程模型和硬件特性，提 出GPU并 行计算优化方法。实验表明，使用存储优化、内联函数优化、算术指令优化后的GPU算法相对于 原 GPU算法提高了硬件的利用效率，相对于CPU程序有明显的性能优势。 关键词： 图形处理器；并行优化；计算流体力学；非结构网格 DesignandOptimizationofanUnstructuredGridCFD Solver BasedonGPU SongShenyi一，W angYangang，LiuBing一，LuZhonghua ． ComputerNetworkInformationCenter,ChineseAcademyofSciences，Beijing100190,China 2．GraduateUniversityofChineseAcademyofSciences，Bejiing100190,China Abstract： ThispaperpresentsanimplementationofaCFD finitevolumeunstructuredgridsolverbasedonGPU andproposesanoptimizationmethodforGPU parallelcomputation，aimingtothebottleneckofRandom AccessMemoryforunstructuredgridsolver,withconsideringGPU programmingmodelandhardware architecture．TheresultindicatesthattheoptimizedGPU algorithm improvestheutilizationofGPU card andshowsanevidentadvantage． Keywords： GPU；Paralleloptimization；CFD：Unsturcturedgrid 基金项 目：国家高技术研究发展计划资助项 目(2009AA01A139)；国家 自然科学基金项 目 30 宋慎义 等：基于GPU的非缱 网格cF 鲢鲞 堡兰± 垡 1．引言 2．非结构网格上的有限体积法 CFD求解 计算流体力学 (ComputationalFluidDynamics， 2．1控制方程 简称 CFD1使用数值方法在计算机中对流体力学的控 就并行计算方法而言，Navier．Stokes方程和 制方程进行求解，从而可预测流场的流动。与流体实 Euler方程 并无实质差别，本文使用的求解器使用 验相比，数值计算 (CFD)周期短，投入少，而且还 Euler方程。积分形式的可压缩二维 Euler方程为 ： 可以模拟许多实验技术无法实现的流动工况。数值计 算与流动实验及理论分析相互配合，是解决流体力学 Wdg2+ 。 问题的很好途径Ⅲ。 流体力学 问题一般使用 Navier．StokeS方程 其中， 为守恒变量组成的一个矢量(密度、密 或 Eular方程来描述，只有在极少数情况下，方程才 度与X方向速度的乘积、密度与Y方向速度的乘积、 有解析解 。而数值模拟是利用计算机对计算区域网 密度与能量的乘积)，F为无粘通量矢量，Q 为控制 体， Q为控制体单元的边界，d~2为体积微元， 为 格化，通过数值求解微分方程来模拟流体运动 。常 面积微元。 用的方法有有限差分法 (FiniteDifferenceMethod1、 有限体积法 (FiniteVolumeMethod)、