人为耗散和通量改正.ppt

9. 差分格式的定性分析 差分算法的差分性质 定性分析方法(1) 定性分析方法(2) 10. 人为耗散和通量改正 人为耗散方法 通量改正法 通量改正方法之小结 * * 精度：差分方程对源方程的逼近误差 相容性：时空步长趋于0时，差分方程的极限为源方程 精度、相容性、稳定性、收敛性、耗散性、色散性、和守恒性等 稳定性：任何初值扰动对差分数值解的影响随时间推移 不再增加（强稳定）或在一段时间内有界（弱稳定） 收敛性：当步长趋于0，差分数值解收敛于源问题的真解 Lax等价定理：对一个适定的初值问题，在满足相容条件的 前提下，稳定性是收敛性的充要条件。 耗散性：差分余项对解产生的耗散效应 色散性：差分余项对解产生的色散效应 守恒性：数值解保持真解所固有的守恒性的程度 1 Taylor展开法 (与求修正方程一样） c0 迎风格式 精度： 满足相容条件： 若c0则恒为负，相应差分格式为逆耗散格式，按Hirt论断，一切逆耗散格式均不稳； 该耗散项为正的条件是 或 耗散性和稳定性： 色散性： 时空步长趋于0时不受任何限制 耗散项： 色散项 2 Fourier法(或称为Von Neumann方法) c0 迎风格式 代入源方程 代入差分方程 其中 当 (长波)，则有 其中 可见，差分解的幅度增长因子为 相应的增长率为 稳定条件 解析解与差分解的辐角差为 由于Ep与 k 有关, 故产生色散效应 Ep0 称为正色散,差分解超前真实解； Ep0 k称为负色散,差分解滞后于真实解 r=1或1/2时无色散 通常只关心稳定条件本身，则直接用下述方法： 代入差分方程，解出放大因子 稳定性条件是放大因子的模小于等于1。 即使满足稳定性条件，当计算进行到一定时间步数时，仍出现不稳定，此源于非线性效应。即使格式自身存在的数值耗散（或叫隐耗散）, 在激波区一般不足以克服波头振荡等计算不稳定现象。故需要在格式中人为加入耗散项。 数值耗散在非激波区会带来过分的平滑，以致于模糊物理量分布，通量改正是为了尽量减小格式的平滑作用，提高精度。 在格式中人为加入耗散项。 目标 所有物理量能平滑地通过激波区 在激波区两侧满足Rankine-Hugoniot关系 激波传播速度不失真，且与跃变条件相容 激波过渡区的厚度越小越好，且与激波强度无关 在激波区外不产生或有很弱的人为耗散 对一维问题， Neumann-Richtmyer人为耗散法应用较广泛，对多维问题，Lapidus人为耗散法效果较好。 1 Neumann-Richtmyer人为耗散法 (1950，JAP, 21, 232) 在动量方程中加入人为粘性压力项 ，其人为粘性系数比例于速度梯度 （若激波在x向传播） 式中?为密度，u为速度，b为可调参数，一般取b=1.5-2 最后激波区的厚度大致为 即4－5个格距 2 Lapidus人为耗散法 (1967, JCP, 2, 154) 对各分量方程同时加入人为耗散项 可由分裂格式差分求解 个别量的方程不该加人为耗散，如磁通量函数 强耗散格式（或引入人为耗散） 过分平滑 弱耗散格式 激波波头的振荡 为兼顾稳定性和精度，有各种办法，其中广泛应用的是 1973, JCP, 11, 38; 1975, JCP, 18, 248; 1976, JCP, 20, 397; 1979, JCP, 31, 335 通量改正法 （Flux-Corrected Transport） 正耗散：从守恒型方程出发，采用强耗散格式，保证无振荡 逆耗散：对差分解在非激波区进行逆耗散，以抵消差分解的耗散误差 方程及其差分格式 它描述通量F的盈亏 导致U改变的输运过程 第一步：用低阶格式 ( 强耗散，如迎风格式或Lax格式 ) 计算通量 第二步：用逆耗散通量 消除耗散 x U 强耗散格式 弱耗散格式 过度的改正将出现波头振荡，故须对逆耗散通量 加以限制 第二步变为 的选取满足 涉及 两步 首先确定 然后确定 方案3 x Ud 方案1 当精确解存在峰或谷时，此法有明显的削波现象 方案2 此法提高 减小 ， 从而降低削波效应 如何改正A 恒等式 i点流入量 i点流出量 均非负数 均非负数 我们可将上式改为充分条件：