推导7阶精度的WENO格式.PPTVIP

* * Copyright by Li Xinliang Copyright by Li Xinliang 计算流体力学讲义 第八讲 差分方法（4） 李新亮 lixl@ ；力学所主楼219； 知识点： 高精度激波捕捉格式——WENO 激波捕捉方法—— Godnov,MUSCL 常用的隐式处理方法——LU-SGS * Copyright by Li Xinliang 讲义、课件上传至 (流体中文网） - “流体论坛” -“ CFD基础理” 讲课录像及课件上传至网盘 /browse.aspx/.Public Copyright by Li Xinliang * 知识回顾 1. Roe 格式—— 准确特征方向的守恒型格式 标量方程的情况： 等价于 标准的1阶迎风差分格式 f(u)的平均增长率 平均斜率 双曲守恒方程组的Roe 格式 平均增长率矩阵 平均斜率 Roe 平均 含义： 点 处的瞬时增长率刚好为区间 上的平均增长率 Copyright by Li Xinliang * 2. TVD 格式 概念： 网格Reynolds数 单调格式、保单调格式及TVD格式 Harten定理： 正系数原则 TVD 保单调 单调 TVD格式= 1阶迎风+ j *（LW格式-1阶迎风） 二阶精度区 TVD区 二阶精度TVD区（二者交集） Copyright by Li Xinliang * § 8.1 WENO格式 —— 高精度的激波捕捉法 1. 基本思路 {j-3, j-2,j-1,j,j+1,j+2} {j-3,j-2,j-1,j}； {j-2,j-1,j,j+1}； {j-1,j,j+1,j+2} 五个基架点被分成三个组 1） 若高精度逼近 ， 必然利用多个基架点 2） 如果该基架点内函数有间断，会导致振荡 3） 间断不可能处处存在 4） 把基架点分成多个组（模板）， 每个模板独立计算j点导数的逼近。 —— 得到多个差分 5）根据每个模板的光滑程度，设定权重 6） 对多个差分结果进行加权平均 。光滑度越高，权重越大。如果某模板存在间断，则权重趋于0； 如果都光滑，则组合成更高阶格式。 Copyright by Li Xinliang * 2. WENO格式的原理描述 考虑线性单波方程： 注： 为了简便，以非守恒型形式为例讲授其思路，实际使用时，请采用下一节介绍的守恒形式 （1） 确定网格基架点： 6个点 {j-3 , j-2,j-1,j,j+1,j+2} 构造出该基架点上的目标差分格式 计算 这6个点可构造5阶迎风差分： 该格式为WENO 的“目标”格式， 即， 光滑区WENO 逼近于该格式。 利用Taylor展开，可唯一确定系数 （可利用小程序coeff-schemes.f ) 实际上，还可利用分辨率优化技术，可构造出新的目标格式（降低精度、提高分辨率，见第4讲）。目前大量WENO的优化版做这种工作。 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越大，表示越不光滑。 光滑