Continuum-based models and concepts for the transport of nanoparticles in saturated porous media: A state-of-the-science review

Environmental applications of nanoparticles (NP) increasingly result in widespread NP distribution within porous media where they are subject to various concurrent transport mechanisms including irreversible deposition, attachment/detachment (equilibrium or kinetic), agglomeration, physical straining, site-blocking, ripening, and size exclusion. Fundamental research in NP transport is typically conducted at small scale, and theoretical mechanistic modeling of particle transport in porous media faces challenges when considering the simultaneous effects of transport mechanisms. Continuum modeling approaches, in contrast, are scalable across various scales ranging from column experiments to aquifer. They have also been able to successfully describe the simultaneous occurrence of various transport mechanisms of NP in porous media such as blocking/straining or agglomeration/deposition/detachment. However, the diversity of model equations developed by different authors and the lack of effective approaches for their validation present obstacles to the successful robust application of these models for describing or predicting NP transport phenomena. This review aims to describe consistently all the important NP transport mechanisms along with their representative mathematical continuum models as found in the current scientific literature. Detailed characterizations of each transport phenomenon in regards to their manifestation in the column experiment outcomes, i.e., breakthrough curve (BTC) and residual concentration profile (RCP), are presented to facilitate future interpretations of BTCs and RCPs. The review highlights two NP transport mechanisms, agglomeration and size exclusion, which are potentially of great importance in controlling the fate and transport of NP in the subsurface media yet have been widely neglected in many existing modeling studies. A critical limitation of the continuum modeling approach is the number of parameters used upon application to larger scales and when a series of transport mechanisms are involved. We investigate the use of simplifying assumptions, such as the equilibrium assumption, in modeling the attachment/detachment mechanisms within a continuum modelling framework. While acknowledging criticisms about the use of this assumption for NP deposition on a mechanistic (process) basis, we found that its use as a description of dynamic deposition behavior in a continuum model yields broadly similar results to those arising from a kinetic model. Furthermore, we show that in two dimensional (2-D) continuum models the modeling efficiency based on the Akaike information criterion (AIC) is enhanced for equilibrium vs kinetic with no significant reduction in model performance. This is because fewer parameters are needed for the equilibrium model compared to the kinetic model. Two major transport regimes are identified in the transport of NP within porous media. The first regime is characterized by higher particle-surface attachment affinity than particle-particle attachment affinity, and operative transport mechanisms of physicochemical filtration, blocking, and physical retention. The second regime is characterized by the domination of particle-particle attachment tendency over particle-surface affinity. In this regime although physicochemical filtration as well as straining may still be operative, ripening is predominant together with agglomeration and further subsequent retention. In both regimes careful assessment of NP fate and transport is necessary since certain combinations of concurrent transport phenomena leading to large migration distances are possible in either case.