双曲抛物型方程.pptVIP

* 第1页，共20页，编辑于2022年，星期五 4 . 1 差分格式 1. 迎风格式 (c0) (c0) 精度： 稳定条件： 2. 蛙跳（leapfrog）格式 精度： 稳定条件： 仅对对流项 * 第2页，共20页，编辑于2022年，星期五 3. 时间前差，空间中心差分的显格式 精度： 稳定条件： 4. Lax-Wendroff格式 精度： 稳定条件： * 第3页，共20页，编辑于2022年，星期五 5. Crank-Nicholson格式 精度： 稳定条件：恒稳 6. 全隐格式 精度： 稳定条件：恒稳 可对对流项及扩散项分别采用不同的格式 * 第4页，共20页，编辑于2022年，星期五 非线性Burgers方程 为N-S方程的一个较好的模型方程 对于特殊选定的初始值和边界条件，及特别的函数f，可得准确解 设 f =u2/2，定解域为 0? x ? L, t ?0, 边界值 u(0, t)=u0, u(L, t)=0 定常问题的准确解： 其中 自行尝试各种格式 * 第5页，共20页，编辑于2022年，星期五 时间微商的差分逼近 上述差分方法中在对时间偏微分时，只分为 1 单步二层格式（前向差分）一阶精度 2 蛙跳格式 （中心差分）： 二阶精度 3 若维持有关欧拉格式的空间微商的差分形式不变（简记为Ld），则我们可引入龙格－库塔(Runge-Kutta)格式，其时间差分精度为4阶（只是内存占用量大） * 第6页，共20页，编辑于2022年，星期五 范例 若对空间微商采用中心差分，即 则龙格－库塔格式为 稳定条件： * 第7页，共20页，编辑于2022年，星期五 守恒型与非守恒型 作为源方程的方程组可以有不同的解析形式，如欧拉形式与拉格朗日形式，守恒形式与非守恒形式，积分形式与特征形式等。不同形式的源方程在解析分析中完全等效，但在数值计算方面则不尽然。 对任一物理量U=U(x,t) ，若能写成 则称U为守恒型变量，F为其通量（密度），形如此方程者为守恒型方程 * 第8页，共20页，编辑于2022年，星期五 质量守恒 动量守恒 能量守恒 磁通守恒 涡旋守恒 * 第9页，共20页，编辑于2022年，星期五 差分格式 从守恒型方程出发设计的格式具有总体守恒的特性，故称为守恒型格式，如该方程可写为 半格点值由插值得到 若体系与外界无交换，则 * 第10页，共20页，编辑于2022年，星期五 差分格式子类 半格点插值公式的选取 若干子类 1. 欧拉显格式 二阶 四阶 2. 欧拉全隐格式 精度： * 第11页，共20页，编辑于2022年，星期五 3. Lax格式 4. 蛙跳－梯形格式 蛙跳步 梯形步 相当于 * 第12页，共20页，编辑于2022年，星期五 5. Lax-Wendroff格式 (等价于半步长Lax＋半步长蛙跳) 格式4、5的时间精度均为二级。 特别是碰到激波时，守恒格式能使激波关系较为精确地满足，因此在激波的计算中应首先考虑使用守恒型格式。 * 第13页，共20页，编辑于2022年，星期五 5. 说明 对流方程的差分格式 耗散性和色散性对于对流方程的差分格式的稳定性具有重要的影响，而理想（磁）流体方程具有与对流方程相似的形式，同属于双曲型方程，故设计格式时须注意耗散性及色散性 在稳定条件满足的情况下，迎风格式、Lax格式及全隐格式具有一阶耗散，属强耗散格式，L-W格式具有三阶耗散，龙格库塔格式具有五阶耗散，属弱耗散格式，色散效应起主导作用。蛙跳、跳点和C-N格式的耗散余项为0，属无耗散格式，在碰到激波等强间断时，弱（无）耗散格式的色散性会导致波头振荡和计算失稳，须引入人为耗散。 * 第14页，共20页，编辑于2022年，星期五