颗粒轨迹全靠猜?拉格朗日多相流没设对,粒子数量和壁面磨损全是错的
做含颗粒、液滴、气泡这类仿真的工程师,基本都被多相流折磨过:管道里明明有固体颗粒冲蚀弯头,喷嘴里明明有液滴在飞,结果一跑,粒子要么凭空消失,要么数量对不上,壁面磨损更是差得离谱。问题大多出在拉格朗日模型没设对。
多相流有两条路:欧拉法把各相当成连续介质,拉格朗日法把颗粒当成一个个离散的"包裹"去追踪。涉及颗粒轨迹、停留时间、冲蚀、喷雾、分离的场景,拉格朗日往往才是对的工具。今天就用 5 个点把 STAR-CCM+ 里拉格朗日多相流的基础和设置讲清楚。
先搞懂概念:粒子、包裹,以及拉格朗日到底在追踪什么
拉格朗日方法的核心,是把离散相看成一个个独立运动的粒子,逐个求解它们的运动轨迹。这里有两个关键概念:
Particle(粒子):真实的单个颗粒;
Parcel(包裹):为了节省计算量,把一群尺寸、速度、状态相同的粒子打包成一个"计算单元",一次追踪代表一整包。
粒子有尺寸和形状,而阻力系数、升力系数都和形状、尺寸直接相关。理解"包裹"是打包代表,才能明白为什么粒子数量能既对上物理、又不至于算到天荒地老。
控制方程:粒子的运动,靠的是这几个守恒
每个粒子的运动,本质上是在求解守恒方程:
动量守恒:粒子受到的力决定它的加速度和轨迹,力包括曳力(阻力)、重力、升力、压力梯度力等;
能量守恒:粒子和连续相之间有热交换时,决定粒子温度;
质量守恒:涉及蒸发、凝结、反应时,粒子质量会变化。
此外还有两个常被忽略的效应:湍流弥散——湍流的随机脉动会把粒子"抖"散开,不开启时粒子轨迹会过于整齐;重力与浮力——粒子和连续相之间的密度差会产生沉降或上浮。这些子模型按实际工况勾选,缺了就会和物理对不上。
参数设置:粒子模型、入射器,一个都不能少
软件里要先在拉格朗日模型中选对粒子类型,再定义粒子从哪、以什么方式进入:
有质量粒子/ 无质量粒子:无质量粒子只跟随流场(用于可视化流线),有质量粒子才受惯性、重力影响,冲蚀、分离必须用有质量;
Point injector(点入射器):从一个点喷出;
Coneinjector(锥入射器):以锥形喷雾方式喷出,适合喷嘴;
Surface injector(面入射器):从整个面均匀释放粒子。
入射器里要给粒子的尺寸分布、初始速度、质量流量或粒子数量。这一步给错,后面粒子数量、浓度全是错的。
易错细节:壁面边界条件,决定粒子是弹开还是消失
粒子撞到壁面后怎么办,是冲蚀和分离仿真的关键,STAR-CCM+ 提供了几种壁面边界:
Escape(逃逸):粒子到达边界就离开计算域,常用于出口;
Rebound(反弹):粒子撞壁后弹回,冲蚀、弯头必用;
Stick(粘附):粒子粘在壁面上;
Splash(飞溅):液滴撞壁后溅出更小的液滴;
Composite(组合):多种方式组合。
最常见的坑,就是弯头、壁面默认成了 Escape,粒子一撞就"消失",自然算不出冲蚀和反弹轨迹。涉及颗粒冲蚀的壁面,必须设成 Rebound 并配合冲蚀模型。
结果解读:看轨迹、看收敛,更要看粒子数对不对
跑完后,先看粒子轨迹是否符合物理预期——该沉降的沉降、该反弹的反弹、该分离的分离;再核对粒子数量、质量流量是否和入口设定守恒,有没有粒子莫名其妙地丢失。
同时要确认流场本身已经收敛,残差稳定下降到合理量级,否则粒子是在一个错误的流场里运动,轨迹再漂亮也没有意义。最后再结合冲蚀速率、停留时间等结果,去回答弯头哪里最容易磨、分离效率是多少这类真正的工程问题。
小结
拉格朗日多相流把离散相当成"包裹"逐个追踪,适合轨迹、冲蚀、喷雾、分离。要理解粒子与包裹、动量/能量/质量守恒和湍流弥散、重力浮力;按有无质量选粒子模型,用点/锥/面入射器定义入口;壁面按逃逸、反弹、粘附、飞溅正确设置;最后核对轨迹、粒子数守恒和流场收敛。把这些设对,颗粒仿真才敢和实测对标。
【配套资料领取】
扫码添加小助理,私信回复关键词【拉格朗日】,领取本篇配套资料包:①拉格朗日粒子模型选型表 ②粒子边界条件设置速查 ③实操清单










评论0