Impact of saturation on evaporation-driven density instabilities in porous media: mathematical and numerical analysis

TL;DR

This paper investigates how water saturation levels affect the formation of density-driven salt concentration instabilities during evaporation in porous media, using mathematical analysis and numerical simulations to predict instability onset.

Contribution

It provides new criteria for the onset of density instabilities based on linear stability analysis considering saturation effects, supported by numerical simulations.

Findings

Density instabilities depend on evaporation rate, permeability, and saturation.

Decreased water saturation reduces salt transport and delays instability formation.

The study offers predictive tools for instability development in saline porous media.

Abstract

Evaporation from a porous medium partially saturated with saline water, causes the salinity (salt concentration) to increase near the top of the porous medium as water leaves while salt stays behind. As the density of the water increases with increased salt concentration, the evaporation leads to a gravitational unstable setting, where density instabilities can form. Whether density instabilities form, depends on a large range of parameters like the evaporation rate and intrinsic permeability of the porous medium, but also on the water saturation. As water saturation decreases, the storage, convection and diffusion of salt also decrease, which all influence the onset of instabilities. By performing a linear stability analysis on the governing equations, we give criteria for onset of instabilities. Numerical simulations give information about the further development of these instabilities. With this knowledge we can predict whether and when density instabilities occur, and how they will influence the further development of salt concentration in the porous medium.