一种基于块雅可比迭代的高阶FR格式隐式方法

摘要：

最近,基于非结构网格的高阶通量重构格式(flux reconstruction,FR)因其构造简单且通用性强而受到越来越多人的关注.但将FR格式应用于大规模复杂流动的模拟时仍面临计算开销大、求解时间长等问题.因此,亟需发展与之相适应的高效隐式求解方法和并行计算技术.本文提出了一种基于块Jacobi迭代的高阶FR格式求解定常二维欧拉方程的单GPU隐式时间推进方法.由于直接求解FR格式空间和隐式时间离散后的全局线性方程组效率低下并且内存占用很大.而通过块雅可比迭代的方式,能够改变全局线性方程组左端矩阵的特征,克服影响求解并行性的相邻单元依赖问题,使得只需要存储和计算对角块矩阵.最终将求解全局线性方程组转化为求解一系列局部单元线性方程组,进而又可利用LU分解法在GPU上并行求解这些小型局部线性方程组.通过二维无黏Bump流动和NACA0012无黏绕流两个数值实验表明,该隐式方法计算收敛所用的迭代步数和计算时间均远小于使用多重网格加速的显式Runge-Kutta格式,且在计算效率方面至少有一个量级的提升.

关键词：高精度格式;通量重构格式;隐式时间推进方法;块雅可比迭代;GPU计算

所属期刊栏目：53

分类号：O355(流体力学)

资助基金：国家重点研发计划资助项目2018YFB1501102

在线出版日期：2021-11-12（万方平台首次上网日期，不代表论文的发表时间）

页数：共13页

页码：1586-1598