颗粒轨迹全靠猜?拉格朗日多相流没设对,粒子数量和壁面磨损全是错的

做含颗粒、液滴、气泡这类仿真的工程师,基本都被多相流折磨过:管道里明明有固体颗粒冲蚀弯头,喷嘴里明明有液滴在飞,结果一跑,粒子要么凭空消失,要么数量对不上,壁面磨损更是差得离谱。问题大多出在拉格朗日模型没设对。

多相流有两条路:欧拉法把各相当成连续介质,拉格朗日法把颗粒当成一个个离散的"包裹"去追踪。涉及颗粒轨迹、停留时间、冲蚀、喷雾、分离的场景,拉格朗日往往才是对的工具。今天就用 5 个点把 STAR-CCM+ 里拉格朗日多相流的基础和设置讲清楚。


先搞懂概念:粒子、包裹,以及拉格朗日到底在追踪什么


拉格朗日方法的核心,是把离散相看成一个个独立运动的粒子,逐个求解它们的运动轨迹。这里有两个关键概念:

Particle(粒子):真实的单个颗粒;

Parcel(包裹):为了节省计算量,把一群尺寸、速度、状态相同的粒子打包成一个"计算单元",一次追踪代表一整包。

粒子有尺寸和形状,而阻力系数、升力系数都和形状、尺寸直接相关。理解"包裹"是打包代表,才能明白为什么粒子数量能既对上物理、又不至于算到天荒地老。


控制方程:粒子的运动,靠的是这几个守恒


每个粒子的运动,本质上是在求解守恒方程:

动量守恒:粒子受到的力决定它的加速度和轨迹,力包括曳力(阻力)、重力、升力、压力梯度力等;

能量守恒:粒子和连续相之间有热交换时,决定粒子温度;

质量守恒:涉及蒸发、凝结、反应时,粒子质量会变化。

此外还有两个常被忽略的效应:湍流弥散——湍流的随机脉动会把粒子"抖"散开,不开启时粒子轨迹会过于整齐;重力与浮力——粒子和连续相之间的密度差会产生沉降或上浮。这些子模型按实际工况勾选,缺了就会和物理对不上。


参数设置:粒子模型、入射器,一个都不能少


软件里要先在拉格朗日模型中选对粒子类型,再定义粒子从哪、以什么方式进入:

有质量粒子/ 无质量粒子:无质量粒子只跟随流场(用于可视化流线),有质量粒子才受惯性、重力影响,冲蚀、分离必须用有质量;

Point injector(点入射器):从一个点喷出;

Coneinjector(锥入射器):以锥形喷雾方式喷出,适合喷嘴;

Surface injector(面入射器):从整个面均匀释放粒子。

入射器里要给粒子的尺寸分布、初始速度、质量流量或粒子数量。这一步给错,后面粒子数量、浓度全是错的。

图片

颗粒轨迹全靠猜?拉格朗日多相流没设对,粒子数量和壁面磨损全是错的的图2


易错细节:壁面边界条件,决定粒子是弹开还是消失


粒子撞到壁面后怎么办,是冲蚀和分离仿真的关键,STAR-CCM+ 提供了几种壁面边界:

Escape(逃逸):粒子到达边界就离开计算域,常用于出口;

Rebound(反弹):粒子撞壁后弹回,冲蚀、弯头必用;

Stick(粘附):粒子粘在壁面上;

Splash(飞溅):液滴撞壁后溅出更小的液滴;

Composite(组合):多种方式组合。

最常见的坑,就是弯头、壁面默认成了 Escape,粒子一撞就"消失",自然算不出冲蚀和反弹轨迹。涉及颗粒冲蚀的壁面,必须设成 Rebound 并配合冲蚀模型。

图片

颗粒轨迹全靠猜?拉格朗日多相流没设对,粒子数量和壁面磨损全是错的的图4


结果解读:看轨迹、看收敛,更要看粒子数对不对


跑完后,先看粒子轨迹是否符合物理预期——该沉降的沉降、该反弹的反弹、该分离的分离;再核对粒子数量、质量流量是否和入口设定守恒,有没有粒子莫名其妙地丢失。

同时要确认流场本身已经收敛,残差稳定下降到合理量级,否则粒子是在一个错误的流场里运动,轨迹再漂亮也没有意义。最后再结合冲蚀速率、停留时间等结果,去回答弯头哪里最容易磨、分离效率是多少这类真正的工程问题。

图片

颗粒轨迹全靠猜?拉格朗日多相流没设对,粒子数量和壁面磨损全是错的的图6


小结


拉格朗日多相流把离散相当成"包裹"逐个追踪,适合轨迹、冲蚀、喷雾、分离。要理解粒子与包裹、动量/能量/质量守恒和湍流弥散、重力浮力;按有无质量选粒子模型,用点/锥/面入射器定义入口;壁面按逃逸、反弹、粘附、飞溅正确设置;最后核对轨迹、粒子数守恒和流场收敛。把这些设对,颗粒仿真才敢和实测对标。


【配套资料领取】


扫码添加小助理,私信回复关键词【拉格朗日】,领取本篇配套资料包:①拉格朗日粒子模型选型表 ②粒子边界条件设置速查 ③实操清单

颗粒轨迹全靠猜?拉格朗日多相流没设对,粒子数量和壁面磨损全是错的的图7
查看全文
我的头像
默认 最新
ansys结构交流群