In this work, the influence of pore space geometry on solute transport in porous media is investigated performing computational fluid dynamics pore-scale simulations of fluid flow and solute transport. The three-dimensional periodic domains are obtained from three different pore structure configurations, namely, face-centered-cubic (fcc), body-centered-cubic (bcc), and sphere-in-cube (sic) arrangements of spherical grains. Although transport simulations are performed with media having the same grain size and the same porosity (in fcc and bcc configurations), the resulting breakthrough curves present noteworthy differences, such as enhanced tailing. The cause of such differences is ascribed to the presence of recirculation zones, even at low Reynolds numbers. Various methods to readily identify recirculation zones and quantify their magnitude using pore-scale data are proposed. The information gained from this analysis is then used to define macroscale models able to provide an appropriate description of the observed anomalous transport. A mass transfer model is applied to estimate relevant macroscale parameters (hydrodynamic dispersion above all) and their spatial variation in the medium; a functional relation describing the spatial variation of such macroscale parameters is then proposed.

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 often rely on effective deposition rates defined by simple linear reaction ODEs, neglecting the advection-diffusion interplay, and assuming slow reactions (or large Péclet numbers). Therefore, when these rates are inserted into general macroscopic transport equations, they can lead to conceptual inconsistencies and, therefore, often qualitatively wrong results. In this work, we study the upscaling of Brownian deposition on face-centred cubic (FCC) spherical arrangements using a linear effective reaction rate, defined by volume averaging, and a macroscopic advection-diffusion-reaction equation. The case of partial deposition, defined by an attachment probability, is studied and the limit of ideal deposition is retrieved as a particular case. We make use of 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 Péclet 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 definition of deposition rate is more robust and extendable to the whole range of Péclet numbers; it also is consistent with the classical heat and mass transfer literature.

No abstract

No abstract

Porous and heterogeneous materials are found in many applications from composites, membranes, chemical reactors, and other engineered materials to biological matter and natural subsurface structures. In this work we propose an integrated approach to generate, study and upscale transport equations in random and periodic porous structures. The geometry generation is based on random algorithms or ballistic deposition. In particular, a new algorithm is proposed to generate random packings of ellipsoids with random orientation and tunable porosity and connectivity. The porous structure is then meshed using locally refined Cartesian-based or unstructured strategies. Transport equations are thus solved in a finite-volume formulation with quasi-periodic boundary conditions to simplify the upscaling problem by solving simple closure problems consistent with the classical theory of homogenisation for linear advectiondiffusion-reaction operators. Existing simulation codes are extended with novel developments and integrated to produce a fully open-source simulation pipeline. A showcase of a few interesting three-dimensional applications of these computational approaches is then presented. Firstly, convergence properties and the transport and dispersion properties of a periodic arrangement of spheres are studied. Then, heat transfer problems are considered in a pipe with layers of deposited particles of different heights, and in heterogeneous anisotropic materials.

Transport and particulate processes are ubiquitous in environmental, industrial and biological applications, often involving complex geometries and porous media. In this work we present a general population balance model for particle transport at the pore-scale, including aggregation, breakage and surface deposition. The various terms in the equations are analysed with a dimensional analysis, including a novel collision-induced breakage mechanism, and split into one- and two-particles processes. While the first are linear processes, they might both depend on local flow properties (e.g. shear). This means that the upscaling (via volume averaging and homogenisation) to a macroscopic (Darcy-scale) description requires closures assumptions. We discuss this problem and derive an effective macroscopic term for the shear-induced events, such as breakage caused by shear forces on the transported particles. We focus on breakage events as prototype for linear shear-induced events and derive upscaled breakage frequencies in periodic geometries, starting from nonlinear power-law dependence on the local fluid shear rate. Results are presented for a two-dimensional channel flow and a three dimensional regular arrangement of spheres, for arbitrarily fast (mixing-limited) events. Implications for linearised shear-induced collisions are also discussed. This work lays the foundations of a new general framework for multiscale modelling of particulate flows.

No abstract

No abstract