An extension of the open-source porousMultiphaseFoam toolbox dedicated to groundwater flows solving the Richards' equation

In this note, the existing porousMultiphaseFoam toolbox, developed initially for any two-phase flow in porous media is extended to the specific case of the Richards' equation which neglect the pressure gradient of the non-wetting phase. This model is typically used for saturated and unsaturated groundwater flows. A Picard's algorithm is implemented to linearize and solve the Richards' equation developed in the pressure head based form. This new solver of the porousMultiphaseFoam toolbox is named groundwaterFoam. The validation of thesolver is achieved by a comparison between numerical simulations and results obtained from the literature. Finally, a parallel efficiency test is performed on a large unstructured mesh and exhibits a super-linear behavior as observed for the other solvers of the toolbox.

[1]  R. H. Brooks,et al.  Hydraulic properties of porous media , 1963 .

[2]  Cyprien Soulaine,et al.  An open-source toolbox for multiphase flow in porous media , 2015, Comput. Phys. Commun..

[3]  H. Rusche.  Computational fluid dynamics of dispersed two-phase flows at high phase fractions , 2003 .

[4]  Dmitri Kavetski,et al.  Adaptive time stepping and error control in a mass conservative numerical solution of the mixed form of Richards equation | NOVA. The University of Newcastle's Digital Repository , 2001 .

[5]  Linda M. Abriola,et al.  Mass conservative numerical solutions of the head‐based Richards equation , 1994 .

[6]  Cass T. Miller,et al.  An evaluation of temporally adaptive transformation approaches for solving Richards' equation , 1999 .

[7]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[8]  Michel Quintard,et al.  An open source massively parallel solver for Richards equation: Mechanistic modelling of water fluxes at the watershed scale , 2014, Comput. Phys. Commun..

[9]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .

[10]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[11]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[12]  Morris Muskat,et al.  Physical principles of oil production , 1949 .

[13]  J. W. Sheldon,et al.  One-Dimensional, Incompressible, Noncapillary, Two-Phase Fluid Flow in a Porous Medium , 1959 .

[14]  E. Youngs,et al.  Fundamentals of Soil Physics. , 1982 .