3. 传热相似准则数
在化工放大中,传热过程的相似放大依赖于一组无量纲准则数。它们不是凭空定义的符号,而是从对流传热微分方程组出发,通过量纲分析和方程无量纲化得到的“物理意义标签”。理解每个准则数的来源与含义,才能在放大时判断哪些相似条件必须保持、哪些可以放松。
一、从对流传热微分方程组出发
以不可压缩牛顿流体、常物性、无内热源的二维稳态对流传热为例,控制方程包括连续性方程、动量方程(Navier–Stokes 方程)和能量方程。取特征长度 L、特征速度 u∞、特征温度差 ΔT,定义无量纲量:
X = x/L, Y = y/L, U = u/u∞, V = v/u∞
P = p/(ρu∞²), θ = (T - T∞)/ΔT
代入动量方程后,惯性项系数为 1,黏性项前出现组合 ρu∞L/μ,即雷诺数 Re。代入能量方程后,对流项系数为 1,导热项前出现组合 μcp/λ(即 Pr)与 ρu∞L/μ(即 Re)的乘积。这说明:动量传递与热量传递的无量纲方程由 Re 和 Pr 共同控制。
二、五个准则数的物理意义
雷诺数 Re = ρuL/μ:惯性力与黏性力之比。Re 小表示黏性主导、流动层流;Re 大表示惯性主导、易失稳转捩为湍流。它决定流场结构,从而间接决定传热强弱。
普朗特数 Pr = μcp/λ = ν/a:动量扩散率(运动黏度 ν)与热扩散率(导温系数 a)之比。Pr 反映速度边界层与温度边界层的相对厚度:Pr≈1 时两层近似重合(如气体);Pr≫1 时温度边界层更薄(如油类、高分子熔体);Pr≪1 时温度边界层更厚(如液态金属)。
努塞尔数 Nu = hL/λ:壁面对流换热量与同温差下纯导热量的比值,是无量纲化的对流传热系数。Nu 不是物性数,而是待求的结果量,通常表示为 Nu = f(Re, Pr) 或 Nu = f(Gr, Pr)。
格拉晓夫数 Gr = gβΔT L³/ν²:浮升力与黏性力之比,衡量自然对流的强弱。Gr 越大,浮升驱动越强,自然对流越剧烈。它与 Re 的平方在强制对流中地位类似,因此自然对流常用 Gr 与 Pr 的组合(瑞利数 Ra = Gr·Pr)关联。
斯坦顿数 St = Nu/(Re·Pr) = h/(ρcpu):实际对流传热量与流体热容量输运量之比,直接反映传热效率。在传热与摩擦类比(如 Colburn 类比 jH = St·Pr2/3 = f/2)中,St 是核心量,便于把传热与压降关联起来。
三、推导思路:以 Nu 与 St 的关系为例
由定义 Nu = hL/λ,两边同除以 Re·Pr:
Nu/(Re·Pr) = (hL/λ) / [(ρuL/μ)(μcₚ/λ)]
= hL/λ · μλ/(ρuL·μcₚ)
= h/(ρcₚu) = St
这一步说明 St 并非独立新量,而是 Nu、Re、Pr 的组合。它的价值在于把“传热系数”换算成“无量纲热流”,从而与摩擦因子 f 建立类比关系。工程上常用 Colburn 类比 jH = St·Pr2/3 ≈ f/2 来由压降估算传热,或反之由传热反推阻力,这在换热器放大与板式/翅片结构选型中非常实用。
四、放大中的实际判断
假设某夹套搅拌釜小试数据:桨径 D = 0.1 m,转速 N = 300 r/min,物料密度 ρ = 1000 kg/m³、黏度 μ = 0.01 Pa·s、比热 cp = 4000 J/(kg·K)、导热系数 λ = 0.6 W/(m·K)。可先算 Re = ρND²/μ ≈ 1000×5×0.01/0.01 = 5000,属过渡区;Pr = μcp/λ ≈ 66.7。若放大到 D = 0.5 m,要保持 Re 相同,则 N 须降为原来的 1/25,即 12 r/min。此时单位体积传热面积下降,往往需要改用更大面积的换热结构或提高温差来补偿。这类“先算准则数、再定操作参数、最后校核 h 与压降”的路径,正是相似准则在放大中的落地方式。