《一维畦灌施肥地表水流与溶质运移耦合模型——Ⅱ.模型验证 One-dimensional model for coupling surface water flow with solute transport for border fertigationII.Model verification》.pdfVIP

第22卷第2期 水科学进展 V01．22，No．2 2011年3月 ADVANCESINWATERSCIENCE Mar．．2011 10313．2236．016 DOI：CNKI：32-1309／P．201 一维畦灌施肥地表水流与溶质运移耦合模型 ——Ⅱ．模型验证 章少辉1”，许 迪1’2，李益农1’2，白美健L2 (1．中国水利水电科学研究院水利研究所．北京100048；2．国家节水灌溉北京工程技术研究中心。北京100048) 摘要：基于典型畦灌施肥试验观测结果及其模拟结果．对比分析利用混合数值解法和Roe有限体积法分别求解一 维畦灌施肥地表水流与溶质运移过程控制方程在数值稳定性与收敛性、计算精度与效率上的差异，验证混合数值 解法的计算性能与模拟效果。结果表明．混合数值解法比Roe有限体积法表现出更佳的数值稳定性和收敛性，产生 的水平衡误差和平均相对误差较低，同样度量环境下的计算效率提高2．5倍以上。基于混合数值解法的一维畦灌施 肥地表水流与溶质运移耦合模拟方法可明显增强数值计算的稳定性和收敛性，有效提高计算精度和效率，为开展 畦灌施肥系统设计与评价提供了有效的数值模拟工具。 关键词：畦灌施肥；混合数值解法；稳定性；收敛性；计算精度与效率 中图分类号：0351．2；$275．3文献标志码：A 文章编号：1001J02-0196-07 利用Roe有限体积法同步耦合求解一维全水动力学畦灌模型和守恒型对流-扩散方程，可使数值模拟结 果在局部与整体时空上保持较好的质量与动量守恒性．体现出流速与溶质对流间的数值互动耦合效应…， 但受时间与空间高精度格式构造复杂等因素制约，模拟易导致过高的计算耗时，难以在数值稳定性、收敛 性、计算精度及效率方面达到协调统一【2l。为此，作者基于隐．显混合时间格式，利用有限差分法、AUSM 与中心格式有限体积法和有限单元法分别对一维畦灌施肥地表水流与溶质运移耦合模型中的各矢量项进行空 间离散，对形成的控制方程代数方程组进行数值求解，构建起基于混合数值解法的一维畦灌施肥地表水流与 溶质运移耦合模型¨】。本文基于典型畦灌施肥试验数据，对比分析混合数值解法和Roe有限体积法之间在 数值稳定性与收敛性、计算精度与效率上的差异，验证基于混合数值解法的一维畦灌施肥地表水流与溶质运 移耦合模型的模拟效果。 1 畦灌施肥试验和数值解法 1．1 典型畦灌施肥试验 1．1．1 试验处理 在国家节水灌溉北京工程技术研究中心大兴试验基地内，于2009年4月5日开展冬小麦返青期典型畦 mm。试验区内表土 灌施肥试验。该基地位于北京市南郊，属半干旱温带大陆性季风气候，年均降水量540 为砂质壤土，平均干容重I．54S／era3，灌溉水源为承压地下水。 收稿日期：2010-04-22：圈络出版时问：20ll-03·13 网络出版地址：http：／／www．cnki．net／kcm矗／detail／32．1309．P2236．016．html 国水利水电科学研究院科研专项项目(节集1038) 作者篱介：章少辉(1977一)。男．汉，河北石家庄人．工程师．博士研究生．主要从事田问节水灌溉技术研究o E．mail：ahangsh@iwhr．∞m 万方数据 第2期 章少辉，等：一维畦灌施肥地表水流与溶质运移耦合模型——Ⅱ．模型验证 197 典型畦灌施肥试验分为两个处理，入畦单宽流量分别为3L／(s·m)和6L／(s·m)，灌水时间分别为60 min和40min，畦田规格均为长102m、宽2 em m，平均地面坡度4／10000。畦面微地形标准偏差分别为2．21 和3．13em。 试验前将12kg硫酸铵化肥放入储液罐内，配水使其充分