第7讲差分方法3.pptVIP

Copyright by Li Xinliang * 将这6个基架点分割成3个组（称为模板） 每个组独立计算 的差分逼近 模板1 模板2 模板3 模板1： {j-3,j-2,j-1,j} 模板2： {j-2,j-1,j,j+1} 模板3： {j-1,j,j+1,j+2} 利用这三个模板的基架点，可构造出逼近 的3阶精度差分格式 计算j点的导数u’, 竟然算出了三个不同的值，怎么办？ ENO 方法： 选择最优（最光滑）的，舍弃其余两个 WENO的处理方法： 三个都要，加权平均它们。 利用Taylor展开式，可唯一确定这些系数）（可利用小程序coeff-schemes.f ) 也可运用优化技术，降低精度、提高分辨率 …… Copyright by Li Xinliang * （3） 对这3个差分值进行加权平均，得到总的差分值 原则： A. 模板内函数越光滑，则权重越大； 模板内有间断时，权重趋于0 B. 三个模拟内函数都光滑时，这三个三阶精度的逼近式可组合成一个五阶精度的逼近式。 “理想权重” （3.1） 确定理想权重 令： 5阶精度 容易解出： Copyright by Li Xinliang * （3.2） 度量每个模板内函数的光滑程度 IS越大，表示越不光滑。 光滑区，不同模板上的IS趋近同一值。 具体形式见下一节。 （3.3） 给出实际权重 构造IS方法很多， 例如： ：第k个模板 光滑区逼近 O（1）量级 间断区 量级，很大 特点： 间断区权重很小 光滑区，趋近于理想权重 （3.4） 给出最终的差分逼近 Copyright by Li Xinliang * 3. Jiang Shu 的五阶WENO格式 守恒型；目 前使用的WENO格式均为守恒型 针对方程： 模板1 模板2 模板3 构造差分格式如下： 构造方法与前文相同 （但注意这里构造的是通量，而前文是直接构造差分格式） 针对整个网格基，构造出5阶精度的通量（理想情况下的通量） 并构造出每个模板上的通量，计算出理想权重。 仍利用程序coeff-schemes.f求系数 理想权重 光滑度量因子 实际权重 Copyright by Li Xinliang * 光滑度量因子的计算 （Jiang Shu） k=1 k=2 k=3 其中： j-2 j-1 j j+1 j+2 是使用模板k 得到的插值函数 利用{j-2,j-1,j}点上的值构造的插值函数 ， 特点： 光滑区趋近同一个值 非光滑区值远大于光滑区 O(1) j 点一阶、二阶导数的差分逼近（用模板k计算） 代入 Copyright by Li Xinliang * 最终5阶 WENO 格式为 正通量情况（ a0） 负通量情况 （a0） 注： 正通量差分格式中下标“j+k”改成“j-k”即得到负通量的差分格式(除了第1式 不变） 注意:是j-1/2 而不是j+1/2 k=1 k=2 k=3 Copyright by Li Xinliang * 4. WENO格式的边界处理 （1） 简易的（降阶）处理方法： 如果某模板用到边界外的点，简单将该模板权重设为0即可 如果用到边界点外的点，则该权重设为0 效果不错，但会边界降阶 （推荐） （2） 构造单边差分的WENO格式 优点： 精度高 缺点： 稳定性不易保证 可能会出现负权重，造成不稳定，可用如下文献的方法处理： Shi J, Hu CQ, and Shu CW, 2002, A Technique of Treating Negative Weights in WENO Schemes, Journal of Computational Physics 175, 108–127 j=1 2 3 4