A Mass Conservative Method for Numerical Modeling of Axisymmetric flow

Document Type: Research Paper


Department of Civil Engineering, K.N. Toosi University of Technology, Tehran, Iran


In this paper, the axisymmetric flow toward a pumping well has been numerically solved by the cell-centered finite volume method. The numerical model is descretized over unstructured and triangular-shaped grids which allows to simulate inhomogeneous and complex-shaped domains. Due to the non-orthogonality of the irregular grids, the multipoint flux approximation (MPFA) schemes are used to discretize the flux term. In this work, the diamond scheme as the MPFA method has been employed and the least square method is applied to express the full discrete form of the vertex-values of the hydraulic head. The scheme has been verified via the Theis solution, as a milestone in well hydraulics. The numerical results show the capability of the developed model in evaluating transient drawdown in the confined aquifers. The proposed numerical model leads to the  stable and local conservative solutions contrary to the standard finite element methods. Also this numerical technique has the second order of accuracy.


  1. Theis CV (1935) The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using ground-water storage. Transactions American Geophysical Union, 16, 519–24.
  2. Hantush, MS, Jacob CE (1955) Non-steady radial flow in an infinite leaky aquifer. Eos, Transactions American Geophysical Union, 36(1), 95-100.
  3. Neuman SP (1972) Theory of flow in unconfined aquifers considering delayed gravity response. Water Resources Research 8(4), 1031–1045.
  4. Neuman SP (1973) Supplementary comments on theory of flow in unconfined aquifers considering delayed gravity response. Water Resources Research 9(4), 1102–1103.
  5. Neuman SP (1974) Effect of partial penetration on flow in unconfined aquifers considering delayed gravity response. Water Resources Research 10(2), 303–312.
  6. Neuman SP (1975) Analysis of pumping test data from anisotropic unconfined aquifers. Water Resources Research, 11(2), 329–345.
  7. Vermeer PA, Verruijt A (1981) An accuracy condition for consolidation by finite elements. International Journal of Numeical and Analytical Methods in Geomechanics, 5(1), 1-14.
  8. Murad MA, Loula AFD (1992) Improved accuracy in finite element analysis of Biot’s consolidation problem. Computer Methods in Applied Mechanics and Engineering 95(3), 359–382.
  9. Ferronato M, Castelletto N, Gambolati G (2010) A fully coupled 3-D mixed finite element model of Biot consolidation. Journal of Computational Physics, 229(12), 4813–4830.
  10. Kim J (2010) Sequential methods for coupled geomechanics and multiphase flow. PhD thesis, Stanford University.
  11. Asadi R, Ataie-Ashtiani B (2015) A comparison of finite volume formulations and coupling strategies for two-phase flow in deforming porous media. Computers and Geotechnics, 67, 17-32.
  12. Asadi R, Ataie-Ashtiani B, Simmons CT (2014) Finite volume coupling strategies for the solution of a Biot consolidation model. Computers and Geotechnics, 55, 494-505.
  13. Asadi R, Ataie-Ashtiani B (2016) Numerical modeling of subsidence in saturated porous media: A mass conservative method. Journal of hydrology, 542, 423-436.
  14. Caviedes-Voullième D, Garc─▒ P, Murillo J (2013) Verification, conservation, stability and efficiency of a finite volume method for the 1D Richards equation. Journal of hydrology, 480, 69-84.
  15. Manzini G, Ferraris S (2004) Mass-conservative finite volume methods on 2-D unstructured grids for the Richards’ equation. Advances in Water Resources, 27(12), 1199-1215.
  16. Lee SH, Jenny P, Tchelepi HA (2002) A finite-volume method with hexahedral multiblock grids for modeling flow in porous media. Computers and Geosciences, 6(3-4), 353-379.
  17. Edwards MG (2002) Unstructured, control-volume distributed, full-tensor finite-volume schemes with flow based grids. Computational Geosciences, 6(3-4), 433-452.
  18. Szymkiewicz A (2012) Modelling water flow in unsaturated porous media: accounting for nonlinear permeability and material heterogeneity. Springer Science & Business Media.
  19. Coudière Y, Villedieu P (2000) Convergence rate of a finite volume scheme for the linear convection-diffusion equation on locally refined meshes. Mathematical Modelling and Numerical Analysis, 34(6), 1123-1149.
  20. Coudière Y, Vila JP, Villedieu P (1999) Convergence rate of a finite volume scheme for a two dimensional convection–diffusion problem, Mathematical Modelling and Numerical Analysis, 33, 493–516.
  21. Bertolazzi E, Manzini G (2004) A cell-centered second-order accurate finite volume method for convection–diffusion problems on unstructured meshes. Mathematical Models and Methods in Applied Sciences, 14(8), 1235-1260.
  22. Bertolazzi E, Manzini G (2005) A unified treatment of boundary conditions in least-square based finite-volume methods. Computers & Mathematics with Applications, 49(11), 1755-1765.
  23. Bevilacqua I, Canone D, Ferraris S (2011) Acceleration techniques for the iterative resolution of the Richards equation by the finite volume method. International Journal for Numerical Methods in Biomedical Engineering, 27(8), 1309-1320.
  24. Biot MA (1941) General theory of three-dimensional consolidation. Journal of Applied Physics, 12, 155–164.
  25. Edwards MG, Rogers C (1998) Finite volume discretization with imposed flux continuity for the general tensor pressure equation. Computers & Geosciences, 2(4), 259–290.
  26. Schöberl J (1997) NETGEN An advancing front 2D/3D-mesh generator based on abstract rules. Computing and Visualization in Science, 1(1), 41-52.
  27. Fetter CW (1994) Applied Hydrogeology, third ed. Prentice Hall, New Jersey.
  28. Bedient PB, Rifai HS, Newell CJ (1994). Ground water contamination: transport and remediation. Prentice-Hall International, Inc.