Papers › Three-dimensional phase-field simulations of water freezing and thawing at pore-scale
Three-dimensional phase-field simulations of water freezing and thawing at pore-scale
Pavel Strachota
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
This work deals with numerical simulation of water freezing and thawing in a complex three-dimensional geometry of a porous medium. The porous structure is represented by a virtual container filled with glass beads. Phase transition modeling is approached at both macro-scale and micro-scale, combining heat transfer in a heterogeneous medium and a phase-field approximation of the Gibbs-Thomson relation by means of the Allen-Cahn equation. The formulation of the model contains novel components tailored for the given purpose. At the macro-scale, surface tension effects are negligible and phase transition focusing based on temperature can replace the Allen-Cahn equation. In contrast to that, simulations of equilibrium states at the micro-scale allow to eliminate the heat equation by assuming constant supercooling. For numerical solution, an efficient hybrid parallel algorithm based on the finite volume method and the Runge-Kutta-Merson solver with adaptive time stepping are employed. The results of different model variants at different scales are discussed. In a parametric study, the full phase-field model is demonstrated to deliver consistent results across a wide range of surface tension values, exhibiting curvature-induced premelting if surface tension is artificially exaggerated. As surface tension tends to the realistic values, the results of the phase-field approach those of the simplifed temperature-driven phase transition model. In addition, micro-scale simulations of water freezing at different supercooling values aim to predict the unfrozen water content and compare the results with data from literature. Numerical stability, accuracy, and computational costs are also discussed.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections