坐标、法向和无量纲量
考虑固壁位于 y = 0 y=0 y = 0 y > 0 y>0 y > 0 x x x n ^ = + y ^ \hat{\mathbf n}=+\hat{\mathbf y} n ^ = + y ^ u = ( u , v ) \mathbf u=(u,v) u = ( u , v ) v > 0 v>0 v > 0
u ( x , 0 ) = 0 , v ( x , 0 ) = 0. u(x,0)=0,\qquad v(x,0)=0. u ( x , 0 ) = 0 , v ( x , 0 ) = 0. 边界层外缘 y → ∞ y\to\infty y → ∞ u → U e ( x ) u\to U_e(x) u → U e ( x ) − y ^ -\hat{\mathbf y} − y ^
不可压 Newton 流体满足
∇ ⋅ u = 0 , \nabla\cdot\mathbf u=0, ∇ ⋅ u = 0 , ρ ( ∂ t u + u ⋅ ∇ u ) = − ∇ p + μ ∇ 2 u . \rho(\partial_t\mathbf u+\mathbf u\cdot\nabla\mathbf u)
=-\nabla p+\mu\nabla^2\mathbf u. ρ ( ∂ t u + u ⋅ ∇ u ) = − ∇ p + μ ∇ 2 u . ρ \rho ρ k g m − 3 \mathrm{kg\,m^{-3}} kg m − 3 μ \mu μ P a s \mathrm{Pa\,s} Pa s ν = μ / ρ \nu=\mu/\rho ν = μ / ρ m 2 s − 1 \mathrm{m^2\,s^{-1}} m 2 s − 1 U U U L L L L / U L/U L / U ρ U 2 \rho U^2 ρ U 2
∂ t ∗ u ∗ + u ∗ ⋅ ∇ ∗ u ∗ = − ∇ ∗ p ∗ + 1 R e ∇ ∗ 2 u ∗ , \partial_{t^*}\mathbf u^*+\mathbf u^*\cdot\nabla^*\mathbf u^*
=-\nabla^*p^*+\frac1{Re}\nabla^{*2}\mathbf u^*, ∂ t ∗ u ∗ + u ∗ ⋅ ∇ ∗ u ∗ = − ∇ ∗ p ∗ + R e 1 ∇ ∗ 2 u ∗ , R e = U L ν = ρ U L μ . Re=\frac{UL}{\nu}=\frac{\rho UL}{\mu}. R e = ν UL = μ ρ UL . R e Re R e R e Re R e
薄层尺度与方程
令沿壁变化尺度 L L L δ ≪ L \delta\ll L δ ≪ L v / U ∼ δ / L v/U\sim\delta/L v / U ∼ δ / L u ∂ x u ∼ U 2 / L u\partial_xu\sim U^2/L u ∂ x u ∼ U 2 / L ν ∂ y y u ∼ ν U / δ 2 \nu\partial_{yy}u\sim\nu U/\delta^2 ν ∂ yy u ∼ νU / δ 2
δ L ∼ R e − 1 / 2 . \frac\delta L\sim Re^{-1/2}. L δ ∼ R e − 1/2 . 保留主导项，二维稳态边界层方程为
∂ x u + ∂ y v = 0 , \partial_xu+\partial_yv=0, ∂ x u + ∂ y v = 0 , u ∂ x u + v ∂ y u = − 1 ρ d p e d x + ν ∂ y y u . u\partial_xu+v\partial_yu
=-\frac1\rho\frac{\mathrm dp_e}{\mathrm dx}
+\nu\partial_{yy}u. u ∂ x u + v ∂ y u = − ρ 1 d x d p e + ν ∂ yy u . 法向动量主阶给 ∂ y p = 0 \partial_yp=0 ∂ y p = 0 p e ( x ) p_e(x) p e ( x ) U e d U e / d x = − ( 1 / ρ ) d p e / d x U_e\mathrm dU_e/\mathrm dx=-(1/\rho)\mathrm dp_e/\mathrm dx U e d U e / d x = − ( 1/ ρ ) d p e / d x δ \delta δ
这个尺度结论也可由双尺度无量纲化逐项核验。令
x = L x ∗ , y = δ y ∗ , u = U u ∗ , v = U δ L v ∗ . x=Lx^*,\qquad y=\delta y^*,\qquad
u=Uu^*,\qquad v=U\frac{\delta}{L}v^*. x = L x ∗ , y = δ y ∗ , u = U u ∗ , v = U L δ v ∗ . 连续性两项均为 U / L U/L U / L U 2 / L U^2/L U 2 / L x x x R e − 1 Re^{-1} R e − 1 y y y R e − 1 ( L / δ ) 2 Re^{-1}(L/\delta)^2 R e − 1 ( L / δ ) 2 R e − 1 ( L / δ ) 2 = O ( 1 ) Re^{-1}(L/\delta)^2=O(1) R e − 1 ( L / δ ) 2 = O ( 1 ) δ / L = O ( R e − 1 / 2 ) \delta/L=O(Re^{-1/2}) δ / L = O ( R e − 1/2 ) O ( R e − 1 ) O(Re^{-1}) O ( R e − 1 ) ∂ x x u \partial_{xx}u ∂ xx u x → 0 x\to0 x → 0
局部速度剖面还可压缩为三个可测厚度：
δ ∗ = ∫ 0 ∞ ( 1 − u U e ) d y , θ = ∫ 0 ∞ u U e ( 1 − u U e ) d y , H = δ ∗ θ . \delta^*=\int_0^\infty\left(1-\frac{u}{U_e}\right)\mathrm dy,
\qquad
\theta=\int_0^\infty\frac{u}{U_e}\left(1-\frac{u}{U_e}\right)\mathrm dy,
\qquad H=\frac{\delta^*}{\theta}. δ ∗ = ∫ 0 ∞ ( 1 − U e u ) d y , θ = ∫ 0 ∞ U e u ( 1 − U e u ) d y , H = θ δ ∗ . 位移厚度 δ ∗ \delta^* δ ∗ θ \theta θ H H H
τ w ρ U e 2 = d θ d x + θ U e ( 2 + H ) d U e d x . \frac{\tau_w}{\rho U_e^2}
=\frac{\mathrm d\theta}{\mathrm dx}
+\frac{\theta}{U_e}(2+H)\frac{\mathrm dU_e}{\mathrm dx}. ρ U e 2 τ w = d x d θ + U e θ ( 2 + H ) d x d U e . 等式每项都无量纲；它把剪应力、剖面增长和外流加速联系起来，但只给一条关系，仍需相似剖面或经验闭合来确定 H H H τ w \tau_w τ w u ( x , y ) u(x,y) u ( x , y )
工程数据中的 δ 99 \delta_{99} δ 99 u / U e = 0.99 u/U_e=0.99 u / U e = 0.99 δ ∗ \delta^* δ ∗ θ \theta θ
例 1：层厚的尺度估算
水样流体取 ν = 1.0 × 10 − 6 m 2 s − 1 \nu=1.0\times10^{-6}\,\mathrm{m^2\,s^{-1}} ν = 1.0 × 1 0 − 6 m 2 s − 1 U = 2.0 m s − 1 U=2.0\,\mathrm{m\,s^{-1}} U = 2.0 m s − 1 x = 0.50 m x=0.50\,\mathrm m x = 0.50 m
R e x = U x ν = 1.0 × 10 6 . Re_x=\frac{Ux}{\nu}=1.0\times10^6. R e x = ν Ux = 1.0 × 1 0 6 . 尺度估计 δ / x ∼ R e x − 1 / 2 = 10 − 3 \delta/x\sim Re_x^{-1/2}=10^{-3} δ / x ∼ R e x − 1/2 = 1 0 − 3 δ ∼ 0.50 m m \delta\sim0.50\,\mathrm{mm} δ ∼ 0.50 mm 99 % 99\% 99% δ 99 ≈ 2.5 m m \delta_{99}\approx2.5\,\mathrm{mm} δ 99 ≈ 2.5 mm
Blasius 相似解与摩擦
零压强梯度、均匀 U e = U ∞ U_e=U_\infty U e = U ∞
η = y U ∞ ν x , ψ = ν U ∞ x f ( η ) , \eta=y\sqrt{\frac{U_\infty}{\nu x}},
\qquad
\psi=\sqrt{\nu U_\infty x}\,f(\eta), η = y νx U ∞ , ψ = ν U ∞ x f ( η ) , u = ∂ y ψ u=\partial_y\psi u = ∂ y ψ v = − ∂ x ψ v=-\partial_x\psi v = − ∂ x ψ
f ′ ′ ′ + 1 2 f f ′ ′ = 0 , f'''+\frac12ff''=0, f ′′′ + 2 1 f f ′′ = 0 , 边界 f ( 0 ) = f ′ ( 0 ) = 0 f(0)=f'(0)=0 f ( 0 ) = f ′ ( 0 ) = 0 f ′ ( ∞ ) = 1 f'(\infty)=1 f ′ ( ∞ ) = 1
C f , x = τ w 1 2 ρ U ∞ 2 ≈ 0.664 R e x , C_{f,x}=\frac{\tau_w}{\tfrac12\rho U_\infty^2}
\approx\frac{0.664}{\sqrt{Re_x}}, C f , x = 2 1 ρ U ∞ 2 τ w ≈ R e x 0.664 , τ w = μ ( ∂ y u ) w \tau_w=\mu(\partial_yu)_w τ w = μ ( ∂ y u ) w x x x τ w \tau_w τ w
Blasius 剖面的动量厚度为
θ = 0.664 ν x U ∞ . \theta=0.664\sqrt{\frac{\nu x}{U_\infty}}. θ = 0.664 U ∞ νx . 零压强梯度时积分动量式化为 τ w / ( ρ U ∞ 2 ) = d θ / d x \tau_w/(\rho U_\infty^2)=\mathrm d\theta/\mathrm dx τ w / ( ρ U ∞ 2 ) = d θ / d x d θ / d x = 0.332 / R e x \mathrm d\theta/\mathrm dx=0.332/\sqrt{Re_x} d θ / d x = 0.332/ R e x C f , x = 0.664 / R e x C_{f,x}=0.664/\sqrt{Re_x} C f , x = 0.664/ R e x x 1 / 2 x^{1/2} x 1/2 x − 1 / 2 x^{-1/2} x − 1/2
例 2：平板局部剪应力
沿用 ρ = 1000 k g m − 3 \rho=1000\,\mathrm{kg\,m^{-3}} ρ = 1000 kg m − 3 U = 2.0 m s − 1 U=2.0\,\mathrm{m\,s^{-1}} U = 2.0 m s − 1 R e x = 10 6 Re_x=10^6 R e x = 1 0 6
C f , x = 6.64 × 10 − 4 , C_{f,x}=6.64\times10^{-4}, C f , x = 6.64 × 1 0 − 4 , ∣ τ w ∣ = 1 2 ρ U 2 C f , x = 1.33 P a . |\tau_w|=\frac12\rho U^2C_{f,x}=1.33\,\mathrm{Pa}. ∣ τ w ∣ = 2 1 ρ U 2 C f , x = 1.33 Pa . 这是局部值，不是从前缘到 x x x
Blasius 解之所以存在，是因为平板问题在改变 x x x U e ( x ) = K x m U_e(x)=Kx^m U e ( x ) = K x m
f ′ ′ ′ + f f ′ ′ + β ( 1 − f ′ 2 ) = 0 , β = 2 m m + 1 , f'''+ff''+\beta\bigl(1-f'^2\bigr)=0,
\qquad \beta=\frac{2m}{m+1}, f ′′′ + f f ′′ + β ( 1 − f ′2 ) = 0 , β = m + 1 2 m , 并配合 f ( 0 ) = f ′ ( 0 ) = 0 f(0)=f'(0)=0 f ( 0 ) = f ′ ( 0 ) = 0 f ′ ( ∞ ) = 1 f'(\infty)=1 f ′ ( ∞ ) = 1 m > 0 m>0 m > 0 m < 0 m<0 m < 0 U e ( x ) U_e(x) U e ( x ) x x x
若壁面有均匀吸气或吹气，边界条件改为 v ( x , 0 ) = v w v(x,0)=v_w v ( x , 0 ) = v w + y +y + y v w < 0 v_w<0 v w < 0 v w > 0 v_w>0 v w > 0 S = v w x / ( ν U e ) S=v_w\sqrt{x/(\nu U_e)} S = v w x / ( ν U e ) S S S v w v_w v w U e U_e U e
逆压梯度与分离
d p e / d x > 0 \mathrm dp_e/\mathrm dx>0 d p e / d x > 0 U e U_e U e
τ w = μ ( ∂ y u ) w = 0 \tau_w=\mu(\partial_yu)_w=0 τ w = μ ( ∂ y u ) w = 0 为局部判据，其后可能出现近壁回流。边界层在分离邻域变厚，薄层和定常假设可能失效；零剪应力位置只是渐近理论的预警，不能自动给完整涡脱落区。
例 3：由外流判断压强梯度
若外缘速度 U e ( x ) = U 0 ( 1 − x / L ) U_e(x)=U_0(1-x/L) U e ( x ) = U 0 ( 1 − x / L ) 0 < x < L 0<x<L 0 < x < L
d p e d x = − ρ U e d U e d x = ρ U 0 U e L > 0. \frac{\mathrm dp_e}{\mathrm dx}
=-\rho U_e\frac{\mathrm dU_e}{\mathrm dx}
=\frac{\rho U_0U_e}{L}>0. d x d p e = − ρ U e d x d U e = L ρ U 0 U e > 0. 这是逆压梯度，具有分离倾向。仅由 U e U_e U e τ w \tau_w τ w x → L x\to L x → L
壁面附近不一定由沿流长度控制。若无限平板以角频率 Ω \Omega Ω ∂ t u ∼ Ω U \partial_tu\sim\Omega U ∂ t u ∼ Ω U ν U / δ s 2 \nu U/\delta_s^2 νU / δ s 2
δ s = 2 ν Ω . \delta_s=\sqrt{\frac{2\nu}{\Omega}}. δ s = Ω 2 ν . Ω \Omega Ω s − 1 \mathrm{s^{-1}} s − 1 δ s \delta_s δ s 2 \sqrt2 2 e − y / δ s e^{-y/\delta_s} e − y / δ s δ s \delta_s δ s R e x − 1 / 2 Re_x^{-1/2} R e x − 1/2
定常二维流中，τ w = 0 \tau_w=0 τ w = 0
线性稳定性：先写基流和扰动类
取近似平行基流 U = ( U ( y ) , 0 ) \mathbf U=(U(y),0) U = ( U ( y ) , 0 ) u ′ = ϵ u ~ \mathbf u'=\epsilon\tilde{\mathbf u} u ′ = ϵ u ~
ψ ~ = ϕ ( y ) e i ( α x − ω t ) . \tilde\psi=\phi(y)e^{i(\alpha x-\omega t)}. ψ ~ = ϕ ( y ) e i ( αx − ω t ) . α \alpha α m − 1 \mathrm{m^{-1}} m − 1 ω \omega ω s − 1 \mathrm{s^{-1}} s − 1 α \alpha α Im ω > 0 \operatorname{Im}\omega>0 Im ω > 0 R e Re R e
线性稳定只判断无穷小扰动的初始增长。算符非正规时，即使所有本征模衰减，多个非正交模仍可产生有限时间瞬态增长；有限振幅扰动还可触发非线性转捩。实验“临界 Reynolds 数”依赖来流噪声、粗糙度、压力梯度和判据，不是只由量纲分析给出的普适常数。
把速度、长度分别用 U 0 U_0 U 0 h h h D = d / d y D=\mathrm d/\mathrm dy D = d / d y c = ω / α c=\omega/\alpha c = ω / α
( U − c ) ( D 2 − α 2 ) ϕ − U ′ ′ ϕ = 1 i α R e ( D 2 − α 2 ) 2 ϕ . (U-c)(D^2-\alpha^2)\phi-U''\phi
=\frac{1}{i\alpha Re}(D^2-\alpha^2)^2\phi. ( U − c ) ( D 2 − α 2 ) ϕ − U ′′ ϕ = i α R e 1 ( D 2 − α 2 ) 2 ϕ . 此式中的 y y y α \alpha α U U U c c c R e = U 0 h / ν Re=U_0h/\nu R e = U 0 h / ν ϕ = D ϕ = 0 \phi=D\phi=0 ϕ = D ϕ = 0 ϕ , D ϕ → 0 \phi,D\phi\to0 ϕ , D ϕ → 0 y max y_{\max} y m a x
扰动能量给出比单个本征值更直接的机制判断。对无边界能量通量的无量纲平行流，二维扰动满足
d d t ∫ V ∣ u ′ ∣ 2 2 d V = − ∫ V u ′ v ′ U ′ d V − 1 R e ∫ V ∣ ∇ u ′ ∣ 2 d V . \frac{\mathrm d}{\mathrm dt}\int_V\frac{|\mathbf u'|^2}{2}\,\mathrm dV
=-\int_V u'v' U'\,\mathrm dV
-\frac1{Re}\int_V|\nabla\mathbf u'|^2\,\mathrm dV. d t d ∫ V 2 ∣ u ′ ∣ 2 d V = − ∫ V u ′ v ′ U ′ d V − R e 1 ∫ V ∣∇ u ′ ∣ 2 d V . 第二项恒为非正黏性耗散；第一项可正可负，表示 Reynolds 应力从基流剪切提取或返还能量。若采用有入口、出口的空间增长问题，还要保留边界能量通量。瞬态增长应由演化算子的奇异值或给定初值的能量增益
G ( t ) = E ( t ) / E ( 0 ) G(t)=E(t)/E(0) G ( t ) = E ( t ) / E ( 0 )
局部时间分析假设基流沿 x x x
例 4：振幅增长和能量增长
某时间模的 Im ω = 20 s − 1 \operatorname{Im}\omega=20\,\mathrm{s^{-1}} Im ω = 20 s − 1 A 0 = 10 − 4 U A_0=10^{-4}U A 0 = 1 0 − 4 U 0.10 s 0.10\,\mathrm s 0.10 s
A A 0 = e ( 20 ) ( 0.10 ) = e 2 = 7.39. \frac{A}{A_0}=e^{(20)(0.10)}=e^2=7.39. A 0 A = e ( 20 ) ( 0.10 ) = e 2 = 7.39. 扰动能量与振幅平方成正比，增长因子 e 4 = 54.6 e^4=54.6 e 4 = 54.6
其他不稳定与模型边界
Kelvin–Helmholtz 型不稳定由剪切层速度差供能；Rayleigh–Bénard 由不稳定密度分层和浮力驱动；离心不稳定与曲线流线有关。它们的控制量除 R e Re R e
可压高速流需 Mach 数和能量方程；热边界层涉及 Prandtl 数；非 Newton 流体需本构时间尺度；粗糙壁面引入粗糙度比。边界层近似离开壁面不替代外部 Euler 解，两者必须匹配。
比较不同转捩实验时，最低限度应报告外缘速度和温度、局部 R e x Re_x R e x R e θ = U e θ / ν Re_\theta=U_e\theta/\nu R e θ = U e θ / ν U e U_e U e
数值计算也要区分物理不稳定和离散不稳定。前者在网格、时间步和计算域收敛后保留一致的增长率；后者常随离散参数漂移，或由入口反射、出口回流、人工耗散不足触发。可复核流程是：先验证基流残差和守恒，再施加幅值已知的单频扰动，在仍属线性区间内拟合 log A ( t ) \log A(t) log A ( t ) log A ( x ) \log A(x) log A ( x )
这些方法的适用层级应保持清楚：尺度分析判断主导项，积分方程估算厚度和摩擦，相似解给特殊外流的完整剖面，线性理论诊断无穷小扰动，直接数值模拟才可能解析给定初边值问题的非线性演化。后一层的计算量更大，却不自动修复错误的物性、几何或入口条件。
练习
练习 1：Re 量纲
所属知识 无量纲化
难度 2/5 查看提示 用
ν \nu ν 的
m 2 / s m^{2}/s m 2 / s 单位。
查看解答 U L / ν = ( m / s ⋅ m ) / ( m 2 / s ) = 1 UL/\nu=(m/s\cdot m)/(m^{2}/s)=1 UL / ν = ( m / s ⋅ m ) / ( m 2 / s ) = 1 ；也等于
ρ U L / μ \rho UL/\mu ρ UL / μ 。
练习 2：法向速度尺度
所属知识 边界层
难度 3/5 查看提示 用
∂ x u + ∂ y v = 0 \partial_x u+\partial_y v=0 ∂ x u + ∂ y v = 0 。
查看解答 U / L ∼ V n / δ U/L\sim V_n/\delta U / L ∼ V n / δ ，故
V n / U ∼ δ / L ∼ R e − 1 / 2 V_n/U\sim \delta/L\sim Re^{-1/2} V n / U ∼ δ / L ∼ R e − 1/2 。
练习 3：局部与平均摩擦
所属知识 Blasius
难度 3/5 查看提示 局部
C f ∝ x − 1 / 2 C_f\propto x^{-1/2} C f ∝ x − 1/2 ，对剪力沿 x 积分。
查看解答 平均系数为
1.328 / R e L 1.328/\sqrt{\mathrm{Re}_L} 1.328/ Re L ，是末端局部系数的两倍；只适用于层流零压梯度平板。
练习 4：分离符号
所属知识 压强梯度
难度 3/5 查看提示 外流 Euler 关系给
d p / d x = − ρ U e d U e / d x dp/dx=-\rho U_e dU_e/dx d p / d x = − ρ U e d U e / d x 。
查看解答 U e U_e U e 减小给
d p / d x > 0 dp/dx>0 d p / d x > 0 ，逆压梯度使
τ w \tau_w τ w 降低；
τ w = 0 \tau_w=0 τ w = 0 是经典分离判据。
练习 5：时间稳定判据
所属知识 模态
难度 3/5 查看提示 展开
e i ( α x − ω t ) e^{i(\alpha x-\omega t)} e i ( αx − ω t ) 中复
ω \omega ω 。
查看解答 ω = ω r + i ω i \omega=\omega_r+i\omega_i ω = ω r + i ω i 时振幅为
e ω i t e^{\omega_i t} e ω i t ；
ω i > 0 \omega_i>0 ω i > 0 增长。空间分析不可沿用同一符号。
练习 6：转捩证据
所属知识 稳定性边界
难度 4/5 查看提示 区分本征模、瞬态、有限振幅和非线性饱和。
查看解答 谱稳定不排除非正规瞬态增长；瞬态增长也不保证持续湍流。需报告扰动幅度、噪声、时间窗和非线性结果。
关系与资源
课程 · 2013 Advanced Fluid Mechanics Gareth McKinley
用于核对 P09 的守恒方程符号、应力量纲、无黏与黏性流推导、边界层尺度和流动稳定性例题。
打开官方来源 MIT OpenCourseWare 2.25 覆盖边界层、相似性和稳定性。本文明确了 R e Re R e
课程 · 2015 Numerical Fluid Mechanics Pierre Lermusiaux
用于核对 P11 的离散化误差、场方程网格算法、时间步稳定性、边界条件和数值验证流程。
打开官方来源 MIT OpenCourseWare 2.29《Numerical Fluid Mechanics》支撑本章对薄边界层和不稳定模态的数值分辨要求：壁面法向网格必须解析速度梯度，计算域与入口扰动会影响增长率，时间步和空间离散还需通过加密比较区分物理放大与数值耗散。该来源具体支持把壁面剪切、分离位置和模态增长率纳入验证清单，而不允许用单一彩色流线图确认转捩机制。