流体流动仿真：从层流到湍流的完整分析流程与经验

流体流动仿真(CFD)是工程仿真的重要分支，涵盖层流、湍流、可压与不可压流动等多种形式。结合多个CFD项目经验，聊聊流体仿真的核心思路与实操细节。

1. 流体仿真的分类

流体流动问题按雷诺数可以分成层流(Re < 2300.通常)和湍流(Re > 4000.通常)。层流的数值模拟相对简单，N-S方程可以直接求解;湍流的数值模拟则复杂得多，需要用湍流模型(RANS方法)或者在大涡模拟(LES)的框架下求解滤波后的N-S方程。

我刚接触CFD的时候，以为所有流动问题都需要用湍流模型。后来做了一个低速、小尺寸的微流控通道仿真，雷诺数只有几十，属于典型的层流，用层流模型(Laminar)就足够了。用了湍流模型反而引入了不必要的模型误差。这个经历让我意识到：在打开软件设置模型之前，先估算一下雷诺数，判断流动状态，是很有必要的。

2. 前处理：流体域抽取是第一步

固体力学仿真的建模对象是结构本身;流体流动仿真的建模对象是流体域——也就是流体所占的空间。如果仿真的是一个管道内的流动，需要抽取的是管道内部空间，而不是管道壁面。这个”流体域抽取”(Flow Domain Extraction)是前处理里最容易让新手卡住的一步。

我用ANSYS的SpaceClaim做流体域抽取，步骤是：先导入包含壁面的几何模型 → 用”外壳”工具生成一个包围流体域的实体 → 用”布尔运算”从外壳中减去固体区域 → 剩下的就是流体域。这个过程说起来简单，但如果几何模型里有小缝隙或者面不闭合，布尔运算会失败。我的经验是：做流体域抽取之前，一定先用几何修复工具把所有缝隙补好，面与面之间要完全闭合。

3. 边界层网格：近壁区流动解析的关键

湍流问题中，速度梯度在壁面附近最大。如果壁面附近的网格太疏，近壁区的速度分布会被严重抹平，摩擦阻力和传热系数的计算就会出现较大误差。解决办法是：在壁面附近生成边界层网格(Boundary Layer Mesh)，也就是在紧贴壁面的区域布置多层很薄的网格，厚度方向按一定增长率向外扩展。

边界层网格的设置里，最关键的两个参数是：第一层网格高度(用y+值来控制)和边界层层数。对于k-omega SST这类低雷诺数湍流模型，y+应该接近1(即第一层网格的中心距离壁面的无量纲距离约为1)，这样才能解析黏性底层。我做外流场仿真时，先用经验公式估算一个y+值，生成网格后在Fluent里检查实际的y+分布，如果偏离目标太远，就重新调整第一层网格高度——这个”预估→计算→检查→调整”的迭代过程，是每个CFD工程师的必经之路。

4. 压力基 vs 密度基求解器

Fluent提供了两大类求解器：压力基(Pressure-Based)和密度基(Density-Based)。压力基求解器适用于不可压和低马赫数流动(Ma < 0.3)，是目前使用最广泛的求解器;密度基求解器适用于可压缩流动(高亚音速、跨音速、超音速)，通过同时求解质量、动量和能量守恒方程来处理密度变化。

我的选择标准：马赫数低于0.3 → 压力基，没错。马赫数超过0.5 → 密度基，更合适。介于0.3和0.5之间 → 两个都可以试试，看哪个收敛更快、结果更合理。我在做一个高速列车的气动仿真时，列车时速380 km/h，对应马赫数约0.31.属于那个模糊区间。分别用两种求解器跑了一次，结果差异很小，但密度基的收敛速度明显更快——这主要是因为密度基求解器对可压缩效应的处理更直接。

5. 后处理：不只是看速度云图

流体仿真的结果后处理，很多人停留在”看速度云图和压力云图”的层面。但实际上，CFD可以输出的信息远不止这些。流线(Streamlines)可以直观展示流动路径和回流区;矢量图(Vector Plot)可以展示局部的速度方向和大小;Q准则(Q-Criterion)可以识别和可视化湍流中的涡结构——这对于理解复杂流动中的能量耗散和混合过程非常有帮助。

我在做一个建筑群风环境分析时，速度云图显示某区域风速较低，初步判断是”无风区”。但画了流线之后发现，这个区域其实有显著的涡旋流动，只是速度幅值较低。这个发现改变了对该区域通风效果的判断——涡旋虽然速度低，但能有效促进污染物扩散，通风效果不一定差。这个案例让我意识到：单一视角的结果呈现可能会误导判断，多角度分析才是严谨的做法。