Nonlinear numerical simulation of rainwater infiltration through road embankments by FEM

Abstract This paper presents the results of nonlinear numerical simulations carried out to analyze rainwater infiltration through a road embankment structure using the finite element method (FEM) technique. The analysis of this paper includes a simplification of a cross-sectional road embankment located in Santander (Northern Spain) and considers different realistic hydraulic material parameters, according to the relevant literature, in order to study their effect on water movement throughout the section over time. The model solved the nonlinear Richards’ equation using the constitutive relationships of van Genuchten–Mualem. Finally, the analysis of the numerical results reveals four main conclusions: (1) models based on FEM are suitable for analyzing rainwater infiltration in road structures, (2) the hard-shoulder area of the highway represents an important infiltration point for rainwater in road structures, (3) the embankment material can greatly affect the soil’s water content in this part of the road, and (4) the analysis of the soil water retention curve of the embankment material provides relevant information about its hydraulic behavior under unsaturated conditions.

[1]  Defne Apul,et al.  A REVIEW OF WATER MOVEMENT IN THE HIGHWAY ENVIRONMENT: IMPLICATIONS FOR RECYCLED MATERIALS USE , 2002 .

[2]  D. Fredlund,et al.  Soil Mechanics for Unsaturated Soils , 1993 .

[3]  S. P. Neuman,et al.  Finite element simulation of flow in saturated - unsaturated soils considering water uptake by plants , 1974 .

[4]  J. Z. Zhu,et al.  The finite element method , 1977 .

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

[6]  Richard H. Gallagher Viscous flow and hydrodynamics , 1975 .

[7]  S. P. Neuman,et al.  Field test of a modified numerical model for water uptake by root systems , 1974 .

[8]  Kok-Kwang Phoon,et al.  NUMERICAL STUDY OF FINITE ELEMENT METHOD BASED SOLUTIONS FOR PROPAGATION OF WETTING FRONTS IN UNSATURATED SOIL , 2004 .

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

[10]  H. Rügner,et al.  Model-based prediction of long-term leaching of contaminants from secondary materials in road constructions and noise protection dams. , 2009, Waste management.

[11]  L. A. Richards Capillary conduction of liquids through porous mediums , 1931 .

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

[13]  Lars-Christer Lundin,et al.  Modeling Water Flow Patterns in Flexible Pavements , 2005 .

[14]  R. Carsel,et al.  Developing joint probability distributions of soil water retention characteristics , 1988 .

[15]  Defne S Apul,et al.  Modeling hydrology and reactive transport in roads: the effect of cracks, the edge, and contaminant properties. , 2007, Waste management.

[16]  Michal Kuraz,et al.  Solving the Nonstationary Richards Equation With Adaptive hp-FEM , 2011 .

[17]  Andrew Dawson,et al.  Water in Road Structures , 2009 .

[18]  W. G. Gray,et al.  Finite Element Simulation in Surface and Subsurface Hydrology , 1977 .

[19]  S. P. Neuman,et al.  SATURATED-UNSATURATED SEEPAGE BY FINITE ELEMENTS , 1973 .

[20]  Harianto Rahardjo,et al.  Response parameters for characterization of infiltration , 2010 .