Peter Alt-Epping

Benchmark reactive transport simulations of a column experiment in compacted bentonite with multispecies diffusion and explicit treatment of electrostatic effects

Bentonite clay is considered as a potential buffer and backfill material in subsurface repositories for high-level nuclear waste. As a result of its low permeability, transport of water and solutes in compacted bentonite is driven primarily by diffusion. Developing models for species transport in bentonite is complicated, because of the interaction of charged species and the negative surface charge of clay mineral surfaces. The effective diffusion coefficient of an ion in bentonite depends on the ion’s polarity and valence, on the ionic strength of the solution, and on the bulk dry density of the bentonite. These dependencies need to be understood and incorporated into models if one wants to predict the effectiveness of bentonite as a barrier to radionuclides in a nuclear repository. In this work, we present a benchmark problem for reactive transport simulators based on a flow-through experiment carried out on a saturated bentonite core. The measured effluent composition shows the complex interplay of species transport in a charged medium in combination with sorption and mineral precipitation/dissolution reactions. The codes compared in this study are PHREEQC, CrunchFlow, FLOTRAN, and MIN3P. The benchmark problem is divided into four component problems of increasing complexity, leading up to the main problem which addresses the effects of advective and diffusive transport of ions through bentonite with explicit treatment of electrostatic effects. All codes show excellent agreement between results provided that the activity model, Debye-Hückel parameters, and thermodynamic data used in the simulations are consistent. A comparison of results using species-specific diffusion and uniform species diffusion reveals that simulated species concentrations in the effluent differ by less than 8 %, and that these differences vanish as the system approaches steady state.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

Benchmarking the simulation of Cr isotope fractionation

A benchmark problem set consisting of four problem levels was developed for the simulation of Cr isotope fractionation in 1D and 2D domains. The benchmark is based on a recent field study where Cr(VI) reduction and accompanying Cr isotope fractionation occurs abiotically by an aqueous reaction with dissolved Fe 2+ (Wanner et al., 2012., Appl. Geochem., 27, 644–662). The problem set includes simulation of the major processes affecting the Cr isotopic composition such as the dissolution of various Cr(VI) bearing minerals, fractionation during abiotic aqueous Cr(VI) reduction, and non-fractionating precipitation of Cr(III) as sparingly soluble Cr-hydroxide.Accuracy of the presented solutions was ensured by running the problems with four well-established reactive transport modeling codes: TOUGHREACT, MIN3P, CRUNCHFLOW, and FLOTRAN. Results were also compared with an analytical Rayleigh-type fractionation model. An additional constraint on the correctness of the results was obtained by comparing output from the problem levels simulating Cr isotope fractionation with the corresponding ones only simulating bulk concentrations. For all problem levels, model to model comparisons showed excellent agreement, suggesting that for the tested geochemical processes any code is capable of accurately simulating the fate of individual Cr isotopes.

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions.

Transient convective circulation in a sedimented ridge setting: Implications for applying geochemical data to estimate permeability

[1] For an open, sedimented ridge setting we investigate by means of transient numerical simulations hydrologic factors that can influence chemical signatures of mineral alteration, and how these features relate to the permeability of the upper oceanic crust. The transient state we consider involves three stages: (1) starting from a static flow field and a conductive thermal regime, flow in basement rocks is initiated as unstable convection; (2) with the inflow of seawater to the basement, the convection pattern evolves into a stable convection cell linking the recharge and discharge sites on the seafloor; and (3) convection weakens through time as the system offsets the initial thermal condition to attain a steady state. The duration of each stage, and the extent of the region in the basement undergoing fluid exchange with seawater varies with permeability. By the end of stage 2, for basement permeabilities between 1 × 10 -14 and 1 × 10 -13 m 2 , most of the basement and the deeper sediments surrounding the discharge zone is hydrologically isolated, but with active recirculation inside this zone. Through-flow of seawater is restricted to a relatively narrow boundary layer. The unstable flow geometry during stages one and two is reflected in fluctuations of recharge and discharge rates into and out of the basement, and periodic variations in temperature of the discharging fluids. Fluid discharge at steady state is approximately one-half of the average occurring during stage one. Changing flow patterns through time lead to a strongly heterogeneous evolution of water/rock ratios in the basement. Average water/ rock ratios for the basement approach a log-log linear correlation with basement permeability as steady state is reached.