Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.     

A novel two-dimensional model for colloid transport in physically and geochemically heterogeneous porous media

A two-dimensional model for colloid transport in geochemically and physically heterogeneous porous media is presented. The model considers patchwise geochemical heterogeneity, which is suitable to describe the chemical variability of many surficial aquifers with ferric oxyhydroxide-coated porous matrix, as well as spatial variability of hydraulic conductivity, which results in heterogeneous flow field. The model is comprised of a transient fluid flow equation, a transient colloid transport equation, and an equation for the dynamics of colloid deposition and release. Numerical simulations were carried out with the model to investigate the colloid transport behavior in layered and randomly heterogeneous porous media. Results demonstrate that physical and geochemical heterogeneities markedly affect the colloid transport behavior. Layered physical or geochemical heterogeneity can result in distinct preferential flow paths of colloidal particles. Furthermore, the combined effect of layered physical and geochemical heterogeneity may result in enhanced or reduced preferential flow of colloids. Random distribution of physical heterogeneity (hydraulic conductivity) results in a random flow field and an irregularly distributed colloid concentration profile in the porous medium. Contrary to random physical heterogeneity, the effect of random patchwise geochemical heterogeneity on colloid transport behavior is not significant. It is mostly the mean value of geochemical heterogeneity rather than its distribution that governs the colloid transport behavior.     

Transport of a non-volatile contaminant in unsaturated porous media in the presence of colloids

Stable colloidal particles can travel long distances in subsurface environments and carry particle-reactive contaminants with them to locations further than predicted by the conventional advective-dispersive transport equation. When such carriers exist in a saturated porous medium, the system can be idealized as consisting of three phases: an aqueous phase, a carrier phase, and a stationary solid matrix phase. However, when colloids are present in an unsaturated porous medium, the system representation should include one more phase, i.e. the air phase. In the work reported, a mathematical model was developed to describe the transport and fate of the colloidal particles and a non-volatile contaminant in unsaturated porous media. The model is based on mass balance equations in a four-phase porous medium. Colloid mass transfer mechanisms among aqueous, solid matrix, and air phases, and contaminant mass transfer between aqueous and colloid phases are represented by kinetic expressions. Governing equations are non-dimensionalized and solved to investigate colloid and contaminant transport in an unsaturated porous medium. A sensitivity analysis of the transport model was utilized to assess the effects of several parameters on model behavior. The colloid transport model matches successfully with experimental data of Wan and Wilson. The presence of air-water interface retards the colloid transport significantly counterbalancing the facilitating effect of colloids. However, the retardation of contaminant transport by colloids is highly dependent on the properties of the contaminant and the colloidal surface.     

Virus retention and transport through Al-oxide coated sand columns: effects of ionic strength and composition

Knowledge of the factors that influence the fate and transport of viruses in porous media is very important for accurately determining groundwater vulnerability and for developing protective regulations. In this study, six saturated sand column experiments were performed to examine the effects of a positively charged Al-oxide, which was coated on sand particles, on the retention and transport of viruses (phiX174 and MS-2) in background solutions of different ionic strength and composition. We found that the Al-oxide coating on sand significantly removed viruses during their transport in a phosphate buffered saline (PBS) solution. Mass balance calculations showed that 34% of the input MS-2 was inactivated/irreversibly sorbed on the surface of Al-oxide coated sand whereas 100% of phiX174 was recovered. Results from this study also indicated that higher ionic strength facilitated the transport of both phiX174 and MS-2 through the Al-oxide coated sand. This was attributed to the effect of ion shielding, which at higher ionic strength decreased the electrostatic attraction between the viral particles and the sand surface and consequently decreased virus sorption. Strong effect of the ionic strength indicates that an outer-sphere complexation mechanism was responsible for the virus sorption on the Al-oxide coated sand. Ion composition of the background solutions was also found to be a significant factor in influencing virus retention and transport. Virus transport was enhanced in the presence of phosphate (HPO(4)(2-)) as compared to bicarbonate (HCO(3)(-)), and the effect of HPO(4)(2-) was more significant on MS-2 than on phiX174. The presence of bivalent cations (Ca(2+) and Mg(2+)) increased virus transport because the cations partially screened the negative charges on the viruses therefore decreased the electrostatic attraction between the positively charged sand surface and the negatively charged viruses. Mass recovery data indicated that bivalent cations gave rise to a certain degree of inactivation/irreversibly sorption of phiX174 on the surface of Al-oxide coated sand. On the contrary, the bivalent cations appeared to have protected MS-2 from inactivation/irreversibly sorption. This study provides some insights into the mechanisms responsible for virus retention and transport in porous media.     

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.     

Contaminant transport in dual-porosity media with dissolved organic matter and bacteria present as mobile colloids

In riverbank filtration, contaminant transport is affected by colloidal particles such as dissolved organic matter (DOM) and bacterial particles. In addition, the subsurface heterogeneity influences the behavior of contaminant transport in riverbank filtration. A mathematical model is developed to describe the contaminant transport in dual-porosity media in the presence of DOM and bacteria as mobile colloids. In the model development, a porous medium is divided into the mobile and immobile regions to consider the presence of ineffective micropores in physically heterogeneous riverbanks. We assume that the contaminant transport in the mobile region is controlled by the advection and dispersion while the contaminant transport in the immobile region occurs due to the molecular diffusion. The contaminant transfer between the mobile and immobile regions takes place by diffusive mass transfer. The mobile region is conceptualized as a four-phase system: two mobile colloidal phases, an aqueous phase, and a solid matrix. The complete set of governing equations is solved numerically with a fully implicit finite difference method. The model results show that in riverbank filtration, the contaminant can migrate further than expected due to the presence of DOM and bacteria. In addition, the contaminant mobility increases further in the presence of the immobile region in aquifers. A sensitivity analysis shows that in dual-porosity media, earlier breakthrough of the contaminant takes place as the volumetric fraction of the mobile region decreases. It is also demonstrated that as the contaminant mass transfer rate coefficient between the mobile and immobile regions increases, the contaminant concentration gradient between the two regions reverses at earlier pore volumes. The contaminant mass transfer coefficient between the mobile and immobile regions mainly controls the tailing effect of the contaminant breakthrough. The contaminant breakthrough curves are sensitive to changes in contaminant adsorption and desorption rate coefficients on DOM and bacteria. In situations where the contaminant is released in the presence of DOM and bacteria in dual-porosity media, the early breakthrough and tailing occur due to the colloidal facilitation and presence of immobile regions.     

Modeling colloid transport for performance assessment

The natural system is expected to contribute to isolation at the proposed high-level nuclear waste (HLW) geologic repository at Yucca Mountain, NV (YM). In developing performance assessment (PA) computer models to simulate long-term behavior at YM, colloidal transport of radionuclides has been proposed as a critical factor because of the possible reduced interaction with the geologic media. Site-specific information on the chemistry and natural colloid concentration of saturated zone groundwaters in the vicinity of YM is combined with a surface complexation sorption model to evaluate the impact of natural colloids on calculated retardation factors (RF) for several radioelements of concern in PA. Inclusion of colloids into the conceptual model can reduce the calculated effective retardation significantly. Strongly sorbed radionuclides such as americium and thorium are most affected by pseudocolloid formation and transport, with a potential reduction in RF of several orders of magnitude. Radioelements that are less strongly sorbed under YM conditions, such as uranium and neptunium, are not affected significantly by colloid transport, and transport of plutonium in the valence state is only moderately enhanced. Model results showed no increase in the peak mean annual total effective dose equivalent (TEDE) within a compliance period of 10,000 years, although this is strongly dependent on container life in the base case scenario. At longer times, simulated container failures increase and the TEDE from the colloidal models increased by a factor of 60 from the base case. By using mechanistic models and sensitivity analyses to determine what parameters and transport processes affect the TEDE, colloidal transport in future versions of the TPA code can be represented more accurately.     

Determination of spatially-resolved porosity, tracer distributions and diffusion coefficients in porous media using MRI measurements and numerical simulations

This work is focused on measuring the concentration distribution of a conservative tracer in a homogeneous synthetic porous material and in heterogeneous natural sandstone using MRI techniques, and on the use of spatially resolved porosity data to define spatially variable diffusion coefficients in heterogeneous media. The measurements are made by employing SPRITE, a fast MRI method that yields quantitative, spatially-resolved tracer concentrations in porous media. Diffusion experiments involving the migration of H(2)O into D(2)O-saturated porous media are conducted. One-dimensional spatial distributions of H(2)O-tracer concentrations acquired from experiments with the homogeneous synthetic calcium silicate are fitted with the one-dimensional analytical solution of Fick's second law to confirm that the experimental method provides results that are consistent with expectations for Fickian diffusion in porous media. The MRI-measured concentration profiles match well with the solution for Fick's second law and provide a pore-water diffusion coefficient of 1.75×10(-9)m(2)s(-1). The experimental approach was then extended to evaluate diffusion in a heterogeneous natural sandstone in three dimensions. The relatively high hydraulic conductivity of the sandstone, and the contrast in fluid density between the H(2)O tracer and the D(2)O pore fluid, lead to solute transport by a combination of diffusion and density-driven advection. The MRI measurements of spatially distributed tracer concentration, combined with numerical simulations allow for the identification of the respective influences of advection and diffusion. The experimental data are interpreted with the aid of MIN3P-D - a multicomponent reactive transport code that includes the coupled processes of diffusion and density-driven advection. The model defines local diffusion coefficients as a function of spatially resolved porosity measurements. The D(e) values calculated for the heterogeneous sandstone and used to simulate diffusive and advective transport range from 5.4×10(-12) to 1.0×10(-10)m(2)s(-1). These methods have broad applicability to studies of contaminant migration in geological materials.     

Influence of soluble copper on the electrokinetic properties and transport of copper oxychloride-based fungicide particles

This article describes the influence of dissolved copper on the electrokinetic properties and transport of a copper oxychloride-based fungicide (COF) in porous media. The Zeta potential (ζ) of COF particles increases (viz. becomes less negative) with increasing concentration of Cu(2+) in the bulk solution. ζ decreases for COF when the electrolyte (NaNO(3)) concentration is raised from 1 to 10mM. This can be ascribed to ion correlation of Cu(2+) in the electrical double layer (EDL). COF transport tests in quartz sand columns showed the addition of Cu(2+) to the bulk solution to result in increased retention of the metal. Modelling particle deposition dynamics provided results consistent with kinetic attachment. Based on the effect of soluble Cu on colloid mobility, the transport of particulate and soluble forms of copper is coupled via the chemistry of pore water and colloid interactions. Mutual effects between cations and colloids should thus be considered when determining the environmental fate of particulate and soluble forms of copper in soil and groundwater (especially at copper-contaminated sites).     

The concentration, apparent molecular weight and chemical reactivity of silica from groundwater in southern Nevada

Sorption of radionuclides, metals and organic compounds to colloidal particles has been suggested to increase the mobility of these pollutants in groundwater. Because silicates and alumino-silicates can be important components of groundwater colloids, we have conducted a study to characterize the nature of silica in various springs and wells in Southern Nevada and to determine the extent that silica may be associated with colloidal particles that can participate in pollutant transport. The total silica content was measured using inductively coupled plasma emission spectroscopy (ICP). In addition, reactive silica was measured using the silica molybdate colorimetric technique. The apparent molecular weight of the silica was investigated using split-flow-lateral-transport-thin-cell-fractionation (SPLITT) which can readily distinguish between colloidal and low molecular weight associations. This study indicates that silica does not tend to form stable inorganic colloids in Southern Nevada groundwaters but exists as low molecular weight species. However, water from one of the test facilities on the Nevada Test Site (NTS) did contain stable siliceous colloids that could have important implications for the modeling the transport of radionuclides at this site.     

Gel barrier formation in unsaturated porous media

The gel barrier formation by a gelling liquid (Colloidal Silica) injection in an unsaturated porous medium is investigated by developing a mathematical model and conducting numerical simulations. Gelation process is initiated by adding electrolytes such as NaCl, and the gel phase consisting of cross-linked colloidal silica particles grows as the gelation process proceeds. The mathematical model describing the transport and gelation of Colloidal Silica (CS) is based on coupled mass balance equations for the gel mixture (the sol phase plus the gel phase), gel phase (cross-linked colloidal silica particles plus water captured between cross-linked particles), and colloidal silica particles (discrete and cross-linked) and NaCl in the sol (suspension of discrete colloidal silica particles in water) and gel phases. The solutions in terms of volumetric fraction of the gel phase yield the gel mixture viscosity via the dependency on the volumetric fraction of gel phase. This dependency is determined from a kinetic gelation model with time-normalized viscosity curves. The proposed model is verified by comparing experimentally and numerically determined hydraulic conductivities of gel-treated soil columns at different CS injection volumes. The numerical experiments indicate that an impermeable gel layer is formed within the time period twice the gel-point in a one-dimensional flow system. At the same normalized time corresponding to twice the gel-point, the CS solutions with lower NaCl concentrations result in further migration and poor performance in plugging the pore space. The viscosity computation proposed in this study is compared with another method available in the literature. It is observed that the other method estimates the viscosity at the mixing zone higher than the one proposed by the authors. The proposed model can simulate realistic injection scenarios with various combinations of operating parameters such as NaCl concentration and NaCl mixing time, and thus providing guidelines in performing this technology on site.     

Particle transport in unsaturated fractured chalk under arid conditions

A series of field and laboratory experiments were conducted to study the mechanisms of particle detachment and transport from fractures in vadose chalk. Experiments of intermittent flow events along fracture surfaces were carried out in the laboratory. In the field, water was percolated from land surface via a discrete fracture into a compartmental sampler installed inside a horizontal corehole located I m below the surface. The mass, size distribution, and composition of the particles drained from the fracture voids were examined along with flow rates and salt dissolution. Two boreholes penetrating the underlying saturated zone were sampled and analyzed for colloidal concentration and composition. Most of the particle and solute release at the drained effluents occurred during the first several hours of flow, but erratic pulses of particles were still observed after long periods of time. Most of the detached particles had a mean diameter of >2 microm, while the mobile colloidal phase in the groundwater had a mean diameter of approximately 1 microm. Mineralogical composition of the groundwater colloids and the particles detached from the upper vadose fracture were similar. Laboratory observations demonstrated the importance of the existence of a coating layer, made of weathered particles and salts, on particle detachment. The results of this study suggest that: (1) particle detachment causes flow-rate variability in the unsaturated fracture; (2) the mechanisms of particle detachment and salt dissolution within the fracture are linked: and (3) although most of the detached particles are large and likely to accumulate inside fractures, some colloidal particles also eroded from the fracture void and are likely to be transported to the groundwater.     

Suspended particulate matter determines physical speciation of Fe,Mn, and trace metals in surface waters of Loire watershed

This study investigates the spatiotemporal variability of major and trace elements, dissolved organic carbon (DOC), total dissolved solids (TDS), and suspended particulate matter (SPM) in surface waters of several hydrosystems of the Loire River watershed in France. In particular, this study aims to delineate the impact of the abovementioned water physicochemical parameters on natural iron and manganese physical speciation (homoaggregation/heteroaggregation) among fine colloidal and dissolved (<?10 nm), colloidal (10–450 nm) and particulate (>?450 nm) phases in Loire River watershed. Results show that the chemistry of the Loire River watershed is controlled by two end members: magmatic and metamorphic petrographic context on the upper part of the watershed; and sedimentary rocks for the middle and low part of the Loire. The percentage of particulate Fe and Mn increased downstream concurrent with the increase in SPM and major cations concentration, whereas the percentage of colloidal Fe and Mn decreased downstream. Transmission electron microscopy analyses of the colloidal and particulate fractions (from the non-filtered water sample) revealed that heteroaggregation of Fe and Mn rich natural nanoparticles and natural organic matter to the particulate phase is the dominant mechanism. The heteroaggregation controls the partitioning of Fe and Mn in the different fractions, potentially due to the increase in the ionic strength, and divalent cations concentration downstream, and SPM concentration. These findings imply that SPM concentration plays an important role in controlling the fate and behavior of Fe and Mn in various sized fractions.

Physical speciation by heteroaggregation of (Fe-Mn) compounds: high [SPM]?→?[Fe-Mn] particulate faction; low {SPM]?→?[Fe-Mn] colloid-dissolved fraction.

     

Straining phenomena in bacteria transport through natural porous media

Background, aim, and scope  

Transport of bacteria through natural porous media is an issue of increasing concern arising in several very important environmental processes. These include the percolation of bacteria from fecal waste to drinking water reservoirs, thus leading to a risk for human health, or the bioremediation of contaminated soils in which the bacteria are expected to travel long distances underground in order to reach contaminated areas and degrade chemicals originating from accidental spills. An understanding of bacterial retention and transport mechanisms in porous media would be of great help in the development of models able to predict the distance covered by bacterial suspensions in these situations.     

Influences of dissolved and colloidal organic carbon on the uptake of Ag, Cd, and Cr by the marine mussel Perna viridis

The cross-flow ultrafiltration and radiotracer techniques were used to study the influences of natural dissolved organic carbon (DOC) and colloidal organic carbon (COC) on the bioavailability of Ag, Cd, and Cr to the green mussel Perna viridis. We examined the uptake of these metals by the mussels at different concentrations of DOC and COC from different origins (estuarine, coastal, and diatom decomposed). Using the DOC originating from the decomposed diatom (Thalassiosira pseudonana), we demonstrated that Cd and Cr uptake, quantified by the concentration factor (DCF), increased linearly with increasing DOC concentration. There was, however, no consistent influence of natural DOC concentration on the metal uptake when the DOC was obtained from different sources of seawater (coastal and estuarine). The influences of COC on metal bioavailability were metal-specific and dependent on the geochemical properties of colloids and colloid-metal complexation. Cd uptake rate was not influenced by the COC concentrations. Uptake of diatom-decomposed colloidal Cr was enhanced by 3.4x, whereas the uptake of diatom-decomposed colloidal Ag was decreased by 8.2x compared with the uptake of low molecular weight Cr and Ag (<1 kDa). The uptake of diatom-decomposed colloidal Cr and Ag was generally lower than the uptake of metals bound with the same type of colloids for 2 days. Further aging of the colloid-metal binding reduced metal bioavailability to the mussels. In the presence of different sizes of colloidal particles where there was no major binding of colloids with the metals, metal uptake by the mussels was not influenced by different COC concentrations. Overall, our study suggests that although metal dissociation from colloids may be an important step for the uptake of colloidal metals, other mechanisms such as pinocytosis and co-transport may also be involved in the uptake of these metals, especially in aquatic environments with high DOC and COC concentrations.     

Tracer diffusion coefficients in sedimentary rocks: correlation to porosity and hydraulic conductivity

Matrix diffusion is an important transport process in geologic materials of low hydraulic conductivity. For predicting the fate and transport of contaminants, a detailed understanding of the diffusion processes in natural porous media is essential. In this study, diffusive tracer transport (iodide) was investigated in a variety of geologically different limestone and sandstone rocks. Porosity, structural and mineralogical composition, hydraulic conductivity, and other rock properties were determined. The effective diffusion coefficients were measured using the time-lag method. The results of the diffusion experiments indicate that there is a close relationship between total porosity and the effective diffusion coefficient of a rock (analogous to Archie's Law). Consequently, the tortousity factor can be expressed as a function of total porosity. The relationship fits best for thicker samples (> 1.0 cm) with high porosities (> 20%), because of the reduced influence of heterogeneity in larger samples. In general, these correlations appear to be a simple way to determine tortuosity and the effective diffusion coefficient from easy to determine rock porosity values.     

Estimating the relevance of engineered carbonaceous nanoparticle facilitated transport of hydrophobic organic contaminants in porous media

Naturally occurring nanoparticles (NP) enhance the transport of hydrophobic organic contaminants (HOCs) in porous media. In addition, the debate on the environmental impact of engineered nanoparticles (ENP) has become increasingly important. HOC bind strongly to carbonaceous ENP. Thus, carbonaceous ENP may also act as carriers for contaminant transport and might be important when compared to existing transport processes. ENP bound transport is strongly linked to the sorption behavior, and other carbonaceous ENP-specific properties. In our analysis the HOC-ENP sorption mechanism, as well as ENP size and ENP residence time, was of major importance. Our results show that depending on ENP size, sorption kinetics and residence time in the system, the ENP bound transport can be estimated either as (1) negligible, (2) enhancing contaminant transport, or (3) should be assessed by reactive transport modeling. One major challenge to this field is the current lack of data for HOC-ENP desorption kinetics.     

Virus transport in physically and geochemically heterogeneous subsurface porous media

A two-dimensional model for virus transport in physically and geochemically heterogeneous subsurface porous media is presented. The model involves solution of the advection-dispersion equation, which additionally considers virus inactivation in the solution, as well as virus removal at the solid matrix surface due to attachment (deposition), release, and inactivation. Two surface inactivation models for the fate of attached inactive viruses and their subsequent role on virus attachment and release were considered. Geochemical heterogeneity, portrayed as patches of positively charged metal oxyhydroxide coatings on collector grain surfaces, and physical heterogeneity, portrayed as spatial variability of hydraulic conductivity, were incorporated in the model. Both layered and randomly (log-normally) distributed physical and geochemical heterogeneities were considered. The upstream weighted multiple cell balance method was employed to numerically solve the governing equations of groundwater flow and virus transport. Model predictions show that the presence of subsurface layered geochemical and physical heterogeneity results in preferential flow paths and thus significantly affect virus mobility. Random distributions of physical and geochemical heterogeneity have also notable influence on the virus transport behavior. While the solution inactivation rate was found to significantly influence the virus transport behavior, surface inactivation under realistic field conditions has probably a negligible influence on the overall virus transport. It was further demonstrated that large virus release rates result in extended periods of virus breakthrough over significant distances downstream from the injection sites. This behavior suggests that simpler models that account for virus adsorption through a retardation factor may yield a misleading assessment of virus transport in "hydrogeologically sensitive" subsurface environments.     

Pore-scale modeling of dispersion in disordered porous media

We employ a direct pore-level model of incompressible flow that uses the modified moving particle semi-implicit (MMPS) method. The model is capable of simulating both unsteady- and steady-state flow directly in microtomography images of naturally-occurring porous media. We further develop this model to simulate solute transport in disordered porous media. The governing equations of flow and transport at the pore level, i.e., Navier-Stokes and convection-diffusion, are solved directly in the pore space mapped by microtomography techniques. Three naturally-occurring sandstones are studied in this work. We verify the accuracy of the model by comparing the computed longitudinal dispersion coefficients against the experimental data for a wide range of Peclet numbers, i.e., 5×10(-2)相似文献