Coupled large strain consolidation and solute transport. I: Model development

The development of a numerical model, CST1, for coupled large strain consolidation and solute transport in saturated porous media is presented. The consolidation algorithm is one-dimensional and includes the capabilities of a previous code, CS2, with the addition of time-dependent loading, unload/reload effects, and an externally-applied hydraulic gradient. The solute transport algorithm is two-dimensional and accounts for advection, longitudinal and transverse dispersion, first-order decay reactions, and linear equilibrium sorption. Solute transport is consistent with temporal and spatial variations of porosity and seepage velocity in the consolidating layer. The key to the transport model is the definition of two Lagrangian fields of elements that follow the motions of fluid and solid phases separately. This reduces numerical dispersion and simplifies transport calculations to that of dispersion mass flow between contiguous fluid elements. The effect of relative numerical resolution of fluid and solid elements on the accuracy of sorption/desorption is also discussed. This paper presents the theoretical and numerical development of the CST1 model. A companion paper presents verification checks of CSTI and the results of simulations that illustrate the significance of consolidation-induced solute transport for some interesting numeric examples.