| 二维各向同性湍流直接数值模拟的六边形谱方法及GPU实现和优化 |
| 其他题名 | hexagonal spectral methods for direct numerical simulation of two-dimensional homogeneous isotropic turbulence and their gpu implementation and optimization
|
| 乔海军; 李会元
|
| 2013
|
| 发表期刊 | 数值计算与计算机应用
 |
| ISSN | 1000-3266
|
| 卷号 | 34期号:2页码:147-160 |
| 摘要 | 本文在六边形傅里叶分析及六边形快速傅里叶变换的基础上,提出了二维各向同性湍流直接数值模拟的对偶六边形傅里叶谱方法和六边形傅里叶谱方法,基于二维Navier-Stokes方程的涡度-速度形式,构造了两种六边形傅里叶谱方法的离散格式,设计了其快速求解算法,并且在GPU高性能平台上研制并优化了相应的数值模拟程序.根据方程的具体形式和六边形傅里叶谱方法的特点,从算法层面对方程的求解过程,尤其是非线性Jacobian项快速计算进行优化,经过优化之后,方程求解算法的计算复杂度减少了约30%;根据GPU的体系结构和数值模拟程序的功能要求,将计算模块全部设计为在GPU上运行的kernel函数,尽量避免内存与显存之间的数据拷贝,并在软件工程层面上对代码进行性能优化.优化后的GPU程序获得了高达50倍的加速比.在此基础上,我们对二维各向同性湍流进行了初步的数值模拟,并考察了在不同初始雷诺数条件下,能量和拟涡能随着时间的演变曲线.计算结果表明六边形傅里叶谱方法与传统的傅里叶谱方法一样高效精确. |
| 收录类别 | CNKI
; WANFANG
; CSCD
|
| 其他摘要 | In this paper, a hexagonal spectral method and a dual hexagonal spectral method are proposed for the direct numerical simulation (DNS) of the tow-dimensional homogeneous isotropic turbulences on the basis of Fourier analysis and fast Fourier transform (FFT) on hexagons. Taking advantage of the vorticity-velocity form of the two-dimensional Navier-Stokes equations, we develop a discrete scheme of the two hexagonal spectral methods, and design their fast implementation algorithm on the GPU platform. A collection of optimizations are then made in the algorithm level for seeking the numerical solution, in particular, for evaluation the nonlinear Jacobian term. These efforts yield a reduction of 30 percent for the total computational complexity. On the other hand, a number of optimization efforts are made in the software engineering level to improve performance of our programs, such as avoiding data copy between host memory and GPU memory, increasing the parallel granularity, optimizing the storage hierarchy, merging the kernel functions. These steps finally obtain a 50-times speedup or our GPU programs in comparison to the CPU programs. With the help of our high performance GPU programs, a series of experiments are carried out for the DNS of two-dimensional homogeneous isotropic turbulences. The evolution of energy and enstrophy of turbulences with different initial Reynolds numbers are studied. Our experimental results show that the DNS results of both the hexagonal and the dual hexagonal Fourier spectral methods possess the same accuracy and efficiency as the classical rectangular Fourier spectral method. |
| 关键词 | Gpu
Cuda
六边形傅里叶谱方法
二维各向同性湍流
直接数值模拟
|
| 部门归属 | 中国科学院软件研究所;中国科学院研究生院;
|
| 学科领域 | Mechanics
; Computer Science (Provided By Thomson Reuters)
|
| 资助者 | 国家自然科学基金项目(No.91130014)资助|广东省计算科学重点实验室和广东省引进创新科研团队计算科学科研团队支持
|
| 语种 | 中文
|
| CSCD记录号 | CSCD:4852089
|
| 内容类型 | 期刊论文
|
| URI标识 | http://ir.iscas.ac.cn/handle/311060/15605
|
| 专题 | 中国科学院软件研究所
|
推荐引用方式
GB/T 7714 |
乔海军,李会元. 二维各向同性湍流直接数值模拟的六边形谱方法及GPU实现和优化[J]. 数值计算与计算机应用,2013,34(2):147-160.
|
| APA |
乔海军,&李会元.(2013).二维各向同性湍流直接数值模拟的六边形谱方法及GPU实现和优化.数值计算与计算机应用,34(2),147-160.
|
| MLA |
乔海军,et al."二维各向同性湍流直接数值模拟的六边形谱方法及GPU实现和优化".数值计算与计算机应用 34.2(2013):147-160.
|
修改评论