Numerical Simulation of Groundwater Table Falling in Horizontal and Sloping Aquifers by Differential Quadrature Method (DQM)
Publication: Journal of Hydrologic Engineering
Volume 17, Issue 8
Abstract
Since a nonlinear partial differential equation was developed by Boussinesq based on Darcy’s law and the Dupuit-Forchheimer assumption, it has played an essential role in the development of simulation models for the solution of various flow problems in porous media. To solve Boussinesq’s equation, both analytical and numerical solutions, such as finite difference (FD) or finite element (FE), have been sought. The differential quadrature method (DQM) is a new numerical method frequently used by researchers to solve partial differential equations. In this paper, DQMs have been coupled with explicit, implicit, and Crank-Nicholson FD and applied to solve Boussinesq’s equation for a drainage problem. This case has been solved previously both analytically and numerically, including by DQM, by other researchers. Those who had employed DQM had linearized Boussinesq’s equation first and solved the linearized form implicitly. In this work, the nonlinear form of Boussinesq’s equation has been solved using three models based on DQM for solving the dimensionless form and on one approach for solving the dimensional form of Boussinesq’s equation. The obtained results and their degree of accuracy are compared with available experimental and numerical data found in the literature, and on that basis, it is concluded that DQM generates accurate results, is very easy to formulate and operate, does not need large mesh size, and is very time efficient.
Get full access to this article
View all available purchase options and get full access to this article.
References
Bear, J., and Cheng, A. H.-D. (2010). Modeling groundwater flow and contaminant transport, Springer, London.
Bellman, R. E., and Casti, J. (1971). “Differential quadrature and long-term integration.” J. Math., 34, 235–238.
Bellman, R. E., Kashef, B. G., and Casti, J. (1972). “Differential quadrature: A technique for the rapid solution of nonlinear partial differential equations.” J. Comput. Phys.JCTPAH, 10(1), 40–52.
Bert, C. W., and Malik, M. (1996). “Differential quadrature method in computational mechanics, a review.” Appl. Mech. Rev.AMREAD, 49(1), 1–27.
Boussinesq, J. (1877). “Du mouvement non permanent des eaux souterraines.” Essai sur le theorie des courantes, Imprimerie Nationale, Paris, 680 (in French).
Boussinesq, J. (1904). “Recherches theoretiques sur l’ecoulement des nappes d’eau infiltrees dans le sol et sur le debit des sources.” J. de Math. Pures Appl.JMPAAM, 5, 5–78 (in French).
Cabrera, G., and Matthey, G. (1984). “Space-time finite element model for two stream unconfined aquifer systems.” J. Hydrol. (Amsterdam)JHYDA7, 73(3–4), 299–314.
Chen, C. N. (2004). “DQEM and DQFDM irregular elements for analyses of 2-D heat conduction in orthotropic media.” Appl. Math. Modell.AMMODL, 28(7), 617–638.
Chen, Z. X., Bodvarsson, G. S., Witherspoon, P. A., and Yortsos, Y. C. (1995). “An integral equation formulation for the unconfined flow of groundwater with variable inlet conditions.” J. Transp. Porous MediaTPMEEI, 18(1), 15–36.
Chu, S. T., and DeBoer, D. W. (1976). “Field and laboratory evaluation of bi-level drainage theory.” Trans. ASABETARSBB, 19(3), 478–481.
Ding, H., Shu, C., Yeo, K. S., and Lu, Z. L. (2005). “Simulation of natural convection in eccentric annuli between a square outer cylinder and a circular inner cylinder using local MQ-DQ method.” Numer. Heat Transfer, Part A, 47(3), 291–313.
Ding, H., Shu, C., Yeo, K. S., and Xu, D. (2004). “Numerical computation of three-dimensional incompressible viscous flows in the primitive variable form by local multiquadric differential quadrature method.” J. Comput. Methods Appl. Mech. Eng., 195, 516–533.
Dogrul, E. C., and Kadir, T. N. (2006). “Flow computation and mass balance in Galerkin finite-element groundwater models.” J. Hydrol. Eng.JHYEFF, 132(11), 1206–1214.
Erduran, K. S., Erduran, V. C., and Macalister, C. R. (2005). “Finite volume solution to integrated shallow surface-saturated groundwater flow.” Int. J. Numer. Methods FluidsIJNFDW, 49(7), 763–783.
Guo, W. (1997). “Transient groundwater flow between reservoirs and water-table aquifers.” J. Hydrol. (Amsterdam)JHYDA7, 195(1–4), 370–384.
Hashemi, M. R., Abedini, M. J., and Malekzadeh, P. (2006). “Numerical modeling of long waves in shallow water using incremental differential quadrature method.” J. Ocean Eng. OCENBQ, 33(13), 1749–1764.
Hashemi, M. R., Abedini, M. J., and Malekzadeh, P. (2007). “A differential quadrature analysis of unsteady open channel flow.” Appl. Math. Modell., 31(8), 1594–1608.
Hashemi, M. R., Abedini, M. J., Neill, S. P., and Malekzadeh, P. (2008). “Tidal and surge modelling using differential quadrature: A case study in the Bristol Channel.” Coastal Eng.COENDE, 55(10), 811–819.
He, X., and Ren, L. (2005). “Finite volume multiscale finite element method for solving the groundwater flow problems in heterogeneous porous media.” Water Resour. Res.WRERAQ, 41, W10417.
Hornberger, G. M., Ebert, J., and Remson, I. (1970). “Numerical solution of the Boussinesq equation for aquifer-stream interaction.” Water Resour. Res.WRERAQ, 6(2), 601–608.
Jain, M. K., Iyenger, S. R. K., and Jain, R. K. (1994). Computational methods for partial differential equations, Wiley Eastern, New Delhi, India.
Kalaidzidou-Paikou, N., Karamouzis, D., and Moraitis, D. (1997). “A finite element model for the unsteady groundwater flow over sloping beds.” Water Resour. Manage.WRMAEJ, 11(1), 69–81.
Karakostas, Z. C., and Manolis, G. D. (1998). “A stochastic boundary element solution applied to groundwater flow.” Eng. Anal. Boundary Elem.EABAEL, 21(1), 9–21.
Karami, G., and Malekzadeh, P. (2002). “A new differential quadrature methodology for beam analysis and the associated differential quadrature element method.” Comput. Methods Appl. Mech. Eng.CMMECC, 191(32), 3509–3536.
Karami, G., and Malekzadeh, P. (2003). “Application of a new differential quadrature methodology for free vibration analysis of plates.” Int. J. Numer. Methods Eng.IJNMBH, 56(6), 847–868.
Kaya, B., and Arisoy, Y. (2010). “Differential quadrature method for linear long wave propagation in open channels.” Wave propagation in materials for modern applications, Petrin, A., ed., InTech, Rijeka, Croatia, 253–266.
Kaya, B., and Arisoy, Y. (2011). “Differential quadrature solution for one-dimensional aquifer flow.” Math. Comput. Appl., 16(2), 524–534.
Kaya, B., Arisoy, Y., and Ulke, A. (2010). “Differential quadrature method (DQM) for numerical solution of the diffusion wave model.” J. Flood Eng., 1(2) 133–147.
Khebchareon, M., and Saenton, S. (2005). “Finite element solution for 1-D groundwater flow, advection-dispersion and interphase mass transfer. I: Model development.” Thai J. Math., 3(2), 223–240.
Kim, D. J., and Ann, M. J. (2001). “Analytical solutions of water table variation in a horizontal unconfined aquifer: Constant recharge and bounded by parallel streams.” Hydrol. ProcessesHYPRE3, 15(13), 2691–2699.
Lesnic, D., Elliott, L., and Ingham, D. B. (1998). “A boundary element method for the determination of the transmissivity of a heterogeneous aquifer in groundwater flow systems.” Eng. Anal. Boundary Elem.EABAEL, 21(3), 223–234.
Li, L., Zhou, H., and Gomez-Hernandez, J. J. (2010). “Steady-state saturated groundwater flow modeling with full tensor conductivities using finite differences.” Comput. Geosci.CGEODT, 36(10), 1211–1223.
Lin, C. L. (1972). “Digital simulation of the Boussinesq equation for a water table aquifer.” Water Resour. Res.WRERAQ, 8(3), 691–698.
Lockington, D. A. (1997). “Response of unconfined aquifer to sudden change in boundary head.” J. Irrig. Drain. Eng.JIDEDH, 123(1), 24–27.
Lockington, D. A., Parlange, J. Y., Parlange, M. B., and Selker, J. (2000). “Similarity solution of the Boussinesq equation.” Adv. Water Resour.AWREDI, 23(7), 725–729.
Malekzadeh, P., and Karami, G. (2003). “Out-of-plane static analysis of circular arches by DQM.” Int. J. Solids Struct.IJSOAD, 40(23), 6527–6545.
Malekzadeh, P., and Karami, G. (2004). “Vibration of non-uniform thick plates on elastic foundation by differential quadrature method.” J. Eng. Struct., 26(10), 1473–1482.
Malekzadeh, P., and Karami, G. (2005). “Polynomial and harmonic differential quadrature method for free vibration of variable thickness thick skew plates.” J. Eng. Struct., 27(10), 1563–1575.
Mehl, S., and Hill, M. C. (2004). “Three-dimensional local grid refinement for block-centered finite-difference groundwater models using iteratively coupled shared nodes: A new method of interpolation and analysis of errors.” Adv. Water Resour.AWREDI, 27(9), 899–912.
Moutsopoulos, K. N. (2010). “The analytical solution of the Boussinesq equation for flow induced by a step change of the water table elevation revisited.” J. Transp. Porous MediaTPMEEI, 85(3), 919–940.
Onyejekwe, O. O. (1999). “Subsurface drainage of sloping land.” Eng. Anal. Boundary Elem., 23(7), 619–624.
Ostfeld, A., Muzaffar, E., and Lansey, K. E. (1999). “Analytical ground-water flow solutions for channel-aquiferinteraction.” J. Irrig. Drain. Eng.JIDEDH, 125(4), 196–202.
Parlange, J. Y., Hogarth, W. L., Govindaraju, R. S., Parlange, M. B., and Lockington, D. (2000). “On an exact analytical solution of the Boussinesq equation.” J. Transp. Porous MediaTPMEEI, 39(3), 339–345.
Polubarinova-Kochina, P. Y.-A. (1962). Theory of groundwater movement, Princeton University Press, Princeton, NJ.
Robati, A., and Barani, G. A. (2009). “Modeling of water surface profile in subterranean channel by differential quadrature method (DQM).” Appl. Math. Modell., 33(3), 1295–1305.
Rushton, K. R., and Chan, Y. (1977). “Numerical pumping test analysis in unconfined aquifers.” J. Irrig. Drain Div.JRCEA4, 103(IR1), 1–12.
Sabti, N. A. (1989). “Linear and nonlinear solution of the Boussinesq equation for the bi-level drainage problem.” Agric. Water Manage.AWMADF, 16(4), 269–278.
Shu, C. (1991). “Generalized differential quadrature-integral quadrature and application to the simulation of incompressible viscous flows including parallel computation.” Ph.D. thesis, Univ. of Glasgow, Glasgow, Scotland.
Shu, C. (1996). “Free vibration analysis of composite laminated conical shells by generalized differential quadrature.” J. Sound Vib.JSVIAG, 194(4), 587–604.
Shu, C. (1999). “Application of differential quadrature method to simulate natural convection in a concentric annulus.” Int. J. Numer. Methods FluidsIJNFDW, 30(8), 977–993.
Shu, C. (2000). Differential quadrature and its application in engineering, Springer, London.
Shu, C., and Chew, Y. T. (1997). “Fourier expansion-based differential quadrature and its application to Helmholtz eigenvalue problems.” Commun. Numer. Methods Eng., 13(8), 643–653.
Shu, C., Chew, Y. T., Khoo, B. C., and Yeo, K. S. (1995). “Application of GDQ scheme to simulate incompressible viscous flows around complex geometries.” J. Mech. Res. Commun.MRCOD2, 22(4), 319–325.
Shu, C., Ding, H., Chen, H. Q., and Wang, T. G. (2005a). “An upwind local RBF-DQ method for simulation of inviscid compressible flows.” Comput. Methods Appl. Mech. Eng.CMMECC, 194(18–20), 2001–2017.
Shu, C., Ding, H., and Yeo, K. S. (2003). “Local radial basis function-based differential quadrature method and its application to solve two-dimensional incompressible Navier-Stokes equations.” J. Comput. Methods Appl. Mech. Eng., 192(7–8), 941–954.
Shu, C., Ding, H., and Yeo, K. S. (2005b). “Computation of incompressible Navier-Stokes equations by local RBF-based differential quadrature method.” Comput. Model. Eng. Sci.CMESCE, 7(2), 195–205.
Shu, C., Khoo, B. C., and Yeo, K. S. (1994a). “Numerical solution of incompressible Navier-Stokes equations generalized differential quadrature.” Finite Elem. Anal. Des., 18(1–3), 83–97.
Shu, C., Khoo, B. C., Yeo, K. S., and Chew, Y. T. (1994b). “Application of GDQ scheme to simulate natural convection in a square cavity.” Int. Commun. Heat Mass Transf.IHMTDL, 21(6), 809–817.
Shu, C., and Richards, B. E. (1990). “High resolution of natural convection in a square cavity by generalized differential quadrature.” Proc., 3rd Conf. on Advances in Numerical Methods in Engineering: Theory and Applications, Vol. 2, Elsevier Applied Science, New York, 978–985.
Shu, C., and Richards, B. E. (1992). “Application of generalized differential quadrature to solve 2-dimensional incompressible Navier-Stokes equation.” Int. J. Numer. Methods FluidIJNFDW, 15(7), 791–798.
Shu, C., and Xue, H. (1997). “Explicit computation of weighting coefficients in the harmonic differential quadrature.” J. Sound Vib.JSVIAG, 204(3), 549–555.
Shukla, K. N., Chauhan, H. S., and Srivastava, V. K. (1990). “Finite difference solution of Boussinesq in unsteady-state equation for highly sloping lands.” J. Irrig. Drain. Eng.JIDEDH, 116(1), 107–113.
Singh, S. R., and Jacob, C. M. (1976). “Numerical solution of Boussinesq equation.” J. Eng. Mech. Div.JMCEA3, 102(5), 807–824.
Singh, R. M., Singh, K. K., and Singh, S. R. (2006). “Falling water tables in a sloping/nonsloping aquifer under various initial water table profiles.” Agric. Water Manage.AWMADF, 82(1-2), 210–222.
Szekely, F. (2008). “Three-dimensional mesh resolution control in finite difference groundwater flow models through boxed spatial zooming.” J. Hydrol. (Amsterdam)JHYDA7, 351(3–4), 261–267.
Taigbenu, A. E. (1986). “An integral formulation applied to the diffusion and Boussinesq equations.” Int. J. Numer. Methods Eng.IJNMBH, 23(6), 1057–1079.
Taigbenu, A. E. (1991). “A simplified finite element treatment of the nonlinear Boussinesq equation.” Adv. Water Resour.AWREDI, 14(1), 42–49.
Taigbenu, A. E. (2003). “A time-dependent Green’s function-based model for stream-unconfined aquifer flows.” Water SAWASADV, 29(3), 273–280.
Telyakovskiy, A. S., Bragab, G. A., and Furtadoa, F. (2002). “Approximate similarity solutions to the Boussinesq equation.” Adv. Water Resour.AWREDI, 25(2), 191–194.
Tolikas, P. K., Sidiropoulos, E. G., and Tzimopolous, C. D. (1984). “A simple analytical solution for the Boussinesq one-dimensional groundwater flow equation.” Water Resour. Res.WRERAQ, 20(1), 24–28.
Upadhyaya, A., and Chauhan, H. S. (1998). “Solutions of Boussinesq equation in semiinfinite flow region.” J. Irrig. Drain. Eng.JIDEDH, 124(5), 265–270.
Upadhyaya, A., and Chauhan, H. S. (2000). “An analytical solution for bi-level drainage design in the presence of evapotranspiration.” Agric. Water Manage.AWMADF, 45(2), 169–184.
Upadhyaya, A., and Chauhan, H. S. (2001a). “Interaction of stream and sloping aquifer receiving constant recharge.” J. Irrig. Drain Eng.JRCEA4, 127(5), 295–301.
Upadhyaya, A., and Chauhan, H. S. (2001b). “Falling water tables in non-sloping/sloping aquifer.” J. Irrig. Drain. Eng.JIDEDH, 127(6), 378–384.
Verhoest, N. E. C., and Troch, P. A. (2000). “Some analytical solutions of the linearized Boussinesq equation with recharge for a sloping aquifer.” Water Resour. Res.WRERAQ, 36(3), 793–800.
Verma, A. K., Gupta, S. K., Singh, K. K., and Chauhan, H. S. (1998). “An analytical solution for design of bi-level drainage systems.” Agric. Water Manage.AWMADF, 37(1), 75–92.
Wang, X., Gan, L., and Zhang, Y. (2008). “Differential quadrature analysis of the buckling of thin rectangular plates with cosine-distributed compressive loads on two opposite sides.” J. Adv. Eng. Software, 39(6), 497–504.
Zissis, T. S., Teloglou, I. S., and Terzidis, G. A. (2001). “Response of a sloping aquifer to constant replenishment and to stream varying water level.” J. Hydrol. (Amsterdam)JHYDA7, 243(3–4), 180–191.
Zucker, M. B., Remson, I., Ebert, J., and Aguado, E. (1973). “Hydrologic studies using the Boussinesq equation with a recharge term.” Water Resour. Res.WRERAQ, 9(3), 586–592.
Information & Authors
Information
Published In
Journal of Hydrologic Engineering
Volume 17 • Issue 8 • August 2012
Pages: 869 - 879
Copyright
© 2012. American Society of Civil Engineers.
History
Received: Jan 16, 2011
Accepted: Sep 30, 2011
Published online: Oct 6, 2011
Published in print: Aug 1, 2012
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.