Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

IF value: 4.713
IF4.713
IF 5-year value: 4.927
IF 5-year
4.927
CiteScore value: 8.0
CiteScore
8.0
SNIP value: 1.425
SNIP1.425
IPP value: 4.65
IPP4.65
SJR value: 2.353
SJR2.353
Scimago H <br class='widget-line-break'>index value: 71
Scimago H
index
71
h5-index value: 53
h5-index53
Preprints
https://doi.org/10.5194/tc-2020-293
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/tc-2020-293
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

  19 Oct 2020

19 Oct 2020

Review status
This preprint is currently under review for the journal TC.

A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation

Niccolò Tubini1, Stephan Gruber2, and Riccardo Rigon1 Niccolò Tubini et al.
  • 1Department of Civil, Environmental and Mechanical Engineering, Universtity of Trento, Trento, Italy
  • 2Department of Geography and Environmental Studies, Carleton University, Ottawa, ON, K1S 5B6, Canada

Abstract. The accurate simulation of heat transfer with phase change is a central problem in cryosphere studies. This is because the nonlinear behaviour of enthalpy as function of temperature can prevent thermal models of snow, ice and frozen soil from converging to the correct solution. Existing numerical techniques rely on increased temporal resolution in trying to keep corresponding errors withing acceptable bounds. Here, we propose an algorithm, originally applied to solve water flow in soils, as a method to solve these integration issues with guaranteed convergence and conservation of energy for any time step size.

We review common modeling approaches, focusing on the fixed-grid method and on frozen soil. Based on this, we develop a conservative formulation of the governing equation and outline problems of alternative formulations in discretized form. Then, we apply the nested Newton-Casulli-Zanolli (NCZ) algorithm to a one-dimensional finite-volume discretization of the energy-enthalpy formulation.

Model performance is demonstrated against the Neumann and Lunardini analytical solutions and by comparing results from numerical experiments with integration time steps of one hour, one day, and ten days. Using our formulation and the NCZ algorithm, the convergence of the solver is guaranteed for any time step size. With this approach, the integration time step can be chosen to match the time scale of the processes investigated.

Niccolò Tubini et al.

Interactive discussion

Status: open (until 14 Dec 2020)
Status: open (until 14 Dec 2020)
AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment
[Subscribe to comment alert] Printer-friendly Version - Printer-friendly version Supplement - Supplement

Niccolò Tubini et al.

Data sets

OSF project Niccolò Tubini, Stephan Gruber, and Riccardo Rigon https://doi.org/10.17605/OSF.IO/8RDWS

Niccolò Tubini et al.

Viewed

Total article views: 345 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
239 102 4 345 2 3
  • HTML: 239
  • PDF: 102
  • XML: 4
  • Total: 345
  • BibTeX: 2
  • EndNote: 3
Views and downloads (calculated since 19 Oct 2020)
Cumulative views and downloads (calculated since 19 Oct 2020)

Viewed (geographical distribution)

Total article views: 307 (including HTML, PDF, and XML) Thereof 304 with geography defined and 3 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 

Cited

Saved

No saved metrics found.

Discussed

No discussed metrics found.
Latest update: 27 Nov 2020
Publications Copernicus
Download
Short summary
We present a new method to compute temperature changes with melting and freezing, a fundamental challenge in cryosphere research, extremely efficiently and with guaranteed correctness of the energy balance for any time step size. This is a key feature since the integration time step can then be chosen according to the timescale of the processes to study, from seconds to days.
We present a new method to compute temperature changes with melting and freezing, a fundamental...
Citation