中国科学院机构知识库网格
Chinese Academy of Sciences Institutional Repositories Grid
二维各向同性湍流直接数值模拟的六边形谱方法及GPU实现和优化

文献类型:期刊论文

作者乔海军 ; 李会元
刊名数值计算与计算机应用
出版日期2013
卷号34期号:2页码:147-160
ISSN号1000-3266
关键词GPU CUDA 六边形傅里叶谱方法 二维各向同性湍流 直接数值模拟
其他题名hexagonal spectral methods for direct numerical simulation of two-dimensional homogeneous isotropic turbulence and their gpu implementation and optimization
中文摘要本文在六边形傅里叶分析及六边形快速傅里叶变换的基础上,提出了二维各向同性湍流直接数值模拟的对偶六边形傅里叶谱方法和六边形傅里叶谱方法,基于二维Navier-Stokes方程的涡度-速度形式,构造了两种六边形傅里叶谱方法的离散格式,设计了其快速求解算法,并且在GPU高性能平台上研制并优化了相应的数值模拟程序.根据方程的具体形式和六边形傅里叶谱方法的特点,从算法层面对方程的求解过程,尤其是非线性Jacobian项快速计算进行优化,经过优化之后,方程求解算法的计算复杂度减少了约30%;根据GPU的体系结构和数值模拟程序的功能要求,将计算模块全部设计为在GPU上运行的kernel函数,尽量避免内存与显存之间的数据拷贝,并在软件工程层面上对代码进行性能优化.优化后的GPU程序获得了高达50倍的加速比.在此基础上,我们对二维各向同性湍流进行了初步的数值模拟,并考察了在不同初始雷诺数条件下,能量和拟涡能随着时间的演变曲线.计算结果表明六边形傅里叶谱方法与传统的傅里叶谱方法一样高效精确.
英文摘要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.
学科主题Mechanics ; Computer Science (provided by Thomson Reuters)
资助信息国家自然科学基金项目(No.91130014)资助|广东省计算科学重点实验室和广东省引进创新科研团队计算科学科研团队支持
收录类别CNKI ; WANFANG ; CSCD
语种中文
CSCD记录号CSCD:4852089
公开日期2013-09-17
源URL[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.

入库方式: OAI收割

来源:软件研究所

浏览0
下载0
收藏0
其他版本

除非特别说明,本系统中所有内容都受版权保护,并保留所有权利。