摘要

本文开发了Krylov子空间法中的Bi-CGSTAB、GMRES(m)、CGS、TFQMR及QMR方法的计算程序,并将其实施于SIMPLER算法作为其内迭代方法,针对CFD/NHT领域的问题,研究了它们的求解特性;发现: Bi-CGSTAB方法有着高效的收敛速度和良好的稳定性;N-S方程求解中不同方程不同m值的协调选取是GMRES(m)方法在CFD/NHT领域推广应用的关键;CGS和QMR方法易于中断;TFQMR方法收敛速度慢于其他方法,但能适用于更广泛问题的求解。