Fritjof Fagerlund

Empirical correlations to estimate agglomerate size and deposition during injection of a polyelectrolyte-modified Fe0 nanoparticle at high particle concentration in saturated sand

Controlled emplacement of polyelectrolyte-modified nanoscale zerovalent iron (NZVI) particles at high particle concentration (1-10 g/L) is needed for effective in situ subsurface remediation using NZVI. Deep bed filtration theory cannot be used to estimate the transport and deposition of concentrated polyelectrolyte-modified NZVI dispersions (>0.03 g/L) because particles agglomerate during transport which violates a fundamental assumption of the theory. Here we develop two empirical correlations for estimating the deposition and transport of concentrated polyelectrolyte-modified NZVI dispersions in saturated porous media when NZVI agglomeration in porous media is assumed to reach steady state quickly. The first correlation determines the apparent stable agglomerate size formed during NZVI transport in porous media for a fixed hydrogeochemical condition. The second correlation estimates the attachment efficiency (sticking coefficient) of the stable agglomerates. Both correlations are described using dimensionless numbers derived from parameters affecting deposition and agglomeration in porous media. The exponents for the dimensionless numbers are determined from statistical analysis of breakthrough data for polyelectrolyte-modified NZVI dispersions collected in laboratory scale column experiments for a range of ionic strength (1, 10, and 50mM Na(+) and 0.25, 1, and 1.25 mM Ca(2+)), approach velocity (0.8 to 55 × 10(-4)m/s), average collector sizes (d(50)=99 μm, 300 μm, and 880 μm), and polyelectrolyte surface modifier properties. Attachment efficiency depended on approach velocity and was inversely related to collector size, which is contrary to that predicted from classic filtration models. High ionic strength, the presence of divalent cations, lower extended adsorbed polyelectrolyte layer thickness, decreased approach velocity, and a larger collector size promoted NZVI agglomeration and deposition and thus limited its mobility in porous media. These effects are captured quantitatively in the two correlations developed. The application and limitations of using the correlations for preliminary design of in situ NZVI emplacement strategies is discussed.

Polymer-modified Fe0 nanoparticles target entrapped NAPL in two dimensional porous media: effect of particle concentration, NAPL saturation, and injection strategy.

Polymer-modified nanoscale zerovalent iron (NZVI) particles are delivered into porous media for in situ remediation of nonaqueous phase liquid (NAPL) source zones. A systematic and quantitative evaluation of NAPL targeting by polymer-modified NZVI in two-dimensional (2-D) porous media under field-relevant conditions has not been reported. This work evaluated the importance of NZVI particle concentration, NAPL saturation, and injection strategy on the ability of polymer-modified NZVI (MRNIP2) to target the NAPL/water interface in situ in a 2-D porous media model. Dodecane was used as a NAPL model compound for this first demonstration of source zone targeting in 2-D. A driving force for NAPL targeting, the surface activity of MRNIP2 at the NAPL/water interface was verified ex situ by its ability to emulsify NAPL in water. MRNIP2 at low particle concentration (0.5 g/L) did not accumulate in or near entrapped NAPL, however, MRNIP2 at moderate and high particle concentrations (3 and 15 g/L) did accumulate preferentially at entrapped NAPL, i.e., it was capable of in situ targeting. The amount of MRNIP2 that targets a NAPL source depends on NAPL saturation (S(n)), presumably because the saturation controls the available NAPL/water interfacial area and the flow field through the NAPL source. At effective S(n) close or equal to 100%, MRNIP2 bypassed NAPL and accumulated only at the periphery of the entrapped NAPL region. At lower S(n), flow also carries MRNIP2 to NAPL/water interfaces internal to the entrapped NAPL region. However, the mass of accumulated MRNIP2 per unit available NAPL/water interfacial area is relatively constant (∼0.8 g/m(2) for MRNIP2 = 3 g/L) from S(n) = 13 to ∼100%, suggesting that NAPL targeting is mostly controlled by MRNIP2 sorption onto the NAPL/water interface.

PCE dissolution and simultaneous dechlorination by nanoscale zero-valent iron particles in a DNAPL source zone.

While the capability of nanoscale zero-valent iron (NZVI) to dechlorinate organic compounds in aqueous solutions has been demonstrated, the ability of NZVI to remove dense non-aqueous phase liquid (DNAPL) from source zones under flow-through conditions similar to a field scale application has not yet been thoroughly investigated. To gain insight on simultaneous DNAPL dissolution and NZVI-mediated dechlorination reactions after direct placement of NZVI into a DNAPL source zone, a combined experimental and modeling study was performed. First, a DNAPL tetrachloroethene (PCE) source zone with emplaced NZVI was built inside a small custom-made flow cell and the effluent PCE and dechlorination byproducts were monitored over time. Second, a model for rate-limited DNAPL dissolution and NZVI-mediated dechlorination of PCE to its three main reaction byproducts with a possibility for partitioning of these byproducts back into the DNAPL was formulated. The coupled processes occurring in the flow cell were simulated and analyzed using a detailed three-dimensional numerical model. It was found that subsurface emplacement of NZVI did not markedly accelerate DNAPL dissolution or the DNAPL mass-depletion rate, when NZVI at a particle concentration of 10g/L was directly emplaced in the DNAPL source zone. To react with NZVI the DNAPL PCE must first dissolve into the groundwater and the rate of dissolution controls the longevity of the DNAPL source. The modeling study further indicated that faster reacting particles would decrease aqueous contaminant concentrations but there is a limit to how much the mass removal rate can be increased by increasing the dechlorination reaction rate. To ensure reduction of aqueous contaminant concentrations, remediation of DNAPL contaminants with NZVI should include emplacement in a capture zone down-gradient of the DNAPL source.

Effects of single-fracture aperture statistics on entrapment, dissolution and source depletion behavior of dense non-aqueous phase liquids.

Understanding of the entrapment and dissolution behavior of dense non-aqueous phase liquids (DNAPLs) in single fractures is important for modeling contaminant flux generation from fractured sites. Here a systematic numerical study is presented to investigate the effect of fracture aperture statistics on DNAPL migration, entrapment and dissolution within individual, variable-aperture fractures. Both fractures with open and closed bottom boundaries were considered. For the simulation a continuum-based two-phase model was used with a capillary pressure function which calculates the entry pressure based on the local aperture. Prior to application the model was compared against the invasion percolation approach and found more suitable for the present study, in particular as it allows a more versatile presentation of boundary conditions. The results showed that increasing aperture standard deviation and/or decreasing correlation length lead to larger amounts of entrapped DNAPL (due to the fact that larger standard deviation produces more distinct contrast between small and large aperture regions and the fact that longer correlation length provides more possible channels through the fracture) as well as larger maximum and average sizes of DNAPL blobs, and subsequently lead to longer times for complete dissolution. To understand the relationship between the solute flux and the remaining mass, a simplified source depletion function which links the outflow concentration to the DNAPL saturation was found adequate to describe the dissolution process for the case where the bottom boundary is open for DNAPL migration and thus the DNAPL does not accumulate to form a pool. The parameters in this function were not very sensitive to variations in correlation length but were sensitive to aperture standard deviation. The same average entrapped DNAPL saturation produced considerably smaller solute concentrations in cases with larger aperture variability due to the larger average size of DNAPL blobs (i.e., smaller contact area for DNAPL dissolution). Boundary conditions had a significant impact on DNAPL entrapment and dissolution. A closed boundary at the bottom led to DNAPL pooling (i.e., large continuous blobs) which causes significant tailing in the dissolution breakthrough curve due to water bypassing.

Modeling of nonaqueous phase liquid (NAPL) migration in heterogeneous saturated media: Effects of hysteresis and fluid immobility in constitutive relations

The confidence in model predictions for nonaqueous phase liquid (NAPL) transport in stochastically heterogeneous systems is limited. The fundamental approaches as well as the constitutive models ha ...

Dissolution of dense non-aqueous phase liquids in vertical fractures : Effect of finger residuals and dead-end pools

Understanding the dissolution behavior of dense non-aqueous phase liquids (DNAPLs) in rock fractures under different entrapment conditions is important for remediation activities and any related predictive modeling. This study investigates DNAPL dissolution in variable aperture fractures under two important entrapment configurations, namely, entrapped residual blobs from gravity fingering and pooling in a dead-end fracture. We performed a physical dissolution experiment of residual DNAPL blobs in a vertical analog fracture using light transmission techniques. A high-resolution mechanistic (physically-based) numerical model has been developed which is shown to excellently reproduce the experimentally observed DNAPL dissolution. We subsequently applied the model to simulate dissolution of the residual blobs under different water flushing velocities. The simulated relationship between the Sherwood number Sh and Peclet number Pe could be well fitted with a simple power-law function (Sh=1.43Pe⁰·⁴³). To investigate mass transfer from dead-end pools, another type of trapping in rock fractures, entrapment and dissolution of DNAPL in a vertical dead-end fracture was simulated. As the entrapped pool dissolves, the depth of the interface between the DNAPL and the flowing water increases linearly with decreasing DNAPL saturation. The interfacial area remains more or less constant as DNAPL saturation decreases, unlike in the case of residual DNAPL blobs. The decreasing depth of the contact interface changes the flow field and causes decreasing water flow velocity above the top of the DNAPL pool, suggesting the dependence of the mass transfer rate on the depth of the interface, or alternatively, the remaining mass percentage in the fracture. Simulation results show that the resultant Sherwood number Sh is significantly smaller than in the case of residual blobs for any given Peclet number, indicating slower mass transfer. The results also show that the Sh can be well fitted with a power-law function of Pe and remaining mass percentage. The obtained relationships of dimensionless groups concerning the mass transfer characteristics at the level of individual fractures can be further used in predictive modeling of dissolution at a larger (fracture network) scale.

Refractive‐Light‐Transmission Technique Applied to Density‐Driven Convective Mixing in Porous Media With Implications for Geological CO2 Storage

Density-driven convection has been identified to accelerate the rate of CO2 solubility trapping during geological CO2 storage in deep saline aquifers. In this paper, we present an experimental meth ...

Spectral induced polarization of Na-montmorillonite dispersions

Montmorillonite (Mt) clays have a high specific surface area and surface charge, which confer them remarkable adsorption properties. Nevertheless, their electrochemical and aggregation behavior are not completely elucidated because of the complexity of their microstructural and interfacial properties. In this work, the conductive and dispersive properties of Na-Mt suspensions of weight fractions 0.5-5.2% were investigated for the first time using the spectral induced polarization method. A four-electrode system was used to reduce errors introduced by electrode polarization and contact resistances. Complex conductivity spectra in the low-frequency range of 0.1Hz to 45kHz were successfully described using a triple layer model of the basal surface of Mt and a complex conductivity model that considers conduction of the diffuse layer and polarization of the Stern layer. Aggregate size distributions were inferred from inverted relaxation time distributions. We found that the negative and permanent surface charge of the basal plane of Na-Mt controls its quadrature (imaginary) conductivity, which is not very sensitive to pH and salinity (NaCl) in the 100Hz to 45kHz frequency range. For lower frequencies, the sudden increase of the quadrature conductivity at the highest salinities was explained by considering coagulation of Na-Mt particles.

Transport of surface-modified carbon nanotubes through a soil column.

Carbon nanotubes (CNTs) are widely manufactured nanoparticles, which are being utilized in a number of consumer products, such as sporting goods, electronics and biomedical applications. Due to their accelerating production and use, CNTs constitute a potential environmental risk if they are released to soil and groundwater systems. It is therefore essential to improve the current understanding of environmental fate and transport of CNTs. The transport and retention of CNTs in both natural and artificial media have been reported in literature, but the findings widely vary and are thus not conclusive. There are a number of physical and chemical parameters responsible for variation in retention and transport. In this study, a complete procedure of selected multiwalled carbon nanotubes (MWCNTs) is presented starting from their surface modification to a complete set of laboratory column experiments at critical physical and chemical scenarios. Results indicate that the stability of the commercially available MWCNTs are critical with their attached surface functional group which can also influence the transport and retention of MWCNT through the surrounding medium.

Detecting Sub-Surface Groundwater Flow in Fractured Rock Using Self- Potential (SP) Methods

In the wine district of Clare Valley in South Australia, the natural-voltage SP-signal generated by fluid flow in fractured rock during a pumping test was carefully monitored over time. From ten electrode-locations surrounding the pumping well logged every 15 s, the drawdown cone produced by pumping was determined on basis of the SP measurements together with laboratory measurements of the grain-boundary ζ (zeta) potential. Such measurements allowed calculation of the fractured-rock aquifer’s transmissivity and average permeability. Results were confirmed by piezometer measurements to the extent that data were available. The study has revealed that SP-signals generated during pumping tests are of complex nature. However, if the pumping test is sufficiently long to allow the signal to stabilise, and careful field procedures are in place, then SP-measurements have the potential to add significant hydrogeological information. SP-measurements are relatively easy and cheap, and are, contrary to traditional hydrogeological methods, not restricted to the locations of existing piezometers, which is particularly useful in fractured-rock aquifers.