Solution of the one-dimensional Stefan problem with two transitions for modelling of the water freezing in a glacial crevasse

This article presents a numerical solution of the one-dimensional Stefan problem with two phase transitions, which is implemented on a non-uniform grid. The system of equations is written in a general form, i.e. it includes not only conductive, but also convective and dissipative terms. The problem...

Full description

Saved in:
Bibliographic Details
Main Author: S. V. Popov
Format: Article
Language:Russian
Published: Nauka 2023-04-01
Series:Лëд и снег
Subjects:
Online Access:https://ice-snow.igras.ru/jour/article/view/1157
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1850025762101395456
author S. V. Popov
author_facet S. V. Popov
author_sort S. V. Popov
collection DOAJ
description This article presents a numerical solution of the one-dimensional Stefan problem with two phase transitions, which is implemented on a non-uniform grid. The system of equations is written in a general form, i.e. it includes not only conductive, but also convective and dissipative terms. The problem is solved numerically by the front-fixing method on a non-uniform grid using an implicit finite-difference scheme, which is implemented by the sweep method. This algorithm can also be used to create more complex mathematical models of heat and mass transfer, as well as to describe glacial and subglacial processes. The mathematical apparatus proposed in the article was used to solve a specific problem of water freezing in a glacial crevasse. The presence and progression of crevasses, in turn, is a demonstrative factor indicating the dynamic activity of the glacier. Crevasses formed in one way or another can not only expand, but also decrease in size until they completely disappear. One of the reasons for their closure is the freezing of near-surface meltwater in the crevasse. Such a process was observed on glaciers near Mirny and Novolazarevskaya stations (East Antarctica). This process is modeled as an example of solving the Stefan problem. It is believed that all media are homogeneous and isotropic. The temperature of the water in the crevasse corresponds to the melting temperature of the ice. Modeling has shown that for the coastal part of the cold Antarctic glacier with an average temperature of –10°C and below, crevasses 5–10 cm of width freeze in less than a week. Wider ones freeze a little longer. 30 cm wide crevasses close in about two to three weeks, depending on the temperature of the glacier.
format Article
id doaj-art-ab5251f5cd67455ca063beabf1c23ba6
institution DOAJ
issn 2076-6734
2412-3765
language Russian
publishDate 2023-04-01
publisher Nauka
record_format Article
series Лëд и снег
spelling doaj-art-ab5251f5cd67455ca063beabf1c23ba62025-08-20T03:00:45ZrusNaukaЛëд и снег2076-67342412-37652023-04-0163113014010.31857/S2076673423010131789Solution of the one-dimensional Stefan problem with two transitions for modelling of the water freezing in a glacial crevasseS. V. Popov0Polar Marine Geosurvey Expedition; Saint-Petersburg State University; Melnikov Permafrost InstituteThis article presents a numerical solution of the one-dimensional Stefan problem with two phase transitions, which is implemented on a non-uniform grid. The system of equations is written in a general form, i.e. it includes not only conductive, but also convective and dissipative terms. The problem is solved numerically by the front-fixing method on a non-uniform grid using an implicit finite-difference scheme, which is implemented by the sweep method. This algorithm can also be used to create more complex mathematical models of heat and mass transfer, as well as to describe glacial and subglacial processes. The mathematical apparatus proposed in the article was used to solve a specific problem of water freezing in a glacial crevasse. The presence and progression of crevasses, in turn, is a demonstrative factor indicating the dynamic activity of the glacier. Crevasses formed in one way or another can not only expand, but also decrease in size until they completely disappear. One of the reasons for their closure is the freezing of near-surface meltwater in the crevasse. Such a process was observed on glaciers near Mirny and Novolazarevskaya stations (East Antarctica). This process is modeled as an example of solving the Stefan problem. It is believed that all media are homogeneous and isotropic. The temperature of the water in the crevasse corresponds to the melting temperature of the ice. Modeling has shown that for the coastal part of the cold Antarctic glacier with an average temperature of –10°C and below, crevasses 5–10 cm of width freeze in less than a week. Wider ones freeze a little longer. 30 cm wide crevasses close in about two to three weeks, depending on the temperature of the glacier.https://ice-snow.igras.ru/jour/article/view/1157mathematical modellingnumerical solutionstefan problemfinite difference schemesice crevassesantarctica
spellingShingle S. V. Popov
Solution of the one-dimensional Stefan problem with two transitions for modelling of the water freezing in a glacial crevasse
Лëд и снег
mathematical modelling
numerical solution
stefan problem
finite difference schemes
ice crevasses
antarctica
title Solution of the one-dimensional Stefan problem with two transitions for modelling of the water freezing in a glacial crevasse
title_full Solution of the one-dimensional Stefan problem with two transitions for modelling of the water freezing in a glacial crevasse
title_fullStr Solution of the one-dimensional Stefan problem with two transitions for modelling of the water freezing in a glacial crevasse
title_full_unstemmed Solution of the one-dimensional Stefan problem with two transitions for modelling of the water freezing in a glacial crevasse
title_short Solution of the one-dimensional Stefan problem with two transitions for modelling of the water freezing in a glacial crevasse
title_sort solution of the one dimensional stefan problem with two transitions for modelling of the water freezing in a glacial crevasse
topic mathematical modelling
numerical solution
stefan problem
finite difference schemes
ice crevasses
antarctica
url https://ice-snow.igras.ru/jour/article/view/1157
work_keys_str_mv AT svpopov solutionoftheonedimensionalstefanproblemwithtwotransitionsformodellingofthewaterfreezinginaglacialcrevasse