饱和-非饱和土中渗流Richards方程有限元算法

发展了一种求解理查森方程的一般有限元算法.该方法采用积分法处理孔隙水压力对时间的导数项,采用集中质量技术处理有限元方程中质量矩阵来保证数值稳定.所采用的质量守恒迭代方法不须改变迭代方式,采用一般的Picard迭代方法.该方法能求解入渗、地下水位瞬变和排水等范围广泛的饱和-非饱和渗流问题.对3个已公开发表具有详细试验数据的算例的模拟表明,该方法对入渗锋、稳定渗流地下水位和非稳定渗流溢出面都模拟很好.Pieard迭代方法效率很高,且无数值振荡发生.