ISSN ONLINE(2319-8753)PRINT(2347-6710)

Homotopy Analysis Method Applied To Determine Pressure Head In Unsaturated Soil During Infiltration Phenomenon

Kajal K. Patel1, M. N. Mehta2, Twinkle R. Singh3
  1. Research Scholar, Department of Applied Mathematics & Humanities, S. V. National Institute of Technology, Surat, Gujarat, India
  2. Professor, Department of Applied Mathematics & Humanities, S. V. National Institute of Technology, Surat, Gujarat, India
  3. Assistant Professor, Department of Applied Mathematics & Humanities, S. V. National Institute of Technology, Surat, Gujarat, India
Related article at Pubmed, Scholar Google

Visit for more related articles at International Journal of Innovative Research in Science, Engineering and Technology


The Homotopy analysis method has been applied to Richard's equation under considering the value of diffusivity coefficient D  and hydraulic conductivity K  under certain assumption. The solution of this equation has been obtained by considering the guess value of pressure head in an unsaturated soil. It is concluded that the pressure head in an unsaturated soil is decreasing as depth Z is increasing for given time T > 0. The graph and a numerical value are given by using Maple coding


Richard's equation, Unsaturated porous media, Homotopy Analysis Method, Infiltration Phenomenon.


From last few years, groundwater flow and mass transport in the unsaturated porous media has significantly increasingly and it has a great important for hydrologist, agriculturists and people related with water resources sciences. In hydrology, volumetric water content is very important phenomenon in infiltration of groundwater in an unsaturated porous media. The groundwater flow and mass transport in the unsaturated porous media is a region of underground as its top by the soils (porous media) and below by the groundwater table, where physical phenomena such as infiltration, groundwater recharge, moisture content and others are happened. Ground water recharge problem has been discussed by many researchers with different viewpoints. Swartzendruber uses Philip’s [31] method to get graphical illustration of the mathematical solution for horizontal water function. Verma and Mishra [29] have obtained solution by similarity transformation of a unidimensional vertical ground water recharges through porous media. Mehta [27] has obtained an approximate solution by the method of singular perturbation technique. Hari Prasad et al. [30] had developed a numerical model for unsaturated zones for the process of gravity drainage and infiltration. Almost all unsaturated flow simulations use the Richard's equation. Richards [15] investigated the influence capillarity has on liquid infiltration in soils, and equation derived is now known as Richards' equation. He derived a governing nonlinear equation for water flow in porous media based on a continuum approach. In the derivation of Richards' equation, it is assumed the liquid infiltration is driven by capillarity and gravity, while the second fluid phase occupying the unsaturated void-spaces is largely inert and does not affect the evolution. The sophistication of the analytical and numerical methods that are available for solving the governing nonlinear equation of unsaturated flow in porous media (Richard's Equation) require specialized forms of the conductivity (hydraulic conductivity) and water diffusivity (diffusivity functions). Several investigators have described different relation between the diffusivity coefficient and hydraulic conductivity in groundwater recharge problem.
In 1958, Gradner [32] model provided a relationship between the pressure head 'h' and the volumetric water content  as,   b h  a     where a & b are empirical constants. The exponential hydraulic conductivity function has been widely used, but it is known to have a limited range of application to many real soils. Other functions developed by Brooks and Corey [33] and Van Genuchten [31] are firmly established for practical applications. Such special forms of the hydraulic functions make it possible to linearize the governing flow equations, and hence solve them analytically. Solutions to the linearized unsaturated flow equations are generally limited to the steady flow in semi-infinite, homogenous soils (Broadbridge and White, 1988; Warrick and Yeh, 1990) and to transient flow in homogeneous and layered soils (Srivastava and Yeh) [7]. In 1988,The Broadbridge and White Model [5] has adopted a function form for the diffusivity given by Philip and Knight (1974)[6] which allows for the transformation of the soil water diffusivity D  as function has the form:
conditions, the Hopf-Cole transformations are applied to reduce a nonlinear equation to a linear form that possesses an exact parametric solution [5]. In this paper, Broadbridge and White Model [5] was applied to solve infiltration phenomenon in unsaturated porous media. In 1992, Liao[12] employed the basic idea of the homotopy in topology to propose a general analytical method for nonlinear problems, namely Homotopy Analysis method. Homotopy analysis method has been applied to nonlinear fluid dynamics problems by Liao [9], [11]. The important factor of this method is the generalized Taylor expansion. Liao point out that the generalized Taylor series provided a way to control and adjust the convergence region through an auxiliary parameter such that this method is particularly suitable for strongly nonlinear problems. This method does not dependence on small or large parameter and is easy to adjust the convergence region and rate of approximation of the series. In this paper, ground water recharge phenomenon is solved by Homotopy analysis method.
In the present model, Darcian-based unsaturated flow equation described over the large basin L (i.e. L is the length of a large basin) of homogeneous porous media, and there is no air resistance to the flow (i.e. the porous medium contains only the flowing liquid water and empty voids of air). The air in the void space is approximately at atmospheric pressure. The flowing liquid (water) is considered continuous at a microscopic level, incompressible and isothermal, where the pressure head in unsaturated soil is considered as time depended function. In this case consider the flow of infiltrated water in vertical downward direction up to length L by neglecting spreading in other directions. We obtain one dimensional nonlinear partial differential equation by combining Darcy's law [1] for unsaturated flow with the continuity equation.
In this paper, the Homotopy analysis method is implemented to obtain the approximate solutions of the governing nonlinear partial second order differential equation for infiltration phenomenon, which is able to solve the onedimensional h-based forms of Richard's equation. Its solution provides the infiltrated water into the porous media (soil) at any depth Z at time T > 0. The graph of pressure head versus time T is given at depth Z > 0, which shows at any depth Z pressure head is decreases as time T increases.




[1] Bear J. and Bachmat Y., "Introduction to Modeling of Transport Phenomena in Porous Media", Kluwer, Academic Publishers, Dordrecht, The Netherlands, 1990.

[2] Bear J.," Dynamics of Fluids in Porous Media", American Elsevier, New York, 1972.

[3] Broadbridge P. and White I., "Constant rate rainfall infiltration: A versatile nonlinear model", I, Analytical solution, Water Resources Research, Vol. 24, No. 1, 145-154, 1988.

[4] Brooks R. H. and Corey A. T., "Hydraulic Properties of Porous Media", Hydrology Paper 3, Coloardo State Uni., Fort Collins, CO., 1964.

[5] Darcy H., "Les Fontaines publiques" de loville de Dijon, Dalmont, Paris, 1856.

[6] Gardner W. R., "Some steady state solutions of unsaturated moisture flow equations with application to evaporation from a water table", Soil sci., 85, 228-232, 1958.

[7] Hari Prasad, K. S., Mohan Kumar and M. Sekhar, "Modelling Flow through unsaturated Zone: Sensitivity to unsaturated Soil Properties", Sadhana, Indian Academy of Sciences, 26 (6),517-528, 2001.

[8] Liao S. J., "A kind of approximate solution technique which does not depend upon small parameters (II): an application in fluid mechanics", Int. J. Nonlinear. Mech, 32-815, 1997.

[9] Liao S. J., "An approximate solution technique which does not depend upon small parameters: a special example", Int. J. Non-linear Mech. 30, 371-380, 1995.

[10] Liao S. J.," Beyond Perturbation: introduction to Homotopy Analysis Method", CRC Press: Chapman & Hall; 2003.

[11] Liao S. J., "On the homotopy analysis method for nonlinear problems", appl. Math. Comput. 147, 499-513, 2004.

[12] Liao S. J.," The proposed homotopy analysis technique for the solution of nonlinear problems", Ph. D. Thesis. Shanghaihao tong University, 1992.

[13] M. N. Mehta, "A singular perturbation solution of one-dimensional flow in unsaturated porous media with small diffusivity coefficient", Proc. FMFP E1 to E4, 1975.

[14] Meher R., Mehta M. N. and. Meher S. K., "Adomian decomposition method for moisture content in one dimensional fluid flow through unsaturated porous media", International Journal of Applied Mathematics and mechanics, 6 (7), 13-23, 2010.

[15] Philips J. R. and Knight J. H.," On solving the unsaturated flow equation, 3, New quasi-analytical technique", soil science, Vol. 117, 1-13, 1974.

[16] Richards L. A., "Physics 1", 318-333, 1931.

[17] Scheidegger A. E.," The physics of flow through Porous media", Third Edition, University of Toronto Press, Toronto, 1974.

[18] Srivastava R. and Yeh J. T.-C., "Analytical solutions for one-dimensional", Transient Infiltration toward the water table in homogenous and Layered soils, Wat. Resour. Res., 27, 753-762, 1991.

[19] Van Genuchten, M Th., "A closed-form Equation for orediciting the HyraulicConducitivity of unsaturated soils", soil sci. Soc., Am. J., 44, 892-898, 1980.

[20] Verma A. P. and Mishra S. K., "A similarity solution of a one-dimensional vertical groundwater recharge through porous media", Rev. Roum. Scien. Techn. Ser. Mec.Appl. 18 (2), 345-351, 1973.

[21] Warrick A. W., Islas A., Lomen D. O., "An analytical solution to Richards equation for time varying infiltration", Water Resour., Res. 27, 763-766, 1991.