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.
翻译:将离散速度(EVP)的计算方法可以是将离散速度(离散速度)放置在三角贝壳、细胞或边缘。由于细胞和边缘离散离异的特性比脊椎多,因此细胞和边缘离散的特性模拟出与脊离异的相同网格更多的线性运动特征。但是,基于细胞和边缘速度的离散作用受到压力率或压力差操作器内核的分解作用,需要如此处为细胞速度提议的特殊压力率计算,或如更早为边缘速度提议的稳定。一个简单的四重分析说明了内核是如何被清除的,还表明细胞和边缘速度的放置导致压力分解操作器的毛骨悚然的分流,产生很大的负向性叶值。尽管有与快速衰变和预期不会扭曲海洋冰动态的微分,但它们需要更小的内部时间步骤或更稳定的参数,以明确的EVP方法为准。