机构地区: 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室
出 处: 《气候与环境研究》 1998年第1期74-82,共9页
摘 要: 利用Galerkin有限元集中质量法,讨论了非饱和土壤水流动的数值计算模型,用以模拟均质土壤,地下水埋藏很深且在不同的初始和边界条件下的水分运动。通过将第一类边界条件齐次化,而对于第二类非线性边界条件采用变分的办法将其直接转化为已知边界通量的计算,从而容易且较好地处理了带通量约束边界条件。对于入渗或蒸发问题,用通常的离散方法容易产生振荡非物理解。我们采用有限元集中质量法,将质量矩阵的非对角元集中到主对角元,避免了这种非物理解的产生。利用该方案统一地模拟如下情况:稳定入渗和蒸发、蒸发一段时间后达到风干率,以及波动入渗蒸发时土壤水势的变化状况。模拟表明,在已知有关土壤参数的情况下,该模型能较好地模拟非饱和流水势及含水量分布。 A numerical model for an unsaturated flow problem by using finite element method is established in order to simulate liquid moisture flow in an unsaturated zone with homogeneous soil and deep subsurface water, and with different initial and boundary conditions. For infiltration or evaporation problems, a traditional method usually yields oscillatory non-physics profiles. However, we obtain nonoscillatory solutions and evade non-physics solutions for these problems by using the mass-lumped finite element method. Moreover, the kind of boundary conditions is handled very well.