Improving numerical accuracy for the viscous-plastic formulation of sea ice
Accurate modeling of sea ice dynamics is critical for predicting environmental variables and is important in applications such as navigating ice breaker ships. Research for both modeling and simulating sea ice dynamics is ongoing, with the most widely accepted model based on the viscous-plastic (VP)...
Published in: | Journal of Computational Physics |
---|---|
Main Authors: | , , |
Format: | Text |
Language: | unknown |
Published: |
2022
|
Subjects: | |
Online Access: | http://arxiv.org/abs/2206.10061 https://doi.org/10.1016/j.jcp.2023.112184 |
id |
ftarxivpreprints:oai:arXiv.org:2206.10061 |
---|---|
record_format |
openpolar |
spelling |
ftarxivpreprints:oai:arXiv.org:2206.10061 2023-09-05T13:23:00+02:00 Improving numerical accuracy for the viscous-plastic formulation of sea ice Li, Tongtong Gelb, Anne Lee, Yoonsang 2022-06-20 http://arxiv.org/abs/2206.10061 https://doi.org/10.1016/j.jcp.2023.112184 unknown http://arxiv.org/abs/2206.10061 doi:10.1016/j.jcp.2023.112184 Mathematics - Numerical Analysis Physics - Atmospheric and Oceanic Physics text 2022 ftarxivpreprints https://doi.org/10.1016/j.jcp.2023.112184 2023-08-16T17:08:32Z Accurate modeling of sea ice dynamics is critical for predicting environmental variables and is important in applications such as navigating ice breaker ships. Research for both modeling and simulating sea ice dynamics is ongoing, with the most widely accepted model based on the viscous-plastic (VP) formulation introduced by Hibler in 1979. Due to its highly nonlinear features, this model is intrinsically challenging for computational solvers. In particular, sea ice simulations often significantly differ from satellite observations. This study therefore focuses on improving the numerical accuracy of the VP sea ice model. Since the poor convergence observed in existing numerical simulations stems from the nonlinear nature of the VP formulation, this investigation proposes using the celebrated weighted essentially non-oscillatory (WENO) scheme -- as opposed to the frequently employed centered difference (CD) scheme -- for the spatial derivatives in the VP sea ice model. We then proceed to numerically demonstrate that WENO yields higher-order convergence for smooth solutions, and that furthermore it is able to resolve the discontinuities in the sharp features of sea ice covers -- something that is not possible using CD methods. Finally, our proposed framework integrates a potential function method that utilizes the phase field method to naturally incorporates the physical restrictions of ice thickness and ice concentration in transport equations, resulting in a modified transport equations which includes additional forcing terms. Our method does not require post-processing, thereby avoiding the possible introduction of discontinuities and corresponding negative impact on the solution behavior. Numerical experiments are provided to demonstrate the efficacy of our new methodology. Text Sea ice ArXiv.org (Cornell University Library) Breaker ENVELOPE(-67.257,-67.257,-67.874,-67.874) Journal of Computational Physics 487 112184 |
institution |
Open Polar |
collection |
ArXiv.org (Cornell University Library) |
op_collection_id |
ftarxivpreprints |
language |
unknown |
topic |
Mathematics - Numerical Analysis Physics - Atmospheric and Oceanic Physics |
spellingShingle |
Mathematics - Numerical Analysis Physics - Atmospheric and Oceanic Physics Li, Tongtong Gelb, Anne Lee, Yoonsang Improving numerical accuracy for the viscous-plastic formulation of sea ice |
topic_facet |
Mathematics - Numerical Analysis Physics - Atmospheric and Oceanic Physics |
description |
Accurate modeling of sea ice dynamics is critical for predicting environmental variables and is important in applications such as navigating ice breaker ships. Research for both modeling and simulating sea ice dynamics is ongoing, with the most widely accepted model based on the viscous-plastic (VP) formulation introduced by Hibler in 1979. Due to its highly nonlinear features, this model is intrinsically challenging for computational solvers. In particular, sea ice simulations often significantly differ from satellite observations. This study therefore focuses on improving the numerical accuracy of the VP sea ice model. Since the poor convergence observed in existing numerical simulations stems from the nonlinear nature of the VP formulation, this investigation proposes using the celebrated weighted essentially non-oscillatory (WENO) scheme -- as opposed to the frequently employed centered difference (CD) scheme -- for the spatial derivatives in the VP sea ice model. We then proceed to numerically demonstrate that WENO yields higher-order convergence for smooth solutions, and that furthermore it is able to resolve the discontinuities in the sharp features of sea ice covers -- something that is not possible using CD methods. Finally, our proposed framework integrates a potential function method that utilizes the phase field method to naturally incorporates the physical restrictions of ice thickness and ice concentration in transport equations, resulting in a modified transport equations which includes additional forcing terms. Our method does not require post-processing, thereby avoiding the possible introduction of discontinuities and corresponding negative impact on the solution behavior. Numerical experiments are provided to demonstrate the efficacy of our new methodology. |
format |
Text |
author |
Li, Tongtong Gelb, Anne Lee, Yoonsang |
author_facet |
Li, Tongtong Gelb, Anne Lee, Yoonsang |
author_sort |
Li, Tongtong |
title |
Improving numerical accuracy for the viscous-plastic formulation of sea ice |
title_short |
Improving numerical accuracy for the viscous-plastic formulation of sea ice |
title_full |
Improving numerical accuracy for the viscous-plastic formulation of sea ice |
title_fullStr |
Improving numerical accuracy for the viscous-plastic formulation of sea ice |
title_full_unstemmed |
Improving numerical accuracy for the viscous-plastic formulation of sea ice |
title_sort |
improving numerical accuracy for the viscous-plastic formulation of sea ice |
publishDate |
2022 |
url |
http://arxiv.org/abs/2206.10061 https://doi.org/10.1016/j.jcp.2023.112184 |
long_lat |
ENVELOPE(-67.257,-67.257,-67.874,-67.874) |
geographic |
Breaker |
geographic_facet |
Breaker |
genre |
Sea ice |
genre_facet |
Sea ice |
op_relation |
http://arxiv.org/abs/2206.10061 doi:10.1016/j.jcp.2023.112184 |
op_doi |
https://doi.org/10.1016/j.jcp.2023.112184 |
container_title |
Journal of Computational Physics |
container_volume |
487 |
container_start_page |
112184 |
_version_ |
1776203573320744960 |