A mass-conservative control volume-finite element method for solving Richards’ equation in heterogeneous porous media

Cumming, Benjamin, Moroney, Timothy J., & Turner, Ian (2011) A mass-conservative control volume-finite element method for solving Richards’ equation in heterogeneous porous media. BIT Numerical Mathematics.

View at publisher


We present a mass-conservative vertex-centred finite volume method for efficiently solving the mixed form of Richards’ equation in heterogeneous porous media. The spatial discretisation is particularly well-suited to heterogeneous media because it produces consistent flux approximations at quadrature points where material properties are continuous. Combined with the method of lines, the spatial discretisation gives a set of differential algebraic equations amenable to solution using higher-order implicit solvers. We investigate the solution of the mixed form using a Jacobian-free inexact Newton solver, which requires the solution of an extra variable for each node in the mesh compared to the pressure-head form. By exploiting the structure of the Jacobian for the mixed form, the size of the preconditioner is reduced to that for the pressure-head form, and there is minimal computational overhead for solving the mixed form.

The proposed formulation is tested on two challenging test problems. The solutions from the new formulation offer conservation of mass at least one order of magnitude more accurate than a pressure head formulation, and the higher-order temporal integration significantly improves both the mass balance and computational efficiency of the solution.

Impact and interest:

4 citations in Scopus
3 citations in Web of Science®
Search Google Scholar™

Citation counts are sourced monthly from Scopus and Web of Science® citation databases.

These databases contain citations from different subsets of available publications and different time periods and thus the citation count from each is usually different. Some works are not in either database and no count is displayed. Scopus includes citations from articles published in 1996 onwards, and Web of Science® generally from 1980 onwards.

Citations counts from the Google Scholar™ indexing service can be viewed at the linked Google Scholar™ search.

ID Code: 43401
Item Type: Journal Article
Refereed: Yes
Additional Information: Springer First Online Publication.
Additional URLs:
Keywords: Finite volume, Richards' equation, DAE
DOI: 10.1007/s10543-011-0335-3
ISSN: 1572-9125
Subjects: Australian and New Zealand Standard Research Classification > MATHEMATICAL SCIENCES (010000) > NUMERICAL AND COMPUTATIONAL MATHEMATICS (010300) > Numerical Solution of Differential and Integral Equations (010302)
Divisions: Past > QUT Faculties & Divisions > Faculty of Science and Technology
Past > Schools > Mathematical Sciences
Copyright Owner: Copyright 2011 Springer Science + Business Media B.V.
Copyright Statement: The original publication is available at SpringerLink
Deposited On: 20 Jul 2011 22:09
Last Modified: 18 Mar 2016 06:23

Export: EndNote | Dublin Core | BibTeX

Repository Staff Only: item control page