第四章离散化的基本方法.ppt

计算流体动力学Computational fluid mechanics 南京工业大学机械与动力工程学院 凌祥 第四章 离散化的基本方法 本章主线路 有限差分基础 差分方程 显式方法和隐式方法 误差与稳定性分析 显式方法和隐式方法 例如在右图中，在网格点3处 可以据此依次求解！ 以上过程就是显式的求解过程 定义如下： 显式方法中每一个差分方程只包含一个未知数，从而这个未知数可以用直接计算的方式显示求解 显式方法和隐式方法 事实上，热传导方程不仅仅只能写出一种差分格式，如果把刚才右边的空间差分写出第n个时间层的量和第n+1个时间层的量的平均值，即： 克兰克－尼科尔森公式（广泛用于求解抛物型方程问题） 隐式方法 显式方法和隐式方法 在右图中，说明隐式格式： 整理上式（17），得： 为简化，设： 于是，方程转化为： 显式方法和隐式方法 注意式（19）Ki是由第n个时间层上的量构成的，而这些量是已知的，所以Ki也是已知的。 根据式（19），可构造网格点2～6的式子 网格点2 由于T1已知 点3 点4 点6 点5 由于T7已知 显式方法和隐式方法 上面5个方程，有5个未知数，写成矩阵形式： 这是一个三对角矩阵，可以追赶法求解。 显式方法和隐式方法 根据上面的例子可以看出，隐式方法比显示方法复杂的多。 现在考虑非线性方程，比如： 热传导系数是温度的函数。 以上方程对显示法几乎没有影响，此时的差分方程为： 此时方程仍然是线性的。因为只有一个未知数 显式方法和隐式方法 但是，如果对上式采用克兰克－尼科尔森公式 相应的差分方程仍然可以写成式（18）的形式，显然，新得到的差分方程含有未知函数T的乘积。变成了非线性方程组的求解。这项工作极为困难！！！ 这是隐式法一个很大的缺点。通常要采用“线性化”处理，第6章详细讨论这个问题。 显式方法和隐式方法 显式方法和隐式方法区别： 显式方法： 取定Δx， 那Δt就不是独立的，不能任意取值，它受到稳定性条件的限制 而隐式方法： 比显式方法大的多的Δt仍然能保持稳定性，甚至有些隐式方法是无条件稳定的 显式方法和隐式方法 CFD领域，时间推进的方法一般用于下面两类计算： 1、由给定的初始条件得到流场的定常解 时间推进并不要求是时间精确的，只要得到正常的流场。 2、对真正的非定常流，求其时间精确解 时间推进法的时间精度是绝对必要的。 也就是说：在某些情况下，显式方法是最合理的，而在某些情况下，隐式方法才是正确的选则。 为了搞清楚这一点，下面讨论两种方法的主要优缺点 显式方法和隐式方法 两种方法的优缺点： 1、显式方法： 优点：方程建立及编程相对简单 缺点：对取定的Δx， Δt必须小于稳定性条件对他提 出的限制，某些情况下， Δt可能很小才能保持 稳定，这就需要很长的计算时间。 显式方法和隐式方法 显式方法和隐式方法 两种方法的优缺点： 2、隐式方法： 优点：用大的多的Δt也能保持稳定性。 缺点：方法的建立和编程更复杂，每一步计算时间都 比显式法长的多。另外Δt取的越大，截断误差 也越大，此方法在跟踪严格瞬态变化不如显式 法精确 误差与稳定性分析 误差与稳定性分析主要针对显示方法引出。 本节的稳定性分析是针对特定的差分方程的，其结果也仅限于这些特定的方程。 仍然用一维的热传导方程为例： 这个方程的显式差分格式： 那么，稳定性到底是什么？计算为什么会不稳定？ 误差与稳定性分析 很大程度上由数值误差引起的。 计算推进过程中→误差放大→不稳定 计算推进过程中→误差不增加或者衰减→稳定 下一步进行数值误差的讨论。 误差与稳定性分析 考虑刚才的一维的热传导方程，这个方程的数值解 收到两种误差的影响。 （1）离散误差 偏微分方程的精确解（解析解）与相应差分方程的精确解（无舍入误差的解）之间的差。 （2）舍入误差 对数值多次重复产生的数值误差。 误差与稳定性分析 记A=偏微分方程的精确解 D=差分方程的精确解 N=在某个具有有限精度的计算机上实际计算 出的解 离散误差=A-D 舍入误差=ε=A-D 那么N=D+ε ε是舍入误差。本节简称误差。 误差与稳定性分析 将带有舍入误差的解带入方程，得到 而根据定义，D是差分方程的精确解，所以： （20）式-（21）式得到： 误差ε也满足差分方程 误差与稳定性分析 现在，考虑热传导方程的稳定性 假定在求解方程的过程中，从第n步推进