块三对角线性代数方程组的gpu加速求解-nvidia.pdfVIP

块三对角线性代数方程组的GPU 加速求解 一、研究背景 计算流体力学（Computational Fluid Dynamics，简称CFD ）是21 世纪流体力学领域的 重要技术之一，CFD 相当于 “虚拟”地在计算机做实验，用以模拟仿真实际的流体流动情 况。而其基本原理则是数值求解控制流体流动的微分方程，得出流体流动的流场在连续区域 上的离散分布，从而近似模拟流体流动情况。 随着计算机技术的推广普及和计算方法的新发展，几十年来CFD 技术取得了蓬勃的发 展。由于数值模拟相对于实验研究有很独特的优点，比如成本低，周期短，能获得完整的数 据，能模拟出实际运行过程中各种所测数据状态，对于设计、改造等商业或实验室应用起到 重要的指导作用。近年来，随着计算机容量的大大提高、先进的数值计算技术的出现以及各 种湍流模型的提出，CFD 技术已经广泛应用于航空航天、化工、能源、环境、水利工程、 气象预测等领域。美国海空军下一代F-35 战斗机所使用的附面层分离进气道是CFD 的成果 之一。附面层分离进气道通过特殊设计形状的突起分离流速较慢的附面层以改善涡轮风扇发 动机的进气流场。此设计比传统的附面层隔板方法可以减轻数百公斤重量，同时在一定速度 范围内能够维持很好的分离效率。 二、数学问题描述 基于连续介质假设的计算流体力学控制方程是非线性偏微分方程组，必须进行数值求 解。其本质是偏微分方程经由数值方法离散成一个线性代数系统的求解。在实际的工程背景 下，为了避免离散得到的大型稀疏矩阵所增加的存储开销，往往通过对矩阵元素的重新排列， 将其重组为块状或带状矩阵。以二维系统为例，在如图 1 所示的网格化计算区域上（假设x 方向有m 个节点，y 方向有n 个节点） N W E P S 图1 计算区域网格化 运用偏微分方程离散化方法，通用变量 在任一点P 处将满足如下离散方程  c  c  c  c  c  d P P E E W W S S N N P 其中，E 、W、S 、N 分别表示与P 点相邻的东、西、南、北计算节点，c 和d 分别表示系数 和常数项。 将待求变量沿y 方向逐行存储，对应的未知向量可表示为 X T (x T , x T , , x T ) 1 2 n 其中 x T ( , , , ) j 1,j 2,j m ,j 对应的离散方程组系数矩阵 B1 C1 0 0 0 0 0    A B C 0 0  2 2 2 