Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique
The momentum equations that describe sea ice drift for aviscous-plastic (VP) ice rheology are difficult to solve numerically,because the associated bulk and shear viscosities can be verylarge. Traditionally, implicit solution techniques for the VP rheologyare thought to be expensive; the explicit el...
Main Authors: | , |
---|---|
Format: | Conference Object |
Language: | unknown |
Published: |
2009
|
Subjects: | |
Online Access: | https://epic.awi.de/id/eprint/20436/ https://hdl.handle.net/10013/epic.32562 |
id |
ftawi:oai:epic.awi.de:20436 |
---|---|
record_format |
openpolar |
spelling |
ftawi:oai:epic.awi.de:20436 2023-05-15T18:17:59+02:00 Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique Losch, Martin Danilov, Sergey 2009 https://epic.awi.de/id/eprint/20436/ https://hdl.handle.net/10013/epic.32562 unknown Losch, M. orcid:0000-0002-3824-5244 and Danilov, S. orcid:0000-0001-8098-182X (2009) Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique , EGU General Assembly, Vienna, April . hdl:10013/epic.32562 EPIC3EGU General Assembly, Vienna, April, pp. 19-24 Conference notRev 2009 ftawi 2021-12-24T15:33:25Z The momentum equations that describe sea ice drift for aviscous-plastic (VP) ice rheology are difficult to solve numerically,because the associated bulk and shear viscosities can be verylarge. Traditionally, implicit solution techniques for the VP rheologyare thought to be expensive; the explicit elastic-viscous-plastic(EVP) method was designed to be more efficient and accurate. In orderto assess their relative performance, experiments with idealizedgeometry are used to compare model solutions of implicit VP- andexplicit EVP-solvers in two very different ice-ocean codes: theregular-grid, finite-volume Massachusetts Institute of Technologygeneral circulation model (MITgcm) and the Alfred Wegener InstituteFinite Element Ocean Model (FEOM). For both codes the obtainedsolutions of implicit VP- and EVP-solvers can differ significantly,because the EVP solutions tend to have smaller ice viscosities(weaker ice). EVP solutions tend to converge to implicit VPsolutions for very small sub-cycling time steps. A limiting scheme forEVP viscosities, that addresses a noise problem, reduces the viscosityeven further and, especially in the case of the variable-resolutionunstructured grids of FEOM, can lead to unexpected ice distributionsthat are dramatically different from solutions without thisscheme. Implicit VP-solvers are found to be generally faster than theEVP-solvers, most likely because the ice distribution does not changemuch within the short time steps of this study. Short time steps arethought to be typical of present day high resolution ice-ocean models,so that previous timing results for long time steps may no longer berepresentative. Conference Object Sea ice Alfred Wegener Institute for Polar- and Marine Research (AWI): ePIC (electronic Publication Information Center) |
institution |
Open Polar |
collection |
Alfred Wegener Institute for Polar- and Marine Research (AWI): ePIC (electronic Publication Information Center) |
op_collection_id |
ftawi |
language |
unknown |
description |
The momentum equations that describe sea ice drift for aviscous-plastic (VP) ice rheology are difficult to solve numerically,because the associated bulk and shear viscosities can be verylarge. Traditionally, implicit solution techniques for the VP rheologyare thought to be expensive; the explicit elastic-viscous-plastic(EVP) method was designed to be more efficient and accurate. In orderto assess their relative performance, experiments with idealizedgeometry are used to compare model solutions of implicit VP- andexplicit EVP-solvers in two very different ice-ocean codes: theregular-grid, finite-volume Massachusetts Institute of Technologygeneral circulation model (MITgcm) and the Alfred Wegener InstituteFinite Element Ocean Model (FEOM). For both codes the obtainedsolutions of implicit VP- and EVP-solvers can differ significantly,because the EVP solutions tend to have smaller ice viscosities(weaker ice). EVP solutions tend to converge to implicit VPsolutions for very small sub-cycling time steps. A limiting scheme forEVP viscosities, that addresses a noise problem, reduces the viscosityeven further and, especially in the case of the variable-resolutionunstructured grids of FEOM, can lead to unexpected ice distributionsthat are dramatically different from solutions without thisscheme. Implicit VP-solvers are found to be generally faster than theEVP-solvers, most likely because the ice distribution does not changemuch within the short time steps of this study. Short time steps arethought to be typical of present day high resolution ice-ocean models,so that previous timing results for long time steps may no longer berepresentative. |
format |
Conference Object |
author |
Losch, Martin Danilov, Sergey |
spellingShingle |
Losch, Martin Danilov, Sergey Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
author_facet |
Losch, Martin Danilov, Sergey |
author_sort |
Losch, Martin |
title |
Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_short |
Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_full |
Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_fullStr |
Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_full_unstemmed |
Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_sort |
solving the momentum equations of dynamic sea ice models with implicit solvers and the elastic-viscous-plastic technique |
publishDate |
2009 |
url |
https://epic.awi.de/id/eprint/20436/ https://hdl.handle.net/10013/epic.32562 |
genre |
Sea ice |
genre_facet |
Sea ice |
op_source |
EPIC3EGU General Assembly, Vienna, April, pp. 19-24 |
op_relation |
Losch, M. orcid:0000-0002-3824-5244 and Danilov, S. orcid:0000-0001-8098-182X (2009) Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique , EGU General Assembly, Vienna, April . hdl:10013/epic.32562 |
_version_ |
1766193839304867840 |