首页/文章/ 详情

LPT简介及曳力模型探讨

1年前浏览3091
拉格朗日颗粒追踪(Lagrangian Particle Tracking)是计算流体力学CFD仿真中处理气固或液固流的常用方法,典型应用场景包括燃油雾化、悬浮颗粒、流化床、旋风分离器、化学反应器等。

图1 典型应用场景

当我们做纯流体仿真的结果后处理时,经常会使用流线streamline功能用于分析流体的流动特性,如图2所示。流线以拉格朗日法的角度可以理解为在追踪无质量的细小颗粒,颗粒运动完全依赖于流场运动。

图2 流体仿真的流线后处理

然而,现实中颗粒在重力、曳力以及其它相关力的作用下通常会表现出如图3所示的运动特性,即颗粒在重力作用下向下运动,在曳力作用下滞后于流体运动

图3 有质量颗粒的运动特性

以单一球形颗粒举例来看,颗粒质量对颗粒运动轨迹的影响如图4所示。  

图4 颗粒质量对颗粒运动的影响

颗粒控制方程为:

在不考虑颗粒质量情况下,颗粒速度等于流场速度,即: 

 

而在考虑颗粒质量时,颗粒速度不等于流场速度,即:

 
拉格朗日法利用牛顿第二定律作为控制方程求解颗粒速度Up   

其中曳力通常为最大的力,因此在接下来的算例中主要讨论曳力的影响


图5 颗粒受力示意图



不同曳力模型对仿真结果的影响

Fluent、STAR-CCM+以及VirtualFlow均提供不同的曳力模型供用户选择。以下模型根据横风颗粒算例中颗粒属性选取,但受限于不同软件植入了不同的曳力模型,无法进行同种曳力模型不同软件的对比。另外,在实际工程算例计算中,往往需要计算大量的颗粒,不同颗粒具有不同的属性,因此需要选取一个通用模型。

  • Fluent曳力模型

Fluent对于球形颗粒采用Morsi-Alexander曳力模型,其中各系数根据粒子雷诺数的不同有不同的取值,表达式为:   

  • STAR-CCM+曳力模型
STAR-CCM+采用Schiller-Naumann曳力模型,同样根据雷诺数的不同有不同的取值,表达式为:  
 
  • VirtualFlow曳力模型

VirtualFlow采用了Clift-Gauvin曳力模型,其表达式为:


接下来,采用不同CFD软件模拟单个油滴颗粒在横风作用下的运动情况,计算域与边界条件信息如图6所示。  

图6 横风颗粒算例计算域与边界信息

其它相关设定如下:

• 瞬态仿真

• 单相耦合

• 仅考虑曳力影响

• k-epsilon湍流模型

• 气相为空气属性

• 颗粒为油滴属性

• 颗粒直径0.5mm

Fluent、STAR-CCM+及VirtualFlow三款软件对单个油滴颗粒在横风作用下的运动轨迹的预测结果如图7所示。

图7 不同软件预测结果

可以看出,对单个油滴颗粒在横风作用下的运动轨迹的预测,不同软件预测结果整体趋势相同,但颗粒轨迹横向偏移不同。其中,VirtualFlow预测的颗粒轨迹横向偏移最少,STAR-CCM+结果较大,Fluent结果介于两者之间。计算偏差主要源自颗粒曳力模型选取的不同。曳力模型及其他相关受力模型的选取要综合考虑流场属性,颗粒与流场耦合作用等多方面决定,并通过试验数据进行校核,才能最大程度上保证仿真的准确性。

本篇作为拉格朗日颗粒追踪方法的介绍篇,希望先帮助大家建立一个初步认知。曳力模型对比部分由于缺乏试验对标数据,未作结论性总结,下次有机会再做定量对比分析,抛砖引玉,欢迎大家在评论区留言讨论。


来源:多相流在线
Fluent多相流燃烧化学湍流通用航空航天船舶轨道交通控制积鼎 CFD
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2023-06-23
最近编辑:1年前
积鼎科技
联系我们13162025768
获赞 108粉丝 110文章 300课程 0
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈