参数设置
空气 20°C: ν ≈ 1.5×10⁻⁵ m²/s, ρ ≈ 1.2 kg/m³。水 20°C: ν ≈ 1.0×10⁻⁶ m²/s, ρ ≈ 998 kg/m³。
边界层速度剖面 u(y)/U
平板上方流动 / 绿实线＝Blasius 剖面 u/U / 蓝虚线＝δ_99 / 橙＝δ* / 紫＝θ
理论与主要公式
对于零压力梯度下平板上的层流边界层，Blasius 相似解给出如下代表性结果。
局部雷诺数。U 为主流速度，x 为距前缘的距离，ν 为运动粘度：
$$\mathrm{Re}_x = \frac{U\,x}{\nu}$$
99% 厚度、排挤厚度与动量厚度：
$$\delta_{99} \approx \frac{5.0\,x}{\sqrt{\mathrm{Re}_x}},\quad \delta^{*} \approx \frac{1.721\,x}{\sqrt{\mathrm{Re}_x}},\quad \theta \approx \frac{0.664\,x}{\sqrt{\mathrm{Re}_x}}$$
局部摩擦系数与壁面切应力。ρ 为流体密度：
$$C_f = \frac{0.664}{\sqrt{\mathrm{Re}_x}},\qquad \tau_w = \tfrac{1}{2}\,C_f\,\rho\,U^2$$
形状因子 H = δ*/θ ≈ 2.59。当 Re_x ≳ 5×10⁵ 时层流过渡到湍流，上述公式不再适用。
关于 Blasius 边界层模拟器
🙋
边界层就是物体表面附近流速变慢的那一薄层对吧？Blasius 解为什么这么重要？
🎓
通俗地说，它是历史上第一次把 Prandtl 提出的边界层方程在最简单情形——平板、层流、零压力梯度——下严格解出来的。从此之后翼型设计、涡轮机械都把它当作基本参照。在模拟器里拖动"位置 x"，你会看到 δ_99 按 √x 的规律增长，这就是 Blasius 解的精髓。
🎓
问得好。边界层在数学上是渐近趋向主流的，没有一个明确的"上边界"。所以工程上约定取流速达到主流 99% 的高度作为边界层厚度。Blasius 解给出 δ_99 ≈ 5x/√Re_x。在默认参数（U=1 m/s, x=100 mm, ν=1.5×10⁻⁵ m²/s）下应该是 6.1 mm 左右。
🙋
卡片里还有排挤厚度和动量厚度，它们跟 δ_99 区别是什么？
🎓
可以把 δ_99 当作"看得见的厚度"，把 $\delta^{*}$ 和 $\theta$ 当作"效应厚度"。$\delta^{*}$ 表示边界层把主流挤出去的等效距离——在管道里相当于把有效流通截面减少了 $\delta^{*}$。$\theta$ 是损失的动量被换算成的长度，本质就是壁面摩擦的积分。Blasius 解中 $\delta^{*}/\theta \approx 2.59$，这是判别层流边界层的标志性数值。
🙋
C_f 随 1/√Re_x 下降，是说雷诺数大的时候摩擦就变小了？
🎓
下降的只是"无量纲"的摩擦系数。真实壁面切应力 τ_w = ½ C_f ρ U²。如果把 U 从 1 m/s 提到 5 m/s，C_f 会减少约 √5 倍，但 τ_w 会增大约 11 倍——因为 U² 项远比 C_f 的下降快。评估实际阻力时这两点都要考虑。
🙋
Re_x 一旦超过 5×10⁵，红色标签就亮起来了。Blasius 解会失效吗？
🎓
是的。光滑平板上的自然过渡大约在 Re_x ≈ 5×10⁵ 开始。再往后边界层变成湍流，厚度按 x^(4/5) 增长，而 C_f 远高于层流值（典型为 0.0592/Re_x^(1/5)）。真实机翼、船体等大多在前缘后不远就转为湍流，工程上把前段的 Blasius 公式与后段的湍流相关公式合起来用。模拟器会在超出适用范围时给出警告。
常见问题
在相同位置 x 处比较，层流边界层比湍流薄、摩擦系数也较小，但更容易发生分离。湍流由于强烈混合把主流的高动量带到壁面附近，因此抗分离能力强——这正是高尔夫球凹坑和飞机扰流片人为触发湍流的原因。设计上的取舍是："想减小摩擦的部位保持层流，必须避免分离的部位强制湍流"。
把局部摩擦系数 C_f 从 0 积分到板长 L，得到全层流情形下的平均摩擦系数 C_F = 1.328/√Re_L。摩擦阻力为 D = C_F · ½ρU² · (b·L)，其中 b 为板宽。若 Re_L 超过 5×10⁵，应将平板分为前段层流区和后段湍流区，分别用对应公式再合计。本模拟器只给出局部值，全板积分需要另行计算。
首先是近壁网格：首层中心建议 y+ < 1，δ_99 范围内至少要有 10–20 层网格。其次是入口与前缘几何：Blasius 假设锐边前缘，因此应让前缘位于计算域内部，前面留出一段滑移壁面"助流区"，而不是把入口直接放在前缘处。再次是必须选择层流求解——若启用 k-ε 等湍流模型则结果不会与 Blasius 一致。
不能。Blasius 是 dP/dx = 0 的特殊解。对楔形流动（主流速度 U ∝ x^m）应使用 Falkner–Skan 解；对更复杂的压力分布则用积分形式的边界层方程（如 Thwaites 法）。逆压力梯度（dP/dx > 0）会使层流边界层迅速分离，这往往是过渡或失速的诱因。Blasius 解是上述扩展理论的"零梯度参考"。
实际应用
机翼、船体与涡轮叶片的摩擦阻力估算：大多数气动表面前缘附近会保留一小段层流区，工程上对该段用 Blasius 与平均化公式 C_F = 1.328/√Re_L 计算，再与下游的湍流相关式合计，得到总摩擦阻力。NACA 翼型数据与船舶阻力图表都建立在这种"层流-湍流混合边界层"模型之上。
强制对流换热估算：热边界层与动量边界层形态高度相似（Pr = 1 时几乎重合）。Blasius 解直接导出局部 Nusselt 数 Nu_x ≈ 0.332 Re_x^(1/2) Pr^(1/3)，这是翅片和冷板入口段广为使用的换热相关式。模拟器中 δ_99 越薄，热边界层也越薄，意味着对流换热系数越高。
CFD 软件的基准验证：由于 Blasius 是严格相似解，常被用来检验 CFD 软件的精度。测试新湍流模型或近壁网格策略时，工程师通常先看层流基线在多个 x 处是否能复现 δ_99 和 C_f。本模拟器输出的值可直接作为 CFD 后处理的比较基准。
微流体与微型飞行器：微通道与 MAV 翼运行在很低的雷诺数下，整层边界层都保持层流。此时 Blasius 解几乎可以直接使用：用 δ* 修正通道有效截面，用 C_F 估算微翼阻力。在这些尺度下层流分析远比湍流重要，Blasius 解成了日常使用的工作公式。
常见误解与注意点
最常见的误解是认为边界层厚度"与主流速度无关"。公式 δ_99 = 5x/√Re_x 突出 x，但 Re_x = Ux/ν 中含有 U，因此实际上 δ_99 = 5√(νx/U)，与 √U 成反比。把流速增大 4 倍，边界层厚度就减半。把模拟器中的流速从 1 m/s 拉到 4 m/s，δ_99 会从 6 mm 左右降到 3 mm 左右。
第二个常见错误是把排挤厚度 $\delta^{*}$ 与 99% 厚度 $\delta_{99}$ 混为一谈。Blasius 解中 $\delta^{*}/\delta_{99} \approx 0.344$，$\delta^{*}$ 只有 $\delta_{99}$ 的约 1/3。如果用 $\delta_{99}$ 修正喷嘴的有效流通截面会严重过修正，导致流量与压力损失预测偏差很大。经验规则是：内部流动的有效截面用 $\delta^{*}$，外部流动的"可视厚度"用 $\delta_{99}$，阻力积分用 $\theta$。模拟器同时绘出三者，便于直观比较。
第三个误用是在过渡之后仍套用 Blasius 解。当 Re_x 超过约 5×10⁵ 后边界层变为湍流，厚度按 x^(4/5) 增长，C_f 的量级约为 0.0592/Re_x^(1/5)，远高于层流值。例如长 1 m 的平板在 10 m/s 气流中，板尾的 Re_L ≈ 6.7×10⁵，若把 Blasius 套用到全长会严重低估摩擦。模拟器在 Re_x ≥ 5×10⁵ 时给出红色警告，请务必先检查再把结果用于设计。