第18章：传热与流动耦合

本章大纲

1. 开篇与导言 - 传热与流动的相互作用 - 工程应用背景

2. 强制对流换热 - 物理机制与特征 - 平板、圆管、绕流体的换热 - 雷诺类比与普朗特类比

3. 自然对流与浮力效应 - 浮力驱动流动 - 瑞利数与格拉晓夫数 - 垂直平板与水平平板的自然对流

4. 热边界层理论 - 速度边界层与温度边界层 - 普朗特数的物理意义 - 热边界层厚度估算

5. 努塞尔数关联式 - 无量纲传热系数 - 典型几何的经验关联式 - 工程设计中的应用

6. 历史人物：让-巴蒂斯特·傅里叶 - 傅里叶与热传导理论 - 从军事工程师到数学物理学家

7. 高级话题：共轭传热与热管理优化 - 流固耦合传热 - 现代热管理技术

1. 开篇与导言

当你用手触摸刚从冰箱拿出的易拉罐时，为什么会感觉到冰凉？当你站在电风扇前，为什么会感觉凉爽？这些日常现象背后，都涉及到传热与流动的耦合作用。在自然界和工程应用中，热量传递与流体运动几乎总是相伴而生，它们相互影响、相互制约，形成了复杂而有趣的物理现象。

想象一下笔记本电脑的散热系统：热量从CPU芯片产生，通过导热管传递到散热片，再由风扇带动的气流将热量带走。这个过程中，固体导热、流体对流、甚至辐射传热都在同时发生。理解这些机制，对于设计高效的热管理系统至关重要。

本章将探讨传热与流动耦合的基本规律，重点关注对流换热现象。我们将从物理直觉出发，建立简单而实用的工程估算方法，让你能够快速评估各种传热问题。

2. 强制对流换热

2.1 物理机制

强制对流换热是指流体在外力（如泵、风扇）驱动下流过固体表面时发生的热交换。这种换热方式在工程中极为常见：汽车散热器、空调换热器、电子设备冷却等都属于强制对流。

关键物理过程：

• 热量从壁面传递到流体：通过分子扩散（导热）
• 热量在流体中输运：通过对流（宏观流动）
• 边界层内的温度变化：从壁面温度过渡到主流温度

壁面 ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
      ↑ 热流 q"
    ╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱╱  ← 热边界层
  ─────────────────────────────→  主流 V∞, T∞

2.2 平板强制对流

考虑均匀来流掠过恒温平板的情况。这是最基本的强制对流问题，许多复杂几何都可以简化为平板问题。

层流边界层（Re_x < 5×10^5）：

• 局部传热系数：$h_x = 0.332 \frac{k}{x} Re_x^{1/2} Pr^{1/3}$
• 平均传热系数：$\bar{h} = 0.664 \frac{k}{L} Re_L^{1/2} Pr^{1/3}$

湍流边界层（Re_x > 5×10^5）：

• 局部传热系数：$h_x = 0.0296 \frac{k}{x} Re_x^{4/5} Pr^{1/3}$

工程估算技巧：

• 空气在常温下：Pr ≈ 0.7
• 水在常温下：Pr ≈ 7
• 传热系数量级：空气 10-100 W/(m²·K)，水 500-10000 W/(m²·K)

2.3 圆管内强制对流

圆管内的流动换热是工程中最常见的情况之一。换热器、管道系统都涉及这类问题。

充分发展流动：

• 层流（Re_D < 2300）：Nu_D = 3.66（恒温壁面）或 Nu_D = 4.36（恒热流）
• 湍流（Re_D > 10^4）：Nu_D = 0.023 Re_D^{4/5} Pr^n
• 加热流体：n = 0.4
• 冷却流体：n = 0.3

入口段效应：

• 热入口长度：L_t ≈ 0.05 Re_D Pr · D（层流）
• 在入口段，传热系数较高（可达充分发展值的2-3倍）

2.4 绕流体强制对流

圆柱横向绕流：

     来流 →  ○  ← 圆柱
           涡街

• Nu_D = C Re_D^m Pr^{1/3}
• C和m取决于Re_D范围
• 驻点处传热最强，分离点附近传热最弱

球体绕流：

• Nu_D = 2 + (0.4 Re_D^{1/2} + 0.06 Re_D^{2/3}) Pr^{0.4}
• 适用范围：0.71 < Pr < 380，3.5 < Re_D < 7.6×10^4

2.5 雷诺类比

雷诺类比揭示了动量传递与热量传递的相似性：

$$St = \frac{C_f}{2}$$ 其中St是斯坦顿数（$St = Nu/(Re·Pr)$），$C_f$是摩擦系数。

物理意义：

• 动量传递越快（摩擦系数大），热量传递也越快
• 可以通过测量阻力来推算传热系数
• 对光滑表面、Pr≈1的流体最准确

修正的雷诺类比（普朗特-泰勒类比）： $$St = \frac{C_f/2}{1 + 5\sqrt{C_f/2}(Pr - 1)}$$ 适用于0.6 < Pr < 60的范围。

3. 自然对流与浮力效应

3.1 浮力驱动的流动

自然对流是由温度差引起的密度差驱动的流动。没有外部强制，流体完全靠浮力自发运动。这种现象在日常生活中随处可见：

• 室内暖气的热循环：热空气上升，冷空气下沉
• 海陆风：白天陆地升温快，形成海风；夜晚陆地降温快，形成陆风
• 烟囱效应：热烟气的浮力产生抽吸作用
• 冰箱内的对流：冷空气下沉，形成自然循环

3.2 控制参数

格拉晓夫数（Grashof number）： $$Gr = \frac{g\beta\Delta T L^3}{\nu^2}$$ 物理意义：浮力与粘性力之比的平方

• g：重力加速度
• β：体积膨胀系数（理想气体：β = 1/T）
• ΔT：温差
• L：特征长度
• ν：运动粘度

瑞利数（Rayleigh number）： $$Ra = Gr \cdot Pr = \frac{g\beta\Delta T L^3}{\nu\alpha}$$ 物理意义：浮力驱动与热扩散的竞争

• Ra < 10^9：层流自然对流
• Ra > 10^9：湍流自然对流

3.3 垂直平板自然对流

考虑垂直热板在静止空气中的自然对流：

    ↑ 热羽流
    │
  ╱╱│╱╱  边界层
  ╱╱│╱╱  
热板 │    静止空气
  ╱╱│╱╱
  ╱╱│╱╱
    │

层流区（Ra < 10^9）：

• 局部努塞尔数：$Nu_x = 0.68 + \frac{0.670 Ra_x^{1/4}}{[1 + (0.492/Pr)^{9/16}]^{4/9}}$
• 平均努塞尔数：$\bar{Nu} = 0.68 + \frac{0.670 Ra^{1/4}}{[1 + (0.492/Pr)^{9/16}]^{4/9}}$

湍流区（Ra > 10^9）：

• $Nu = 0.825 + \frac{0.387 Ra^{1/6}}{[1 + (0.492/Pr)^{9/16}]^{8/27}}$

工程简化（空气，Pr ≈ 0.7）：

• 层流：Nu ≈ 0.59 Ra^{1/4}
• 湍流：Nu ≈ 0.10 Ra^{1/3}

3.4 水平平板自然对流

热面朝上（不稳定）：

  ↑↑↑↑↑ 热羽流（贝纳德对流胞）
━━━━━━━━━━━━
  热板

• Ra < 10^7：Nu = 0.54 Ra^{1/4}
• Ra > 10^7：Nu = 0.15 Ra^{1/3}

热面朝下（稳定）：

  热板
━━━━━━━━━━━━
  稳定分层

• Nu = 0.27 Ra^{1/4}（仅适用于10^5 < Ra < 10^10）

3.5 封闭空间自然对流

垂直夹层（两平行板间）：

热│  ←循环  │冷
板│    ↓    │板
  │    ↑    │

纵横比影响很大：

• 窄缝（H/L >> 1）：单胞环流
• 方腔（H/L ≈ 1）：主环流+角涡
• 宽腔（H/L << 1）：多胞结构

经验关联式（空气）： $$Nu = 0.18 \left(\frac{Pr}{0.2 + Pr} Ra\right)^{0.29} \left(\frac{H}{L}\right)^{-1/4}$$ 适用范围：10^3 < Ra(L/H) < 10^7

3.6 混合对流

当强制对流和自然对流同时存在时：

判据： $$\frac{Gr}{Re^2} = \frac{\text{浮力}}{\text{惯性力}}$$

• Gr/Re² << 1：强制对流主导，忽略浮力
• Gr/Re² ≈ 1：混合对流，两者都重要
• Gr/Re² >> 1：自然对流主导，忽略强制流动

协同与对抗：

• 协同流（浮力与强制流同向）：增强传热
• 对抗流（浮力与强制流反向）：可能出现流动反转
• 横向流（浮力与强制流垂直）：产生二次流

4. 热边界层理论

4.1 速度边界层与温度边界层

当流体流过不同温度的固体表面时，会同时形成速度边界层和温度边界层：

壁面 ━━━━━━━━━━━━━━━━━━━━━━━━━━
     ╱╱╱╱╱╱╱╱╱╱╱╱╱  速度边界层 δ
   ════════════════  温度边界层 δt
 ─────────────────→ 主流

两个边界层的相对厚度取决于普朗特数：

• Pr < 1（液态金属）：δt > δ，热扩散快于动量扩散
• Pr ≈ 1（气体）：δt ≈ δ，两者相当
• Pr > 1（油、水）：δt < δ，动量扩散快于热扩散

4.2 普朗特数的物理意义

$$Pr = \frac{\nu}{\alpha} = \frac{\mu c_p}{k} = \frac{\text{动量扩散率}}{\text{热扩散率}}$$ 典型值：

• 液态金属：Pr ≈ 0.01-0.05（导热极好）
• 气体：Pr ≈ 0.7-1.0（动量与热量传递相当）
• 水：Pr ≈ 1-10（随温度变化）
• 油类：Pr ≈ 50-100,000（导热很差）

工程意义：

• Pr << 1：热边界层很厚，温度影响范围大
• Pr >> 1：热边界层很薄，温度变化集中在壁面附近

4.3 热边界层厚度估算

层流平板： $$\frac{\delta_t}{\delta} \approx Pr^{-1/3}$$ 具体厚度：

• 速度边界层：$\delta \approx 5\sqrt{\frac{\nu x}{V_\infty}}$
• 温度边界层：$\delta_t \approx \delta \cdot Pr^{-1/3}$

圆管入口段：

• 速度入口长度：L_h/D ≈ 0.05 Re_D（层流）
• 热入口长度：L_t/D ≈ 0.05 Re_D Pr（层流）

4.4 积分方法估算

利用积分动量方程和能量方程，可以得到简单的估算：

动量积分方程： $$\frac{d}{dx}\int_0^\delta u(V_\infty - u)dy = \nu\left(\frac{\partial u}{\partial y}\right)_{y=0}$$ 能量积分方程： $$\frac{d}{dx}\int_0^{\delta_t} u(T_\infty - T)dy = \alpha\left(\frac{\partial T}{\partial y}\right)_{y=0}$$ 通过假设速度和温度分布（如多项式分布），可以求解这些方程得到传热系数。

4.5 相似解

对于某些特殊情况，存在相似解：

楔形流（Falkner-Skan流动）：

• 主流速度：$V_\infty = Cx^m$
• m = 0：平板流动（Blasius解）
• m = 1：驻点流动
• m < 0：逆压梯度流动

5. 努塞尔数关联式与工程应用

5.1 努塞尔数的物理意义

$$Nu = \frac{hL}{k} = \frac{\text{对流传热}}{\text{导热}}$$ 努塞尔数表示对流换热相对于纯导热的增强程度：

• Nu = 1：纯导热（静止流体）
• Nu > 1：对流增强传热
• Nu >> 1：强对流，导热可忽略

与其他无量纲数的关系：

• 强制对流：Nu = f(Re, Pr)
• 自然对流：Nu = f(Ra) 或 Nu = f(Gr, Pr)
• 混合对流：Nu = f(Re, Gr, Pr)

5.2 典型几何的经验关联式汇总

5.2.1 外部流动

平板： | 流态 | 关联式 | 适用范围 |

圆柱横流： | Re 范围 | C | m |

通用形式：Nu = C Re^m Pr^{1/3}

球体：

• 白色公式：Nu = 2 + (0.4 Re^{1/2} + 0.06 Re^{2/3}) Pr^{0.4} (μ/μ_s)^{1/4}

5.2.2 内部流动

圆管层流（充分发展）：

• 恒温壁：Nu = 3.66
• 恒热流：Nu = 4.36
• 入口段修正：Nu_x/Nu_fd = 1 + (D/x)^{0.7}

圆管湍流：

• Dittus-Boelter：Nu = 0.023 Re^{0.8} Pr^n
• 加热：n = 0.4
• 冷却：n = 0.3
• Gnielinski（更精确）： $$Nu = \frac{(f/8)(Re-1000)Pr}{1 + 12.7(f/8)^{1/2}(Pr^{2/3}-1)}$$ 其中 f = (0.79 ln Re - 1.64)^{-2}

非圆形管道： 使用水力直径 D_h = 4A/P

• 矩形：D_h = 2ab/(a+b)
• 环形：D_h = D_o - D_i
• 三角形：D_h = 4A/P

5.3 工程应用实例

5.3.1 换热器设计

管壳式换热器：

壳侧流体 →  ═══════════  ← 管侧流体
            ║ ○ ○ ○ ○ ║
            ║ ○ ○ ○ ○ ║
            ═══════════

总传热系数： $$\frac{1}{UA} = \frac{1}{h_i A_i} + \frac{\ln(r_o/r_i)}{2\pi kL} + \frac{1}{h_o A_o} + R_f$$ 其中 R_f 是污垢热阻。

板式换热器：

• 紧凑高效，Nu 可达 200-400
• 典型关联式：Nu = C Re^m Pr^{1/3}
• C = 0.15-0.40，m = 0.65-0.85（取决于板型）

5.3.2 电子冷却

芯片散热器：

     ┌─┬─┬─┬─┬─┐ 散热片
     │ │ │ │ │ │
━━━━━━━━━━━━━━━━━ 基板
     ▓▓▓▓▓▓▓▓▓   芯片

关键参数：

• 肋效率：η = tanh(mL)/(mL)，其中 m = √(hP/kA_c)
• 总热阻：R_total = R_芯片-基板 + R_基板-散热片 + R_对流
• 优化目标：最小化热阻，同时控制压降

热管： 利用相变传热，有效导热系数可达铜的 100-1000 倍

• 蒸发段：q" = h_evap ΔT_sat
• 冷凝段：q" = h_cond ΔT_sat
• 毛细极限：ΔP_cap > ΔP_液 + ΔP_蒸 + ΔP_重力

5.3.3 建筑节能

双层玻璃窗：

室内 │  空气层  │ 室外
     │          │
  T_i │    Ra    │ T_o
     │          │

夹层传热：

• 导热为主（Ra < 10^3）：Nu ≈ 1
• 自然对流（10^3 < Ra < 10^5）：Nu = 0.18 Ra^{0.29}
• 优化间距：使 Ra ≈ 10^3，最小化传热

通风空腔：

• 夏季：利用烟囱效应排热
• 冬季：预热新风
• 设计流速：0.5-2 m/s（平衡传热与压降）

5.4 相变传热简介

5.4.1 沸腾传热

池沸腾曲线：

log q"
  │     ╱CHF
  │   ╱╲
  │ ╱核沸腾
  │╱
  └────────── log ΔT
自然对流│过渡│膜沸腾

临界热流密度（CHF）：

• Zuber关联式：q"_CHF = 0.131 h_fg ρ_v^{1/2} [gσ(ρ_l - ρ_v)]^{1/4}

流动沸腾：

• 强制对流区：Nu = 0.023 Re^{0.8} Pr^{0.4}
• 核沸腾区：q" = f(ΔT_sat)
• 两相流型：泡状流→弹状流→环状流→雾状流

5.4.2 凝结传热

膜状凝结（Nusselt理论）： 垂直平板： $$\bar{h} = 0.943 \left[\frac{g\rho_l(\rho_l - \rho_v)h_{fg}k_l^3}{\mu_l L(T_{sat} - T_w)}\right]^{1/4}$$ 水平管： $$\bar{h} = 0.725 \left[\frac{g\rho_l(\rho_l - \rho_v)h_{fg}k_l^3}{\mu_l D(T_{sat} - T_w)}\right]^{1/4}$$ 滴状凝结：

• 传热系数比膜状凝结高 5-10 倍
• 需要疏水表面
• 工程应用：添加促进剂或表面改性

5.5 强化传热技术

5.5.1 被动强化

扩展表面：

• 肋片：增加传热面积
• 针肋：适用于气体冷却
• 内肋管：增强管内传热

粗糙表面：

• 人工粗糙度：e/D = 0.01-0.05
• 增强湍流混合
• 传热增强 2-3 倍，压降增加 4-9 倍

涡流发生器：

  ↗↘ 纵向涡
━━━△━━━ 三角翼

• 产生二次流
• 减薄边界层
• 延迟分离

5.5.2 主动强化

振动：

• 机械振动：f = 10-100 Hz
• 声学振动：产生声流
• 脉动流：周期性扰动

电场/磁场：

• 电流体动力学（EHD）
• 磁流体动力学（MHD）
• 适用于特殊流体

射流冲击：

• 驻点传热系数最高
• 典型 Nu = 2-3 倍于平行流
• 应用：钢板急冷、电子冷却

6. 历史人物：让-巴蒂斯特·傅里叶

6.1 从军事工程师到数学物理学家

让-巴蒂斯特·约瑟夫·傅里叶（Jean-Baptiste Joseph Fourier, 1768-1830）的一生充满传奇色彩。作为一个裁缝的儿子，他9岁成为孤儿，却最终成为19世纪最伟大的数学物理学家之一。

早年经历：

• 1768年生于法国欧塞尔
• 在本笃会学校接受教育，展现数学天赋
• 21岁成为巴黎高等师范学校首批学生
• 1798年随拿破仑远征埃及，担任科学顾问

埃及岁月的影响： 在埃及的三年（1798-1801）深刻影响了傅里叶的科学生涯。沙漠的极端温差让他对热现象产生浓厚兴趣。他注意到：

• 白天沙漠表面温度可达70°C
• 夜晚温度急剧下降至接近冰点
• 金字塔内部却保持恒温

这些观察激发了他对热传导问题的思考。

6.2 《热的解析理论》

1822年，傅里叶发表了划时代的著作《Théorie Analytique de la Chaleur》（热的解析理论），奠定了传热学的数学基础。

核心贡献：

1. 傅里叶定律： $$q = -k\nabla T$$ 首次将热流与温度梯度联系起来

2. 热传导方程： $$\frac{\partial T}{\partial t} = \alpha \nabla^2 T$$ 描述温度场的时空演化

3. 傅里叶级数： $$f(x) = \frac{a_0}{2} + \sum_{n=1}^{\infty}(a_n\cos(nx) + b_n\sin(nx))$$ 革命性的数学工具，影响远超传热学

物理洞察：

• 热量是一种运动形式（当时革命性观点）
• 温度场的连续性
• 边界条件的重要性
• 相似性原理的雏形

6.3 傅里叶的实验精神

傅里叶不仅是理论家，也是实验家。他设计了多个巧妙的实验：

金属棒实验：

热源 ━━━━━━━━━━━━━━━━━━ 冷端
     T₁  T₂  T₃  T₄  T₅

通过测量不同位置的温度，验证了指数衰减规律。

环形实验： 研究环形金属片的温度分布，发现了周期性边界条件的重要性。

多层材料实验： 研究复合材料的热阻，为现代保温材料设计奠定基础。

6.4 遗产与影响

傅里叶的工作影响深远：

数学领域：

• 傅里叶变换成为信号处理基础
• 谱分析方法
• 偏微分方程理论

物理应用：

• 量子力学（波函数展开）
• 固体物理（晶格振动）
• 光学（衍射理论）

工程应用：

• 传热设备设计
• 建筑节能
• 电子冷却
• 材料加工

名言：

"数学的主要目标应该是公众的利益和自然现象的解释。"

7. 高级话题：共轭传热与热管理优化

7.1 共轭传热问题

共轭传热是指固体导热与流体对流同时求解的问题。在固-液界面上，温度和热流连续：

界面条件：

• 温度连续：$T_{solid} = T_{fluid}$
• 热流连续：$-k_s\frac{\partial T}{\partial n}\bigg|_s = -k_f\frac{\partial T}{\partial n}\bigg|_f$

典型应用：

• 涡轮叶片冷却
• 微电子散热
• 热交换器优化
• 核反应堆安全分析

7.2 现代热管理技术

7.2.1 微通道冷却

 ═══════════════
 ├─┼─┼─┼─┼─┼─┼─┤ 微通道 (D < 1mm)
 ═══════════════
    热源

特点：

• 极高的表面积体积比
• 层流传热系数高
• 压降大，需要优化设计
• 典型传热系数：10^4-10^5 W/(m²·K)

7.2.2 相变材料（PCM）

利用相变潜热进行热缓冲：

工作原理：

• 吸热：固→液，储存能量
• 放热：液→固，释放能量
• 恒温特性：相变期间温度不变

应用实例：

• 建筑墙体（石蜡，T_m ≈ 25°C）
• 电池热管理（共晶盐）
• 航天器热控（正十八烷）

7.2.3 喷雾冷却

  ∘∘∘∘ 雾滴
   ↓↓↓
 ════════ 液膜
 ▓▓▓▓▓▓▓▓ 热表面

传热机制：

• 液滴冲击传热
• 薄膜蒸发
• 核沸腾（高热流时）
• 可达热流密度：> 1000 W/cm²

7.2.4 热电冷却

基于Peltier效应的固态制冷：

优点：

• 无运动部件
• 精确温控（±0.1°C）
• 可逆操作（制冷/加热）

缺点：

• 效率低（COP < 1）
• 成本高
• 适用于小功率场合

7.3 拓扑优化与仿生设计

7.3.1 拓扑优化

通过数学优化确定最佳的材料分布：

目标函数：

• 最小化热阻
• 最大化散热量
• 均匀化温度分布

约束条件：

• 材料用量
• 压降限制
• 制造可行性

优化结果： 树状、分形等仿生结构

7.3.2 仿生热管理

自然界的启发：

叶脉网络：

• 分级结构
• Murray定律：$d_0^3 = d_1^3 + d_2^3$
• 最小化流动阻力

大象耳朵：

• 大面积散热
• 血管网络优化
• 主动控制（扇动）

白蚁巢穴：

• 被动通风系统
• 烟囱效应利用
• 恒温环境维持

7.4 人工智能在热管理中的应用

7.4.1 机器学习预测

传热关联式开发：

• 神经网络拟合复杂几何
• 超越传统经验公式
• 实时性能预测

故障诊断：

• 异常热模式识别
• 预测性维护
• 优化运行策略

7.4.2 强化学习优化

动态热管理：

• 实时调节冷却策略
• 多目标优化（性能vs能耗）
• 自适应控制

案例：数据中心冷却：

• Google使用DeepMind减少40%冷却能耗
• 动态调节风扇转速、冷水温度
• 预测负载变化

7.5 未来展望

新材料：

• 石墨烯（k > 5000 W/m·K）
• 金属泡沫
• 纳米流体

新技术：

• 量子点冷却
• 磁制冷
• 声子工程

系统集成：

• 多尺度耦合
• 跨学科优化
• 智能自适应系统

本章小结

本章探讨了传热与流动耦合的基本规律和工程应用。关键要点：

1. 强制对流换热： - 由外部驱动力（泵、风扇）引起 - Nu = f(Re, Pr) - 平板、圆管、绕流体各有特定关联式 - 雷诺类比揭示动量与热量传递的相似性

2. 自然对流换热： - 由浮力驱动，无需外部动力 - Nu = f(Ra) 或 f(Gr, Pr) - 垂直板、水平板、封闭空间有不同流动形态 - Gr/Re² 判断混合对流的主导机制

3. 热边界层理论： - 速度边界层与温度边界层的相对厚度取决于Pr - Pr < 1：热边界层厚；Pr > 1：热边界层薄 - 边界层厚度估算对工程设计至关重要

4. 努塞尔数关联式： - Nu 表示对流相对于导热的增强程度 - 不同几何和流动条件有相应的经验关联式 - 工程应用：换热器、电子冷却、建筑节能

5. 现代热管理技术： - 微通道、相变材料、喷雾冷却等新技术 - 拓扑优化和仿生设计提供新思路 - AI技术正在革新热管理策略

核心公式汇总：

• 傅里叶定律：$q = -k\nabla T$
• 牛顿冷却定律：$q = h\Delta T$
• 努塞尔数：$Nu = hL/k$
• 雷诺数：$Re = \rho VL/\mu$
• 普朗特数：$Pr = \nu/\alpha = \mu c_p/k$
• 格拉晓夫数：$Gr = g\beta\Delta T L^3/\nu^2$
• 瑞利数：$Ra = Gr \cdot Pr$

练习题

基础题

习题18.1 平板强制对流 一块长度L=0.5m的平板置于20°C、速度3m/s的空气流中，平板表面温度保持80°C。估算： (a) 临界雷诺数对应的位置 (b) 平板末端的局部传热系数 (c) 整板的平均传热系数

提示：空气性质在50°C时查取，Re_cr = 5×10^5

习题18.2 圆管内流动换热 水以0.5 kg/s的流量流过内径25mm、长度3m的圆管，管壁温度恒定为90°C，水的入口温度为20°C。求出口温度。

提示：使用对数平均温差，水的性质：ρ=1000 kg/m³，c_p=4180 J/(kg·K)，k=0.6 W/(m·K)，μ=10^-3 Pa·s

习题18.3 垂直板自然对流 一块高0.3m、宽0.5m的垂直热板在20°C的静止空气中，板温80°C。计算自然对流散热量。

提示：空气在50°C时β=1/T_f，其他性质查表

挑战题

习题18.4 混合对流判断 水平圆管（D=50mm）内有60°C的水以0.1m/s流动，管壁温度20°C。判断是强制对流、自然对流还是混合对流主导？

提示：计算Gr/Re²，考虑重力方向

习题18.5 肋片优化设计 设计CPU散热器的铝制肋片阵列。芯片功率100W，尺寸30×30mm，环境温度25°C，芯片最高温度85°C。空气强制对流h=50 W/(m²·K)。求最优肋片间距和高度。

提示：考虑肋效率和压降平衡

习题18.6 热管设计 设计一根热管传输200W热量，蒸发段和冷凝段长度各100mm，工作温度60°C。选择工作流体并估算热管直径。

提示：考虑毛细极限和携带极限

习题18.7 双层玻璃优化 设计双层玻璃窗以最小化传热。玻璃间距可调（5-20mm），室内20°C，室外-10°C。何时导热主导？何时对流主导？最优间距是多少？

提示：计算不同间距下的Ra数

习题18.8 数据中心冷却策略 某数据中心机柜功率密度10kW/m²，采用地板送风。设计气流组织使芯片温度不超过70°C，环境温度25°C。比较冷/热通道封闭的效果。

提示：考虑气流短路和热回流

常见陷阱与错误 (Gotchas)

1. 混淆局部与平均传热系数 - 错误：直接用局部关联式计算总传热量 - 正确：积分或使用平均值关联式

2. 忽略入口段效应 - 错误：短管用充分发展流动关联式 - 正确：L/D < 10时必须考虑入口段修正

3. Pr数的温度依赖 - 错误：用错误温度下的物性 - 正确：用膜温度或特定参考温度

4. 自然对流的方向性 - 错误：水平板上下表面用相同关联式 - 正确：热面朝上和朝下完全不同

5. 混合对流的错误叠加 - 错误：Nu_total = Nu_forced + Nu_natural - 正确：Nu_total = (Nu_forced^n ± Nu_natural^n)^(1/n)，n≈3

6. 忽略辐射传热 - 错误：高温时只考虑对流 - 正确：T > 100°C时辐射可能占主导

7. 肋片效率误用 - 错误：认为肋片总是增强传热 - 正确：低导热材料或过长肋片可能η < 0.5

8. 相变传热的不稳定性 - 错误：认为沸腾总是增强传热 - 正确：膜态沸腾传热恶化，须避免CHF

最佳实践检查清单

设计阶段

• [ ] 明确传热模式（强制/自然/混合/相变）
• [ ] 计算关键无量纲数（Re, Pr, Gr, Ra）
• [ ] 选择合适的关联式和适用范围
• [ ] 考虑入口效应、发展段长度
• [ ] 评估辐射传热的重要性
• [ ] 检查流动状态（层流/过渡/湍流）

计算验证

• [ ] 物性参数用正确的参考温度
• [ ] 关联式适用范围（Re, Pr范围）
• [ ] 单位一致性检查
• [ ] 量级估算合理性
• [ ] 考虑不确定性（±20-30%）

优化考虑

• [ ] 传热增强vs压降代价权衡
• [ ] 材料成本vs性能提升平衡
• [ ] 可制造性和维护性
• [ ] 污垢和老化的影响
• [ ] 控制策略和运行灵活性

实验验证

• [ ] 测点布置避开入口段
• [ ] 考虑热损失和边缘效应
• [ ] 稳态条件确认
• [ ] 重复性测试
• [ ] 与理论预测对比偏差分析