Skip to main content
Log in

Numerical Solution and Analysis of Three-Dimensional Transient Darcy Flow

  • Published:
Transport in Porous Media Aims and scope Submit manuscript

Abstract

This paper presents a detailed analysis of a numerical solution of three-dimensional transient Darcy flow. The numerical solution of the governing parabolic partial differential equations is obtained by using stabilized mixed Galerkin method and backward Euler method for the discretization of space and time, respectively. The resulting well-posed system of algebraic equations is subsequently solved using conjugate gradient method. The proposed model is validated against Mongan’s analytical model for underground water flow using a set of hexahedral and tetrahedral meshes. The model is used to analyze the transient behavior by simulating the Darcy flow through homogeneous and heterogeneous as well as isotropic and anisotropic media. For large meshes, a parallel algorithm of the transient Darcy flow is also developed for shared memory architecture using OpenMP library. For structured meshes, a speedup of over 22 is obtained on dual AMD Opteron processors. The proposed numerical method for transient Darcy flow offers stability, ease of implementation in higher dimensions and parallel solution for large and complex geometry using standard finite element spaces.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10

Similar content being viewed by others

References

  • Ansari, S.U., Hussain, M., Mazhar, S., Rashid, A. Ahmad, S.M.: Parallel stabilized mixed Galerkin method for three-dimensional Darcy flow using OpenMP. In: National Software Engineering Conference (NSEC 2015), Rawalpindi, Pakistan (2015a)

  • Ansari, S.U., Hussain, M., Mazhar, S., Rashid, A., Ahmad, S.M.: Stabilized mixed Galerkin method for transient analysis of Darcy flow. In: International Conference on Modeling, Simulation and Applied Optimization (ICMSAO’15), Istanbul, Turkey (2015b)

  • Ansari, S.U., Hussain, M., Mazhar, S., Rashid, A., Ahmad, S.M.: Three-dimensional stabilized mixed Galerkin method for Darcy flow. In: 13th International Conference on Frontiers of Information Technology (FIT2015), Islamabad, pp. 104–108 (2015c). https://doi.org/10.1109/fit.2015.67

  • Ansari, S.U., Hussain, M., Ahmad, S. M., Rashid, A., Mazhar, S.: Stabilized mixed finite element method for transient Darcy flow (to appear in Transaction of the Canadian Society for Mechanical Engineering) (2017)

  • Brezzi, F., Russo, A.: Choosing bubbles for advection–diffusion problems. Math. Models Methods Appl. Sci. 4, 571–587 (1994)

    Article  Google Scholar 

  • Brooks, A.N., Hughes, T.J.R.: Streamline upwind/Petrov–Galerkin methods for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations. Comput. Methods Appl. Mech. Eng. 32, 199–259 (1982)

    Article  Google Scholar 

  • D’Angelo, C., Zunino, P.: Numerical approximation with Nitsche’s coupling of Stokes’/Darcy’s flow problems applied to hemodynamics. Appl. Numer. Math. 62, 378–395 (2012)

    Article  Google Scholar 

  • Darcy, H.: Les Fontaines Publiques de la Ville de Dijon. Dalmont, Paris (1856)

    Google Scholar 

  • El-Amin, M.F., Salama, A., Sun, S.: A conditionally stable scheme for a transient flow of a non-Newtonian fluid saturating a porous medium. In: International Conference on Computer Science (ICCS2012), Procedia Computer Science, vol. 9, pp. 651–660 (2012)

  • Franca, L.P., Stenberg, R.: Error analysis for some Galerkin-least-squares methods for the elasticity equations. SIAM J. Numer. Anal. 28(6), 1680–1797 (1991)

    Article  Google Scholar 

  • Furman, A., Neuman, S.P.: Laplace-transform analytic element solution of transient flow in porous media. Adv. Water Resour. 26, 1229–1237 (2003)

    Article  Google Scholar 

  • Huang, T., Guo, X., Chen, F.: Modeling transient pressure behavior of a fractured well for shale gas reservoirs based on the properties of nanopores. J. Nat. Gas Sci. Eng. 23, 387–398 (2015)

    Article  Google Scholar 

  • Hughes, T.J.R.: The Finite Element Method: Linear Static and Dynamic Finite Element Analysis. Prentice-Hall, Englewoods Cliffs, NJ (1987). (Dover edition, 2000)

    Google Scholar 

  • Hughes, T.J.R., Brooks, A.N.: An algorithm for solving the Navier–Stokes equations based upon the streamline-upwind/Petrov–Galerkin formulation. In: Lewis, R.W., et al. (eds.) Numerical Methods in Coupled Systems, pp. 387–404. Wiley, London (1984)

    Google Scholar 

  • Hughes, T.J.R., Franca, L.P., Balestra, M.: A new finite element formulation for computational fluid dynamics. V. Circumventing the Babuska-Brezzi condition: a stable Petrov–Galerkin formulation of the stokes problem accommodating equal-order interpolations. Comput. Methods Appl. Mech. Eng. 59, 85–99 (1986)

    Article  Google Scholar 

  • Hughes, T.J.R., Franca, L.P., Hulbert, G.M.: A new finite element formulation for computational fluid dynamics. VIII. The Galerkin/least-squares method for advective–diffusive equations. Comput. Methods Appl. Mech. Eng. 73, 173–189 (1989)

    Article  Google Scholar 

  • Hussain, M., Kavokin, A.: A calculation of 3D model of ground water flux at evaporation from water table using parallel algorithm—MPICH. Int. J. Math. Phys. 3(2), 128–132 (2012)

    Google Scholar 

  • Hussain, M., Kavokin, A.: A 2D parallel algorithm using MPICH for calculation of ground water flux at evaporation from water table. In proceedings of FIT’09, Abbottabad, Pakistan (2009)

  • Hussain, M., Abid, M., Ahmad, M.: Stabilized mixed finite elements for Darcy’s law on distributed memory systems. In: Proceedings of International Symposium on Frontiers of Computational Sciences, Islamabad, Pakistan, pp. 39–47 (2012)

  • Hussain, M., Abid, M., Ahmad, M., Hussain, S.F.: A parallel 2D stabilized finite element method for Darcy flow on distributed systems. World Appl. Sci. J. 27(9), 1119–1125 (2013)

    Google Scholar 

  • Istok, J.: Groundwater Modeling by the Finite Element Method. Wiley (2013)

  • Jia, H., Li, K., Liu, S.: Characteristic stabilized finite element method for transient Navier-Stokes equations. Comput. Methods Appl. Mech. Eng 199, 2996–3004 (2010)

    Article  Google Scholar 

  • Jiao, J., Zhang, Y.: A method based on local approximate solutions (LAS) for inverting transient flow in heterogeneous aquifers. J. Hydrol. 514, 145–149 (2014)

    Article  Google Scholar 

  • Lewis, R.W., Nithiarasu, P., Seetharamu, K.: Fundamentals of the Finite Element Method for Heat and Fluid Flow. Wiley, Hoboken (2004)

    Book  Google Scholar 

  • Loula, A.F.D, Correa, M.R.: Numerical analysis of stabilized mixed finite element methods for darcy flow. In: III European Conference on Computational Mechanics Solids, Structures and Coupled Problems in Engineering, Lisbon, Portugal (2006)

  • Masud, A., Hughes, T.J.R.: A stabilized mixed finite element method for Darcy flow. Comput. Methods Appl. Mech. Eng. 191, 4341–4370 (2002)

    Article  Google Scholar 

  • Mongan, C.E.: Validity of Darcy’s law under transient conditions. U.S. Geological Survey Professional Paper, vol. 1331 (1985)

  • Noetinger, B.: A quasi steady state method for solving transient Darcy flow in complex 3D fractured networks accounting for matrix to fracture flow. J. Comput. Phys. 283, 204–223 (2015)

    Article  Google Scholar 

  • Philip, J.R.: Transient Fluid Motions in Saturated Porous Media. CSIRO, Canberra (1952)

    Google Scholar 

  • Silva, L., Valette, R., Laure, P., Coupez, T.: A new three-dimensional mixed finite element for direct numerical simulation of compressible viscoelastic flows with moving free surfaces. Int.J. Mater. Form. 5(1), 55–72 (2012)

    Article  Google Scholar 

  • Singh, A.K., Agnihotri, P., Singh, N.P., Singh, A.K.: Transient and non-Darcian effects on natural convection flow in a vertical channel partially filled with porous medium: analysis with Forchheimer–Brinkman extended Darcy model. Int. J. Heat Mass Transf. 54, 1111–1120 (2011)

    Article  Google Scholar 

  • Tezduyar, T.E.: Stabilized finite element formulations for incompressible flow computations. Adv. Appl. Mech. 28, 1–44 (1992)

    Google Scholar 

  • Wu, Y.S., Pruess, K.: Integral solutions for transient fluid flow through a porous medium with pressure-dependent permeability. Int. J. Rock Mech. Min. Sci. 37, 51–61 (2000)

    Article  Google Scholar 

  • Zheng, C., Bennett, G.D.: Applied Contaminant Transport Modeling, vol. 2. Wiley-Interscience, New York (2002)

    Google Scholar 

Download references

Acknowledgements

The authors are very grateful to Dr. Arif Masud in University of Illinois at Urbana-Champaign for valuable guidance and suggestions for this research work. The authors are also obliged to Ghulam Ishaq Khan (GIK) Institute for their financial support in this research work.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shahab U. Ansari.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ansari, S.U., Hussain, M., Rashid, A. et al. Numerical Solution and Analysis of Three-Dimensional Transient Darcy Flow. Transp Porous Med 123, 289–305 (2018). https://doi.org/10.1007/s11242-018-1041-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11242-018-1041-2

Keywords

Navigation