Dittus-Boelter关联式 — 管内湍流传热基本公式
理论与物理
概述 — 什么是 Dittus-Boelter 公式
老师，Dittus-Boelter 公式在哪里使用呢？虽然每章教科书都有出现，但说实话我还没太理解…
这是估算管内湍流换热系数 $h$ 最简便的公式。公式形式如下：
流体被加热时使用 $n = 0.4$，被冷却时使用 $n = 0.3$。冷却水管道、散热片间流道、化工厂的热交换器——在这些“流体在管内湍流流动”的场景的初步设计中，这是最先使用的公式。
诶，仅凭这么简单的公式就能知道换热系数吗？只需要输入 Re 和 Pr？
是的，所以它才是“初步估算之王”。例如设计汽车散热器管道时，只要确定了冷却水流速和管道内径，计算出 Re，再代入水的 Pr（常温下约 7），就能立刻得到 $h$ 的估算值。这样就能对“这个管径是否足够冷却？”做出最初的判断。
但是您说“估算”，意思是也有靠不住的情况对吧？
问得好。这个公式有前提条件。仅适用于 Re > 10,000 的充分发展湍流，且管道足够长（$L/D > 60$）的情况。对于入口效应强的短管，或 Re = 2,300〜10,000 的过渡流区域，Gnielinski 公式要准确得多。如果不了解这个“适用范围”就盲目使用，可能会轻易引入 2〜3 成的误差。
核心公式
整理构成 Dittus-Boelter 公式的两个公式。
努塞尔数关联式：
此处各无量纲数的定义如下：
向换热系数的转换：
| 符号 | 物理量 | SI单位 |
|---|---|---|
| $D$ | 管道内径 | m |
| $u$ | 管内平均流速 | m/s |
| $\rho$ | 流体密度 | kg/m³ |
| $\mu$ | 动力粘度（粘性系数） | Pa·s |
| $c_p$ | 定压比热 | J/(kg·K) |
| $k$ | 流体的热导率 | W/(m·K) |
| $h$ | 对流换热系数 | W/(m²·K) |
| $n$ | Pr 指数（加热 0.4／冷却 0.3） | — |
适用条件与局限
适用条件，能再具体点告诉我吗？除了“Re > 10,000”还有其他限制对吧？
嗯，总结起来是这样的：
| 条件 | 范围 | 违反时的后果 |
|---|---|---|
| 雷诺数 | $\mathrm{Re}_D > 10{,}000$ | 过渡流/层流时会被高估。应使用 Gnielinski 公式 |
| 普朗特数 | $0.7 < \mathrm{Pr} < 160$ | 液态金属（Pr≪1）或高粘度油（Pr≫160）时误差很大 |
| 管长/直径比 | $L/D > 60$ | 入口效应会使实际 h 偏高。需要入口修正系数 |
| 壁面温度与流体温度差 | 中等程度（$\Delta T$ 较小） | 大温差时粘度变化大→应使用 Sieder-Tate 公式 |
$L/D > 60$ 就是说，内径 20mm 的管道需要 1.2m 以上才行吗？
就是这个意思。电子设备的冷板或小型散热片的翅片间流道，很多时候 $L/D = 10〜30$ 左右。对于这种短管，入口附近局部换热会非常高，如果假设 Dittus-Boelter 公式给出均匀的 $h$，就会导致低估。实际工作中，乘以入口修正系数 $(1 + (D/L)^{0.7})$ 或切换到 Gnielinski 公式会更安全。
各项的物理意义
Re0.8 — 为什么是 0.8 次方
雷诺数是“惯性力/粘性力”之比，是湍流强度的指标。指数不是 1 而是 0.8，是因为湍流区域的 Nu 数随 Re 的增加并非线性，而是稍显平缓。这源于湍流边界层的发展模式——壁面附近的粘性底层随流速增加而变薄但并未完全消失。许多实验已证实，与实验数据的最佳吻合点在 0.8 附近。
Prn — 加热与冷却为什么指数不同
普朗特数是“动量扩散/热扩散”之比。流体被加热时，壁面附近的粘度降低，速度边界层变薄→换热增强→采用较大的指数（0.4）。被冷却时，壁面附近的粘度增加，速度边界层变厚→换热抑制→采用较小的指数（0.3）。这种不对称性，是简易地纳入了温度依赖的粘度变化影响。
系数 0.023 的意义
这个系数并非理论推导得出，而是通过大量圆管流动实验数据的回归分析确定的经验常数。根据 Winterton (1998) 的文献，原论文的实验数据约有 ±25% 的离散，0.023 是其中间值。
历史背景
公式的真正提出者问题
被称为“Dittus-Boelter 公式”的 $\mathrm{Nu} = 0.023\,\mathrm{Re}^{0.8}\,\mathrm{Pr}^n$ 源自 F.W. Dittus 和 L.M.K. Boelter 于 1930 年发表的论文。但近年来的文献调查（Winterton, 1998）发现，指数 $n$ 的加热/冷却区别使用（0.4/0.3）是 Boelter 在 1942 年单独修正的，在 Dittus 的原论文中并不存在。此外，原论文并非正式的同行评审期刊，而是加州大学的内部报告，原本的获取本身就困难。尽管如此，由于被引用过于广泛，“Dittus-Boelter 公式”的名称已经固定下来——这是工程史上一个有趣的轶事。
数值解法与实现
换热系数的计算步骤
请一步步教我实际求 $h$ 的步骤。
步骤如下：
- 确定膜温度：$T_f = (T_w + T_b)/2$（壁面温度 $T_w$ 与整体温度 $T_b$ 的平均值）
- 获取物性值：从物性表中读取 $T_f$ 下的 $\rho$, $\mu$, $c_p$, $k$
- 计算 Re：$\mathrm{Re}_D = \rho u D / \mu$
- 计算 Pr：$\mathrm{Pr} = \mu c_p / k$
- 确认适用条件：Re > 10,000 且 0.7 < Pr < 160 且 L/D > 60
- 计算 Nu：$\mathrm{Nu}_D = 0.023\,\mathrm{Re}_D^{0.8}\,\mathrm{Pr}^n$
- 计算 $h$：$h = \mathrm{Nu}_D \cdot k / D$
如果不知道膜温度——比如设计初期壁面温度还没确定时怎么办？
问得好。可以先假设壁面温度进行计算，然后用得到的 $h$ 重新计算壁面温度→再次更新物性值→重新计算 $h$… 进行迭代。实际工作中，2〜3 次迭代就足够收敛了。对于像水这样物性随温度变化不大的流体，即使使用整体温度来取物性值，误差也在百分之几以内。
计算示例 — 冷却水管道
我想看具体的数字。比如发动机的冷却水系统之类的。
那我们用典型的发动机冷却水条件来试试看。
条件：内径 $D = 0.02$ m，流速 $u = 1.5$ m/s，水温 80°C（加热状态）
80°C 水的物性值：
- $\rho = 972$ kg/m³
- $\mu = 3.55 \times 10^{-4}$ Pa·s
- $c_p = 4,197$ J/(kg·K)
- $k = 0.670$ W/(m·K)
步骤 1 — Re：
$$ \mathrm{Re} = \frac{972 \times 1.5 \times 0.02}{3.55 \times 10^{-4}} = 82{,}028 $$→ 完全湍流。OK。
步骤 2 — Pr：
$$ \mathrm{Pr} = \frac{3.55 \times 10^{-4} \times 4{,}197}{0.670} = 2.22 $$→ 在范围内。OK。
步骤 3 — Nu（因为是加热，所以 $n = 0.4$）：
$$ \mathrm{Nu} = 0.023 \times 82{,}028^{0.8} \times 2.22^{0.4} = 0.023 \times 10{,}173 \times 1.37 = 320.6 $$步骤 4 — $h$：
$$ h = \frac{320.6 \times 0.670}{0.02} = 10{,}740 \;\text{W/(m²·K)} $$也就是说约 $h \approx 10{,}700$ W/(m²·K)。这是汽车冷却水系统中典型的值。
哦，相当大的值呢。空气的强制对流大概是多少？
空气的话 Pr ≈ 0.71，$k$ ≈ 0.03 W/(m·K)，比水小一个数量级，所以即使相同的流速和管径，$h$ 也只有 50〜200 W/(m²·K) 左右。水的 $h$ 是空气的 50〜100 倍，是因为水的热导率和热容量要大得多。所以发动机冷却才用水。
物性值的评估温度
在哪个温度下取物性值——实际工作中最容易被忽视的点
Dittus-Boelter 公式原本是从壁面温度与整体温度差较小的实验数据推导出来的。因此，物性值应在膜温度 $T_f = (T_w + T_b)/2$ 下评估是标准做法。不过，在壁温未知的阶段，通常先用整体温度代替，求得 $h$ 后再重新计算壁温并更新物性值，采用这种迭代法。
当壁面与流体温差较大时（例如：高温油在管道壁面被急冷的情况），壁面附近的粘度会发生很大变化。此时 Dittus-Boelter 公式无法充分应对，应切换到包含粘度修正项 $(\mu_b/\mu_w)^{0.14}$ 的 Sieder-Tate 公式。
实践指南
CAE 中的边界条件设置
用 CAE 软件进行热分析时，Dittus-Boelter 公式怎么用呢？是把公式输入软件吗？
主要有两种用法：
- 作为对流边界条件手动设置 $h$：不模拟流体侧，只解析固体部分的情况。将 Dittus-Boelter 公式计算出的 $h$ 和流体整体温度 $T_\infty$ 指定给壁面。这是 Ansys Mechanical 或 Abaqus 热分析中最常见的方法。
- CFD 求解器内部自动计算：在 Fluent 或 STAR-CCM+ 中直接求解管内流动时，求解器会通过壁面函数自动计算 $h$。这种情况下不需要手动使用 Dittus-Boelter 公式。相反，其用法变成了作为 CFD 结果合理性验证，与 Dittus-Boelter 公式的估算值进行比较。
也就是说，“不解流体的热分析”时 Dittus-Boelter 公式才有用武之地。反过来，如果用 CFD 直接求解，就作为验证用？
正是如此。实际工作中常见的是“设计初期用 Dittus-Boelter 公式进行概算→详细设计阶段用 CFD”这样的流程。如果初期概算的 $h$ 与 CFD 结果差异很大，很可能是因为超出了适用条件。
与 Gnielinski 公式的选用区分
请再详细讲讲与 Gnielinski 公式的区别。我经常不知道该用哪个…
我们来看看对比表：
| 项目 | Dittus-Boelter 公式 | Gnielinski 公式 |
|---|---|---|
| 公式形式 | $0.023\,\mathrm{Re}^{0.8}\,\mathrm{Pr}^n$ | $\frac{(f/8)(\mathrm{Re}-1000)\mathrm{Pr}}{1+12.7\sqrt{f/8}(\mathrm{Pr}^{2/3}-1)}$ |
| Re 适用范围 | > 10,000 | 2,300 〜 5×10⁶ |
| Pr 适用范围 | 0.7 〜 160 | 0.5 〜 2,000 |
| 入口效应 | 未考虑（以 L/D > 60 为前提） | 可通过修正项应对 |
| 精度（与实验吻合度） | ±25% | ±10% |
| 计算工作量 | 计算器即可立即计算 | 需要额外计算摩擦系数 $f$ |
| 最佳用途 | 初期概算、灵敏度分析 | 详细设计、报告依据 |
那么需要精度时就选 Gnielinski，Dittus-Boelter 是“想粗略看看的时候”用的感觉吗？
基本是这样。不过，在 Re > 10,000 的完全湍流区域，两者的差异大多在百分之几以内。也就是说，对于充分发展湍流的长管，Dittus-Boelter 就足够了。相反，Gnielinski 公式真正发挥价值是在过渡流区域（Re = 2,300〜10,000）或入口效应大的短管。
与 Sieder-Tate 公式的选用区分
我也想了解与 Sieder-Tate 公式的区别。
Sieder-Tate 公式的形式是 $\mathrm{Nu} = 0.027\,\mathrm{Re}^{0.8}\,\mathrm{Pr}^{1/3}\,(\mu_b/\mu_w)^{0.14}$，包含了壁面与流体整体的粘度比 $(\mu_b/\mu_w)^{0.14}$。
也就是说，这是用于壁面温度与流体温度差很大，粘度变化显著的情况的公式。例如原油管道壁面被加热器加热的情况，壁面附近的粘度会急剧下降。这种情况下 Dittus-Boelter 公式就不够用了，Sieder-Tate 公式的粘度修正项就会起作用。
反过来，对于水或空气这样粘度随温度变化很小的流体，Dittus-Boelter 公式和 Sieder-Tate 公式的结果几乎一致。
工业应用实例
| 领域 | 典型条件 | 推算的 $h$ [W/(m²·K)] | 备注 |
|---|---|---|---|
| 汽车散热器管道 | 水，$u$ = 1〜2 m/s，$D$ = 15〜25 mm | 5,000〜15,000 | 在 80°C 附近使用 |
| 化工厂热交换器 | 有机溶剂，$u$ = 0.5〜3 m/s | 500〜5,000 | Pr 较大，因此 h 相对较高 |
| 空调风管 | 空气，$u$ = 5〜15 m/s，$D$ = 100〜300 mm | 30〜100 | 空气的 k 小，因此 h 也小 |
| 数据中心液冷 | 惰性液体，$u$ = 0.5〜2 m/s | 2,000〜8,000 | 微通道需注意 L/D |
| 燃气轮机冷却通道 | 空气，$u$ = 30〜100 m/s | 200〜1,000 | 高 Re、旋转效应需要修正 |
在热水管道设计中的常用
在日本的热水系统设计（依据 JIS B 8417）中，广泛使用 Dittus-Boelter 公式进行管内强制对流的概算，流速 0.5〜2 m/s 的热水管道可得到 $h$ = 3,000〜8,000 W/(m²·K) 左右。在计算管道保温设计“通过包裹保温材料能减少多少散热”时，首先需要管内侧的 $h$。
相关主题
この記事の評価ご回答ありがとうございます！参考に
なったもっと
詳しく誤りを
報告
なった
詳しく
報告