球的径向热传导模拟器 返回
热工程

球的径向热传导模拟器

对通过球壳(中空球壳)径向流动的定常热传导进行可视化的工具。改变内半径、外半径、热导率和内外表面温度时,热阻、热流量、热流密度和径向温度分布实时显示,可用于隔热罐和极低温杜瓦瓶的保温设计。

参数设置
内半径 r₁
m
球壳内侧半径(高温侧)
外半径 r₂
m
球壳外侧半径(低温侧)。必须大于 r₁
热导率 k
W/(m·K)
保温材 0.03~0.05、玻璃约 1、铜约 400
内表面温度 T₁
°C
外表面温度 T₂
°C
计算结果
热阻 R_th (K/W)
热流量 Q (W)
内表面热流密度 (W/m²)
外表面热流密度 (W/m²)
中间半径温度 (°C)
隔热性能判定
球壳断面 — 径向热流动画

球壳的横截面图。热从内侧(高温)向外侧(低温)呈放射状流动。同一热流量 Q 通过扩展的球面,外侧箭头密度更低。

温度分布 T(r)
热流量 Q vs 内半径 r₁
理论与主要公式

$$Q=\frac{T_1-T_2}{R_{th}},\qquad R_{th}=\frac{r_2-r_1}{4\pi k\,r_1\,r_2}$$

球壳定常热传导。Q:热流量 [W]、R_th:球壳热阻 [K/W]、T₁,T₂:内外表面温度、r₁,r₂:内外半径、k:热导率。热流量沿传导路径恒定,但热流密度按 q(r) = Q/(4πr²) 与 1/r² 成正比减小。

$$T(r)=T_1-(T_1-T_2)\,\frac{\dfrac{1}{r_1}-\dfrac{1}{r}}{\dfrac{1}{r_1}-\dfrac{1}{r_2}}$$

球壳内温度分布。T 是 r 而非 1/r 的一次函数,温度降幅集中在内侧。

球的径向热传导基础

🙋
"球的径向热传导"在什么情况下会出现?我在课上学过平板热传导,但球形的情况不太理解……
🎓
实际应用很常见啊,比如保温瓶内液、液化天然气巨大球形储罐、血液保存杜瓦瓶等。当需要把内部温度与外界隔离时,球壳壁就会有热通过径向传导,这就是模型。压力容器和航天燃料罐经常用球形,一是强度好,二是同体积外表面积最小,所以保温角度也占便宜。
🙋
外表面积小就是散热少,对吧。那球形无条件地更有利于隔热吗?
🎓
这是个有趣的地方,形状本身有利和同厚度壁面有利是两回事。同厚度比较的话,球外侧面积按 r² 展开,反而比平板和圆筒更容易透热。下面的感度图中,把内半径 r₁ 减小,热流量 Q 反而增大,你能看到这一点。理由是"小内表面承受整个温差",而外表面大但散热缓慢,内侧成了限速步骤。
🙋
哦,确实 Q 随着 r₁ 减小而增大呢。那加保温层时,外半径变大,反而散热更快?
🎓
正好说到"保温临界半径"这个经典问题,圆筒电线绝缘设计里经常出现。加保温层会"增加热阻",但"外表面积也增大",对流散热更容易。小球外加低导热率保温层时,可能出现加层反而热损失增加的反直觉现象。圆筒临界半径 r_c = k/h,球约为 r_c = 2k/h(h 是外表面对流系数)。本工具采纯传导边界条件,不含外侧对流效应,实际工程设计必须超过临界半径。
🙋
所以保温层不是越厚越好啊。还有,T(r) 曲线不是直线,内侧突然下降,这也是球的特点吗?
🎓
完全对。定常时,同一热流量 Q 要通过逐渐扩大的球面 4πr²,温度梯度 dT/dr 按 1/r² 迅速变小。积分后 T(r) 就成了 1/r 的一次函数,内侧陡降、外侧缓降。本工具的图就显示这个"内快外缓"。与平板(直线)和圆筒(对数)比,球最集中在内侧。比如 20 cm 直径热球用 1 cm 薄壳包裹,温度降幅有 3~4 mm 集中在壳的最内层,意思是保温层内侧品质和贴合度决定整体性能。
🙋
很有意思。在有限元模拟球体传热时,用这个解析解能快速验证,不用特别细致的网格。
🎓
完全对,这就是"理智检验"的核心。FEM 模拟球壳出的热流量跟这个式子偏几个数量级,说明边界条件或单位肯定哪里错了。如果吻合,再加对流和辐射看偏差有多大,就有了量化的出发点。极低温、隔热、压力容器领域,这个公式是最先写上去的基本式,设计现场人人都用。

常见问题

定常情况下,内表面 r₁、外表面 r₂、热导率 k 的球壳径向热传导热阻为 R_th = (r₂ − r₁) / (4π k r₁ r₂)。这与平板 R = L/(kA) 和圆筒 R = ln(r₂/r₁)/(2π k L) 一样,遵循热版本的欧姆定律 ΔT = Q·R_th。公式中的 (r₁·r₂) 对应内外半径的几何平均面积,是对径向展开的传热面积的平均化表示。
同样厚度对比,球外侧的面积按 r² 展开,因此热流量比平板和圆筒更大,单纯意义上不是"有利于隔热"。但球形状内表面积与外表面积的比值最小,所以单位容积的表面积(即热损失的入口)最少。极低温球形杜瓦瓶和液化天然气球形罐选择球形正是这个原因,只有配合保温层厚度才能实现低热损失。
不一定。这是著名的"保温临界半径"效应。给小球外部加低导热率保温层时,外表面积增大,对流散热更容易,可能出现加保温层反而热损失增加的反向现象。临界半径圆筒为 r_c = k/h,球为 r_c = 2k/h(h 为外表面热传递系数)。本工具采用内外表面温度边界条件的纯传导模型,不直接包含外侧对流效应。实际工程必须选择超过临界半径的厚度。
定常径向热流的球壳中,相同热流量 Q 依次通过更大的球面 4πr²,导致温度梯度 dT/dr 按 1/r² 快速减小。积分后温度是 1/r 的一次函数。即 T(r) 在内侧快速下降,外侧缓和下降。本工具的 T(r) 图表显示"内侧快、外侧缓"的曲线。与平板(直线)和圆筒(对数)比较,球的温度降幅最集中在内侧。

实际应用案例

极低温储罐、液化天然气球形罐:液氢、液氮储存杜瓦瓶、-162 ℃ 液化天然气 Moss 型球形罐都是这个球壳热传导问题的直接应用。罐内壁低温,外壁外气温,中间有多层真空隔热(MLI)或珍珠岩填充层。设计现场的关键指标是"单位时间蒸发百分率(煤气化率)",通过 Q = ΔT/R_th 反推保温材的热导率和厚度。本工具的内半径感度图能大致显示增加容量(增大 r₁)如何影响单位体积热损失。

压力容器、核反应堆安全壳:球形同体积外表面积最小,内压应力分布最均匀,高压气体罐和某些核反应堆安全壳采用球形。事故工况传热分析时,内侧高温高压流体、外侧空气或冷却水,边界条件用同样的热阻网络(各层热阻串联求和)。多层结构(钢、隔热材、混凝土)时,分别计算 R_th 然后加总,这是基本做法。

球形保温材、催化剂粒子、电子器件发热分析:催化剂球形颗粒、蓄热球、半导体封装内焊珠、甚至冰淇淋勺子这样的"球形物体"温度响应,首先用这个定常解做初步评估。从半径和热导率的组合推断热时常 τ ≈ ρcr²/k 的量级,判断"秒级升温还是时间级",再进瞬态分析,这是标准流程。

地球和行星科学:地球内部热输运、行星冷却模型、月球和冰质天体地下海温度评估都用球壳传导公式作出发点。实际上涉及对流、辐射、潜热、内部发热,但首先得估"假如纯传导,行星表面要放散多少瓦特",然后判断对流是否需要。研究人员的电子表格第一行往往就是 R_th = (r₂−r₁)/(4πkr₁r₂) 这个式子。

常见误区与注意

最大陷阱是"球的热阻公式与平板、圆筒混淆"。R = (r₂−r₁)/(4πkr₁r₂) 只对球形,平板 R = L/(kA) 和圆筒 R = ln(r₂/r₁)/(2πkL) 形式完全不同。同厚度下,球因为 r₁r₂ 在分母,大半径时阻值更小。多层伝热网络设计时,若各层都用平板"面积恒定"会双向出错:薄层用平板会低估、厚层会高估。必须按形状选对应的阻值公式。

其次,"内外表面温度边界条件不含实际对流和辐射"。本工具是纯传导模型,用户直接指定内外壁温。现实中知道"内液温"和"外气温",但壁面温度由对流系数 h_in, h_out 决定,是未知量。实际热流要把内对流阻 1/(h_in·A_in) 和外对流阻 1/(h_out·A_out) 与本式的 R_th 串联。忽视这点,本工具值会偏大(只算壁内传导)。用作"理智检验"时要心里有数。

最后,"热导率 k 误认为常数"。实际保温材温度依赖性很强,聚氨酯泡沫和真空多层在低温时性能快速变化,高温时辐射成分按温度三次方增长,使得"表观 k"剧烈变化。极低温设备设计时用"温度范围积分有效 k",不用中点温度单值。本工具固定 k 时,应取内外温度中点附近的 k 值作第一近似,精细设计务必查材料温度曲线。

使用指南

  1. 输入内球半径(mm)和外球半径(mm)以设置球壳尺寸。例如:内球半径50mm、外球半径100mm的聚氨酯泡沫隔热罐
  2. 指定材料热导率(W/m·K)。聚氨酯0.025、玻璃棉0.04、岩棉0.05是典型值
  3. 输入内外表面温度(°C)以确定边界条件,热阻、热流量和各层温度分布实时计算
  4. 查看中间半径的温度分布曲线,观察从内侧到外侧的温度变化趋势

计算示例

液氮储存杜瓦瓶球壳设计:内球半径 r₁=40mm、外球半径 r₂=80mm、聚氨酯填充(λ=0.026 W/m·K)、内表面温度 Tin=77K(-196°C)、外表面温度 Tout=293K(20°C) 时,热阻 Rth=(r₂-r₁)/(4πλr₁r₂)=(0.04)/(4π×0.026×0.04×0.08)≈3.04 K/W,热流量 Q=(293-77)/3.04≈71W。内表面热流密度约 710 W/m²,外表面约 355 W/m²。

工程实务注意

  1. 空气层热导率约 0.026 W/m·K,真空抽至 1×10⁻⁴Pa 以下可降至 0.001 以下。二层壁间压力越低,辐射屏蔽越重要
  2. 聚氨酯泡沫施工不均会局部密度下降,热导率可升至 0.04。务必进行样品测试(JIS A1412 标准)
  3. 大温差时材料热导率有温度依赖,应用平均温度修正值。极低温环境数值会显著改变,推荐用中温值重算
  4. 外表面结露环境下考虑加隔热层或电加热带,计算中需包含