公开(公告)号：CN110275733B

公开(公告)日：2022-11-22

申请号：CN201910564493.4

申请日：2019-06-27

Applicant: 上海交通大学

Inventor： 文敏华 ,  刘永志 ,  林新华 ,  鲍华 ,  沈泳星 ,  胡跃 ,  王一超 ,  韦建文

IPC: G06F9/38 ,  G06F17/12 ,  G06F17/16

Abstract: 一种基于有限体积法求解声子玻尔兹曼方程的GPU并行加速方法，通过划分非结构网格，确定边界条件和计算参数并初始化能量密度分布后，从CPU内存向GPU显存传输每个网格单元之间的影响系数；然后计算声子散射项，并使用稳定双共轭梯度法(BiCGSTAB)求解线性方程组，对每个网格单元的能量密度分布进行更新并通过GPU对声子模式温度分布和平衡态分布函数进行更新，最后通过比较每个网格单元的能量密度分布更新前后的变化，并当满足收敛条件时停止计算并输出结果。本发明在GPU上并行计算求解过程中的主要迭代部分，CPU负责整个计算过程的数据读取、数据输出以及计算流程控制，从而显著提高了计算效率。

公开(公告)号：CN110275733A

公开(公告)日：2019-09-24

申请号：CN201910564493.4

申请日：2019-06-27

Applicant: 上海交通大学

Inventor： 文敏华 ,  刘永志 ,  林新华 ,  鲍华 ,  沈泳星 ,  胡跃 ,  王一超 ,  韦建文

IPC: G06F9/38 ,  G06F17/12 ,  G06F17/16

Abstract: 一种基于有限体积法求解声子玻尔兹曼方程的GPU并行加速方法，通过划分非结构网格，确定边界条件和计算参数并初始化能量密度分布后，从CPU内存向GPU显存传输每个网格单元之间的影响系数；然后计算声子散射项，并使用稳定双共轭梯度法(BiCGSTAB)求解线性方程组，对每个网格单元的能量密度分布进行更新并通过GPU对声子模式温度分布和平衡态分布函数进行更新，最后通过比较每个网格单元的能量密度分布更新前后的变化，并当满足收敛条件时停止计算并输出结果。本发明在GPU上并行计算求解过程中的主要迭代部分，CPU负责整个计算过程的数据读取、数据输出以及计算流程控制，从而显著提高了计算效率。