基于离散元方法的极地浮碎冰区船舶冰阻力

采用离散元方法建立浮/碎冰模型,结合欧拉多相流VOF方法数值模拟船舶与浮/碎冰之间的接触碰撞过程。重点对比分析船、冰作用过程中Hertz-Mindlin模型、Linear Spring模型和Walton Braun模型3种接触力计算模型对船-冰作用形式和冰阻力数值的影响。结果表明,在同一工况下,3种接触碰撞模型对船体周围碎冰运动模式变化的影响一致。对于船体所受冰阻力,通过分析船长方向和船宽方向的冰阻力数值,并将计算结果与采用Du Brovin经验公式计算得到的结果相对比,得出采用Linear Spring接触力计算模型计算得到的阻力幅值发散低,数值相对稳定,计算误差较小,为极地冰区船舶与浮碎冰相互作用的阻力计算提供参考。