计算流体力学课件：差分方法3.pptxVIP

差分方法（3）知识点： 激波捕捉格式—— TVD、WENO、MUSCL、NND 1 知识回顾1. 差分格式的分辨率有效网格点数： 一个波长里面的网格点数 （PPW: Point per Wavelength)修正波数2. 群速度控制(GVC) 33. 流通矢量分裂“在二阶迎风与二阶中心格式中选一个” “选接近一阶迎风的”GVC2 格式 4§ 5.1 非物理振荡及TVD格式1. 数值解中的非物理振荡 间断附近非物理振荡的根源理论1：色散误差导致各波传播速度不同 （第4讲）理论2：粘性耗散不足思路： 物理问题—— 有粘； 物理粘性足以克服本身振荡 数值方法错误计算了物理粘性——不足以克服振荡物理问题本身也可能振荡。但如果错误计算物理粘性，则会错误地加剧（或衰减）振荡。1） 非物理振荡的原因分析理论3： 格式不能保单调 5数值实验 二阶中心差分计算域[0,1],网格点201 （Dx=0.005 ）时间步长Dx=0.0005 T=0.1时刻的u分布Re=200Dx=0.005 现象： Dx一定时，减小Reynolds数可抑制振荡 Reynolds数一定时，减小Dx可抑制振荡暗示是某一特征量Re=2000Dx=0.005 Re=2000Dx=0.0005 相同 6对流-扩散方程的特性：nn+1（线性）差分方程：某点的值是上一时刻周围几个点上值的线性组合物理上要求系数 ak 均非负含义： 某处浓度的增加对下一时刻周围浓度的影响为正。j-2 j-1 j j+1 j+2差分方程单调性（无振荡）条件： 差分方程 （1）中的系数非负网格Reynolds数 72） 重要概念： 网格 Reynolds数 以网格尺度度量的Reynolds数含义： 数值振荡—— 流动尺度为网格尺度 网格 Reynolds数小，该尺度的能量被耗散掉—— 不发生振荡jj+1j-1过于苛刻的条件单方向网格点数106， 三维1018 单纯靠物理粘性抑制振荡，网格间距必须足够小，通常难以实现网格足够小：不会发生振荡； 网格小于激波的实际厚度，则不会振荡网格Reynolds数足够小时，物理粘性发挥作用，抑制振荡 83） 人工粘性 物理粘性 足够小才发挥作用， Reynolds数很高时很难做到 思路： 人为增加粘性系数 （添加人工粘性） 抑制振荡优点：方法简便，有抑制振荡效果缺点：改变了物理问题，带来误差湍流、分离流等——对粘性敏感： 非物理解分离流—— 对粘性敏感转捩——对粘性敏感很难计算对粘性敏感的问题改进措施： A: 局部施加人工粘性 B: 高阶人工粘性Von NeumannMacCormack人工粘性系数大梯度区，加大人工粘性光滑区为二阶小量人工粘性项 94) 数值振荡的定量描述—— 总变差对于离散函数{uj} 定义总变差:单调函数振荡函数j=1j=N含义： 反映了振荡的剧烈程度双曲型守恒方程特点： 沿特征线 ， u不变特征线未相交——总变差不变特征线相交—— 总变差减小结论： 单个双曲型方程，总变差不增（Total Variation Diminishing： TVD） 102 概念： 单调格式、保单调格式与TVD格式n时刻： 单调函数j=1j=Nn+1时刻： 仍是单调函数j=1j=N设n时刻 是单调的，如果n+1时刻的解 仍保证单调，则称该格式为保单调格式。保单调格式基本结论： 常系数的单调格式只能是一阶 单调格式必是保单调的； 线性格式，单调与保单调等价格式：如果满足 则称其为单调格式。 单调格式： 单调格式保单调格式：TVD格式 总变差不增TVD保单调单调 113. TVD格式的理论基础—— Harten定理Harten定理：如果差分格式可写成如下形式：且则格式（1）是TVD格式（1）可验证： Roe格式是TVD格式保证“系数非负”含义：“单调格式必是TVD格式” 12例7.2.1：考虑线性单波方程：试讨论如下 Lax-Wendroff 格式二阶中心人工粘性是否满足Harten条件常系数的单调格式—— 只有一阶精度对比条件：不满足Harten 条件 13知识回顾： Lax-Wendroff 格式Taylor展开，写出修正方程时-空二阶精度巧妙添加人工粘性，不但克服了不稳定