No abstract

The use of zerovalent iron micro- and nanoparticles (MZVI and NZVI) for groundwater remediation is hindered by colloidal instability, causing aggregation (for NZVI) and sedimentation (for MZVI) of the particles. Transportability of MZVI and NZVI in porous media was previously shown to be significantly increased if viscous shear-thinning fluids (xanthan gum solutions) are used as carrier fluids. In this work, a novel modeling approach is proposed and applied for the simulation of 1D flow and transport of highly concentrated (20 g/L) non-newtonian suspensions of MZVI and NZVI, amended with xanthan gum (3 g/L). The coupled model is able to simulate the flow of a shear thinning fluid including the variable apparent viscosity arising from changes in xanthan and suspended iron particle concentrations. The transport of iron particles is modeled using a dual-site approach accounting for straining and physicochemical deposition/release phenomena. A general formulation for reversible deposition is herein proposed, that includes all commonly applied dynamics (linear attachment, blocking, ripening). Clogging of the porous medium due to deposition of iron particles is modeled by tying porosity and permeability to deposited iron particles. The numerical model proved to adequately fit the transport tests conducted using both MZVI and NZVI and can develop into a powerful tool for the design and the implementation of full scale zerovalent iron applications.

The use of nanoscale ferrihydrite particles, which are known to effectively enhance microbial degradation of a wide range of contaminants, represents a promising technology for in situ remediation of contaminated aquifers. Thanks to their small size, ferrihydrite nanoparticles can be dispersed in water and directly injected into the subsurface to create reactive zones where contaminant biodegradation is promoted. Field applications would require a detailed knowledge of ferrihydrite transport mechanisms in the subsurface, but such studies are lacking in the literature. The present study is intended to fill this gap, focusing in particular on the influence of flow rate and ionic strength on particle mobility. Column tests were performed under constant or transient ionic strength, including injection of ferrihydrite colloidal dispersions, followed by flushing with particle-free electrolyte solutions. Particle mobility was greatly affected by the salt concentration, and particle retention was almost irreversible under typical salt content in groundwater. Experimental results indicate that, for usual ionic strength in European aquifers (2 to 5 mM), under natural flow condition ferrihydrite nanoparticles are likely to be transported for 5 to 30 m. For higher ionic strength, corresponding to contaminated aquifers, (e.g., 10 mM) the travel distance decreases to few meters. A simple relationship is proposed for the estimation of travel distance with changing flow rate and ionic strength. For future applications to aquifer remediation, ionic strength and injection rate can be used as tuning parameters to control ferrihydrite mobility in the subsurface and therefore the radius of influence during field injections.

Flow of non-Newtonian fluids through porous media at high Reynolds numbers is often encountered in chemical, pharmaceutical and food as well as petroleum and groundwater engineering and in many other industrial applications. Under the majority of operating conditions typically explored, the dependence of pressure drops on flow rate is non-linear and the development of models capable of describing accurately this dependence, in conjunction with non trivial rheological behaviours, is of paramount importance. In this work pore-scale single-phase flow simulations conducted on synthetic two-dimensional porous media are performed via computational fluid dynamics for both Newtonian and non-Newtonian fluids, and the results are used for the extension and validation of the Darcy-Forchheimer law, herein proposed for shear thinning fluid models of Cross, Ellis and Carreau. The inertial parameter β is demonstrated to be independent of the viscous properties of the fluids. The results of flow simulations show the superposition of two contributions to pressure drops: one, strictly related to the non-Newtonian properties of the fluid, dominates at low Reynolds numbers, while a quadratic one, arising at higher Reynolds numbers, is dependent on the porous medium properties. The use of pore-scale flow simulations on limited portions of the porous medium is here proposed for the determination of the macroscale averaged parameters (permeability K, inertial coefficient β and shift factor α), which are required for the estimation of pressure drops via the extended Darcy-Forchheimer law. The method can be applied for those fluids which would lead to critical conditions (high pressures for low permeability media and/or high flow rates) in laboratory tests.

Pore-scale simulation of fluid flow and solute dispersion in three-dimensional porous media In the present work fluid flow and solute transport through porous media are described by solving the governing equations at the pore scale with finite-volume discretization. Instead of solving the simplified Stokes equation (very often employed in this context) the full Navier-Stokes equation is used here. The realistic three-dimensional porous medium is created in this work by packing together, with standard ballistic physics, irregular and polydisperse objects. Emphasis is placed on numerical issues related to mesh generation and spatial discretization, which play an important role in determining the final accuracy of the finite-volume scheme and are often overlooked. The simulations performed are then analyzed in terms of velocity distributions and dispersion rates in a wider range of operating conditions, when compared with other works carried out by solving the Stokes equation. Results show that dispersion within the analyzed porous medium is adequately described by classical power laws obtained by analytic homogenization. Eventually the validity of Fickian diffusion to treat dispersion in porous media is also assessed. KAUST Repository

No abstract

In the present work column transport tests were performed in order to study the mobility of guar-gum suspensions of microscale zero-valent iron particles (MZVI) in porous media. The results were analyzed with the purpose of implementing a radial model for the design of full scale interventions. The transport tests were performed using several concentrations of shear thinning guar gum solutions as stabilizer (1.5, 3 and 4g/l) and applying different flow rates (Darcy velocity in the range 1·10(-4) to 2·10(-3)m/s), representative of different distances from the injection point in the radial domain. Empirical relationships, expressing the dependence of the deposition and release parameters on the flow velocity, were derived by inverse fitting of the column transport tests using a modified version of E-MNM1D (Tosco and Sethi, 2010) and the user interface MNMs (www.polito.it/groundwater/software). They were used to develop a comprehensive transport model of MZVI suspensions in radial coordinates, called E-MNM1R, which takes into account the non Newtonian (shear thinning) rheological properties of the dispersant fluid and the porous medium clogging associated with filtration and sedimentation in the porous medium of both MZVI and guar gum residual undissolved particles. The radial model was run in forward mode to simulate the injection of MZVI dispersed in guar gum in conditions similar to those applied in the column transport tests. In a second stage, we demonstrated how the model can be used as a valid tool for the design and the optimization of a full scale intervention. The simulation results indicated that several concurrent aspects are to be taken into account for the design of a successful delivery of MZVI/guar gum slurries via permeation injection, and a compromise is necessary between maximizing the radius of influence of the injection and minimizing the injection pressure, to guarantee a sufficiently homogeneous distribution of the particles around the injection point and to prevent preferential flow paths.

No abstract