Mixing and spreading of gravity currents in heterogeneous porous media

This study utilizes high-fidelity numerical simulations to demonstrate that while permeability heterogeneity generally accelerates gravity current migration and enhances dissolution through dispersive effects, the interplay between instability size, correlation length, and density stratification ultimately determines dissolution efficiency, with homogeneous media often outperforming heterogeneous ones except at low Rayleigh numbers.

The Gist

Imagine you are pouring a thick, heavy syrup into a glass of water. In a perfect, uniform glass, the syrup would slide down the side in a smooth, predictable sheet. But what if the glass wasn't empty? What if it was packed with a sponge that had some holes the size of a pin and others the size of a marble?

This is the world of gravity currents in porous media (like underground rock or soil) that Albert Jiménez-Ramos and Juan J. Hidalgo explored in their paper. They used powerful computer simulations to watch how fluids mix and move through these "spongy" underground layers, specifically looking at how the unevenness of the rock changes the game.

Here is the story of their findings, broken down into everyday concepts:

The Setup: Two Types of Fluids

The researchers studied two scenarios, like two different types of parties:

1. The Stable Party: Imagine pouring oil into water. The oil is lighter and stays on top, or if it's saltwater, it sinks smoothly. The fluids don't fight each other; they just slowly mix at the boundary.
2. The Unstable Party: This is like pouring a fluid that is heavy in the middle but light at the edges. It's a chaotic situation where the heavy parts want to sink and the light parts want to float, creating "fingers" of fluid that dive down or shoot up, mixing everything violently.

The "Sponge" Effect: Heterogeneity

In the real world, underground rock isn't uniform. It's a heterogeneous mix of hard, tight spots (low permeability) and loose, open spots (high permeability). The researchers treated this like a sponge with random holes.

What they found:

• The Barrier Effect: When the fluid hits a tight, hard spot in the rock, it gets stuck. It's like trying to run through a crowd; if there's a wall (a low-permeability zone), you have to go around it. This "barrier" usually slows down the mixing process because the fluid can't get through easily.
• The Trap: Sometimes, the light fluid gets trapped in a tight pocket surrounded by heavy fluid. It's like a bubble trapped in a net. Eventually, this trapped bubble dissolves rapidly, creating a little explosion of mixing.

The Big Surprise: Chaos vs. Order

The most interesting discovery was how the "sponge" (the rock) interacted with the "fingers" (the unstable mixing).

• In the Stable Case: The uneven rock acted like a disperser. It spread the fluid out, making the mixing zone wider and slower. It was like running through a forest; you get scattered, and you don't get very far very fast.
• In the Unstable Case: You might think the rock would slow down the chaotic fingers, but it didn't. The chaotic "fingering" was so strong that it overpowered the rock's tendency to scatter things. The fingers punched through the rock's barriers.
  • The Result: The mixing became more efficient in the unstable case than in the stable one. The fingers made the interface between fluids narrower and sharper, allowing them to dissolve into each other faster than if the rock were perfectly smooth.

The "Speed vs. Mixing" Trade-off

The paper highlights a tug-of-war between how fast the fluid moves and how well it mixes:

• High Speed (High Rayleigh Number): When the fluid is very dense and moves fast, it tends to stay in a tight stream. In a uniform rock, it mixes well. But in a bumpy rock, the "barrier effect" wins. The fluid gets blocked, moves faster along the easy paths, but mixes less overall.
• Low Speed (Low Rayleigh Number): When the fluid moves slowly, diffusion (the natural tendency to spread out) does the work. Here, the uneven rock actually helps. The early chaos caused by the rock's bumps makes the fluid mix better than it would in a smooth, uniform rock.

The "Anisotropy" Factor: Direction Matters

The researchers also looked at the direction of the rock's holes.

• Horizontal Layers (Like a Layer Cake): If the rock has horizontal layers of hard and soft spots, it acts like a series of shelves. The sinking fingers hit a shelf and stop. This shuts down the mixing quickly.
• Vertical Layers (Like a Stack of Papers): If the layers are vertical, the fingers can slide down them easily, but the whole current moves slower because it has to navigate the vertical walls.

The Bottom Line

The paper concludes that the efficiency of mixing depends on a delicate balance:

1. The Size of the Chaos: How big are the "fingers" of mixing?
2. The Size of the Rock's Bumps: How big are the holes in the sponge?

If the fingers are small and the rock's bumps are big (high speed, high variance), the rock acts as a barrier, slowing down mixing and letting the fluid travel further.
If the fingers are large and the rock's bumps are small (low speed), the rock's bumps actually help kick-start the mixing, making it more efficient than a smooth rock would.

In short: Nature's "spongy" underground doesn't just slow things down; it changes the rules of the game. Sometimes it blocks the flow, and sometimes, if the fluid is chaotic enough, it helps the fluids mix together faster than they would in a perfect, smooth world.

Technical Summary

1. Problem Statement

The study investigates the complex interplay between gravity currents (fluids driven by density differences) and porous media heterogeneity. While previous research has examined gravity currents in homogeneous media or the effects of heterogeneity on simple mixing, this work specifically targets the combined impact of:

• Permeability Heterogeneity: Represented by log-normal stochastic fields with varying variances and correlation lengths.
• Density Stratification: Comparing stable scenarios (linear density-concentration relationship, e.g., saltwater intrusion) against unstable scenarios (non-monotonic density relationship, e.g., CO₂ injection in brine where fingering instabilities occur).
• Dissolution and Mixing: Understanding how heterogeneity alters the dissolution flux, migration speed, and spreading of the current, which are critical for applications like aquifer remediation and geological carbon storage.

2. Methodology

The authors employed high-fidelity numerical simulations to solve the governing equations for two-dimensional, miscible flow in porous media.

• Governing Equations: The system is modeled using the Darcy equation for flow and the advection-diffusion equation for concentration transport. The equations are non-dimensionalized using the Rayleigh number ($Ra$).
• Permeability Modeling: Heterogeneity is modeled as a log-normal stochastic process. Key parameters varied include:
  • Variance ($\sigma^2_{\log k}$): 1, 3, and 5.
  • Correlation Lengths ($\lambda_x, \lambda_z$): Used to create isotropic and anisotropic fields (horizontal vs. vertical stratification).
  • Anisotropy Ratio ($\gamma = \lambda_x / \lambda_z$): Tested values of 0.5, 1, 2, and 10.
• Density Laws:
  • Stable Case: Linear relationship ($\rho = -c$).
  • Unstable Case: Non-monotonic relationship mimicking CO₂-brine systems, featuring an intermediate density maximum that triggers fingering instabilities.
• Numerical Setup:
  • Solver: SECUReFOAM (OpenFOAM package) using the finite volume method.
  • Mesh: High-resolution structured grid ($16,000 \times 800$ cells) to resolve thin fingers and concentration gradients.
  • Metrics: The study quantified buoyant mass, scalar dissipation rate ($\chi$), dissolution flux, interface length, and interface width.

3. Key Contributions

The paper provides a comprehensive analysis of how heterogeneity modifies the physics of gravity currents, specifically challenging the assumption that heterogeneity always enhances mixing. Key contributions include:

• Counteracting Mechanisms: Demonstrating that in unstable (fingering) regimes, the dispersive effect of heterogeneity is counteracted by the fingering instability, leading to different outcomes than in stable regimes.
• Anisotropy Effects: Identifying a "barrier effect" in horizontally stratified media (high $\gamma$) that suppresses vertical convection and reduces dissolution, contrasting with the behavior in vertically stratified media.
• Variance Thresholds: Establishing that the variance of the permeability field only significantly enhances dissolution efficiency compared to homogeneous cases when the Rayleigh number is low.
• Blob Trapping: Detailed the mechanism where buoyant fluid gets trapped in low-permeability zones, forming "blobs" that dissolve rapidly and create localized upward fingering, even in stable stratification.

4. Key Results

A. Homogeneous vs. Heterogeneous Media

• Stable Case: Heterogeneity acts primarily as a dispersive agent. It widens the mixing interface, reduces the concentration gradient, and consequently decreases dissolution. However, it increases the migration speed of the current nose.
• Unstable Case: Heterogeneity accelerates the onset of convection. While this initially increases dissolution, the dispersive effect of the heterogeneity eventually narrows the interface compared to the stable case but creates a complex interaction where the homogeneous case often outperforms the heterogeneous case in total dissolved mass for high $Ra$.

B. Impact of Anisotropy ($\gamma$)

• Horizontal Stratification ($\gamma > 1$): Creates a barrier effect. Low-permeability horizontal layers trap buoyant fluid, suppressing finger formation and shutting down convection locally. This leads to a reduction in total dissolution compared to homogeneous media, despite an initial spike in dissolution flux.
• Vertical Stratification ($\gamma < 1$): Acts as a barrier to horizontal migration but favors vertical finger descent. This results in an intermittent dissolution behavior as the current crosses permeability barriers, with an earlier onset of the "shutdown" regime (where dense fluid accumulation at the bottom halts mixing).

C. Impact of Variance ($\sigma^2_{\log k}$)

• High Variance: Increases the dispersive effect, widening the interface and reducing the concentration gradient driving dissolution.
• Efficiency Trade-off: For high Rayleigh numbers (strong convection), heterogeneous cases generally dissolve less mass than homogeneous cases because the barrier and dispersive effects dominate.
• Low Rayleigh Numbers: Heterogeneity can enhance dissolution compared to the homogeneous case due to the early onset of convection and interface elongation, as diffusion plays a more dominant role relative to the instability size.

D. Interface Dynamics

• In unstable cases, the interface is generally narrower than in stable cases due to the dominance of fingering over diffusion.
• The presence of trapped "blobs" in low-permeability regions causes oscillatory behavior in interface length and dissolution flux, particularly in anisotropic media.

5. Significance and Implications

This research is critical for optimizing Geological Carbon Storage (GCS) and aquifer remediation:

• Carbon Storage: The findings suggest that in high-permeability, high-convection scenarios (high $Ra$), natural heterogeneity might actually hinder the solubility trapping of CO₂ compared to idealized homogeneous models. Engineers must account for the "barrier effect" of horizontal stratification which can limit vertical mixing.
• Remediation: For DNAPL remediation, understanding that heterogeneity can trap fluids in "blobs" that dissolve rapidly helps in predicting the longevity of remediation efforts.
• Modeling Guidelines: The study highlights that simple upscaling of heterogeneous media is insufficient. The interaction between the size of convective instabilities and the correlation length of the permeability field is the governing factor for dissolution efficiency. Optimal dissolution requires a specific balance between the Rayleigh number and the strength of heterogeneity.

In conclusion, the paper establishes that heterogeneity does not universally enhance mixing; rather, its impact is conditional on the stability of the density stratification and the relative scales of the flow instabilities versus the geological heterogeneity.