Journal of Computational Physics ( IF 3.8 ) Pub Date : 2020-12-29 , DOI: 10.1016/j.jcp.2020.110086 Carolin Mehlmann , Peter Korn
We present a discretization of the dynamics of sea-ice on triangular grids. Our numerical approach is based on the nonconforming Crouzeix-Raviart finite element. An advantage of this element is that it facilitates the coupling to an ocean model that employs an Arakawa C-type staggering of variables. We show that the Crouzeix-Raviart element implements a discretization of the viscous-plastic and elastic-viscous-plastic stress tensor that suffers from unacceptable small scale noise in the velocity field. To resolve this issue we introduce an edge-based stabilization of the Crouzeix-Raviart element. Through a blend of theoretical considerations, based on the Korn inequality, and numerical experiments we show that the stabilized Crouzeix-Raviart element provides a stable discretization of sea-ice dynamics on triangular grids that is relevant for sea-ice modelling in ocean and climate science.
中文翻译:
三角网格上的海冰动力学
我们提出了三角网格上海冰动力学的离散化。我们的数值方法基于不合格的Crouzeix-Raviart有限元。此元素的一个优点是,它有助于与使用Arakawa C型变量交错的海洋模型进行耦合。我们表明,Crouzeix-Raviart元素实现了粘塑性和弹性粘塑性应力张量的离散化,该应力张量在速度场中遭受不可接受的小规模噪声。为了解决此问题,我们引入了Crouzeix-Raviart元素的基于边缘的稳定性。结合基于科恩不等式的理论考虑,