加速收敛方法概述.docVIP

加速收敛方法概述 1 当地时间步长法 原理： 根据稳定性条件，对于方程的显示时间推进必须遵循Courant条件。为了增加时间步长，提高收敛效率，采用当地时间步长。当地时间步长方法就是在时间推进求解每个网格上的数值解时，采用该网格单元满足稳定性条件的最小时间步长，而不是整个计算域内的所有结点都满足稳定性条件的最小时间步长，这可以大大减少计算量。 UIJt=R 原来时间步长：?t~minCFL?xλ* 受制于最小空间步长。边界层近壁空间网格 当地时间步长：每点采用不同的时间步长推进 U 数值方法： S表示该网格单元面积，n表示该网格单元的边法向，c表示当地声速，CFL为Courant数。 适用性： 对定常问题，收敛后不影响计算精度，可大幅加速收敛。 2多重网格方法 多重网格时近几十年来发展起来的一种加速收敛方法。它先被用于加速收敛椭圆型问题，该方法能够使得迭代矩阵的谱半径与网格间距无关。随着计算流体力学的发展，多重网格在求解欧拉方程、N-S方程过程中得到了应用：Jamson等人首先将其运用到中心差分格式中，并结合runge-kutta法加速了收敛：D.J Marviplis等人将其运用到非结构网格中，并取得了较好的效果。 多重网格算法的基本思想是引入一系列连续变粗的网格，并将其计算流场发展的部分任务转移到粗网格上进行。细网格上的低频误差在粗网格上相当于高频误差，因此用一种消除高频误差的有效方法，在各自的网格上消除相对于该网格的高频误差，但对细网格而言，就消除了一系列频率的误差。这样做的目的有两个好处： （1）在粗网格上推进一步所需要的时间要少得多，工作量小，提高计算效率； （2）在粗网格上空间步长大，推动了解的快速发展，从而使得迭代较少的步数就可能将误差推到计算域外。 这两个优点皆能加速流场的发展，使得迭代的步数较少，减少了工作量。具体的执行是在粗网格与密网格上交替进行。 多重网格的缺点之一是难以表达粗网格与密网格之间的几何关系。尤其对非结构网格，由于非结构网格本身的无序性，这个缺点尤为突出。D.J Marviplis等人提出了叠加式及嵌套式的多重网格，粗网格与密网格无关，各自用非结构网格生成法生成，然后把它们叠加起来，找出它们的投影关系，用投影关系来进行流场参数的传递。该方法需各自生成连续变粗的网格，对网格的生成方法要求较高：构造能充分表达其几何结构特征的粗网格将是很困难的，粗网格往往不能表达复制几何结构的某些关键特征。D.J Marviplis 提出了聚合式的多重网格法，粗网格的生成是对密网格单元进行聚合而得到，该方法的显著优点是只需要一套网格，粗网格是在此基础上通过聚合而得到。其缺点是对解算器的要求较高。