Modeling coupled nanoparticle aggregation and transport in porous media: a Lagrangian approach

Changes in nanoparticle size and shape due to particle-particle interactions (i.e., aggregation or agglomeration) may significantly alter particle mobility and retention in porous media. To date, however, few modeling studies have considered the coupling of transport and particle aggregation processes. The majority of particle transport models employ an Eulerian modeling framework and are, consequently, limited in the types of collisions and aggregate sizes that can be considered. In this work, a more general Lagrangian modeling framework is developed and implemented to explore coupled nanoparticle aggregation and transport processes. The model was verified through comparison of model simulations to published results of an experimental and Eulerian modeling study (Raychoudhury et al., 2012) of carboxymethyl cellulose (CMC)-modified nano-sized zero-valent iron particle (nZVI) transport and retention in water-saturated sand columns. A model sensitivity analysis reveals the influence of influent particle concentration (ca. 70 to 700 mg/L), primary particle size (10-100 nm) and pore water velocity (ca. 1-6 m/day) on particle-particle, and, consequently, particle-collector interactions. Model simulations demonstrate that, when environmental conditions promote particle-particle interactions, neglecting aggregation effects can lead to under- or over-estimation of nanoparticle mobility. Results also suggest that the extent to which higher order particle-particle collisions influence aggregation kinetics will increase with the fraction of primary particles. This work demonstrates the potential importance of time-dependent aggregation processes on nanoparticle mobility and provides a numerical model capable of capturing/describing these interactions in water-saturated porous media.