A Simple Line-Element Model for Three-Dimensional Analysis of Steady Free Surface Flow through Porous Media

Abstract

Considering the fact that only pores can transport water, pores in the homogeneous control volume are conceptualized as a three-dimensional orthogonal network of line elements, which is in contrast to the continuum hypothesis in traditional numerical approaches. The related flow velocity, hydraulic conductivity and continuity equation equivalent to the continuum model are formulated based on the principle of flow balance. Subsequently, the unified form for flow velocity and continuity equation is established based on the local coordinate system, and a finite line-element method is developed, in which three-dimensional steady free surface flow is reduced to one-dimensional form, and the numerical difficulty is greatly decreased. The proposed line-element model is validated by the good agreements of free surface locations with other methods through steady flow in a rectangular dam and a right trapezoidal dam, respectively. It is found that the proposed line-element model is not heavily dependent on the mesh size and penalty parameter. Steady free surface flow on the left bank abutment slope of the Kajiwa Dam in Southwestern China is further evaluated, and a parabolic variational inequality algorithm based on the continuum model is also employed for comparison. The consistent results indicate that the proposed line-element model can capture the steady free surface flow behavior as well as the continuum-based method. Moreover, the proposed line-element model can rapidly achieve accurate solutions whether for simple examples or for complicated engineering applications.