2.3 双曲型方程的差分方法.ppt

2.3.1 一阶线性常系数双曲型方程 　1） 迎风格式（upwind scheme） 　　 迎风格式就是将解沿特征线传播看成是“风”，当风从右边(a0)吹来，迎风作向前差商；当风从左边吹来(a0)，迎风作向后差商，这样构造得差风格式是稳定的。 讨论： 1 )稳定性条件为 ； 2 )截断误差为 ，一般取 所以该格式为一阶精度。 3 )该格式可以不考虑对应的微分方程的特征线的走向，而迎风格式需考虑特征线的走向。为此，我们也可以把迎风格式写成统一的形式： 上两式是左边相同，他们都以 趋近于对流方程。因此这两种格式的截断误差取决于两式的右端项的大小。 由稳定性限制要求有 ， 如果 ，上两式相等。但在实际中总是取 ， 因此Lax-Friedrichs格式的截断误差（Trunction error）TE比迎风格式TE大。 由： 这样，上式就可写为： * 2.3 双曲型方程的差分方法 计 算 力 学 基 础 第二章 有限差分方法 它的解为： 推导： 令： 　　从生活常识上来说，当你突然闻到飘来的花香，你自然会迎风去搜寻，从数学上讲，这就是符合特征线的走向差分离散处理。 该格式是绝对不稳定的差分格式。 1954年，Lax和Friedrichs为克服上述格式的不稳定性，提出了如下的差分格式： 2）Lax-Friedrichs格式 该格式称为Lax-Friedrichs格式，也称为Lax格式。 而Lax-Friedrichs可写为： 实际上，这两种格式还是有很大的区别，现仅从截断误差来考虑，现假设a0，此时迎风格式可写为： 可得： 把上两式代入（B）式，得 (C) （B） 3）Lax-wendroff格式 迎风格式和Lax格式是一阶精度的差分格式，1960年Lax和Wendroff构造出一个二阶精度的格式。 而： 把此两式代入(C)式，得： 略去高阶项，可以得到如下的差分格式： （D） 差分格式（D）称为Lax-wendroff格式。 讨论：1 该格式的稳定性条件为 2 二阶精度显示格式。 我们也可以建两步Lax-wendroff格式，对（D）式进行改造： （E） 差分格式（E）称为两步Lax-wendroff显示格式。 由（C）式 有： 而： 4）MacCormack二步二阶精度的显示格式 把 　和 　　代入（C）式，得： （F） 略去高阶小量，可得如下的差分格式： （G） 若引入： 则有： 将向前、向后差商代入上两式，得： 或者： 5）蛙跳格式 利用双曲型方程的解在特征线为常数这一事实也可以构造出以上差分格式。 为确定起见，不妨设 a0 。 设在 时间层上网格点 和 上的 值给定。现 计算 时间层上的网格点 点的 值。假定 条件成立，过 点特征线与 交于 ，由解的性质可知， ，当 不是 网格点时， 是未知的，但可以用插值方法给出 近似值。 6）利用特征线构造差分格式 两点线性插值： 三点抛物线插值： (1) 利用 两点线性插值 由此可推导出差分格式： 迎风格式 (2) 利用 两点进行线性插值 由此得到： 改写为： Lax-Friedrichs格式 (3) 如果使用 和 三个点进行抛物插值，则得到 由此可得差分格式： Lax-wendroff格式 *