地下水数值模拟03讲解.ppt

第三章 有限差分法 主要内容 有限差分法的基本原理 几种主要的差分格式 二维渗流问题的差分方程 一般差分方程组的解法 一、有限差分法的基本思想 1、基本原理 一、有限差分法的基本思想 2、求解步骤 差分网格 一、有限差分法的基本思想 3、优缺点 二、导数的有限差分近似表示 1、差分的概念 二、导数的有限差分近似表示 2、导数的有限差分形式 差分公式对比 思考与练习 高阶导数的有限差分形式 思考与练习 高阶导数的有限差分形式 三、简单水文地质模型的有限差分方程组 1、水文地质模型 三、简单水文地质模型的有限差分方程组 2、数学模型 三、简单水文地质模型的有限差分方程组 3、有限差分方程 三、简单水文地质模型的有限差分方程组 3、有限差分方程 三、简单水文地质模型的有限差分方程组 3、有限差分方程 水文地质模型描述 网格剖分 以等距剖分为例 将研究区域[0，L]用直线等分为n份，把时间段[0,T]用直线等分成m份 以 表示结点（i,k）处的水头 一、一维显式有限差分格式 显式差分方程的求解 显式差分方程的收敛性和稳定性 二、一维隐式有限差分格式 隐式差分方程的求解 三、一维中心式有限差分格式 小结——三种有限差分格式的对比 * * 第一节 有限差分法的 基本原理 从物理现象引出相应微分方程（方程+边界条件）； 用差分网格离散求解域； 用差分公式将基本方程转化为差分方程（代数方程）； 用差分方程的解作为微分方程的近似解。 选取网格； 对微分方程及定解条件选择差分近似，列出差分格式； 求解差分格式； 讨论差分格式解对于微分方程解的收敛性及误差估计。 Δx，Δy——空间步长 Δt——时间步长 结点 格点 1、对自然边界处理的灵活性较差。 2、对溶质运移等问题，精度受限 1.简单问题的数学表达式和计算的执行过程比较直观、易懂； 2.算法效率比较高，易编程； 缺 点 优 点 设有x的解析函数y(x)，函数y对x的导数为： —— 是函数对自变量的导数，又称微商。 —— 、 分别称为函数及其自变量的差分 ——dy、dx分别是函数及自变量的微分, —— 为函数对自变量的差商。 由导数（微商）和差商的定义可知，当自变量的差分（增量）趋近于零时，就可以由差商得到导数。因此在数值计算中常用差商近似代替导数。 用泰勒级数展开可以推导出导数的有限差分形式。 中心差 后差 前差 截断误差 公式 名称 二阶导数 一阶导数 ? ? ? ? 0 1 5 2 3 4 6 7 8 9 10 11 12 以一维河间地块承压含水层中的水流问题为例。 考虑一个以通过x＝0和x＝L处的长且直的河流为界的承压含水层，该含水层均质各向同性，顶底板水平，上覆弱透水层，垂向补给强度为W，两河流边界的水位分别φ0、 φL为且不随时间变化。 试研究含水层的水头分布。 解析解 （1）网格剖分 沿河流的方向取单宽[0,L]作为研究区域，将L等分为n份，空间步长 对每个结点进行编号，结点编号由左向右依次为0，1，2，…，i，…，n。共有n+1个结点，其中2个边界结点，n-1个内结点。 （2）建立有限差分方程 先任取一结点i进行分析。 移项整理，得： 对于所有内结点1、2、…、n-1，建立结点相应的差分方程组 n-1个线性代数方程，未知量共n-1个 故方程可解。 （3）求解 将方程表示成矩阵形式，则有： 线性代数方程组！ 系数矩阵中的元素仅位于三条对角线上，系数矩阵对称且正定，故称为三对角线性代数方程组。 最有效的求解方法 ——追赶法 第二节 几种主要的差分格式 以一维河间地块承压含水层中的水流问题为例。 含水层均质各向同性，不考虑垂向补给，两河流边界的水位随时间变化，分别φ0(t)、 φL(t) 。 试研究含水层的水头分布。 该问题属于一维承压非稳定流的定解问题。求解该问题，需要对空间和时间进行离散，形成的差分网格称为时空网格。 导数可以利用一阶、二阶导数的差商代替，由于一阶导数可以有三种差商表示，因此分别对水头关于时间的导数项分别运用前差、后差、中心差将得到三种差分格式。 显式有限差分 ← 前差 隐式有限差分 ← 后差 中心式有限差分 ← 中心差 (i,k) (i-1,k) (i+1,k) (i,k+1) 向前差分 = 整理得： 定义 截断误差为：O([Δx]2)+O(Δt) 只要知道了k时段初始时刻tk各结点的hik值，便可计算出k时段末了时刻tk+1的hik+1值，各方程可独立求解，因此，这种方程称为显式有限差分方程。 差分方程