Abstract:In the study of sea ice rheological behavior under different temporal and spatial scales, a series of numerical methods have been developed in the past several decades. Nowadays, there are mainly four methods applied commonly, which are Finite Different (FD) method, Particle-In-Cell (PIC), Smoothed Particle Hydrodynamics (SPH) and granular flow (GF) method. The Eulerian FD method is the most widely applied method for its high computational efficiency and stability in the polar and Marginal Ice Zone (MIZ) at large scale. It was also applied into other seas at meso-scale, such as Bohai Sea, Baltic Sea. Some new schemes, such as Line Successive Over-Relaxation (LSOR) and Alternative Direction Implicit (ADI), were adopted into the FD method to improve its computational precision. The most shortcoming of FD method is the obvious numerical diffusion in solving momentum and continuity equations, especially at the ice edge. To remedy this problem, the coupled Lagrangian and Eulerian PIC approach was established for sea ice dynamics at large and meso scales. In the PIC method, the sea ice in fixed cells is divided into a series ice particle. The ice mass in cells is adjusted with the drifting of Lagrangian particles, and the particle velocity is interpolated from Eulerian cells. In the Lagrangian SPH method, the Gaussian kernel function is used to integrate the ice parameters from discrete particles to continuous field, and the sea ice rehology can be described precisely with the drifting, deformation of ice particles. In the three methods above, Hibler's Viscous plastic constitutive law was used generally. In the GF method, the sea ice is simulated as discrete medium instead of the continuous medium assumed in other methods. The viscous-elastic-plastic law was established to model the interaction among ice particles, and the dynamics processes of ice ridging, rafting and breakup can be simulated at small scales. But the biggest cost of this increased accuracy is a significant increase in computational time when compared with other methods, especially in its application at large and meso scales. Thus, different numerical methods for the different demands for scale, precision or efficiency accordingly. Meanwhile, with the modification of existing methods, other new numerical methods, such as Arbitrary-Lagrangian-Eulerian (ALE), should be developed. Moreover, the study of numerical methods for sea ice dynamics should be coupled with other sea ice problems, such as constitutive law and thermodynamics, to improve the computational precision and efficiency comprehensively.

An unstructured CD-grid variational formulation for sea ice dynamics

Simulating sea-ice deformation in viscous-plastic sea-ice models with CD-grids

Multi-Scale Sea Ice Kinematics Modeling with a Grid Hierarchy in Community Earth System Model (version 1.2.1)

Comparison of sea ice kinematics at different resolutions modeled with a grid hierarchy in the Community Earth System Model (version 1.2.1)

Discretization of sea ice dynamics in the tangent plane to the sphere by a CD‐grid‐type finite element

Improving numerical accuracy for the viscous-plastic formulation of sea ice

ADVANCES IN NUMERICAL METHODS FOR SEA ICE DYNAMICS

CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5

Modeling sea ice in the marginal ice zone as a dense granular flow with rheology inferred from discrete element model data

CD-type discretization for sea ice dynamics in FESOM version 2

A Computational Fluid Dynamics Model for the Small-Scale Dynamics of Wave, Ice Floe and Interstitial Grease Ice Interaction

Mechanical and numerical models for sea ice dynamics on smallmeso scale

Modeling sea ice in the marginal ice zone as a dense granular flow with rheology inferred from a discrete element model

A dynamical core based on a discontinuous Galerkin method for higher-order finite-element sea ice modeling

Mechanical and Numerical Models for Sea Ice Dynamics on Small-Meso Scale

A Comparative Study of Viscoelastic Models for Ocean Wave Dissipation in Ice-Covered Regions of the Eastern Canadian Shelf

A nested high-resolution unstructured grid 3-D ocean-sea ice-ice shelf setup for numerical investigations of the Petermann ice shelf and fjord

Multi-scale Sea Ice Kinematics Modeling Dataset with a Tripolar Grid Hierarchy (TS Grids) in CESM

Dynamically Coupling Full Stokes and Shallow Shelf Approximation for Marine Ice Sheet Flow Using Elmer/Ice (V8.3)

SPH Modelling of Sea-ice Pack Dynamics

RESEARCH ADVANCES IN THE STUDY OF CONSTITUTIVE MODELS FOR SEA ICE DYNAMICS