9992 字
50 分钟
FluidMechanics-Chapter6:流体阻力及能量损失

概述#

这一章研究一个核心问题:

真实流体为什么会损失机械能,以及如何计算这些损失。

整章的逻辑链条是:

  1. 黏性与湍动产生阻力;
  2. 阻力做功,使机械能不可逆地转化为内能,表现为水头损失;
  3. 先判断流态,再确定沿程阻力系数;
  4. 将沿程损失与各处局部损失相加;
  5. 对物体绕流,还要分析边界层、分离、尾流和绕流阻力。

本章常见英文术语:

  • flow resistance:流动阻力
  • energy loss / head loss:能量损失 / 水头损失
  • frictional head loss:沿程水头损失
  • minor/local head loss:局部水头损失
  • laminar flow:层流
  • turbulent flow:湍流
  • viscous sublayer:黏性底层
  • boundary layer:边界层
  • drag force:绕流阻力
TIP

教学范围判断

四份课件依次覆盖本章 Section 1–9:阻力概述、两种流态、圆管层流、湍流、圆管湍流沿程损失、非圆管损失、局部损失、边界层、绕流阻力。课件中没有发现明确的“not required”或“not examined”标注,因此本笔记全部纳入。

课件中的外部论文链接属于 extended reading;雷诺应力和混合长度的完整推导考试优先级较低,但其定义、物理意义和最终公式仍应掌握。

WARNING

本章最容易混淆的三组量

  1. uu:某点速度;uˉ\bar u:时均速度;vv:断面平均速度;
  2. τ0\tau_0:壁面切应力;uu_*:摩阻流速,它有速度量纲,但不是实际流速;
  3. λ\lambda:沿程阻力系数;ζ\zeta:局部阻力系数。二者对应的损失形式不同。

目录#


第一部分:课程笔记#

1 流动阻力与水头损失#

1.1 阻力与能量损失的根源#

真实流体运动时会产生机械能损失,其根源可概括为:

  • fluid viscosity(流体黏性):相邻流层存在相对运动,产生黏性切应力;
  • turbulence(湍动):流体微团发生脉动、横向交换与涡旋混合,进一步耗散机械能。

机械能并没有消失,而是不可逆地转化为内能。工程上通常用单位重量流体损失的机械能表示,即 水头损失 hwh_w,单位为 m。

1.2 两类流动阻力#

沿程阻力与沿程水头损失#

Frictional resistance / frictional head loss 指流体沿着形状和尺寸基本不变的长直流段运动时,由壁面切应力持续作用而产生的阻力和能量损失。

记为:

hfh_f

主要特点:

  • 分布在整段流路上;
  • 与管长有关;
  • 在均匀流中,断面速度分布沿程不变,但总水头仍持续下降。

局部阻力与局部水头损失#

Local resistance / minor head loss 指管道形状、尺寸或流动方向突然改变时,由速度分布剧烈重组、流动分离和涡旋产生的阻力与能量损失。

记为:

hjh_j

典型位置:

  • 突然扩大、突然缩小;
  • 弯头、阀门、闸门;
  • 管道进口与出口;
  • 三通、分流、汇流等。
NOTE

“minor loss”中的 minor 只表示损失发生在局部构件处,不代表数值一定很小。短管、阀门很多的管路中,局部损失可能占主要部分。

1.3 水头损失叠加原理#

两断面之间的总水头损失等于所有沿程损失和局部损失之和:

hw=i=1nhf,i+k=1mhj,k\boxed{ h_w=\sum_{i=1}^{n}h_{f,i}+\sum_{k=1}^{m}h_{j,k} }

其中:

  • hwh_w:两断面之间的总水头损失,m;
  • hf,ih_{f,i}:第 ii 段等直径或等断面流段的沿程损失,m;
  • hj,kh_{j,k}:第 kk 个局部构件造成的局部损失,m;
  • nn:沿程流段数量;
  • mm:局部构件数量。

插图占位|图 6-1:插入课件《Flow resistance-fundamentals》第 7 页“总水头线、测压管水头线及弯头、突扩、突缩、闸门造成的损失”示意图。重点标出总水头线始终下降,局部构件处下降更明显。

1.4 Darcy–Weisbach 公式#

圆管恒定均匀流的沿程水头损失:

hf=λLdv22g\boxed{ h_f=\lambda\frac{L}{d}\frac{v^2}{2g} }

其中:

  • λ\lambda:Darcy friction factor,沿程阻力系数,无量纲;
  • LL:管长,m;
  • dd:圆管内径,m;
  • vv:断面平均流速,m/s;
  • gg:重力加速度,m/s2^2

对非圆形断面,引入水力半径:

R=AχR=\frac{A}{\chi}

于是:

hf=λL4Rv22g\boxed{ h_f=\lambda\frac{L}{4R}\frac{v^2}{2g} }

其中:

  • AA:过流断面面积;
  • χ\chi:湿周,即流体与固体边界接触的周长;
  • RR:水力半径。

圆管有:

R=πd2/4πd=d4R=\frac{\pi d^2/4}{\pi d}=\frac d4

因此 4R=d4R=d,可回到圆管公式。

1.5 局部损失的一般式#

hj=ζv22g\boxed{ h_j=\zeta\frac{v^2}{2g} }

其中:

  • ζ\zeta:local loss coefficient,局部阻力系数;
  • vv:与该 ζ\zeta 配套定义的特征流速。
WARNING

查表得到 ζ\zeta 后,必须确认它以哪一段流速为基准。对于突扩、突缩、三通等不同直径构件,选错流速会直接造成数量级错误。

1.6 均匀流基本方程#

取长度为 LL 的均匀流段,对其作沿流向受力平衡。由于均匀流前后断面速度分布相同,沿程损失等于测压管水头下降:

hf=JLh_f=JL

其中 JJ 为水力坡度:

J=hfLJ=\frac{h_f}{L}

对任意以原流线为边界的子流管:

τ=ρgRJ\boxed{ \tau=\rho gR'J }

对整个断面壁面:

τ0=ρgRJ\boxed{ \tau_0=\rho gRJ }

其中:

  • τ\tau:子流管边界上的切应力;
  • τ0\tau_0:壁面切应力;
  • R=A/χR'=A'/\chi':子流管的水力半径;
  • R=A/χR=A/\chi:整个断面的水力半径;
  • ρ\rho:流体密度。

这组公式只来自整体受力平衡,因此:

对圆管、非圆管、明渠,层流与湍流均成立。

1.7 摩阻流速#

定义 friction velocity(摩阻流速)

u=τ0ρ\boxed{ u_*=\sqrt{\frac{\tau_0}{\rho}} }

τ0=ρgRJ\tau_0=\rho gRJ 得:

u=gRJ\boxed{ u_*=\sqrt{gRJ} }

再由 Darcy–Weisbach 公式得:

u=vλ8\boxed{ u_*=v\sqrt{\frac{\lambda}{8}} }

uu_* 反映壁面摩擦强弱,是湍流近壁区无量纲化的重要速度尺度。

1.8 圆管断面切应力分布#

对半径为 rr 的同心流管:

R=πr22πr=r2R'=\frac{\pi r^2}{2\pi r}=\frac r2

于是:

τ=ρgr2J\tau=\rho g\frac r2J

在管壁 r=r0r=r_0

τ0=ρgr02J\tau_0=\rho g\frac{r_0}{2}J

两式相除:

ττ0=rr0\boxed{ \frac{\tau}{\tau_0}=\frac r{r_0} }

因此圆管均匀流中:

  • 中心 r=0r=0τ=0\tau=0
  • 向管壁线性增大;
  • 管壁 r=r0r=r_0τ=τ0\tau=\tau_0,达到最大值。

这个线性分布对圆管层流和湍流都成立,因为它来自整体动量平衡。


2 层流与湍流#

2.1 两种流态#

层流#

Laminar flow 中,流体质点分层、有序运动,各层之间横向混合很弱。

特点:

  • 运动规则;
  • 脉动很弱;
  • 黏性力占主导;
  • 水头损失与平均速度一次方成正比。

湍流#

Turbulent flow 中,速度、压力等物理量在时间和空间上发生不规则脉动,流体微团强烈掺混。

特点:

  • 随机性和无序性;
  • 强烈横向动量交换;
  • 涡旋众多;
  • 能量耗散较强;
  • 工程流动中非常常见。

2.2 Reynolds experiment#

雷诺实验通过向玻璃管中注入细染色液束观察流态:

  • 低流速:染色线保持清晰,说明为层流;
  • 增大流速:染色线摆动并扩散,进入过渡状态;
  • 高流速:染色液迅速混合,成为湍流。

实验还发现,升速和降速时发生转捩的速度不同:

  • upper critical Reynolds number:升速过程中层流失稳的上临界值,受扰动影响较大;
  • lower critical Reynolds number:降速过程中湍流恢复层流的下临界值,较稳定。

工程上使用下临界雷诺数作为判别标准。

插图占位|图 6-2:插入课件《Flow resistance-fundamentals》第 17 页雷诺实验装置与染色线形态,并插入第 19 页 lghflgv\lg h_f-\lg v 图。

2.3 水头损失与流速的关系#

实验可写成:

hf=kvmh_f=kv^m

取对数:

lghf=lgk+mlgv\lg h_f=\lg k+m\lg v

其中斜率 mm 反映流态:

  • 层流:m=1m=1,即 hfvh_f\propto v
  • 湍流:通常 m1.752.0m\approx1.75\sim2.0
  • 水力粗糙区:λ\lambdaReRe 无关,故 hfv2h_f\propto v^2

2.4 Reynolds number#

圆管流动:

Re=vdν=ρvdμ\boxed{ Re=\frac{vd}{\nu}=\frac{\rho vd}{\mu} }

其中:

  • vv:断面平均速度;
  • dd:圆管直径;
  • ν\nu:运动黏度;
  • μ\mu:动力黏度。

物理意义:

Re惯性力黏性力\boxed{ Re\sim\frac{\text{惯性力}}{\text{黏性力}} }

圆管判别:

Re<2300层流\boxed{ Re<2300\quad\Rightarrow\quad\text{层流} }2300Re4000过渡区2300\lesssim Re\lesssim4000\quad\Rightarrow\quad\text{过渡区}Re>4000通常按湍流处理Re>4000\quad\Rightarrow\quad\text{通常按湍流处理}

对于非圆管有压流,可使用水力直径:

de=4Rd_e=4R

并写为:

Re=vdeνRe=\frac{vd_e}{\nu}

课件对明渠采用:

ReR=vRνRe_R=\frac{vR}{\nu}

对应下临界值约为:

ReR,c=575Re_{R,c}=575

因为 de=4Rd_e=4R,所以 2300/4=5752300/4=575,两种写法本质一致。

2.5 为什么用下临界雷诺数#

上临界值受入口扰动、管壁振动、来流脉动等影响很大;下临界值较稳定,主要由断面形状决定。

因此判断时使用下临界值:

  • Re<RecRe<Re_c:一定为层流;
  • Re>RecRe>Re_c:层流已可能失稳,继续结合数值范围判断过渡或湍流。

2.6 课堂例题:自来水管流态#

已知:

d=0.1 m,v=1.0 m/s,ν=1.31×106 m2/sd=0.1\ \text{m},\qquad v=1.0\ \text{m/s},\qquad \nu=1.31\times10^{-6}\ \text{m}^2/\text{s}

则:

Re=vdν=1.0×0.11.31×106=7.63×104Re=\frac{vd}{\nu} =\frac{1.0\times0.1}{1.31\times10^{-6}} =7.63\times10^4

故为湍流。

保持层流的最大平均速度:

vmax=Recνd=2300×1.31×1060.10.030 m/sv_{\max}=\frac{Re_c\nu}{d} =\frac{2300\times1.31\times10^{-6}}{0.1} \approx0.030\ \text{m/s}

2.7 课堂例题:矩形明渠流态#

已知:

b=0.2 m,h=0.1 m,v=0.12 m/sb=0.2\ \text{m},\quad h=0.1\ \text{m},\quad v=0.12\ \text{m/s}

过流面积与湿周:

A=bh=0.02 m2A=bh=0.02\ \text{m}^2χ=b+2h=0.4 m\chi=b+2h=0.4\ \text{m}

水力半径:

R=Aχ=0.05 mR=\frac A\chi=0.05\ \text{m}

20C20^\circ\mathrm C 水的 ν1.0×106 m2/s\nu\approx1.0\times10^{-6}\ \text{m}^2/\text{s}

ReR=vRν=0.12×0.051.0×106=6000>575Re_R=\frac{vR}{\nu} =\frac{0.12\times0.05}{1.0\times10^{-6}} =6000>575

故为湍流。


3 圆管层流#

3.1 推导前提#

圆管层流速度分布的推导基于:

  • 恒定流;
  • 不可压缩流体;
  • 圆管内充分发展均匀层流;
  • 牛顿流体;
  • 轴对称流动;
  • 管壁满足无滑移条件。
WARNING

入口段速度分布仍在发展,不满足充分发展条件,因此本节抛物线速度分布与 λ=64/Re\lambda=64/Re 不能直接用于入口附近。

3.2 速度分布#

设:

  • 管半径为 r0r_0
  • 距管轴半径为 rr
  • 距壁面距离为 y=r0ry=r_0-r
  • 轴向速度为 u(r)u(r)

切应力分布:

τ=12ρgJr\tau=\frac12\rho gJr

牛顿内摩擦定律:

τ=μdudy=μdudr\tau=\mu\frac{du}{dy}=-\mu\frac{du}{dr}

联立:

μdudr=12ρgJr-\mu\frac{du}{dr}=\frac12\rho gJr

积分:

u=ρgJ4μr2+Cu=-\frac{\rho gJ}{4\mu}r^2+C

管壁无滑移:

r=r0,u=0r=r_0,\qquad u=0

得到:

u(r)=ρgJ4μ(r02r2)\boxed{ u(r)=\frac{\rho gJ}{4\mu}\left(r_0^2-r^2\right) }

速度沿半径呈抛物线分布;在三维中形成旋转抛物面。

插图占位|图 6-3:插入课件《Flow resistance-fundamentals》第 28–29 页圆管层流速度抛物线与切应力线性分布图。标出中心速度最大、壁面速度为零。

3.3 最大速度与断面平均速度#

中心 r=0r=0

umax=ρgJ4μr02\boxed{ u_{\max}=\frac{\rho gJ}{4\mu}r_0^2 }

因此:

u=umax[1(rr0)2]\boxed{ u=u_{\max}\left[1-\left(\frac r{r_0}\right)^2\right] }

断面平均速度:

v=1AAudAv=\frac1A\int_Au\,dA

积分可得:

v=12umax\boxed{ v=\frac12u_{\max} }

即:

umax=2v\boxed{ u_{\max}=2v }

也可将速度分布写成:

u=2v[1(rr0)2]\boxed{ u=2v\left[1-\left(\frac r{r_0}\right)^2\right] }

3.4 动能与动量修正系数#

圆管层流速度分布很不均匀:

α=2\boxed{ \alpha=2 }β=43\boxed{ \beta=\frac43 }

其中:

  • α\alpha:kinetic-energy correction factor,动能修正系数;
  • β\beta:momentum correction factor,动量修正系数。

3.5 Hagen–Poiseuille 关系#

由平均速度:

v=ρgJr028μv=\frac{\rho gJr_0^2}{8\mu}

J=hf/LJ=h_f/Lr0=d/2r_0=d/2

hf=32μLvρgd2h_f=\frac{32\mu Lv}{\rho gd^2}

压强降形式:

Δp=32μLvd2\boxed{ \Delta p=\frac{32\mu Lv}{d^2} }

流量形式:

Q=vπd24Q=v\frac{\pi d^2}{4}

所以:

Q=πd4128μΔpL\boxed{ Q=\frac{\pi d^4}{128\mu}\frac{\Delta p}{L} }

这说明在其他条件不变时:

Qd4Q\propto d^4

因此管径增大一倍,层流流量增大为原来的 24=162^4=16 倍。

3.6 层流沿程阻力系数#

将层流损失与 Darcy–Weisbach 公式对比:

32μLvρgd2=λLdv22g\frac{32\mu Lv}{\rho gd^2} =\lambda\frac Ld\frac{v^2}{2g}

得到:

λ=64Re\boxed{ \lambda=\frac{64}{Re} }

因此层流中:

  • λ\lambda 只与 ReRe 有关;
  • 与相对粗糙度 Δ/d\Delta/d 无关;
  • hfvh_f\propto v

3.7 壁面切应力#

由摩阻系数关系:

τ0=λ8ρv2\boxed{ \tau_0=\frac{\lambda}{8}\rho v^2 }

层流中代入 λ=64/Re\lambda=64/Re

τ0=8μvd\tau_0=\frac{8\mu v}{d}

3.8 课堂例题:油在圆管内层流#

已知:

ρ=850 kg/m3,ν=0.18×104 m2/s\rho=850\ \text{kg/m}^3,\quad \nu=0.18\times10^{-4}\ \text{m}^2/\text{s}d=0.1 m,v=0.0635 m/sd=0.1\ \text{m},\quad v=0.0635\ \text{m/s}
  1. 中心最大速度:
umax=2v=0.127 m/su_{\max}=2v=0.127\ \text{m/s}
  1. 距中心 r=20 mmr=20\ \text{mm}r0=50 mmr_0=50\ \text{mm}
u=0.127[1(0.020.05)2]=0.107 m/su=0.127\left[1-\left(\frac{0.02}{0.05}\right)^2\right] =0.107\ \text{m/s}
  1. 雷诺数与阻力系数:
Re=vdν=0.0635×0.10.18×104353Re=\frac{vd}{\nu} =\frac{0.0635\times0.1}{0.18\times10^{-4}} \approx353λ=643530.181\lambda=\frac{64}{353}\approx0.181
  1. 壁面切应力:
τ0=λ8ρv20.077 Pa\tau_0=\frac{\lambda}{8}\rho v^2 \approx0.077\ \text{Pa}
  1. 每千米管长损失:
hf=λ10000.10.063522g0.37 mh_f=\lambda\frac{1000}{0.1}\frac{0.0635^2}{2g} \approx0.37\ \text{m}

3.9 课堂例题:毛细管黏度计#

已知:

d=6 mm,L=2 m,Q=77×106 m3/sd=6\ \text{mm},\quad L=2\ \text{m},\quad Q=77\times10^{-6}\ \text{m}^3/\text{s}

水银压差计读数:

hp=0.3 m,ρoil=900 kg/m3h_p=0.3\ \text{m},\quad \rho_{\mathrm{oil}}=900\ \text{kg/m}^3

水银密度取 13600 kg/m313600\ \text{kg/m}^3

压差对应油柱水头:

hf=(ρHgρ1)hp4.23 mh_f=\left(\frac{\rho_{Hg}}{\rho}-1\right)h_p \approx4.23\ \text{m}

平均速度:

v=4Qπd22.72 m/sv=\frac{4Q}{\pi d^2}\approx2.72\ \text{m/s}

层流关系:

hf=64ReLdv22g=32νLvgd2h_f=\frac{64}{Re}\frac Ld\frac{v^2}{2g} =\frac{32\nu Lv}{gd^2}

解得:

ν8.6×106 m2/s\boxed{ \nu\approx8.6\times10^{-6}\ \text{m}^2/\text{s} }μ=ρν7.7×103 Pas\boxed{ \mu=\rho\nu\approx7.7\times10^{-3}\ \text{Pa}\cdot\text{s} }

校核:

Re=vdν1.91×103<2300Re=\frac{vd}{\nu}\approx1.91\times10^3<2300

假设层流成立。


4 湍流基本理论#

4.1 湍流断面结构#

圆管湍流断面可分为:

  1. viscous sublayer(黏性底层):紧贴壁面,流体横向脉动受到壁面抑制,黏性切应力占主导;
  2. turbulent core region(湍流核心区):远离壁面,湍动交换强烈,湍流切应力占主导。
NOTE

“黏性底层”是湍流内部的一小层,并不意味着整个管流属于层流。判定整个管流流态要看断面平均速度和雷诺数。

插图占位|图 6-4:插入课件《Flow resistance-Turbulence》第 3–5 页“黏性底层与湍流核心区”示意图。

4.2 瞬时量、时均量与脉动量#

以速度为例:

ui=uˉi+ui\boxed{ u_i=\bar u_i+u_i' }

其中:

  • uiu_i:instantaneous velocity,瞬时速度;
  • uˉi\bar u_i:time-averaged velocity,时均速度;
  • uiu_i':fluctuating velocity,脉动速度。

时均速度:

uˉi=1T0Tui(t)dt\boxed{ \bar u_i=\frac1T\int_0^Tu_i(t)\,dt }

脉动量的时均值为零:

ui=0\overline{u_i'}=0

压力同样可分解:

p=pˉ+pp=\bar p+p'

4.3 湍流中的恒定流#

湍流瞬时速度始终在波动,因此不能用“瞬时速度是否变化”判断恒定流。

若时均量不随时间改变:

uˉt=0\frac{\partial\bar u}{\partial t}=0

则称为恒定湍流。

4.4 湍流强度#

常用湍流强度定义为脉动速度均方根与平均速度之比:

I=ux2+uy2+uz23uˉ\boxed{ I=\frac{\sqrt{\dfrac{\overline{u_x'^2}+\overline{u_y'^2}+\overline{u_z'^2}}{3}}}{\bar u} }

它反映速度脉动相对于平均运动的强弱。

4.5 湍流总切应力#

湍流中总切应力由两部分组成:

τ=τν+τt\boxed{ \tau=\tau_\nu+\tau_t }

黏性切应力#

τν=μduˉdy\boxed{ \tau_\nu=\mu\frac{d\bar u}{dy} }

它与动力黏度和时均速度梯度有关,近壁区占主导。

Reynolds stress#

二维情况下:

τt=ρuv\boxed{ \tau_t=-\rho\overline{u'v'} }

它来自脉动造成的横向质量和动量交换。

物理意义:

  • 与密度和脉动强度有关;
  • 不直接由分子黏度决定;
  • 在湍流核心区通常占主导。

负号来自典型脉动相关性:高动量流体向低速区运动、低动量流体向高速区运动,使 uu'vv' 往往异号。

4.6 Prandtl mixing-length theory#

Prandtl 假设流体微团在横向移动一段距离后,才与周围流体充分混合。这一特征距离称为 mixing length ll

近壁区取:

l=κy\boxed{ l=\kappa y }

其中:

  • yy:距壁面距离;
  • κ0.4\kappa\approx0.4:Kármán 常数。

湍流切应力可近似写成:

τt=ρl2duˉdyduˉdy\boxed{ \tau_t=\rho l^2\left|\frac{d\bar u}{dy}\right|\frac{d\bar u}{dy} }

若速度梯度为正,可简写为:

τt=ρl2(duˉdy)2\tau_t=\rho l^2\left(\frac{d\bar u}{dy}\right)^2

定义涡黏度:

η=ρl2duˉdy\boxed{ \eta=\rho l^2\left|\frac{d\bar u}{dy}\right| }

于是:

τ=(μ+η)duˉdy\boxed{ \tau=(\mu+\eta)\frac{d\bar u}{dy} }

η\eta 不是流体本身的物性参数,它依赖流态、位置和湍动强度。

4.7 黏性底层速度分布#

黏性底层中湍流切应力可忽略,且切应力近似等于壁面切应力:

τ0μduˉdy\tau_0\approx\mu\frac{d\bar u}{dy}

积分并用壁面无滑移条件:

uˉ=τ0μy\bar u=\frac{\tau_0}{\mu}y

引入摩阻流速 uu_* 和运动黏度 ν\nu

uˉu=yuν\boxed{ \frac{\bar u}{u_*}=\frac{yu_*}{\nu} }

令:

u+=uˉu,y+=yuνu^+=\frac{\bar u}{u_*},\qquad y^+=\frac{yu_*}{\nu}

则:

u+=y+\boxed{ u^+=y^+}

即黏性底层中速度沿 yy 近似线性分布。

4.8 黏性底层厚度#

实验取黏性底层边缘:

y+11.6y^+\approx11.6

所以:

δl=11.6νu\boxed{ \delta_l=\frac{11.6\nu}{u_*} }

u=vλ/8u_*=v\sqrt{\lambda/8}

δl=32.8dReλ\boxed{ \delta_l=\frac{32.8d}{Re\sqrt\lambda} }

结论:

  • 黏性底层通常很薄;
  • 在管径和流体性质不变时,流速增大、ReRe 增大,δl\delta_l 变薄;
  • 粗糙凸起是否伸出黏性底层,决定管壁在水力学上表现为光滑或粗糙。

4.9 湍流核心区速度分布#

在近壁湍流核心区,取:

ττ0,l=κy\tau\approx\tau_0,\qquad l=\kappa y

可得对数律:

uˉu=1κlny+C\boxed{ \frac{\bar u}{u_*}=\frac1\kappa\ln y+C }

严格的无量纲形式为:

uˉu=1κln(yuν)+B\boxed{ \frac{\bar u}{u_*}=\frac1\kappa\ln\left(\frac{yu_*}{\nu}\right)+B }

光滑管经验式:

uˉu=5.75lg(yuν)+5.5\boxed{ \frac{\bar u}{u_*}=5.75\lg\left(\frac{yu_*}{\nu}\right)+5.5 }

粗糙管经验式:

uˉu=5.75lg(yΔ)+8.48\boxed{ \frac{\bar u}{u_*}=5.75\lg\left(\frac y\Delta\right)+8.48 }

其中 Δ\Delta 为等效粗糙高度。

4.10 指数速度分布#

圆管湍流也常用幂律近似:

uumax=(yr0)1/n\boxed{ \frac{u}{u_{\max}}=\left(\frac y{r_0}\right)^{1/n} }

经典的 one-seventh power law

uumax=(yr0)1/7\frac{u}{u_{\max}}=\left(\frac y{r_0}\right)^{1/7}

湍流速度分布比层流更均匀:

  • 层流:umax=2vu_{\max}=2v
  • 湍流:umaxu_{\max} 通常约为 1.15v1.26v1.15v\sim1.26v

5 圆管湍流沿程损失#

5.1 绝对粗糙度与相对粗糙度#

  • Δ\Delta:absolute roughness,绝对粗糙度;
  • Δ/d\Delta/d:relative roughness,相对粗糙度。

相同粗糙高度对小管影响更明显,因此实际判断要使用 Δ/d\Delta/d

5.2 Nikuradse experiment#

Nikuradse 在人工均匀砂粒粗糙管中测量 λ\lambda,得到不同 ReReΔ/d\Delta/d 下的阻力规律。

λ\lambda 的控制因素可分为五个区域:

  1. 层流区: λ=64Re\lambda=\frac{64}{Re}
  2. 层流向湍流过渡区:λ\lambda 不稳定;
  3. 水力光滑湍流区: λ=f(Re)\lambda=f(Re)
  4. 光滑向粗糙的过渡湍流区: λ=f(Re,Δd)\lambda=f\left(Re,\frac\Delta d\right)
  5. 水力粗糙湍流区: λ=f(Δd)\lambda=f\left(\frac\Delta d\right)

插图占位|图 6-5:插入课件《Flow resistance-Turbulence》第 21–22 页 Nikuradse 图及粗糙凸起与黏性底层的三种相对关系。

5.3 三个湍流阻力区的物理解释#

水力光滑区#

粗糙凸起完全埋在黏性底层中,主流感受不到粗糙度:

λ=f(Re)\lambda=f(Re)

过渡粗糙区#

一部分粗糙凸起伸出黏性底层:

λ=f(Re,Δd)\lambda=f\left(Re,\frac\Delta d\right)

水力粗糙区#

粗糙凸起明显穿出黏性底层,形状阻力占主导:

λ=f(Δd)\lambda=f\left(\frac\Delta d\right)

此时 λ\lambdaReRe 无关,所以:

hfv2h_f\propto v^2

5.4 光滑管与完全粗糙管公式#

光滑湍流区:

1λ=2lg(2.51Reλ)\boxed{ \frac1{\sqrt\lambda} =-2\lg\left(\frac{2.51}{Re\sqrt\lambda}\right) }

完全粗糙区:

1λ=2lg(Δ3.7d)\boxed{ \frac1{\sqrt\lambda} =-2\lg\left(\frac{\Delta}{3.7d}\right) }

5.5 Colebrook–White equation#

适用于圆管湍流的光滑区、过渡粗糙区和粗糙区:

1λ=2lg(2.51Reλ+Δ3.7d)\boxed{ \frac1{\sqrt\lambda} =-2\lg\left( \frac{2.51}{Re\sqrt\lambda} +\frac{\Delta}{3.7d} \right) }

它是隐式方程,需要迭代或借助 Moody 图。

5.6 Barr 显式公式#

课件采用:

1λ=2lg(5.1286Re0.89+Δ3.7d)\boxed{ \frac1{\sqrt\lambda} =-2\lg\left( \frac{5.1286}{Re^{0.89}} +\frac{\Delta}{3.7d} \right) }

它可直接计算,适合考试手算。

5.7 Moody diagram#

Moody 图以:

  • 横轴 ReRe
  • 纵轴 λ\lambda
  • 曲线参数 Δ/d\Delta/d

表示圆管阻力规律。

读图步骤:

  1. 计算 ReRe
  2. 计算 Δ/d\Delta/d
  3. 找到对应相对粗糙度曲线;
  4. 在给定 ReRe 处读取 λ\lambda
  5. 代入 Darcy–Weisbach 公式。

插图占位|图 6-6:插入课件《Flow resistance-Turbulence》第 28 页 Moody diagram。建议在图旁标注层流直线 64/Re64/Re、光滑管曲线和完全粗糙区水平趋势。

5.8 课件给出的阻力分区判别#

课件依据 Moody 图给出一组工程判别边界:

  1. Re4000Re\le 4000 则处于层流区或层流向湍流的过渡区;
  2. Δd0.0008,4000<Re10(Δd)1\frac\Delta d\le0.0008, \qquad 4000<Re\le10\left(\frac\Delta d\right)^{-1} 可按水力光滑湍流区处理;
  3. Re>576.12(Δd)1.119Re>576.12\left(\frac\Delta d\right)^{-1.119} 可按水力粗糙湍流区处理;
  4. 其余情况属于光滑区向粗糙区的过渡区。

这些阈值是根据 Moody 图整理的经验边界。计算时直接使用 Moody 图或 Colebrook–White 公式通常更稳妥,不建议脱离课件强行死记。

5.9 圆管湍流计算流程#

给定 Q,d,L,ν,ΔQ,d,L,\nu,\Delta

  1. 平均速度: v=4Qπd2v=\frac{4Q}{\pi d^2}
  2. 雷诺数: Re=vdνRe=\frac{vd}{\nu}
  3. 判断流态;
  4. 若为层流:λ=64/Re\lambda=64/Re
  5. 若为湍流:用 Moody、Colebrook 或 Barr 求 λ\lambda
  6. 计算: hf=λLdv22gh_f=\lambda\frac Ld\frac{v^2}{2g}

5.10 课堂例题:粗糙圆管沿程损失#

已知:

L=500 m,d=0.2 m,Δ=0.1 mmL=500\ \text{m},\quad d=0.2\ \text{m},\quad \Delta=0.1\ \text{mm}Q=10×103 m3/s,ν=1.31×106 m2/sQ=10\times10^{-3}\ \text{m}^3/\text{s},\quad \nu=1.31\times10^{-6}\ \text{m}^2/\text{s}

平均速度:

v=Qπd2/4=0.318 m/sv=\frac{Q}{\pi d^2/4}=0.318\ \text{m/s}

雷诺数:

Re=vdν4.86×104Re=\frac{vd}{\nu}\approx4.86\times10^4

相对粗糙度:

Δd=0.00010.2=5.0×104\frac\Delta d=\frac{0.0001}{0.2}=5.0\times10^{-4}

用 Barr 公式:

λ0.0227\lambda\approx0.0227

沿程损失:

hf=0.02275000.20.31822g0.293 mh_f=0.0227\frac{500}{0.2}\frac{0.318^2}{2g} \approx0.293\ \text{m}

6 非圆管与明渠沿程损失#

6.1 非圆形有压管#

定义等效直径:

de=4R=4Aχ\boxed{ d_e=4R=\frac{4A}{\chi} }

湍流沿程损失:

hf=λLdev22g\boxed{ h_f=\lambda\frac{L}{d_e}\frac{v^2}{2g} }

课件给出的近似适用范围:

  • 矩形管:长宽比小于约 8;
  • 椭圆管:长短轴比小于约 3。
WARNING

直接用 de=4Rd_e=4R 代替圆管直径主要用于湍流。非圆管层流的阻力系数还与断面形状有关,不能一概使用 64/Re64/Re

6.2 Chezy formula#

明渠均匀流:

v=CRJ\boxed{ v=C\sqrt{RJ} }

其中:

  • CC:Chezy coefficient,谢才系数,量纲为 m1/2^{1/2}/s;
  • R=A/χR=A/\chi:水力半径;
  • JJ:水力坡度,均匀流中等于底坡。

与 Darcy 阻力系数的关系:

C=8gλ\boxed{ C=\sqrt{\frac{8g}{\lambda}} }

6.3 Manning formula#

C=1nR1/6C=\frac1nR^{1/6}

代入 Chezy 公式:

v=1nR2/3J1/2\boxed{ v=\frac1nR^{2/3}J^{1/2} }

流量:

Q=Av=AnR2/3J1/2\boxed{ Q=Av=\frac A nR^{2/3}J^{1/2} }

其中 nn 为 Manning roughness coefficient,曼宁糙率。

NOTE

Manning 公式是经验公式,nn 的数值与所采用的单位制配套。本课程公式和表值均按 SI 制使用。

课件给出的常见范围:

  • 铸铁管、钢管:n0.0110.013n\approx0.011\sim0.013
  • 混凝土管:n0.0110.014n\approx0.011\sim0.014

6.4 课堂例题:新供水管流量#

已知:

d=0.4 m,L=100 m,n=0.011,hf=0.4 md=0.4\ \text{m},\quad L=100\ \text{m},\quad n=0.011,\quad h_f=0.4\ \text{m}

圆管满流:

A=πd24=0.126 m2A=\frac{\pi d^2}{4}=0.126\ \text{m}^2R=d4=0.1 mR=\frac d4=0.1\ \text{m}J=hfL=0.004J=\frac{h_f}{L}=0.004

Manning 公式:

v=10.011(0.1)2/3(0.004)1/21.24 m/sv=\frac1{0.011}(0.1)^{2/3}(0.004)^{1/2} \approx1.24\ \text{m/s}Q=Av0.156 m3/sQ=Av\approx0.156\ \text{m}^3/\text{s}

7 局部水头损失#

7.1 一般表达式#

hj=ζv22g\boxed{ h_j=\zeta\frac{v^2}{2g} }

湍流中,局部损失通常与速度平方成正比。

ζ\zeta 主要受以下因素影响:

  • 局部构件几何形状;
  • 面积比;
  • 转弯角度;
  • 阀门开度;
  • 选取的基准流速;
  • 某些情况下还受 ReRe 影响。

7.2 突然扩大损失#

设水平管道从断面 1 突然扩大到断面 2,且:

A2>A1,v1>v2A_2>A_1,\qquad v_1>v_2

动量方程、能量方程和连续方程联立得到 Borda–Carnot 公式:

hj=(v1v2)22g\boxed{ h_j=\frac{(v_1-v_2)^2}{2g} }

利用连续方程:

v1A1=v2A2v_1A_1=v_2A_2

若以 v1v_1 为基准:

hj=(1A1A2)2v122g\boxed{ h_j=\left(1-\frac{A_1}{A_2}\right)^2\frac{v_1^2}{2g} }ζ1=(1A1A2)2\boxed{ \zeta_1=\left(1-\frac{A_1}{A_2}\right)^2 }

若以 v2v_2 为基准:

hj=(A2A11)2v222g\boxed{ h_j=\left(\frac{A_2}{A_1}-1\right)^2\frac{v_2^2}{2g} }ζ2=(A2A11)2\boxed{ \zeta_2=\left(\frac{A_2}{A_1}-1\right)^2 }

插图占位|图 6-7:插入课件《Flow resistance-local loss》第 3 页突扩控制体图,标出 p1A1p_1A_1p2A2p_2A_2、回流区及 v1,v2v_1,v_2

7.3 突然缩小损失#

课件给出的近似式,以缩小后速度 v2v_2 为基准:

hj=ζv222g\boxed{ h_j=\zeta\frac{v_2^2}{2g} }ζ0.5(1A2A1)\boxed{ \zeta\approx0.5\left(1-\frac{A_2}{A_1}\right) }

突然缩小后先形成收缩断面,再发生扩散和混合,因此损失机理与反向突扩不同。

7.4 管道出口#

管道排入很大的水池或容器:

ApipeApool0\frac{A_{pipe}}{A_{pool}}\approx0

可看成突然扩大:

ζout=1.0\boxed{ \zeta_{out}=1.0 }hout=v22g\boxed{ h_{out}=\frac{v^2}{2g} }

物理意义:管内速度水头几乎全部耗散。

7.5 管道进口#

课件给出的典型值:

  • 直角进口: ζ0.5\zeta\approx0.5
  • 圆角进口: ζ0.050.25\zeta\approx0.05\sim0.25
  • 外伸进口: ζ1.0\zeta\approx1.0

进口越接近流线型,流动分离越弱,损失越小。

7.6 多级突扩的优化#

速度由 v1v_1 降到 v2v_2,中间加入速度为 vv 的管段,两次突扩总损失:

hj=(v1v)22g+(vv2)22gh_j=\frac{(v_1-v)^2}{2g}+\frac{(v-v_2)^2}{2g}

令:

dhjdv=0\frac{dh_j}{dv}=0

得到:

v=v1+v22\boxed{ v=\frac{v_1+v_2}{2} }

此时:

hj,min=(v1v2)24g\boxed{ h_{j,\min}=\frac{(v_1-v_2)^2}{4g} }

单次突扩:

hj,1=(v1v2)22gh_{j,1}=\frac{(v_1-v_2)^2}{2g}

所以两级等速度差突扩的最小损失为单次突扩的:

12\boxed{ \frac12 }

7.7 减小局部损失的方法#

核心原则:

让固体边界变化更平缓,使流线转弯和断面变化更连续。

例如:

  • 突扩改为渐扩;
  • 突缩改为圆滑收缩;
  • 弯头增大曲率半径;
  • 进口做成圆角;
  • 避免阀门局部强烈节流。

8 边界层#

8.1 边界层概念#

当大雷诺数流体沿固体表面流动时,黏性影响主要集中在紧邻壁面的薄层内,该薄层称为 boundary layer(边界层)

  • 边界层内:速度梯度大,黏性不可忽略;
  • 边界层外:速度梯度很小,可近似按理想流体或势流处理。

这使真实高 ReRe 流动可分区求解。

8.2 边界层厚度#

通常定义边界层边缘速度达到外部来流速度的约 99%99\%

u(y=δ)0.99Uu(y=\delta)\approx0.99U_\infty

其中:

  • UU_\infty:边界层外来流速度;
  • δ\delta:边界层厚度。

边界层沿流向逐渐增厚。

8.3 层流边界层与湍流边界层#

边界层可经历:

  1. 层流边界层;
  2. 转捩区;
  3. 湍流边界层。

局部雷诺数:

Rex=Uxν\boxed{ Re_x=\frac{U_\infty x}{\nu} }

课件给出的光滑平板临界范围:

3×105<Rex,cr<3×106\boxed{ 3\times10^5<Re_{x,cr}<3\times10^6 }

临界值受来流扰动强度影响,并非固定常数。

湍流边界层最靠近壁面处仍存在黏性底层。

插图占位|图 6-8:插入课件《Flow resistance-local loss》第 15–17 页光滑平板上层流边界层、转捩区、湍流边界层和黏性底层示意图。

8.4 光滑平板边界层厚度近似式#

零压强梯度光滑平板:

层流边界层:

δ5xRex\boxed{ \delta\approx\frac{5x}{\sqrt{Re_x}} }

湍流边界层:

δ0.377xRex1/5\boxed{ \delta\approx\frac{0.377x}{Re_x^{1/5}} }

这些是经验近似式,使用前要确认流态和适用条件。

8.5 压强梯度与边界层分离#

若沿流向压强不变:

dpdx=0\frac{dp}{dx}=0

平板边界层通常只会持续增厚,不发生分离。

若沿流向压强升高:

dpdx>0\boxed{ \frac{dp}{dx}>0 }

称为 adverse pressure gradient(逆压梯度)

流体沿流向减速,近壁低动量流体受到黏性阻力后可能停滞并反向运动,导致边界层脱离壁面。

分离点满足:

τ0=μ(uy)y=0=0\boxed{ \tau_0=\mu\left(\frac{\partial u}{\partial y}\right)_{y=0}=0 }

分离点后出现回流和涡旋。

插图占位|图 6-9:插入课件《Flow resistance-local loss》第 19–20 页逆压梯度下边界层分离图。标出分离点、τ0=0\tau_0=0、回流区。

8.6 Wake#

边界层分离后,分离流线与物体后部包围的区域称为 wake(尾流区)

尾流的影响:

  • 增大绕流阻力;
  • 消耗机械能;
  • 产生冲刷或淤积;
  • 交替脱落的涡旋可能诱发结构振动。

减小尾流的主要方法:

  • 将物体外形设计得更流线化;
  • 减缓逆压梯度;
  • 推迟边界层分离。

9 绕流阻力与球体沉降#

9.1 Drag 与 lift#

固体与流体存在相对运动时,流体对固体的作用力可分解为:

  • drag force FDF_D:平行于来流方向;
  • lift force FLF_L:垂直于来流方向。

9.2 Kármán vortex street#

钝体绕流发生边界层分离后,两侧涡旋可能交替脱落,形成 Kármán vortex street(卡门涡街)

若涡脱落频率接近结构固有频率,可能发生共振,造成:

  • 立管振动;
  • 桥梁拉索风致振动;
  • 构件疲劳破坏。

插图占位|图 6-10:插入课件《Flow resistance-Cd》第 4 页卡门涡街图。

9.3 绕流阻力公式#

FD=CDρu022A\boxed{ F_D=C_D\frac{\rho u_0^2}{2}A }

其中:

  • FDF_D:绕流阻力,N;
  • CDC_D:drag coefficient,阻力系数;
  • ρ\rho:流体密度;
  • u0u_0:未受扰动来流与物体的相对速度;
  • AA:物体迎流投影面积。

绕流阻力可理解为两部分的合力:

  • skin-friction drag:物体表面切应力形成的摩擦阻力;
  • pressure drag / form drag:前后压强不对称形成的压差阻力。

流线型物体的摩擦阻力相对重要;钝体发生分离并形成大尾流后,压差阻力通常占主导。

CDC_D 主要取决于:

  • 物体形状;
  • 雷诺数;
  • 表面粗糙度;
  • 分离与尾流状态。

9.4 小球阻力系数#

球体:

Re=u0dνRe=\frac{u_0d}{\nu}

课件给出的分段近似:

CD={24Re,Re<1,13Re,1<Re<103,0.44,103<Re<2×105.\boxed{ C_D= \begin{cases} \dfrac{24}{Re}, & Re<1,\\[6pt] \dfrac{13}{\sqrt{Re}}, & 1<Re<10^3,\\[6pt] 0.44, & 10^3<Re<2\times10^5. \end{cases} }

9.5 Stokes drag#

Re<1Re<1

CD=24ReC_D=\frac{24}{Re}

代入一般阻力公式,球的迎流面积 A=πd2/4A=\pi d^2/4

FD=3πμdu0\boxed{ F_D=3\pi\mu du_0 }

9.6 球体终端沉降速度#

球体下沉达到稳定时:

G=FB+FDG=F_B+F_D

即:

ρsgπd36=hogπd36+3πμdu0\rho_sg\frac{\pi d^3}{6} = ho g\frac{\pi d^3}{6}+3\pi\mu du_0

解得:

u0=d2g(ρsρ)18μ\boxed{ u_0=\frac{d^2g(\rho_s-\rho)}{18\mu} }

其中:

  • ρs\rho_s:球体密度;
  • ρ\rho:流体密度;
  • dd:球径;
  • u0u_0:球相对于流体的终端沉降速度。
WARNING

该式由 Stokes 阻力推导,必须校核:

Re=u0dν<1Re=\frac{u_0d}{\nu}<1

若不满足,应改用相应 CD(Re)C_D(Re) 关系重新求解。


10 本章综合整理#

10.1 计算题总流程#

遇到管路能量损失问题:

  1. 画流动系统,划分等直径流段与局部构件;
  2. 由连续方程求各段平均速度;
  3. 计算各段 ReRe
  4. 判断层流或湍流;
  5. 求各段 λ\lambda
  6. 计算沿程损失: hf=λLdv22gh_f=\lambda\frac Ld\frac{v^2}{2g}
  7. 查取各局部构件 ζ\zeta,确认基准速度;
  8. 计算局部损失: hj=ζv22gh_j=\zeta\frac{v^2}{2g}
  9. 叠加: hw=hf+hjh_w=\sum h_f+\sum h_j
  10. 与 Bernoulli 方程联立求未知流量、压强或所需作用水头。

10.2 圆管层流与湍流对比#

项目圆管层流圆管湍流
判别Re<2300Re<2300通常 Re>4000Re>4000
主要动量传递分子黏性湍动交换 + 分子黏性
速度分布抛物线,较不均匀对数律/幂律,较均匀
umax/vu_{\max}/v221.151.261.15\sim1.26
切应力总分布由中心向壁面线性增大同样线性增大
阻力系数λ=64/Re\lambda=64/Reλ=f(Re,Δ/d)\lambda=f(Re,\Delta/d)
粗糙度影响视阻力区而定
水头损失与速度hfvh_f\propto vhfv1.752h_f\propto v^{1.75\sim2}
修正系数α=2,β=4/3\alpha=2,\beta=4/3工程中常近似 αβ1\alpha\approx\beta\approx1

10.3 高频易错点#

  1. 圆管切应力在中心为零,壁面最大;并非断面均匀分布。
  2. 湍流断面靠壁处存在黏性底层,仍由黏性切应力主导。
  3. 水力光滑不等于几何表面绝对光滑,它取决于粗糙凸起与黏性底层厚度的相对关系。
  4. 水力粗糙区 λ\lambdaReRe 无关,但水头损失仍随 v2v^2 增大。
  5. λ=64/Re\lambda=64/Re 只适用于充分发展圆管层流。
  6. 突扩的局部阻力系数随选用 v1v_1v2v_2 而变化,实际 hjh_j 不变。
  7. 出口损失 ζ=1\zeta=1 表示速度水头在大水池中耗散。
  8. 边界层分离的关键条件是逆压梯度与近壁动量不足;分离点壁面切应力为零。
  9. Stokes 沉降速度算完后必须校核 Re<1Re<1


11 最后速记#

11.1 三个总公式#

hw=hf+hj\boxed{h_w=\sum h_f+\sum h_j}hf=λLdv22g\boxed{h_f=\lambda\frac Ld\frac{v^2}{2g}}hj=ζv22g\boxed{h_j=\zeta\frac{v^2}{2g}}

11.2 均匀流与圆管切应力#

τ0=ρgRJ\boxed{\tau_0=\rho gRJ}ττ0=rr0\boxed{\frac\tau{\tau_0}=\frac r{r_0}}u=τ0ρ=vλ8\boxed{u_*=\sqrt{\frac{\tau_0}{\rho}}=v\sqrt{\frac\lambda8}}

11.3 圆管层流#

Re<2300\boxed{Re<2300}u=2v[1(rr0)2]\boxed{u=2v\left[1-\left(\frac r{r_0}\right)^2\right]}umax=2v\boxed{u_{\max}=2v}λ=64Re\boxed{\lambda=\frac{64}{Re}}Q=πd4128μΔpL\boxed{Q=\frac{\pi d^4}{128\mu}\frac{\Delta p}{L}}

11.4 湍流#

ui=uˉi+ui\boxed{u_i=\bar u_i+u_i'}τ=μduˉdyρuv\boxed{\tau=\mu\frac{d\bar u}{dy}-\rho\overline{u'v'}}δl=11.6νu\boxed{\delta_l=\frac{11.6\nu}{u_*}}1λ=2lg(2.51Reλ+Δ3.7d)\boxed{ \frac1{\sqrt\lambda} =-2\lg\left( \frac{2.51}{Re\sqrt\lambda}+\frac\Delta{3.7d} \right) }

11.5 非圆管与明渠#

de=4R=4Aχ\boxed{d_e=4R=\frac{4A}{\chi}}v=CRJ\boxed{v=C\sqrt{RJ}}v=1nR2/3J1/2\boxed{v=\frac1nR^{2/3}J^{1/2}}

11.6 局部损失#

突扩:

hj=(v1v2)22g\boxed{h_j=\frac{(v_1-v_2)^2}{2g}}

出口:

ζ=1\boxed{\zeta=1}

11.7 边界层与绕流#

Rex=Uxν\boxed{Re_x=\frac{U_\infty x}{\nu}}

分离点:

τ0=0\boxed{\tau_0=0}

阻力:

FD=CDρu022A\boxed{F_D=C_D\frac{\rho u_0^2}{2}A}

Stokes 沉降:

u0=d2g(ρsρ)18μ,Re<1\boxed{u_0=\frac{d^2g(\rho_s-\rho)}{18\mu},\quad Re<1}

资料对应说明#

  • 《English_Flow resistance-fundamentals》:Section 1–3,阻力概述、流态判别、圆管层流;
  • 《English_Flow resistance-Turbulence》:Section 4–6,湍流理论、圆管湍流阻力、非圆管与明渠;
  • 《English_Flow resistance-local loss》:Section 7–8,局部损失与边界层;
  • 《English_Flow resistance-Cd》:Section 9 与全章总结,绕流阻力和球体沉降;
  • 教材第 6 章:用于补充推导、适用条件与习题;
  • 四份历年卷:用于筛选 Reynolds 数、圆管切应力、黏性底层、粗糙区、层流流量和突扩损失等高频考点。

第二部分:练习题#

本部分按照“概念判断 → 管内流动 → 能量损失 → 边界层与绕流”的顺序安排,共 16 题

题目来源标记:

  • 【24–25 真题】:直接整理自 2024–2025 春夏期末卷;
  • 【历年卷改编】:依据 2021–2024 历年卷中的同类考点整理;
  • 【课件例题】:来自本章四份课堂课件;
  • 【教材习题】:来自教材第 6 章习题;
  • 【综合变式】:围绕课堂公式设计的少量综合题。
TIP

建议先合上答案独立作答:

  • 概念题:每题 2–4 分钟;
  • 单公式计算题:每题 6–8 分钟;
  • 综合损失题:每题 10–15 分钟。

A 基础概念与判断#

A1 Reynolds number 的物理意义#

【24–25 真题】

English

The Reynolds number represents the ratio of ______ force to ______ force in a fluid flow.

中文

雷诺数表示流动中 ______ 力与 ______ 力之比。

答案与讲解inertial force / viscous force\boxed{\text{inertial force / viscous force}}

即:

Re惯性力黏性力\boxed{Re\sim\frac{\text{惯性力}}{\text{黏性力}}}

ReRe 小,黏性作用相对显著,流动趋于有序;ReRe 大,惯性作用强,扰动更容易发展为湍流。

A2 管径变化与层流流量#

【24–25 真题与 23–24 真题改编】

English

A fully developed laminar flow occurs in a circular pipe. The pipe length and the pressure difference between its ends remain unchanged. If the pipe diameter is doubled, the flow rate becomes eight times the original value. True or false?

中文

圆管内为充分发展层流,管长和两端压强差保持不变。若管径增大为原来的 2 倍,流量变为原来的 8 倍。判断正误。

答案与讲解False\boxed{\text{False}}

Hagen–Poiseuille 公式:

Q=πd4128μΔpLQ=\frac{\pi d^4}{128\mu}\frac{\Delta p}{L}

其余条件不变时:

Qd4Q\propto d^4

所以:

Q2Q1=(d2d1)4=24=16\frac{Q_2}{Q_1}=\left(\frac{d_2}{d_1}\right)^4=2^4=16

正确结论为流量增大到原来的 16 倍

A3 圆管切应力与粗糙区阻力#

【24–25 真题】

English

Judge the following statements.

  1. The shear stress is uniformly distributed over the cross-section of a uniform circular-pipe flow.
  2. In the fully rough turbulent regime, the friction factor is independent of Reynolds number.

中文

判断下列说法:

  1. 圆管均匀流过流断面上的切应力均匀分布。
  2. 在水力粗糙湍流区,沿程阻力系数与雷诺数无关。
答案与讲解
  1. False\boxed{\text{False}}

    圆管均匀流:

    ττ0=rr0\frac\tau{\tau_0}=\frac r{r_0}

    切应力从中心的 0 线性增加到壁面的最大值 τ0\tau_0

  2. True\boxed{\text{True}}

    水力粗糙区中粗糙凸起穿出黏性底层,阻力主要由粗糙形状决定:

    λ=f(Δd)\lambda=f\left(\frac\Delta d\right)

    因而与 ReRe 无关。

A4 黏性底层厚度#

【24–25 真题】

English

For water flowing in a pipe, the pipe diameter, water temperature and friction factor remain constant. As the flow rate increases, the thickness of the viscous sublayer will:

A. increase
B. decrease
C. remain unchanged
D. be uncertain

中文

管径、水温和沿程阻力系数保持不变。随流量增大,黏性底层厚度将:

A. 增大
B. 减小
C. 不变
D. 无法判断

答案与讲解B. decrease\boxed{\text{B. decrease}}δl=32.8dReλ\delta_l=\frac{32.8d}{Re\sqrt\lambda}

d,ν,λd,\nu,\lambda 不变时,流量增大意味着 vv 增大,从而 Re=vd/νRe=vd/\nu 增大,所以 δl\delta_l 变小。

A5 为什么使用下临界雷诺数#

【课件思考题】

English

Why is the lower critical Reynolds number, rather than the upper critical Reynolds number, used as the criterion for distinguishing laminar and turbulent flows?

中文

为什么判别层流和湍流时使用下临界雷诺数,而不使用上临界雷诺数?

答案与讲解

上临界雷诺数对应升速过程中层流失稳,其数值对入口扰动、设备振动和来流脉动很敏感,不稳定。

下临界雷诺数对应降速过程中湍流恢复为层流,其数值较稳定,主要由过流断面形状决定。因此工程判别采用下临界值。

圆管中通常取:

Rec2300Re_c\approx2300

A6 圆管湍流的速度与切应力分布#

【21–22 真题改编】

English

Describe the velocity distribution and shear-stress distribution of a fully developed turbulent flow in a circular pipe, and sketch them qualitatively.

中文

说明圆管充分发展湍流的速度分布和切应力分布特点,并定性画图。

答案与讲解

速度分布:

  • 管壁无滑移,u=0u=0
  • 黏性底层中近似线性分布;
  • 湍流核心区近似满足对数律或幂律;
  • 断面速度较层流均匀,中心速度约为 1.15v1.26v1.15v\sim1.26v

总切应力分布:

ττ0=rr0\frac\tau{\tau_0}=\frac r{r_0}
  • 管轴处为 0;
  • 向壁面线性增大;
  • 管壁处达到 τ0\tau_0

组成变化:

  • 近壁区:黏性切应力占主导;
  • 核心区:Reynolds 湍流切应力占主导。

作图时,速度曲线应比层流抛物线更“平坦”;切应力画成从中心 0 到壁面最大值的直线。


B 圆管层流与湍流计算#

B1 通风管流态与最大层流流量#

【教材习题 6-2】

English

Air flows through a ventilation pipe of diameter d=0.25 md=0.25\ \mathrm m at an average velocity of v=3 m/sv=3\ \mathrm{m/s}. At 30C30^\circ\mathrm C, the kinematic viscosity is ν=16.07×106 m2/s\nu=16.07\times10^{-6}\ \mathrm{m^2/s}.

  1. Determine the flow regime.
  2. Determine the maximum flow rate for which laminar flow can be maintained, using Rec=2300Re_c=2300.

中文

通风管直径 d=0.25 md=0.25\ \mathrm m,平均风速 v=3 m/sv=3\ \mathrm{m/s}30C30^\circ\mathrm C 时运动黏度 ν=16.07×106 m2/s\nu=16.07\times10^{-6}\ \mathrm{m^2/s}

  1. 判断流态;
  2. Rec=2300Re_c=2300,求保持层流的最大流量。
答案与讲解
  1. 雷诺数:
Re=vdν=3×0.2516.07×1064.67×104Re=\frac{vd}{\nu} =\frac{3\times0.25}{16.07\times10^{-6}} \approx4.67\times10^4

因此:

turbulent flow\boxed{\text{turbulent flow}}
  1. 最大层流平均速度:
vmax=Recνd=2300×16.07×1060.250.148 m/sv_{\max}=\frac{Re_c\nu}{d} =\frac{2300\times16.07\times10^{-6}}{0.25} \approx0.148\ \mathrm{m/s}

断面面积:

A=πd24=0.0491 m2A=\frac{\pi d^2}{4}=0.0491\ \mathrm{m^2}

最大流量:

Qmax=Avmax0.0491×0.148=7.26×103 m3/sQ_{\max}=Av_{\max} \approx0.0491\times0.148 =7.26\times10^{-3}\ \mathrm{m^3/s}Qmax7.3×103 m3/s\boxed{Q_{\max}\approx7.3\times10^{-3}\ \mathrm{m^3/s}}

B2 粗糙圆管沿程损失#

【课件例题】

English

A pipe has length L=500 mL=500\ \mathrm m, diameter d=0.20 md=0.20\ \mathrm m, roughness height Δ=0.10 mm\Delta=0.10\ \mathrm{mm} and carries water at 10C10^\circ\mathrm C at a flow rate Q=0.010 m3/sQ=0.010\ \mathrm{m^3/s}. Take ν=1.31×106 m2/s\nu=1.31\times10^{-6}\ \mathrm{m^2/s}. Use the Barr equation to determine the frictional head loss.

中文

管长 500 m500\ \mathrm m、直径 0.20 m0.20\ \mathrm m、粗糙高度 0.10 mm0.10\ \mathrm{mm},输送 10C10^\circ\mathrm C 水,流量 0.010 m3/s0.010\ \mathrm{m^3/s},取 ν=1.31×106 m2/s\nu=1.31\times10^{-6}\ \mathrm{m^2/s}。用 Barr 公式求沿程水头损失。

答案与讲解

平均速度:

v=4Qπd2=0.318 m/sv=\frac{4Q}{\pi d^2}=0.318\ \mathrm{m/s}

雷诺数:

Re=vdν4.86×104Re=\frac{vd}{\nu} \approx4.86\times10^4

相对粗糙度:

Δd=0.00010.2=5.0×104\frac\Delta d=\frac{0.0001}{0.2}=5.0\times10^{-4}

Barr 公式:

1λ=2lg(5.1286Re0.89+Δ3.7d)\frac1{\sqrt\lambda} =-2\lg\left( \frac{5.1286}{Re^{0.89}}+\frac\Delta{3.7d} \right)

解得:

λ0.0227\lambda\approx0.0227

沿程损失:

hf=0.02275000.20.31822g0.293 mh_f=0.0227\frac{500}{0.2}\frac{0.318^2}{2g} \approx0.293\ \mathrm mhf0.29 m\boxed{h_f\approx0.29\ \mathrm m}

B3 均匀圆管流的切应力#

【教材习题 6-6】

English

A uniform pipe flow has L=100 mL=100\ \mathrm m, d=0.020 md=0.020\ \mathrm m and hydraulic gradient J=0.008J=0.008. For water with ρ=1000 kg/m3\rho=1000\ \mathrm{kg/m^3}, determine:

  1. the wall shear stress;
  2. the shear stress at r=0.005 mr=0.005\ \mathrm m;
  3. the head loss.

中文

均匀圆管流长 100 m100\ \mathrm m、直径 0.020 m0.020\ \mathrm m、水力坡度 J=0.008J=0.008。取水密度 1000 kg/m31000\ \mathrm{kg/m^3},求壁面切应力、r=0.005 mr=0.005\ \mathrm m 处切应力和水头损失。

答案与讲解

圆管水力半径:

R=d4=0.005 mR=\frac d4=0.005\ \mathrm m
  1. 壁面切应力:
τ0=ρgRJ=1000×9.81×0.005×0.008\tau_0=\rho gRJ =1000\times9.81\times0.005\times0.008τ0=0.392 Pa\boxed{\tau_0=0.392\ \mathrm{Pa}}
  1. 半径 rr 处的同心流管水力半径为 R=r/2R'=r/2
τ=ρgr2J\tau=\rho g\frac r2Jτ=1000×9.81×0.0052×0.008\tau=1000\times9.81\times\frac{0.005}{2}\times0.008τ=0.196 Pa\boxed{\tau=0.196\ \mathrm{Pa}}

也可用:

τ=τ0rr0=0.392×0.0050.010=0.196 Pa\tau=\tau_0\frac r{r_0} =0.392\times\frac{0.005}{0.010} =0.196\ \mathrm{Pa}
  1. 水头损失:
hf=JL=0.008×100=0.80 m\boxed{h_f=JL=0.008\times100=0.80\ \mathrm m}

C 沿程损失、局部损失与明渠#

C1 推导突然扩大损失#

【21–22 真题改编】

English

For a horizontal pipe with a sudden enlargement from area A1A_1 to A2A_2, derive the head-loss formula

hj=(v1v2)22g.h_j=\frac{(v_1-v_2)^2}{2g}.

Assume steady incompressible flow and take α1=α2=β1=β2=1\alpha_1=\alpha_2=\beta_1=\beta_2=1.

中文

水平管道由 A1A_1 突然扩大到 A2A_2。假设恒定不可压缩流,且 α1=α2=β1=β2=1\alpha_1=\alpha_2=\beta_1=\beta_2=1,推导:

hj=(v1v2)22g.h_j=\frac{(v_1-v_2)^2}{2g}.
答案与讲解

连续方程:

A1v1=A2v2=QA_1v_1=A_2v_2=Q

取突扩区域为控制体。除入口断面上的压力 p1A1p_1A_1 外,突扩台阶环形面上还受近似为 p1p_1 的压力:

p1(A2A1)p_1(A_2-A_1)

因此总的上游压力作用为:

p1A1+p1(A2A1)=p1A2p_1A_1+p_1(A_2-A_1)=p_1A_2

沿流向动量方程:

(p1p2)A2=ρQ(v2v1)(p_1-p_2)A_2=\rho Q(v_2-v_1)

利用 Q=A2v2Q=A_2v_2

p1p2ρ=v2(v2v1)\frac{p_1-p_2}{\rho}=v_2(v_2-v_1)

能量方程:

p1ρg+v122g=p2ρg+v222g+hj\frac{p_1}{\rho g}+\frac{v_1^2}{2g} = \frac{p_2}{\rho g}+\frac{v_2^2}{2g}+h_j

将动量方程给出的压强关系代入,化简得:

hj=(v1v2)22g\boxed{h_j=\frac{(v_1-v_2)^2}{2g}}

进一步利用 v2=v1A1/A2v_2=v_1A_1/A_2

hj=(1A1A2)2v122g\boxed{ h_j=\left(1-\frac{A_1}{A_2}\right)^2\frac{v_1^2}{2g} }

**得分关键:**必须同时列出连续、动量、能量三个方程,并说明速度与动量修正系数取 1。

C2 两级突扩的最小损失#

【课件例题】

English

A flow slows down from velocity v1v_1 to v2v_2 through a sudden enlargement. An intermediate pipe is inserted, producing two successive sudden enlargements. Neglect interaction between the two local regions.

  1. Determine the intermediate velocity vv that minimizes the total head loss.
  2. Compare the minimum loss with that of a single sudden enlargement.

中文

流速由 v1v_1 经突扩降至 v2v_2。在中间加入一段管道,使过程变为两次连续突扩,忽略两局部区相互影响。求使总损失最小的中间速度,并与单次突扩比较。

答案与讲解

两次突扩总损失:

hj=(v1v)22g+(vv2)22gh_j=\frac{(v_1-v)^2}{2g}+\frac{(v-v_2)^2}{2g}

求极小值:

dhjdv=2(v1v)+2(vv2)2g=0\frac{dh_j}{dv} =\frac{-2(v_1-v)+2(v-v_2)}{2g}=0

因此:

v=v1+v22\boxed{v=\frac{v_1+v_2}{2}}

代回:

hj,min=2×[(v1v2)/2]22g=(v1v2)24gh_{j,\min} =2\times\frac{[(v_1-v_2)/2]^2}{2g} =\frac{(v_1-v_2)^2}{4g}

单次突扩:

hj,1=(v1v2)22gh_{j,1}=\frac{(v_1-v_2)^2}{2g}

所以:

hj,min=12hj,1\boxed{h_{j,\min}=\frac12h_{j,1}}

本质:将总速度差均匀分配到两次突扩,可减弱每次混合损失。

C3 梯形明渠均匀流#

【教材习题 6-15】

English

A trapezoidal earth channel has bottom width b=2 mb=2\ \mathrm m, side slope m=1.5m=1.5 (horizontal to vertical), water depth h=1.5 mh=1.5\ \mathrm m, bed slope J=0.0004J=0.0004 and Manning coefficient n=0.0225n=0.0225. Determine the mean velocity and discharge.

中文

梯形土渠底宽 b=2 mb=2\ \mathrm m,边坡系数 m=1.5m=1.5,水深 h=1.5 mh=1.5\ \mathrm m,底坡 J=0.0004J=0.0004,曼宁糙率 n=0.0225n=0.0225。求平均流速与流量。

答案与讲解

梯形断面面积:

A=h(b+mh)=1.5(2+1.5×1.5)=6.375 m2A=h(b+mh) =1.5(2+1.5\times1.5) =6.375\ \mathrm{m^2}

湿周:

χ=b+2h1+m2\chi=b+2h\sqrt{1+m^2}χ=2+2×1.51+1.527.408 m\chi=2+2\times1.5\sqrt{1+1.5^2} \approx7.408\ \mathrm m

水力半径:

R=Aχ0.861 mR=\frac A\chi \approx0.861\ \mathrm m

Manning 公式:

v=1nR2/3J1/2v=\frac1nR^{2/3}J^{1/2}v=10.0225(0.861)2/3(0.0004)1/20.804 m/sv=\frac1{0.0225}(0.861)^{2/3}(0.0004)^{1/2} \approx0.804\ \mathrm{m/s}v0.80 m/s\boxed{v\approx0.80\ \mathrm{m/s}}

流量:

Q=Av=6.375×0.8045.13 m3/sQ=Av=6.375\times0.804 \approx5.13\ \mathrm{m^3/s}Q5.13 m3/s\boxed{Q\approx5.13\ \mathrm{m^3/s}}

C4 引水管所需作用水头#

【教材习题 6-23 改编】

English

Water is discharged from a large reservoir through a pipe of diameter d=0.50 md=0.50\ \mathrm m and total length L=1200 mL=1200\ \mathrm m. The flow rate is Q=0.20 m3/sQ=0.20\ \mathrm{m^3/s} and the Darcy friction factor is λ=0.015\lambda=0.015. The system contains a sharp-edged entrance with ζe=0.5\zeta_e=0.5, three bends each with ζb=0.8\zeta_b=0.8, and a valve with ζv=0.2\zeta_v=0.2. The outlet is a free jet. Determine the required head HH between the reservoir surface and the outlet centreline.

中文

大水池通过直径 0.50 m0.50\ \mathrm m、总长 1200 m1200\ \mathrm m 的管道自由出流。流量 0.20 m3/s0.20\ \mathrm{m^3/s},Darcy 阻力系数 λ=0.015\lambda=0.015。系统有直角进口 ζe=0.5\zeta_e=0.5、3 个弯头且每个 ζb=0.8\zeta_b=0.8、阀门 ζv=0.2\zeta_v=0.2。求水池自由液面到出口中心的所需作用水头 HH

答案与讲解

管内平均速度:

A=πd24=0.1963 m2A=\frac{\pi d^2}{4}=0.1963\ \mathrm{m^2}v=QA=0.200.19631.019 m/sv=\frac QA=\frac{0.20}{0.1963} \approx1.019\ \mathrm{m/s}

速度水头:

v22g0.0529 m\frac{v^2}{2g}\approx0.0529\ \mathrm m

从自由液面到自由射流出口列 Bernoulli 方程:

H=v22g+hf+hjH=\frac{v^2}{2g}+h_f+\sum h_j

沿程损失:

hf=λLdv22gh_f=\lambda\frac Ld\frac{v^2}{2g}

局部损失系数总和:

ζ=0.5+3×0.8+0.2=3.1\sum\zeta=0.5+3\times0.8+0.2=3.1

因此:

H=[1+λLd+ζ]v22gH=\left[1+\lambda\frac Ld+\sum\zeta\right]\frac{v^2}{2g}H=[1+0.01512000.5+3.1]0.0529H=\left[1+0.015\frac{1200}{0.5}+3.1\right]0.0529H=(40.1)(0.0529)2.12 mH=(40.1)(0.0529)\approx2.12\ \mathrm mH2.12 m\boxed{H\approx2.12\ \mathrm m}

**易错点:**出口为自由射流,因此出口速度水头 v2/(2g)v^2/(2g) 要保留;只有排入大水池时才可用出口局部损失 ζ=1\zeta=1 代替该速度耗散过程。


D 边界层与绕流#

D1 光滑平板边界层厚度#

【综合变式】

English

Water with ν=1.0×106 m2/s\nu=1.0\times10^{-6}\ \mathrm{m^2/s} flows over a smooth flat plate at U=1.0 m/sU_\infty=1.0\ \mathrm{m/s}.

  1. At x=0.10 mx=0.10\ \mathrm m, assume a laminar boundary layer and estimate its thickness.
  2. At x=5.0 mx=5.0\ \mathrm m, assume a turbulent boundary layer and estimate its thickness.

Use

δl=5xRex,δt=0.377xRex1/5.\delta_l=\frac{5x}{\sqrt{Re_x}},\qquad \delta_t=\frac{0.377x}{Re_x^{1/5}}.

中文

运动黏度 1.0×106 m2/s1.0\times10^{-6}\ \mathrm{m^2/s} 的水以 1.0 m/s1.0\ \mathrm{m/s} 流过光滑平板。分别估算 x=0.10 mx=0.10\ \mathrm m 处的层流边界层厚度和 x=5.0 mx=5.0\ \mathrm m 处的湍流边界层厚度。

答案与讲解
  1. x=0.10 mx=0.10\ \mathrm m
Rex=Uxν=1.0×0.101.0×106=1.0×105Re_x=\frac{U_\infty x}{\nu} =\frac{1.0\times0.10}{1.0\times10^{-6}} =1.0\times10^5δl=5×0.10105=1.58×103 m\delta_l=\frac{5\times0.10}{\sqrt{10^5}} =1.58\times10^{-3}\ \mathrm mδl1.58 mm\boxed{\delta_l\approx1.58\ \mathrm{mm}}
  1. x=5.0 mx=5.0\ \mathrm m
Rex=1.0×5.01.0×106=5.0×106Re_x=\frac{1.0\times5.0}{1.0\times10^{-6}} =5.0\times10^6δt=0.377×5.0(5.0×106)1/58.62×102 m\delta_t=\frac{0.377\times5.0}{(5.0\times10^6)^{1/5}} \approx8.62\times10^{-2}\ \mathrm mδt0.086 m\boxed{\delta_t\approx0.086\ \mathrm m}

结果反映边界层沿流向不断增厚。

D2 边界层分离#

【课件思考题】

English

Explain how boundary-layer separation is formed. State the condition at the separation point and give one effective method for reducing the wake.

中文

说明边界层分离如何形成,写出分离点条件,并给出一种减小尾流区的方法。

答案与讲解

当流体沿流向进入升压区:

dpdx>0\frac{dp}{dx}>0

外部主流减速。边界层内近壁流体原本动量较小,又持续受到黏性阻力,可能先停滞,再出现反向流动。边界流线随之离开壁面,形成分离与尾流。

分离点:

τ0=μ(uy)y=0=0\boxed{ \tau_0=\mu\left(\frac{\partial u}{\partial y}\right)_{y=0}=0 }

减小尾流的方法:将物体设计为流线型,使断面变化平缓,减弱逆压梯度并推迟分离。

D3 烟气中的颗粒能否沉降#

【教材习题 6-24】

English

Flue gas rises at 0.60 m/s0.60\ \mathrm{m/s}. Its kinematic viscosity is 2.23×104 m2/s2.23\times10^{-4}\ \mathrm{m^2/s} and density is 0.20 kg/m30.20\ \mathrm{kg/m^3}. A spherical particle has diameter d=0.10 mmd=0.10\ \mathrm{mm} and density ρs=1.1×103 kg/m3\rho_s=1.1\times10^3\ \mathrm{kg/m^3}. Determine whether the particle can settle downward. Use Stokes’ law and verify its validity.

中文

烟气以 0.60 m/s0.60\ \mathrm{m/s} 向上运动,运动黏度 2.23×104 m2/s2.23\times10^{-4}\ \mathrm{m^2/s}、密度 0.20 kg/m30.20\ \mathrm{kg/m^3}。球形颗粒直径 0.10 mm0.10\ \mathrm{mm}、密度 1.1×103 kg/m31.1\times10^3\ \mathrm{kg/m^3}。用 Stokes 公式判断颗粒能否向下沉降,并校核适用性。

答案与讲解

动力黏度:

μ=ρν=0.20×2.23×104=4.46×105 Pas\mu=\rho\nu =0.20\times2.23\times10^{-4} =4.46\times10^{-5}\ \mathrm{Pa\cdot s}

颗粒相对于烟气的终端沉降速度:

u0=d2g(ρsρ)18μu_0=\frac{d^2g(\rho_s-\rho)}{18\mu}

代入 d=1.0×104 md=1.0\times10^{-4}\ \mathrm m

u00.134 m/su_0\approx0.134\ \mathrm{m/s}

校核:

Re=u0dν=0.134×1042.23×1040.060<1Re=\frac{u_0d}{\nu} =\frac{0.134\times10^{-4}}{2.23\times10^{-4}} \approx0.060<1

Stokes 公式适用。

烟气向上速度为 0.60 m/s0.60\ \mathrm{m/s},大于颗粒相对向下沉降速度 0.134 m/s0.134\ \mathrm{m/s}。以地面为参考,颗粒速度为:

0.600.134=0.466 m/s向上0.60-0.134=0.466\ \mathrm{m/s}\quad\text{向上}

因此:

颗粒不能向下沉降,将被烟气向上携带。\boxed{\text{颗粒不能向下沉降,将被烟气向上携带。}}

FluidMechanics-Chapter6:流体阻力及能量损失
https://www.sleepyfish2031.top/posts/课程笔记/流体力学/chapter6/
作者
Sleepyfish
发布于
2026-06-16
许可协议
CC BY-NC-SA 4.0