CFD2011-第5讲-差分方法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 Copyright by Li Xinliang * * 计算流体力学讲义2011 第五讲 差分方法（3） 李新亮 lixl@ ；力学所主楼219； 知识点： 激波捕捉格式—— TVD、WENO、MUSCL、NND * Copyright by Li Xinliang 讲义、课件上传至 (流体中文网） - “流体论坛” -“ CFD基础理论 ” 下载地址2： /browse.aspx/.Public 知识回顾 1. 差分格式的分辨率 有效网格点数： 一个波长里面的网格点数 （PPW: Point per Wavelength) 修正波数 2. 群速度控制(GVC) Copyright by Li Xinliang * 3. 流通矢量分裂 “在二阶迎风与二阶中心格式中选一个” “选接近一阶迎风的” GVC2 格式 Copyright by Li Xinliang * § 5.1 非物理振荡及TVD格式 1. 数值解中的非物理振荡 间断附近非物理振荡的根源 理论1：色散误差导致各波传播速度不同 （第4讲） 理论2：粘性耗散不足 思路： 物理问题—— 有粘； 物理粘性足以克服本身振荡 数值方法错误计算了物理粘性——不足以克服振荡 物理问题本身也可能振荡。但如果错误计算物理粘性，则会错误地加剧（或衰减）振荡。 1） 非物理振荡的原因分析 理论3： 格式不能保单调 Copyright by Li Xinliang * 数值实验 二阶中心差分 计算域[0,1], 网格点201 （Dx=0.005 ） 时间步长Dx=0.0005 T=0.1时刻的u分布 Re=200 Dx=0.005 现象： Dx一定时，减小Reynolds数可抑制振荡 Reynolds数一定时，减小Dx可抑制振荡 暗示 是某一特征量 Re=2000 Dx=0.005 Re=2000 Dx=0.0005 相同 Copyright by Li Xinliang * 对流-扩散方程的特性： n n+1 差分方程： 某点的值是上一时刻周围几个点上值的线性组合 物理上要求系数 ak 均非负 含义： 某处浓度的增加对下一时刻周围浓度的影响为正。 j-2 j-1 j j+1 j+2 差分方程单调性（无振荡）条件： 差分方程 （1）中的系数非负 网格Reynolds数 Copyright by Li Xinliang * 2） 重要概念： 网格 Reynolds数 以网格尺度度量的Reynolds数 含义： 数值振荡—— 流动尺度为网格尺度 网格 Reynolds数小，该尺度的能量被耗散掉—— 不发生振荡 j j+1 j-1 过于苛刻的条件 单方向网格点数106， 三维1018 单纯靠物理粘性抑制振荡，网格间距必须足够小，通常难以实现 网格足够小：不会发生振荡； 网格小于激波的实际厚度，则不会振荡 网格Reynolds数足够小时，物理粘性发挥作用，抑制振荡 Copyright by Li Xinliang * 3） 人工粘性 物理粘性 足