intuitive_cfd

第18章：传热与流动耦合

本章大纲

  1. 开篇与导言
    • 传热与流动的相互作用
    • 工程应用背景
  2. 强制对流换热
    • 物理机制与特征
    • 平板、圆管、绕流体的换热
    • 雷诺类比与普朗特类比
  3. 自然对流与浮力效应
    • 浮力驱动流动
    • 瑞利数与格拉晓夫数
    • 垂直平板与水平平板的自然对流
  4. 热边界层理论
    • 速度边界层与温度边界层
    • 普朗特数的物理意义
    • 热边界层厚度估算
  5. 努塞尔数关联式
    • 无量纲传热系数
    • 典型几何的经验关联式
    • 工程设计中的应用
  6. 历史人物：让-巴蒂斯特·傅里叶
    • 傅里叶与热传导理论
    • 从军事工程师到数学物理学家
  7. 高级话题：共轭传热与热管理优化
    • 流固耦合传热
    • 现代热管理技术

1. 开篇与导言

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

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

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

2. 强制对流换热

2.1 物理机制

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

关键物理过程：

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

2.2 平板强制对流

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

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

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

工程估算技巧：

2.3 圆管内强制对流

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

充分发展流动：

入口段效应：

2.4 绕流体强制对流

圆柱横向绕流：

     来流 →  ○  ← 圆柱
           涡街

球体绕流：

2.5 雷诺类比

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

\[St = \frac{C_f}{2}\]

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

物理意义：

修正的雷诺类比（普朗特-泰勒类比）： \(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}\)

物理意义：浮力与粘性力之比的平方

瑞利数（Rayleigh number）： \(Ra = Gr \cdot Pr = \frac{g\beta\Delta T L^3}{\nu\alpha}\)

物理意义：浮力驱动与热扩散的竞争

3.3 垂直平板自然对流

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

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

层流区（Ra < 10^9）：

湍流区（Ra > 10^9）：

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

3.4 水平平板自然对流

热面朝上（不稳定）：

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

热面朝下（稳定）：

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

3.5 封闭空间自然对流

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

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

纵横比影响很大：

经验关联式（空气）： \(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{惯性力}}\)

协同与对抗：

4. 热边界层理论

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

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

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

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

4.2 普朗特数的物理意义

\[Pr = \frac{\nu}{\alpha} = \frac{\mu c_p}{k} = \frac{\text{动量扩散率}}{\text{热扩散率}}\]

典型值：

工程意义：

4.3 热边界层厚度估算

层流平板： \(\frac{\delta_t}{\delta} \approx Pr^{-1/3}\)

具体厚度：

圆管入口段：

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流动）：

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

5.1 努塞尔数的物理意义

\[Nu = \frac{hL}{k} = \frac{\text{对流传热}}{\text{导热}}\]

努塞尔数表示对流换热相对于纯导热的增强程度：

与其他无量纲数的关系：

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

5.2.1 外部流动

平板： | 流态 | 关联式 | 适用范围 | |——|——–|———-| | 层流 | Nu = 0.664 Re^{1/2} Pr^{1/3} | Re < 5×10^5, Pr > 0.6 | | 湍流 | Nu = 0.037 Re^{4/5} Pr^{1/3} | 5×10^5 < Re < 10^7 | | 混合 | Nu = (0.037 Re^{4/5} - 871) Pr^{1/3} | Re_cr = 5×10^5 |

圆柱横流： | Re 范围 | C | m | |———|—|—| | 0.4-4 | 0.989 | 0.330 | | 4-40 | 0.911 | 0.385 | | 40-4000 | 0.683 | 0.466 | | 4000-40000 | 0.193 | 0.618 | | 40000-400000 | 0.027 | 0.805 |

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

球体：

5.2.2 内部流动

圆管层流（充分发展）：

圆管湍流：

非圆形管道： 使用水力直径 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 是污垢热阻。

板式换热器：

5.3.2 电子冷却

芯片散热器：

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

关键参数：

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

5.3.3 建筑节能

双层玻璃窗：

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

夹层传热：

通风空腔：

5.4 相变传热简介

5.4.1 沸腾传热

池沸腾曲线：

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

临界热流密度（CHF）：

流动沸腾：

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.5 强化传热技术

5.5.1 被动强化

扩展表面：

粗糙表面：

涡流发生器：

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

5.5.2 主动强化

振动：

电场/磁场：

射流冲击：

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

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

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

早年经历：

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

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

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 共轭传热问题

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

界面条件：

典型应用：

7.2 现代热管理技术

7.2.1 微通道冷却

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

特点：

7.2.2 相变材料（PCM）

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

工作原理：

应用实例：

7.2.3 喷雾冷却

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

传热机制：

7.2.4 热电冷却

基于Peltier效应的固态制冷：

优点：

缺点：

7.3 拓扑优化与仿生设计

7.3.1 拓扑优化

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

目标函数：

约束条件：

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

7.3.2 仿生热管理

自然界的启发：

叶脉网络：

大象耳朵：

白蚁巢穴：

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

7.4.1 机器学习预测

传热关联式开发：

故障诊断：

7.4.2 强化学习优化

动态热管理：

案例：数据中心冷却：

7.5 未来展望

新材料：

新技术：

系统集成：

本章小结

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

  1. 强制对流换热：
    • 由外部驱动力（泵、风扇）引起
    • Nu = f(Re, Pr)
    • 平板、圆管、绕流体各有特定关联式
    • 雷诺类比揭示动量与热量传递的相似性
  2. 自然对流换热：
    • 由浮力驱动，无需外部动力
    • Nu = f(Ra) 或 f(Gr, Pr)
    • 垂直板、水平板、封闭空间有不同流动形态
    • Gr/Re² 判断混合对流的主导机制
  3. 热边界层理论：
    • 速度边界层与温度边界层的相对厚度取决于Pr
    • Pr < 1：热边界层厚；Pr > 1：热边界层薄
    • 边界层厚度估算对工程设计至关重要
  4. 努塞尔数关联式：
    • Nu 表示对流相对于导热的增强程度
    • 不同几何和流动条件有相应的经验关联式
    • 工程应用：换热器、电子冷却、建筑节能
  5. 现代热管理技术：
    • 微通道、相变材料、喷雾冷却等新技术
    • 拓扑优化和仿生设计提供新思路
    • AI技术正在革新热管理策略

核心公式汇总：

练习题

基础题

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

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

答案 取50°C空气性质：ν = 1.8×10^-5 m²/s，k = 0.028 W/(m·K)，Pr = 0.7 (a) 临界位置： Re_cr = V·x_cr/ν = 5×10^5 x_cr = 5×10^5 × 1.8×10^-5 / 3 = 0.3 m (b) 末端Re_L = 3×0.5/(1.8×10^-5) = 8.33×10^4（层流） h_L = 0.332(k/L)Re_L^0.5·Pr^(1/3) = 11.5 W/(m²·K) (c) 平均传热系数： h̄ = 0.664(k/L)Re_L^0.5·Pr^(1/3) = 23.0 W/(m²·K)

习题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

答案 Re = 4ṁ/(πDμ) = 4×0.5/(π×0.025×10^-3) = 2.55×10^4（湍流） Pr = μc_p/k = 10^-3×4180/0.6 = 6.97 Nu = 0.023Re^0.8Pr^0.4 = 0.023×(2.55×10^4)^0.8×6.97^0.4 = 154 h = Nu·k/D = 154×0.6/0.025 = 3696 W/(m²·K) 使用能量平衡： Q = ṁc_p(T_out - T_in) = hA·LMTD 其中A = πDL = π×0.025×3 = 0.236 m² 迭代求解得：T_out ≈ 58°C

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

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

答案 膜温度T_f = (80+20)/2 = 50°C = 323K 空气性质：ν = 1.8×10^-5 m²/s，α = 2.5×10^-5 m²/s，k = 0.028 W/(m·K) β = 1/323 = 3.1×10^-3 K^-1 Ra = gβΔTL³/(να) = 9.8×3.1×10^-3×60×0.3³/(1.8×10^-5×2.5×10^-5) = 1.1×10^8（层流） Nu = 0.59Ra^0.25 = 0.59×(1.1×10^8)^0.25 = 60.3 h = Nu·k/L = 60.3×0.028/0.3 = 5.6 W/(m²·K) Q = hA·ΔT = 5.6×(0.3×0.5)×60 = 50.4 W

挑战题

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

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

答案 水在40°C的性质：ν = 6.6×10^-7 m²/s，β = 3.85×10^-4 K^-1 Re = VD/ν = 0.1×0.05/(6.6×10^-7) = 7576 Gr = gβΔTD³/ν² = 9.8×3.85×10^-4×40×0.05³/(6.6×10^-7)² = 4.3×10^7 Gr/Re² = 4.3×10^7/7576² = 0.75 由于0.1 < Gr/Re² < 10，属于混合对流区域。水平管内，浮力产生二次流，会显著影响传热。顶部（冷却侧）和底部（加热侧）的传热系数会有很大差异。 建议采用混合对流关联式或分别计算强制和自然对流贡献。

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

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

答案 铝的导热系数k = 200 W/(m·K) 可用温差ΔT = 85 - 25 = 60°C 肋片参数m = √(hP/kA_c) = √(2h/kt)（矩形肋片） 假设肋厚t = 2mm：m = √(2×50/(200×0.002)) = 15.8 m^-1 肋效率η = tanh(mL)/(mL) 对于L = 30mm：mL = 0.47，η = 0.44/0.47 = 0.94 最优间距（基于边界层不重叠）： s_opt ≈ 2δ_t ≈ 2×3√(νx/V) ≈ 3-5 mm 选择：肋片高度30mm，间距4mm，厚度2mm 肋片数量：30/4 ≈ 7片 总散热面积：A_total = 基板面积 + η×肋片面积 ≈ 0.05 m² 验证：Q = hA_total·ΔT = 50×0.05×60 = 150W > 100W ✓

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

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

答案 选择水作为工作流体（60°C）： - 汽化潜热h_fg = 2.36×10^6 J/kg - 蒸汽密度ρ_v = 0.13 kg/m³ - 液体密度ρ_l = 983 kg/m³ - 表面张力σ = 0.066 N/m - 动力粘度μ_l = 4.7×10^-4 Pa·s 携带极限（entrainment limit）： Q_ent = A_v·h_fg·ρ_v·√(σρ_v/2) 毛细极限（capillary limit）： Q_cap = (2σ/r_eff)·(K·A_w/L_eff)·(ρ_l·h_fg/μ_l) 假设烧结粉末芯结构，r_eff ≈ 50μm，K ≈ 10^-10 m² 对于直径D = 10mm的热管： - 蒸汽通道面积A_v ≈ 0.5πD²/4 = 3.9×10^-5 m² - Q_ent ≈ 450W > 200W ✓ - Q_cap ≈ 280W > 200W ✓ 选择：直径10mm，长度200mm的水热管

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

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

答案 空气层平均温度T_m = 5°C，ΔT = 30°C 空气性质：ν = 1.4×10^-5 m²/s，α = 1.9×10^-5 m²/s，k = 0.024 W/(m·K) β = 1/278 = 3.6×10^-3 K^-1 Ra = gβΔTL³/(να)，其中L是间距 不同间距的Ra数： - L = 5mm：Ra = 660（导热主导，Nu ≈ 1） - L = 10mm：Ra = 5300（弱对流，Nu ≈ 1.2） - L = 15mm：Ra = 1.8×10^4（对流增强，Nu ≈ 2.0） - L = 20mm：Ra = 4.2×10^4（强对流，Nu ≈ 2.8） 传热量Q = Nu·k·A·ΔT/L 最小传热量出现在L ≈ 10mm处，此时对流刚开始但还不强烈。 更大间距虽然热阻增加，但对流增强抵消了这个优势。 最优设计：间距10-12mm，可减少70%的传热损失。

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

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

答案 基础参数： - 机柜尺寸：0.6m(宽)×2m(高)×1m(深) - 单机柜功率：10kW/m² × 1.2m² = 12kW - 所需风量（ΔT=20°C）：Q = ṁc_pΔT → ṁ = 0.6 kg/s 开放式布置问题： - 热回流：30-40%的热空气回流 - 实际温升：ΔT_real = ΔT_ideal/(1-γ) = 20/0.7 = 28.6°C - 入口温度需降至：70 - 28.6 = 41.4°C 冷通道封闭： - 消除热回流，γ ≈ 0 - 送风温度可提高至50°C - 节能效果：制冷COP提高20-30% 热通道封闭： - 防止热空气扩散 - 提高回风温度，增强economizer使用时间 - 可采用更高的送风温度（15-18°C） 推荐方案： 1. 热通道封闭 + 行间精密空调 2. 送风温度18°C，回风温度38°C 3. 变频风机根据负载调节 4. 预期PUE：1.3-1.4（vs 传统2.0）

常见陷阱与错误 (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

最佳实践检查清单

设计阶段

计算验证

优化考虑

实验验证