计算流体力学CFD课件.ppt

计算流体力学CFD(1) 计算流体力学CFD(2) 计算流体力学CFD(3) 误差与稳定性分析 由于误差?也满足差分方程，故有 而 则 误差与稳定性分析 解得 放大因子 误差与稳定性分析 要使 放大因子 必须满足 误差与稳定性分析 上式就是差分方程 的稳定性条件。 对于给定的?x，?t的值必须足够小，才能满足上述稳定性条件，以保证计算过程中误差不会放大。 误差与稳定性分析 稳定性条件的具体形式取决于差分方程的形式。 的差分方程 是无条件不稳定的。 比如，一阶波动方程： 误差与稳定性分析 但如果用 则 （Lax方法） 误差与稳定性分析 令误差 则放大因子 式中 误差与稳定性分析 则放大因子 稳定性要求 则 误差与稳定性分析 稳定性要求 式中的C称为柯朗(Courant)数。 误差与稳定性分析 稳定性要求 上式称为柯朗－弗里德里奇－列维(Courant-Friedrichs-Lewy)条件，一般写成CFL条件。 误差与稳定性分析 下面来看CFL条件的物理意义。 CFL条件： 也是二阶波动方程： 的稳定性条件。 误差与稳定性分析 下面来看CFL条件的物理意义。 二阶波动方程： 的特征线为 CFL条件的物理意义：要保证稳定性，数值解的依赖区域必须全部包含解析解的依赖区域。 误差与稳定性分析 CFL条件的物理意义：要保证稳定性，数值解的依赖区域必须全部包含解析解的依赖区域。 误差与稳定性分析 显式方法 边界条件已知 显式方法 显式方法中每一个差分方程只包含一个未知数，从而这个未知数可以用直接计算的方法显式地求解。显式方法是最简单的方法。 隐式方法 隐式方法 克兰克－尼科尔森格式 隐式方法 对于排列在同一时间层所有网格点上的未知量，必须将它们联立起来同时求解，才能求出这些未知量，这种方法就定义为隐式方法。 隐式方法 由于需要求解联立的代数方程组，隐式方法通常涉及大型矩阵的运算。隐式方法比显式方法需要更多、更复杂的计算。 隐式方法 隐式方法 A，B，Ki 均为已知量 隐式方法 A，B，Ki 均为已知量 隐式方法 在网格点2： A，B，Ki 均为已知量 T1 为边界条件，已知量 隐式方法 在网格点3： A，B，Ki 均为已知量 在网格点4： 在网格点5： 隐式方法 A，B，Ki 均为已知量 在网格点6： T7 为边界条件，已知量 隐式方法 于是有关于T2，T3，T4，T5，T6这五个未知数的五个方程 A，B，Ki 均为已知量 隐式方法 写成矩阵形式： 隐式方法 系数矩阵是一个三对角矩阵，仅在三条对角线上有非零元素。 求解线性代数方程组的标准方法是高斯消去法。应用于三对角方程组，通常采用托马斯算法（国内称为追赶法）求解。 显式方法与隐式方法的比较 显式方法与隐式方法的比较 对于显式方法，一旦?x取定，那么?t的取值必须受到稳定性条件的限制，其取值必须小于等于某个值。否则，计算不稳定。因此，?t必须取得很小，才能保持计算稳定，要算到某个给定的时间值，程序要运行很长时间。 显式方法与隐式方法的比较 隐式方法没有稳定性限制，可以取比显式方法大得多的?t，仍能保持计算稳定。要计算某个给定的时间值，隐式方法所用的时间步数比显式方法少很多。 显式方法与隐式方法的比较 对某些应用来说，虽然隐式方法一个时间步的计算会比显式方法花的时间长，但由于时间步数少，总的运行时间可能比显式方法少。 显式方法与隐式方法的比较 另外，当?t取得较大时，截断误差就大，隐式方法在跟踪严格的瞬态变化（未知函数随时间的变化）时，可能不如显式方法精确。 不过，对于以定常态为最终目标的时间相关算法，时间上够不够精确并不重要。 显式方法与隐式方法的比较 当流场中某些局部区域的网格点分布很密，采用显式方法，小的时间步长会导致计算时间特别长。 例如，高雷诺数粘性流，物面附近的流场会产生急剧的变化，因此，物面附近需要更密的空间网格。 在这种情况下，若采用隐式方法，即使对于很密的空间网格，也能采用较大的时间步长，就会减少程序运行时间。 误差与稳定性分析 误差与稳定性分析 在从一个推进步进行到下一步时，如果某个特定的数值误差被放大了，那么计算就变成不稳定。如果误差不增长，甚至在从一个推进步进行到下一步时，误差还在衰减，那么计算通常就是稳定的。 误差与稳定性分析 A=偏微分方程的精确解（解析解） D=差分方程的精确解 离散误差=A-D 误差与稳定性分析 D=差分方程的精确解 舍入误差=?=N-D N=在某个有限精度的计算机上实际计算出来的解 （数值解） N=D+? 误差与稳定性分析 数值解N=精确解D+误差? 数值解N满足差分方程，于是有 误差与稳定性分析 数值解N=精确解D+误差? 精确解D也必然满足差分方程，于是有 误差与稳定性分析 数值解N=精确