spac系统计算机模拟.pptVIP

肖娟 zhangxd626@163.com  太原理工大学 水利科学与工程学院 目 录 1 引言 考察 Richardson 格式的稳定性。 用 表示计算 所产生的误差，如果右端 无误差存在， 则 满足： 取 （14） 假设k - 1层之前无误差存在。即 ，而在第k层产生了 误差。 ，这一层其它点也无误差，而且在计算过程 中不再产生新的误差，利用（14）式算出误差? 的传播如下表： r = ? 时 Richardson 格式的误差传播 j j0 – 4 j0 – 3 j0 – 2 j0 –1 j0 j0 +1 j0 +2 j0 +3 j0 +4 k ?? ? -2? ? ? -4? 7? 4? ? ? -6? 17? -24? 17? -6? ? ? -8? 31? -68? 89? -68? 31? -8? ? -10? 49? -144? 277? -388? 277? -144? 49? -10? 71? -260? 641? -109? 1311? -109? 641? -260? 71? r≤ 1/2 时古典显格式的误差传播 j j0 – 4 j0 – 3 j0 – 2 j0 –1 j0 j0 +1 j0 +2 j0 +3 j0 +4 k ?? 0.5? 0 0.5? 0.25? 0 0.5? 0 0.25 ? 0.125? 0 0. 375? 0 0.375 ? 0 0.125 ? 0.0625? 0 0.25? 0 0.375? 0 0.25? 0 0.0625? 如果选用 r = ? 时的古典显格式，误差方程为： 稳定性概念： 差分格式关于初值稳定的实际含义是：如果其解 在某一层存在误差，则由它引起的以后各层上的误差不 超过原始误差的M倍（M为与? 无关的常数）。 因此，在稳定的条件下，只要初始误差足够小， 以后各层的误差也能足够小。 以上构造的几种差分格式中， 古典格式：r ≤1/2时稳定 古典陷格式：绝对稳定 Richardson格式：绝对不稳定 六点对称格式：绝对稳定。 （1）水流为稳态和无其反应情况下的溶质运移方程 考虑六点对称格式，在点 列方程， 令 ， ，有： 舍去o(h2)，o(? 2)得到六点对称格式： 令 上式变为 令 上式可以写成： j = 1, 2,…, N –1 对时间变量用向后差分，对空间变量用中心差分， 可得到隐格式： 令 整理得： （2）．非稳态方程 非稳态一维垂直流情况下，土壤溶质的基本方程为： 式中容积含水量? 的求法如下： 式中，h为负压水头；t为时间；z为到原点的距离（cm）， 向下为正；C(h)为容水度， ，K(h)为土壤导水率， 可由一些常用的经验公式算出。 先解非饱和垂直水流方程 求出h后，用水分特征曲线换算成相应的? 值。v =q/? ，q为通量。 式中： 将以上方程整理后可写成： 2、有限元法 设有微分方程 定义在由边界? 围成的区域以上，L为微分算子。 设{?j}（j = 1, 2, …，n, …），是一完备的函数系。伽辽金方法是求形如 的近似解，其中aj (j = 1, 2, …, n )为待定常数。un称为试探函数， ?j称为形状函数（或基函数，插值函数，为{?j}(j = 1, 2, …, n,…) 中前n个线性无关的函数）。 若u是方程（*）的精确解，则必有 （1） (2) * SPAC系统计算机模拟 2 土壤水流方程 3 土壤溶质运移方程 4 微分方程数值解 5 HYDRUS 模型应用 土壤溶质运移所研究的是溶于土壤水中的溶质在 土壤中运移的过程、规律和机理。 土壤溶液中溶质不但在组成和形态处于一个永恒的 变化中，而且它还无时不刻在运动中，土壤水（溶液） 受到自然和人为环境因子的作用，如降雨、灌溉、气温、 风速和植物根系吸收等影响而在土壤中不断地进行