Krylov子空间迭代法在大规模线性方程组求解中应用.pdfVIP

Krylov子空间迭代

亨克·A·乌得勒

支大学

1背景

在连续的模拟中，例如流体通过管道的流动或空气绕飞机的流动，通常会在感的区域上

施加一个网格，并将自己限制在计算相关参数，如、流速或温度，在网格点上的值。定

律导致这些参数在相邻网格点之间的近似关系，结合边界网格点的预定行为和给定的源，最终会

导致非常大的线性方程组Ax=b。向量x表示网格点上未知参数的值，b表示给定的输入，矩

阵A描述了网格点之间参数的关系。由于这些关系通常仅限于附近的网格点，因此矩阵的大多数

元素为零。

当网格变得更加精细，即网格点之间的距离减小时，模型会变得更加准确，而在三维模拟中，

这很容易导致非常大的方程组。即使每个坐标方向只有几百个网格点，也已经会导致数百万个未

知数的系统。许多其他问题也会导致非常大的系统：电路模拟、磁场计算、天气预报、化学过程、

器件模拟、核反应堆安全问题、机械结构中的应力等。

这些线性系统的数值解法基于巧妙实现的消元法，其中尽可能地利用了线性系统的

稀疏结构，以避免与零元素进行计算，并避免