第7章 土壤水分运动的数值解.ppt

7.4 土壤水分运移方程的线性化方法 以基质势水头h为因变量的一维垂直向土壤水分运动基本方程为： 7.4 土壤水分运移方程的线性化方法 7.4 土壤水分运移方程的线性化方法 7.4 土壤水分运移方程的线性化方法 （2）线性外推法（已知 时刻，外推 时刻） 假定含水率在相邻时段内是线性变化的，根据前一个时段末的含水率，用直线外推近似求出计算时段末的含水率： 7.5 线性化方法与土壤水分参数的取值 7.5.2.1直接法 先求出两结点的参数，然后 (1)取两结点参数的算术平均值： 7.5 线性化方法与土壤水分参数的取值 7.5.2 2间接法 先求出两结点处含水率的算术、调和、几何平均值，然后再求出相应的参数值。 7.5 线性化方法与土壤水分参数的取值 7.5.2 3其他方法 （1）求出的参数和两结点处的参数和，三点参数取平均，故可称为“三点式”： （a）若取绝对误差，判断条件为： or 7.5.1线性化方法 7.5.1.3迭代法 （b）若取相对误差，判断条件为： or 7.5.1线性化方法 7.5.1.3迭代法 ?1与?2 的取值可根据实际问题的精度要求而定。 此法计算稍复杂，但精度较高，数值计算时经常采用。为避免迭代次数太多，一般可设置一个判断语句：当迭代次数m≥5时，将时间步长  减半，再重新计算。 7.5.1线性化方法 7.5.1.3迭代法 7.5.2参数取值 (2)取两结点参数的调和平均值： (3) 取两结点参数的几何平均值： (1)算术平均值： (2)调和平均值： (3)几何平均值： （2） （3）由邻近结点的参数代替，如 7.6 边界条件的处理及追赶法求解三对角方程 前已叙及，采用显式差分格式可直接求解，采用隐式或中心差分格式时总可以得到如下形式的三对角方程组（以? 方程为例）： 共n - 2个方程，再加上上、下边界条件，即可组成n ? n 阶的n元一次方程组，线性化后求解该方程组可获得原定解问题的解。 7.6.1边界条件的处理 上、下边界不论为第一、第二或第三类边界，总可以化为： (2) (1) 7.6 边界条件的处理及追赶法求解三对角方程 7.6.1边界条件的处理 验 证： ? 当上、下边界为第一类边界时 ? ? ? í ì = = = ? = = = ? T + + ) ( , 0 , 1 ) 2 .( Eqn ) ( , 0 , 1 ) 1 .( Eqn 1 2 n 1 1 1 1 1 k n n k t f H E F t f H G F 7.6 边界条件的处理及追赶法求解三对角方程 7.6.1边界条件的处理 ? 当上、下边界为第二、三类边界时 ? ? ? ? ? í ì = + - - - = + - - - T + - - + - + - + + + + ) ( ) , ( 1 2 2 / 1 1 1 1 1 2 1 1 1 1 1 2 3 1 2 1 1 1 2 2 3 k n n n k n k n n k k k k t q K z z D t q K z z D q q q q q * 第7章 土壤水分运动的数值解 7.1 概述 7.2 差分的概念及分类 7.3 显式差分、隐式差分和中心差分 7.4 土壤水分运移方程的线性化方法 7.5 线性化方法与土壤水分参数的取值 7.6 边界条件的处理及追赶法求解三对角方程 7.7 垂直一维非饱和土壤水流计算流程图 第7章 土壤水分运动的数值解 主要参考书目 Huyakorn, P. S. and G. F. Pinder. Computational methods in subsurface flow. Academic Press, New York, 1983. 薛禹群，谢春红，水文地质学的数值法，煤炭工业出版社，1980. 土壤水分运动方程的求解方法： ? 解析解（Analytical solution） ? 数值解（Numerical solution） (1) 有限差分法（Finite Difference method） (2) 有限单元法（Finite Element method） 7.1 概述 7.2 差分的概念及分类 差分的由来（Taylor展开） if ? = ? (x) exists, then (1) (2) ? 第7章 土壤水分运动的数值解 when ?x ? 0, then ① 向前