Y+该小于5还是大于30？搞懂湍流能量级串，边界层网格就不会白画

Y+ 判据、稳态瞬态怎么选、能量级串三个分区、RANS 为什么算不准——一篇打通。

很多人画边界层网格是靠「抄别人的层数」，算完 y+ 不合格还不自知；也有人明明该用瞬态，硬拿稳态糊过去，结果和实测差一大截。

根子都在于没理解湍流的能量是怎么产生、怎么传递、又怎么死掉的。这篇把 STAR-CCM+ 课程里最硬的一节理论课拆开讲：网格判据怎么验、稳态瞬态怎么选、能量级串三个分区各管什么、以及相干性和间歇性为什么让 RANS 精度打折。

一、Y+ 判据从哪来：缓冲层是「禁区」

湍流边界层分内层和外层。内层再细分为三层，每层的速度规律完全不同：

粘性底层　U⁺ 与 Y⁺ 成正比（线性）关系

→ 粘性力主导，速度梯度极大

缓冲层 ⚠没有一个特别适用的函数能描述它

→ 粘性力和湍流应力都不占绝对主导，是经验公式的空白区

对数律层　U⁺ 与 Y⁺ 成对数关系

→ 湍流应力主导，壁面函数就建在这一层

所以划网格的铁律是

Y⁺ < 5　或　Y⁺ > 30

要么落进粘性底层（直接解析），要么落进对数律层（用壁面函数）——就是别停在缓冲层里

图1　边界层内层三层：粘性底层（线性）、缓冲层（无适用函数）、对数律层（对数）

而且——第一层 Y+ 合格，不代表网格就合格。检查边界层网格必须同时看两件事：

① 第一层网格高度是否满足 Y+ 要求（落在粘性底层或对数律层）

② 边界层总厚度是否满足速度梯度要求——也就是有没有把整个边界层包在里面

只看第一条，很可能出现「Y+ 完美，但棱柱层在第 10 层就截断了，边界层外面还有一大段速度梯度没被解析」的情况。

📐 平板边界层实例（入口 10 m/s）

· 网格：棱柱层 + 多面体

· 边界层控制方式：第一层网格高度

· 边界层总厚度 0.03 m，共 30 层

· 第一层高度用 Y Plus 计算器算出 = 3.5×10⁻⁵ m

· 物理模型：定常 / 气体 / 分离流 / 恒密度 / 湍流 K-Epsilon（自动选 Realizable K-E）

· 壁面处理：All Y+ Wall Treatment——它能自动适配你的边界层网格，不用纠结首层落在哪一层

图2　物理模型选择：Realizable K-Epsilon + All Y+ 壁面处理

两条判据怎么验

验第一条：建一个 XY Plot，域选平板壁面，左轴数据选 Y+。结果是大部分位置 Y+ < 1，只有起始位置 1.5，全线都在 5 以下 → 第一层网格稳稳落在粘性底层，合格。

图3　Y+ 沿壁面分布：绝大部分 < 1，起始处 1.5，全部满足 < 5

验第二条：看出口速度剖面。横轴是「速度 / 自由流速度」，纵轴是无量纲距离。当速度比达到 1（即恢复到自由流）时，对应的距离是 0.02 m；而边界层网格总厚度是 0.03 m。0.02 < 0.03，说明整个速度梯度都被边界层网格包住了，两条判据同时满足。

图4　出口速度剖面：速度比达 1 时距离 0.02 m，小于边界层总厚 0.03 m

顺带一个直觉：边界层是沿平板逐渐发展的，入口处最薄、出口处最厚。所以用一个固定的首层高度配合足够的总厚度，才能全程覆盖。

二、稳态还是瞬态？只看「时均值」变不变

湍流本身是随机的，速度一直在波动。那为什么还能用稳态算？关键在于时间平均值：

稳态（Steady）

定义：流动特征不随时间变化。湍流速度虽然随时间波动，但时间平均下来的时均值始终恒定 → 就可以用稳态模拟，只关注宏观特性。

目的：简化分析、计算效率高。

适用场景：换热器、压缩机、泵、风机、涡轮的稳定运行阶段。

瞬态（Transient）

定义：流动特征随时间变化。速度波动的同时，时均值本身也呈现上升或下降的趋势。

适用场景：非定常动态过程——泵的启动、冲击、周期性流动、发动机活塞缸；储罐连续填充（往水箱冲水）；往复式机械（活塞泵、内燃机活塞上下运动）；阀门快速启闭。

代价：更贴近真实物理世界，但计算成本高得多。

图5　典型瞬态场景：储罐连续填充、往复式机械、阀门快速启闭

一句话判据：你关心的那个物理量，时均值恒定就用稳态，时均值还在漂就必须上瞬态。

三、能量级串：湍流的能量是怎么「死掉」的

能量级串的定义：能量从大尺度涡向小尺度涡逐步传递的规律。一句话概括就是——大涡孕育小涡，小涡靠大涡的动能维持，小涡又衍生出更小的涡，逐级传递，直到能量因为粘性耗散而消失。

主流 / 泵 / 风机 注入动能↓

大尺度涡（含能区）↓ 涡拉伸，分裂

中等尺度涡（惯性子区）↓ 继续传递

小尺度涡（耗散区）↓ 粘性主导

耗散成内能（热能）

图6　能量级串：红色为大尺度涡（强度高），蓝色为小尺度涡，最终粘性耗散

按能量级串的过程，可以把它分成三个区，每个区各由一个长度尺度来描述：

① 含能区　｜　尺度：积分长度尺度

大涡主导，能量最高，能量从外部（主流、泵、风机）注入。

积分长度尺度与湍流动能的 3/2 次方成正比、与能量耗散率成反比。它是能量注入的主要尺度，反映湍流的宏观流动特性和流动几何尺寸——比如管道直径、机翼弦长。

② 惯性子区　｜　尺度：泰勒微尺度

中等尺度涡，是介于大涡和小涡之间的过渡区域。能量只在涡之间传递，既没有能量产生，也没有耗散成热能——就是大涡拆成小涡的纯粹搬运过程。此时涡开始受到粘性效应影响，泰勒微尺度由运动粘度、湍流动能和能量耗散率共同决定。

③ 耗散区　｜　尺度：科尔莫戈洛夫长度尺度

小涡主导，分子粘性效应显著，这是湍流当中最小尺度的涡。粘性耗散完全主导，能量无法再向更小的涡传递，最终转换成热能散发掉，能量归零。小涡的动能通过粘性摩擦变成内能，完成整个能量级串的闭环。

图7　能谱三分区：含能区 → 惯性子区 → 耗散区

湍流动能 k 是什么

定义：三个方向速度脉动的时间平均之和的一半，即 k = ½( u′² + v′² + w′² ) 的时均。

物理意义：代表所有尺度涡能量的总和，是描述湍流强弱的核心量。

工程意义：网格必须解析到你想关心的那个尺度。只关心含能区（工程上大多数情况），网格按积分长度尺度给就够了；想一路解析到耗散区，那就是 DNS——网格要细到科尔莫戈洛夫尺度，这也是 DNS 在工程上算不起的根本原因。RANS/LES/DNS 的取舍，本质就是「你打算在级串的哪一层接手」。

图8　积分长度尺度 / 泰勒微尺度 / 科尔莫戈洛夫长度尺度，分别对应三个能量分区

四、相干性与间歇性：RANS 为什么会失准

相干性：湍流不是完全随机的

定义：流动当中存在清晰、持续且有序的大尺度漩涡结构，这种结构就叫相干结构。

触发机制＝不稳定性：当两股流体存在速度差或剪切作用时，流动会失稳，产生规则的大尺度漩涡——这些漩涡就是相干结构的核心形态。

常见场景：混合层、射流的剪切层、钝体后的尾流（比如圆柱绕流：涡一会儿在上、一会儿在下，看着乱，其实节奏有序）。

结论：湍流虽然随机，但并非完全随机，它包含着有序的大尺度涡结构。

图9　相干结构：圆柱绕流中清晰、持续、有序的大尺度涡脱落

间歇性：湍流也不是「一直湍」的

定义：流动在时间和空间上呈现不规则、非持续的湍流脉动。它有两层含义：

① 空间点上层流/湍流来回切换　某个固定空间点并非完全处于湍流状态，而是在层流和湍流之间随时间切换。这要用转捩模型来描述，其核心参数就叫间歇因子。

② 完全湍流区里 k 和 ε 也在剧烈波动　即使处在完全湍流区域，湍流动能和耗散率也会随时间剧烈波动，而不是一个恒定值。

⚠ 矛盾就在这里

常规的雷诺平均（RANS）湍流模型，假设的是湍流动能 k 和耗散率 ε 在时空上恒定。

但实际流动中，k 和 ε 随空间和时间都会显著波动，流动状态是随机切换的。理想模型 vs 实际流动的这个差距，最终导致传统 RANS 模型的计算精度下降。

结论：想捕捉间歇性的瞬态结构、拿到更高精度，就要用大涡模拟（LES）或转捩模型。间歇性主要出现在层流-湍流转捩区域和复杂流动区域——这两类工况，就是 RANS 最容易翻车的地方。

图10　间歇性：理想化的 RANS 假设与实际波动的差距，就是精度损失的来源

✅ 一页纸记住这四个结论

1. 网格：Y⁺ < 5 或 > 30，避开缓冲层；而且必须同时验「首层 Y+」和「边界层总厚度是否包住整个速度梯度」

2. 稳态/瞬态：时均值恒定用稳态（换热器、泵、风机的稳定运行段）；时均值还在漂就上瞬态（启动、冲击、往复式、阀门启闭）

3. 能量级串：含能区（大涡·积分长度尺度）→ 惯性子区（只传递·泰勒微尺度）→ 耗散区（小涡·科尔莫戈洛夫尺度）；k = ½Σ速度脉动平方时均

4. 精度：RANS 假设 k、ε 时空恒定，遇到相干结构和间歇性就会失准 → 转捩区、复杂流动区上 LES 或转捩模型

扫码添加小助理，回复关键词「STAR CCM+」获取转行学习路径与课程大纲