湍流边界层模拟器 — 平板 1/7 次幂律与 Cf

飞机机翼与摩擦阻力预测：翼面边界层能否保持层流、何时转捩到湍流，是与燃油消耗直接相关的最重要设计项。设计师通过层流翼（NLF 翼）形状优化或人为安装转捩条来控制阻力。1/7 次幂律与 Blasius 解的对比，是定量评估「保持层流的好处」的第一步。

船舶的船体摩擦阻力：船舶总阻力的 50–80% 来自船体表面的摩擦阻力，其中大部分由湍流边界层引起。船级社标准摩擦线（如 ITTC 1957）本质上就是平板湍流摩擦系数 Cf 的实验式加形状修正。Froude 法把模型试验外推到实船性能的起点，正是本工具的公式。

CFD 壁函数与验证：在 RANS 或 LES 中，壁函数选择与计算网格 y+ 分布显著影响近壁预测的精度。用本工具计算教科书的平板 Cf、δ、τ_w 并与 CFD 结果对比，是验证网格质量与湍流模型设置的快速「sanity check」，设计现场不可或缺。

对流传热系数预测：Reynolds-Colburn 类似律将对流传热系数 h 与摩擦系数 Cf 紧密关联（St ≈ Cf/2，Pr ≈ 1）。湍流边界层 Cf 较大，对流传热也比层流高出数倍。换热器与冷却翅片设计中经常故意安装促湍流肋条，正是基于此。

最常见的误解是凭直觉认为「湍流边界层更厚，所以阻力应该更小」。事实恰恰相反：湍流边界层不仅更厚，壁面摩擦也更大。原因是湍流涡剧烈混合动量，将外部高速流体输送到壁面附近，使壁面速度梯度 $\partial u/\partial y$ 比层流（抛物线分布）更陡，于是剪应力 $\tau_w = \mu(\partial u/\partial y)$ 增大。在本工具中将 C_f 与 Blasius 解比较即可一目了然：相同 Re_x 下湍流值大数倍。「更厚但阻力更大」是湍流边界层的本质。

第二个常见误解是把边界层厚度 $\delta_{99}$ 当作「物理上明确的边界」。实际流体中并不存在不连续的边界，$\delta_{99}$ 只是「u 达到来流速度 U 的 99% 处」这一便利位置。根据定义不同也可使用 $\delta_{95}$ 或 $\delta_{90}$，数值会不同。设计中更可靠的是物理意义更明确的位移厚度 $\delta^*$ 和动量厚度 $\theta$，形状因子 H = $\delta^*/\theta$ 在判别转捩与分离时也比 $\delta_{99}$ 本身更有效。

最后注意，不要在 Re_x = 0 或非常小的 Re_x 处套用 1/7 次幂律公式。本式针对完全湍流域 (Re_x ≥ 5×10⁵) 推导，套用到层流域会得到物理上不可能的数值。例如 Re_x = 10⁴ 代入湍流式得 δ/x = 0.37·(10⁴)^(-0.2) ≈ 0.059，即边界层厚度为位置的 6%，明显不合理。本工具在 Re_x < 5×10⁵ 时显式显示警告，请务必确认。实际工程中，应在层流域使用 Blasius，在湍流域使用本式（或对数律式）。