Simulation procedure of unconfined seepage with an inner seepage face in a heterogeneous field

An inner seepage face phenomenon is given and a numerical simulation procedure has been developed. It may appear at the interface of two materials when an unconfined seepage flows from a porous media to a coarser porous media with a higher permeability. Inaccuracy and divergent problems may arise both in a saturated-only and in a variably saturated analysis while an inner seepage face is not simulated with a special procedure. The position of the seepage face is determined during the nonlinear iteration process and the flux of the inner seepage face nodes is transferred to the downstream side nodes. Validity and efficiency of the procedure are illustrated by the simulation of two dimensional steady state seepage examples of heterogeneous zoned dams which is usually used to validate algorithms. An analysis of a three-dimensional earth core rockfill dam is also presented here. The procedure can also be applied to general transient seepage problems.

[1]  J. Prévost,et al.  Flow through porous media: A procedure for locating the free surface , 1987 .

[2]  M. I. Herreros,et al.  Application of level-set approach to moving interfaces and free surface problems in flow through porous media , 2006 .

[3]  Colin B. Brown,et al.  Darcy Flow Solutions with a Free Surface , 1967 .

[4]  R. Borja,et al.  One the solution of elliptic free-boundary problems via Newton's method , 1991 .

[5]  Farhang Daneshmand,et al.  Unconfined seepage analysis in earth dams using smoothed fixed grid finite element method , 2012 .

[6]  Pierre Perrochet,et al.  On the primary variable switching technique for simulating unsaturated–saturated flows , 1999 .

[7]  M. Tamer Ayvaz,et al.  Modeling three-dimensional free-surface flows using multiple spreadsheets , 2007 .

[8]  Carmelo Gentile,et al.  A nonlinear programming analysis of unconfined steady‐state seepage , 1987 .

[9]  Jean-Pierre Bardet,et al.  A practical method for solving free-surface seepage problems , 2002 .

[10]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[11]  K. Bathe,et al.  Finite element free surface seepage analysis without mesh iteration , 1979 .

[12]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[13]  Mengxi Wu,et al.  A finite-element algorithm for modeling variably saturated flows , 2010 .

[14]  Y M Cheng,et al.  An efficient method for the free surface seepage flow problems , 1993 .

[15]  Chandrakant S. Desai.  Seepage Analysis of Earth Banks under Drawdown , 1972 .

[16]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[17]  S. P. Neuman,et al.  Finite Element Method of Analyzing Steady Seepage with a Free Surface , 1970 .

[18]  J. Oden,et al.  Theory of variational inequalities with applications to problems of flow through porous media , 1980 .

[19]  W. D. Liam Finn.  Finite-Element Analysis of Seepage Through Dams , 1967 .

[20]  R. L. Cooley.  Some new procedures for numerical solution of variably saturated flow problems , 1983 .

[21]  Chandrakant S. Desai,et al.  Finite element residual schemes for unconfined flow , 1976 .

[22]  R. Allan Freeze,et al.  Influence of the Unsaturated Flow Domain on Seepage Through Earth Dams , 1971 .