Analysis of Expanded Mixed Finite Element Methods for a Nonlinear Parabolic Equation Modeling Flow into Variably Saturated Porous Media

We present an analysis of expanded mixed finite element methods applied to Richards' equation, a nonlinear parabolic partial differential equation modeling the flow of water into a variably saturated porous medium. We consider the full range of saturated to completely unsaturated media. In the case of the lowest order Raviart--Thomas spaces and the range of all possible saturations, we bound the H-1-norm of the error in capacity in terms of approximation error. This estimate uses a time-integrated scheme and the Kirchhoff transformation to handle a degeneracy in the case of completely unsaturated flow. Optimal convergence is then shown for a nonlinear form related to the error in the capacity for the case of saturated to partially saturated flow. Convergence rates depending on the Holder continuity of the capacity term are derived. Last, optimal convergence of pressures and fluxes is stated for the case of strictly partially saturated flow.

[1]  M. Wheeler,et al.  Mixed Finite Elements for Elliptic Problems with Tensor Coefficients as Cell-Centered Finite Differences , 1997 .

[2]  Todd Arbogast,et al.  Numerical methods for the simulation of flow in root-soil systems , 1993 .

[3]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[4]  M. Vauclin,et al.  A Comparative Study of Three Forms of the Richard Equation Used for Predicting One-Dimensional Infiltration in Unsaturated Soil1 , 1981 .

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

[6]  Keith L. Bristow,et al.  Simulating Water Movement in Layered and Gradational Soils Using the Kirchhoff Transform , 1990 .

[7]  Todd Arbogast,et al.  A Nonlinear Mixed Finite Eelement Method for a Degenerate Parabolic Equation Arising in Flow in Porous Media , 1996 .

[8]  Todd Arbogast.  An Error Analysis for Galerkin Approximations to an Equation of Mixed Elliptic-Parabolic Type , 1990 .

[9]  Mary F. Wheeler,et al.  A Two-Grid Finite Difference Scheme for Nonlinear Parabolic Equations , 1998 .

[10]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

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

[12]  A. Weiser,et al.  On convergence of block-centered finite differences for elliptic-problems , 1988 .

[13]  L. D. Marini,et al.  Two families of mixed finite elements for second order elliptic problems , 1985 .

[14]  M. E. Rose,et al.  Numerical methods for flows through porous media. I , 1983 .

[15]  Stephan Luckhaus,et al.  Quasilinear elliptic-parabolic differential equations , 1983 .

[16]  Ricardo H. Nochetto,et al.  Approximation of Degenerate Parabolic Problems Using Numerical Integration , 1988 .