BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Matteo Icardi

Matteo Icardi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynphysics.comp-phcond-mat.mtrl-scicond-mat.stat-mechInformation Theorymath.APmath.ITphysics.app-phphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Emergence of non-Fickian transport in truncated pluri-Gaussian permeability fields

We present a numerical simulation study of advective-diffusive scalar transport in three-dimensional high-contrast discontinuous permeability fields, generated with a truncated pluri-Gaussian geostatistical approach. A range of permeability contrasts, correlation lengths, and Péclet numbers are studied to characterise the transition to non-Fickian transport behaviour. This is triggered by high permeability contrasts between different zones and is enhanced by the presence of connected higher permeability channels, which are characterised by high advective flow velocities. In this case, the overall transport behaviour deviate from the macroscopic advection-dispersion model based on a Fickian analogy. The numerical experiments are run with an Eulerian approach using a novel unified numerical framework based on the finite-volume library \of, for i) generating random pluri-Gaussian porous media, ii) solving the steady state Darcy-scale flow, iii) solving the advection diffusion equation, iv) computing post-processing quantities such as first order statistics, spatial probability density functions and breakthrough curves. We identify a hierarchy of non-Fickian transport triggering factors, the strength of permebility contrast being the pivotal driver. The Péclet number and the characteristic length at which facies transitions are observed as secondary factors. Transport remains Fickian when the facies conductivities differ by up to one order of magnitude. Greater permeability contrasts act strengthen the emergence of fast flow channels leading to non-Fickian transport.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Heterogeneous Multi-Rate mass transfer models in OpenFOAM

We implement the Multi-Rate Mass Transfer (MRMT) model for mobile-immobile transport in porous media within the open-source finite volume library \textsc{OpenFOAM}\reg \citep{Foundation2014}. Unlike other codes available in the literature [Geiger, S., Dentz, M., Neuweiler, I., SPE Reservoir Characterisation and Simulation Conference and Exhibition (2011); Silva, O., Carrera, J., Dentz, M., Kumar, S., Alcolea, A., Willmann, M., Hydrology and Earth System Sciences 13, (2009)], we propose an implementation that can be applied to complex three-dimensional geometries and highly heterogeneous fields, where the parameters of the MRMT can arbitrarily vary in space. Furthermore, being built over the widely diffused OpenFOAM\reg library, it can be easily extended and included in other models, and run in parallel. We briefly describe the structure of the < multiContinuumModels > library that includes the formulation of the MRMT based on the works of [Haggerty, R., Gorelick, S.M., Water Resources Research 31, (1995)] and [F. Municchi and M. Icardi Phys. Rev. Research 2, 013041, (2020)]. The implementation is verified against benchmark solutions and tested on two- and three-dimensional random permeability fields. The role of various physical and numerical parameters, including the transfer rates, the heterogeneities, and the number of terms in the MRMT expansions is investigated. Finally, we illustrate the significant role played by heterogeneity in the mass transfer when permeability and porosity are represented using Gaussian random fields.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Macroscopic models for filtration and heterogeneous reactions in porous media

Derivation of macroscopic models for advection-diffusion processes in the presence of dominant heterogeneous (e.g., surface) reactions using homogenisation theory or volume averaging is often deemed unfeasible due to the strong coupling between scales that characterise such systems. In this work, we show how the upscaling can be carried out by applying and extending the methods presented in literature. The approach relies on the decomposition of the microscale concentration into a reactive component, given by the eigenfunction of the advection-diffusion operator, the associated eigenvalue which represents the macroscopic effective reaction rate, and a non-reactive component. The latter can be then upscaled with a two-scale asymptotic expansion and the final macroscopic equation is obtained for the leading order. The same method can also be used to overcome another classical assumption, namely of non solenoidal velocity fields, such as the case of deposition of charged colloidal particles driven by electrostatic potential forces. The whole upscaling procedure, which consists in solving three cell problems, is implemented for arbitrarily complex two- and three-dimensional periodic structures using the open-source finite volume library OpenFOAM. We provide details on the implementation and test the methodology for two-dimensional periodic arrays of spheres, and we compare the results against fully resolved numerical simulations, demonstrating the accuracy and generality of the upscaling approach. The effective velocity, dispersion and reaction coefficients are obtained for a wide range of Péclet and surface Damk{ö}hler numbers, and for Coulomb-like forces to the grains.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimisation of crust freezing in meat processing via Computational Fluid Dynamics

In this work a numerical model for two-dimensional axisymmetric continuous freezing by impingement of processed meat or similar products in food industry moving along a conveyor belt is presented. The model represents a more computationally efficient alternative to solve conjugate heat transfer between a fluid and a solid, accompanied by phase change in some constituents of the solid phase. In the model presented here it is assumed that the solid can be represented as an homogeneous medium, with its thermophysical properties depending on the temperature. The impingement freezing model is conceived to be valid for highly processed vegetarians products or meat such as sausages, mince or ham freezing. Furthermore, this approach is much simpler in terms of computational cost whilst it still captures the complexity of continuous freezing under industrial setting. The methodology is implemented as a new solver in the widely used open-source Computational Fluid Dynamics (CFD) library OpenFOAM. Overall, highly non-linear freezing behaviour was found due to the phase change inside the solid and the associated heat of fusion. We studied the effect of high fluid Reynolds numbers as well as investigating the optimal distance between the jet and the solid surface for different speeds of the conveyor. We found that the maximum freezing is obtained positioning the jet at a distance H = 7.2D (where D is the diameter of the impinging jet) and setting the speed of the conveyor such that the Péclet number of the solid is $\text{Pe}_{\text{s}}=8244$. The methodology developed allows to obtain detailed insight on the freezing process for various impingement configurations at a minimum computational cost using a freely available open-source tool.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probability density function (PDF) models for particle transport in porous media

Mathematical models based on probability density functions (PDF) have been extensively used in hydrology and subsurface flow problems, to describe the uncertainty in porous media properties (e.g., permeability modelled as random field). Recently, closer to the spirit of PDF models for turbulent flows, some approaches have used this statistical viewpoint also in pore-scale transport processes (fully resolved porous media models). When a concentration field is transported, by advection and diffusion, in a heterogeneous media, in fact, spatial PDFs can be defined to characterise local fluctuations and improve or better understand the closures performed by classical upscaling methods. In the study of hydrodynamical dispersion, for example, PDE-based PDF approach can replace expensive and noisy Lagrangian simulations (e.g. trajectories of drift-diffusion stochastic processes). In this work we derive a joint position-velocity Fokker-Planck equation to model the motion of particles undergoing advection and diffusion in in deterministic or stochastic heterogeneous velocity fields. After appropriate closure assumptions, this description can help deriving rigorously stochastic models for the statistics of Lagrangian velocities. This is very important to be able to characterise the dispersion properties and can, for example, inform velocity evolution processes in Continuous Time Random Walk (CTRW) dispersion models. The closure problem that arises when averaging the Fokker Planck equation shows also interesting similarities with the mixing problem and can be used to propose alternative closures for anomalous dispersion.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Generalised Multi-Rate Models for conjugate transfer in heterogeneous materials

We propose a novel macroscopic model for conjugate heat and mass transfer between a \emph{mobile region}, where advective transport is significant, and a set of \emph{immobile regions} where diffusive transport is dominant. Applying a spatial averaging operator to the microscopic equations, we obtain a \emph{multi-continuum} model, where an equation for the average concentration in the mobile region is coupled with a set of equations for the average concentrations in the immobile regions. Subsequently, by mean of a spectral decomposition, we derive a set of equations that can be viewed as a generalisation of the multi-rate mass transfer (MRMT) model, originally introduced by Haggerty & Gorelick. This new formulation does not require any assumption on local equilibrium or geometry. We then show that the MRMT can be obtained as the leading order approximation, when the mobile concentration is in local equilibrium. The new Generalised Multi-Rate Mode (GMRM) has the advantage of providing a direct method for calculating the model coefficients for immobile regions of arbitrary shapes, through the solution of appropriate micro-scale cell problems. An important finding is that a simple re-scaling or re-parametrisation of the transfer rate coefficient (and thus, the memory function) is not sufficient to account for the flow field in the mobile region and the resulting non-uniformity of the concentration at the interfaces between mobile and immobile regions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Molecular Information Delivery in Porous Media

Information delivery via molecular signals is abundant in nature and potentially useful for industry sensing. Many propagation channels (e.g., tissue membranes and catalyst beds) contain porous medium materials and the impact this has on communication performance is not well understood. Here, communication through realistic porous channels is investigated for the first time via statistical breakthrough curves. Assuming that the number of arrived molecules can be approximated as a Gaussian random variable and using fully resolved computational fluid dynamics results for the breakthrough curves, the numerical results for the throughput, mutual information, error probability, and information diversity gain are presented. Using these numerical results, the unique characteristics of the porous medium channel are revealed.

0 citations0 reviewsNo code linkedOpen access