A regularization strategy for discontinuity when modelling coupled water and heat flow in freezing unsaturated soil

Result code in IS VaVaI

<a href="https://www.isvavai.cz/riv?ss=detail&h=RIV%2F60460709%3A41330%2F24%3A100795" target="_blank" >RIV/60460709:41330/24:100795 - isvavai.cz</a>

Result on the web

<a href="https://doi.org/10.1016/j.compfluid.2024.106299" target="_blank" >https://doi.org/10.1016/j.compfluid.2024.106299</a>

DOI - Digital Object Identifier

<a href="http://dx.doi.org/10.1016/j.compfluid.2024.106299" target="_blank" >10.1016/j.compfluid.2024.106299</a>

Result language

angličtina

Original language name

A regularization strategy for discontinuity when modelling coupled water and heat flow in freezing unsaturated soil

Original language description

Freezing tightly couples the water and heat flow. In most porous media, the interface between liquid and frozen water is not sharp and a slushy zone is present. There are two distinct approaches for mathematical modelling of freezing. It is the equilibrium approach which allows an instant freezing under given conditions and non -equilibrium approach where a specific timing in the freezing process is considered. In this contribution, we specifically target the equilibrium approach . The key mathematical model for the equilibrium approach is the Clausius-Clapeyron equation, which allows the derivation of a soil freezing curve relating temperature to pressure head. Implementing freezing soil accurately is not a straight -forward. Using the Clausius-Clapeyron equation creates a discontinuity in the freezing rate and latent heat at the freezing point. Little attention has been paid to the adequate description of the numerical treatment of this phenomenon and to the computational challenges that it poses. Numerical approximation of this discontinuous system is prone to spurious oscillations. In this contribution, we show the application of regularization of the discontinuous term. This treatment successfully stabilizes the computation and can remove oscillations. To avoid over -regularization, we present here a minimax strategy to determine optimal regularization parameters. We further compare an over -regularized setup with the non -equilibrium approach , where we show that a successful regularization is equivalent to incorporating a minimal timing to the freezing process. Finally - for validating our implementation and computational approach - experimental laboratory data of the volumetric water content and temperature profiles from previously published soil freezing experiments were represented here with the optimally regulized equilibrium approach .

Czech name

—

Czech description

—

Type

J_imp - Article in a specialist periodical, which is included in the Web of Science database

CEP classification

—

OECD FORD branch

10503 - Water resources

Project

<a href="/en/project/LTE119008" target="_blank" >LTE119008: AgriClima: Adding Value through Piloting of the EU-CELAC Climate Services Market in Agriculture</a><br>

Continuities

P - Projekt vyzkumu a vyvoje financovany z verejnych zdroju (s odkazem do CEP)<br>S - Specificky vyzkum na vysokych skolach

Publication year

2024

Confidentiality

S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů

Name of the periodical

COMPUTERS & FLUIDS

ISSN

0045-7930

e-ISSN

0045-7930

Volume of the periodical

277

Issue of the periodical within the volume

106299

Country of publishing house

CZ - CZECH REPUBLIC

Number of pages

13

Pages from-to

1-13

UT code for WoS article

001243399800001

EID of the result in the Scopus database

2-s2.0-85193294753