用于模拟管道颗粒两相流动的非解析方法及电子设备技术

本公开提供了一种用于模拟管道颗粒两相流动的非解析方法及电子设备。其中所述方法包括：遍历颗粒周边的网格；计算颗粒附近流场的平均孔隙率与速度；计算颗粒所受的拖曳力大小；进行颗粒与流场的耦合，得到颗粒的运动情况和流场的流动情况。本公开对一种用于模拟管道颗粒两相流动的高精度非解析算法能在颗粒大小与网格尺寸相近的情况下很好地描述流场中的颗粒运动，降低了CFD‑DEM耦合算法对网格尺寸的要求，可用于处理液固、气固流化床中的两相流动问题。

【技术实现步骤摘要】
用于模拟管道颗粒两相流动的非解析方法及电子设备
本公开属于计算力学
，更具体的涉及一种用于模拟管道颗粒两相流动的非解析方法及电子设备。
技术介绍
流固两相流动广泛存在与自然环境与工业生产中，例如河流中的沙石沉降、化工领域的喷动床等。采用计算力学的手段获得流场和颗粒的运动情况具有很高的实际应用价值，计算结果可用于预测自然形貌的变化、指导工业设备的设计等。计算流体动力学-离散元(CFD-DEM)耦合算法是一类用于计算流固两相流动的常用方法，采用纳维尔-斯托克斯(Navier-Stokes)方程计算流场的流动，采用离散元法计算颗粒的运动，并通过耦合算法得到流体与固体之间的相互作用。CFD-DEM耦合算法又可细分为两类：解析的耦合算法与非解析的耦合算法。在解析的耦合算法下，流体网格的尺寸需要远小于颗粒的大小，并且在颗粒表面附近需要进行动态加密，计算量较大，无法用于计算成千上万个颗粒的运动情况。而在非解析的耦合算法下，网格尺寸一般是颗粒直径的三倍以上，用拖曳力模型表征颗粒与流场之间的相互作用，计算量远小于解析的耦合算法，所以常用非解析的CFD-DEM耦合算法计算多个颗粒在流场中的运动。然本文档来自技高网...

【技术保护点】
1.一种用于模拟管道颗粒两相流动的非解析方法，包括：设置流场和流场内颗粒的初始条件和边界条件，对所述流场区域划分网格；采用希尔伯特曲线遍历颗粒周边的网格；确定颗粒附近流场的平均孔隙率与速度；确定颗粒所受的拖曳力；根据之前所述的拖曳力，进行颗粒与流场的耦合，将拖曳力施加至颗粒及流体上，确定颗粒的运动情况和流场的流动情况。

【技术特征摘要】
1.一种用于模拟管道颗粒两相流动的非解析方法，包括：设置流场和流场内颗粒的初始条件和边界条件，对所述流场区域划分网格；采用希尔伯特曲线遍历颗粒周边的网格；确定颗粒附近流场的平均孔隙率与速度；确定颗粒所受的拖曳力；根据之前所述的拖曳力，进行颗粒与流场的耦合，将拖曳力施加至颗粒及流体上，确定颗粒的运动情况和流场的流动情况。2.根据权利要求1所述的方法，其中，用希尔伯特曲线对流场中的网格进行重新编号，并沿该曲线对颗粒周边的网格进行遍历。3.根据权利要求1所述的方法，其中，所述确定颗粒附近流场的平均孔隙率与速度，包括：采用加权平均的形式计算平均孔隙率和速度，其表达式如下：式中，是平均孔隙率，i是网格编号，αi是第i个网格的孔隙率，Wα是计算平均孔隙率时所用的权函数，ri是颗粒与第i个网格的中心之间的距离，ΔVi是第i个网格的体积，是颗粒附近流场的平均速度，Uf，i是第i个网格内的流场速度，Wu是计算流场平均速度时所用的权函数。4.根据权利要求3所述的方法，其中，计算平均孔隙率所用的权函数所述Wα为阶跃函数，其表达式如下：式中，r是网格中心与颗粒中心之...

【专利技术属性】
技术研发人员：刘谋斌，滕郁骏，王泽坤，
申请(专利权)人：北京大学，
类型：发明
国别省市：北京,11