简化模型：流体仿真为什么卡顿，如何提速？
1 绪论
从事产品设计开发的一线工程师，日常始终在和「时间」「计算资源」激烈博弈。为快速推出性能完善的产品，需要反复迭代设计方案并完成验证；计算流体力学（CFD）能够大幅削减实体样机制作成本与测试周期，现已深度落地于各行业研发工作。
左侧：精细网格复杂几何模型；右侧：简化模型；示意图提示：若完整建模螺栓、销钉、风机等全部结构，会造成仿真计算严重卡顿
但现场工程师常会陷入矛盾困境：为精准捕捉流动特征，把图纸里所有细小结构（螺栓、销钉、风机叶片）全部纳入仿真模型，最终建模、等待计算完成消耗的时间，不亚于实体物理测试，直接背离了引入 CFD 实现「快速敲定设计方向」的初衷，形成仿真瓶颈。
1.1 高成本流体仿真概述
不做任何数学假设、1:1 完整还原产品全部几何外形与运动状态的直接模拟法，会产生极高数值计算成本，甚至直接导致项目无法推进。本章将识别高资源消耗的物理工况，介绍工程简化方案，实现低成本、高可靠仿真。
1.1.1 数值计算成本
CFD 仿真成本不单纯由计算机硬件性能决定，数值算法与物理工况复杂度叠加，共同决定整体计算开销。
矩阵规模与自由度（DOF）的关联
CFD 通过离散纳维 - 斯托克斯方程求解巨型线性方程组；网格数量翻倍时，矩阵规模呈几何级膨胀，内存占用、CPU 计算量非线性暴涨。尤其数千万网格的模型，计算数据超出 CPU 运算缓存上限，受内存带宽限制，计算效率大幅下滑。 时间尺度约束与 CFL 条件（库朗 - 弗里德里希斯 - 莱维条件）
若模型包含运动结构，必须采用瞬态仿真；为保证数值稳定，时间步长 Δt 必须设置极小，几何越复杂、网格越细密，步长取值越小。仅模拟数秒物理流动，就需要数十万次迭代计算，大幅抬高仿真成本。
1.1.2 高资源消耗物理现象案例
建模前需先识别以下工况，提前制定简化方案：
运动边界与相对运动
叶轮、往复活塞等运动结构，网格需要随结构变形或采用重叠网格，每一步都额外增加物理量插值运算；相比纯静止结构仿真，算力消耗提升 10~50 倍以上。 多尺度几何结构
大型机壳（米级）内部共存微米级过滤微孔；若为捕捉微细流道压降，全域统一加密网格，会造成网格总量爆炸，仿真周期急剧拉长。 流体机械系统仿真
多风机强制向系统输入能量，叶片尖端会产生强涡流与流动分离；完整模拟该现象需要壁面附近极致加密网格，加重湍流模型计算负担，收敛速度显著下降。
1.1.3 简化建模的降本收益
将高消耗物理工况替换为数学模型（MRF、多孔介质、风机边界等），可获得如下直接收益：
仿真类型转换（瞬态→稳态）
最显著的降本手段：瞬态仿真转为稳态。例如采用 MRF 方法，可将旋转件动态行为转化为静态方程描述，整体仿真时长缩短 90% 以上。 提升数值鲁棒性、加速收敛
复杂几何带来的网格扭曲、长宽比失衡，是残差震荡的核心诱因；简化模型从根源消除这类数值不稳定因素，大幅降低仿真发散风险。
*工程师小贴士：简化模型应用判定标准
简化不等于无差别删除结构，通过 3 个问题判断适用范围：
a.是否属于核心关注区域？ 若研究风机自身振动，需精细建模；若仅分析系统整体流量，选用风机边界条件即可。 b.能否保留流动核心特性？ 使用多孔模型时，必须有实验数据支撑，保证模型复现与真实过滤器一致的压降曲线。 c.计算资源是否存在上限？ 结合可用硬件、项目交付周期，找到精度与计算成本平衡的最优解，是一线工程师核心能力。
1.2 用简化模型替代复杂流动现象
简化建模的核心不是单纯删除几何，而是把结构对流动的影响等效为数学函数。计算机无需处理数十亿复杂几何网格，仅求解简化代数方程，仿真速度大幅提升。
1.2.1 完整复现真实物理场景难度大的原因
a.微细流道造成网格激增
仿真精度由网格分辨率决定；若要捕捉 1mm 小孔内流速梯度，孔内至少划分 5~10 层网格，设备整体尺寸越大，总网格数呈几何倍数上涨。 b.网格畸变导致收敛变差
强行对复杂几何划分网格，会大幅提升网格非正交性；劣质网格会降低单节点矩阵收敛性，迭代时残差持续震荡，最终仿真发散。
1.2.2 旋转运动简化：MRF 与滑移网格
旋转件仿真核心思路：是直接让结构转动，还是等效施加相对速度。
相对速度等效模拟旋转流动
MRF（移动参考坐标系）不物理旋转叶片，而是对指定区域引入旋转坐标系。
附加惯性力项
不真实转动几何，而是在动量方程中增加虚拟惯性力：离心力、科里奥利加速度，仅需极低计算成本即可精准复现旋转流动。
1.2.3 多孔介质流动简化：多孔介质模拟方法
过滤器、催化剂等重复微细流道结构，不单独计算每条流道，而是从整体体积阻力角度建模。
体积阻力模型原理
删除复杂内部几何，对该区域流体施加与流速成正比的粘性阻力、与流速平方成正比的惯性阻力（动量损失）。
各向异性调控
仅通过数值参数，无需划分网格，即可实现定向通流、侧向阻流的流动规整效果。
1.2.4 流体机械系统简化：风机边界条件法
风机、叶片是向流体输入能量的部件，简化模型中将其等效为「能量阶跃」数学模型。
精细风机网格模型→简化单一面替代风机结构
压力阶跃面原理
三维叶片结构简化为单个二维面；流体穿过该面时，沿流动方向产生不连续压力抬升。
自动求解工作点
输入风机 P-Q 性能曲线后，CFD 求解器自动匹配系统总压损与风机升压，求解稳态工作流量。
2 旋转流体机械仿真方法
旋转流体机械依靠叶轮、叶片物理旋转，向流体传递能量（泵、风机）或回收流体能量（涡轮）。从 CFD 角度，旋转部件仿真相比常规静止结构（管路、外流）存在独特难点：如何数值处理运动边界与静止域之间的能量、物理量交互。
2.1 移动参考坐标系（MRF）
MRF 不物理旋转叶轮，而是给旋转域流体施加反向速度分量，等效复现旋转效果；原理和风洞固定车辆、吹风模拟行驶工况一致。
2.1.1 数值实现原理与理论
midas NFX CFD 求解器在旋转坐标系下描述流动，将绝对速度
与相对速度
关系式引入控制方程：速度关系式：
修正动量方程
旋转域属于非惯性坐标系，无需真实转动几何，通过增加两项惯性力复现旋转效应：
推动流体远离旋转中心，生成压力梯度；
垂直于流体相对速度施加惯性力，精准模拟旋流2.1.2 midas NFX 单元连接方式：节点合并
MRF 前提：旋转域与静止域相对位置不随时间变化（稳态）。实操设置：midas NFX 前处理中，转子（旋转区）与定子（静止区）交界面使用节点合并，节点完全共享。优势：节点完全重合，跨界面质量、动量守恒完美满足，无插值带来的数值误差。
2.1.3 稳态仿真效率与实操注意事项
MRF 多用于快速求解旋转设备稳态性能（扬程、效率）。
优点：无需每步重构网格，计算速度快，收敛判定清晰； 局限：几何位置全程固定，无法实时捕捉叶片与壳体周期性干涉效应。
工程师实操提示：MRF 叶轮摆放建议
若叶片与壳体凸起（蜗舌）距离过近并固定，会产生局部非物理高压；建模时建议将叶片调整至不会干扰主流场的相位。
2.2 滑移网格（Sliding Mesh）
滑移网格将叶轮周边网格与静止域分离，随时间步真实转动网格，是高精度瞬态仿真方案。
2.2.1 物理量传递机制与控制方程
每一时间步网格发生位移，采用基于移动控制体积的积分型控制方程：
：网格移动速度
：流体速度非共形界面：转子与定子交界面节点不重合；CFD 求解器通过守恒型插值，实现界面间无损耗传递质量、动量。
2.2.2 midas NFX 单元连接方式：邻面匹配
滑移网格网格持续滑动，绝对不能使用固定节点合并。实操设置：划分网格时使用「邻面匹配（接触条件）」生成非共形界面。优势：转子、定子节点相互独立，可精准捕捉叶轮旋转时叶片与壳体的动态干涉（RSI）。
2.2.3 瞬态仿真精度控制
若需分析外部结构耦合影响，必须采用滑移网格：
转子 - 定子干涉（RSI）
：叶轮扫过蜗壳、扩压器等固定结构产生的周期性压力脉动； 叶片通过频率（BPF）噪声分析
3 阻力件流动简化：多孔介质与边界条件
多孔模型删除复杂内部几何，仅通过数学公式复现流体穿过该区域产生的压降 ΔP；根据结构厚度分为两类：三维体域「多孔介质」、二维面域「多孔阶跃边界 Porous Jump」
3.1 多孔介质（三维体域建模）
适用于厚型过滤器、催化层、散热翅片等整体体积持续产生流动阻力的结构。
3.1.1 理论基础：达西 - 福希海默方程
CFD 在动量方程中增加虚拟阻力源项
实现多孔效果：
决定设置路径：【属性】菜单，对流体区域输入多孔介质参数。网格要求：必须节点合并，源项需要在网格节点连续计算，多孔域与周边流体域交界面网格节点完全连通。
3.2 薄型阻力件简化：多孔阶跃边界（Porous Jump）
厚度仅 1~2mm 的穿孔板、细密金属网，三维建模会生成海量网格，效率极低，选用多孔阶跃边界。
3.2.1 二维不连续面模拟压降原理
将阻力件等效为无厚度二维面，流体穿过时上下游产生数值压力阶跃：
𝚫𝒎：阻力件实际厚度
3.2.2 midas NFX 实操设置
设置路径：【边界条件】，输入阻力系数与厚度 Δm。网格要求：禁止节点合并，使用邻面匹配；上下游节点分离，才能正常计算两面压差。
3.3 各向异性多孔模型调控流动
多孔模型最大优势：分方向设置不同阻力系数，实现各向异性流动控制。
Honeycomb蜂窝结构沿孔道方向通流顺畅，侧向几乎无流动；无需精细建模实体，主流向设低阻力，另外两垂直方向阻力放大 100~1000 倍，低成本实现流动整流效果。
工程师避坑指南：多孔模型常见错误
单位换算失误：
不可直接输入整体压降数据，需将实验 ΔP 除以厚度，换算为单位长度压力损失；
孔隙率误用：
孔隙率用于缩小有效流通截面积、提升局部流速，参数错误会严重干扰对流换热计算，需谨慎填写。
3.4 多孔模型选型汇总
4 风机边界条件 Fan Boundary Condition
删除复杂三维叶轮几何，等效为向流体输入动能的虚拟二维面。
4.1 数值原理与工作点自动求解
能量输入模型：在风机所在面添加动量源项，流体穿过该面时，叶片做功转化为数值压力抬升 ΔP。
自动求解工作点：不同于固定流量入口，风机边界每轮迭代自动求解系统阻力曲线与风机性能曲线交点，确定实际稳定流量。
4.2 风机 P-Q 性能曲线耦合逻辑
核心输入参数为厂商提供的流量 - 压力（P-Q）曲线：CFD 实时读取当前截面流量
，匹配曲线对应升压值，施加于截面两侧
4.3 风机安装位置对应的 midas NFX 设置策略
风机安装位置不同，系统边界与网格处理完全不同，需提前规划。
4.3.1 外置风机（入口 / 出口风机）
设置在计算域外边界（与大气接触处），耦合大气压与风机性能曲线计算。
4.3.2 内置风机
安装在系统管路中间，等效为内部薄壁不连续面。
网格要求：必须使用邻面匹配，禁止节点合并，保证截面压力不连续；
流动方向控制：在边界条件中手动指定方向矢量，矢量反向会使风机变为阻流壁面，务必核对。
工程师提速技巧
a.增设虚拟直管段：
风机前后有障碍物易产生数值震荡；在风机面上下游增加 1~2 倍风机直径的直管段规整流场，大幅提升收敛速度；
b.分步初始化流量：
初期用固定流量条件建立稳定流场，迭代中期切换风机边界自动寻优，避免初始逆流发散；
c.区分静压 / 全压：
厂商曲线需确认基准为静压，绝大多数 CFD 求解器仅支持静压曲线，混淆会导致流量严重失真。
5 不同工况最优湍流模型选型策略
简化模型会引入剧烈压力阶跃、强旋流场；简化方法自带数学源项，湍流模型对涡、耗散的近似方式直接决定结果精度，必须针对性选型。
5.1 主流湍流模型特性与选型原则
k-ε 系列：
主流区自由流仿真表现优秀，数值收敛性极强；适合分析系统宏观压降、整体流动趋势，计算成本低。
k-ω SST Shear Stress Transport：
精准捕捉近壁流动分离、逆压梯度；精细设计阶段分析叶片表面流动、阻力件尾迹涡流时推荐使用。
5.2 简化方法配套湍流模型匹配指南
5.2.1 旋转流体机械（MRF / 滑移网格）
推荐模型：k-ω SST
标准 k-ε 会过度预测湍流粘度，抑制旋流强度，产生非物理误差；滑移网格需精准捕捉叶片表面分离、压力脉动，优先选用 SST。
5.2.2 多孔介质 / 多孔阶跃边界
常规系统仿真推荐：k-ε
多孔区域内部能量损失已全部纳入达西 - 福希海默源项，无需精细求解内部湍流；追求全域压力平衡时，低开销、稳定的 k-ε 性价比更高；仅当需要精确分析过滤器下游扩散尾迹时，选用 SST。
5.3 风机边界条件配套湍流模型
推荐模型：k-ε 系列
风机边界多用于大型系统（机房、中央空调）全域流量、工作点预测；k-ε 收敛稳定，能快速求解风机性能与系统阻力平衡点。
实操排障指南
a.残差震荡、数值不稳定：SST 对壁面第一层网格高度 y + 敏感；可先用 k-ε 计算稳定初始流场，再切换 SST 续算； b.湍流粘度比超限警告：多为风机 / MRF 界面网格尺寸突变；平缓调整网格梯度比例即可； c.模型参数优先级高于湍流模型：多孔阻力系数、风机 P-Q 曲线数据存在误差时，再高精度湍流模型也无法得到可靠结果，优先校验输入实验数据。
5.4 湍流模型选型汇总表
6 结论：高效仿真完整方案
CFD 仿真的核心目标不是视觉上 1:1 复刻实体几何，而是用最低算力成本，输出支撑设计决策所需物理量。本文介绍多种几何简化方法，通过数学函数替代复杂结构，平衡精度与计算效率。
6.1 核心技术要点回顾
6.1.1 旋转部件仿真
MRF 移动参考坐标系：瞬态转稳态，节点合并保证物理连续，快速求解设备整体性能； 滑移网格：真实旋转网格，邻面匹配传递物理量，用于压力脉动、叶片通过频率噪声分析。
6.1.2 阻力件多孔仿真
三维多孔介质：厚过滤 / 催化层，属性面板赋值，必须节点合并； 二维多孔阶跃边界：薄板、金属网，边界条件赋值，邻面匹配分隔网格。
6.1.3 风机边界条件
去除叶片几何，输入 P-Q 曲线自动求解工作点；内置风机必须邻面匹配，实现截面压力不连续。
6.1.4 数值稳定性（湍流模型）
全域流量、压降分析选用收敛性更好的 k-ε；叶片分离、涡流精细分析选用 k-ω SST。
6.2 简化方法选型总览表
6.3 工程师寻找最优平衡点三大准则
a.输入数据质量决定仿真结果
遵循「垃圾进、垃圾出」原则；简化模型高度依赖实验压降、风机性能曲线，仿真前优先核验厂商参数、实验数据的物理合理性。 b.分步迭代仿真策略
不要直接使用滑移网格 + SST 高精度模型起步；先用 MRF+k-ε 低负载方案建立稳定初始流场，再按需提升精度，避免仿真发散、缩短项目周期。 c.严守网格连接核心规则
midas NFX CFD 核心实操准则：
需要全域物理连续（MRF、三维多孔介质）：节点合并； 仅需数值传递、允许压力不连续（滑移网格、多孔阶跃、内置风机）：邻面匹配；两类操作混淆会直接导致不收敛或完全错误的仿真结果。
