Two Film Approach to Continuum Scale Mixing and Dispersion with Equilibrium Bimolecular Reaction

Simon A. Mathias,Diogo Bolster,Sergii Veremieiev

DOI: https://doi.org/10.1007/s11242-024-02091-y

IF: 3.6097

2024-05-28

Transport in Porous Media

Abstract:Reliable reactive transport models require careful separation of mixing and dispersion processes. Here we treat displacing and displaced fluids as two separate fluid phases and invoke Whitman's classical two-film theory to model mass transfer between the two phases. We use experimental data from Gramling's bimolecular reaction experiment to assess model performance. Gramling's original model involved just three coupled PDEs. In this context, our new formulation leads to a set of seven coupled PDEs but only requires the specification of two extra parameters, associated with the mass transfer coefficient and its dependence on time. The two film mass transfer model provides a simple and theoretically based method for separating mixing from dispersion in Eulerian continuum-scale methods. The advantage of this approach over existing methods is that it enables the simulation of equilibrium chemical reactions without having to invoke unrealistically small reaction rate coefficients. The comparison with Gramling's experimental data confirms that our proposed method is suitable for simulating realistic and complicated bimolecular reaction behaviour. However, further work is needed to explore alternative methods for avoiding the need of a time-dependent mass transfer rate coefficient.

engineering, chemical

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to more accurately simulate reactive transport phenomena in porous media, especially in the case of bimolecular chemical reactions. Specifically, the author aims to distinguish between mixing and dispersion processes and improve the description of these processes in existing models to avoid over - estimating the reaction rate.

### Background problems
The traditional advection - dispersion equation (ADE) cannot effectively distinguish between mixing and dispersion processes when simulating solute transport in porous media. Since the dispersion coefficient is usually much larger than the molecular diffusion coefficient, this leads to an over - estimation of the degree of mixing, and thus an over - estimation of the reaction rate. For example, in the study by Gramling et al. (2002), the amount of CuEDTA²⁻ generated predicted by the ADE model was approximately 20% higher than the actual experimental results.

### Core problems of the paper
1. **Distinguishing between mixing and dispersion**: Traditional models regard mixing and dispersion as the same process, while in fact they are different physical mechanisms. Mixing is caused by molecular diffusion, while dispersion is caused by flow heterogeneity.
2. **Accuracy of the reaction rate coefficient**: Existing models often need to introduce unrealistically small reaction rate coefficients to fit experimental data, which is not in line with the actual situation.
3. **Time - dependent mass transfer coefficient**: Many models assume that the mass transfer coefficient changes with time, but lack a theoretical basis.

### Solutions
To solve the above problems, this paper proposes a new method based on Whitman's two - film theory. This method regards the displacing fluid and the displaced fluid as two independent phases and distinguishes between mixing and dispersion processes by introducing a mass transfer term. The advantages of this method are:
- It can more accurately simulate equilibrium chemical reactions without using unrealistic small reaction rate coefficients.
- It provides a simple and theoretically - based method to separate mixing and dispersion processes.

### Model characteristics
- **Seven sets of coupled partial differential equations (PDE)**: The new model contains seven sets of coupled PDEs, but only two additional parameters (the mass transfer coefficient and its time - dependence) need to be specified.
- **Mass transfer coefficient**: The mass transfer coefficient \( k \) is defined as:
\[
k=\frac{\alpha D_0}{2\phi b}
\]
where \( \alpha \) is the fluid - fluid interface area, \( D_0 \) is the molecular diffusion coefficient, \( \phi \) is the effective porosity, and \( b \) is the thickness of the two - film system.
- **Time - dependence**: The mass transfer coefficient \( k \) changes with time in the form of:
\[
k = k_0t^{-\beta}
\]
where \( k_0 \) and \( \beta \) are empirical parameters.

### Experimental verification
The author uses the experimental data of Gramling et al. to verify the validity of the model. The experiments include non - reactive and reactive tracer experiments at different flow rates. By comparing the model predictions and experimental results, the applicability of the new method in simulating complex bimolecular reaction behaviors is verified.

In conclusion, this paper aims to improve the existing reactive transport models by introducing the two - film theory, so as to more accurately simulate the mixing, dispersion and reaction processes in porous media.