Upscaling of nanoparticle transport in porous media under unfavorable conditions: Pore scale to Darcy scale.

Transport and deposition of nanoparticles in porous media is a multi-scale problem governed by several pore-scale processes, and hence, it is critical to link the processes at pore scale to the Darcy-scale behavior. In this study, using pore network modeling, we develop correlation equations for deposition rate coefficients for nanoparticle transport under unfavorable conditions at the Darcy scale based on pore-scale mechanisms. The upscaling tool is a multi-directional pore-network model consisting of an interconnected network of pores with variable connectivities. Correlation equations describing the pore-averaged deposition rate coefficients under unfavorable conditions in a cylindrical pore, developed in our earlier studies, are employed for each pore element. Pore-network simulations are performed for a wide range of parameter values to obtain the breakthrough curves of nanoparticle concentration. The latter is fitted with macroscopic 1-D advection-dispersion equation with a two-site linear reversible deposition accounting for both equilibrium and kinetic sorption. This leads to the estimation of three Darcy-scale deposition coefficients: distribution coefficient, kinetic rate constant, and the fraction of equilibrium sites. The correlation equations for the Darcy-scale deposition coefficients, under unfavorable conditions, are provided as a function of measurable Darcy-scale parameters, including: porosity, mean pore throat radius, mean pore water velocity, nanoparticle radius, ionic strength, dielectric constant, viscosity, temperature, and surface potentials of the particle and grain surfaces. The correlation equations are found to be consistent with the available experimental results, and in qualitative agreement with Colloid Filtration Theory for all parameters, except for the mean pore water velocity and nanoparticle radius.

Dissolution kinetics of volatile organic compound vapors in water: An integrated experimental and computational study.

In this study we performed batch experiments to investigate the dissolution kinetics of trichloroethylene (TCE) and toluene vapors in water at room temperature and atmospheric pressure. The batch systems consisted of a water reservoir and a connected headspace, the latter containing a small glass cylinder filled with pure volatile organic compound (VOC). Results showed that air phase concentrations of both TCE and toluene increased relatively quickly to their maximum values and then became constant. We considered subsequent dissolution into both stirred and unstirred water reservoirs. Results of the stirred experiments showed a quick increase in the VOC concentrations with time up to their solubility limit in water. VOC vapor dissolution was found to be independent of pH. In contrast, salinity had a significant effect on the solubility of TCE and toluene vapors. VOC evaporation and vapor dissolution in the stirred water reservoirs followed first-order rate processes. Observed data could be described well using both simplified analytical solutions, which decoupled the VOC dynamics in the air and water phases, as well as using more complete coupled solutions. However, the estimated evaporation (ke) and dissolution (kd) rate constants differed by up to 70% between the coupled and uncoupled formulations. We also numerically investigated the effects of fluid withdrawal from the small water reservoir due to sampling. While decoupling the VOC air and water phase mass transfer processes produced unreliable estimates of kd, the effects of fluid withdrawal on the estimated rate constants were found to be less important. The unstirred experiments showed a much slower increase in the dissolved VOC concentrations versus time. Molecular diffusion of the VOCs within the aqueous phase became then the limiting factor for mass transfer from air to water. Fluid withdrawal during sampling likely caused some minor convection within the reservoir, which was simulated by increasing the apparent liquid diffusion coefficient.

Modeling the co-transport of viruses and colloids in unsaturated porous media.

A mathematical model is developed to simulate the co-transport of viruses and colloids in unsaturated porous media under steady-state flow conditions. The virus attachment to the mobile and immobile colloids is described using a linear reversible kinetic model. Colloid transport is assumed to be decoupled from virus transport; that is, we assume that colloids are not affected by the presence of attached viruses on their surface. The governing equations are solved numerically using an alternating three-step operator splitting approach. The model is verified by fitting three sets of experimental data published in the literature: (1) Syngouna and Chrysikopoulos (2013) and (2) Walshe et al. (2010), both on the co-transport of viruses and clay colloids under saturated conditions, and (3) Syngouna and Chrysikopoulos (2015) for the co-transport of viruses and clay colloids under unsaturated conditions. We found a good agreement between observed and fitted breakthrough curves (BTCs) under both saturated and unsaturated conditions. Then, the developed model was used to simulate the co-transport of viruses and colloids in porous media under unsaturated conditions, with the aim of understanding the relative importance of various processes on the co-transport of viruses and colloids in unsaturated porous media. The virus retention in porous media in the presence of colloids is greater during unsaturated conditions as compared to the saturated conditions due to: (1) virus attachment to the air-water interface (AWI), and (2) co-deposition of colloids with attached viruses on its surface to the AWI. A sensitivity analysis of the model to various parameters showed that the virus attachment to AWI is the most sensitive parameter affecting the BTCs of both free viruses and total mobile viruses and has a significant effect on all parts of the BTC. The free and the total mobile viruses BTCs are mainly influenced by parameters describing virus attachment to the AWI, virus interaction with mobile and immobile colloids, virus attachment to solid-water interface (SWI), and colloid interaction with SWI and AWI. The virus BTC is relatively insensitive to parameters describing the maximum adsorption capacity of the AWI for colloids, inlet colloid concentration, virus detachment rate coefficient from the SWI, maximum adsorption capacity of the AWI for viruses and inlet virus concentration.

Virus-sized colloid transport in a single pore: model development and sensitivity analysis.

A mathematical model is developed to simulate the transport and deposition of virus-sized colloids in a cylindrical pore throat considering various processes such as advection, diffusion, colloid-collector surface interactions and hydrodynamic wall effects. The pore space is divided into three different regions, namely, bulk, diffusion and potential regions, based on the dominant processes acting in each of these regions. In the bulk region, colloid transport is governed by advection and diffusion whereas in the diffusion region, colloid mobility due to diffusion is retarded by hydrodynamic wall effects. Colloid-collector interaction forces dominate the transport in the potential region where colloid deposition occurs. The governing equations are non-dimensionalized and solved numerically. A sensitivity analysis indicates that the virus-sized colloid transport and deposition is significantly affected by various pore-scale parameters such as the surface potentials on colloid and collector, ionic strength of the solution, flow velocity, pore size and colloid size. The adsorbed concentration and hence, the favorability of the surface for adsorption increases with: (i) decreasing magnitude and ratio of surface potentials on colloid and collector, (ii) increasing ionic strength and (iii) increasing pore radius. The adsorbed concentration increases with increasing Pe, reaching a maximum value at Pe=0.1 and then decreases thereafter. Also, the colloid size significantly affects particle deposition with the adsorbed concentration increasing with increasing particle radius, reaching a maximum value at a particle radius of 100nm and then decreasing with increasing radius. System hydrodynamics is found to have a greater effect on larger particles than on smaller ones. The secondary minimum contribution to particle deposition has been found to increase as the favorability of the surface for adsorption decreases. The sensitivity of the model to a given parameter will be high if the conditions are favorable for adsorption. The results agree qualitatively with the column-scale experimental observations available in the literature. The current model forms the building block in upscaling colloid transport from pore scale to Darcy scale using Pore-Network Modeling.

Oxidation of volatile organic vapours in air by solid potassium permanganate.

Volatile organic compounds (VOCs) may frequently contaminate groundwater and pose threat to human health when migrating into the unsaturated soil zone and upward to the indoor air. The kinetic of chemical oxidation has been investigated widely for dissolved VOCs in the saturated zone. But, so far there have been few studies on the use of in situ chemical oxidation (ISCO) of vapour phase contaminants. In this study, batch experiments were carried out to evaluate the oxidation of trichloroethylene (TCE), ethanol, and toluene vapours by solid potassium permanganate. Results revealed that solid potassium permanganate is able to transform the vapour of these compounds into harmless oxidation products. The degradation rates for TCE and ethanol were higher than for toluene. The degradation process was modelled using a kinetic model, linear in the gas concentration of VOC [ML(-3)] and relative surface area of potassium permanganate grains (surface area of potassium permanganate divided by gas volume) [L(-1)]. The second-order reaction rate constants for TCE, ethanol, and toluene were found to be equal to 2.0×10(-6) cm s(-1), 1.7×10(-7) cm s(-1), and 7.0×10(-8) cm s(-1), respectively.