递归投影方法(RPM)和其在流场数值模拟中的加速收敛性能研究.pdf

第 28卷 第 1期 空 气 动 力 学 学 报 Vol_28．No．1 2O1O年 2月 ACTA AERODYNAM ICA SINICA Feb．，2O10 文章编号 ：0258—1825(20l0)Ol一0017—06 递归投影方法 (RPM)及其在流场数值 模拟中的加速收敛性能研究 李孝伟，王正裕 (上海大学 ，上海市应用数学和力学研究所 。上海 200072) 摘 要 ：随着计算流体力学中的物理模型的日益复杂和数据运算量的增大 ．人们越来越重视提高数值方法的计算 效率 。本文在研究 RPM方法的基础上，提出了一种更易于保证实施精度的具有 良好加速收敛性能的改型的RPM 方法。一方面．通过构造的一维迭代格式研究了几个控制参数如Krylov接受比、Krylov向量维数等对 RPM方法 的加速收敛性能的影响。另一方面，研究了在流场数值模拟中即使采用了常规 的加速收敛措施如当地时问步长和 隐式残值光顺 以后 ，改型的RPM方法 的加速收敛性能；结果表明，此时．RPM方法仍然能使总体上的收敛速度提 高近 1／3。 关键词 ：RPM 方法 ；不变子空间；当地时间步长 中圉分类号 ：V2l1。3 文献标识码 ：A 2001年，为了提高 RPM 方法对于不变子空间的近似 0 引 言 精度，V．Janovsky和 Liberda．PrahaL3提出一种改 进的RPM 方法 。该方法使用Cayley变形处理雅可 计算流体力学的发展过程是一个流体动力学模 比矩阵，以修正超出稳定边界特征值所对应的子空间 型的精确化和细致化的过程，同时流场数值模拟的收 的正交基 ；在得到的修正 的正交基上，提出一种耦合 敛性和计算效率也成为了计算流体力学 中的关键 问 题 。为了提高计算效率 ，广大计算流体力学工作者发 的双重迭代的RPM方法 ，最后证 明该修正的RPM 展了多种加速收敛的措施 ，典型的如采用当地时间步 方法的能够发现一些原始 RPM 方法发现不 了的 长和隐式残值光顺技术等。但人们并不满足，因为现 Hopf分叉点。2005年 ，M．Dorobantu和 J．M611er 实期待一些更为有效的加速收敛方法的出现。 将RPM 用于计算流体力学的NSMB程序块中作为 1993年 ，G．M．Shroff与 H．B．Kellerc提 出 加速数值收敛进程 的外壳 。]，计算 了二维超音速喷 一 种递归投影算法 (theRecursiveProjectionMeth— 管流动和二维亚音速翼型绕流，并证 明RPM 方法在 od，RPM)，目的是 为了稳定和加速非线性分叉 和混 二维定常流场的计算中可以达到约 5倍的加速效果。 沌的数值研究中的迭代过程。对于离散 的偏微分方 从以上研究者们的结果可 以看 出，RPM 方法的 程，该方法对迭代变量的残差构造Krylov子空间，并 加速收敛效果和作用非常显著，这激发了本文作者对 从 中分离出影响迭代收敛的强势特征子空间，进而通 它的研究兴趣 。在研究中发现 ，原有 RPM 方法在实 过使用牛顿迭代法消除其对迭代收敛进程带来的负 施过程中不易保证正交基的精度，针对此，本文提 出 面影响 。 了一种改型的 RPM 方法。计算实践表明，该改型的 1999年 ，PhilipLove__2在三维的Kolmog