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

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

1、§2-3、差分格式稳定性分析法、一、离散扰动的稳定性分析理论基础；误差传播特性；归纳法例；FTCS格式设初始扰动（误差）为：则为了使则要求及再计算当，误差分布可能出现的情况如下；m-1mm＋1m-1mm＋1（1）（2）m-1mm＋1（3）以此误差分布，从n时间步计算到时间步，并要求在时间层上：这个最后条件，如第一步计算中，附加误差修正不过冲的条件即。但“不过冲”与误差传播振幅不扩大的含义并不一致。一、矩阵法（谱分析法）一个更严格的关于初值问题差分格式稳定性分析的方法是矩阵法。解是一个解向量，经过后解的改变由A是一个变换矩阵，分析变换矩阵的性质，讨论解的性

2、质变化。（仍讨论上述的例子）但为了满足的最大值，经矩阵相乘后，其幅度不增大，从线性代数的分析可知，其必要条件是：矩阵的谱半经不大于1；（谱半径的定义是，矩阵的所有特征值中的绝对值的最大值。）一致有界的充要条件是A的谱半经设：A的特征值为；相应的特征向量；（m=1,2,···，J-1）由此递推;l线性代数中关于求三对角矩阵特征值的定理：设矩阵A为M阶的三对角矩阵，即：则A的特征值为利用该定理，A的特征值是：三,Von.Neumann稳定性分析性,基本思想:分析差分数值解的耗散特性,判断数值解的是否有界的特性；(或曰差分方程对误差的传播性质)初始解(初始误差

3、)利用Fourier展开成Fourier级数线性问题中,通过分析任意一个Fourier分量解的性质特性;数值解有界（或初始误差在传播过程中不扩大），则要求层间放大因子（放大矩阵）要求解有界：所以，放大矩阵（因子）一致有界性若：G为复数，则：G为矩阵此处引入的层间放大因子故可设等。例1；Lax-Wendroff格式由即：容易解出是的必要条件，即Von-Neumann分析的稳定性条件是。例2、FTCS格式；设:代入整理，并求得到2s10在复平面上是表示实半轴为2S，虚半轴为C的椭圆要求即所有椭圆上的点均应在单位园之内。显然应：另对于椭圆：其顶点的曲率半经分别

4、为；故还应有：若引入网格雷诺数的概念则由：例：则则*此处提请注意而不是有些参考书上所说的在运用此格式求解Burger’s方程时,对于数值解是否出现振荡的问题是关键的，伪振荡显然是不希望;但在，的条件下，FTCS格式将产生伪振荡。伪振荡并非数值计算的不稳定！其实质是格式的非单调性所引起的。若将格式写成：当，格式是单调的（见单调格式一样）而时，格式是非单调的。可能出现数值解的伪振荡。例：定解条件如下：初值：边值：简单地，选用11个网格布局，第一个时间步时，仅有第10个网格点地值为非零。因：则例c=0.4s=0.1继续计算Thewiggleswillevent

5、uallypropagatetotheotherboundarybutwillremainboundedthroughouttheiterationtosteadystate.Theoscillationsthatoccurinthiscasearesimilartotheoscillationswhichappearwhenasecond-order(orhigher)schemeinusedtosolvetheinviscidBurgersequationforapropagatingdiscontinuity.例3,方程组问题的稳定性分析Lax-W

6、endroff格式：写成向量稳定条件：即求G的特征值时若为稳定性条件例4；非线性问题的局部线化稳定性分析非线性方程差分方程可由守恒或非守恒方程出发写出；等价守恒方程为预测步校正步：或直接由非守恒方程出发，采用MacCormack格式；stepP；stepC；由于Von-Neumawm稳定性分析只能针对线性问题，上述格式中或必须将它当成常数才能进行分析；为此取等于中绝对值较大的一个，并保留原符号；并类似地记：（C可正可负），将非线性问题局部线化（或称时间上冻结非线性项），则稳定分性分析的问题成为；（改写成三步形式）stepp:stepc:step3:ste

7、pp用Von-Neumann分析，得stepcstep3无稳定性问题.所以：故若有则若分别考虑step1和step2的稳定要求:由,采用复平面上与单位园的对比条件可知要求(不失一般性，令c>0);实半轴：虚半轴：端点曲率比较：归纳为由复平面上的要求是；综合，不等式右侧部分较严格的要求是，左侧部分较严格的要求是；例5、关于MacCormack格式稳定性条件更精确的讨论在计算条件下满足,为此将G1G2的关系代入G记：进一步计算可得其中的充分必要条件是(1).讨论必要条件；由以上要求，则至少必须在两端点上满足：，分别讨论：（符合要求）上式方括号中，记（并求根）

8、则有上式二次项的系数为正，所以二次抛物线的开口向上，（如图）X1X2从图中可看出