Unified analysis of algorithms for equilibrium, non-equilibrium, and hysteresis models of phase transition in permafrost

TL;DR

This paper introduces a unified analysis framework for heat flow models with phase transitions in permafrost, accommodating equilibrium, non-equilibrium, and hysteresis effects, with new algorithms and convergence proofs.

Contribution

It develops a comprehensive analysis and numerical methods that handle various phase transition regimes in a unified manner, including hysteresis effects.

Findings

Proposed numerical algorithms with proven convergence.

Unified analysis covering equilibrium, non-equilibrium, and hysteresis cases.

Numerical illustrations demonstrating model behavior.

Abstract

In this paper we consider a nonlinear partial differential equation describing heat flow with ice-water phase transition in permafrost soils. Such models and their numerical approximations have been well explored in the applications literature. In this paper we describe a new direction in which the allow relaxation and hysteresis of the phase transition which introduce additional nonlinear terms and complications for the analysis. We present numerical algorithms as well as analysis of the well-posedness and convergence of the fully implicit iterative schemes. The analysis we propose handles the equilibrium, non-equilibrium, and hysteresis cases in a unified way. We also illustrate with numerical examples for a model ODE and PDE.