第八章 地下水模拟.docVIP

数学模型一般由控制方程与约束条件组成，控制方程可分为线性方程和非线性方程，约束条件有边界条件（三类边界）和初始条件（稳定流时无）。求解数学模型的解法可分为解析法和数值法两大类。20世纪60年代以前，解含水层地下水的水头和流量问题，多偏重于解析解法，比如《地下水动力学》课程中所讲，无论是以稳定流为基础的裘布衣公式，还是以非稳定流为基础的泰斯公式，它们的推导都有许多假设，在水文地质条件满足这些假设时，当然没有问题，但要解决大范围的地下水系统计算时，由于水文地质条件的复杂性，解析解法就无能为力了。60年代后期随着电子计算机的发展，人们把数值模拟应用到水文地质计算中来。由于电模拟制作和参数调试都比数值法麻烦，所以目前应用更多的是数值解法。 数值法是以离散化方法求解数学模型的。这里的离散化方法是指一种“离散后集合”的方法，即将研究区域和计算时段分别分割成若干个小区域(单元)与小时段(对静态模型不需要)，把各单元视为区域的等价区，根据数学模型建立各单元特征点物理量之间关系的代数方程，按单元、按时段顺次求解由这些代数方程组成的方程组，从而最终得出全区域变量值的集合。 ①有限差分法是通过以有限个离散点代替连续区域，以差商代换离散点上微商的途径，将偏微分方程转化为线性代数方程来求解数学模型；②有限单元法则依据变分原理、剩余加权法或均衡原理进行剖分插值，得出代数方程组，以求解数学模型；边 界元法是近几年开始广泛应用于工程问题上的一种方法，它是通过将区域水流的基本微分方程变换为边界积分方程后离散成线性代数方程组的途径来求解数学模型。