Math. J. Interdiscip. Sci.

Approximate Analytical Solution of Advection-Dispersion Equation By Means of OHAM

Dipak J Prajapati and N B Desai

  • Download PDF
  • DOI Number

Advection, Dispersion, Convergence-control parameter, Discrete squared residual

PUBLISHED DATE September, 2018
PUBLISHER The Author(s) 2018. This article is published with open access at

This work deals with the analytical solution of advection dispersion equation arising in solute transport along unsteady groundwater flow in finite aquifer. A time dependent input source concentration is considered at the origin of the aquifer and it is assumed that the concentration gradient is zero at the other end of the aquifer. The optimal homotopy analysis method (OHAM) is used to obtain numerical and graphical representation.


Solutions of advection-dispersion equation (ADE) may be used to predict the concentration of solutes in unsteady groundwater flow. Advection causes the contaminant plum to flow in the direction of groundwater water without any change in the shape. Dispersion of plum arises due to the variation in groundwater velocity. The heterogeneity of the porous medium is responsible for dispersion. The solute transport in heterogeneous aquifer is thus the combined process of advection and dispersion. The dispersion in the direction of groundwater flow is called longitudinal dispersion and the transverse dispersion is perpendicular to the groundwater flow direction. The ADE can be derived using Fick’s law and low of conservation of mass.

Analytical solutions in one-dimensional problems through semi-infinite or finite porous media have been presented by several researchers: (Mazaheri et. al. 2013, Kumar et. al. 2010, Marino et al. 1974, Singh et al. 2008) etc. The objective of this work is to derive an approximate analytical solution of ADE with the help of Optimal Homotopy Analysis Method (OHAM). In this work, an approximate analytical solution of one-dimensional ADE in heterogeneous finite aquifer is derived for continuous time dependent input source concentration of increasing nature.

Page(s) 15-20
ISSN Print: 2278-9561, Online: 2278-957X

The approximate analytical solution of advection-dispersion equation with variable coefficients is obtained by optimal homotopy analysis method with time-dependent input source concentration. The contaminant transport behaves as expected i.e. the concentration at a fixed time decreases as the distance increases and the concentration at a fixed position advances with time. The solution is useful as a preliminary predictive tool for simulating the solute migration in aquifer due to the release of a time-dependent source.

  • Baxter M., Van Gorder R. A., Vajravelu K (2014). On the choice of auxiliary linear operator in the optimal homotopy analysis of the Cahn-Hilliard initial value problem, Numer Algor, 66(2) 269–298.
  • Das P., Begam S., Singh M. K. (2017). Mathematical modeling of groundwater contamination with varying velocity field, Journal of Hydrology and Hydromechanics, 65(2) 192–204.
  • Kumar A., Jaiswal D. K., Kumar A. (2010). Analytical solutions to one-dimensional advection diffusion equation with variable coefficients in semi-infinite media, Journal of Hydrology, 380, 330–337.
  • Liao S. J. (2010). An optimal homotopy analysis approach for strongly nonlinear differential equations, Commun. Nonlinear Sci. Numer. Simul., 15, 2003–2016.
  • Marino M. A. (1974). Distribution of contaminants in porous media flow, Water Resources Research, 10(5).
  • Mazaheri M., Samani J. M. V., Samani H. M. V. (2013). Analytical solution to one-dimensional advection-diffusion equation with several point sources through arbitrary time-dependent emission rate patterns, J. Agr. Sci. Tech., 15, 1231–1245.
  • Prajapati D. J., Desai N. B. (2015). The solution of immiscible fluid flow by means of optimal homotopy analysis method, International Journal of Computer and Mathematical Sciences, 4(8).
  • Prajapati D. J., Desai N. B. (2016). Application of the basic optimal homotopy analysis method to fingering phenomenon, Global Journal of Pure and Applied Mathematics, 12(3), 2011–2022.
  • Prajapati D. J., Desai N. B. (2017). Optimal homotopy analysis solution of fingero-imbibition phenomenon in homogeneous porous medium with magnetic fuid effect, Kalpa Publications in Computing, 2, 85–94.
  • Prajapati D. J., Desai N. B. A. (2017). pproximate analytical solution for the fingero-imbibition phenomenon by optimal homotopy analysis method, International Journal of Computational and Applied Mathematics, 12(3), 751–761.
  • Prajapati D. J., Desai N. B. (2017). Analytic analysis for oil Recovery during cocurrent imbibition in inclined homogeneous porous medium, International Journal on Recent and Innovation Trends in Computing and Communication, 5(7), 189–194.
  • Prajapati D. J., Desai N. B. (2017). Mathematical investigation of counter-current imbibition phenomenon in heterogeneous porous medium, PRAJÑ? - Journal of Pure and Applied Sciences, 24-25 38–46.
  • Singh M. K., Singh P., Singh V. P. (2010). Analytical solution for solute transport along and against time dependent source concentration in homogeneous finite aquifers, Advances in Theoretical and Applied Mechanics, 3(3), 99–119.
  • Singh M. K., Mahto N. K., Singh P. (2008). Longitudinal dispersion with time-dependent source concentration in semi-infinite aquifer, Journal of Earth System Science, 117(6) 945–949.
  • Scheidegger A. E. (1961). General theory of dispersion in porous media, Journal of Geophysical Research, 66(10) 3273–3278.
  • Xu M., Eckstein Y. (1997). Statistical analysis of the relationships between dispersivity and other physical properties of porous media. Hydrogeology Journal, 5(4) 4–20.