计算机应用 ›› 2013, Vol. 33 ›› Issue (10): 2783-2786.

• 先进计算 • 上一篇    下一篇

交替方向隐式CFD解法器的GPU并行计算及其优化

邓亮,徐传福,刘巍,张理论   

  1. 国防科学技术大学 计算机学院,长沙 410073
  • 收稿日期:2013-04-24 修回日期:2013-06-13 出版日期:2013-10-01 发布日期:2013-11-01
  • 通讯作者: 邓亮
  • 作者简介:邓亮(1989-),男,湖南娄底人,硕士研究生,CCF会员,主要研究方向:高性能计算、GPU并行加速、计算流体力学; 徐传福(1980-),男,安徽六安人,助理研究员,博士,主要研究方向:高性能计算、并行计算机性能评测;刘巍(1980-),男,吉林人,博士,主要研究方向:计算流体力学、高性能计算;张理论(1975-),男,河南南阳人,研究员,博士,主要研究方向:高性能计算。
  • 基金资助:
    家自然科学基金资助项目;空气动力学国家重点实验室开放课题

Parallelization and optimization of alternating direction implicit CFD solver on GPU

Liang DENG,XU Chuanfu,LIU Wei,ZHANG Lilun   

  1. School of Computer Science, National University of Defense Technology, Changsha Hunan 410073, China
  • Received:2013-04-24 Revised:2013-06-13 Online:2013-11-01 Published:2013-10-01
  • Contact: Liang DENG

摘要: 交替方向隐格式(ADI)是常见的偏微分方程离散格式之一,目前对ADI格式在计算流体力学(CFD)实际应用中的GPU并行工作开展较少。从一个有限体积CFD应用出发,通过分析ADI解法器的特点和计算流程,基于统一计算架构(CUDA)编程模型设计了基于网格点与网格线的两类细粒度GPU并行算法,讨论了若干性能优化方法。在天河-1A系统上,采用128×128×128网格规模的单区结构网格算例,无粘项、粘性项及ADI迭代计算的GPU并行性能相对于单CPU核,分别取得了100.1、40.1和10.3倍的加速比,整体ADI CFD解法器的GPU并行加速比为17.3

关键词: 交替方向隐式, 计算流体力学解法器, 结构网格, GPU并行

Abstract: Alternating Direction Implicit (ADI) scheme is a typical discretization scheme for solving partial differential equations. However, there are few researches on the implementations and optimizations of ADI scheme on GPUs for practical Computational Fluid Dynamics (CFD) applications. In this paper, through analysis of the characteristics and calculation processes of ADI solver in a practical CFD application, the authors implemented fine-grained GPU parallelization algorithm for the ADI solver based on grid points and grid lines by a Compute Unified Device Architecture (CUDA) model. Some performance optimization methods were discussed. The experimental results on the TianHe-1A supercomputer show that the proposed GPU-enabled ADI solver can achieve overall speedup of 17.3 compared to single CPU core when simulating a 128×128×128 grid. The speedups for inviscid flux calculation, viscous flux calculation and ADI iteration are 100.1, 40.1 and 10.3 respectively.

Key words: Alternating Direction Implicit (ADI), Computational Fluid Dynamics (CFD) solver, structured grid, GPU parallelization

中图分类号: