地下水动力学课件.pptVIP

* * 第八章 研究地下水运动的数值模拟方法 求解地下水运动问题的方法有解析法、物理模拟法和数值法。 §8-1 概述 解析解物理概念清楚、计算步骤简便，便于分析各种影响因素之间相互联系又相互制约的内在规律及对地下水运动状态的影响，因此在生产实践中得到广范应用。 解析方法是用数学上的积分方法或积分变换等方法求得数学模型的解析表达式，通常称为解析解或精确解。解析解能够给出所定义的渗流区域内任意点任意时刻的水头值。 1.解析法： 解析法只适用含水层几何形状规则、基本方程简单、定解条件单一的情况。 2.数值法 数值法是把刻划地下水运动的数学模型离散化，把定界问题化成代数方程，解出渗流区域内有限个结点上的数值解（近似解）。 数值法适用性广（复杂的含水层、定解条件等），通用性强（计算机模拟）、并可程序化，修改模型方便。目前，在求解大型地下水流问题时被广泛应用。 数值法模拟计算过程没有物理模拟法逼真、直观，而且计算工作量大，需要借助于计算机进行模拟计算。 数值法中，最常用的是有限差分法（FDM）和有限单元法（FEM)。有限差分法是建立在用差商代替导数的基础上；而有限单元法是建立在直接求函数的近似解的基础上。 占用计算内 存大，计工 作量大。 1.计算程序的通用性强； 2.对不规则边界处理方便； 3.计算单元划分灵活； 4.水流问题、物质输运问题解的精度 一般较FDM的精度高。 有限单 元 法 对自然边界处理的灵活性较差。 1.简单问题的数学表达式和计算的执 行过程比较直观、易懂； 2.算法效率比较高； 3.计算精度高； 4.有可使用的商用软件。 有限差 分 法 缺 点 优 点 数值法 有限差分法和有限单元法优缺点比较 §8-2 有限差分法 有限差分法的基本思想 把渗流区域按一定的方式剖分成许多小区域 （均衡域），用该区域中心点（结点）的集 集合代替连续的渗流区域，在这些点上用差 商近似地代替导数，将描述地下水流问题的 数学模型化为一组以有限个未知函数值为未 知量的差分方程（代数方程）组，通过求 解差分方程组，得到所求解在结点上的近似 值。 用差商近似代替导数 差分的概念 一、承压水一维非稳定流数学模型的概化 图8-1 河间地块承压含水层示意图 0 设在两条平行的河流之间有一均质、各向同性、等厚、无越流补给的承压含水层； 两河流间距离为L，两河水位的变化规律为Φ0(t)和Φl(t) ； 初始时刻含水层内的水头分布为：H0(x)(0≤x≤L)。 在上述条件下，地下水有由高水位一边向低水位一边流动，构成承压水一维非稳定流问题。取如上图所示的坐标系，则该问题的数学模型为： 由于导数有不同的差分格式，所以差分方程也有不同的形式。 0≤x≤L, 0≤t≤Tsum (8-1) 0≤x≤L (8-2) 0≤t≤Tsum (8-3) 0≤t≤Tsum (8-4) Ⅰ 1. 显式差分方程的建立 空间离散：首先将研究区域[0，L]用直线等分为l份，空间步长△x=L/l; 时间离散：将时间段[0，Tsum]用直线等分为m份，时间步长△t= Tsum /m； 二、一维显式有限差分格式 （1） 离散化 图8-2 空间、时间离散示意图 l m 河流 河流 下面以典型结点（i，k) 为例，说明有显式限差分方程建立的思路、过程及求解方法。 如离散网格图所示，任一结点的坐标为： xi=i△x (i=0,1,2…,l)； tk=k△t (k=0,1,2…,m)。 简记为（i，k)，并以Hik表示H(xi， tk）=H( i△x，k△t）,用hik 表示原方程的近似解，即差分方程的解。 在t=tk时刻，（8-1）式左端用二阶导数的近似表达式，右端用时间导数的向前差分近似代替导数，有： (2) 微分方程的差分化 略去0(△t )和 0 (△x)2 ，可得（8-5）式的对应的差分方程： （8-5） （8-6） （8-7）式表明：只要知道了k时段初始时刻tk各结点的hik值，便可计算出k时段末了时刻tk+1的hik+1值（l≤i ≤l-1, 1≤k ≤m-1)，各方程可独立求解，因此，这种方程称为显式有限差分方程。 （8-7） 若定义 ，则（8-6）式可变为： 上式也可变为： (3) 定解条件离散化 在t=0的初始时刻，各结点水头值由初始条件给出： （8-8） 边界结点上各时刻的水头值由边界条件给出： （8-9） 差分方程