В работе представлена математическая модель связанного тепломассопереноса с фазовыми переходами вода–лёд–пар в гетерогенных пористых грунтах. Модель основана на системе уравнений Ричардса, диффузии водяного пара и теплопереноса, связанных через температурно-зависимую влагопроводность, описываемую законом Козени–Кармана, и термодинамическое равновесие фаз, задаваемое соотношением Клапейрона–ван Генухтена. Целью исследования является выявление физического механизма подавления влагопереноса при промерзании и количественная характеристика резкого снижения его интенсивности. Получено аналитическое выражение для коэффициента диффузии почвенной влаги, определяющего скорость перераспределения влаги и связывающего влагопроводность грунта с его удельной влагоёмкостью. Показано, что поведение системы может быть описано в компактной безразмерной форме с управляющим параметром, разделяющим два предельных режима подавления влагопереноса: кинематический, обусловленный геометрической блокировкой порового пространства льдом, и термодинамический, связанный с фазовой инерцией системы вследствие латентной теплоты фазового перехода. Установлено, что при прохождении фронта замерзания происходит резкое падение коэффициента диффузии влаги, вызванный одновременным уменьшением влагопроводности и многократным ростом эффективной влагоёмкости среды. Определены характерный масштаб снижения коэффициента диффузии и критическая температура, соответствующая переходу между кинематическим и термодинамическим режимами подавления влагопереноса. Полученные результаты раскрывают физический механизм критического подавления влагопереноса в замерзающих пористых средах и могут быть использованы при разработке моделей тепловлагопереноса в мерзлых грунтах, прогнозировании процессов сезонного промерзания почв, а также при инженерных расчётах устойчивости грунтовых оснований и инфраструктуры в условиях холодного климата.

Статья подготовлена на основе англоязычного препринта, ранее опубликованного на платформе EarthArXiv (без процедуры рецензирования). Данная статья является его переработанной и существенно расширенной версией: дополнен раздел верификации модели (включая валидацию по натурным данным), проведено сопоставление с современными моделями криолитозоны, а также переработано обсуждение практических приложений.

A mathematical model of coupled heat and mass transfer with water–ice–vapor phase transitions in heterogeneous porous soils is presented. The model comprises the Richards equation, the water vapor diffusion equation, and the heat transfer equation, coupled through a temperature-dependent hydraulic conductivity governed by the Kozeny–Carman relation and thermodynamic phase equilibrium described by the Clapeyron–van Genuchten framework. The study aims to identify the physical mechanism responsible for the suppression of moisture transport during soil freezing and to quantitatively characterize the sharp reduction in the soil moisture diffusion coefficient. An analytical expression for the soil moisture diffusion coefficient is derived, linking the hydraulic conductivity to the specific moisture capacity and determining the rate of moisture redistribution in the soil. The analytical result admits a compact dimensionless form controlled by a single governing parameter that separates two physically distinct limiting regimes of moisture transport degradation: a kinematic regime dominated by geometric pore blockage due to ice formation, and a thermodynamic regime governed by phase inertia associated with the latent heat of the phase transition. Numerical simulations using the finite element method reveal a sharp collapse of the soil moisture diffusion coefficient upon passage of the freezing front, driven by the simultaneous reduction of hydraulic conductivity and a manifold increase in effective moisture capacity. The characteristic magnitude of the collapse and the critical temperature marking the transition between the two regimes are determined analytically. The~results elucidate the physical mechanism of critical moisture transport suppression in freezing porous media and provide a physically based parameterization applicable to thermo-hydrological models of frozen soils, prediction of seasonal freezing processes, and engineering assessment of soil foundation stability under cold-climate conditions.

The manuscript is based on an English-language preprint previously posted on EarthArXiv. The submitted Russian version is substantially revised, expanded, and translated. Full disclosure is provided in the manuscript.