ARTICLE · INTELLIGENCE

战地情报 · 详情页

来自尧图项目组的一线实战观察与深度解析

CFD仿真进阶:多相流模型选型、动网格与y+湍流模型搭配实战

CFD仿真进阶:多相流模型选型、动网格与y+湍流模型搭配实战 简介ANSYS FLUENT 中文培训教材 PPT part2 以“边界条件”为主线面向使用 ANSYS FLUENT 进行 CFD 仿真分析的初学者与工程技术人员帮助读者厘清流体域、固体域、多孔介质等概念掌握确定边界位置、选择边界类型及输入边界数据的基本思路。压缩包内共包含 1 个 PPT 文件容量约 1.98MB图文并茂适合随课程渐进学习也可作为日常设置入口、出口、壁面、对称面及周期边界条件时的速查参考。预览内容覆盖速度进口、压力进口、流量入口、压力出口等常用边界条件的特点与适用场景并给出减少近边界网格扭曲、避免回流影响等实用建议能够有效降低初学者因边界条件设置不当带来的计算发散或结果失真风险。内容还涉及压力入口与流量入口在可压缩/不可压缩流动中的不同处理方式以及壁面热边界条件、多孔介质阻力系数等设定要点。资源已有 92 人学习下载说明其作为 FLUENT 入门配套材料具备一定参考价值。 一口气连做三个CFD案例、软件界面已经玩得挺溜、残差曲线掉到1e-4就心里踏实——这是很多人在学ANSYS FLUENT时读完培训教材Part 1的状态。然后他们翻开《ANSYS_FLUENT中文培训教材PPT part2》翻到多相流那一章突然发现教材不教点哪里了开始跟你讲VOF的界面重构、讲Eulerian模型里相间曳力系数、讲y到底应该取30还是取1。很多人就是在这里放弃的因为他们把第二册当成又一本操作手册来读而Part 2真正的意图是逼你从会点鼠标转向会做判断。这篇内容不打算替你复述PPT里的每一页而是把Part 2阶段最值得死磕的五个硬骨头拆开来讲——分别是模型选型逻辑、多相流的应用边界、运动区域与动网格的处理思路、收敛性判断、以及y和湍流模型的搭配。顺便把我这十多年里在上面踩过的坑、得出的判断经验一起交代清楚给正在啃这套教材第二册的人一条相对顺畅的路径。1. 教程Part 1到Part 2的本质差异从“按钮位置”到“决策依据”很多学员觉得Part 2内容变难了难在教材里突然冒出了大量模型名称和适用范围。这其实是培训体系设计的必然第一册解决的是能不能跑通一个算例第二册解决的是面对一个没有标准答案的工程问题你凭什么选择某条求解路径。1.1 第一册教操作第二册教判断Part 1的内容一般围绕几何导入、网格划分、边界条件设置、求解器参数、一个经典层流或简单湍流案例展开。你跟着做一遍得到一条漂亮的残差曲线会觉得自己已经会CFD了。但到了真正的工程项目里问题不会自己告诉你该用VOF还是Mixture也不会告诉你旋流器这种强旋转流动到底该用哪套湍流模型——这正是Part 2要解决的。我自己的经验是读第二册时不要按章节顺序硬啃而是要带着问题去读。每学一个新模型先问三个问题这个模型的物理假设是什么它把真实流动简化掉了哪些东西它的计算开销处在什么量级有没有成本更低的替代方案如果选错了模型结果会偏在哪个方向是偏保守还是偏危险带着这三个问题去翻教材你会注意到很多以前忽略的细节比如Mixture模型那节教材里一定会写适用于相间滑移速度较小的场合这句话在考试里可能只是划重点在工程里就是选型红线。1.2 学完第二册你该具备哪些能力我习惯把Part 2的学习目标归纳成一张能力清单读者可以拿来自检能力维度学完后的具体表现物理问题归类看到工程问题能快速判断属于单相/多相/传热/运动边界中的哪类模型选型能说清楚为什么选A模型而不选B模型前提假设是什么数值参数调整残差不收敛时知道先动哪些参数而不是随机尝试结果可信度判断能通过质量守恒、y范围、网格无关性验证等途径自证结果可靠这一章内容本身不长但它决定了后面几章你能吸收多少。如果现在让你写一个选型决策表还是觉得无从下手建议先把Part 1的案例重做一遍这次不要照着PPT点试着关掉教程自己完成整个设置流程。点击速度会慢一半但对每个设置项的理解会深一倍。2. 多相流模型选型VOF、Mixture、Eulerian的边界条件多相流是Part 2教材里篇幅最大的板块之一也是最容易让人一头雾水的板块。FLUENT里多相流模型有好几个每个底下还有一堆子模型、相间相互作用参数新手很容易陷进去出不来。我的建议是先把三个主力模型的边界死死划清楚。2.1 三种模型的一句话区分用一句话概括它们各自的定位VOF流体体积法所有相共享同一套动量和压力场用体积分数追踪相界面适合界面清晰、两相不相互贯穿的场合。Mixture混合物模型把多相当作一种混合流体求解混合动量方程再通过滑移速度描述相间相对运动适合相间耦合较强、分布相对均匀的场合。Eulerian欧拉模型每一相单独建立质量、动量方程通过相间交换系数耦合适合相间作用复杂、各相运动差异明显的场合。这里的关键判断点在于界面是否清晰和相间是否有明显滑移。做液舱晃荡、溃坝、射流破碎这类界面分明的自由液面问题优先VOF做气泡在液体里弥散上升、沉降槽颗粒浓度分布这类问题Mixture往往比VOF更实用做流化床、气力输送、高浓度颗粒流这类相间强烈耦合的问题Eulerian是绕不开的选择。2.2 选型决策表与典型应用说得再直白一点可以按下面这张表快速定位工况特征推荐模型选型理由自由液面、液-气分界面明显VOF界面捕捉算法成熟能还原清晰的相界面气泡/液滴弥散在连续相中界面尺度小Mixture 或 EulerianVOF对每个气泡都要解析界面计算量无法承受颗粒浓度高、相间曳力主导Eulerian需要完整的相间动量交换描述只需要总体浓度分布不关心局部界面细节Mixture方程数量少收敛快工程精度够用说一个我自己处理过的搅拌槽气液分散案例。刚开始选型时想用VOF看得清楚结果网格加密到接近八百万还是看不到完整气泡群气含率与实验偏差很大。后来换回Eulerian模型配合合适的曳力系数一小时就算完了气含率误差控制在百分之十以内。这个案例让我记住一条原则VOF不是万能的计算资源有限时模型简化程度要跟问题需要的精度相匹配。2.3 VOF使用中最容易忽略的Courant数问题如果最终确定用VOF有一个参数新手经常忽略Free Surface Courant Number。这个值一般建议控制在1以内最好不超过2否则自由液面会出现破碎、扩散体积分数云图看起来像一团浆糊。这也解释了为什么VOF计算通常需要较小的时间步。瞬态计算时时间步长设太大会导致界面处每个时间步内流体移动超过一个网格尺度界面重构算法Geo-Reconstruct就会失效。具体设置上可以在求解器里勾选Variable Time Step让软件根据Courant数自动调整时间步长效率和安全性能兼顾。另外当液滴或气泡尺寸与网格尺度接近时是否考虑表面张力模型也是个决策点。判断标准可以看Weber数惯性力主导就忽略表面张力表面张力主导就必须开启表面张力模型并指定接触角。教材这部分通常讲得比较简略工程上碰到喷雾、液膜这类问题表面张力的影响往往比想象中大得多。3. 动区域与动网格别把MRF、滑移网格和动态层混为一谈第二册里运动边界相关的章节是另一个容易让人犯迷糊的地方。旋转机械、搅拌器、活塞运动、阀门开闭这些工况在工程里太常见了但处理方式却分好几个层次。3.1 三种运动处理方式的使用边界同样用一句话理清MRF多重参考系把旋转区域放到旋转坐标系里网格本身不转是一种稳态近似用于 rotor-stator 相互作用不强的工况。滑移网格Sliding Mesh旋转域真实转动交界面处网格相对滑动瞬态求解用于必须捕捉动静件瞬时干涉的工况。动网格Dynamic Mesh边界真实运动网格在计算中光顺、层铺或重划用于阀芯运动、活塞压缩、结构变形等场景。这三者的计算成本递增物理真实性也递增但并不意味着每次都选最真实的方案。有个做离心风机性能曲线的实际项目用MRF稳态算出来的风压、效率跟实验值非常贴近而一上来就用滑移网格瞬态计算单算例耗时翻了近十倍最终结果精度提升却有限。这类工程取舍教材不会替你决定但教材中关于MRF适用于不考虑转子与定子瞬态干涉的稳态性能计算这句话就是判断依据。3.2 交界面设置里的暗坑使用MRF或滑移网格时旋转域与静止域之间必须建立interface交界面。这个步骤看似简单实际上出错率很高。常见问题有三种。第一种是两个交界面网格节点不匹配软件会做插值传递但插值会带来数值耗散可能影响涡结构传播第二种是把interface设成了wall导致流场完全被隔断结果出口流量为零这种错误一般通过查看速度矢量图或质量流量报告就能抓出来第三种是滑移网格计算时没有在交界面上选择正确的帧运动选项导致旋转域转了一阵后界面上的通量出现振荡。我的排查习惯是每次设置完interface先跑一个很短的迭代步然后在后处理里画交界面两侧的速度云图观察是否有不连续的跳变。分界面不光滑就先停下手里的计算把交界面设置检查一遍而不是继续往下算。3.3 动网格参数经验真正使用动网格时三个方法的选择也要心里有数弹簧光顺Spring Smoothing适合小变形动态层Dynamic Layering适合拉伸压缩方向明确的运动局部重划Local Remeshing适合大幅位移、网格质量急剧恶化的情况。动态层方法对网格拓扑有要求运动方向上的网格最好是六面体或棱柱层弹簧光顺则对三角形、四面体网格更友好。我在做阀门开启过程的动网格模拟时起初用弹簧光顺算到开启行程过半时出现负体积导致迭代发散后来改成局部重划加弹簧光顺的组合方案设置好最小/最大网格尺度后问题就解决了。负体积是动网格最典型的杀手锏一旦出现计算几乎必发散。预防措施就两条时间步长宁小勿大让边界每次移动的距离不超过网格尺度的几分之一网格光顺系数保持默认到中等偏小给网格变形留出缓冲空间。4. 收敛不等于结果正确从残差曲线到全局质量平衡很多新手把残差曲线低于1e-3或1e-4当作成功标准这是Part 2阶段需要纠正的重大误区。无论培训教材里有没有强调我都要说残差只是迭代层面的收敛信号结果可信还需要满足物理层面的自洽。4.1 三个层次的收敛缺一不可我的经验是可以把收敛拆成三个层次理解迭代收敛残差曲线下降到设定标准并且不再明显波动。物理收敛监控点出入口流速、壁面平均温度、压差等不再单调漂移呈稳定值或稳定的周期振荡。网格收敛加密网格后关键结果指标的变化幅度在可接受范围内。三个层次都满足才算一份经得起推敲的仿真结果。实际操作中物理收敛往往比残差指标更能反映问题——有时候残差卡在1e-2上不去但你要的出口流量已经稳定不变了这照样可以判断结果可用反之残差漂亮地掉到1e-5但监测的壁面温度持续走高那说明物理过程还没达到稳态继续硬算没有意义。4.2 发散问题的一条完整排查链路真遇到发散时不要慌把发散当成信息来读。下面的排查链是我多年以来一直沿用的顺序看残差曲线中哪个变量先爆。如果是速度残差先爆问题大概率在动量方程稳定性上如果是能量残差先爆多半跟温度场或边界条件设置有关如果是湍流量残差爆往往跟近壁面网格质量或者湍流参数初值不对有关。看监控点曲线是高频振荡还是单调发散。高频振荡常见于时间步长过大或松弛因子过高单调发散通常是网格质量或边界设置出了问题比如存在负体积网格、出口回流设置不当。查看净质量流量失衡imbalance数值。如果进出口质量流量差超过总流量的1%说明连续性方程没有守恒好问题可能出在interface或者边界条件层面。降低欠松弛因子。压力取0.20.3动量取0.40.5先让流场稳定下来再逐步调回来这是成本最低的尝试。检查网格质量。偏斜度Skewness超过0.9的单元在高速流动区域就是隐患需要局部加密或重新划分。这套流程走完绝大多数发散问题都能定位到原因。别一发散就清空设置从头再来那样只会浪费时间并且不知道自己错在哪儿。4.3 提高稳定性的惯用手段除了欠松弛因子几个常用的提高稳定性手段也可以试用伪瞬态方法Pseudo Transient配合耦合式求解器收敛性通常比分离式好初始化时先给一个合理的全域初值比如均匀速度场而不是完全从零开始计算初期先关闭高精度格式用一阶迎风格式把大尺度流场稳定下来再切换到二阶格式。这里要补一句很多教材在讲到格式切换时会用有利于收敛这种模糊表述实际上二阶格式在旋流、激波这类大梯度区域能明显减少数值耗散但也更容易导致残差震荡。我的做法是先用一阶跑两三百步确认全场没有发散迹象再切到二阶继续收敛。这套流程虽然保守但稳定可靠。5. y与湍流模型搭配近壁面网格的尺子怎么用Part 2里涉及湍流模型的内容一定会反复出现y。很多新手对它的理解停留在要很小或者要落在30到300之间却没有真正理解y是一把衡量近壁面网格密度的尺子而且这把尺子的目标区间取决于你选的是哪套湍流模型。5.1 y到底在衡量什么y是一个无量纲量本质上是第一层网格质心到壁面的距离基于壁面摩擦速度做的无量纲化。它衡量的是你的第一层网格抓住了粘性底层的哪个位置。湍流边界层从壁面往外可以大致划分为粘性底层y 5、过渡层5 y 30和对数律层y 30。如果你打算用壁面函数来模拟近壁流动第一层网格要放在对数律层里也就是y处于30到300之间如果你打算解析粘性底层第一层网格就要薄到y接近1。这就是为什么y要取多少没有统一答案它跟湍流模型的选择是绑定的。5.2 第一层网格高度估算从目标y反推实际建模时我们只能通过估算来设定第一层网格高度再用计算结果反查实际y是否达标。估算公式是壁面剪切应力τ_w 0.5 × ρ × U² × C_f摩擦速度u_τ sqrt(τ_w / ρ)第一层网格高度Δy y × μ / (ρ × u_τ)其中C_f可以用平板湍流经验公式估算C_f 0.058 × Re^(-0.2)。举个实际例子。空气外掠平板来流速度U10 m/s特征长度L0.5 m空气密度ρ1.225 kg/m³动力黏度μ1.789×10⁻⁵ Pa·s。雷诺数Re ≈ 3.4×10⁵。代入经验公式C_f 0.058 × (3.4×10⁵)^(-0.2) ≈ 0.0048τ_w 0.5 × 1.225 × 10² × 0.0048 ≈ 0.29 Pau_τ sqrt(0.29 / 1.225) ≈ 0.49 m/s如果目标是y ≈ 1Δy 1 × 1.789×10⁻⁵ / (1.225 × 0.49) ≈ 3.0×10⁻⁵ m也就是约0.03毫米。如果目标改成y ≈ 30Δy ≈ 0.9毫米。两者差了三十倍网格层数可能差出一大截计算量差异非常明显。这就是为什么工程上做高雷诺数外流场时很少有人真的把全局第一层网格做到y1——成本太高。通常的做法是用标准壁面函数配合k-epsilon模型让y落在30到300之间只有在对壁面传热或分离流动预测要求高的地方才局部加密到y≈1并用SST k-omega这类模型。5.3 不同湍流模型与y要求的对应关系下面的对应关系建议保存下来湍流模型近壁处理方式推荐y范围k-epsilon 标准/RNG/Realizable标准壁面函数30300k-epsilon 配合增强壁面处理增强壁面处理局部y≈1k-omega / SST k-omega直接积分到壁面y≈1大涡模拟LES须解析粘性底层y≈1壁面模型LES可放宽需要提醒的是SST k-omega虽然对分离流动预测得比较准但在壁面网格达不到y≈1时计算精度可能比k-epsilon加壁面函数还差。曾有个管道内流动案例我用粗糙网格配合SST k-omega算摩擦阻力结果跟实验差了百分之二十以上后来老老实实加密近壁网格到y1计算结果才回落正常范围。网格和模型是配套的拆开谈任何一个都没有意义。5.4 后处理中怎么验证y计算完成后在后处理里创建Wall类型的面显示y云图检查实际y是否落在自己设定的目标区间内。如果目标区间是30300却发现大面积的y只有十几或几百上千就要调整边界层网格高度重新计算。这里有一条个人心得y的检查不要只看全局最大值和最小值要看主流区域和关键壁面如叶片表面、管道壁面的分布是否均匀。若只在局部角落出现极高的y往往对应的是流动滞止或强加速区域这是物理现象而非网格问题若整个壁面都超出目标范围那就要大刀阔斧地改网格。写在最后的一点体会教材Part 2这一册信息密度比第一册大不少心态上不要指望一遍就全部消化。我自己当年学到这里也是反反复复奔跑在残差曲线和网格质量之间真正把这些模型用顺是在做了三四个实际工程项目之后的事。如果只能给出一条建议那就是每做完一个案例把物理问题—模型选择—关键参数—结果表现四条信息记在一张小表里一个月后回看你会发现自己对FLUENT这套软件的理解上升了整整一个层次。建模永远从物理判断开始点击鼠标只是最后一步。本文还有配套的精品资源点击获取
RELATED READING

延伸阅读

更多一线实战笔记与深度复盘,助您持续精进