Abstract:The effective mixing behavior of solutes in porous media is fundamentally connected to the development of a local mixing interface between the two initial solutions, which is characterized by a complex lamellar structure. The deformation of the interface is controlled by the interplay of advection and diffusion, which generate the mechanisms of lamella stretching and shrinking, respectively. Based on the results of pore-scale numerical simulations, we develop a mechanistic single parabolic lamella model (SPLM) to capture the interface evolution across various temporal and Péclet number scales. The model shows near-perfect agreement with a 2D parallel plates scenario and promising results for a 3D porous medium. The SPLM model also establishes Péclet regimes for the equilibrium area and temporal regimes for the transient behavior of the interface. These findings represent a step forward towards eventually incorporating mixing limitation into general macroscopic reactive transport models.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to solve the problem of effective solute mixing behavior in porous media. Specifically, the paper focuses on the development of the local mixing interface between two initial solutions, which has a complex layered structure. The deformation of the interface is controlled by the interaction between convection and diffusion, which respectively lead to the stretching and shrinking of the layered structure. Based on the results of pore - scale numerical simulations, the researchers developed a mechanical single parabolic lamella model (Single Parabolic Lamella Model, SPLM) to capture the interface evolution at different time and Péclet number (Pe) scales. This model shows an almost perfect fit in the two - dimensional parallel - plate scenario and also achieved encouraging results in three - dimensional porous media. In addition, the SPLM model also established the time scales of the equilibrium area and the interface transient behavior, and these findings will help to finally incorporate mixing limitations into general macroscopic reaction - transport models.
### Key point summary
1. **Research background**:
- Reaction - transport with mixing limitations is widespread in porous media and involves multiple fields such as groundwater pollution and remediation, CO2 sequestration, enhanced oil recovery, and river sediment ecology.
- As a key process for reactant combination, the development of a prediction model for mixing is crucial.
- The complex structure in porous media leads to velocity variations at the pore scale, which cannot be resolved at the Darcy scale and further affect the mixing behavior at larger scales.
2. **Research methods**:
- Use high - resolution numerical simulation data to extract concentration isosurfaces and analyze the evolution of their surface areas.
- Propose the single parabolic lamella model (SPLM), considering the influence of convection and diffusion on the evolution of the mixing - layer interface.
- Verify the validity of the model by comparing simulation results under different Péclet numbers.
3. **Main findings**:
- The interface area increases over time but eventually reaches a stagnant plateau.
- In the early time scale, the growth of the interface area exhibits super - linear behavior.
- The area value in the plateau period increases with the increase of the Péclet number, but not linearly.
4. **Model application**:
- The SPLM model has been verified in two - dimensional parallel - plate and three - dimensional porous media.
- The model performs excellently in two - dimensional parallel - plate and also achieved relatively good results in three - dimensional porous media, although there are some deviations in the early - time behavior.
### Formula summary
- **Stretching rate**:
\[
\frac{\partial s}{\partial t}\big|_{\text{stretch}} = c_A\bar{v}\frac{1}{1 + \frac{s}{\lambda}}
\]
where \(c_A\) is a dimensionless stretching proportionality constant, \(\bar{v}\) is the average flow velocity, and \(\lambda\) is the velocity - related length.
- **Shrinking rate**:
\[
\frac{\partial s}{\partial t}\big|_{\text{shrink}} = - c_DD\frac{s}{h^2}
\]
where \(c_D\) is a dimensionless shrinking proportionality constant, \(D\) is the diffusion coefficient, and \(h\) is the layer width.
- **Equilibrium condition**:
\[
s^*_{\text{equil}}=\frac{\lambda^*}{2}\left(\sqrt{1 + \frac{4c\text{Pe}}{\lambda^*}} - 1\right)
\]
where \(s^*=\frac{s}{h}\), \(\lambda^*=\frac{\lambda}{h}\), \(c = \frac{c_A}{c_D}\), \(\text{Pe}=\frac{\bar{v}h}{D}\).
- **Interface area growth**:
\[
G=\frac{A}{...}