引言
流体力学的核心问题是求解 Navier–Stokes 方程，而这组非线性偏微分方程的精确解屈指可数。1904 年 Prandtl 提出的边界层理论，通过”黏性作用只集中在壁面附近薄层内”的物理洞察，把大雷诺数流动分解为外部无黏势流与内部黏性薄层两个可分别求解的问题，成为近代流体力学与空气动力学的奠基性思想。本文从控制方程出发，推导边界层近似与平板层流的 Blasius 相似解，用数值方法复现这一经典结果，并以平板摩阻估算为例说明其工程用法。
Navier–Stokes 方程
对流体微团应用质量守恒与动量守恒，并采用牛顿黏性本构（应力与应变率线性相关），得到不可压缩常黏度流动的控制方程。连续性方程为：
$$\nabla\cdot\boldsymbol{u} = \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} = 0$$
动量方程（Navier–Stokes 方程）为：
$$\rho\left(\frac{\partial \boldsymbol{u}}{\partial t} + (\boldsymbol{u}\cdot\nabla)\boldsymbol{u}\right) = -\nabla p + \mu\nabla^2\boldsymbol{u} + \rho\boldsymbol{f}$$
其中 $\rho$ 为密度，$p$ 为压强，$\mu$ 为动力黏度，$\boldsymbol{f}$ 为单位质量体力。将方程无量纲化后，惯性项与黏性项的比值归结为唯一的无量纲参数——雷诺数 $\mathrm{\Re} = \rho U L/\mu$，它控制流动的相似性：几何相似且雷诺数相等的流动动力学相似。当 $\mathrm{\Re}$ 很大时黏性项在流场主体中可忽略，但若直接略去黏性，方程降阶导致壁面无滑移条件无法满足——这一矛盾正是边界层理论的出发点。
边界层理论
Prandtl 的关键观察是：大雷诺数下，速度从壁面零值过渡到来流值的区域厚度 $\delta$ 远小于特征长度 $L$。在边界层内对二维定常方程做量级分析：法向动量方程退化为 $\partial p/\partial y \approx 0$，即压强沿法向不变，由外部势流”印加”到壁面；流向动量方程中只保留法向黏性项，得到边界层方程：
$$u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = -\frac{1}{\rho}\frac{\mathrm{d}p}{\mathrm{d}x} + \nu\frac{\partial^2 u}{\partial y^2}$$
其中 $\nu = \mu/\rho$ 为运动黏度。方程虽被简化为抛物型（可沿流向步进求解），但对零压强梯度平板仍存在精确相似解。Blasius 引入相似变量 $\eta = y\sqrt{U_\infty/(\nu x)}$ 与流函数 $f(\eta)$，使 $u/U_\infty = f'(\eta)$，偏微分方程化为常微分方程边值问题：
$$f''' + \frac{1}{2}f,f'' = 0, \qquad f(0) = f'(0) = 0,\quad f'(\infty) = 1$$
该问题无初等函数解，数值结果为 $f”(0) = 0.332057$，由此导出全部工程关心的特征量，汇总如下表（$\mathrm{\Re}_x = U_\infty x/\nu$）。
| 特征量 | 定义 | Blasius 结果 |
|---|---|---|
| 名义厚度 $\delta_{99}$ | $u = 0.99,U_\infty$ 处 | $5.0,x/\sqrt{\mathrm{\Re}_x}$ |
| 位移厚度 $\delta^*$ | 流量亏损等效 | $1.721,x/\sqrt{\mathrm{\Re}_x}$ |
| 动量厚度 $\theta$ | 动量亏损等效 | $0.664,x/\sqrt{\mathrm{\Re}_x}$ |
| 局部摩阻系数 $c_f$ | $2\tau_w/(\rho U_\infty^2)$ | $0.664/\sqrt{\mathrm{\Re}_x}$ |
| 形状因子 $H$ | $\delta^*/\theta$ | 2.591 |
对更一般的有压强梯度流动，将边界层方程沿法向积分并利用连续性方程消去法向速度，可得 von Kármán 动量积分关系式：
$$\frac{\mathrm{d}\theta}{\mathrm{d}x} + (2+H),\frac{\theta}{U}\frac{\mathrm{d}U}{\mathrm{d}x} = \frac{c_f}{2}$$
该方程把壁面摩阻与位移厚度、动量厚度两个积分量直接联系起来，是边界层近似解法（如 Thwaites 法）的基础：给定外部势流速度分布 $U(x)$，即可沿流向积分求得 $\theta(x)$ 与摩阻分布，而无需求解偏微分方程，并可根据形状因子的增长预估分离位置（层流分离大致对应 $H \approx 3.5$ 到 4.0）。
边界层理论的物理价值不止于平板解：它解释了逆压梯度（$\mathrm{d}p/\mathrm{d}x > 0$）如何使近壁低动量流体减速直至壁面剪应力归零——流动分离，进而导致压差阻力骤增；这也解释了流线型设计为何能大幅减阻。工程上判断边界层状态的经验转捩准则为 $\mathrm{\Re}_x \approx 5\times10^5$，超过该值层流解不再适用。
数值方法与实现
Blasius 方程是三阶非线性两点边值问题，可用打靶法（把 $f”(0)$ 作为待定参数迭代）或直接用边值问题求解器求解。下面用 scipy.integrate.solve_bvp 在截断区间 $\eta\in[0,8]$ 上求解，并输出与文献值的对照：
import numpy as np
from scipy.integrate import solve_bvp
# Blasius 方程 f''' + f f''/2 = 0 化为一阶方程组，y = [f, f', f'']
def ode(eta, y):
return np.vstack((y[1], y[2], -0.5 * y[0] * y[2]))
def bc(ya, yb):
return np.array([ya[0], ya[1], yb[1] - 1.0]) # f(0)=f'(0)=0, f'(∞)=1
eta = np.linspace(0.0, 8.0, 161)
guess = np.zeros((3, eta.size))
guess[1] = np.clip(eta / 3.0, 0.0, 1.0) # f' 的初始猜测
sol = solve_bvp(ode, bc, eta, guess)
f, fp, fpp = sol.sol(eta)
print(f"f''(0) = {fpp[0]:.6f} (文献值 0.332057)")
i99 = np.argmax(fp >= 0.99)
print(f"η_99 = {eta[i99]:.3f} (对应名义厚度系数)")
print(" η f f' f''")
for e in [0, 1, 2, 3, 4, 5, 6]:
fv = sol.sol(float(e))
print(f"{e:4.1f} {fv[0]:7.4f} {fv[1]:7.4f} {fv[2]:7.4f}")程序输出 $f”(0) = 0.332059$，与文献值 $0.332057$ 的偏差约 $2\times10^{-6}$，源于无穷远边界截断在 $\eta = 8$；$\eta_{99} \approx 4.95$，对应工程中常用的名义厚度系数 5.0；$f'(2) \approx 0.63$、$f'(3) \approx 0.85$ 等数值与经典 Blasius 表逐项吻合，验证了数值实现的正确性。
工程案例分析：平板边界层与摩阻估算
某风洞试验段内气流以 $U_\infty = 10$ m/s 流过一块长 $L = 0.3$ m 的平板，常温空气取 $\rho = 1.225$ kg/m³、$\nu = 1.5\times10^{-5}$ m²/s。板尾雷诺数为：
$$\mathrm{\Re}_L = \frac{U_\infty L}{\nu} = \frac{10\times0.3}{1.5\times10^{-5}} = 2.0\times10^{5}$$
低于转捩临界值 $5\times10^5$，全板可按层流处理。代入上表公式：板尾名义厚度 $\delta_{99} = 5.0\times0.3/\sqrt{2.0\times10^5} \approx 3.35$ mm——不足板长的 1.2%，印证了 $\delta \ll L$ 的前提；板尾局部摩阻系数 $c_f \approx 1.49\times10^{-3}$，壁面剪应力约 0.091 Pa；整板平均摩阻系数 $C_f = 1.328/\sqrt{\mathrm{\Re}_L} \approx 2.97\times10^{-3}$，单面单位展宽的摩擦阻力约 0.055 N。这组计算展示了工程估算的标准流程：先由雷诺数判断流态，再选择层流或湍流关联式；若同一平板速度提高到 40 m/s，$\mathrm{\Re}_L$ 将达 $8\times10^5$，必须改用湍流边界层公式（如 $c_f = 0.0592,\mathrm{\Re}_x^{-1/5}$），否则摩阻将被低估近一半。
结论与展望
Navier–Stokes 方程完整描述了黏性流动，但其非线性使解析求解极为困难；边界层理论利用大雷诺数下的尺度分离，把问题分解为外部势流与内部薄层，并以 Blasius 相似解给出了层流摩阻的基准结果。该理论同时揭示了逆压梯度导致分离的机理，是理解翼型失速、扩压器喘振等工程问题的钥匙。进一步的方向包括湍流边界层的半经验理论、可压缩与热边界层，以及与 Navier–Stokes 直接数值模拟相结合的现代壁面建模方法。
发表回复