《一维畦灌施肥地表水流与溶质运移耦合模型——Ⅰ.模型建立 One-dimensional model for coupling surface water flow with solute transport for border fertigationⅠ.Model solution》.pdfVIP

第22卷第2期 水科学进展 V01．22。No．2 2011年3月 ADVANCESINWATERSCIENCE Mar．．201I 10313．2236．015 DOl：CNKI：32·1309／P．201 一维畦灌施肥地表水流与溶质运移耦合模型 ．——I．模型建立 章少辉1’2，许 迪k2，李益农1”，白美健1·2 (1．中国水利水电科学研究院水利研究所．北京100048；2．国家节水灌溉北京工程技术研究中心．北京100048) 摘要：提高畦灌施肥地表水流与溶质运移数值模型的稳定性、收敛性及计算精度．有利于改善地面畦灌施肥系统的 设计与评价工作。在利用隐-显混合时间格式对一维畦灌地表水流与溶质运移耦合模型中包含的各矢量项进行时间 离散基础上，借助有限差分法、有限体积法和有限单元法分别对由隐时间格式生成的物理矢量线性近似式空间导 数、物理矢量空问导数、溶质扩散矢量和地形矢量项进行空间离散。对最终形成的控制方程代数方程组进行数值求 解，构建起基于混合数值解法的一维畦灌施肥地表水流与溶质运移耦合模型。 关键词：畦灌施肥；混合数值解法；耦合；水流运动；溶质运移；时空离散 中圈分类号：0351．2；S275．3文献标志码：A 文章编号：1001-679112011)02m189．07 畦灌施肥是将预先溶解的化肥随灌溉水流一并施用的地面施肥方式，与化肥撒施后再进行灌溉的传统方 法相比，具有高效利用化肥、减少肥分流失、降低劳力投入等突出特点¨剖。开展畦灌施肥地表水流与溶质 运移过程模拟研究，可为比较不同施肥方式下的畦灌系统性能提供对比依据，对于推广先进的施肥方式，减 轻由农业施肥不当带来的面源污染等问题具有重要的理论与现实意义¨J。 在数值模拟畦灌施肥过程中，通常是在利用全水动力学、零惯量等一维畦灌模型模拟地表水流运动过程 基础上，采用守恒型对流一扩散方程模拟地表溶质运移过程¨引。近年来，伴随着计算机技术与数值解法的快 速发展，国外一些学者构建起一维畦灌施肥地表水流与溶质运移耦合模型，并利用Roe有限体积法同步求解 一维全水动力学畦灌模型和守恒型对流一扩散方程”…。与以往顺次非耦合求解地表水流和溶质运移控制方程 的传统数值解法相比，同步耦合解法能使数值模拟结果在局部与整体时空上保持较好的质量与动量守恒性。 并体现出流速与溶质对流间的数值互动耦合效应冲1。然而，受时间与空间高精度格式构造复杂等因素制约， Roe有限体积法在模拟非定常态的畦灌施肥地表水流与溶质运移过程中。易导致过高的计算耗时，模拟过程 难以在数值稳定性、收敛性、计算精度及效率方面达到协调统一。为此，有必要针对一维畦灌施肥地表水流 与溶质运移耦合模型中包含的各矢量项特征，采用不同的数值解法分别对其进行解算，以期达到提高计算精 度、增强数值计算稳定性与收敛性、缩短模拟计算耗时的目的o¨，为畦灌施肥系统的设计、管理与性能评 价提供更有效的数值分析工具。 本文基于隐-显混合时间格式，利用有限差分法、有限体积法和有限单元法分别对一维畦灌施肥地表水 流与溶质运移耦合模型中的各矢量项进行空间离散。对形成的控制方程代数方程组进行数值求解，构建起基 于混合数值解法的一维畦灌施肥地表水流与溶质运移耦合模型。 l旬3-13 收稿日期：2010拼．22；网络出版时问：201 网络出版地址：http：／／www．cnki．net／kcms／detail／32．1309．P2236．015．htmi 国水利水电科学研究院科研专项项目(节集1038) 作者简介：章少辉(1977一)．男．河北石家庄人，工程师．博士研究生。主要从事田间节水灌溉技术研究。 E-mail：zhangsh@iwhr．corn 万方数据 190 水科学进展 第22卷 1 一维畦灌施肥地表水流与溶质运移耦合模型 1．1耦合控制方程 采用基于Saint．Venant方程的全水动力学模型描述畦灌施肥地表水