Heat and super-diffusive melting fronts in unsaturated porous media

TL;DR

This paper develops a nonlinear diffusion model for melting fronts in unsaturated porous media, revealing super-diffusive spreading influenced by pore size distribution and geological heating effects.

Contribution

It introduces an analytical solution for super-diffusive melting fronts considering pore size distribution and energy conservation, with applications to geological subsurface melting.

Findings

Melting fronts exhibit super-diffusive spreading with a specific diffusion exponent.

The model accurately predicts the depth of melting fronts in geological systems.

Energy conservation limits the melting front dynamics in the derived solution.

Abstract

When water is present in a medium with pore sizes in a range around 10nm the corresponding freezing point depression will cause long range broadening of a melting front. Describing the freezing-point depression by the Gibbs-Thomson equation and the pore size distribution by a power law, we derive a non-linear diffusion equation for the fraction of melted water. This equation yields super-diffusive spreading of the melting front with a diffusion exponent which is given by the spatial dimension and the exponent describing the pore size distribution. We derive this solution analytically from energy conservation in the limit where all the energy is consumed by the melting and explore the validity of this approximation numerically. Finally, we explore a geological application of the theory to the case of one-dimensional sub-surface melting fronts in granular or soil systems. These fronts, which are produced by heating of the surface, spread at a super-diffusive rate and affect the subsurface to significantly larger depths than would a system without the effects of freezing point depression.