ARTICLE DETAIL

资讯详情

深耕网站视觉设计与运营推广的一线实战洞察。

Ansys Fluent DPM实战:旋风分离器颗粒轨迹与分离效率仿真全攻略

Ansys Fluent DPM实战:旋风分离器颗粒轨迹与分离效率仿真全攻略 前阵子帮客户调一个旋风分离器的分离效率问题几何并不复杂但DPM参数设得不对仿真结果和实验差了十几个百分点。后来一圈圈排查最核心的坑居然出在颗粒追踪步数没给够壁面反弹系数设得太理想。今天就把这套完整流程写出来基于Ansys Fluent的DPM模型做旋风分离器颗粒轨迹仿真从几何清理、网格划分、物理模型选型到粒子注入、轨迹计算、分离效率统计最后把踩过的坑、排查思路都列出来案例文件的组织方式和使用顺序也会一并说明。这篇内容适合三类人刚接触DPM、想拿真实工业案例练手的Fluent新手已经做过简单管道颗粒模拟、但没碰过强旋流设备的朋友以及做气固分离设备设计、需要把分离效率仿真结果当成设计依据的工程师。文章不会只贴参数表每一个关键设置我都会解释“为什么这么设”以及设错了会出现什么诡异现象。1. 先搞清楚旋风分离器仿真到底在算什么东西1.1 设备原理与DPM模型的对应关系旋风分离器属于典型的惯性分离设备核心原理并不复杂含尘气流沿切向进入圆筒段后被迫做旋转运动颗粒因密度远大于气体受离心力作用被甩向器壁在壁面附近随下行气流滑落到锥体底部的排尘口净化后的气体则从中心排气管排出。这个过程中气体连续相一直在做强旋流运动而颗粒相当于被“扔”进流场里的离散个体。DPM模型全称Discrete Phase Model中文叫离散相模型它做的事情就是把这些颗粒当成一个个独立的质点在已求解好的连续相流场里做受力积分从而得到每条颗粒的运动轨迹。Fluent实现这个轨迹追踪时会逐个求解颗粒的运动方程把重力、拖曳力、Saffman升力、虚拟质量力等逐项叠加进去再通过时间步长积分得到位置。你可以把DPM理解成“在已经算好的风中扔羽毛”风是连续相羽毛是颗粒风不会因为羽毛的存在而改变方向这对应的是单向耦合。这个思路和旋风分离器特别契合因为工业条件下旋风分离器内颗粒体积浓度通常很低颗粒对气流的反作用可以忽略不计单向耦合已经足够。不过如果你的工况颗粒浓度很高比如入口浓度超过每立方米几十克甚至上百克单向耦合就不够用了得考虑双向耦合、DDPM甚至DEM。工程上怎么判断用不用双向耦合一个经验指标是颗粒质量加载率当颗粒质量流量与气体质量流量之比小于10%时单向耦合基本可靠超过这个值就要关注气固耦合效应。1.2 这个仿真能解决什么问题不能解决什么问题用DPM做旋风分离器仿真最核心的产出是两条一是颗粒轨迹可视化能直观看到不同粒径的颗粒在分离器内部怎么运动、最终被分离还是逃逸二是分离效率曲线也就是不同粒径颗粒的收集和逃逸比例。对于设备设计人员来说这两样东西直接决定选型、优化筒体尺寸、调整入口形式等方向。但我要泼一盆冷水DPM给出的分离效率是“单颗粒在时均流场中的统计结果”它没有考虑颗粒之间的碰撞、团聚、壁面再生夹带这些真实工业现象。这意味着仿真出的效率曲线在粗颗粒段比如50微米以上通常比较准但在细颗粒段比如1到5微米往往偏乐观。颗粒越小布朗运动、近壁湍流扩散的影响越大而普通DPM模型对这些效应的刻画非常粗糙。所以做这套仿真之前先想清楚自己要回答什么问题。如果是要看结构改动的相对变化趋势比如入口宽度从100毫米改成80毫米分离效率是变好还是变坏DPM完全够用如果是要出一个绝对效率值作为设备验收依据那就要谨慎。把DPM结果当趋势参考配合实验标定才符合工程实际。2. 前处理准备几何清理与网格划分的实操思路2.1 几何建模避坑与可参考的尺寸体系很多人在几何阶段就栽了跟头。旋风分离器看起来就是一个圆筒加一个锥体但真正决定仿真质量的细节都在“接口”和“插入深度”上。比如入口是矩形截面切向接入这个矩形入口的宽度、位置、与筒体相切的方式直接影响初始旋流强度又如排气管内筒的插入深度它决定含尘气流从排气管下端翻越进入内筒的路程长短这个长度对细颗粒分离影响非常大。我在实际项目中用过一个典型几何尺度体系筒体直径300毫米圆筒段高度400毫米锥体段高度400毫米矩形入口截面尺寸是100毫米乘50毫米排气管直径110毫米排气管插入深度150毫米排尘口直径90毫米。这套尺寸出处是常见的中试级旋风分离器文献数据非常适合作仿真复现。你完全可以按比例缩放但要注意缩放之后入口速度和筒体直径的匹配关系会变颗粒的斯托克斯数也会变分离特性不是线性变化的。用SpaceClaim或DesignModeler建模时有几个操作细节值得注意。第一流体域一定要建“气体实际流动区域”也就是筒体内部扣除排气管占用的体积如果直接把整个几何画出来再切分后面设置边界条件会很痛苦第二入口段最好留出一段平直过渡管长度约为入口水力直径的两到三倍这样能让入口气流充分发展避免入口端面直接顶着强湍流扰动第三所有倒角、圆角、螺栓孔这类结构细节在初版仿真里一律删掉它们是网格杀手对分离效率的贡献可以忽略。2.2 网格策略边界层、多面体与无关性验证旋风分离器是强旋流设备网格质量对结果的影响极其敏感这个敏感程度远超普通管道流动。原因是强旋流中湍流各向异性非常显著中心涡核区的速度梯度极大如果网格太粗涡核强度会在迭代过程中被人为耗散掉表现出来的结果就是旋流衰减过快、分离效率异常偏低。网格策略上我推荐用多面体网格或者六面体网格尽量避免纯四面体。多面体网格在复杂几何上生成效率高而且数值耗散比四面体低不少。边界层网格至少布置五到八层第一层高度对应的y值怎么控制关键看壁面函数的选择。用标准壁面函数时y最好落在30到300之间用增强壁面处理时y得小于1。旋风分离器内部流场结构复杂想让全壁面都满足y1代价很高所以我通常用标准壁面函数配合y约30到60算下来既稳又准。体网格规模方面中试尺寸的旋风分离器初始网格做到80万到150万量级通常能得到比较合理的流场。但这不是终点做完初版网格后一定要做网格无关性验证把网格加密到两三百万对比出口速度、压降和分离效率的变化幅度。一般规律是如果加密后分离效率变化超过几个百分点说明原网格还没收敛到网格无关状态需要继续加密。不要嫌麻烦这一步不做后面所有参数分析都是空中楼阁。3. 物理模型与DPM核心参数设置3.1 连续相用RSM为什么我不敢用k-epsilon湍流模型选择是这套仿真最重要的一步。旋风分离器内部是强旋流加强弯曲流线湍流表现出明显的各向异性。如果你用标准的k-epsilon模型它默认湍流粘性是各向同性的会严重高估涡核区域的湍流粘性导致旋流强度被过量耗散分离效率被急剧低估。实际仿真中会出现一个标志性异常颗粒本该被离心力甩到壁面结果几乎全部跟着气流从排气管跑了。解决这个问题有两种主流方案。第一种是雷诺应力模型RSM它直接求解六个雷诺应力分量组成的输运方程能够捕捉旋流中的各向异性计算代价比k-epsilon高但对旋风分离器这类强旋流设备来说RSM是工程精度和计算成本的平衡点第二种是大涡模拟LES精度更高能解析涡核的瞬态摆动代价是必须非稳态计算对网格和时间步长要求极高一台高性能工作站也要算好几天工程上一般不作为首选。我在工程算例中通常用RSM加标准壁面函数。压力插值格式建议用PRESTO!因为PRESTO!在处理强旋流和大幅压力梯度时比默认的Standard格式稳定得多动量、湍动能、雷诺应力都用二阶迎风格式一阶格式的数值耗散在旋流中过于严重基本不用考虑。操作压力设置为101325帕温度293开尔文气体密度用理想气体或者常密度都行入口速度按实际工况设定一般旋风分离器入口截面平均速度在15到25米每秒。3.2 DPM注入方式、粒径分布与颗粒力的取舍连续相流场收敛之后再看DPM设置。颗粒注入方式最贴近实际的是面注射选中入口截面作为颗粒发射面定义颗粒数量、速度分布、粒径分布和初始方向。颗粒数量我建议设置为1000到5000个之间太少统计效率时样本不足太多后处理时轨迹线密密麻麻看不清。速度分布可以直接用对应入口面连续相的速度分量这样颗粒一进流场就和当地气体同速更符合实际。粒径分布的设置要特别留意。如果你只看一个粒径比如20微米得到的是单粒径分离效率画不出分级效率曲线。更合理的做法是用Rosin-Rammler分布模拟实际粉尘的粒度分布需要给定中位径、分布指数和最小最大粒径。比如粉尘中位径20微米分布指数1.5粒径范围2到80微米Fluent会自动在范围内按分布规律生成多个粒径组。颗粒材料密度按实际粉尘设置石英砂大约2600千克每立方米煤粉大约1200到1500千克每立方米密度值直接影响离心力大小千万别随手填一个默认的1000。颗粒受力项怎么取舍在气固分离模拟中拖曳力是必须的它主导颗粒跟随气流的能力重力也要开着对粗颗粒的沉降趋势有影响。Saffman升力在近壁区剪切流中作用明显建议打开。虚拟质量力、压力梯度力、热泳力这些在常温常压气固分离场景下影响微乎其微开着只会增加计算负担和数值不稳定风险我一般直接关闭。3.3 边界条件语义Escape、Trap、Reflect别搞混这个坑我必须重点讲因为太多人栽在边界条件语义的歧义上。DPM模型的颗粒边界条件有几种类型escape表示颗粒穿过边界后离开计算域并停止追踪trap表示颗粒在该边界上被捕集也被定格在这里reflect表示颗粒与边界发生反弹反弹后的速度由恢复系数决定。这三种语义在旋风分离器里对应的是物理上的“逃逸出口”“被捕集壁面”和“反射壁面”。典型的设置方案是入口截面设置为escape表示颗粒如果从入口反着飞出去就算逃逸排气管出口设置为escape表示没被分离下来的颗粒从这里跟着净化气体跑出去了筒体和锥体壁面设置为reflect或trap。问题恰恰出在这里如果你用reflect模拟颗粒撞壁后的真实反弹那颗粒可能会反复弹跳模型中很难完全停住如果你用trap颗粒撞壁就被“没收”相当于假设壁面捕集效率100%这会高估分离效率。我在工程中常用的做法是壁面用reflect法向恢复系数设0.3切向恢复系数设0.7这组参数接近干燥颗粒与不锈钢壁面的典型碰撞行为排尘口底端设置成escape表示颗粒落到集尘箱里就算分离成功排气管出口设置成escape表示逃逸颗粒。这样统计分离效率时陷阱出在“排尘口如果设成reflect颗粒会在锥底部弹来弹去既不算逃逸也不算被捕集”导致最后粒子账目对不上。颗粒追踪步数也叫Max Number of Steps这个值别省旋风分离器里颗粒尤其是细颗粒会在内部转很多圈步数不够的话颗粒会没算完就被强行截断效率统计直接失真我一般设到50000以上。4. 求解流程两步法、收敛控制与中途暂停续算4.1 先连续相后颗粒的两步法DPM求解有两种典型路径一种是非耦合求解先把连续相流场算到收敛然后一次性注入颗粒、追踪轨迹另一种是耦合求解在每个迭代步中都更新连续相与颗粒之间的动量、质量和能量交换。对于旋风分离器这个场景我强烈建议先用非耦合两步法理由很直接颗粒浓度低单向耦合足够先拿到稳定收敛的旋流场再往里“放颗粒”逻辑清楚、调试方便。第一步关闭DPM只计算连续相入口设成速度入口出口设成压力出口用RSM做稳态迭代残差目标设到10的负4次方以下同时监控出口速度、进出口压差这两个物理量确认它们不再随迭代变化。第二步开启DPM保持连续相流场不变打开Discrete Phase Model面板中的Interaction with Continuous Phase选项不勾选表示只做单向追踪。然后做稳态颗粒追踪Fluent会基于已经“冻结”的流场一次性对所有注入颗粒进行轨迹积分。这里有个细节容易忽略在两步法里连续相的“冻结”状态是否合理取决于分离器内部是否存在明显的非稳态流动。旋风分离器内涡核天然存在进动现象也就是涡核会绕着轴线像摆锤一样晃动严格来说流场不是稳态的。RSM稳态算出来的是一个时间平均意义上的旋流场用它做单次轨迹追踪会丢失颗粒对涡核摆动的时间响应。工程上这个影响对粗颗粒不明显对细颗粒会有偏差。想更精确可以用冻结流场加非稳态粒子追踪也就是Unsteady Particle Tracking让颗粒在静止流场中按非稳态方式积分能部分补偿涡核摆动的影响。4.2 求解器与松弛因子怎么调求解器设置上压力速度耦合我推荐Coupled比SIMPLEC在强旋流中更稳健尽管内存占用稍高。空间离散格式前面提过压力用PRESTO!动量、湍动能、雷诺应力全用二阶迎风。亚松弛因子是另一个容易翻车的地方RSM加Coupled耦合时如果保持默认松弛因子初场质量不好很容易直接发散。我常用的调法是开始阶段把动量松弛因子降到0.3湍动能和湍流耗散率降到0.4雷诺应力松弛因子保持0.2到0.3之间先按低松弛迭代几百步等残差曲线起落平稳后再逐步调回默认值附近。残差曲线能不能作为收敛的唯一判据不能。强旋流工况下残差很难像普通管内流动那样掉到1e-6经常在1e-3到1e-4之间震荡。这时候你要看的是工程物理量出口面平均流速、进出口压降、排气管出口的质量流量不平衡率这几个量稳定了残差小幅波动可以接受。我最常用的监控方式是创建排气管出口面的面积加权平均速度、压力以及整个计算域的进出口质量流量差每十步记录一次当这些曲线都基本变成水平直线时连续相就算收敛了。收敛判定这块还有一个经验值可以参考进出口质量流量不平衡率要小于0.1%。如果这个值一直下不去大概率是局部回流或者压力出口边界离分离区太近可以把出口延伸段的长度加长让压力出口处在更稳定的区域。4.3 计算中途能关电脑吗保存、暂停与续算这个问题被问过太多次这里一次性讲透。先说结论Fluent计算中途完全可以暂停、保存甚至关掉电脑停机下次开机读入保存的case和data文件就能从断点继续迭代。核心前提是之前已经把计算进度成功保存成了文件。Fluent中保存完整状态涉及两个文件case文件存网格、物理模型、边界条件、求解器设置data文件存当前迭代步的流场、压力场、DPM轨迹结果等所有数据。任何想中断的时刻点File菜单下的Write Case和Write Data或者直接Save Case Data就会把当前状态落盘。那2024版本有没有更省事的办法有Auto Save功能。在求解计算页面里可以设置自动保存每迭代多少步保存一次case和data还可以开启只保存data、不保存case腾出磁盘空间。我习惯设置每500步保存一次data文件配合每两千步保存一次case文件这样即便计算结果爆炸或者软件崩溃最多丢失五百步的进度不至于从头再来。手工暂停也有两种方式一是直接在计算控制面板点Pause当前步的求解会停下来但不保存到磁盘二是先点Pause或Stop Calculation再手动保存文件。这里有个关键区别Pause只是驻留在内存里如果电脑断电或软件崩溃内存中的计算进展照样丢。稳妥的做法是先保存再退出顺序不要搞反。至于要不要升级到并行模式如果你用的是单机多核并行并行进程下保存文件时记得使用h5格式File偏好设置里选择Save Data as .dat.h5这种格式在大模型、多进程场景下保存更快、更不容易损坏。保存好后关机、重启、重新启动Fluent、读入刚才的case和data点Run Calculation继续迭代即可。5. 后处理与分离效率统计5.1 颗粒轨迹可视化的正确姿势颗粒轨迹的可视化不只是为了截图放到报告里更是诊断问题的工具。Fluent后处理中创建颗粒轨迹的方式是加载DPM结果后在图形区显示Particle Tracks通过输入粒子ID或者按释放面、粒径段来筛选不同类别的颗粒。我常用的做法是先看小颗粒再看大颗粒逐粒径段观察运动轨迹。小颗粒和大颗粒在旋风分离器中的典型轨迹差异非常明显。大颗粒的轨迹会在筒体上部很快向壁面靠拢然后螺旋下行到锥体底部基本不进入排气管下方的中心区域小颗粒则会明显跟随气流向中心涡核区运动多次在筒体内回旋出现“逃逸”趋势。如果在可视化中看到大颗粒也直奔中心排气管那就说明流场根本没转起来或者颗粒受力设置有误回到前面章节排查。轨迹显示时有个参数容易被忽略就是粒子显示密度。后处理软件中Particle Tracks面板里有个Particle Scale或Point Size设置如果粒径很小比如几微米默认显示尺寸在屏幕上就是一个点几乎看不见。这时候要统筹调整显示尺寸让轨迹视觉上清晰可见但这只是显示问题不影响实际数据。我可以给你一个排查技巧如果轨迹显示正常但看不到颗粒检查Media Definition里颗粒是否真的被注入了可以在DPM报告的Summary面板里查看注入颗粒总数和已追踪完成的粒子数是否一致。5.2 分离效率的计算口径分离效率的计算看起来简单其实口径很容易糊涂。物理上旋风分离器的总分离效率定义为被捕集颗粒质量占总注入颗粒质量的百分比。在DPM稳态追踪中Fluent会把每个终止于escape或trap边界的颗粒记录下来统计哪个边界捕集了多少颗粒。按我前面推荐的边界设置方案颗粒最终被统计为三类从排气管出口escape的算逃逸颗粒、从排尘口出口escape的算被捕集颗粒、从入口方向escape的算反常逃逸。总分离效率的算式就是N排尘口/(N排尘口N排气管出口)。这里N是颗粒数量。但工程上更关心的往往是分级效率也就是按粒径区间分别统计通常做法是在DPM设置中按Rosin-Rammler分布分档或者在注入面创建多个不同粒径的面注射组每组设一个固定粒径分别做一次DPM统计最后把每个粒径组的捕集率串成一条分级效率曲线。我用的最多的情况是设置六到八个粒径点从2微米、5微米、10微米、20微米、35微米、50微米到80微米每个粒径组单独追踪、单独统计这样既能快速诊断问题又能给设备选型提供直观依据。统计时要注意粒子数权重的口径问题。如果每组粒径释放相同数量颗粒统计出来的效率是“数量效率”不是“质量效率”。处理实际粉尘时单位体积内粗颗粒的质量远大于细颗粒做质量加权或者按入口粉尘的粒径分布权重计算得到的才是和工程对应的效率。所以表格中我会同时列出数量统计和质量统计建议你在正式报告中采用后者。5.3 后处理里的坑后处理阶段有几个坑我一次说全。第一个坑是颗粒轨迹的结果显示为一条直线仿佛颗粒完全不受流场影响。出现这种问题最可能的原因是DPM在稳态追踪时没有正确读取冻结流场或者连续相根本没有算收敛就开启了DPM。把连续相迭代充分收敛后重新追踪轨迹就会恢复了。第二个坑是DPM报告中的Escape和Trap总数和注入总数对不上少了一部分。这是Max Number of Steps给少了颗粒还在计算域内打转就被截断对应术语是Incomplete。遇到这种问题提升最大步数到十万甚至更大同时检查会不会有颗粒卡在壁面附近不断reflect一直无法离开这种情况下可以给超时步数设置一个外部边界条件让它最终被判定为逃逸或者捕集。第三个坑是分离效率算出来异常高比如超过99%。如果你用trap壁面且恢复系数设得很低那基本是设定出来的虚假结果而不是真实性能。我经常提醒团队的一点是如果结果好到超出直觉先别急着写报告花十分钟检查边界条件看是不是把入口设成了trap或者把透气性差的壁面设成了全捕集。6. 常见问题与排查技巧实录6.1 计算发散从残差暴涨到流场失稳连续相计算发散是RSM模拟强旋流时最常遇到的麻烦。现象一般是残差曲线在几百步后突然暴涨或者Coupled求解器报告压力残差异常增大紧接着出口流量变成无意义的值。排查顺序我从低到高列一下。第一步检查网格质量重点关注偏斜度超过0.85的单元用Mesh面板里的Quality指标查一下如果有大量高偏斜度单元旋流剧烈环境下极易发散第二步降低松弛因子把动量调到0.2到0.3湍流、雷诺应力也都降到0.3以下先稳住再慢慢加第三步把初始流场从简单均匀场换成分步初始化也就是先用一阶格式算几百步再切二阶格式第四步如果仍然发散看看入口速度是不是太高强旋流下Re数过高稳态RSM本身收敛就困难可以在入口段加长延伸段。我遇到最顽固的一次是RSM加Coupled怎么调都发散最后改用SIMPLEC加低松弛反而稳了。所以遇到疑难杂症别只钉在一种组合里换一换求解算法和松弛策略往往就破了局。6.2 颗粒行为异常穿透、反弹、卡在计算域颗粒行为异常是DPM仿真的第二大类问题。现象之一是大颗粒直接从排气管逃逸、完全看不到壁面聚集轨迹。这种基本可以确定是连续相流场旋流强度不够排查湍流模型是不是被k-epsilon或过粗网格耗散掉了也有可能是入口速度太低没有获得足够的离心力。将入口速度提上来或者把RSM开起来问题会立刻改观。现象之二是颗粒在壁面附近反弹个不停永远走不到排尘口。这种往往是恢复系数设得太保守法向恢复系数设成了0.1甚至更低颗粒撞墙后速度几乎全被吸收然后就贴着壁面慢慢蠕动。建议法向恢复系数放到0.3到0.5切向恢复系数放到0.7到0.9让颗粒保持向下滑动的趋势。现象之三是颗粒卡在某个角落既不逃逸也不被捕集粒子报告里大量显示Incomplete。这种基本是颗粒陷阱通常在锥体底部或者入口角落的局部回流区。处理方案有几个增大最大追踪步数适当增大颗粒尺寸显示但不改变物理粒径更有效地方式是修改几何消除过尖锐的局部微小区域让气流顺畅通过。6.3 仿真与实验偏差过大先检查这五个地方最后汇总一下做旋风分离器DPM仿真时和实验值对不上该怎么排查。按影响权重排序第一入口速度是否和实验一致强旋流对入口速度非常敏感入口速度差20%分离效率能差好几个点第二粒径分布是否和实验粉尘一致特别是中位径和分布宽度差一两个微米对细颗粒段效率影响巨大第三颗粒密度是不是填错了很多仿真把颗粒密度默认成水的密度而实际粉尘密度是它的两倍以上第四湍流模型和网格是否达到网格无关没做网格无关性验证时效率偏差可能完全掩盖真实规律第五分离效率统计口径是否一致实验往往用质量效率仿真如果用数量效率细颗粒占比高时差距会很大。这五个检查点我建议做成一张标准排查清单每次做完仿真和实验对标时逐条过一遍。多数情况下偏差问题都能在这五项里找到答案而不是去调整一些边边角角的模型参数。7. 案例文件内容与使用顺序说明7.1 文件清单标题里说过“附案例文件”这里就该说明白文件里到底有什么、怎么用。我按个人习惯把案例落盘成一个标准目录方便完整复现这套流程。geometry文件夹存放旋风分离器流体域的step几何文件命名如cyclone_fluid_domain.step可用SpaceClaim或DesignModeler直接打开。mesh文件夹存放划分好的网格文件cyclone.msh已经做过边界命名主流Fluent版本都能直接读入。case文件夹存放设置完成的Fluent case文件和初始data文件命名如dpm_cyclone.cas.gz和dpm_cyclone.dat.gz这是最关键的资产把边界条件、物理模型、DPM设置全部固化好了。post文件夹存放后处理可用的示例轨迹数据、切面云图文件和CFD-Post工程文件比如particle_tracks.cdb用来核对你的后处理效果。docs文件夹一份简单的README记录操作顺序、版本要求和一个已经跑通的设计工况表。按这套目录整理案例文件既方便自己回顾也方便发给同事或学生照着操作不用再从零开始摸索。7.2 推荐的操作路径拿到案例文件后我推荐按三步走。第一步用Fluent直接读入dpm_cyclone.cas.gz和dpm_cyclone.dat.gz打开后先别急着重算在当前界面中逐项检查连续相的湍流模型、入口边界条件、DPM模型的注入设置和边界条件对照本文第三章的内容核对一遍熟悉各项参数在面板中的位置和含义。第二步在已经收敛的流场基础上重新开启DPM稳态追踪按本文第四章和第五章的方法完成颗粒轨迹计算和分离效率统计并将结果与post文件夹里的示例数据对比。第三步是刻意练习把入口截面宽度、排气管直径或插入深度改掉重新计算连续相和DPM比较分离效率的变化趋势。这种改动练习比单纯跑通一个案例更能训练工程思维。案例文件里我放了一个已经跑通的基准工况入口速度20米每秒颗粒密度2600千克每立方米中位径20微米RSM加标准壁面函数壁面恢复系数法向0.3、切向0.7。你拿这套基准工况先跑如果数据能和post文件夹里的示例吻合说明环境设置一切正常再放手去做自己的变参数研究。我个人的体会是DPM模型在旋风分离器这个场景里真正决定结果质量的是这么三件事网格和湍流模型能不能把旋流强度保住、颗粒追踪步数给得够不够、各出口边界条件语义是否清楚。第一点决定流场骨架是否可靠第二点决定颗粒轨迹积分是否完整第三点决定效率统计是否真实。尤其是最后一点很多人误把排尘口设成reflect导致颗粒沉降到排尘口后又弹回去效率统计从一开始就是废的。最后再分享一个小技巧如果你要做一系列变参数对比比如不同入口速度下的效率变化建议在Geometry阶段就把入口段做成参数化尺寸配合Fluent的Parametric Study功能批量跑完再自动汇总效率会高很多。这个方向后续可以继续做双向耦合、DDPM甚至和DEM联用那些是更高阶的玩法但先把这套DPM全流程跑通、跑准是最值得投入的一步。
返回列表