Saeed Torkzaban

A Theoretical Analysis of Colloid Attachment and Straining in Chemically Heterogeneous Porous Media

A balance of applied hydrodynamic (T(H)) and resisting adhesive (T(A)) torques was conducted over a chemically heterogeneous porous medium that contained random roughness of height h(r) to determine the fraction of the solid surface area that contributes to colloid immobilization (S(f)*) under unfavorable attachment conditions. This model considers resistance due to deformation and the horizontal component of the adhesive force (F(AT)), spatial variations in the pore scale velocity distribution, and the influence of hr on lever arms for T(H) and T(A). Values of S(f)* were calculated for a wide range of physicochemical properties to gain insight into mechanisms and factors influencing colloid immobilization. Colloid attachment processes were demonstrated to depend on solution ionic strength (IS), the colloid radius (r(c)), the Youngs modulus (K), the amount of chemical heterogeneity (P+), and the Darcy velocity (q). Colloid immobilization was also demonstrated to occur on a rough surface in the absence of attachment. In this case, S(f)* depended on IS, r(c), the roughness fraction (f0), h(r), and q. Roughness tended to enhance T(A) and diminish T(H). Consequently, the effect of IS on S(f)* was enhanced by h(r) relative to attachment. In contrast, the effects of r(c) and q on S(f)* were diminished by hr in comparison to attachment. Colloid immobilization adjacent to macroscopic roughness locations shares many similarities to grain-grain contact points and may be viewed as a type of straining process. In general, attachment was more important for higher IS and variance in the secondary minimum, and for smaller r(c), q, and K, but diffusion decreased these values. Conversely, straining was dominant for the opposite conditions. Discrepancies in the literature on mechanisms of colloid retention are likely due to a lack of consideration of all of these factors.

Colloid interaction energies for physically and chemically heterogeneous porous media.

The mean and variance of the colloid interaction energy (Φ*) as a function of separation distance (h) were calculated on physically and/or chemically heterogeneous solid surfaces at the representative elementary area (REA) scale. Nanoscale roughness was demonstrated to have a significant influence on the colloid interaction energy for different ionic strengths. Increasing the roughness height reduced the magnitude of the energy barrier (Φmax*) and the secondary minimum (Φ2min*). Conversely, increasing the fraction of the solid surface with roughness increased the magnitude of Φmax* and Φ2min*. Our results suggest that primary minimum interactions tend to occur in cases where only a portion of the solid surface was covered with roughness (i.e., isolated roughness pillars), but their depths were shallow as a result of Born repulsion. The secondary minimum was strongest on smooth surfaces. The variance in the interaction energy was also a strong function of roughness parameters and h. In particular, the variance tended to increase with the colloid size, the magnitude of Φ*, the height of the roughness, and especially the size (cross-sectional area) of the heterogeneity. Nonzero values of the variance for Φ2min* implied the presence of a tangential component of the adhesive force and a resisting torque that controls immobilization and release for colloids at this location. Heterogeneity reduced the magnitude of Φ* in comparison to the corresponding homogeneous situation. Physical heterogeneity had a greater influence on mean properties of Φ* than similar amounts of chemical heterogeneity, but the largest reduction occurred on surfaces with both physical and chemical heterogeneity. The variance in Φ* tended to be higher for a chemically heterogeneous solid.

Hysteresis of colloid retention and release in saturated porous media during transients in solution chemistry.

Saturated packed column and micromodel transport studies were conducted to gain insight on mechanisms of colloid retention and release under unfavorable attachment conditions. The initial deposition of colloids in porous media was found to be a strongly coupled process that depended on solution chemistry and pore space geometry. During steady state chemical conditions, colloid deposition was not a readily reversible process, and micromodel photos indicated that colloids were immobilized in the presence of fluid drag. Upon stepwise reduction in eluting solution ionic strength (IS), a sharp release of colloids occurred in each step which indicates that colloid retention depends on a balance of applied (hydrodynamic) and resisting (adhesive) torques which varied with pore space geometry, surface roughness, and interaction energy. When the eluting fluid IS was reduced to deionized water, the final retention locations occurred near grain-grain contacts, and colloid aggregation was sometimes observed in micromodel experiments. Significant amounts of colloid retention hysteresis with IS were observed in the column experiments, and it depended on the porous medium (glass beads compared with sand), the colloid size (1.1 and 0.5 mum), and on the initial deposition IS. These observations were attributed to weak adhesive interactions that depended on the double layer thickness (e.g., the depth of the secondary minimum and/or nanoscale heterogeneity), colloid mass transfer on the solid phase to regions where the torque and force balances were favorable for retention, the number and extent of grain-grain contacts, and surface roughness.

Modeling Microorganism Transport and Survival in the Subsurface

An understanding of microbial transport and survival in the subsurface is needed for public health, environmental applications, and industrial processes. Much research has therefore been directed to quantify mechanisms influencing microbial fate, and the results demonstrate a complex coupling among many physical, chemical, and biological factors. Mathematical models can be used to help understand and predict the complexities of microbial transport and survival in the subsurface under given assumptions and conditions. This review highlights existing model formulations that can be used for this purpose. In particular, we discuss models based on the advection-dispersion equation, with terms for kinetic retention to solid-water and/or air-water interfaces; blocking and ripening; release that is dependent on the resident time, diffusion, and transients in solution chemistry, water velocity, and water saturation; and microbial decay (first-order and Weibull) and growth (logistic and Monod) that is dependent on temperature, nutrient concentration, and/or microbial concentration. We highlight a two-region model to account for microbe migration in the vicinity of a solid phase and use it to simulate the coupled transport and survival of species under a variety of environmentally relevant scenarios. This review identifies challenges and limitations of models to describe and predict microbial transport and survival. In particular, many model parameters have to be optimized to simulate a diversity of observed transport, retention, and survival behavior at the laboratory scale. Improved theory and models are needed to predict the fate of microorganisms in natural subsurface systems that are highly dynamic and heterogeneous.

Impacts of bridging complexation on the transport of surface-modified nanoparticles in saturated sand

The transport of polyacrylic acid capped cadmium telluride (CdTe) quantum dots (QDs), carboxylate-modified latex (CML), and bare silica nanoparticles (NPs) was studied in packed columns at various electrolyte concentrations and cation types. The breakthrough curves (BTCs) of QDs and CML particles in acid-treated Accusand showed significant amounts of increasing deposition with 0.5, 1, and 2 mM Ca(2+), but only minute deposition at 50 and 100 mM Na(+). Negligible QD and CML deposition occurred at 2mM Ca(2+) in columns packed with ultrapure quartz sand that was similar in size to the Accusand. These observations are not consistent with interpretations based on Derjaguin-Landau-Verwey-Overbeek (DLVO) calculations of interaction energies. Scanning electron microscopy (SEM) and energy-dispersive X-ray (EDX) analysis demonstrated that there were regions on the acid-treated Accusand covered with small amounts of clay that were absent on the ultrapure quartz sand. A salt cleaning method was therefore used to remove the clay from the acid-treated Accusand. The BTCs of QDs and CML in this acid+salt treated Accusand exhibited much less deposition at any given Ca(2+) concentration compared to those obtained from the acid-treated sand. SEM images showed that most of the QD deposited in acid-treated Accusand occurred on clay surfaces. Unlike our results with QDs and CML, negligible deposition of bare silica NPs occurred at 5 and 10 mM Ca(2+) in acid-treated Accusand. The high deposition of QDs and CML particles was therefore attributed to bridging complexation in which Ca(2+) serves as a bridge between the cation exchange locations on the clay and carboxyl functional groups on the QD and CML particles, which were absent on the bare silica NPs. Our results suggest that the transport of carboxylic ligand-modified NPs may be limited in subsurface environments because of the ubiquitous presence of clay and divalent cations.

Determining Parameters and Mechanisms of Colloid Retention and Release in Porous Media

A modeling framework is presented to determine fundamental parameters and controlling mechanisms of colloid (microbes, clays, and nanoparticles) retention and release on surfaces of porous media that exhibit wide distributions of nanoscale chemical heterogeneity, nano- to microscale roughness, and pore water velocity. Primary and/or secondary minimum interactions in the zone of electrostatic influence were determined over the heterogeneous solid surface. The Maxwellian kinetic energy model was subsequently employed to determine the probability of immobilization and diffusive release of colloids from each of these minima. In addition, a balance of applied hydrodynamic and resisting adhesive torques was conducted to determine locations of immobilization and hydrodynamic release in the presence of spatially variable water flow and microscopic roughness. Locations for retention had to satisfy both energy and torque balance conditions for immobilization, whereas release could occur either due to diffusion or hydrodynamics. Summation of energy and torque balance results over the elementary surface area of the porous medium provided estimates for colloid retention and release parameters that are critical to predicting environmental fate, including the sticking and release efficiencies and the maximum concentration of retained colloids on the solid phase. Nanoscale roughness and chemical heterogeneity produced localized primary minimum interactions that controlled long-term retention, even when mean chemical conditions were unfavorable. Microscopic roughness played a dominant role in colloid retention under low ionic strength and high hydrodynamic conditions, especially for larger colloids.

Release of Quantum Dot Nanoparticles in Porous Media: Role of Cation Exchange and Aging Time

Understanding the fate and transport of engineered nanoparticles (ENPs) in subsurface environments is required for developing the best strategy for waste management and disposal of these materials. In this study, the deposition and release of quantum dot (QD) nanoparticles were studied in saturated sand columns. The QDs were first deposited in columns using 100 mM NaCl or 2 mM CaC12 solutions. Deposited QDs were then contacted with deionized (DI) water and/or varying Na(+) concentrations to induce release. QDs deposited in 100 mM Na(+) were easily reversible when the column was rinsed with DI water. Conversely, QDs deposited in the presence of Ca(2+) exhibited resistance to release with DI water. However, significant release occurred when the columns were flushed with NaCl solutions. This release behavior was explained by cation exchange (Ca(2+) in exchange sites were replaced by Na(+)) which resulted in the breakdown of calcium bridging. We also studied the effect of aging time on the QD release. As the aging time increased, smaller amounts of QDs were released following cation exchange. However, deposited QDs were subsequently released when the column was flushed with DI water. The release behavior was modeled using a single first-order kinetic release process and changes in the maximum solid phase concentration of deposited QDs with transition in solution chemistry. The results of this study demonstrate that the presence of carboxyl groups on ENPs and divalent ions in the solution plays a key role in controlling ENP mobility in the subsurface environment.

An explanation for differences in the process of colloid adsorption in batch and column studies

It is essential to understand the mechanisms that control virus and bacteria removal in the subsurface environment to assess the risk of groundwater contamination with fecal microorganisms. This study was conducted to explicitly provide a critical and systematic comparison between batch and column experiments. The aim was to investigate the underlying factors causing the commonly observed discrepancies in colloid adsorption process in column and batch systems. We examined the colloid adsorption behavior of four different sizes of carboxylate-modified latex (CML) microspheres, as surrogates for viruses and bacteria, on quartz sand in batch and column experiments over a wide range of solution ionic strengths (IS). Our results show that adsorption of colloids in batch systems should be considered as an irreversible attachment because the attachment/detachment model was found to be inadequate in describing the batch results. An irreversible attachment-blocking model was found to accurately describe the results of both batch and column experiments. The rate of attachment was found to depend highly on colloid size, solution IS and the fraction of the sand surface area favorable for attachment (Sf). The rate of attachment and Sf values were different in batch and column experiments due to differences in the hydrodynamic of the system, and the role of surface roughness and pore structure on colloid attachment. Results from column and batch experiments were generally not comparable, especially for larger colloids (≥0.5μm). Predictions based on classical DLVO theory were found to inadequately describe interaction energies between colloids and sand surfaces.

Colloid release and clogging in porous media: Effects of solution ionic strength and flow velocity

The release and retention of in-situ colloids in aquifers play an important role in the sustainable operation of managed aquifer recharge (MAR) schemes. The processes of colloid release, retention, and associated permeability changes in consolidated aquifer sediments were studied by displacing native groundwater with reverse osmosis-treated (RO) water at various flow velocities. Significant amounts of colloid release occurred when: (i) the native groundwater was displaced by RO-water with a low ionic strength (IS), and (ii) the flow velocity was increased in a stepwise manner. The amount of colloid release and associated permeability reduction upon RO-water injection depended on the initial clay content of the core. The concentration of released colloids was relatively low and the permeability reduction was negligible for the core sample with a low clay content of about 1.3%. In contrast, core samples with about 6 and 7.5% clay content exhibited: (i) close to two orders of magnitude increase in effluent colloid concentration and (ii) more than 65% permeability reduction. Incremental improvement in the core permeability was achieved when the flow velocity increased, whereas a short flow interruption provided a considerable increase in the core permeability. This dependence of colloid release and permeability changes on flow velocity and colloid concentration was consistent with colloid retention and release at pore constrictions due to the mechanism of hydrodynamic bridging. A mathematical model was formulated to describe the processes of colloid release, transport, retention at pore constrictions, and subsequent permeability changes. Our experimental and modeling results indicated that only a small fraction of the in-situ colloids was released for any given change in the IS or flow velocity. Comparison of the fitted and experimentally measured effluent colloid concentrations and associated changes in the core permeability showed good agreement, indicating that the essential physics were accurately captured by the model.

Equilibrium and kinetic models for colloid release under transient solution chemistry conditions.

We present continuum models to describe colloid release in the subsurface during transient physicochemical conditions. Our modeling approach relates the amount of colloid release to changes in the fraction of the solid surface area that contributes to retention. Equilibrium, kinetic, equilibrium and kinetic, and two-site kinetic models were developed to describe various rates of colloid release. These models were subsequently applied to experimental colloid release datasets to investigate the influence of variations in ionic strength (IS), pH, cation exchange, colloid size, and water velocity on release. Various combinations of equilibrium and/or kinetic release models were needed to describe the experimental data depending on the transient conditions and colloid type. Release of Escherichia coli D21g was promoted by a decrease in solution IS and an increase in pH, similar to expected trends for a reduction in the secondary minimum and nanoscale chemical heterogeneity. The retention and release of 20nm carboxyl modified latex nanoparticles (NPs) were demonstrated to be more sensitive to the presence of Ca(2+) than D21g. Specifically, retention of NPs was greater than D21g in the presence of 2mM CaCl2 solution, and release of NPs only occurred after exchange of Ca(2+) by Na(+) and then a reduction in the solution IS. These findings highlight the limitations of conventional interaction energy calculations to describe colloid retention and release, and point to the need to consider other interactions (e.g., Born, steric, and/or hydration forces) and/or nanoscale heterogeneity. Temporal changes in the water velocity did not have a large influence on the release of D21g for the examined conditions. This insensitivity was likely due to factors that reduce the applied hydrodynamic torque and/or increase the resisting adhesive torque; e.g., macroscopic roughness and grain-grain contacts. Our analysis and models improve our understanding and ability to describe the amounts and rates of colloid release and indicate that episodic colloid transport is expected under transient physicochemical conditions.