Transport of Neutrally Buoyant and Dense Variably Sized Colloids in a Two-Dimensional Fracture with Anisotropic Aperture

The transport of monodisperse as well as polydisperse colloid suspensions in a two-dimensional, water saturated fracture with spatially variable and anisotropic aperture is investigated with a particle tracking model. Both neutrally buoyant and dense colloid suspensions are considered. Although flow and transport in fractured subsurface formations have been studied extensively by numerous investigators, the transport of dense, polydisperse colloid suspensions in a fracture with spatially variable and anisotropic aperture has not been previously explored. Simulated snapshots and breakthrough curves of ensemble averages of several realizations of a log-normally distributed aperture field show that polydisperse colloids exhibit greater spreading than monodisperse colloids, and dense colloids show greater retardation than neutrally buoyant colloids. Moreover, it is demonstrated that aperture anisotropy oriented along the flow direction substantially increases colloid spreading; whereas, aperture anisotropy oriented transverse to the flow direction retards colloid movement.     

Three-Dimensional Analytical Models for Virus Transport in Saturated Porous Media

Analytical models for virus transport in saturated, homogeneous porous media are developed. The models account for three-dimensional dispersion in a uniform flow field, and first-order inactivation of suspended and deposited viruses with different inactivation rate coefficients. Virus deposition onto solid particles is described by two different processes: nonequilibrium adsorption which is applicable to viruses behaving as solutes; and colloid filtration which is applicable to viruses behaving as colloids. The governing virus transport equations are solved analytically by employing Laplace/Fourier transform techniques. Instantaneous and continuous/periodic virus loadings from a point source are examined.     

Transport of Pseudomonas putida in a 3-D Bench Scale Experimental Aquifer

This study is focused on the transport of Pseudomonas (P.) putida bacterial cells in a 3-D model aquifer. The pilot-scale aquifer consisted of a rectangular glass tank with internal dimensions: 120?cm length, 48?cm width, and 50?cm height, carefully packed with well-characterized quartz sand. The P. putida decay was adequately represented by a first-order model. Transport experiments with a conservative tracer and P. putida were conducted to characterize the aquifer and to investigate the bacterial behavior during transport in water saturated porous media. A 3-D, finite-difference numerical model for bacterial transport in saturated, homogeneous porous media was developed and was used to successfully fit the experimental data. Furthermore, theoretical interaction energy calculations suggested that the extended-DLVO theory seems to predict bacteria attachment onto the aquifer sand better than the classical DLVO theory.     

Transport of Viruses Through Saturated and Unsaturated Columns Packed with Sand

Laboratory-scale virus transport experiments were conducted in columns packed with sand under saturated and unsaturated conditions. The viruses employed were the male-specific RNA coliphage, MS2, and the Salmonella typhimurium phage, PRD1. The mathematical model developed by Sim and Chrysikopoulos (Water Resour Res 36:173–179, 2000) that accounts for processes responsible for removal of viruses during vertical transport in one-dimensional, unsaturated porous media was used to fit the data collected from the laboratory experiments. The liquid to liquid–solid and liquid to air–liquid interface mass transfer rate coefficients were shown to increase for both bacteriophage as saturation levels were reduced. The experimental results indicate that even for unfavorable attachment conditions within a sand column (e.g., phosphate-buffered saline solution; pH = 7.5; ionic strength = 2 mM), saturation levels can affect virus transport through porous media.     

Non-aqueous phase liquid drop formation within a water saturated fracture

This work focuses on the mechanisms of non-aqueous phase liquid (NAPL) drop formation within a single fracture fed from a NAPL reservoir by way of a circular orifice, such as a pore. The fracture is assumed to be fully saturated, the relative wettability of the system is assumed water-wet, and the water velocity profile within the fracture is described by a Poiseuille flow. The size of the NAPL drops is investigated for various water flow velocities and NAPL entrance diameters. A force balancing method was used to determine the radii of detached drops. The drop sizes calculated from the model developed here are shown to be in agreement with available experimental drop size data. It is shown that at low Reynolds numbers the buoyancy force is the dominant force acting on the drop during the formation process and at high Reynolds numbers the viscous forces dominate. A simplified expression relating the geometry of the fractured system to the drop radii is developed from the model equations, and it is shown to predict drop radii that match well with both the model simulations and the available experimental data.     

Generalized Taylor-Aris moment analysis of the transport of sorbing solutes through porous media with spatially-periodic retardation factor

Taylor-Aris dispersion theory, as generalized by Brenner, is employed to investigate the macroscopic behavior of sorbing solute transport in a three-dimensional, hydraulically homogeneous porous medium under steady, unidirectional flow. The porous medium is considered to possess spatially periodic geochemical characteristics in all three directions, where the spatial periods define a rectangular parallelepiped or a unit-element. The spatially-variable geochemical parameters of the solid matrix are incorporated into the transport equation by a spatially-periodic distribution coefficient and consequently a spatially-periodic retardation factor. Expressions for the effective or large-time coefficients governing the macroscopic solute transport are derived for solute sorbing according to a linear equilibrium isotherm as well as for the case of a first-order kinetic sorption relationship. The results indicate that for the case of a chemical equilibrium sorption isotherm the longitudinal macrodispersion incorporates a second term that accounts for the eflect of averaging the distribution coefficient over the volume of a unit element. Furthermore, for the case of a kinetic sorption relation, the longitudinal macrodispersion expression includes a third term that accounts for the effect of the first-order sorption rate. Therefore, increased solute spreading is expected if the local chemical equilibrium assumption is not valid. The derived expressions of the apparent parameters governing the macroscopic solute transport under local equilibrium conditions agreed reasonably with the results of numerical computations using particle tracking techniques.     

Experimental investigation of acoustically enhanced colloid transport in water-saturated packed columns

The effects of acoustic wave propagation on the transport of colloids in saturated porous media were investigated by injecting Uranine (conservative tracer) as well as blue and red polystyrene microspheres (colloids of different diameters; 0.10 and 0.028 mum, respectively) into a column packed with glass beads. Experiments were conducted by maintaining the acoustic pressure at the influent at 23.0 kPa with acoustic frequencies ranging from 30 to 150 Hz. The experimental results suggested that colloid size did not affect the forward and reverse attachment rate coefficients. The acoustic pressure caused an increase in the effective interstitial velocity at all frequencies for the conservative tracer and colloids of both sizes, with maximum increase at 30 Hz. Furthermore, acoustics enhanced the dispersion process at all frequencies, with a maximum at 30 Hz.     

Effective velocity and effective dispersion coefficient for finite-sized particles flowing in a uniform fracture

In this work we derive expressions for the effective velocity and effective dispersion coefficient for finite-sized spherical particles with neutral buoyancy flowing within a water saturated fracture. We considered the miscible displacement of a fluid initially free of particles by another fluid containing particles of finite size in suspension within a fracture formed by two semi-infinite parallel plates. Particle spreading occurs due to the combined actions of molecular diffusion and the dispersive effect of the Poiseuille velocity profile. Unlike Taylor dispersion, here the finite size of the particles is taken into account. It is shown that because the finite size of a particle excludes it from the slowest moving portion of the velocity profile, the effective particle velocity is increased, while the overall particle dispersion is reduced. A similar derivation applied to particles flowing in uniform tubes yields analogous results. The effective velocity and dispersion coefficient derived in this work for particle transport in fractures with uniform aperture are unique and ideally suited for use in particle tracking models.     

Modeling the Transport of Contaminants Originating from the Dissolution of DNAPL Pools in Aquifers in the Presence of Dissolved Humic Substances

A twodimensional finite difference numerical model was developed to describe the transport of dissolved organics originating from nonaqueous phase liquid (NAPL) pool dissolution in saturated porous media in the presence of dissolved humic substances. A rectangular NAPL pool was considered in a homogeneous porous medium with unidirectional interstitial groundwater velocity. It was assumed that dissolved humic substances and aqueous phase contaminants may sorb onto the solid matrix under local equilibrium conditions. The contaminant in the aqueous phase may undergo firstorder decay. Also, the dissolved contaminant may sorb onto humic substances. The transport properties of dissolved humic substances are assumed to be unaffected by sorbing contaminants, because dissolved humic macromolecules are much larger than dissolved contaminants and sorption of nonpolar contaminants onto humic substances do not affect the overall surface charge of humic substances. The sorption characteristics of dissolved humic substances onto clean sand were determined from column experiments. An effective local mass transfer rate coefficient accounting for the presence of dissolved humic substances was developed. Model simulations indicate that dissolved humic substances substantially increase NAPL pool dissolution, and consequently reduce the required pumpandtreat aquifer remediation time.