A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling
<p>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 mo...
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2023-07-01
|
| Series: | Geoscientific Model Development |
| Online Access: | https://gmd.copernicus.org/articles/16/3907/2023/gmd-16-3907-2023.pdf |
| _version_ | 1797781818843332608 |
|---|---|
| author | T. Richter V. Dansereau C. Lessig P. Minakowski |
| author_facet | T. Richter V. Dansereau C. Lessig P. Minakowski |
| author_sort | T. Richter |
| collection | DOAJ |
| description | <p>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.</p> |
| first_indexed | 2024-03-13T00:02:21Z |
| format | Article |
| id | doaj.art-d92a60fc0bdd4a08b2b5024634d396c7 |
| institution | Directory Open Access Journal |
| issn | 1991-959X 1991-9603 |
| language | English |
| last_indexed | 2024-03-13T00:02:21Z |
| publishDate | 2023-07-01 |
| publisher | Copernicus Publications |
| record_format | Article |
| series | Geoscientific Model Development |
| spelling | doaj.art-d92a60fc0bdd4a08b2b5024634d396c72023-07-13T06:06:10ZengCopernicus PublicationsGeoscientific Model Development1991-959X1991-96032023-07-01163907392610.5194/gmd-16-3907-2023A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modelingT. Richter0V. Dansereau1C. Lessig2P. Minakowski3Institute of Analysis and Numerics, Otto von Guericke University Magdeburg, Magdeburg, GermanyInstitut des Sciences de la Terre, Université Grenoble Alpes, CNRS (UMR5275), Gières, FranceInstitute of Simulation and Graphics, Otto von Guericke University Magdeburg, Magdeburg, GermanyInstitute of Analysis and Numerics, Otto von Guericke University Magdeburg, Magdeburg, Germany<p>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.</p>https://gmd.copernicus.org/articles/16/3907/2023/gmd-16-3907-2023.pdf |
| spellingShingle | T. Richter V. Dansereau C. Lessig P. Minakowski A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling Geoscientific Model Development |
| title | 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_short | 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 |
| url | https://gmd.copernicus.org/articles/16/3907/2023/gmd-16-3907-2023.pdf |
| work_keys_str_mv | AT trichter adynamicalcorebasedonadiscontinuousgalerkinmethodforhigherorderfiniteelementseaicemodeling AT vdansereau adynamicalcorebasedonadiscontinuousgalerkinmethodforhigherorderfiniteelementseaicemodeling AT clessig adynamicalcorebasedonadiscontinuousgalerkinmethodforhigherorderfiniteelementseaicemodeling AT pminakowski adynamicalcorebasedonadiscontinuousgalerkinmethodforhigherorderfiniteelementseaicemodeling AT trichter dynamicalcorebasedonadiscontinuousgalerkinmethodforhigherorderfiniteelementseaicemodeling AT vdansereau dynamicalcorebasedonadiscontinuousgalerkinmethodforhigherorderfiniteelementseaicemodeling AT clessig dynamicalcorebasedonadiscontinuousgalerkinmethodforhigherorderfiniteelementseaicemodeling AT pminakowski dynamicalcorebasedonadiscontinuousgalerkinmethodforhigherorderfiniteelementseaicemodeling |