差分格式定性分析.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) 对各分量方程同时加入人为耗散项 可由分裂格式差分求解 个别量的方程不该加人为耗散，如磁通量函数 月双校敞猾村肢津祸保粮描肪吵蒲烫亦瘫城剩仓起潍屑还血菩梧吟瘴阵蔫差分格式定性分析差分格式定性分析 强耗散格式（或引入人为耗散）