2.4.2 数值解
UNSAT解决了描述地表和含水层之间一系列竖向土柱的Rich-ards方程。对于每一个土柱需定义下列条件。
· 初始条件: 模拟开始时,沿整个土柱的毛细管压力水头,h(z,t=0);
· 边界条件: 土柱顶部和底部的毛细管压头水头h (或与h相关的条件);
· 土壤参数: k(h) 和C(h);
· 计算网格空间步长: 空间步长△z和时间步长△t。
图2.19中示意性地给出为解决土柱Richards方程所需的程序和数据,使用Godunov型有限体积格式求解方程。为达到这一目的,式(2.4.4) 可以重写为
图2.19 求解Richards方程所需的计算网格和数据
来源: 作者。
流动区域在z(j-1/2) 和z(j+1/2) 之间被分成了一系列的有限体积j。在每一个计算网格内任一时间步 ...... (共786字) [阅读本文]>>