A HYBRID LAGRANGIAN-EULERIAN NUMERICAL MODEL FOR SEA ICE DYNAMICS

摘要

This paper presents a Hybrid Lagrangian-Eulerian (HLE) method for sea ice dynamics. The method combines the computational efficiency of Finite Difference Method (FDM) with the higher numerical accuracy of Lagragian Smoothed Particle Hydrodynamics (SPH). In this method, the sea ice cover is represented by a group of Lagrangian ice parcels with their own thicknesses and concentrations. These ice variables are interpolated to the Eulerian gird nodes using the SPH kernel interpolation. The FDM is used to calculate the ice velocities at Eulerian grid nodes. These velocities are interpolated, with the Gaussian interpolation function to give the velocity of each Lagrangian ice parcel for calculating their displacements. The model is validated with the analytical solution of ice ridging in a rectangular basin. It is also used to simulate ice dynamics in a vortex wind field, and compared with the results of PIC and SPH methods. The model is applied to simulate the sea ice dynamics in Bohai Sea, and compared well with the field data.