有限元集中质量法在非饱和土壤水流中的应用

利用galerkin有限元集中质量法，讨论了非饱和土壤水流动的数值计算模型，用以模拟均质土壤，地下水埋藏很深且在不同的初始和边界条件下的水分运动。通过将第一类边界条件齐次化，而对于第二类非线性边界条件采用变分的办法将其直接转化为已知边界通量的计算，从而容易且较好地处理了带通量约束边界条件。对于入渗或蒸发问题，用通常的离散方法容易产生振荡非物理解。我们采用有限元集中质量法，将质量矩阵的非对角元集中到主对角元，避免了这种非物理解的产生。利用该方案统一地模拟如下情况：稳定入渗和蒸发、蒸发一段时间后达到风干率，以及波动入渗蒸发时土壤水势的变化状况。模拟表明，在已知有关土壤参数的情况下，该模型能较好地模拟非饱和流水势及含水量分布。