On the use of large time steps with ELLAM for transport with kinetic reactions over heterogeneous domains

Abstract

An Eulerian Lagrangian localized adjoint method (ELLAM) is considered for the resolution of advection‐dominated transport problems in porous media. Contrary to standard Eulerian methods, ELLAM can use large time steps because the advection term is approximated accurately without any CFL restriction. However, it is shown in this article that special care must be taken for the approximation of the dispersive and reactive terms when large time steps are used over heterogeneous domains. An alternative procedure is proposed. It is based on an equivalent dispersion coefficient or an equivalent reaction rate when different zones are encountered during the tracking. Numerical experiments are performed with variable dispersion or variable reaction rates over space (including nonlinearity). When classical ELLAM require numerous time steps to handle heterogeneity, the alternative procedure is shown to perform with the same accuracy in a single time step. © 2009 American Institute of Chemical Engineers AIChE J, 2009

Estimation of macrodispersion in 2‐D highly heterogeneous porous media using the Eulerian‐Lagrangian localized adjoint method

“…The scheme has been shown to be accurate and efficient for transport in highly heterogeneous media with sharp concentration fronts [ Ramasomanana and Younes , ; Ramasomanana et al ., ]. For heterogeneous dispersion, the scheme uses an equivalent dispersion tensor , which is calculated as the arithmetic mean of local dispersion tensors weighted by the mean time spent under the influence of these local values [ Fahs et al ., ; Ramasomanana and Younes , ].…”

Effect of temperature on convective-reactive transport of CO2 in geological formations

“…In addition, in existing codes, reactive processes are usually included via the operator splitting approach. For kinetic reactions, this approach introduces intrinsic splitting errors that are proportional to the time step used in the numerical solution (Fahs et al, 2009). These numerical artifacts can be avoided by using dense computational grids with small time steps, which increase the computational requirements and the CPU time of simulations.…”

Acid/base front propagation in saturated porous media: 2D laboratory experiments and modeling

“…Several numerical techniques and codes have been developed to simulate reactive transport processes in natural groundwater systems (Fahs et al, 2009;Kinzelbach et al, 1991;Saaltink et al, 2001;Van der Lee et al, 2003;Walter et al, 1994;Yeh, and Tripathi, 1991 among others). The mechanisms involved in reactive transport processes are often strongly non-linear and occur on a multiplicity of time scales, thus rendering the solution of a mathematical model very challenging.…”