应用间断有限元方法求解一维对流占优土壤水流问题

摘要

在某些特定条件下,Richards方程的解在时空上呈现陡峭的锋面。为能有效地模拟具有对流占优特性的非饱和多孔介质中的水流问题,推广一种内部惩罚间断有限元(Interior penalty discontinuous Galerkin,IPDG)方法应用于一维非饱和土壤水入渗问题的模拟。针对具有van Genuchten-Mualem模型和Dirichlet入渗边界条件的Richards方程,分别采用间断有限元法和标准有限元方法求解。借助于相对L2模和相对最大模误差进行讨论。几种不同质地的均质土壤水入渗的数值算例结果表明:相比标准有限元方法,间断有限元方法在选取的4种不同网格剖分单元结点上能够有效地模拟非饱和对流占优土壤水流问题,并且能够获得准确的全局质量守恒。