清华大学 计算流体力学讲义 第二章 理论基础(5)

《清华大学 计算流体力学讲义 第二章 理论基础(5)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、§2.4经典差分式介绍。一、双曲型方程经典格式：1、迎风格式，(upstream,windward)M.Eq.:稳定条件耗散及频散特性.2,蛙跳格式(LeapFrog)这是一个无耗散的格式。3.Lax格式M.Eg.:4.Lax-Wendroff格式5、两步Lax-Wendroff格式step1:step2:6,MacCormack显式格式step1:(Predictstep)step2:(Correctstep)作业*分别针对标准的模型方程与非线性方程（例如Eueler方程）讨论Lax-Wendroff,两步L-W格式和MacCormack格式的异同,7,Beam-W

2、arming两步迎风格式(a>0)P:C:该格式合并一步的形式是：;或M,E：稳定条件8,隐式格式;分析与思考：隐式格式的稳定性如何。9,时间中心隐式格式(梯形差分格式)Time-CenteredimplicitMethod(TrapezoidalDifferencingMethod)由Taylor’展开;同时考虑(1)-(2)：对于若采用中心差无条件稳定；但当C很大时解是振荡的10，介绍一种高阶格式，Rusanov(Burstein-Mirin)Mefhgd这是一种三步格式稳定条件：当时无4阶粘性项当时，无5阶频散项。二，抛物型方程经典格式（通过模型方程来描述）1、

3、显式FTCS格式方程：（注意由上述M。Eq。不能直接看出稳定条件,）2．混合(显－隐)格式.时为FTCS稳定性分析;按Von-Neumam稳定性分析方法,可得上式中分母恒大于零解上述不等式；右边；讨论；上式无条件满足无条件稳定稳定条件综合格式之图示如下：1/21/61/2s3,DuFort-Frankel格式、格式；这是绝对不稳定的格式！修正为；D-F格式仍为显式格式讨论稳定性；若则：若则稳定性分析结论是――无条件稳定的。但是对DuFot-Frankel格式的相容性分析,(或推导M.Eg.)有:或由此可见DnFort-Franew格式的相容性是有条件的，即必须有才是相

4、容的，即应是的高阶小量，这对计算是附加了一个严格的要求！为什么？固而对于使用该格式而言，对时间步长仍有限制，虽然不是从稳定性的要求提出的！4,Box(盒式)格式不等距网格划分的Box格式：方程：令由（1）式得；上式两组关系分别代入（2）式，分别得到（3）和（3）’式：将（3）’式之下标增加1、即改为，余类推，得：其中若，x,t两方向为等变长网格划分则5,ADE格式；奇数步：计算中j从小到大偶数步：计算中j从大到小条件：第一类边界条件相容性条件隐式格式生成的三对角短阵的追赶法求解。令代入（*）式或与（）式对比可知；设为第一类边界条件,即已知由此可由式求得若已知；即故由对

5、于为矩阵的方程组问题，可类似导出，（不做除法，改为求逆运算）二，多维问题的差分格式1，ADI方法(Alterating-Direction-ImplicitMethod)基本思想;保持隐式格式的无条件稳定特点的同时,将多维问题仍退化为只求解三对角(块)矩阵的代数方程组问题.i)*Peaceman&RachfordADI格式step1;step2;合并一步考虑,用VoN-Neumanm方法分析,得在二维条件下,是无条件稳定的;但此格式简单推广到三维时,则并非无条件稳定;在三维条件下,Peaceman&Rachford格式应修正为；Step123每一步都只在一个方向上有未

6、知量，仍保持了只求解三对角方程i)DouglasADI方法Douglas格式的相容性讨论由(2)式得;代入(1)式稳定性分析