论文部分内容阅读
介绍Krylov子空间迭代算法及预处理方法。采用有限体积方法对不可压缩流动方程进行离散。对离散形成的大型代数方程组，采用ILUT为预处理的BiCGstab算法进行求解。给出ILUT预处理算法及BiCGstab算法求解代数方程组的步骤。将该算法应用于旋转叶轮内流计算，计算结果与ERCOFTAC叶轮已有的结果较为符合。数值计算表明采用ILUT配合BiCGstab算法比选择标准ILU预处理速度更快，稳定性也较好。