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

Abstract Adaptive time stepping with embedded error control is applied to the mixed form of Richards equation. It is the first mathematically based adaptive scheme applied to this form of Richards equation. The key to the method is the approximation of the local truncation error of the scheme in terms of the pressure head, although, to enforce mass conservation, the principal time approximation is based on the moisture content. The time stepping scheme is closely related to an implicit Thomas–Gladwell approximation and is unconditionally stable and second-order accurate. Numerical trials demonstrate that the new algorithm fully automates stepsize selection and robustly constrains temporal discretisation errors given a user tolerance. The adaptive mechanism is shown to improve the performance of the non-linear solver, providing accurate initial solution estimates for the iterative process. Furthermore, the stepsize variation patterns reflect the adequacy of the spatial discretisation, here accomplished by linear finite elements. When sufficiently dense spatial grids are used, the time step varies smoothly, while excessively coarse grids induce stepsize oscillations.

[1]  W. Press,et al.  Numerical Recipes: The Art of Scientific Computing , 1987 .

[2]  François Lehmann,et al.  Comparison of Iterative Methods for Improved Solutions of the Fluid Flow Equation in Partially Saturated Porous Media , 1998 .

[3]  S. Sloan,et al.  BIOT CONSOLIDATION ANALYSIS WITH AUTOMATIC TIME STEPPING AND ERROR CONTROL PART 1: THEORY AND IMPLEMENTATION , 1999 .

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

[5]  Emil O. Frind,et al.  Comparative error analysis in finite element formulations of the advection-dispersion equation , 1985 .

[6]  Chin-Fu Tsang,et al.  A perspective on the numerical solution of convection‐dominated transport problems: A price to pay for the easy way out , 1992 .

[7]  Philip John Binning,et al.  Practical implementation of the fractional flow approach to multi-phase flow simulation , 1999 .

[8]  Eric F. Wood,et al.  NUMERICAL EVALUATION OF ITERATIVE AND NONITERATIVE METHODS FOR THE SOLUTION OF THE NONLINEAR RICHARDS EQUATION , 1991 .

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

[10]  Frank Uhlig,et al.  Numerical Algorithms with Fortran , 1996 .

[11]  Cass T. Miller,et al.  Accurate and economical solution of the pressure-head form of Richards' equation by the method of lines , 1997 .

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

[13]  S. Sloan,et al.  Adaptive backward Euler time stepping with truncation error control for numerical modelling of unsaturated fluid flow , 2002 .

[14]  Cass T. Miller,et al.  Robust solution of Richards' equation for nonuniform porous media , 1998 .

[15]  Luca Bergamaschi,et al.  MIXED FINITE ELEMENTS AND NEWTON-TYPE LINEARIZATIONS FOR THE SOLUTION OF RICHARDS' EQUATION , 1999 .

[16]  I. Gladwell,et al.  Variable-order variable-step algorithms for second-order systems. Part 2: The codes , 1988 .

[17]  L. Shampine,et al.  Numerical Solution of Ordinary Differential Equations. , 1995 .

[18]  G. Smith,et al.  Numerical Solution of Partial Differential Equations: Finite Difference Methods , 1978 .

[19]  W. L. Wood.  Practical Time-Stepping Schemes , 1990 .

[20]  Shen-Haw Ju,et al.  Mass types, element orders and solution schemes for the Richards equation , 1997 .

[21]  G. Pinder,et al.  Computational Methods in Subsurface Flow , 1983 .

[22]  P. G. Thomsen,et al.  Local error control inSDIRK-methods , 1986 .