Second-order accurate monotone finite volume scheme for Richards' equation

In this work we perform a theoretical and numerical analysis of Richards' equation. For certain types of nonlinearities we provide explicit analytical solutions. These solutions are used to show that conventional unconditionally monotone finite volume schemes have only first-order accuracy. We derive necessary and sufficient conditions for the monotonicity of finite volume discretizations and use these conditions to construct a monotone finite volume discretization accurate to second-order.

[1]  Rajesh Srivastava,et al.  Analytical solutions for one-dimensional, transient infiltration towards the water table in homogeneous and layered soils , 1991 .

[2]  D. A. Barry,et al.  A class of exact solutions for Richards' equation , 1993 .

[3]  S. Spekreijse.  Multigrid solution of monotone second-order discretizations of hyperbolic conservation laws , 1986 .

[4]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[5]  Bram van Leer,et al.  Upwind and High-Resolution Methods for Compressible Flow: From Donor Cell to Residual-Distribution Schemes , 2003 .

[6]  Heiko Berninger,et al.  Domain decomposition methods for elliptic problems with jumping nonlinearities and application to the Richards equation , 2009 .

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

[8]  Y. Stepanyants,et al.  Front Solutions of Richards’ Equation , 2008 .

[9]  Philip Broadbridge,et al.  Constant rate rainfall infiltration: A versatile nonlinear model, I. Analytic solution , 1988 .

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

[11]  Robert Eymard,et al.  The finite volume method for Richards equation , 1999 .

[12]  Dmitri Kuzmin,et al.  A vertex-based hierarchical slope limiter for p-adaptive discontinuous Galerkin methods , 2010, J. Comput. Appl. Math..

[13]  Felix Otto,et al.  L1-Contraction and Uniqueness for Quasilinear Elliptic–Parabolic Equations , 1996 .

[14]  A. W. Warrick,et al.  Numerical approximations of darcian flow through unsaturated soil , 1991 .

[15]  T.-C. Jim Yeh,et al.  One‐dimensional steady state infiltration in heterogeneous soils , 1989 .

[16]  Y. Cohen,et al.  Estimation of internodal permeabilities for numerical simulation of unsaturated flows , 2002 .

[17]  A. Pullan.  The quasilinear approximation for unsaturated porous media flow. , 1990 .

[18]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

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

[20]  Hans-Görg Roos.  A second order monotone upwind scheme , 2005, Computing.

[21]  T.-C. Jim Yeh,et al.  One-dimensional, steady vertical flow in a layered soil profile , 1990 .

[22]  R. Wallach.  RUNOFF CONTAMINATION BY SOIL CHEMICALS : TIME SCALE APPROACH , 1991 .

[23]  M. Vanclooster,et al.  A Set of Analytical Benchmarks to Test Numerical Models of Flow and Transport in Soils , 2005 .

[24]  M. Sepúlveda,et al.  On an upwind difference scheme for strongly degenerate parabolic equations modelling the settling of suspensions in centrifuges and non-cylindrical vessels , 2006 .

[25]  David Russo,et al.  Estimation of finite difference interblock conductivities for simulation of infiltration into initially dry soils , 1992 .

[26]  Gianmarco Manzini,et al.  Mass-conservative finite volume methods on 2-D unstructured grids for the Richards’ equation , 2004 .

[27]  Piotr K. Smolarkiewicz,et al.  Multidimensional positive definite advection transport algorithm: an overview , 2006 .