数值模拟课件043.2差分方程组的建立.ppt

差分方程组的建立 主讲教师：李亚军 《油藏数值模拟方法与应用》第三章 * 第一节 基本有限差分 第二节 差分方程组的建立☆ 第三节 网格排列方式及其系数矩阵 第四节 差分方程组的稳定性分析 第五节 边界条件的处理 第三章 差分方程组的建立 显式差分格式 隐式差分格式 Crank-Nicolson格式 * 小结 显式差分格式：计算过程简单 隐式差分格式：计算过程复杂 显式、隐式差分的截断误差相同： 提示：能否提高时间方向上的精度？ 一、一维单相等距网格的差分格式 * 3. Crank-Nicolson差分格式 一、一维单相等距网格的差分格式 一维单相问题 第 i 个网格在第 n＋1/2 时刻, 即(i, n+1/2)点处： 二阶差商 一阶中心差商 * 一、一维单相等距网格的差分格式 3. Crank-Nicolson差分格式 * 图形表示： 六点差分格式 求解精度高 求解工作量与隐式格式相似 一、一维单相等距网格的差分格式 3. Crank-Nicolson差分格式 * 第二节 差分方程组的建立 * 二、二维单相等距网格的差分格式 一、一维单相等距网格的差分格式 三、其它类型的差分格式 1. 显式差分格式 对二维单相问题 ： 对第 ( i , j ) 个网格，在第 n 时刻进行 的二阶差商和 的一阶前差商。 二、二维单相等距网格的差分格式 * 二、二维单相等距网格的差分格式 1. 显式差分格式 * 图形表示： n+1时刻 n时刻 i-1,j i,j i+1,j 已知值 未知值 i,j-1 i,j+1 i,j 二、二维单相等距网格的差分格式 1. 显式差分格式 * 对二维单相问题 ： 对第 ( i , j ) 个网格，在第 n + 1 时刻进行 的二阶差商和 的一阶后差商。 二、二维单相等距网格的差分格式 2. 隐式差分格式 * 二、二维单相等距网格的差分格式 2. 隐式差分格式 * 二、二维单相等距网格的差分格式 2. 隐式差分格式 系数位置 * 图形表示： 已知值 未知值 二、二维单相等距网格的差分格式 2. 隐式差分格式 n+1时刻 n时刻 i-1,j i,j i+1,j i,j-1 i,j+1 i,j * 求解过程 若有n个网格，可写出n个方程，每个方程中有五个未知量，须将方程联立求解。 系数矩阵 五对角带状矩阵 对角占优 稀疏矩阵 对称矩阵 二、二维单相等距网格的差分格式 2. 隐式差分格式 * 3. Crank-Nicolson差分格式 二、二维单相等距网格的差分格式 对二维单相问题 ： 对第 ( i , j ) 个网格，在第 n + 1/2 时刻进行 的二阶差商和 的一阶中心差商。 * 二、二维单相等距网格的差分格式 3. Crank-Nicolson差分格式 * 图形表示： 已知值 未知值 二、二维单相等距网格的差分格式 3. Crank-Nicolson差分格式 n+1时刻 n时刻 i-1,j i,j i+1,j i,j-1 i,j+1 i,j * 第二节 差分方程组的建立 * 二、二维单相等距网格的差分格式 一、一维单相等距网格的差分格式 三、其它类型的差分格式 1、时间中心显式差分格式 图形表示： n-1时刻 n时刻 n+1时刻 i-1 i i+1 三、其它差分格式 * 2、Dufort-Frankel差分格式 图形表示： n-1时刻 n时刻 n+1时刻 i-1 i i+1 三、其它差分格式 * 总结： 1、显式差分格式：一个方程求出一个未知数 2、隐式、C-N差分格式：N个方程组成的方程组来同时求出N个未知量，均属隐式格式 3、隐式差分格式，一维：3对角系数矩阵 二维：5对角 三维：7对角 4、显式差分格式计算简单，但稳定性差；隐式差分计算复杂，但稳定性好 第二节 差分方程组的建立 * * * * * * * * *