非饱和带有限分析数值模拟的误差分析.pdf

第２７卷 第 １期 水 科 学 进 展 Ｖｏｌ．２７，Ｎｏ．１ ２０１６年１月　 　 ＡＤＶＡＮＣＥＳ ＩＮ ＷＡＴＥＲ ＳＣＩＥＮＣＥ Ｊａｎ．，２０１６ ＤＯＩ： １０１４０４２／ ｊｃｎｋｉ３２１３０９２０１６０１００８ 非饱和带有限分析数值模拟的误差分析 张在勇， 王文科， 陈　 立， 王周锋， 段　 磊， 安可栋 （长安大学旱区地下水文与生态效应教育部重点实验室， 陕西 西安　 ７１００５４ ） 摘要： 由于描述包气带水分运移的Ｒｉｃｈａｒｄ方程具有高度非线性， 常用数值方法求解。 有限分析法可以较好地保持 原有问题的物理特性， 近年来被广泛应用于环境及农业工程领域。 基于不同的原理， 目前有两种有限分析法模拟 非饱和带水分运移问题， 即基于局部线性化有限分析法和基于Ｋｉｒｃｈｈｏｆｆ 变换有限分析法。 为探讨两者在求解非饱 和带水分运移问题时数值表现的差异性， 采用已有解析解对两种方法进行初步验证； 使用两个数值实验比较两种 方法的计算精度； 应用实测数据验证Ｋｉｒｃｈｈｏｆｆ变换有限分析法。 相比于局部线性化有限分析法， Ｋｉｒｃｈｈｏｆｆ变换有限 分析法能够更好地控制误差及获得更高精度的数值解。 该研究对于完善非饱和带水分运移有限分析数值模拟的理 论与方法有重要意义。 关键词： 非饱和流； Ｒｉｃｈａｒｄ方程； 有限分析法； Ｋｉｒｃｈｈｏｆｆ变换 中图分类号： Ｐ１２２　 　 　 文献标志码： Ａ　 　 　 文章编号： １００１⁃６７９１（２０１６）０１⁃００７０⁃１１ 包气带是连接地表、 地下水重要的纽带， 研究包气带水分运移对生态环境的保护、 地下水补给量与排泄 ［１⁃２］ 量的正确估算等具有非常重要的意义 。 描述非饱和带水分运移的Ｒｉｃｈａｒｄ方程具有高度的非线性， 因此 ［３］ ［４⁃６］ 难以获取其解析解 。 常利用数值方法求解， 较为常用的方法为有限差分法和有限元法 。 因为Ｒｉｃｈａｒｄ 方程同时具有双曲线型和抛物线型方程的双重性质， 所以合适的时间、 空间步长选取对于数值解的稳定性和 ［７⁃９］ 精确性极其重要 。 当选取不适当的时间、 空间步长时， 有限差分法及有限元法等数值方法在求解Ｒｉｃｈａｒｄ 方程时容易造成数值弥散、 数值振荡及质量不守恒等一系列问题。 ［１０］ ２０世纪８０年代初期， Ｃｈｅｎ等 提出了一种新的数值计算方法———有限分析法。 有限分析法的基本思 ［１０］ 想 是将控制方程的局部解析解组成整体的数值解。 由于有限分析数值解可以较好地保持原有问题的物理 特性， 格式通过自动调节有限分析系数来反映对流与扩散效应， 因此可以得到单调无振荡解， 数值稳定性 好。 自从该方法问世以来， 众多学者将其应用于探讨解决某些水文地质问题［１１⁃１６］。 为了求解Ｒｉｃｈａｒｄ方程时 获得精度高、 稳定性好的数值解， 一些学者将有限分析法推广于模拟包气带水分运移问题［１７⁃１９］。 由于有限 分析法在解析解的基础上建立， 为了获得局部解析解， 必须对Ｒｉｃｈａｒｄ方程进行特殊处理， 以便获得局部解 ［１７］ ［１８］ 析解。 目前有两种处理方法： ① 如任理 和Ｔｓａｉ等 对Ｒｉｃｈａｒｄ方程进行局部线性化， 在局部单元上推求 ［１９］ 局部解析解， 称为局部线性化有限分析法。 ② 如Ｚｈａｎｇ 等 使用 Ｋｉｒｃｈｈｏｆｆ 变换对 Ｒｉｃｈａｒｄ方程进行变换，