欧拉法融合拉格朗日法高效模拟灌溉二维地表水运动规律.pdfVIP

第 35 卷 第 12 期 农 业 工 程 学 报 Vol.35 No.12 150 2019 年 6 月 Transactions of the Chinese Society of Agricultural Engineering Jun. 2019 欧拉法融合拉格朗日法高效模拟灌溉二维地表水运动规律 ※ 史 源，章少辉 ，白美健，李益农 （中国水利水电科学研究院流域水循环模拟与调控国家重点实验室，北京 100038 ） 摘 要：大规模现代化农业灌溉管理下，为实现快速高效地获知灌溉水运动及其分布的目的，该文基于二维浅水方程组 的欧拉-拉格朗日混合型表达形式，提出了一种高效简洁的欧拉-拉格朗日混合解法。该解法的基本物理变量被严格地定 义在欧拉型非结构化三角形有限体积单元格上，且变量在单元格之间呈现出阶梯分布状态，以精准地捕捉各类地表浅水 波动并有效地保持质量守恒性；由于控制方程中不存在水运动的对流梯度项（或位置加速度项），仅通过拉格朗日迹线追 踪的形式获得未知与已知时间步之间的变量关系，故与广泛应用的欧拉解法相比，离散格式表达式极为简洁易用；在地 表水运动的干湿边界处，地表水位梯度项被做了修正，以严格地保证各物理量之间的数值平衡，进而能高精度的模拟整 个畦田内的地表水流推进/消退全过程。为验证模型的模拟性能，选取一种高效的欧拉解法（非迭代型全隐式标量耗散有 限体积法）求解二维浅水方程组做为对比模型，基于 3 个典型畦灌试验的实测数据，从模拟精度、质量守恒性和计算效 率 3 个方面，对比分析了 2 种数值解法的性能。结果表明，2 种解法在模拟精度方面相差无几，且欧拉-拉格朗日混合解 法比欧拉解法具有更好的质量守恒性；在计算效率方面，欧拉-拉格朗日混合解法比欧拉解法的效率提高了约 5.3 倍。故 该文提出的二维浅水方程组的欧拉-拉格朗日混合解法，更适用于二维灌溉地表水运动的模拟分析。 关键词：灌溉；非线性方程；数值分析；欧拉解法；拉格朗日解法；二维；浅水方程组 doi ：10.11975/j.issn.1002-6819.2019.12.018 中图分类号：S11+4; S275.3 文献标志码：A 文章编号：1002-6819(2019)-12-0150-07 史 源，章少辉，白美健，李益农. 欧拉法融合拉格朗日法高效模拟灌溉二维地表水运动规律[J]. 农业工程学报，2019， 35(12)：150－156. doi ：10.11975/j.issn.1002-6819.2019.12.018 Shi Yuan, Zhang Shaohui, Bai Meijian, Li Yinong. Euler-Lagrange hybrid numerical simulation for two-dimensional surface water flow in irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(12): 150－156. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.12.018 的取值范围，但仍属于有条件稳定。为此，（迭代或非  0 引 言