A robust upscaling of the effective particle deposition rate in porous media

In the upscaling from pore- to continuum (Darcy) scale, reaction and deposition phenomena at the solid-liquid interface of a porous medium have to be represented by macroscopic reaction source terms. The effective rates can be computed, in the case of periodic media, from three-dimensional microscopic simulations of the periodic cell. Several computational and semi-analytical models have been studied in the field of colloid filtration to describe this problem. They typically rely on effective deposition rates defined by complex fitting procedures, neglecting the advection-diffusion interplay, the pore-scale ow complexity, and assuming slow reactions (or large Peclet numbers). Therefore, when these rates are inserted into general macroscopic transport equations, they can lead to several conceptual inconsistencies and significant errors. To study more accurately the dependence of the deposition on the flow parameters, in this work, we advocate a clear distinction between the surface processes (that altogether defines the so-called attachment efficiency), and the pore-scale processes. With this approach, valid when colloidal particles are small enough, we study Brownian and gravity-driven deposition on face-centred cubic (FCC) spherical arrangements, and de ne a robust upscaling based on a linear effective reaction rate. The case of partial deposition, defined by an attachment probability, is studied and the limit of perfect sink is retrieved as a particular case. We introduce a novel upscaling approach and a particularly convenient computational setup that allows the direct computation of the asymptotic stationary value of effective rates. This allows to drastically reduce the computational domain down to the scale of the single repeating periodic unit: the savings are ever more noticeable in the case of higher Peclet numbers, when larger physical times are needed to reach the asymptotic regime, and thus, equivalently, a much larger computational domain and simulation time would be needed in a traditional setup. We show how this new de nition of deposition rate is more robust and extendable to the whole range of Peclet numbers; it also is consistent with the classical heat and mass transfer literature