抛物型方程的有限差分方法II.ppt

基本思想：将每个时间步分为两个小时间步，即分数时间步，每个空间维对应一个小时间步，在第i个小时间步，采用第i个空间变量用隐式，另一个空间变量为隐式的格式。 * * * 将n维问题的每个时间步分为n个子步，在第i个子步中将问题视为第i维空间上的一维问来处理。 直观上这可理解为在每个时间步用n个顺序的一维离散扩过程模拟一个n维的扩散过程。 * Douglas格式也是交替方向隐格式，精度及稳定性完全同P-R格式， 同时出现 ，存储量大， 但是可以推广到三维。 隐格式：新的时间层具有多于一个节点 3维Douglas格式 优点: 追赶法+无条件稳定的隐格式 是三层格式，比PR格式多一层存贮， 但是易于推广到高维 Douglas格式 特点: 引入过渡层+替换 二维：利用PR格式存储量少 三维: D格式 3.局部一维格式 利用最简显式的LOD格式 CN型LOD格式 基本思想：将一个时间区间 分成n个子区间，用一种差分 格式来近似代替一维方程 将n维问题的每个时间步分为n个子步，在第i个子步中将问题视为第i维空间上的一维问来处理。 直观上这可理解为在每个时间步用n个顺序的一维离散扩过程模拟一个n维的扩散过程。 CN型LOD 格式改写为 与P-R格式等价 但LOD格式可推广至三维 4 预测-校正法（PC） 分为两步： 时的近似解; b) 再在 用 二阶精度格式计算 a) 先用一阶精度格式给出 预测，LOD隐式 校正，显式 与D格式 等价 5.跳点格式（无条件稳定，二阶） 一维跳点格式可以直接推广到二维 显格式 隐格式 为奇数 为偶数 引入差分算子： 令 两式相减 相邻 时间层 已知 未知 6.三维问题 显 隐 C-N， 无条件稳定， 二阶精度， 但必须使用ADI方法计算 无条件稳定 条件稳定 运用较少 PR 格式 本章介绍了经典的有限差分方法求解抛物型问题的 显格式、隐格式（包括ADI格式） 显格式的优点是格式构造简单，每个分量 可以独立求解； 隐式构造比较复杂，各分量 需要联立求解，优点 是稳定性好，一维问题追赶法有效，而对于高维 问题的ADI格式，可以通过求解一系列具有主对角 占优三对角系数矩阵的线性方程组来高效求解。 课堂练习 P维扩散方程显格式稳定性条件是什么？ 交替方向隐式格式的引入是为了解决什么问题？ 海潮引起初始水位倾斜的滨海 承压含水层地下水位变化的数值模拟 用隐格式离散为： 线性代数方程组采用传统的求解三对角方程组的追赶法，最终将解随时间的变化曲线画出来，可以分析出滨海水位与海潮的一样具有相同周期，但时间上滞后，其水位变化幅度与海岸线的距离具有衰减性，越远水位变化幅度越小。 用向前差分格式，向后差分格式Crank-Nicolson格式来求解，取 为0.1与0.5进行计算。 精确解 显格式数值解， 显格式数值解在t=0.1时的误差 当采用显格式 时，网格比不满足稳定条件， 数值解不收敛 t=0.1时显示格式，隐式格式， CN格式误差大小的比较。 显格式当时间步长减少到原来的1/4， 空间步长减少到原来的1/2， 误差减少到原来的1/4， 收敛速度为 。 。 CN格式当时间步长和空间步长同时减少到原来的1/2， 误差减少到原来的1/4，收敛速度为 对高维问题的计算，不能不考虑交替方向隐式（ADI） 算法， ADI格式具有无条件稳定，求解快速的优点。本例对二维问题， 利用D’yakonov交替方向隐式方法求数值解，并画出误差图。 D’yakonov交替方向隐式方法 记 ，则 当时间和空间步长 同时缩小到原来的1/2， 误差缩小到原来的1/4， 收敛速度为 . （五）多维问题 截断误差 Fourier方法分析稳定性： 引入记号 向前差分格式 P维显格式稳定性条件为 二维或高维情况显格式不合适 转向考虑隐格式 无条件稳定（绝对稳定） 截断误差： 无条件稳定（绝对稳定） 为提高精度，采用Grank-Nicolson格式 一维隐格式：绝对稳定，系数矩阵为三对角矩阵 并可用追赶法求解 高维隐格式：绝对稳定，但系数矩阵非 三对角矩 阵，不能用追赶法求解，计算量大。 显格式：稳定性限制严格 Alternating direct implicit method 隐格式将 问题: 一维模型的隐格式—可用追赶法计算 多维模型的隐格式, 系数矩阵是块三角矩阵, 追赶法失效, 需要采用迭代算法 ? 多维模型