基于DEM-LBM耦合算法的冰水混合体数值模型研究.pdfVIP

基于DEM-LBM耦合算法的冰水混合体数值模型研究

摘要

冰水混合物为极地船舶最常见的物质形态，主要分布在碎冰航道内以及冰脊的龙骨

部分。为了探究冰水混合物的力学特性，本文开发了DEM-IMB-LBM（离散元—浸没移

动边界—格子玻尔兹曼方法）耦合算法，建立了冰水混合体数值模型进行研究。

数值模拟主要以多边形为研究对象，建立了其与流体相互作用的流固耦合算法，针

对冰脊龙骨内部碎冰间隙的流场影响因素进行了研究，其中格子玻尔兹曼方法来模拟流

体，固体颗粒之间以及与壁面的相互作用采用离散元方法求解，多边形颗粒的复杂动边

界使用浸没移动边界法进行处理。

本文首先采用提出的DEM-IMB-LBM耦合算法对近似球形的正32多边形进行单颗

粒沉降数值模拟，改变格子尺寸、距容器壁的远近等影响因素对沉降模拟进行分析，并

对比前人工作成果；同时，还验证双正32边形的DKT现象（牵引—接触—翻滚），进

一步验证耦合算法DEM模块的准确性。

大量学者探究冰-水相关的流固耦合问题，大都以球形单元为主流，以方形单元的研

究文献明显少的多，因此在完成DEM-IMB-LBM耦合算法的验证之后，对方形颗粒进

行沉降数值模拟，同时研究了雷诺数，颗粒长宽比，容器宽度等因素对方形颗粒沉降的

影响，也验证了双方形颗粒的“倒T形”结构现象以及不同位置的双方形颗粒沉降特性，

并结合前面的工作进一步验证程序。

自然界海冰形状多种多样，其中不乏有大量凹多边形形态的碎冰存在。因此，对前

面的多边形沉降数值模拟进行拓展，模拟了以两种凹多边形形态为代表的凹多边形颗粒

群的沉降数值模拟以及混合多边形的沉降数值模拟，模拟效果良好。以冰脊水线以下龙

骨部分作为研究对象，构建冰水混合体数值模型，分别探究了加载板速度、颗粒密集度

以及颗粒形态的影响规律。

关键词：格子玻尔兹曼方法；离散元方法；浸没移动边界法；流固耦合；冰水混合

物

基于DEM-LBM耦合算法的冰水混合体数值模型研究

Abstract

Ice-watermixturesarethemostcommonmaterialformsofpolarships,mainlydistributed

inthechannelofcrushediceandthekeelpartoftheiceridge.Inordertoexplorethemechanical

propertiesofice-watermixtures,aDEM-IMB-LBM(discreteelement-immersionmoving

boundary-latticeBoltzmannmethod)couplingalgorithmwasdevelopedandanumericalmodel

ofice-watermixturewasestablishedforstudy.

Thenumericalsimulationmainlytakespolygonsastheresearchobject,andestablishesa

fluid-structureinteractionalgorithmforitsinteractionwithfluids,andtheflowfieldinfluencing

factorsoftheicecrushedicegapinsidetheiceridgekeelarestudied,inwhichthelattice

Boltzmannmethodsimulatesthefluid,theinteractionbetweensolidparticlesandwiththewall

issolvedbythediscreteelementmethod,andthecomplexmovingboundaryofpolygon

particlesisprocessedby