湍流边界层模拟器 返回
流体力学模拟器

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

平板上湍流边界层1/7幂律可视化。改变流速、位置、动粘度、密度,实时计算边界层厚度δ、摩擦系数C_f、壁面剪切应力τ_w,学习与Blasius层流解的异同。

参数设置
流速 U
m/s
位置 x
m
运动粘度 ν
m²/s
密度 ρ
kg/m³

湍流公式在 Re_x ≥ 5×10⁵ 时适用;低于该值则按层流区给出警告。

由于已启用减弱动效设置,动画以静止帧显示。

层流区 — 湍流公式不适用
计算结果(实时)
雷诺数 Re_x
湍流边界层厚度 δ_99
局部摩擦系数 C_f
壁面切应力 τ_w
层流 δ_99 (Blasius)
δ湍流 / δ层流
平板上的湍流边界层(实时)

青色箭头=主流 U/灰条=平板/黄色虚线=增长的 δ(x)(湍流 ~x^0.8)/绿色虚线=层流 δ ~√x/旋涡=壁面附近的湍流脉动/圆点=流动示踪。右端为速度分布(蓝=湍流1/7次方律,橙=层流)。红线=转捩点 Re_x=5×10⁵。

速度分布 u(y)/U — 湍流 vs 层流

横轴=无量纲速度 u/U/纵轴=距壁面的无量纲距离 y/δ(蓝=湍流 1/7次方律在壁面附近更"饱满",壁面速度梯度更陡→切应力更高,橙=层流抛物线近似)

理论与主要公式

平板上的湍流边界层由 1/7 次方律速度分布导出的经验关系式描述。与层流(Blasius)解对比,可清晰看出湍流的特性。

雷诺数。U 为主流速度,x 为距前缘的距离,ν 为运动粘度:

$$Re_x = \frac{U\,x}{\nu}, \qquad \text{转捩:}\; Re_x \ge 5\times 10^5$$

湍流边界层厚度、位移厚度、动量厚度(1/7 次方律):

$$\frac{\delta_{99}}{x} = 0.37\,Re_x^{-1/5}, \quad \frac{\delta^{*}}{x} = 0.046\,Re_x^{-1/5}, \quad \frac{\theta}{x} = 0.036\,Re_x^{-1/5}$$

局部摩擦系数与壁面切应力:

$$C_f = 0.059\,Re_x^{-1/5}, \qquad \tau_w = \tfrac{1}{2}\,C_f\,\rho\,U^2$$

与层流(Blasius)解的对比:

$$\frac{\delta_{99}}{x} = 5.0\,Re_x^{-1/2}, \qquad C_f = 0.664\,Re_x^{-1/2}$$

形状因子 H = δ*/θ ≈ 1.28(湍流)、≈ 2.59(层流 Blasius)。在相同 Re_x 下湍流的 Cf 大数倍,这正是湍流摩擦阻力的本质。

湍流边界层模拟器简介

🙋
流体力学课本里提到的"边界层",我只知道平板上有一个薄层。湍流边界层怎么不一样?
🎓
简单说,平板前缘下游的流体被壁面摩擦减速,形成的就是边界层。起初是层流,光滑稳定,但当Re_x超过约5×10⁵时,会转变为湍流。湍流一旦形成,涡流会剧烈搅动运动量,导致边界层比层流厚得多,而且壁面摩擦也大得多。上面的模拟器里,试着增大流速U,会看到Re_x跳升,进入湍流区。
🙋
$\delta/x = 0.37\,Re_x^{-1/5}$ 这个式子是怎么来的?
🎓
这是用"1/7幂律"$u/U = (y/\delta)^{1/7}$这样的经验速度分布,代入动量积分方程推出来的。对比层流的Blasius解$\delta/x = 5.0\,Re_x^{-1/2}$,指数完全不同。湍流的-1/5指数增长得慢,所以Re_x再大,湍流边界层也不会像层流那么薄。模拟器的理论卡里两个式子并排对比。
🙋
摩擦系数Cf湍流比层流大,这为什么是个问题?
🎓
因为摩擦阻力直接吃掉燃油!飞机翼面、船体,大部分阻力来自Cf。标准条件(空气U=20 m/s、x=1 m、ν=1.5×10⁻⁵)下Re_x ≈ 1.33×10⁶,完全湍流,Cf约3.5×10⁻³。同一Re的Blasius层流值(0.664·Re^(-1/2))才5.8×10⁻⁴,只有湍流的1/6!所以飞机设计师费尽心思要"延迟层流转变"、"设计层流翼型"。试试在模拟器里改流速,看Cf的变化。
🙋
壁面剪切应力τ_w有什么用?
🎓
τ_w在平板上积分,就是粘性阻力。CFD仿真算边界层时,也是直接算τ_w或用壁函数逼近。实务中记住"标准空气平板,U=20 m/s时τ_w≈0.84 Pa = 844 mPa"这样的数,能快速判断数值结果是否合理。这就是所谓的"量级核查"。

常见问题

教科书里写的是5×10⁵,但实际会变化很大。表面粗糙、上游湍度、压力分布、振动都会提前转变,有时1×10⁵就转了。反过来,光滑平板在低湍度风洞可能到3×10⁶才转。本模拟器为了实用,设定5×10⁵为参考,低于此值时警告用户。
在Re_x = 5×10⁵~10⁷范围内,与实验的一致性在几%以内。超过10⁷的高Re域,Schlichting对数律式(如Cf = (2·log₁₀Re_x − 0.65)⁻²·³)更准。而Re_x < 5×10⁵就必须用Blasius Cf = 0.664·Re_x^(-1/2),层流公式。
H可以一眼判断流动状态。平板层流(Blasius)约H≈2.59,湍流(1/7幂律)约H≈1.28,接近分离时H会跳到3~4。所以在CFD后处理中,画出H分布就能秒速找出层流-湍流转变、分离前兆等位置,比看原始速度场效率高多了。
本式是"压力梯度=0的平板"专用。逆压力梯度(减速流)会让边界层加厚、易分离,顺压力梯度(加速流)会让边界层变薄。翼型、弯曲体面必须用Thwaites法、Head法或更高级的积分边界层法、CFD。这个模拟器只是"平板基准"。

实际应用

飞机翼的燃效设计:翼面保持层流区域越长,燃耗越低。设计师通过翼型优化、自适应凹槽、吸气等手段延迟转变。本模拟器的Blasius vs 1/7幂律对比,直观展示了"层流有利"的定量差异。

船舶阻力推估:船体总阻力的50~80%是摩擦。航海标准(ITTC 1957)的摩擦线就是以平板乱流Cf实验式为基础,加船型系数。设计流速、吃水长,输入Re,查表or算式得Cf,再乘以表面积、密度、速度平方,估出粘性抗力。

CFD的网格与壁函数:RANS/LES求解边界层,y+分布、壁函数选择至关重要。用本模拟器算出标准平板的δ、τ_w、Cf参考值,与CFD结果对标,能快速检验网格质量和湍流模型合理性。

热交换器与强化对流:对流传热系数h与摩擦Cf有Reynolds-Colburn相似律(St ≈ Cf/2, Pr≈1)。湍流Cf大→h也大,所以热交换器设计中故意用凸起制造湍流来强化传热。

易犯的错误

误区1:"湍流边界层厚,所以阻力小"。完全反了!湍流虽然δ比层流大,但Cf更大,两者相乘后总阻力(τ_w)远比层流高。原因是湍流里涡的剧烈混合导致壁面速度梯度∂u/∂y急峻,剪切应力τ_w = μ·(∂u/∂y)增大。本模拟器对比Cf值一目了然。

误区2:"δ_99就是边界层的物理边界"。流体没有不连续边界。δ_99只是方便定义的"u=0.99U的位置"。按δ_95或δ_90定义值又不一样。学设计用δ*和θ更有物理意义,形状系数H判据也更清晰。

Do not apply the 1/7 power-law correlation at Re_x = 0 or very small Re_x. The model is intended for fully turbulent flow, roughly Re_x >= 5e5. Below transition, use a laminar Blasius estimate; above transition, use this turbulent correlation.

How to use

  1. Set free-stream velocity U in m/s.
  2. Set x as the distance from the plate leading edge in meters.
  3. Set kinematic viscosity nu and density rho for the fluid. For 20 C air, use nu about 1.5e-5 m2/s and rho about 1.2 kg/m3.
  4. Check Re_x, delta99, Cf, and wall shear stress tau_w. Treat Re_x below 5e5 as a laminar/transitional warning.

Measured example

For air with U=10 m/s, x=5 m, nu=1.5e-5 m2/s, and rho=1.2 kg/m3, Re_x=3.333e6. The 1/7 power law gives delta99=0.0917 m (91.7 mm), Cf=2.93e-3, and tau_w=0.1756 Pa, so the turbulent formula is in range.

Notes

  1. Below Re_x=5e5, use Blasius or a transition model rather than the turbulent 1/7 relation as a design value.
  2. For transition location, combine this estimate with real methods such as the Michel criterion or e^N method.
  3. tau_w displays in Pa above 1 Pa and in mPa below 1 Pa.