Empirical versus time stepping with embedded error control for density‐driven flow in porous media

[1] Modeling density-driven flow in porous media may require very long computational time due to the nonlinear coupling between flow and transport equations. Time stepping schemes are often used to adapt the time step size in order to reduce the computational cost of the simulation. In this work, the empirical time stepping scheme which adapts the time step size according to the performance of the iterative nonlinear solver is compared to an adaptive time stepping scheme where the time step length is controlled by the temporal truncation error. Results of the simulations of the Elder problem show that (1) the empirical time stepping scheme can lead to inaccurate results even with a small convergence criterion, (2) accurate results are obtained when the time step size selection is based on the truncation error control, (3) a non iterative scheme with proper time step management can be faster and leads to more accurate solution than the standard iterative procedure with the empirical time stepping and (4) the temporal truncation error can have a significant effect on the results and can be considered as one of the reasons for the differences observed in the Elder numerical results.

[1]  Annamaria Mazzia,et al.  High order Godunov mixed methods on tetrahedral meshes for density driven flow simulations in porous media , 2005 .

[2]  O. Kolditz,et al.  Coupled groundwater flow and transport : 1. Verification of variable density flow and transport models , 1998 .

[3]  M. Buès,et al.  Coupled Groundwater Flow and Transport in Porous Media. A Conservative or Non-conservative Form? , 2001 .

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

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

[6]  G. Wittum,et al.  The saltpool benchmark problem – numerical simulation of saltwater upconing in a porous medium , 2002 .

[7]  Jesús Carrera,et al.  Mass balance errors when solving the convective form of the transport equation in transient flow problems , 2004 .

[8]  Jesús Carrera,et al.  Role of salt sources in density‐dependent flow , 2009 .

[9]  A new coupling algorithm for density‐driven flow in porous media , 2004 .

[10]  I. Aavatsmark,et al.  An Introduction to Multipoint Flux Approximations for Quadrilateral Grids , 2002 .

[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]  Luca Bergamaschi,et al.  On the Reliability of Numerical Solutions of Brine Transport in Groundwater: Analysis of Infiltration from a Salt Lake , 2001 .

[13]  Dmitri Kavetski,et al.  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 , 2001 .

[14]  P. Ackerer,et al.  Solving the advection–dispersion equation with discontinuous Galerkin and multipoint flux approximation methods on unstructured meshes , 2008 .

[15]  S. Majid Hassanizadeh,et al.  On the modeling of brine transport in porous media , 1988 .

[16]  S. Turek.  a Comparative Study of Time-Stepping Techniques for the Incompressible Navier-Stokes Equations: from Fully Implicit Non-Linear Schemes to Semi-Implicit Projection Methods , 1996 .

[17]  Michel C. Boufadel,et al.  Numerical modeling of water flow below dry salt lakes: effect of capillarity and viscosity , 1999 .

[18]  Olaf Kolditz,et al.  Variable-density flow and transport in porous media: approaches and challenges , 2002 .

[19]  P. F. Andersen,et al.  Saltwater intrusion in aquifers: Development and testing of a three‐dimensional finite element model , 1987 .

[20]  Upwelling and downwelling behavior in the Elder‐Voss‐Souza benchmark , 2007 .

[21]  C. Oldenburg,et al.  Dispersive Transport Dynamics in a Strongly Coupled Groundwater‐Brine Flow System , 1995 .

[22]  Using quantitative indicators to evaluate results from variable-density groundwater flow models , 2005 .

[23]  J. W. Elder.  Numerical experiments with free convection in a vertical slot , 1966, Journal of Fluid Mechanics.

[24]  P. Ackerer,et al.  Modeling Variable Density Flow and Solute Transport in Porous Medium: 1. Numerical Model and Verification , 1999 .

[25]  Olaf Kolditz,et al.  Coupled groundwater flow and transport: 2. Thermohaline and 3D convection systems , 1998 .

[26]  M. Putti,et al.  Mixed-finite element and finite volume discretization for heavy brine simulations in groundwater , 2002 .

[27]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[28]  Philippe Ackerer,et al.  Efficient approximations for the simulation of density driven flow in porous media , 2008 .

[29]  P. Frolkovic,et al.  Numerical modelling of convection dominated transport coupled with density driven flow in porous media , 2000 .

[30]  Niko E. C. Verhoest,et al.  Assessment of adaptive and heuristic time stepping for variably saturated flow , 2007 .

[31]  H.-J. Diersch,et al.  Finite element modelling of recirculating density-driven saltwater intrusion processes in groundwater , 1988 .

[32]  C. Paniconi,et al.  Picard and Newton linearization for the coupled model for saltwater intrusion in aquifers , 1995 .

[33]  W. R. Souza,et al.  Variable density flow and solute transport simulation of regional aquifers containing a narrow freshwater‐saltwater transition zone , 1987 .

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