Premium
Mass conservative numerical solutions of the head‐based Richards equation
Author(s) -
Rathfelder Klaus,
Abriola Linda M.
Publication year - 1994
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/94wr01302
Subject(s) - discretization , mathematics , finite element method , chord (peer to peer) , equivalence (formal languages) , term (time) , mathematical analysis , computer science , physics , thermodynamics , distributed computing , discrete mathematics , quantum mechanics
Numerical procedures for efficient mass conservative solutions of the head‐based form of the Richards equation are presented. Mass conservative solutions are shown to result when the capacity coefficient, C , is formulated by equating the storage term and its chain rule expansion in their discretized forms. Equivalence in the storage term expansion is maintained in finite difference models when C is evaluated with a standard chord slope approximation. This scheme is shown to produce excellent global mass balance accuracy in simulations of vertical moisture infiltration. An analogous approach to the expansion of the storage term using finite elements results in element dependent expressions of C . Application of this approach produces mass balance accuracy with errors less than 1%, but also exhibits slow convergence in the consistent form. A nontraditional finite element procedure is presented which maintains equivalence in the storage term expansion when C is evaluated with the standard chord slope approximation. This scheme exhibits excellent mass balance accuracy, in either the consistent or lumped forms, without significant loss in computational efficiency.