同伦分析法推求greenrampt方程显格式精确解.pdfVIP

第４８卷 第１６期 人　民　长　江 Ｖｏｌ．４８，Ｎｏ．１６ ２０１７年 ８月　 　 Ｙａｎｇｔｚｅ　Ｒｉｖｅｒ　 　 Ａｕｇ．， ２０１７ 　　文章编号：１００１－４１７９（２０１７）１６－００８８－０５ 同伦分析法推求 Ｇｒｅｅｎ－Ａｍｐｔ方程显格式精确解 １ ２ 唐岳 灏 ，姜 清 辉 （１．长江工程职业技术学院，湖北 武汉４３０２１２；　２．武汉大学 土木建筑工程学院，湖北 武汉４３００７２） 摘要：针对 Ｇｒｅｅｎ－Ａｍｐｔ（ＧＡ）积水入渗模型只有隐格式解的问题，利用同伦分析法（ＨＡＭ）给出ＧＡ方程的显 式解析解。首先构造一组基底为幂函数的显示解表达，通过选择合适的线性和非线性算子，定义零阶形变方 程。为保证 ＨＡＭ的系数遍历原则，给出了零阶形变方程中辅助函数的选择形式，还探讨了辅助参数对不同 级数解收敛性的影响。通过与其他近似解的比较，验证了所提方法具有理论严谨完备、计算稳定性好、结果精 度高的特点，能有效模拟降雨条件下积水入渗的过程。 关　键　词：Ｇｒｅｅｎ－Ａｍｐｔ模型；同伦分析；累积入渗量；显格式；积水入渗 中图法分类号：ＴＶ１３５　　　文献标志码：Ａ ＤＯＩ：１０．１６２３２／ｊ．ｃｎｋｉ．１００１－４１７９．２０１７．１６．０１９ 显式积分算法得到了一个非迭代型的计算格式，但在 １　研究背景 ［８］ 本质上仍属于数值算法 。Ａｌｉ针对湿润锋深度与积 Ｇｒｅｅｎ－Ａｍｐｔ（ＧＡ）入渗方程是由 Ｇｒｅｅｎ和 Ａｍｐｔ 水深度和湿润锋吸力之和比值的不同范围，用二次函 在 １９１１年提出的基于毛管理论的积水入渗模型，与描 ［９］ 数逼近形式对 ＧＡ模型进行化简 ，但由于湿锋深度 述非饱和土壤水分运动复杂形式的 Ｒｉｃｈａｒｄｓ方程相 本身为未知量，因此需要试算才能确定比值的具体区 比，该方程形式简单，具有明确的物理意义。ＧＡ方程 间。Ｎａｓｓｅｒｉ采用改进 Ａｄｏｍｉａｎ分解和变分迭代法，分 所需土壤的物理参数也比较容易获得，这使得其在非 别推求出两种不同形式的无穷级数的近似解，但关于 饱和土的入渗过程研究中有着广泛的应用。 收敛性则没有讨论［１０］。Ｂａｒｒｙ在忽略积水深度的前提 ＧＡ入渗方程在实际应用中存在的主要问题是模 下，利用 ＬａｍｂｅｒｔＷ 函数给出一个精确解，但 Ｌａｍｂｅｒｔ 型所给出的累积入渗量 Ｆ、入渗速率ｆ与时间ｔ为隐函 Ｗ函数本身为隐函数，因此还必须借助近似解才能显 数关系，在计算时需要借助数值方法中的迭代计算， ［１１］ 式计算 。 应用不是很方便。常见的数值求解方法有 Ｒｕｎｇｅ－ 若对上述近似解按照数学形式进行分类，可以分 ［１］ ［２］ Ｋｕｔｔａ 和 Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ等 。为了避免迭代计 为无理式解（例如 Ｌｉ，Ｓｔｏｎｅ，Ａｌｍｅｄｅｉｊ，Ａｌｉ等）和级数解 算，学者们对 Ｇｒｅｅｎ－Ａｍｐｔ方程解的显式计算方法做 （例如 Ｓｅｒｒａｎｏ，Ｎａｓｓｅｒｉ等）两大类。这些方法都存在一 了很多研究。例如 Ｌｉ得到一个具有二次逼近的近似 些不足。对于无理式解，其受推求方法的影响很大，一 ［３］