Alteration in molecular structure of alkali activated slag with various water to binder ratios under accelerated carbonation.

Carbonation of alkali activated materials is one of the main deteriorations affecting their durability. However, current understanding of the structural alteration of these materials exposed to an environment inducing carbonation at the nano/micro scale remains limited. This study examined the evolution of phase assemblages of alkali activated slag mortars subjected to accelerated carbonation (1% CO2, 60% relative humidity, up to 28 day carbonation) using XRD, FTIR and 29Si, 27Al, and 23Na MAS NMR. Samples with three water to binder (w/b) ratios (0.35, 0.45, and 0.55) were investigated. The results show that the phase assemblages mainly consisted of C-A-S-H, a disordered remnant aluminosilicate binder, and a minor hydrotalcite as a secondary product. Upon carbonation, calcium carbonate is mainly formed as the vaterite polymorph, while no sodium carbonate is found after carbonation as commonly reported. Sodium acts primarily as a charge balancing ion without producing sodium carbonate as a final carbonation product in the 28-day carbonated materials. The C-A-S-H structure becomes more cross-linked due to the decalcification of this phase as evidenced by the appearance of Q4 groups, which replace the Q1 and Q2 groups as observed in the 29Si MAS NMR spectra, and the dominance of Al(IV) in 27Al MAS NMR. Especially, unlike cementitious materials, the influence of w/b ratio on the crystalline phase formation and structure of C-A-S-H in the alkali activated mortars before and after carbonation is limited.

Coupling Flow, Heat, and Reactive Transport Modeling to Reproduce In Situ Redox Potential Evolution: Application to an Infiltration Pond.

Redox potential (Eh) measurements are widely used as indicators of the dominant reduction-oxidation reactions occurring underground. Yet, Eh data are mostly used in qualitative terms, as actual values cannot be used to distinguish uniquely the dominant redox processes at a sampling point and should therefore be combined with a detailed geochemical characterization of water samples. In this work, we have intensively characterized the redox potential of the first meter of soil in an infiltration pond recharged with river water using a set of in situ sensors measuring every 12 min during a 1 year period. This large amount of data combined with hydrogeochemical campaigns allowed developing a reactive transport model capable of reproducing the redox potential in space and time together with the site hydrochemistry. Our results showed that redox processes were mainly driven by the amount of sedimentary organic matter in the system as well as by seasonal variation of temperature. As a subsidiary result, our work emphasizes the need to use a fully coupled model of flow, heat transport, solute transport, and the geochemical reaction network to fully reproduce the Eh observations in the topsoil.

Nested multiresolution hierarchical simulated annealing algorithm for porous media reconstruction.

Microstructure strongly influences flow and transport properties of porous media. Flow and transport simulations within porous media, therefore, requires accurate three-dimensional (3D) models of the pore and solid phase structure. To date, no imaging method can resolve all relevant heterogeneities from the nano- to the centimeter scale within complex heterogeneous materials such as clay, reservoir rocks (e.g., travertine, chalk, ...), hardened cement paste, and concrete. To reconstruct these porous materials it is thus necessary to merge information from different 2D and potentially 3D imaging methods. One porous media reconstruction methodology that has been around for at least two decades is simulated annealing (SA). However, realizations with SA typically suffer an artificially reduced long-range connectivity, while multiphase reconstructions are not feasible in most cases because of a prohibitive computational burden. To solve these problems we propose a hierarchical multiresolution and multiphase simulated annealing algorithm. To decrease the computational cost of multiphase simulation, our algorithm sequentially simulates one phase after another, in a hierarchical way, which enables handling multimodal distributions and topological relations. Building upon recent work, our algorithm improves long-range connectivity and CPU efficiency by simulating larger particles using a coarser resolution that is subsequently refined compared to standard SA; our proposed extension not only offers the possibility to perform multiphase reconstruction but also allows us (i) to improve binary reconstruction quality, as quantified, e.g., by multiple-point histograms by up to one order of magnitude and (ii) to achieve an overall speed-up. The proposed algorithm is also shown to outperform the direct sampling multiple-point statistics method for the generation of cement paste microstructure with respect to both generation time and quality.

Numerically accelerated pore-scale equilibrium dissolution.

Simulation of dissolution processes with a pore-scale reactive transport model increases insight in coupled chemical-physical-transport processes. However, modelling of dissolution process often requires a large number of time steps especially when the buffering capacity of solid phases is high. In this work we analyze the interplay between solid buffering on one hand and transport on the other. Based on this analysis we propose an approach to reduce the number of required time steps for simulating equilibrium dissolution processes. The underlying idea is that the number of time step iterations can be reduced if the buffering is sufficient to bring the system to a steady state, i.e. that the concentration field around solid is time-invariant. If this condition is satisfied, then it is possible to reduce the physical (and thus also computational) time by adjusting the chemical system appropriately. First we derived a dimensionless value - called buffering number - to determine under which conditions reduction in time can be made. Several examples illustrate that below a certain buffering number, the physical time can be reduced without significant effect on result (e.g. dissolution front) as long as the solid volume fraction is sufficient. This means that for a given solid-liquid system, the calculation time can be reduced either by the reduction of mass in solid or by the increase of equilibrium concentration (solubility). We also show that the calculation time for calcium leaching in cementitious systems can be reduced by 50 times with a negligible error.

Experimental and Numerical Investigations of Silver Nanoparticle Transport under Variable Flow and Ionic Strength in Soil.

Unsaturated column experiments were conducted with an undisturbed loamy sand soil to investigate the influence of flow interruption (FI) and ionic strength (IS) on the transport and retention of surfactant-stabilized silver nanoparticles (AgNP) and the results were compared to those obtained under continuous flow conditions. AgNP concentrations for breakthrough curves (BTCs) and retention profiles (RPs) were analyzed by ICP-MS. Experimental results were simulated by the numerical code HP1 (Hydrus-PhreeqC) with the DLVO theory, extended colloid filtration theory and colloid release model. BTCs of AgNP showed a dramatic drop after FI compared to continuous flow conditions. Evaporation increased due to FI, resulting in increased electrical conductivity of the soil solution, which led to a totally reduced mobility of AgNP. A reduction of IS after FI enhanced AgNP mobility slightly. Here the strongly increased Al and Fe concentration in the effluent suggested that soil colloids facilitated the release of AgNP (cotransport). The numerical model reproduced the measured AgNP BTCs and indicated that attachment to the air-water interface (AWI) occurring during FI was the key process for AgNP retention.

Impact of manure-related DOM on sulfonamide transport in arable soils.

Field application of livestock manure introduces colloids and veterinary antibiotics, e.g. sulfonamides (SAs), into farmland. The presence of manure colloids may potentially intensify the SAs-pollution to soils and groundwater by colloid-facilitated transport. Transport of three SAs, sulfadiazine (SDZ), sulfamethoxypyridazine (SMPD), and sulfamoxole (SMOX), was investigated in saturated soil columns with and without manure colloids from sows and farrows, weaners, and fattening pigs. Experimental results showed that colloid-facilitated transport of SMOX was significant in the presence of manure colloids from fattening pigs with low C/N ratio, high SUVA280nm and protein C, while manure colloids from sows and farrows and weaners had little effect on SMOX transport. In contrast, only retardation was observed for SDZ and SMPD when manure colloids were present. Breakthrough curves (BTCs) of colloids and SAs were replicated well by a newly developed numerical model that considers colloid-filtration theory, competitive kinetic sorption, and co-transport processes. Model results demonstrate that mobile colloids act as carriers for SMOX, while immobile colloids block SMOX from sorbing onto the soil. The low affinity of SMOX to sorb on immobile colloids prevents aggregation and also promotes SMOX's colloid-facilitated transport. Conversely, the high affinity of SDZ and SMPD to sorb on all types of immobile colloids retarded their transport. Thus, manure properties play a fundamental role in increasing the leaching risk of hydrophobic sulfonamides.

A reactive transport model for mercury fate in contaminated soil--sensitivity analysis.

We present a sensitivity analysis of a reactive transport model of mercury (Hg) fate in contaminated soil systems. The one-dimensional model, presented in Leterme et al. (2014), couples water flow in variably saturated conditions with Hg physico-chemical reactions. The sensitivity of Hg leaching and volatilisation to parameter uncertainty is examined using the elementary effect method. A test case is built using a hypothetical 1-m depth sandy soil and a 50-year time series of daily precipitation and evapotranspiration. Hg anthropogenic contamination is simulated in the topsoil by separately considering three different sources: cinnabar, non-aqueous phase liquid and aqueous mercuric chloride. The model sensitivity to a set of 13 input parameters is assessed, using three different model outputs (volatilized Hg, leached Hg, Hg still present in the contaminated soil horizon). Results show that dissolved organic matter (DOM) concentration in soil solution and the binding constant to DOM thiol groups are critical parameters, as well as parameters related to Hg sorption to humic and fulvic acids in solid organic matter. Initial Hg concentration is also identified as a sensitive parameter. The sensitivity analysis also brings out non-monotonic model behaviour for certain parameters.

A reactive transport model for mercury fate in soil--application to different anthropogenic pollution sources.

Soil systems are a common receptor of anthropogenic mercury (Hg) contamination. Soils play an important role in the containment or dispersion of pollution to surface water, groundwater or the atmosphere. A one-dimensional model for simulating Hg fate and transport for variably saturated and transient flow conditions is presented. The model is developed using the HP1 code, which couples HYDRUS-1D for the water flow and solute transport to PHREEQC for geochemical reactions. The main processes included are Hg aqueous speciation and complexation, sorption to soil organic matter, dissolution of cinnabar and liquid Hg, and Hg reduction and volatilization. Processes such as atmospheric wet and dry deposition, vegetation litter fall and uptake are neglected because they are less relevant in the case of high Hg concentrations resulting from anthropogenic activities. A test case is presented, assuming a hypothetical sandy soil profile and a simulation time frame of 50 years of daily atmospheric inputs. Mercury fate and transport are simulated for three different sources of Hg (cinnabar, residual liquid mercury or aqueous mercuric chloride), as well as for combinations of these sources. Results are presented and discussed with focus on Hg volatilization to the atmosphere, Hg leaching at the bottom of the soil profile and the remaining Hg in or below the initially contaminated soil layer. In the test case, Hg volatilization was negligible because the reduction of Hg(2+) to Hg(0) was inhibited by the low concentration of dissolved Hg. Hg leaching was mainly caused by complexation of Hg(2+) with thiol groups of dissolved organic matter, because in the geochemical model used, this reaction only had a higher equilibrium constant than the sorption reactions. Immobilization of Hg in the initially polluted horizon was enhanced by Hg(2+) sorption onto humic and fulvic acids (which are more abundant than thiols). Potential benefits of the model for risk management and remediation of contaminated sites are discussed.

Determination of 137Cs contamination depth distribution in building structures using geostatistical modeling of ISOCS measurements.

Decommissioning of nuclear building structures usually leads to large amounts of low level radioactive waste. Using a reliable method to determine the contamination depth is indispensable prior to the start of decontamination works and also for minimizing the radioactive waste volume and the total workload. The method described in this paper is based on geostatistical modeling of in situ gamma-ray spectroscopy measurements using the multiple photo peak method. The method has been tested on the floor of the waste gas surge tank room within the BR3 (Belgian Reactor 3) decommissioning project and has delivered adequate results.

Inverse optimization of hydraulic, solute transport, and cation exchange parameters using HP1 and UCODE to simulate cation exchange.

Reactive transport modeling is a powerful tool to evaluate systems with complex geochemical relations. However, parameters are not always directly measurable. This study represents one of the first attempts to obtain hydrologic, transport and geochemical parameters from an experimental dataset involving transient unsaturated water flow and solute transport, using an automatic inverse optimization (or calibration) algorithm. The data come from previously published, controlled laboratory experiments on the transport of major cations (Na, K, Mg, Ca) during water absorption into horizontal soil columns that were terminated at different times. Experimental data consisted of the depth profiles of water contents (θ), Cl concentrations, and total aqueous and sorbed concentrations of major cations. The dataset was used to optimize several parameters using the reactive transport model, HP1 and the generic optimization code, UCODE. Although the soil hydraulic and solute transport parameters were also optimized, the study focused mainly on the geochemical parameters because the soil columns were constructed from disturbed soil. The cation exchange capacity and the cation exchange coefficients for two exchange models (Gapon and Rothmund-Kornfeld) were optimized. The results suggest that both calibrated models satisfactorily described the experimental data, although the Rothmund-Kornfeld model fit was slightly better. However, information content and surface response analyses indicated that parameters of the Gapon model are well identifiable, whereas those of the Rothmund-Kornfeld model were strongly correlated. The calibrated geochemical parameters were validated using an independent dataset. In agreement with the identifiability analysis, the Gapon approach was better than the Rothmund-Kornfeld model at calculating the observed concentrations of major cations in the soil solution and on the exchange sites.