什么是湍流模型与y⁺计算
🙋
CFD模拟里总说要选湍流模型，还要算y⁺，这到底是什么呀？
🎓
简单来说，湍流模型就是一套数学公式，用来描述像狂风、漩涡这种混乱流动的平均效果。而y⁺（y-plus）是一个关键的无量纲数，它决定了你靠近壁面的第一层网格该画多密。在实际工程中，比如设计汽车后视镜时，如果网格太粗，就预测不准风噪和阻力。你试着在模拟器里拖动“参考VelocityU”的滑块，看看y⁺和第一层网格高度Δy是怎么跟着变化的，就能直观感受到它的重要性了。
🙋
诶，真的吗？那为什么有的模型要求y⁺≈1，有的却要30到300呢？差别也太大了吧！
🎓
这取决于模型如何处理壁面附近的流动。像k-ω SST这类“低雷诺数模型”，需要把网格画得足够细，直接去解析粘性底层（就是紧贴壁面那层几乎不动的流体），所以要求y⁺≈1。而k-ε模型配合“标准壁面函数”时，它用一个经验公式跳过了粘性底层，直接从对数律区开始算，所以网格可以粗很多（y⁺=30~300）。你可以在模拟器里选择不同的模型，然后观察“推荐y⁺范围”的变化，再对比一下计算出来的Δy，能差出几十倍呢！
🙋
原来是这样！那“湍流强度”和“长度尺度”这两个输入参数，如果我不知道准确值，随便填一个影响大吗？
🎓
影响非常大！这两个参数定义了入口处湍流的“混乱程度”和“漩涡的平均大小”。比如模拟建筑风环境，来流湍流强度可能高达10%以上；而模拟飞机巡航，可能只有0.1%。如果你乱填，仿真结果可能完全失真。工程现场常见的是，根据风洞实验数据或经验公式来估算。你不妨在模拟器里，先把湍流强度从1%调到10%，看看它如何影响计算出的湍动能k和耗散率ε，你就会明白为什么它被称为关键的“入口边界条件”了。
物理模型与关键公式
最核心的概念是无量纲壁面距离y⁺，它由摩擦速度和第一层网格高度决定：
$$y^+ = \frac{u^* \Delta y}{\nu}$$
其中，$u^* = \sqrt{\tau_w / \rho}$ 是摩擦速度（表征壁面剪切力的大小），$\Delta y$ 是第一层网格高度，$\nu$ 是流体的运动粘度。本工具通过雷诺数$Re$和参考速度$U$，利用平板摩擦系数经验公式估算出$u^*$，从而反推出所需的$\Delta y$。
壁面附近的平均速度分布遵循“壁面定律”，这是设置壁面函数的基础：
$$ \text{粘性底层}(y^+<5): \quad u^+ = y^+ $$
$$ \text{对数律区}(y^+>30): \quad u^+ = \frac{1}{\kappa}\ln(y^+) + B \quad (\kappa \approx 0.41, B \approx 5.2)$$
这里 $u^+ = U / u^*$ 是无量纲速度。标准壁面函数就是基于对数律公式，将网格中心点（位于对数律区）的速度与壁面剪切力联系起来，从而避免直接求解粘性底层。
现实世界中的应用
汽车空气动力学：在汽车风阻和升力仿真中，k-ω SST模型配合y⁺≈1的精细网格是行业标准，它能更准确地预测车身后部的气流分离和涡旋，这对优化燃油经济性和高速稳定性至关重要。
航空航天：对于飞机机翼绕流，同样广泛使用SST模型。而在模拟发动机内部高温高压的复杂流动时，可能会采用更高级的尺度自适应模拟（SAS）或分离涡模拟（DES），以捕捉非定常的涡脱落现象。
建筑与风工程：评估摩天大楼的风荷载和行人区风环境时，常使用k-ε Realizable模型配合标准壁面函数（y⁺=30~300）。这种方法在保证精度的同时，能大幅降低网格量和计算成本，适合处理复杂的城市几何。
叶轮机械：设计水泵、风机或涡轮叶片时，流动分离和效率预测是关键。k-ω SST模型因其在逆压力梯度流动中的优异表现而被优先选用，帮助工程师优化叶片型线，提升效率并抑制空化。
常见误解与注意事项
首先，存在一种“只要y+合适就万事大吉”的误解。确实y+是壁面网格划分的第一步，但仅此是不够的。例如，即使将y+设为≈1，如果壁面法向方向的网格层数过少（例如少于5层），就无法解析整个边界层的速度梯度，导致对流动分离的预测错误。通常建议使用10至15层以上的网格来覆盖整个边界层。
其次，误认为“推荐的y+值是绝对的金科玉律”。工具提示的“k-ε模型建议y+≈30”仅仅是一个起点。如果实际流场复杂（存在强压力梯度、流动分离与再附着等情况），有时采用更小的y+值（例如10~20）可能更合适。务必检查模拟结果的壁面剪切应力与压力分布，并根据需要灵活调整网格，这种灵活性至关重要。
最后，随意设置入口边界条件中的“湍流强度”与“长度尺度”。这些参数直接关系到计算的初始值，尤其会影响流动的发展。例如，若要模拟风洞实验，就应使用实验值。若无数据，对于内部流动（如管道）常以百分之几的强度开始，外部流动（如汽车绕流）则常从0.5%~1%左右起步，但切忌“直接采用默认值”，通过敏感性分析来确认其影响才是专业态度的体现。