On 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 a viscous-plastic (VP) ice rheology are difficult to solve numerically, because the associated bulk and shear viscosities can be very large. Traditionally, implicit solution techniques for the VP rheology are thought to be expensive; the explici...
Published in: | Ocean Modelling |
---|---|
Main Authors: | , |
Format: | Article in Journal/Newspaper |
Language: | unknown |
Published: |
2012
|
Subjects: | |
Online Access: | https://epic.awi.de/id/eprint/20437/ https://epic.awi.de/id/eprint/20437/4/icetest_resubmitted2011.pdf https://hdl.handle.net/10013/epic.38521 https://hdl.handle.net/10013/epic.38521.d004 |
id |
ftawi:oai:epic.awi.de:20437 |
---|---|
record_format |
openpolar |
spelling |
ftawi:oai:epic.awi.de:20437 2024-09-15T18:35:21+00:00 On Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique Losch, Martin Danilov, Sergey 2012 application/pdf https://epic.awi.de/id/eprint/20437/ https://epic.awi.de/id/eprint/20437/4/icetest_resubmitted2011.pdf https://hdl.handle.net/10013/epic.38521 https://hdl.handle.net/10013/epic.38521.d004 unknown https://epic.awi.de/id/eprint/20437/4/icetest_resubmitted2011.pdf https://hdl.handle.net/10013/epic.38521.d004 Losch, M. orcid:0000-0002-3824-5244 and Danilov, S. orcid:0000-0001-8098-182X (2012) On Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique , Ocean Modelling, 41 , pp. 42-52 . doi:10.1016/j.ocemod.2011.10.002 <https://doi.org/10.1016/j.ocemod.2011.10.002> , hdl:10013/epic.38521 EPIC3Ocean Modelling, 41, pp. 42-52 Article isiRev 2012 ftawi https://doi.org/10.1016/j.ocemod.2011.10.002 2024-06-24T04:01:33Z The momentum equations that describe sea ice drift for a viscous-plastic (VP) ice rheology are difficult to solve numerically, because the associated bulk and shear viscosities can be very large. Traditionally, implicit solution techniques for the VP rheology are thought to be expensive; the explicit elastic-viscous-plastic (EVP) method was designed to be more efficient and accurate. In order to assess their relative performance, experiments with idealized geometry are used to compare model solutions of implicit VP- and explicit EVP-solvers in two very different ice-ocean codes: the regular-grid, finite-volume Massachusetts Institute of Technology general circulation model (MITgcm) and the Alfred Wegener Institute Finite Element Ocean Model (FEOM). For both codes the obtained solutions 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 VP solutions for very small sub-cycling time steps. A limiting scheme for EVP viscosities, that addresses a noise problem, reduces the viscosity even further and, especially in the case of the variable-resolution unstructured grids of FEOM, can lead to unexpected ice distributions that are dramatically different from solutions without this scheme. Implicit VP-solvers are found to be generally faster than the EVP-solvers, most likely because the ice distribution does not change much within the short time steps of this study. Short time steps are thought to be typical of present day high resolution ice-ocean models, so that previous timing results for long time steps may no longer be representative. Article in Journal/Newspaper Sea ice Alfred Wegener Institute for Polar- and Marine Research (AWI): ePIC (electronic Publication Information Center) Ocean Modelling 41 42 52 |
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 a viscous-plastic (VP) ice rheology are difficult to solve numerically, because the associated bulk and shear viscosities can be very large. Traditionally, implicit solution techniques for the VP rheology are thought to be expensive; the explicit elastic-viscous-plastic (EVP) method was designed to be more efficient and accurate. In order to assess their relative performance, experiments with idealized geometry are used to compare model solutions of implicit VP- and explicit EVP-solvers in two very different ice-ocean codes: the regular-grid, finite-volume Massachusetts Institute of Technology general circulation model (MITgcm) and the Alfred Wegener Institute Finite Element Ocean Model (FEOM). For both codes the obtained solutions 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 VP solutions for very small sub-cycling time steps. A limiting scheme for EVP viscosities, that addresses a noise problem, reduces the viscosity even further and, especially in the case of the variable-resolution unstructured grids of FEOM, can lead to unexpected ice distributions that are dramatically different from solutions without this scheme. Implicit VP-solvers are found to be generally faster than the EVP-solvers, most likely because the ice distribution does not change much within the short time steps of this study. Short time steps are thought to be typical of present day high resolution ice-ocean models, so that previous timing results for long time steps may no longer be representative. |
format |
Article in Journal/Newspaper |
author |
Losch, Martin Danilov, Sergey |
spellingShingle |
Losch, Martin Danilov, Sergey On 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 |
On Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_short |
On Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_full |
On Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_fullStr |
On Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_full_unstemmed |
On Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique |
title_sort |
on solving the momentum equations of dynamic sea ice models with implicit solvers and the elastic-viscous-plastic technique |
publishDate |
2012 |
url |
https://epic.awi.de/id/eprint/20437/ https://epic.awi.de/id/eprint/20437/4/icetest_resubmitted2011.pdf https://hdl.handle.net/10013/epic.38521 https://hdl.handle.net/10013/epic.38521.d004 |
genre |
Sea ice |
genre_facet |
Sea ice |
op_source |
EPIC3Ocean Modelling, 41, pp. 42-52 |
op_relation |
https://epic.awi.de/id/eprint/20437/4/icetest_resubmitted2011.pdf https://hdl.handle.net/10013/epic.38521.d004 Losch, M. orcid:0000-0002-3824-5244 and Danilov, S. orcid:0000-0001-8098-182X (2012) On Solving the Momentum Equations of Dynamic Sea Ice Models with Implicit Solvers and the Elastic-Viscous-Plastic Technique , Ocean Modelling, 41 , pp. 42-52 . doi:10.1016/j.ocemod.2011.10.002 <https://doi.org/10.1016/j.ocemod.2011.10.002> , hdl:10013/epic.38521 |
op_doi |
https://doi.org/10.1016/j.ocemod.2011.10.002 |
container_title |
Ocean Modelling |
container_volume |
41 |
container_start_page |
42 |
op_container_end_page |
52 |
_version_ |
1810478481165778944 |