基于IS–FEM的颗粒曳力系数计算与试验验证  

PARTICLE DRAG COEFFICIENT CALCULATION BASED ON IS-FEM AND EXPERIMENTAL VERIFICATION

作　　者：苏杰 周长江[1] 姜琛 SU Jie;ZHOU Changjiang;JIANG Chen(State Key Laboratory of Advanced Design and Manufacture for Vehicle Body,Hunan University,Changsha 410082,China;Key Laboratory of Traffic Safety on Track of Ministry of Education,Central South University,Changsha 410075,China)

机构地区：[1]湖南大学汽车车身先进设计与制造国家重点实验室,长沙410082 [2]中南大学轨道交通安全教育部重点实验室,长沙410075

出　　处：《力学与实践》2023年第1期119-129,共11页Mechanics in Engineering

基　　金：国家自然科学基金(52075153,12002395);湖南省重点研发基金(2020WK2032)资助项目。

摘　　要：基于浸没光滑有限元模型(immersed smooth finite element model,IS-FEM),计算球形与非球形颗粒曳力系数。设计颗粒曳力图像测量试验,验证IS-FEM模拟精度。颗粒相的运动行为基于连续介质理论的光滑有限元法求解;流体控制方程通过特征分解的半隐式有限元法求解;颗粒与流体相通过非贴体网格交换数据。结果表明,球形颗粒流场特征对称分布,非球形颗粒稳定沉降时长轴与重力方向垂直。不同雷诺数下球形颗粒的曳力系数计算值与Stokes曳力系数一致,非球形颗粒曳力系数高于等效球形颗粒,IS-FEM计算值与沉降试验吻合良好。The drag coefficients of spherical and non-spherical particles are calculated based on the immersed smooth finite element model(IS-FEM). A particle drag image measurement experiment is designed to verify the numerical simulation result. The motion of particle is solved by smooth finite element method with continuum theory. The fluid governing equation is solved by semi-implicit characteristic-based split method. Data exchange between particle and fluid phase is carried out in non-fit mesh. The results show that the velocity field around spherical particles are symmetrically distributed, and the long axis of the non-spherical particles is perpendicular to the gravity direction when the particles are in a stable settlement. The calculated drag coefficients of spherical particles are consistent with the Stokes theory under different Reynolds numbers, and the drag coefficient of non-spherical particles is higher than that of equivalent spherical particles. The predicted value of IS-FEM is in good agreement with the experimental measurement.

关 键 词：曳力系数 浸没光滑有限元法 颗粒沉降 图像处理 

分 类 号：O35[理学—流体力学]