A polyhedral discrete element method for sea ice dynamic process in polar regions
Received:January 18, 2024  Revised:March 12, 2024
View Full Text  View/Add Comment  Download reader
DOI:10.7511/jslx20240118002
KeyWord:sea ice dynamics  discrete element method  ridge model  sea ice fracture
           
AuthorInstitution
李继 大连理工大学 工业装备结构分析优化与CAE软件全国重点实验室, 大连
王嗣强 大连理工大学 工业装备结构分析优化与CAE软件全国重点实验室, 大连
王安良 国家海洋环境预报中心, 北京
季顺迎 大连理工大学 工业装备结构分析优化与CAE软件全国重点实验室, 大连
Hits: 436
Download times: 128
Abstract:
      Sea ice in polar regions has a significant impact on navigational conditions and even climate predictions at geological scales.Addressing the shortcomings of traditional continuum-based grid models for sea ice simulation,the discrete element method (DEM) has emerged as a popular tool in floe dynamics research.In this study,we establish a discrete element method suitable for simulating polar-scale sea ice dynamics.This method employs polyhedrons based on surface triangulation to describe floating ice.The Gilbert-Johnson-Keerthi(GJK) algorithm is employed to enhance computational efficiency.The contact forces between elements are determined based on the energy conservation theory.A bond-fracture model describes the formation of level ice from freezing floes,incorporating plastic forces during the ridge formation process as proposed by Hopkins.Firstly,the accuracy of the contact force calculations is verified through the simulations of collisions between two idealized floes.Subsequently,the capability to simulate ridge formation is validated through the floe field compression under a unidirectional uniform wind.Finally,the influence of bond strength parameters on sea ice fracturing and ridge formation under shear wind conditions is studied.The results indicate that this model effectively captures the dynamics of sea ice on a geological scale and accurately depicts linear kinematic features,such as crevices and ridges,caused by linear motion in ice covers.