讲稿多重网格法及平均现象解释.pptVIP

数值模拟中特殊现象的解释及应用 ;多重网格算法; 对于所以的CFD问题最后都归结为求解线性方程组⑴ AX=d ⑴ 求解线性方程组的方法一般有： 直接方法 迭代法 都比较慢，尤其是在网格数较多时，即使是用松弛因子加速，速度仍较慢。 ;与其他名词的区别;多重网格方法的思想基础是来自用Fourier方法对误差的分析。以一维扩散方程为例来分析，其通用差分格式： ;由Fourier稳定性分析方法可得误差的放大因子G ;如 KmΔx=2π，网格加密一倍 KmΔx=π 所以G 减小 如 KmΔx=π/2，网格放大一倍 KmΔx=π所以G 减小;三、多重网格方法的计算步骤;两层“V”循环多重网格方法的计算步骤 设网格步长分布为H和h 且 H=2h 求解的方程为： Lhuh=fh ⑸ ;2、粗网格修正 Lhuh=fh ⑸ ;多层“V”循环多重网格???法;数值模拟中特殊现象的解释及应用 ;1 、一维波动方程 ⑴ 众所周知他的FTCS差分格式是不稳定的， ⑵ ;2 、类似地一维扩散方程 ⑶ 我们熟知的Richardson格式也不稳定， ⑷ ;二、数学解释;结论1、若数值耗散余项的主项大于零，则该差分格式为正耗散格式，而且格式稳定；反之，为逆耗散格式，且格式不稳定。 ;因此可得与差分格式等价的偏微分方程: ⑹ 因此其耗散余项的主项小于零，所以离散格式为逆耗散格式，不稳定;同样的方法可得与LAX格式等价的偏微分方程： ⑺ 由于其耗散主项大于零，所以格式稳定。 所以实质上平均值改变了差分格式的数值耗散余项： ⑻ ;L-W