讲稿多重网格算法及平均现象的解释.ppt

多重网格算法 对于所以的CFD问题最后都归结为求解线性方程组⑴ AX=d ⑴ 求解线性方程组的方法一般有： 直接方法 迭代法 都比较慢，尤其是在网格数较多时，即使是用松弛因子加速，速度仍较慢。 与其他名词的区别 两层“V”循环多重网格方法的计算步骤 设网格步长分布为H和h 且 H=2h 求解的方程为： Lhuh=fh ⑸ 2、粗网格修正 Lhuh=fh ⑸ 数值模拟中特殊现象的解释及应用 1 、一维波动方程 ⑴ 众所周知他的FTCS差分格式是不稳定的， ⑵ 2 、类似地一维扩散方程 ⑶ 我们熟知的Richardson格式也不稳定， ⑷ 由泰勒展开可得与差分方程等价的偏微分方程，称之为修正的偏微分方程（MPDE） 设原微分方程为 差分方程为 由泰勒展开可得： 所以修正的偏微分方程为： ⑸ 数值耗散余项 数值色散余项 同样的方法可得与LAX格式等价的偏微分方程： ⑺ 由于其耗散主项大于零，所以格式稳定。 所以实质上平均值改变了差分格式的数值耗散余项： ⑻ 刚才我们已经解释FTCS格式通过平均值的替换，就能将逆耗散格式改变为正耗散格式，从而达到稳定，这种做法有时能奏效，但存在很大的盲目性。 * 数值模拟中特殊现象的解释及应用 多重网格算法 一、简介 是求解线性方程组的一种加速收敛的手段 多重网格方法的实质 结构化网格、 非结构化网格 同位网格 交错网格 自适应网格 组合对接、搭接网格、重叠网格 这些术语都是用来定义不同网格特性的术语 多重网格方法的思想基础是来自用Fourier方法对误差的分析。以一维扩散方程为例来分析，其通用差分格式： 二、多重网格算法的原理 可知：误差 也满足方程： ⑵ ⑶ 由Fourier稳定性分析方法可得误差的放大因子G ⑷ 其中： 由 ⑷式可知，放大因子G 在θ∈[0，π]，随θ递减。由此可知不同频率的误差分量在同一网格下的收敛速度不一样，同一频率误差分量在不同网格下的收敛速度也不一样。 如 KmΔx=2π，网格加密一倍 KmΔx=π 所以G 减小 如 KmΔx=π/2，网格放大一倍 KmΔx=π所以G 减小 所以通过这样不断的加密或稀疏网格，可以达到提高收敛速度的目的。 但Δx的改变必须满足稳定性条件 G1 三、多重网格方法的计算步骤 多层“V”循环多重网格方法 两层“V”循环多重网格方法 M M-1 M-2 M-3 M-4 完全多重网格方法 H h 1、设定初值 对Lhuh=fh作1～2次迭代，，得近似解： 两层“V”循环多重网格方法 h H b.从细网格到粗网格转移亏损量： c.在粗网格上精确求解修正量： d.由粗网格到细网格转移修正量： e.计算细网格修正后的量： 多次重复1、2过程直至结果收敛。 a.计算细网格上的亏损量： 多层“V”循环多重网格方法 两层“V”循环多重网格方法 M M-1 M-2 M-3 M-4 完全多重网格方法 一、平均值现象 但当用平均值来代替时， 则变成LAX格式时，该格式却是稳定的。 i-1 i+1 n+1 n-1 同样用算术平均值来代替时，则变成Du