ISRN International Scholarly Research Notices 2356-7872 Hindawi Publishing Corporation 10.1155/2014/189369 189369 Research Article An Analysis on Groundwater Recharge by Mathematical Model in Inclined Porous Media Pathak Shreekant P. Singh Twinkle Fonseca José L. C. Department of Applied Mathematics & Humanities S. V. National Institute of Technology Surat 395007 India svnit.ac.in 2014 2682014 2014 09 04 2014 15 05 2014 26 8 2014 2014 Copyright © 2014 Shreekant P. Pathak and Twinkle Singh. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The present paper discusses the analysis of solution of groundwater flow in inclined porous media. The problem related to groundwater flow in inclined aquifers is usually common in geotechnical and hydrogeology engineering activities. The governing partial differential equation of one-dimensional groundwater recharge problem has been formed by Dupuit's assumption. Three cases have been discussed with suitable boundary conditions and different slopes of impervious incline boundary. The numerical as well as graphical interpretation has been given and its coding is done in MATLAB.

1. Introduction

In nature, groundwater is a key element, in many geological and hydrochemical processes, which sustains spring discharge, river base flow, lakes, and wetlands. Groundwater, extracted from deep geological formations (called aquifers) through pumping wells, constitutes an important component of many water resource systems. Groundwater recharge problem has been discussed by many researchers with different viewpoints. Patel et al.  have obtained a series solution of moisture content in vertical groundwater flow through unsaturated heterogeneous porous media. Swaroop and Mehta  have obtained a solution to the problem of one-dimensional flow in unsaturated porous media taking finite element approach. Swartzendruber has used Philips  method to get graphical illustration of the mathematical solution for horizontal water function. Verma and Mishra  have obtained solution by similarity transformation of one-dimensional vertical ground water recharges through porous media. Mehta  has obtained an approximate solution by the method of singular perturbation technique. This paper is mathematically formulated by Dupuit’s assumption. Dupuit  based his assumptions on the observation that, in most groundwater flows, the slope of the phreatic surface is very small. He assumes that, for the small inclinations of free surface, the velocity is proportional to the slope of free surface but independent of the depth. The proportionality coefficient K (hydraulic conductivity) is a property of vascular plants, soil, or rock that describes the ease with which water can move through pore spaces or fractures. There are different values of K according to different types of porous media; that is, for homogeneous porous media it is constant and for heterogeneous porous media it is the function of distance x . For this mathematical model we considered K = 1 0 - 4  cm/sec. If one of the boundaries of aquifer is flow line or flow surface, the flow analysis is difficult especially for the analytical solution. The water table is flow line or flow surface when there is no recharge from unsaturated zone. Water table is neither flow line nor equipotential line in case of having input flow from vadose zone, and it crosses flow lines according to Todd and Mays . In the analytical method Dupuit’s assumption  has been applied. Hantush and Cruz  have used (1) for flow of water on inclined impermeable boundary: (1) q = - K ( y - m x ) d y d x , where m = tan θ .

Bear and Alexander  have presented the analytical solutions to the one-dimensional fluid flow and steady state on inclined boundary. They considered a constant value for flux Q 0 for the above condition. Their solution also handled the infiltration from the unsaturated zone to water table.

2. Formation of Mathematical Model of the Groundwater Flow

In nature, groundwater is in very large land scale; to understand this problem, we form a mathematical model such that the homogenous porous media are laid on impervious inclined layer with θ angle as shown in Figure 1. If the channels 1 and 2 on boundaries are filled with water to heights h 1 and h 2 , respectively, with respect to the datum line, then the system will reach steady state for such fluid flow problem.

This figure represents unconfined aquifer on incline impervious boundary with positive slope.

If the flow rates in the direction of x and y are denoted by q x and q y , respectively, in the mathematical model then the continuity equation can be written as (2) q x x + q y y = 0 .

But here we have considered only one-dimensional fluid flow in mathematical model (considering fluid flow in x direction only); then, (3) q x x = 0 q x = c 1 .

By using Dupuit’s assumption for the model we get (4) - K ( h - m x ) d h d x = c 1 .

By rewriting the above equation, (5) c 1 d x + K ( h - m x ) d h = 0 .

Equation (5) represents a mathematical model in which ordinary differential equation arises during one-dimensional groundwater recharge phenomena. It is an ordinary differential equation which is not exact; now to convert it into exact, we multiply (5) with integrating factor e ( K m / c 1 ) h : (6) c 1 e ( K m / c 1 ) h d x + K e ( K m / c 1 ) h ( h - m x ) d h = 0 .

As per Figure 1, we choose appropriate boundary conditions: (7) h ( x = 0 ) = h 1 , h ( x = l ) = h 2 .

The solution of (6) is given by (8) ( m x - h + c 1 K m ) e ( K m / c 1 ) h = c 2 .

To determine constants of integration c 1 and c 2 by using boundary conditions equation (7), we get (9) c 1 = h 2 K m - K m 2 + ( c 1 - K m h 1 ) e ( K m / c 1 ) ( h 1 - h 2 ) , (10) c 2 = ( c 1 K m - h 1 ) e ( K m / c 1 ) h 1 .

Three cases have been considered in analysis and consistency of the mathematical model and its solution with suitable boundary conditions and different slopes of impervious incline boundary.

Case 1.

In first case, let the total heads of channels 1 and 2 be 600  cm and 300  cm, respectively, for the mathematical model. The distance between two channels is l = 1000  cm and the hydraulic conductivity of the soil is K = 1 0 - 4  cm/sec. Now, if we take positive slope of impervious boundary m = 0.2 , we can get constants and c 1 = - 0.00973 and c 2 = - 316.768 . In the same way for negative slope m = - 0.2 we can get c 1 = - 0.01659 and c 2 = 473.082 . After substituting the values of and in (8) we get the two curves of free surfaces for the top flow lines as shown in Figure 2, such that the upper curve represents free surface for positive slope m = 0.2 and the lower curve represents free surface of water for negative slope m = - 0.2 .

This figure represents the free surface of water table: h 1 > h 2 for m > 0 and m < 0 .

Case 2.

In the second case, let the total heads of channels 1 and 2 be 200  cm and 600  cm, respectively, for this mathematical model. The other parameters are same as Case 1. Now for the positive slope we get c 1 = 0.01242 and c 2 = 581.305 and for negative slope c 1 = 0.01917 and c 2 = - 940.203 . Using these constants in (8) gives us the upper flow curves of free surface of water as shown in Figure 3, such that the lower curve represents free surface for positive slope m = 0.2 and the upper curve represents free surface of water for negative slope m = - 0.2 .

This figure represents the free surface of water table: h 1 < h 2 for m > 0 and m < 0 .

Case 3.

This case is considered for the same pressure head of both channels 1 and 2 for this mathematical model. Therefore, h 2 = h 1 + m l . Substituting this value in (9) we get (11) c 1 = K m h 1 + ( c 1 - K m h 1 ) e - ( K m 2 l / c 1 ) .

The numerical value of c 1 for positive slope m = 0.2 in (11) is 0.006 , which is same as K m h 1 ; that is, c 1 = K m h 1 = 0.006 and c 2 = 0 . Using the values of c 1 and c 2 in (8) we get (12) h = m x + h 1 .

Equation (12) is an equation of straight line having h intercept as h 1 with positive slope m = 0.2 , and for negative slope m = - 0.2 , c 1 = 0.00972 and c 2 = - 424.33 .

Figure 4 shows that the upper flow line represents free surface of water for positive slope and is parallel to the inclined impervious boundary.

This figure represents the free surface of water table for equal pressure head and m > 0 and m < 0 .

3. Mathematical Analysis of Flow Rate

Differentiates (8) with respect to x , which gives change in free surface with respect to distance x . Hence, (13) d h d x = m 1 - ( c 2 K m / c 1 ) e - ( K m / c 1 ) h .

Also, Darcy’s law gives the flow rate of the system as (14) q = - K ( h - m x ) d h d x , q = K m ( m x - h ) 1 - ( c 2 K m / c 1 ) e - ( K m / c 1 ) h .

Thus, (14) gives flow rate expressed in the terms of distance x and the negative exponential function of height h of free surface of water.

4. Numerical and Graphical Interpretation of Mathematical Model

In the first case, as distance x increases, height h decreases; it shows that the water level of groundwater decreases more speedily when impervious boundary has negative slope compared to positive slope of the impervious boundary, which is consistent physically as well as mathematically and also holds natural area for recharge. The second case is the complete opposite of the first case. In the second case, according to boundary condition, as distance x increases, height h also increases, but the increase rate of free water surface for negative slope is faster than positive slope of impervious boundary. The boundary condition in the third case is set as the heights of pressure heads of both channels are the same. Therefore, for the positive slope the height of free surface of water increase is parallel to the impervious incline boundary.

5. Conclusions

The present mathematical model has referenced the method of Hantush and Cruz  and Zamani  because in their approach the equation of upper flow line has a directional relationship to the average depth of free water surface with respect to the impervious boundary. Equation (14) is the solution of the one-dimensional flow rate of infiltrate water in incline porous media at any distance x and height h . Looking at all the cases, it can be concluded that the level of free surface of water table can be measured at different distances.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are thankful to Professor M. N. Mehta for his valuable suggestions and the Department of Applied Mathematics and Humanities, S. V. National Institute of Technology, Surat, for providing research facilities.

Patel K. Mehta M. Patel T. A series solution of moisture content in vertical groundwater flow through unsaturated heterogeneous porous media International eJournal of Mathematics and Engineering 2012 159 1 1467 1477 Swaroop A. Mehta M. N. A solution to the problem of one dimensional flow in unsaturated porous media taking Finite Element Approach Proceedings of the International Conference on Mathematical Modeling January 2001 Roorkee, India 141 143 Philips R. Advances in Hydro Science 1970 New York, NY, USA Academic Press edited by V. Chow Verma A. Mishra S. A similarity solution of a unidimensional vertical groundwater recharge through porous media Revue Roumaine des Sciences Techniques Série de Mécanique Appliquée 1973 18 2 345 351 2-s2.0-0015569486 Mehta M. N. A singular perturbation solution of one-dimensional flow in unsaturated porous media with small diffusivity coefficient Proceedings of the National Conference on Fluid Mechanics and Fluid Power (FMFP '75) 1975 E1 E4 Bear J. The Hydraulics of Groundwater 1979 New York, NY, USA McGraw-Hill Todd D. K. Mays L. W. Groundwater Hydrology 2005 5th New York, NY, USA John Wiley & Sons Hantush M. M. Cruz J. Hydrogeologic foundations in support of ecosystem restoration United States Enviromental Protection Agency 1999 EPA/600/R-99/104 Bear J. Alexander H. Modeling Groundwater Flow and Contaminant Transport Dynamics of Fluids in Porous Media 2010 New York, NY, USA Springer Zamani M. An analysis on groundwater flow in inclined porous media International Journal of Science and Advanced Technology 2012 147 151