地下水数值模拟03-有限差分法ppt课件.ppt

《地下水数值模拟03-有限差分法ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第三章有限差分法主要内容有限差分法的基本原理几种主要的差分格式二维渗流问题的差分方程一般差分方程组的解法第一节有限差分法的基本原理一、有限差分法的基本思想1、基本原理从物理现象引出相应微分方程（方程+边界条件）；用差分网格离散求解域；用差分公式将基本方程转化为差分方程（代数方程）；用差分方程的解作为微分方程的近似解。一、有限差分法的基本思想2、求解步骤选取网格；对微分方程及定解条件选择差分近似，列出差分格式；求解差分格式；讨论差分格式解对于微分方程解的收敛性及误差估计。差分网格Δx，Δy——空间步长Δt——时间步长结点格点一、有限差分法的基本思想3、优缺点优点缺点1.简单问题的数学表达式和计算

2、的执行过程比较直观、易懂；2.算法效率比较高，易编程；1、对自然边界处理的灵活性较差。2、对溶质运移等问题，精度受限二、导数的有限差分近似表示1、差分的概念设有x的解析函数y(x)，函数y对x的导数为：——是函数对自变量的导数，又称微商。——、分别称为函数及其自变量的差分——dy、dx分别是函数及自变量的微分,——为函数对自变量的差商。由导数（微商）和差商的定义可知，当自变量的差分（增量）趋近于零时，就可以由差商得到导数。因此在数值计算中常用差商近似代替导数。二、导数的有限差分近似表示2、导数的有限差分形式用泰勒级数展开可以推导出导数的有限差分形式。差分公式对比名称公式截断误差一阶导数前差后差

3、中心差二阶导数思考与练习高阶导数的有限差分形式???思考与练习高阶导数的有限差分形式?0152346789101112三、简单水文地质模型的有限差分方程组1、水文地质模型以一维河间地块承压含水层中的水流问题为例。考虑一个以通过x＝0和x＝L处的长且直的河流为界的承压含水层，该含水层均质各向同性，顶底板水平，上覆弱透水层，垂向补给强度为W，两河流边界的水位分别φ0、φL为且不随时间变化。试研究含水层的水头分布。三、简单水文地质模型的有限差分方程组2、数学模型解析解三、简单水文地质模型的有限差分方程组3、有限差分方程（1）网格剖分沿河流的方向取单宽[0,L]作为研究区域，将L等分为n份，空间步长对

4、每个结点进行编号，结点编号由左向右依次为0，1，2，…，i，…，n。共有n+1个结点，其中2个边界结点，n-1个内结点。三、简单水文地质模型的有限差分方程组3、有限差分方程（2）建立有限差分方程先任取一结点i进行分析。移项整理，得：对于所有内结点1、2、…、n-1，建立结点相应的差分方程组n-1个线性代数方程，未知量共n-1个故方程可解。三、简单水文地质模型的有限差分方程组3、有限差分方程（3）求解将方程表示成矩阵形式，则有：线性代数方程组！系数矩阵中的元素仅位于三条对角线上，系数矩阵对称且正定，故称为三对角线性代数方程组。最有效的求解方法——追赶法第二节几种主要的差分格式水文地质模型描述以一

5、维河间地块承压含水层中的水流问题为例。含水层均质各向同性，不考虑垂向补给，两河流边界的水位随时间变化，分别φ0(t)、φL(t)。试研究含水层的水头分布。该问题属于一维承压非稳定流的定解问题。求解该问题，需要对空间和时间进行离散，形成的差分网格称为时空网格。网格剖分以等距剖分为例将研究区域[0，L]用直线等分为n份，把时间段[0,T]用直线等分成m份以表示结点（i,k）处的水头导数可以利用一阶、二阶导数的差商代替，由于一阶导数可以有三种差商表示，因此分别对水头关于时间的导数项分别运用前差、后差、中心差将得到三种差分格式。显式有限差分←前差隐式有限差分←后差中心式有限差分←中心差一、一维显式有限

6、差分格式(i,k)(i-1,k)(i+1,k)(i,k+1)向前差分=整理得：定义截断误差为：O([Δx]2)+O(Δt)只要知道了k时段初始时刻tk各结点的hik值，便可计算出k时段末了时刻tk+1的hik+1值，各方程可独立求解，因此，这种方程称为显式有限差分方程。差分方程与离散后的定界条件构成了数学模型的显式差分方程问题。其求解步骤如下：由以上计算的h11，h21…hn-11值及由边界条件计算的h01和hl1,再次利用差分公式（取k=1，i=1，i=2,…i=n-1，便可计算得t2时刻各结点的水头值。如此重复，便可计算出t3，t4，…各时刻的水头分布值。由初始条件给出t0时刻各结点的水头

7、值h00，h10，…hl0；再根据差分方程，在k=0时，分别取i=1，i=2,…i=n-1,便可求得t1时刻各内结点的水头值h11，…hn-10。显式差分方程的求解开始输入初始参数t=0对i=0,1,…,n循环执行t=t+Δt对i=1,…,n-1循环执行输出结果：t，t