A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling
The ability of numerical sea ice models to reproduce localized deformation features associated with fracture processes is key for an accurate representation of the ice dynamics and of dynamically coupled physical processes in the Arctic and Antarctic. Equally key is the capacity of these models to m...
Published in: | Geoscientific Model Development |
---|---|
Main Authors: | , , , |
Format: | Text |
Language: | English |
Published: |
2023
|
Subjects: | |
Online Access: | https://doi.org/10.5194/gmd-16-3907-2023 https://gmd.copernicus.org/articles/16/3907/2023/ |
id |
ftcopernicus:oai:publications.copernicus.org:gmd109989 |
---|---|
record_format |
openpolar |
spelling |
ftcopernicus:oai:publications.copernicus.org:gmd109989 2023-07-30T03:58:19+02:00 A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling Richter, Thomas Dansereau, Véronique Lessig, Christian Minakowski, Piotr 2023-07-13 application/pdf https://doi.org/10.5194/gmd-16-3907-2023 https://gmd.copernicus.org/articles/16/3907/2023/ eng eng doi:10.5194/gmd-16-3907-2023 https://gmd.copernicus.org/articles/16/3907/2023/ eISSN: 1991-9603 Text 2023 ftcopernicus https://doi.org/10.5194/gmd-16-3907-2023 2023-07-17T16:24:17Z The ability of numerical sea ice models to reproduce localized deformation features associated with fracture processes is key for an accurate representation of the ice dynamics and of dynamically coupled physical processes in the Arctic and Antarctic. Equally key is the capacity of these models to minimize the numerical diffusion stemming from the advection of these features to ensure that the associated strong gradients persist in time, without the need to unphysically re-inject energy for re-localization. To control diffusion and improve the approximation quality, we present a new numerical core for the dynamics of sea ice that is based on higher-order finite-element discretizations for the momentum equation and higher-order discontinuous Galerkin methods for the advection. The mathematical properties of this core are discussed, and a detailed description of an efficient shared-memory parallel implementation is given. In addition, we present different numerical tests and apply the new framework to a benchmark problem to quantify the advantages of the higher-order discretization. These tests are based on Hibler's viscous–plastic sea ice model, but the implementation of the developed framework in the context of other physical models reproducing a strong localization of the deformation is possible. Text Antarc* Antarctic Arctic Sea ice Copernicus Publications: E-Journals Arctic Antarctic Geoscientific Model Development 16 13 3907 3926 |
institution |
Open Polar |
collection |
Copernicus Publications: E-Journals |
op_collection_id |
ftcopernicus |
language |
English |
description |
The ability of numerical sea ice models to reproduce localized deformation features associated with fracture processes is key for an accurate representation of the ice dynamics and of dynamically coupled physical processes in the Arctic and Antarctic. Equally key is the capacity of these models to minimize the numerical diffusion stemming from the advection of these features to ensure that the associated strong gradients persist in time, without the need to unphysically re-inject energy for re-localization. To control diffusion and improve the approximation quality, we present a new numerical core for the dynamics of sea ice that is based on higher-order finite-element discretizations for the momentum equation and higher-order discontinuous Galerkin methods for the advection. The mathematical properties of this core are discussed, and a detailed description of an efficient shared-memory parallel implementation is given. In addition, we present different numerical tests and apply the new framework to a benchmark problem to quantify the advantages of the higher-order discretization. These tests are based on Hibler's viscous–plastic sea ice model, but the implementation of the developed framework in the context of other physical models reproducing a strong localization of the deformation is possible. |
format |
Text |
author |
Richter, Thomas Dansereau, Véronique Lessig, Christian Minakowski, Piotr |
spellingShingle |
Richter, Thomas Dansereau, Véronique Lessig, Christian Minakowski, Piotr A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling |
author_facet |
Richter, Thomas Dansereau, Véronique Lessig, Christian Minakowski, Piotr |
author_sort |
Richter, Thomas |
title |
A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling |
title_short |
A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling |
title_full |
A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling |
title_fullStr |
A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling |
title_full_unstemmed |
A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling |
title_sort |
dynamical core based on a discontinuous galerkin method for higher-order finite-element sea ice modeling |
publishDate |
2023 |
url |
https://doi.org/10.5194/gmd-16-3907-2023 https://gmd.copernicus.org/articles/16/3907/2023/ |
geographic |
Arctic Antarctic |
geographic_facet |
Arctic Antarctic |
genre |
Antarc* Antarctic Arctic Sea ice |
genre_facet |
Antarc* Antarctic Arctic Sea ice |
op_source |
eISSN: 1991-9603 |
op_relation |
doi:10.5194/gmd-16-3907-2023 https://gmd.copernicus.org/articles/16/3907/2023/ |
op_doi |
https://doi.org/10.5194/gmd-16-3907-2023 |
container_title |
Geoscientific Model Development |
container_volume |
16 |
container_issue |
13 |
container_start_page |
3907 |
op_container_end_page |
3926 |
_version_ |
1772821149427695616 |