# Generalised Multi-Rate Models for conjugate transfer in heterogeneous materials

## TL;DR

This paper introduces a generalized multi-rate model for conjugate heat and mass transfer in heterogeneous materials, capturing complex interactions without assuming local equilibrium or specific geometries, and enabling direct coefficient calculation.

## Contribution

It develops a new GMRM framework that generalizes MRMT, allowing for arbitrary shapes and micro-scale computations, improving modeling flexibility and accuracy.

## Key findings

- Derivation of a spectral decomposition-based GMRM model.
- Demonstration that MRMT is a leading order approximation of GMRM.
- Highlighting the limitations of simple re-scaling of transfer coefficients.

## Abstract

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.