A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation
The accurate simulation of heat transfer with phase change is a central problem in cryosphere studies. This is because the non-linear 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...
Published in: | The Cryosphere |
---|---|
Main Authors: | , , |
Format: | Article in Journal/Newspaper |
Language: | English |
Published: |
Copernicus Publications
2021
|
Subjects: | |
Online Access: | https://doi.org/10.5194/tc-15-2541-2021 https://noa.gwlb.de/receive/cop_mods_00056863 https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00056513/tc-15-2541-2021.pdf https://tc.copernicus.org/articles/15/2541/2021/tc-15-2541-2021.pdf |
id |
ftnonlinearchiv:oai:noa.gwlb.de:cop_mods_00056863 |
---|---|
record_format |
openpolar |
spelling |
ftnonlinearchiv:oai:noa.gwlb.de:cop_mods_00056863 2023-05-15T18:32:32+02:00 A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation Tubini, Niccolò Gruber, Stephan Rigon, Riccardo 2021-06 electronic https://doi.org/10.5194/tc-15-2541-2021 https://noa.gwlb.de/receive/cop_mods_00056863 https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00056513/tc-15-2541-2021.pdf https://tc.copernicus.org/articles/15/2541/2021/tc-15-2541-2021.pdf eng eng Copernicus Publications The Cryosphere -- ˜Theœ Cryosphere -- http://www.bibliothek.uni-regensburg.de/ezeit/?2393169 -- http://www.the-cryosphere.net/ -- 1994-0424 https://doi.org/10.5194/tc-15-2541-2021 https://noa.gwlb.de/receive/cop_mods_00056863 https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00056513/tc-15-2541-2021.pdf https://tc.copernicus.org/articles/15/2541/2021/tc-15-2541-2021.pdf https://creativecommons.org/licenses/by/4.0/ uneingeschränkt info:eu-repo/semantics/openAccess CC-BY article Verlagsveröffentlichung article Text doc-type:article 2021 ftnonlinearchiv https://doi.org/10.5194/tc-15-2541-2021 2022-02-08T22:33:52Z The accurate simulation of heat transfer with phase change is a central problem in cryosphere studies. This is because the non-linear 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 within 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 modelling 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 1 h, 1 d, and 10 d. 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 timescale of the processes investigated. Article in Journal/Newspaper The Cryosphere Niedersächsisches Online-Archiv NOA The Cryosphere 15 6 2541 2568 |
institution |
Open Polar |
collection |
Niedersächsisches Online-Archiv NOA |
op_collection_id |
ftnonlinearchiv |
language |
English |
topic |
article Verlagsveröffentlichung |
spellingShingle |
article Verlagsveröffentlichung Tubini, Niccolò Gruber, Stephan Rigon, Riccardo A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation |
topic_facet |
article Verlagsveröffentlichung |
description |
The accurate simulation of heat transfer with phase change is a central problem in cryosphere studies. This is because the non-linear 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 within 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 modelling 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 1 h, 1 d, and 10 d. 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 timescale of the processes investigated. |
format |
Article in Journal/Newspaper |
author |
Tubini, Niccolò Gruber, Stephan Rigon, Riccardo |
author_facet |
Tubini, Niccolò Gruber, Stephan Rigon, Riccardo |
author_sort |
Tubini, Niccolò |
title |
A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation |
title_short |
A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation |
title_full |
A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation |
title_fullStr |
A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation |
title_full_unstemmed |
A method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation |
title_sort |
method for solving heat transfer with phase change in ice or soil that allows for large time steps while guaranteeing energy conservation |
publisher |
Copernicus Publications |
publishDate |
2021 |
url |
https://doi.org/10.5194/tc-15-2541-2021 https://noa.gwlb.de/receive/cop_mods_00056863 https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00056513/tc-15-2541-2021.pdf https://tc.copernicus.org/articles/15/2541/2021/tc-15-2541-2021.pdf |
genre |
The Cryosphere |
genre_facet |
The Cryosphere |
op_relation |
The Cryosphere -- ˜Theœ Cryosphere -- http://www.bibliothek.uni-regensburg.de/ezeit/?2393169 -- http://www.the-cryosphere.net/ -- 1994-0424 https://doi.org/10.5194/tc-15-2541-2021 https://noa.gwlb.de/receive/cop_mods_00056863 https://noa.gwlb.de/servlets/MCRFileNodeServlet/cop_derivate_00056513/tc-15-2541-2021.pdf https://tc.copernicus.org/articles/15/2541/2021/tc-15-2541-2021.pdf |
op_rights |
https://creativecommons.org/licenses/by/4.0/ uneingeschränkt info:eu-repo/semantics/openAccess |
op_rightsnorm |
CC-BY |
op_doi |
https://doi.org/10.5194/tc-15-2541-2021 |
container_title |
The Cryosphere |
container_volume |
15 |
container_issue |
6 |
container_start_page |
2541 |
op_container_end_page |
2568 |
_version_ |
1766216683165319168 |