On discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities
Discretization of the equations of Viscous Plastic and Elastic Viscous Plastic (EVP) sea ice dynamics on triangular meshes can be done by placing discrete velocities at vertices, cells or edges. Since there are more cells and edges than vertices, the cell- and edge-based discretizations simulate mor...
Published in: | Ocean Modelling |
---|---|
Main Authors: | , , |
Format: | Text |
Language: | unknown |
Published: |
2021
|
Subjects: | |
Online Access: | http://arxiv.org/abs/2106.13641 https://doi.org/10.1016/j.ocemod.2021.101937 |
id |
ftarxivpreprints:oai:arXiv.org:2106.13641 |
---|---|
record_format |
openpolar |
spelling |
ftarxivpreprints:oai:arXiv.org:2106.13641 2023-09-05T13:23:02+02:00 On discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities Danilov, S. Mehlmann, C. Fofonova, V. 2021-06-25 http://arxiv.org/abs/2106.13641 https://doi.org/10.1016/j.ocemod.2021.101937 unknown http://arxiv.org/abs/2106.13641 doi:10.1016/j.ocemod.2021.101937 Mathematics - Numerical Analysis text 2021 ftarxivpreprints https://doi.org/10.1016/j.ocemod.2021.101937 2023-08-16T16:33:42Z Discretization of the equations of Viscous Plastic and Elastic Viscous Plastic (EVP) sea ice dynamics on triangular meshes can be done by placing discrete velocities at vertices, cells or edges. Since there are more cells and edges than vertices, the cell- and edge-based discretizations simulate more linear kinematic features at the same mesh than the vertex discretization. However, the discretization based on cell and edge velocities suffer from kernels in the strain rate or stress divergence operators and need either special strain rate computations as proposed here for cell velocities, or stabilization as proposed earlier for edge velocities. An elementary Fourier analysis clarifies how kernels are removed, and also shows that cell and edge velocity placement leads to spurious branches of stress divergence operator with large negative eigenvalues. Although spurious branches correspond to fast decay and are not expected to distort sea ice dynamics, they demand either smaller internal time steps or higher stability parameters in explicit EVP-like methods. Text Sea ice ArXiv.org (Cornell University Library) Ocean Modelling 170 101937 |
institution |
Open Polar |
collection |
ArXiv.org (Cornell University Library) |
op_collection_id |
ftarxivpreprints |
language |
unknown |
topic |
Mathematics - Numerical Analysis |
spellingShingle |
Mathematics - Numerical Analysis Danilov, S. Mehlmann, C. Fofonova, V. On discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities |
topic_facet |
Mathematics - Numerical Analysis |
description |
Discretization of the equations of Viscous Plastic and Elastic Viscous Plastic (EVP) sea ice dynamics on triangular meshes can be done by placing discrete velocities at vertices, cells or edges. Since there are more cells and edges than vertices, the cell- and edge-based discretizations simulate more linear kinematic features at the same mesh than the vertex discretization. However, the discretization based on cell and edge velocities suffer from kernels in the strain rate or stress divergence operators and need either special strain rate computations as proposed here for cell velocities, or stabilization as proposed earlier for edge velocities. An elementary Fourier analysis clarifies how kernels are removed, and also shows that cell and edge velocity placement leads to spurious branches of stress divergence operator with large negative eigenvalues. Although spurious branches correspond to fast decay and are not expected to distort sea ice dynamics, they demand either smaller internal time steps or higher stability parameters in explicit EVP-like methods. |
format |
Text |
author |
Danilov, S. Mehlmann, C. Fofonova, V. |
author_facet |
Danilov, S. Mehlmann, C. Fofonova, V. |
author_sort |
Danilov, S. |
title |
On discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities |
title_short |
On discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities |
title_full |
On discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities |
title_fullStr |
On discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities |
title_full_unstemmed |
On discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities |
title_sort |
on discretizing sea-ice dynamics on triangular meshes using vertex, cell or edge velocities |
publishDate |
2021 |
url |
http://arxiv.org/abs/2106.13641 https://doi.org/10.1016/j.ocemod.2021.101937 |
genre |
Sea ice |
genre_facet |
Sea ice |
op_relation |
http://arxiv.org/abs/2106.13641 doi:10.1016/j.ocemod.2021.101937 |
op_doi |
https://doi.org/10.1016/j.ocemod.2021.101937 |
container_title |
Ocean Modelling |
container_volume |
170 |
container_start_page |
101937 |
_version_ |
1776203619316531200 |