第二讲数值模拟基础及方法理论.pptVIP

3）显示差分方程的求解计算各结点初始时刻水头值利用差分方程计算各结点t1时刻水头值利用边界条件计算边界结点水头值重复2、3步，直到计算出拟计算的各个时刻的水头值显式格式（续3）第26页，共49页，星期日，2025年，2月5日算例（续4）在上述模型中，设L=1000米取空间步长为200米，时间步长为0.25天，分别计算各节点各时刻的水头值。第27页，共49页，星期日，2025年，2月5日Time/dayx=0mx=200mx=400mx=600mx=800mx=1000m02010101010100.252012.5101010100.502013.7510.6251010100.752014.53111.25010.15610101.002015.07811.79710.39110.039101.252015.48812.26610.65410.11710算例（续5）第28页，共49页，星期日，2025年，2月5日Time/dayx=0mx=200mx=400mx=600mx=800mx=1000m02010101010100.252020101010100.502010201010100.75203002010101.0020-10101.252010算例（续6）如果⊿t=1,则第29页，共49页，星期日，2025年，2月5日(6)算例：隐式格式在上述模型中，设L=1000米取空间步长为200米，时间步长为0.25天，用隐式差分格式计算各节点个时刻的水头值。第30页，共49页，星期日，2025年，2月5日在这个例子中，解：隐式格式一般方程为于是有根据初始条件得根据边界条件得第31页，共49页，星期日，2025年，2月5日由初始条件和边界条件由此解得t1时刻的水头值为在上述方程中取n=0,可以得到计算t1时刻水头值的方程所以上述方程变成同理，可计算t2时刻的水头值第32页，共49页，星期日，2025年，2月5日2.地下水数值模型建模步骤第33页，共49页，星期日，2025年，2月5日模拟步骤建立概念模型建立数学模型数值方法及软件（编程）参数拟合模拟：模型校正与检验参数敏感性分析预测模拟软件第34页，共49页，星期日，2025年，2月5日一、概念模型（模型概化）根据详细的地形地貌、地质、水文地质、构造地质、水文地球化学、岩石矿物、水文、气象、工农业利用情况等模拟的区域:含水层类型:潜水（无压）、承压、混合、多层维数:一维、二维、三维水流状态：稳定流/非稳定流、饱和流/非饱和流介质状况:均质和非均质/各向同性和各向异性孔隙/裂隙/双重介质流体的密度差边界条件和初始条件必要时需进行一系列的室内试验与野外试验,以获取有关参数,如渗透系数、弥散系数、分配系数、反应速率常数等。第35页，共49页，星期日，2025年，2月5日第1页，共49页，星期日，2025年，2月5日地下水建模方法和步骤1.求解地下水运动方程的数值方法2.地下水数值模型建模步骤3.建模所需要的基本资料第2页，共49页，星期日，2025年，2月5日绝大部分数学模型是无法用解析法求解的，数值化就是将数学模型转化为可解的数值模型。1.数值方法1.1有限差分法1.2有限单元法1.3分有限差分法1.4半解析半数值法1.5边界元法第3页，共49页，星期日，2025年，2月5日(1)有限差分法原理(2)两种方法建立有限差分方程(3)求解有限差分方程(4)收敛性和稳定性概念(5)算例1.1有限差分法第4页，共49页，星期日，2025年，2月5日(1)有限差分法的基本原理将连续的问题离散后求解：方法一．以地下水流基本微分方程及其定解条件为基础，在渗流区剖分基础上，用差商代替微商，将地下水流微分方程的求解转化为差分方程（代数方程）求解。方法二．在渗流区剖分的基础上，直接由达西定律和水均衡原理，建立各个均衡区的水均衡方程，即差分方程。矩形网格多边形网格第5页，共49页，星期日，2025年，2月5日网格划分的基本类型（1）先划格线，格点位于网格中心均衡网格节点网格（2）先规定格点位置，再垂直平分两相邻结点的连线作格线，形成的网格即为水均衡区第6页，共49页，星期日，2025年，2月5日MODFLOW网格系统第7页，共49页，星期日，