Three-dimensional density-dependent flow and multicomponent reactive transport modeling of chlorinated solvent oxidation by potassium permanganate

A popular method for the treatment of aquifers contaminated with chlorinated solvents is chemical oxidation based on the injection of potassium permanganate (KMnO₄). Both the high density (1025 gL^− 1) and reactivity of the treatment solution influence the fate of permanganate (MnO₄) in the subsurface and affect the degree of contaminant treatment. The MIN3P multicomponent reactive transport code was enhanced to simulate permanganate-based remediation, to evaluate the pathways of MnO₄ utilization, and to assess the role of density contrasts for the delivery of the treatment solution. The modified code (MIN3P-D) provides a direct coupling between density-dependent fluid flow, solute transport, contaminant treatment, and geochemical reactions. The model is used to simulate a field trial of TCE oxidation in a sandy aquifer that is underlain by an aquitard. Three-dimensional simulations are conducted for a coupled reactive system comprised of ten aqueous components, two mineral phases, TCE (dissolved, adsorbed, and NAPL), reactive organic matter, and including ion exchange reactions. Model parameters are constrained by literature data and a detailed data set from the field site under investigation. The general spatial and transient evolution in observed concentrations of the oxidant, dissolved TCE, and reaction products are adequately reproduced by the simulations. The model elucidates the important role of density-induced flow and transport on the distribution of the treatment solution into NAPL containing regions located at the aquifer–aquitard interface. Model results further suggest that reactions that do not directly affect the stability of MnO₄ have a negligible effect on solution density and MnO₄ delivery.     

A review on reactive transport model and porosity evolution in the porous media

Environmental Science and Pollution Research - This work comprehensively reviews the equations governing multicomponent flow and reactive transport in porous media on the pore-scale, mesoscale and...     

Variably saturated flow and multicomponent biogeochemical reactive transport modeling of a uranium bioremediation field experiment

Three-dimensional, coupled variably saturated flow and biogeochemical reactive transport modeling of a 2008 in situ uranium bioremediation field experiment is used to better understand the interplay of transport and biogeochemical reactions controlling uranium behavior under pulsed acetate amendment, seasonal water table variation, spatially variable physical (hydraulic conductivity, porosity) and geochemical (reactive surface area) material properties. While the simulation of the 2008 Big Rusty acetate biostimulation field experiment in Rifle, Colorado was generally consistent with behaviors identified in previous field experiments at the Rifle IFRC site, the additional process and property detail provided several new insights. A principal conclusion from this work is that uranium bioreduction is most effective when acetate, in excess of the sulfate-reducing bacteria demand, is available to the metal-reducing bacteria. The inclusion of an initially small population of slow growing sulfate-reducing bacteria identified in proteomic analyses led to an additional source of Fe(II) from the dissolution of Fe(III) minerals promoted by biogenic sulfide. The falling water table during the experiment significantly reduced the saturated thickness of the aquifer and resulted in reactants and products, as well as unmitigated uranium, in the newly unsaturated vadose zone. High permeability sandy gravel structures resulted in locally high flow rates in the vicinity of injection wells that increased acetate dilution. In downgradient locations, these structures created preferential flow paths for acetate delivery that enhanced local zones of TEAP reactivity and subsidiary reactions. Conversely, smaller transport rates associated with the lower permeability lithofacies (e.g., fine) and vadose zone were shown to limit acetate access and reaction. Once accessed by acetate, however, these same zones limited subsequent acetate dilution and provided longer residence times that resulted in higher concentrations of TEAP reaction products when terminal electron donors and acceptors were not limiting. Finally, facies-based porosity and reactive surface area variations were shown to affect aqueous uranium concentration distributions with localized effects of the fine lithofacies having the largest impact on U(VI) surface complexation. The ability to model the comprehensive biogeochemical reaction network, and spatially and temporally variable processes, properties, and conditions controlling uranium behavior during engineered bioremediation in the naturally complex Rifle IFRC subsurface system required a subsurface simulator that could use the large memory and computational performance of a massively parallel computer. In this case, the eSTOMP simulator, operating on 128 processor cores for 12h, was used to simulate the 110-day field experiment and 50 days of post-biostimulation behavior.     

A triple-continuum approach for modeling flow and transport processes in fractured rock

This paper presents a triple-continuum conceptual model for simulating flow and transport processes in fractured rock. Field data collected from the unsaturated zone of Yucca Mountain, a repository site of high-level nuclear waste, show a large number of small-scale fractures. The effect of these small fractures has not been considered in previous modeling investigations within the context of a continuum approach. A new triple-continuum model (consisting of matrix, small-fracture, and large-fracture continua) has been developed to investigate the effect of these small fractures. This paper derives the model formulation and discusses the basic triple-continuum behavior of flow and transport processes under different conditions, using both analytical solutions and numerical approaches. The simulation results from the site-scale model of the unsaturated zone of Yucca Mountain indicate that these small fractures may have an important effect on radionuclide transport within the mountain.     

Multiphase flow and transport through fractured heterogeneous porous media

The migration of Dense, Non-Aqueous Phase Liquid (DNAPL) and dissolved phase contamination through a fractured heterogeneous porous medium has been investigated through the use of a multiphase compositional model. The sensitivity of the timescales of migration and the distribution of contaminant in the subsurface to the mean permeability, the variance of the permeability, and the degree of fracturing of the domain were examined. It was found that increasing the mean permeability of the domain allowed the DNAPL to penetrate deeper into the subsurface, while decreasing the mean permeability caused the DNAPL to pool at shallower depths. The presence of fractures within the system was found to control the infiltration only in the most fractured domain. Moment analysis of the nonwetting phase showed that large-scale movement had ceased after approximately 9 years (maximum duration of the source-on condition was approximately 4.5 years). This tended to be due to a redistribution of the DNAPL towards a residual configuration, as was evidenced by the gradual trending of average nonwetting phase saturations within the domain to a static value. The dissolved phase plume was found to migrate at essentially the same rate as the nonwetting phase, due to the reduced relative permeability of lenses containing DNAPL, and due to diffusive losses of mass to the matrix of fractured clay and silty-clay lenses. Some exceptions to this were found when the DNAPL could not overcome the displacement pressure of a lens, and could not by-pass the lens due to the lack of available driving force after the source had been shut off.     

Process-based reactive transport modeling of a permeable reactive barrier for the treatment of mine drainage

Reactive transport modeling of a permeable reactive barrier for the treatment of mine drainage was used to integrate a comprehensive data set including pore water chemistry and solid phase data from several sampling events over a >3-year time period. The simulations consider the reduction of sulfate by the organic carbon-based treatment material and the removal of sulfate and iron by precipitation of reduced mineral phases including iron monosulfides and siderite. Additional parameters constraining the model include dissolved H2S, alkalinity and pH, as well as a suite of solid phase S-fractions identified by extractions. Influences of spatial heterogeneity necessitated the use of a 2-dimensional modeling approach. Simulating observed seasonal fluctuations and long-term changes in barrier reactivity required the use of temperature dependent rate coefficients and a multimodal Monod-type rate expression accounting for the variable reactivity of different organic carbon fractions. Simulated dissolved concentrations of SO4, Fe, H2S, alkalinity and pH, as well as solid phase accumulations of reduced sulfur phases generally compare well to observed trends over 23 months. Spatial variations, seasonal fluctuations and the time-dependent decline in reactivity were also captured. The modeling results generally confirm, and further strengthen, the existing conceptual model for the site. Overall sulfate reduction and S-accumulation rates are constrained with confidence within a factor of 1.5.     

Influence factors of multicomponent mixtures containing reactive chemicals and their joint effects

Organic chemicals usually coexist as a mixture in the environment, and the mixture toxicity of organic chemicals has received increased attention. However, research regarding the joint effects of reactive chemicals is lacking. In this study, we examined two kinds of reactive chemicals, cyanogenic toxicants and aldehydes and determined their joint effects on Photobacterium phosphoreum. Three factors were found to influence the joint effects of multicomponent mixtures containing reactive chemicals, including the number of components, the dominating components and the toxic ratios. With an increased number of components, the synergistic or antagonistic effects (interactions) will weaken to the additive effects (non-interactions) if the added component cannot yield a much stronger joint effect with an existing component. Contrarily, the joint effect of the mixture may become stronger instead of weaker if the added components can yield a much stronger joint effect than the existing joint effect of the multicomponent mixture. The components that yield the strongest interactions in their binary mixture can be considered the dominating components. These components contribute more to the interactions of multicomponent mixtures than other components. Moreover, the toxic ratios also influence the joint effects of the mixtures. This study provides an insight into what are the main factors and how they influence the joint effects of multicomponent mixtures containing reactive chemicals, and thus, the findings are beneficial to the study of mixture toxicology.     

Operator-splitting errors in coupled reactive transport codes for transient variably saturated flow and contaminant transport in layered soil profiles

One possible way of integrating subsurface flow and transport processes with (bio)geochemical reactions is to couple by means of an operator-splitting approach two completely separate codes, one for variably-saturated flow and solute transport and one for equilibrium and kinetic biogeochemical reactions. This paper evaluates the accuracy of the operator-splitting approach for multicomponent systems for typical soil environmental problems involving transient atmospheric boundary conditions (precipitation, evapotranspiration) and layered soil profiles. The recently developed HP1 code was used to solve the coupled transport and chemical equations. For steady-state flow conditions, the accuracy was found to be mainly a function of the adopted spatial discretization and to a lesser extent of the temporal discretization. For transient flow situations, the accuracy depended in a complex manner on grid discretization, time stepping and the main flow conditions (infiltration versus evaporation). Whereas a finer grid size reduced the numerical errors during steady-state flow or the main infiltration periods, the errors sometimes slightly increased (generally less than 50%) when a finer grid size was used during periods with a high evapotranspiration demand (leading to high pressure head gradients near the soil surface). This indicates that operator-splitting errors are most significant during periods with high evaporative boundary conditions. The operator-splitting errors could be decreased by constraining the time step using the performance index (the product of the grid Peclet and Courant numbers) during infiltration, or the maximum time step during evapotranspiration. Several test problems were used to provide guidance for optimal spatial and temporal discretization.     

Theory for reactive solute transport through clay membrane barriers

The theoretical development for one-dimensional, coupled migration of solutes with different ionic mobilities through clay soils that behave as ion-restrictive membranes, referred to as clay membrane barriers (CMBs), is presented. The transport formulation is based on principles of irreversible thermodynamics and accounts explicitly for coupling effects of hyperfiltration (ultrafiltration) and chemico-osmotic counter-advection associated with clay membrane behavior in the absence of electrical current. Since, by definition, no solute can enter a "perfect" or "ideal" membrane, the concept of an implicit coupling effect, such that the effective salt-diffusion coefficient, Ds* approaches zero as the chemico-osmotic efficiency coefficient, omega approaches unity is introduced. The theoretical development also illustrates that, even in the absence of membrane behavior, traditional advective-dispersive transport theory based on a constant value of Ds* for the solutes may not be appropriate for simulating transient transport in reactive (ion exchanging) systems. This potential limitation is illustrated through simulations for solute mass flux involving the migration of a binary salt solution (KCl) through a clay barrier with exchange sites saturated with a single exchangeable cation (e.g., Na+) that enters the pore solution upon ion exchange with the salt cation (K+).     

Numerical modeling of coupled variably saturated fluid flow and reactive transport with fast and slow chemical reactions

The couplings among chemical reaction rates, advective and diffusive transport in fractured media or soils, and changes in hydraulic properties due to precipitation and dissolution within fractures and in rock matrix are important for both nuclear waste disposal and remediation of contaminated sites. This paper describes the development and application of LEHGC2.0, a mechanistically based numerical model for simulation of coupled fluid flow and reactive chemical transport, including both fast and slow reactions in variably saturated media. Theoretical bases and numerical implementations are summarized, and two example problems are demonstrated. The first example deals with the effect of precipitation/dissolution on fluid flow and matrix diffusion in a two-dimensional fractured media. Because of the precipitation and decreased diffusion of solute from the fracture into the matrix, retardation in the fractured medium is not as large as the case wherein interactions between chemical reactions and transport are not considered. The second example focuses on a complicated but realistic advective-dispersive-reactive transport problem. This example exemplifies the need for innovative numerical algorithms to solve problems involving stiff geochemical reactions.     

Large-scale modeling of reactive solute transport in fracture zones of granitic bedrocks

Final disposal of high-level radioactive waste in deep repositories located in fractured granite formations is being considered by several countries. The assessment of the safety of such repositories requires using numerical models of groundwater flow, solute transport and chemical processes. These models are being developed from data and knowledge gained from in situ experiments such as the Redox Zone Experiment carried out at the underground laboratory of Äspö in Sweden. This experiment aimed at evaluating the effects of the construction of the access tunnel on the hydrogeological and hydrochemical conditions of a fracture zone intersected by the tunnel. Most chemical species showed dilution trends except for bicarbonate and sulphate which unexpectedly increased with time. Molinero and Samper [Molinero, J. and Samper, J. Groundwater flow and solute transport in fracture zones: an improved model for a large-scale field experiment at Äspö (Sweden). J. Hydraul. Res., 42, Extra Issue, 157–172] presented a two-dimensional water flow and solute transport finite element model which reproduced measured drawdowns and dilution curves of conservative species. Here we extend their model by using a reactive transport which accounts for aqueous complexation, acid–base, redox processes, dissolution–precipitation of calcite, quartz, hematite and pyrite, and cation exchange between Na⁺ and Ca²⁺. The model provides field-scale estimates of cation exchange capacity of the fracture zone and redox potential of groundwater recharge. It serves also to identify the mineral phases controlling the solubility of iron. In addition, the model is useful to test the relevance of several geochemical processes. Model results rule out calcite dissolution as the process causing the increase in bicarbonate concentration and reject the following possible sources of sulphate: (1) pyrite dissolution, (2) leaching of alkaline sulphate-rich waters from a nearby rock landfill and (3) dissolution of iron monosulphides contained in Baltic seafloor sediments. Based on these results, microbially mediated processes are postulated as the most likely hypothesis to explain the measured increase of dissolved bicarbonates and sulphates after tunnel construction.     

A sensibility analysis of model selection in modeling the reactive transport of cesium in crushed granite

We performed a sensibility analysis of model selection in modeling the reactive transport of cesium in crushed granite through model calibration and validation. Based on some solid phase analysis data and kinetic batch experimental results, we hypothesized three two-site sorption models in the LEHGC reactive transport model to fit the breakthrough curves (BTCs) from the corresponding column experiments. The analysis of breakthrough curves shows that both the empirical two-site kinetic linear sorption model and the semi-mechanistic/semi-empirical two-site kinetic surface complexation model, regardless of their complexity, can match our experimental data fairly well under given test conditions. A numerical experiment to further compare the two models shows that they behave differently when the pore velocity is not of the same order of magnitude as our test velocities. This result indicates that further investigations to help determine a better model are needed. We suggest that a multistage column experiment, which tests over the whole range of practical flow velocities, should be conducted to help alleviate inadequate hypothesized models.     

Visual modflow在模拟实验室含水砂槽中柴油运移特征中的应用

Visual modflow是一个可以对三维地下水水流和溶质运移进行数值模拟评价的标准可视化专业软件.建立了砂箱物理实验模型来研究柴油在含水砂槽中的迁移特征.通过模型检验,各个监测时期观测值和预测值相关系数r值在0.564 ～0.669之间,证明这种建模方法是合理的和有效的.利用校正的模型对实验室含水砂槽中柴油运移特征进行模拟,发现所建模型可以较为准确地反映出含水砂槽中柴油污染物的分布特征,拟合、验证和预测结果显示该模型可作为地下水管理的有效工具,这为深入研究柴油污染地下水提供理论依据.     

Two-site kinetic modeling of bacteriophages transport through columns of saturated dune sand

Breakthrough curves, on a semi-log scale, from tests in porous media with block-input of viruses, bacteria, protozoa and colloidal particles often exhibit a typical skewness: a rather slowly rising limb and a smooth transition of a declining limb to a very long tail. One-site kinetic models fail to fit the rising and declining limbs together with the tail satisfactorily. Inclusion of an equilibrium adsorption site does not seem to improve simulation results. This was encountered in the simulation of breakthrough curves from a recent field study on the removal of bacteriophages MS2 and PRD1 by passage through dune sand. In the present study, results of laboratory experiments for the study of this issue are presented. Breakthrough curves of salt and bacteriophages MS2, PRDI, and phiX174 in 1 D column experiments have been measured. One- and two-site kinetic models have been applied to fit and predict breakthrough curves from column experiments. The two-site model fitted all breakthrough curves very satisfactorily, accounting for the skewness of the rising limb as well as for the smooth transition of the declining limb to the tail of the breakthrough curve. The one-site model does not follow the curvature of the breakthrough tail, leading to an overestimation of the inactivation rate coefficient for attached viruses. Interaction with kinetic site 1 is characterized by relatively fast attachment and slow detachment, whereas attachment to and detachment from kinetic site 2 is fast. Inactivation of viruses and interaction with kinetic site 2 provide only a minor contribution to removal. Virus removal is mainly determined by the attachment to site 1. Bacteriophage phiX174 attached more than MS2 and PRD1, which can be explained by the greater electrostatic repulsion that MS2 and PRD1 experience compared to the less negatively charged phiX174.     

Spatial and temporal features of density-dependent contaminant transport: Experimental investigation and numerical modeling

We investigate the spatial and temporal features of variable-density contaminant plumes migration in porous materials. Our analysis is supported by novel experimental results concerning concentration profiles inside a vertical column setup that has been conceived at CEA to this aim. The experimental method relies on X-ray spectrometry, which allows determining solute profiles as a function of time at several positions along the column. The salient outcomes of the measurements are elucidated, with focus on miscible fluids in homogeneous saturated media. The role of the injected solution molarity is evidenced. As molarity increases, the solutes plume transport progressively deviates from the usual Fickian behavior, and pollutants distribution becomes skewed in the direction dictated by gravity. By resorting to a finite elements approach, we numerically solve the nonlinear equations that rule the pollutants migration: a good agreement is found between the simulated profiles and the experimental data. At high molarity, a strong dependence on initial conditions is found. Finally, we qualitatively explore the (unstable) interfacial dynamics between the dense contaminant plume and the lighter resident fluid that saturates the column, and detail its evolution for finite-duration contaminant injections.     

Evaluating the sensitivity of a subsurface multicomponent reactive transport model with respect to transport and reaction parameters.

The input variables for a numerical model of reactive solute transport in groundwater include both transport parameters, such as hydraulic conductivity and infiltration, and reaction parameters that describe the important chemical and biological processes in the system. These parameters are subject to uncertainty due to measurement error and due to the spatial variability of properties in the subsurface environment. This paper compares the relative effects of uncertainty in the transport and reaction parameters on the results of a solute transport model. This question is addressed by comparing the magnitudes of the local sensitivity coefficients for transport and reaction parameters. General sensitivity equations are presented for transport parameters, reaction parameters, and the initial (background) concentrations in the problem domain. Parameter sensitivity coefficients are then calculated for an example problem in which uranium(VI) hydrolysis species are transported through a two-dimensional domain with a spatially variable pattern of surface complexation sites. In this example, the reaction model includes equilibrium speciation reactions and mass transfer-limited non-electrostatic surface complexation reactions. The set of parameters to which the model is most sensitive includes the initial concentration of one of the surface sites, the formation constant (Kf) of one of the surface complexes and the hydraulic conductivity within the reactive zone. For this example problem, the sensitivity analysis demonstrates that transport and reaction parameters are equally important in terms of how their variability affects the model results.     

Microbial growth and transport in porous media under denitrification conditions: experiments and simulations

Soil column experiments were conducted to study bacterial growth and transport in porous media under denitrifying conditions. The study used a denitrifying microbial consortium isolated from aquifer sediments sampled at the U.S. Department of Energy's Hanford site. One-dimensional, packed-column transport studies were conducted under two substrate loading conditions. A detailed numerical model was developed to predict the measured effluent cell and substrate concentration profiles. First-order attachment and detachment models described the interphase exchange processes between suspended and attached biomass. Insignificantly different detachment coefficient values of 0.32 and 0.43 day⁻¹, respectively, were estimated for the high and low nitrate loading conditions (48 and 5 mg l⁻¹ NO₃, respectively). Comparison of these values with those calculated from published data for aerobically growing organisms shows that the denitrifying consortium had lower detachment rate coefficients. This suggests that, similar to detachment rates in reactor-grown biofilms, detachment in porous media may increase with microbial growth rate. However, available literature data are not sufficient to confirm a specific analytical model for predicting this growth dependence.     

Nonideal transport of reactive contaminants in heterogeneous porous media: 7. distributed-domain model incorporating immiscible-liquid dissolution and rate-limited sorption/desorption

The purpose of this work is to present a distributed-domain mathematical model incorporating the primary mass-transfer processes that mediate the transport of immiscible organic liquid constituents in water-saturated, locally heterogeneous porous media. Specifically, the impact of grain/pore-scale heterogeneity on immiscible-liquid dissolution and sorption/desorption is represented in the model by describing the system as comprising a continuous distribution of mass-transfer domains. With this conceptualization, the distributions of the initial dissolution rate coefficient and the sorption/desorption rate coefficient are represented as probability density functions. Several sets of numerical experiments are conducted to examine the effects of heterogeneous dissolution and sorption/desorption on contaminant transport and elution. Four scenarios with different combinations of uniform/heterogeneous rate-limited dissolution and uniform/heterogeneous rate-limited sorption/desorption are evaluated. The results show that both heterogeneous rate-limited sorption/desorption and heterogeneous rate-limited dissolution can significantly increase the time or pore volumes required to elute immiscible-liquid constituents from a contaminated porous medium. However, sorption/desorption has minimal influence on elution behavior until essentially all of the immiscible liquid has been removed. For typical immiscible-liquid constituents that have relatively low sorption, the asymptotic elution tailing produced by heterogeneous rate-limited sorption/desorption begins at effluent concentrations that are several orders of magnitude below the initial steady-state concentrations associated with dissolution of the immiscible liquid. Conversely, the enhanced elution tailing associated with heterogeneous rate-limited dissolution begins at concentrations that are approximately one-tenth of the initial steady-state concentrations. Hence, dissolution may generally control elution behavior of immiscible-liquid constituents in cases wherein grain/pore-scale heterogeneity significantly influences both dissolution and sorption/desorption.     

Operator-splitting procedures for reactive transport and comparison of mass balance errors

Operator-splitting (OS) techniques are very attractive for numerical modelling of reactive transport, but they induce some errors. Considering reactive mass transport with reversible and irreversible reactions governed by a first-order rate law, we develop analytical solutions of the mass balance for the following operator-splitting schemes: standard sequential non-iterative (SNI), Strang-splitting SNI, standard sequential iterative (SI), extrapolating SI, and symmetric SI approaches. From these analytical solutions, the operator-splitting methods are compared with respect to mass balance errors and convergence rates independently of the techniques used for solving each operator. Dimensionless times, NOS, are defined. They control mass balance errors and convergence rates. The following order in terms of decreasing efficiency is proposed: symmetric SI, Strang-splitting SNI, standard SNI, extrapolating SI and standard SI schemes. The symmetric SI scheme does not induce any operator-splitting errors, the Strang-splitting SNI appears to be O(N2OS) accurate, and the other schemes are first-order accurate.     

Influence of indoor transport and mixing time scales on the performance of sensor systems for characterizing contaminant releases

Optimizing real-time sensor systems to detect and identify relevant characteristics of an indoor contaminant event is a challenging task. The interpretation of incoming sensor data is confounded by uncertainties in building operation, in the forces driving contaminant transport, and in the physical parameters governing transport. In addition, simulation tools used by the sensor interpretation algorithm introduce modeling uncertainties. This paper explores how the time scales inherent in contaminant transport influence the information that can be extracted from real-time sensor data. In particular, we identify three time scales (within room mixing, room-to-room transport, and removal from the building) and study how they affect the ability of a Bayesian Monte Carlo (BMC) sensor interpretation algorithm to identify the release location and release mass from a set of experimental data, recorded in a multi-floor building. The research shows that some limitations in the BMC approach do not depend on details of the models or the algorithmic implementation, but rather on the physics of contaminant transport. These inherent constraints have implications for the design of sensor systems.