边界层 (Boundary Layer) — CAE术语解释
边界层是什么
我经常听到"边界层"这个术语，但到底是什么意思呢？应该是指壁面附近的某种流动现象吧？
简单地说，边界层是指壁面附近流体流速从零（贴在壁面上）急剧变化到主流速度的薄层区域。在这一层内，粘性的影响非常大，壁面的阻力和热传递几乎完全由这个区域决定。比如飞机机翼的表面摩擦阻力，就是由边界层内的速度梯度产生的。
为什么壁面处的流速是零呢？
这就是无滑移边界条件（no-slip condition），是粘性流体的基本原理。从分子层面看，贴在壁面上的流体拖拽相邻的流体层——这种链式作用导致流速随离壁距离逐步恢复。这个恢复过程所经历的区域就是边界层。这个概念是由Prandtl在1904年提出的，是流体力学教科书中最重要的主题之一。
边界层的厚度 $\delta$ 通常定义为主流速度 $U_\infty$ 的99%所对应的壁面距离。对于平板上的层流边界层，可以用Blasius解求得：
$$\delta \approx \frac{5.0\,x}{\sqrt{Re_x}}, \qquad Re_x = \frac{U_\infty\,x}{\nu}$$其中 $x$ 是从前缘的距离，$\nu$ 是运动粘性系数。
层流边界层和乱流边界层
边界层也分层流和乱流两种吗？它们有什么区别？
区别很大。层流边界层中流体以有序的层状形式流动，速度分布呈光滑的抛物线形。厚度增长缓慢，所以壁面摩擦也较小。而乱流边界层则不同——除了靠近壁面的粘性底层外，其他地方都是激烈的乱流混合，速度分布要"胖"得多——在靠近壁面的地方速度上升得很陡峭。
那乱流边界层的壁面摩擦更大吗？在实际工程中哪种更常见？
没错，乱流的摩擦阻力更大。但工程应用中的流动Reynolds数都很高，所以大多数情况下我们处理的是乱流边界层。比如汽车表面的Re在10⁶~10⁷的量级，除了前缘非常短的距离外，几乎全部是乱流边界层。反之，对于微流道或MEMS器件，Re在1~100的量级，可能始终保持层流。
平板上乱流边界层厚度的经验公式（基于1/7次幂律）为：
$$\delta \approx \frac{0.37\,x}{Re_x^{1/5}}$$乱流边界层的内部结构根据无次元距离 $y^+$ 分为三个区域：
- 粘性底层（Viscous sublayer）: $y^+ < 5$ — 粘性应力占主导。速度呈线性分布 $u^+ = y^+$。
- 缓冲层（Buffer layer）: $5 < y^+ < 30$ — 粘性应力和乱流应力相互竞争的过渡区域。
- 对数律区（Log-law region）: $30 < y^+ < 300\text{--}500$ — 乱流应力占主导。速度分布遵循对数律。
y⁺（y-plus）的基础
y⁺在CFD设置中经常出现，能准确解释一下它的定义吗？
y⁺是用摩擦速度对壁面距离进行无次元化的量。定义如下：
$y$ 是壁面距离，$u_\tau$ 是摩擦速度，$\tau_w$ 是壁面剪切应力，$\rho$ 是密度，$\nu$ 是运动粘性系数。在CFD中，"第一层单元中心的 $y^+$ 值是多少"是最重要的参数，直接决定了乱流模型的选择和网格设计。
具体y⁺应该设置为多少呢？我经常听说"瞄准y⁺=1"……
这取决于所用的乱流模型。总结一下就是这样：
- 使用壁函数的情况（Standard k-ε等）: $y^+ = 30\text{--}300$。因为要跳过粘性底层，反过来说 $y^+$ 太小的话壁函数的假设就会破裂，导致误差。
- 壁面分解的情况（SST k-ω、Spalart-Allmaras等）: $y^+ \approx 1$。用网格直接求解粘性底层。
实际工作中最常用的SST k-ω模型推荐 $y^+ \approx 1$，但很多求解器有自动切换壁函数的功能，只要 $y^+ = 1\text{--}5$ 左右就基本没问题。
壁函数和壁面分解方法
壁函数就是说可以在壁面附近用粗网格吗？
理解得有对有不对。壁函数是通过对数律"跳过"粘性底层的手法，所以把第一层单元放在对数律区域（$y^+ = 30\text{--}300$）就能大幅降低计算成本。但这里的"粗"只是指壁法线方向的第一层，沿着壁面的方向（流向和周向）还是需要充分分辨的。
有什么情况是壁函数不能用的吗？
有的，而且情况不少。壁函数的前提是"壁面附近的流动遵循对数律"，但在以下情况这个假设就不成立，需要壁面分解：
- 流动分离和再附着：逆压力梯度很强的区域，对数律不再适用。比如汽车后窗或机翼失速区域。
- 强加速或减速：喷嘴内的加速流或扩散器的减速流。
- 需要精确计算热传递：如电子器件冷却或燃气轮机叶片的温度预测，壁面附近的温度梯度很关键。
- 包含层流向乱流的过渡：需要捕捉过渡区域时。
在这些情况下，还是老老实实地做 $y^+ \approx 1$ 的网格吧。虽然计算成本上升，但总比用壁函数得到错误结果好得多。
边界层网格的实践
膨胀层应该怎么设置？第一层厚度的计算方法有点搞不懂。
首先根据目标 $y^+$ 反算第一层厚度 $\Delta s$。需要先从流动条件估算壁面剪切应力，用平板近似式比较实用：
比如空气（$\rho=1.225$ kg/m³、$\mu=1.8\times10^{-5}$ Pa·s）、$U_\infty=30$ m/s、代表长度 $L=1$ m、目标 $y^+=1$ 的话，$\Delta s \approx 0.015$ mm 左右。
增长率和层数怎么确定？我平时都是随便设成1.2、10层……
虽然这样做的工程师不少（笑），但还是有规范的。增长率（Growth Rate）推荐 1.1~1.2。超过1.3的话单元质量会急速下降，导致数值扩散增大、收敛变差。层数的话，膨胀层最外层应该超过边界层厚度 $\delta$，大概 15~25层 是目安。
分析完成后有办法确认y⁺是否符合预期吗？
相关术语
- Reynolds数 — 惯性力与粘性力的比值。与边界层发展和过渡直接相关。
- y⁺ — 壁面无次元距离。边界层网格设计的最重要参数。
- 壁函数 — 用经验公式近似粘性底层的壁面建模手法。
- 乱流模型 — 用于闭合RANS方程的模型组。与边界层处理密切相关。
- 膨胀层 — 沿壁面法线方向铺设薄层状网格的手法。
- 壁律（Law of the Wall） — 壁面附近速度分布的通用对数律。
- 分离 — 逆压力梯度导致边界层离开壁面的现象。
CAE术语的准确理解是团队沟通的基础。— Project NovaSolver也着眼于实务者的学习支援。
在边界层网格设置上有困扰吗？
Project NovaSolver以解决CAE工程师日常面临的课题为目标——网格设计的反复试验、y⁺的调整、壁函数选择的判断——我们的目标是实现这些课题的解决。你的实务经验将成为更好工具开发的原动力。
咨询（准备中）相关主题
价值
详细
错误