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

Federico Municchi

Federico Municchi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynphysics.comp-phcond-mat.mtrl-scimath.APphysics.app-ph

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
5topics
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

5 published item(s)

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
preprint2019arXiv

Computational Analysis of Interfacial Dynamics in Angled Hele-Shaw Cells: Instability Regimes

We present a theoretical and numerical study on the (in)stability of the interface between two immiscible liquids, i.e., viscous fingering, in angled Hele-Shaw cells across a range of capillary numbers ($Ca$). We consider two types of angled Hele-Shaw cells: diverging cells with a positive depth gradient and converging cells with a negative depth gradient, and compare those against parallel cells without a depth gradient. A modified linear stability analysis is employed to derive an expression for the growth rate of perturbations on the interface and for the critical capillary number ($Ca_c$) for such tapered Hele-Shaw cells with small gap gradients. Based on this new expression for $Ca_c$, a three-regime theory is formulated to describe the interface (in)stability: (i) in Regime I, the growth rate is always negative, thus the interface is stable; (ii) in Regime II, the growth rate remains zero (parallel cells), changes from negative to positive (converging cells), or from positive to negative (diverging cells), thus the interface (in)stability possibly changes type at some location in the cell; (iii) in Regime III, the growth rate is always positive, thus the interface is unstable. We conduct three-dimensional direct numerical simulations of the full Navier--Stokes equations, using a phase field method to enforce surface tension at the interface, to verify the theory and explore the effect of depth gradient on the interface (in)stability. We demonstrate that the depth gradient has only a slight influence in Regime I, and its effect is most pronounced in Regime III. Finally, we provide a critical discussion of the stability diagram derived from theoretical considerations versus the one obtained from direct numerical simulations.

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