适应复杂上表面边界条件的一维土壤水运动数值模拟.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 适应复杂上表面边界条件的一维土壤水运动数值模拟 陆垂裕 裴源生 中国水利水电科学研究院水资源研究所 北京 摘要 本文利用有限元法建立了一维土壤水运动数值模型 提出了适应降雨?灌溉 蒸发等复杂的上边界条件的处 理方法 可以用来统一计算土壤剖面降雨 灌溉入渗 地表积水 地表产流 蒸发 蒸腾以及当这些现象交替出现时 的水分运动过程 在数值计算过程中 采用了质量守恒和特殊的迭代求解格式以加快计算收敛速度和减少水量平 衡误差 模拟验证表明 本文模型提出的思路可行 计算方法合理并具有通用性 关键词 土壤水运动 数值模拟 有限元法 上边界条件 质量守恒 中图分类号 文献标识码 土壤水运动是多孔介质流体运动的一种重要形式 地表和地下的各种水文过程 如大气降水 地表 水的入渗 地表产流 潜水蒸发和植物的蒸腾等 都与土壤水的运动密切相关 对土壤水运动的数值模 拟研究在大气科学 土壤学 农业工程 环境工程和地下水动力学等方面均具有重要意义 本文通 过有限元法建立一维土壤水的数值模拟模型 针对土壤水运动数值模拟中对含水量的时间项比较敏感 的问题 利用质量集中法进行处理 并采用特殊的迭代格式以减少模拟过程中的水量平衡误差 一维土壤水数值模型 控制方程 以垂直向上为正 一维土壤水的运动方程 方程 可表达如下 式中 为土壤体积含水率 为土壤水负压 为根系吸水率或其它源汇项 为 饱和?非饱和水力传导度 本文模型采用 的修正模型 来描述三者之间的联系 控制方程的离散 用伽辽金有限元法对式 进行空间离散 则对于每一个计算节点 式中 为计算节点编号 为线性插值基函数 为一维模拟空间 将式 展开并利用分部积分法 有 式中 和 分别为模拟区域上 下边界的垂向位置 将模拟区域内的土壤水负压 用下式进行近似 收稿日期 基金项目 科技部西部开发重大项目 作者简介 陆垂裕 男 湖北黄石人 博士 高级工程师 主要从事地下水资源及水环境研究 式中 为待求系数 其解为各离散节点处的土壤水负压值 为区域内的计算节点总数 在式 中用 替换 并将积分化在各计算单元上进行可得 对上式的时间项用质量集中法 进行处理 即 对计算区域内的 个计算节点 应用上式 和式 进行离散处理可以得到矩阵方程组 式中 为节点边界流量项 对于非边界节点 由于插值基函数 在上下边界处的值 为 因此 时 时 对于时间的离散 为保持计算的稳定性 模型中采用隐式差分的格式 即 式中 代表当前的时间层 代表前一时间层 模型的求解通过 进行 由于方程系数 都是 的函数 因此以上方程组具有高度的非线 性 在每个 计算时段都需要通过迭代法进行求解 迭代过程对如何处理上述矩阵方程中的含水率 项十分敏感 模型中采用了一种称为 质量守恒 的方法 以减少计算过程中的水量平衡误差 这种 方法在迭代过程中将式 中的第一项分成两部分 即 式中 表示当前迭代 表示上一次迭代 式 中右端第二项对于当前迭代过程来说是已知的 将式 中右端第一项转化为用水头表示得 式中 为对角矩阵 其矩阵元素为 为计算节点处的容水度 当迭代结果满足精度要求时 式 中右端第一项基本上等于 这种特性可以有效减少求解过程 中的水量平衡误差 经过上述处理之后 暂不考虑边界条件的处理 迭代求解过程中矩阵方程可以表达如下 当前第 次迭代时 式 右端