3. 集总参数法及其适用性
当物体内部导热热阻远小于表面对流换热热阻时,物体内部温度趋于均匀,此时可以忽略空间温度梯度,把整个物体视为一个温度均匀的“集总热容”。这就是集总参数法的物理前提。它把原本依赖空间坐标的偏微分导热方程,简化为只随时间变化的一阶常微分方程,使瞬态温度响应可以用解析式直接求解。
一、适用性判据:Biot 数
判断能否采用集总参数法,核心是无量纲数 Biot 数。其定义为物体内部导热热阻与表面对流换热热阻之比:
Bi = h·Lc / k
其中 h 为表面换热系数(W/(m²·K)),k 为材料导热系数(W/(m·K)),Lc 为特征长度(m)。对于厚度为 2δ 的平板、半径为 R 的圆柱或球体,若两侧或整个表面参与换热,特征长度取 Lc = V/A(体积与换热表面积之比):平板 Lc = δ,长圆柱 Lc = R/2,球体 Lc = R/3。注意不要把几何半径直接当作特征长度,这是工程计算中最常见的错误之一。
工程上普遍采用的判据是 Bi ≤ 0.1。此时物体内部最大温差相对于表面与流体之间的总温差不超过约 5%,把内部温度视为均匀带来的误差在多数工程问题中可以接受。Bi 越小,集总假设越精确;当 Bi 接近或超过 0.1 时,内部温度梯度不可忽略,应改用一维非稳态导热的近似解或数值方法。
二、集总热容法的温度响应
对控制体做能量平衡:表面换热量等于内能变化率,即
h·A·(T∞ − T) = ρ·V·c·(dT/dt)
引入过余温度 θ = T − T∞,初始条件 θ(0) = θ₀ = T₀ − T∞,积分得:
θ/θ₀ = exp(−t/τ), τ = ρ·V·c / (h·A) = ρ·c·Lc / h
式中 ρ 为密度,c 为比热容,τ 为时间常数(s)。温度按指数规律单调趋近流体温度 T∞,不会越过平衡值。若已知某一时刻温度,也可反解时间:t = τ·ln(θ₀/θ)。
下面用一张图概括从判据到响应的完整逻辑,以及时间常数对曲线形态的影响。
三、时间常数的物理含义与工程用法
时间常数 τ = ρ·c·Lc/h 集中体现了三个因素:材料的热惯性 ρc 越大,温度变化越慢;特征长度 Lc 越大,单位表面积对应的热容越大,响应越慢;表面换热系数 h 越大,热量进出越容易,响应越快。因此强化冷却或加热的途径,要么提高 h(如强制风冷、喷淋、搅拌),要么减小 Lc(如把大块物料切成薄片、采用小直径颗粒或细丝)。
τ 的工程意义在于给出“进度刻度”:经过 1τ,过余温度降为初始值的 36.8%;2τ 约 13.5%;3τ 约 5.0%;4τ 约 1.8%;5τ 约 0.7%。工程上常把 3τ~5τ 视为基本达到热平衡的时间。反过来,若要求某时刻温度不超过给定值,可由 t = τ·ln(θ₀/θ) 估算所需时间,或由 τ = t/ln(θ₀/θ) 反推允许的最大时间常数,再据此选择 h 或尺寸。
四、示例计算与结果判读
设一直径 10 mm 的钢球(k = 45 W/(m·K),ρ = 7800 kg/m³,c = 460 J/(kg·K)),初始温度 500 ℃,突然置于 25 ℃ 的空气流中,h = 80 W/(m²·K)。球体 Lc = R/3 = 0.005/3 ≈ 1.667×10⁻³ m。
Bi = 80 × 1.667e-3 / 45 ≈ 0.00296 ≪ 0.1 → 集总假设成立
τ = 7800 × 460 × 1.667e-3 / 80 ≈ 74.7 s
t(θ/θ₀=0.05) = τ·ln(20) ≈ 74.7 × 3.00 ≈ 224 s
即约 3.7 min 后球温降至约 25 + 0.05×475 ≈ 49 ℃。若把 h 提高到 400 W/(m²·K)(强制风冷),τ 降为约 14.9 s,同样冷却到 5% 过余温度只需约 45 s,说明 h 与冷却时间近似成反比。若把钢球换成直径 50 mm 的球,Lc 增大 5 倍,Bi 仍远小于 0.1,但 τ 也增大 5 倍,冷却明显变慢——这正是“减小尺寸”在热处理和颗粒冷却中常被优先采用的原因。
综上,集总参数法以 Bi ≤ 0.1 为适用边界,用一阶指数响应描述温度随时间的变化,时间常数 τ 则把材料热惯性、几何尺度与表面换热强度统一到一个可比较的尺度上,是瞬态换热估算与方案筛选中最简洁有效的工具之一。