一维畦灌施肥地表水流与溶质运移耦合模型-水科学进展.PDFVIP

第卷第期 水科学进展 ， ２２ ２ Ｖｏｌ２２ Ｎｏ２ 　 　 ２０１１年３ 月 ＡＤＶＡＮＣＥＳ ＩＮ ＷＡＴＥＲ ＳＣＩＥＮＣＥ 　 Ｍａｒ ，２０１１　 ： ： ＤＯＩ ＣＮＫＩ ３２１３０９／ Ｐ２２３６０１６ 一维畦灌施肥地表水流与溶质运移耦合模型 ——— 模型验证 Ⅱ １，２ １，２ １，２ １，２ 章少辉，许　 迪，李益农，白美健 （１ 中国水利水电科学研究院水利研究所，北京　 １０００４８；２ 国家节水灌溉北京工程技术研究中心，北京　 １０００４８） 摘要：基于典型畦灌施肥试验观测结果及其模拟结果，对比分析利用混合数值解法和Ｒｏｅ 有限体积法分别求解一 维畦灌施肥地表水流与溶质运移过程控制方程在数值稳定性与收敛性、计算精度与效率上的差异，验证混合数值 解法的计算性能与模拟效果。结果表明，混合数值解法比Ｒｏｅ有限体积法表现出更佳的数值稳定性和收敛性，产生 的水平衡误差和平均相对误差较低，同样度量环境下的计算效率提高２５倍以上。基于混合数值解法的一维畦灌施 肥地表水流与溶质运移耦合模拟方法可明显增强数值计算的稳定性和收敛性，有效提高计算精度和效率，为开展 畦灌施肥系统设计与评价提供了有效的数值模拟工具。 关键词：畦灌施肥；混合数值解法；稳定性；收敛性；计算精度与效率 中图分类号： ； 文献标志码： 文章编号： （ ） Ｏ３５１２ Ｓ２７５３　 　 　 Ａ　 　 　 １００１６７９１ ２０１１ ０２０１９６０７ 利用Ｒｏｅ有限体积法同步耦合求解一维全水动力学畦灌模型和守恒型对流扩散方程，可使数值模拟结 ［］１ 果在局部与整体时空上保持较好的质量与动量守恒性，体现出流速与溶质对流间的数值互动耦合效应， 但受时间与空间高精度格式构造复杂等因素制约，模拟易导致过高的计算耗时，难以在数值稳定性、收敛 ［２］  ＡＵＳＭ 性、计算精度及效率方面达到协调统一。为此，作者基于隐显混合时间格式，利用有限差分法、 与中心格式有限体积法和有限单元法分别对一维畦灌施肥地表水流与溶质运移耦合模型中的各矢量项进行空 间离散，对形成的控制方程代数方程组进行数值求解，构建起基于混合数值解法的一维畦灌施肥地表水流与 ［３］ 溶质运移耦合模型。本文基于典型畦灌施肥试验数据，对比分析混合数值解法和 有限体积法之间在 Ｒｏｅ 数值稳定性与收敛性、计算精度与效率上的差异，验证基于混合数值解法的一维畦灌施肥地表水流与溶质运 移耦合模型的模拟效果。 １　 畦灌施肥试验和数值解法 １ １　 典型畦灌施肥试验 １ １ １　 试验处理 在国家节水灌溉北京工程技术研究中心大兴试验基地内，于 年月日开展冬小麦返青期典型畦 ２００９ ４ ５ 灌施肥试验。该基地位于北京市南郊，属半干旱温带大陆性季风气候，年均降水量５４０ ｍｍ。试验区内表土 为砂质壤土，平均干容重１５４ ｇ／ ｃｍ ，灌溉水源为承压地下水。３ 收稿日期：２０１００４２２；网络出版时间：２０１１０３１３ 网络出版地址： ：