Modeling calving process of glacier with dilated polyhedral discrete element method  

作　　者：Lu Liu Ji Li Qizhen Sun Chunhua Li Sue Cook Shunying Ji 

机构地区：[1]State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian University of Technology,Dalian 116024,China [2]National Marine Environmental Forecasting Center,Beijing 100081,China [3]Institute for Marine and Antarctic Studies,University of Tasmania,Hobart 7001,Australia

出　　处：《Acta Oceanologica Sinica》2021年第7期159-169,共11页海洋学报（英文版）

基　　金：The National Key R&D Program of China under contract Nos 2018YFA0605902,2016YFC1402705,2016YFC1402706 and 2016YFC1401505;the National Natural Science Foundation of China under contract Nos 41576179 and 51639004;the fund of Australian Research Council’s Special Research Initiative for Antarctic Gateway Partnership under contract No.SR140300001;the China Postdoctoral Science Foundation under contract No.2020M670746.

摘　　要：Mass loss caused by glacier calving is one of the direct contributors to global sea level rise.Reliable calving laws are required for accurate modelling of ice sheet mass balance.Both continuous and discontinuous methods have been used for glacial calving simulations.In this study,the discrete element method(DEM)based on dilated polyhedral elements is introduced to simulate the calving process of a tidewater glacier.Dilated polyhedrons can be obtained from the Minkowski sum of a sphere and a core polyhedron.These elements can be utilized to generate a continuum ice material,where the interaction force between adjacent elements is modeled by constructing bonds at the joints of the common faces.A hybrid fracture model considering fracture energy is introduced.The viscous creep behavior of glaciers on long-term scales is not considered.By applying buoyancy and gravity to the modelled glacier,DEM results show that the calving process is caused by cracks which are initialized at the top of the glacier and spread to the bottom.The results demonstrate the feasibility of using the dilated polyhedral DEM method in glacier simulations,additionally allowing the fragment size of the breaking fragments to be counted.The relationship between crack propagation and internal stress in the glacier is analyzed during calving process.Through the analysis of the Mises stress and the normal stress between the elements,it is found that geometric changes caused by the glacier calving lead to the redistribution of the stress.The tensile stress between the elements is the main influencing factor of glacier ice failure.In addition,the element shape,glacier base friction and buoyancy are studied,the results show that the glacier model based on the dilated polyhedral DEM is sensitive to the above conditions.

关 键 词：glacier calving discrete element method dilated polyhedral element bond and fracture model 

分 类 号：P343.6[天文地球—水文科学]