摘要

光滑粒子流体动力学(SPH)方法采用离散分布的粒子作为离散工具,借助粒子和核函数将场变量及其导数以某种加权求和的方式进行计算,使得该方法更适于处理包含大位移、动边界和自由表面的问题。本文在分析动床条件下的溃坝问题时,以Drunker-Prager准则作为底床颗粒物的侵蚀启动判断条件。满足该条件而发生启动的颗粒物将以非牛顿流体的特征发生运动,即以Herschel-Bulkley-Papanastasiou模型来描述;否则,认为底床颗粒物未启动而保持静止状态。为进一步改进计算结果,采用基于体积比的两相SPH方法描述固液混合体的动力学特征,并基于固液混合颗粒流的动力学特征对混合黏度的计算进行改进。基于上述原理,分析动床条件下的溃坝问题,并将数值计算结果与已有的实验数据进行对比,证明本文提出的混合黏度计算公式改进了数值计算结果。