双曲型方程的有限差分法III.ppt

1 一阶双曲型方程 (1) Lax-Friedrichs格式： （六）二维问题 此格式是一阶精度的。 下面讨论稳定性： 那么Von Neumann 条件满足，格式稳定。 (2) Lax-Wendroff格式： 又由Taylor 展开，有 (3) 分数步长法 例如：以Lax-Wendroff格式来完成二步法 显式格式稳定性有条件限制，多维的更加严格，因此考虑隐格式。 （4）隐式格式 全隐格式 Crank-Nicolson格式 由于隐格式求解二维问题得到的线性方程组 其系数矩阵为宽带状，因此求解不甚顺利， 解决方案：交替方向隐式（ADI ）格式 （Alternate Direction Implicit ） （5）ADI（Alternate direction implicit ） 交替方向隐式格式 ADI-1： X方向隐格式 Y方向隐格式 等价于 ADI-2： 二维Beam-Warming格式 X方向隐格式 Y方向隐格式 ＣＮ格式变形去掉高阶项 * 无条件稳定 ADI-2： * 每一个子步只对一个方向是隐式的，系数矩阵是主对角占优的三对角矩阵，可以利用追赶法，减少了计算量； (2)每个单步都是各方向隐式差分算子的乘积，保证 了格式具有无条件稳定的性质； (3)单步格式与相应CN格式之差是一个局部截断误差 不低于CN格式的差分算子，从而保证了格式的局 部截断误差仍为 优点： * 则称方程组是双曲型方程组。 如果A,B为实对称阵，则方程组是双曲型方程组， 也称为对称双曲型方程组。 2. 一阶双曲型方程组 下面以Lax-Wendroff格式为例，讨论差分方程： 利用多元Taylor展开，有 故有： 利用Fourier方法可讨论上式的稳定性： 可得增长矩阵： 如果A,B是对称阵， 可以证明Lax-Wendroff格式的稳定性条件是： 为放宽稳定性条件，有许多改进技术。 如：Strang方法 * 将n维问题的每个时间步分为n个子步处理，将第i个子步视为第i维空间的一维问题来处理。 * * 每一个子步只对一个方向是隐式的，系数矩阵是主对角占优的三对角矩阵，减少了计算量，且使算法具有内在的可并行性； 每个单步都是各方向隐式差分算子的乘积，比CN格式的增长因子有模更大的分母，保证了格式在有适度扰动的情况下仍具有无条件稳定的性质； 单步格式与相应CN格式之差是一个局部截断误差不低于CN格式的差分算子，从而保证了格式的局部截断误差仍为Ｏ（ｈ＾２） * 算子分裂：把一个复杂的算子表示或近似表示为若干指定形式的简单算子乘积的组合 * （五）双曲型方程及方程组的初边值问题 构造二阶精度的边界条件 确定a,b,c ？ 得方程组 确定a,b,c 得方程组 右边界 二阶精度的边界条件 2.一阶双曲型方程及方程组的边界条件 对流方程 怎样给边界条件使方程适定，区域为 X=1不能给边界条件 X=0不能给边界条件 （初始条件） 为对角线元素为负的对角阵 为对角线元素为零的对角阵 为对角线元素为正的对角阵 S为A的特征向量的列所构成的矩阵 - 0 + 处边界条件数目等于 负特征值数目 处边界条件数目等于 正特征值数目 零特征值不需给出边界条件 - 0 + 3.一阶双曲型方程及方程组的数值边界处理 三层格式 需增补 * n+1 n p Q 实际上是迎风格式 数值边界条件（人工边界条件） 注：采用插值法构造边界条件要用内插公式， 使用外推方法往往是不行。即要用稳 定的格式构造边界条件. 例如：下面的两个不可用的边界条件 再如 注：改进 例２：考虑微分方程组（半无界问题） 定解条件为： 思考：不能给出v(x,t)在x=0处的边界条件， 否则问题不适定。 采用Lax-wendroff 格式 左边界需要附加边界条件 计算 方法一、从特征形式出发 特征型: 采用迎风格式 方法二、从方程本身出发 已知边界条件 有： 利用第一个方程： 利用第二个方程： 见72页 * 将n维问题的每个时间步分为n个子步处理，将第i个子步视为第i维空间的一维问题来处理。 * * 每一个子步只对一个方向是隐式的，系数矩阵是主对角占优的三对角矩阵，减少了计算量，且使算法具有内在的可并行性； 每个单步都是各方向隐式差分算子的乘积，比CN格式的增长因子有模更大的分母，保证了格式在有适度扰动的情况下仍具有无条件稳定的性质； 单步格式与相应CN格式之差是一个局部截断误差不低于CN格式的差分算子，从而保证了格式的局部截断误差仍为Ｏ（ｈ＾２） * 算子分裂：把一个复杂的算子表示或近似表示为若干指定形式的简单算子乘积的组合 *