MATERIALS NOTES阅读地图
8个主要章节
4张图片与图示
1组对照与手算表

把概念、模型计算与实测证据分开,阅读时留意条件、单位和比较口径。

章节速览8 节 · 展开选择

稳态导热里,流入与流出相等,材料内部不再积蓄能量。如果把温热的小物体放进较冷的流体,它的温度仍在变化,此时流出的热量来自物体储存的内能,不能继续套用“各处热流率都一样”的稳态结论。

本文用一个随时间变化的温度T(t)近似整个物体,推导集总热容模型。所有数字都是教学合成参数,不对应真实材料牌号、人体降温、食品安全或设备认证。这个模型最重要的工作,不是按计算器上的指数键,而是判断“一个温度代表整体”是否有依据。

先检查近似内部扩散快,表面对流慢

用Biot数比较内部导热与边界换热,而不是仅凭“物体很小”决定能否使用集总模型。

再计算过程热阻 × 热容 = 时间常数

从能量守恒得到温度轨迹,核对初始值、长期极限、单位和累计散热。

一套假设对应一个温度变量

设物体起始温度均匀为T₀,体积V、密度ρ、比热c、导热系数k均为正常数。整个暴露面积Aₛ以相同且恒定的换热系数h与温度恒定为T∞的充分大流体环境换热;忽略辐射、内部热源、相变和物性随温度的变化,几何尺寸固定。

这些条件分工不同。大环境使T∞近似不随吸热改变;常h让边界传热规律保持线性;内部近似等温才允许用单个T(t)描述物体。若只有一部分表面暴露,Aₛ只计算该部分的有效换热面积,不能把绝热底面也一并加上。若不同表面的h或环境温度不同,需要分区求和,而不是悄悄复用同一个hAₛ。

C = ρVc|热容
单位J/K;温度降低1K对应内能减少C焦耳。这里是总热容,不是比热c。
Rₕ = 1/(hAₛ)|对流热阻
单位K/W;只描述物体与流体边界的换热,不代表物体内部没有任何温差。
Lc = V/Aₛ|特征长度
单位m;本文Biot数与Fourier数统一采用这一定义,不能随意换成直径而不改数值。

从能量守恒推导指数变化

取离开物体的热流率为正,Q̇出=hAₛ(T−T∞)。物体内能变化率为C·dT/dt,守恒要求C·dT/dt=−Q̇出。冷却时温差为正,导数为负;若环境反而更热,两个符号同时反转,模型也能描述加热。

公式C1|集总能量方程与解
C dT/dt = −hAₛ(T − T∞)
τ = C/(hAₛ) = RₕC
T(t) = T∞ + (T₀ − T∞) exp(−t/τ)

恒定物性、h及T∞,内部近似等温,无内部热源、辐射与相变;t≥0。τ单位s,指数无量纲,温差用K或°C时数值相同。来源与推导对照文献1。

令θ=T−T∞,就有dθ/dt=−θ/τ。分离变量积分得到ln|θ(t)|−ln|θ(0)|=−t/τ,结合初始条件得到公式C1。T₀=T∞时直接得到常温解,不需要对零取对数。τ增大,温差衰减更慢;“密度更大散热更慢”只是在其余量固定时由热容增加得到的结论,不能跨不同几何和表面条件直接比较。

当t=τ,剩余温差是初始温差的1/e,约36.8%,并不是“降到初温的36.8%”。例如环境20°C、初温80°C,一个时间常数后的温度是20+60/e≈42.073°C。直接算80/e≈29.43°C,错在把摄氏温标原点当作物理零温差。

先算Biot数,再信任这一条曲线

内部导热并不是瞬间发生的。用k衡量内部传热能力,h衡量边界换热能力,构造Bi=hLc/k。把同一量级的热流估算式hAₛΔT外与(kAₛ/Lc)ΔT内比较,可以理解为ΔT内/ΔT外约随Bi变化。Bi很小,内部温差相对于物体与环境的温差较小,才有理由忽略内部不均匀。

这是量级判断,不是对任何形状和任何时间都严格成立的误差上界。本文不把一个通用阈值宣称为“所有温度误差一定小于某百分比”的认证。若Bi不小、初始温度不均或局部热源显著,应转向空间分布的传热模型;表面测温也不必然等于内部平均温度。

公式C2|两个无量纲数承担不同工作
Bi = hLc/k;α = k/(ρc);Fo = αt/Lc²
采用同一个Lc = V/Aₛ时:t/τ = Bi · Fo

Bi比较内外传热阻力的量级;Fo衡量扩散经历的无量纲时间。α单位m²/s。恒定正参数与相同特征长度是乘积关系的前提;Fo大本身不证明整个过程可用集总近似。参考文献1。

一组完整手算:多久到30°C

设V=10⁻⁶m³,Aₛ=6×10⁻⁴m²,ρ=2500kg/m³,c=800J/(kg·K),k=100W/(m·K),h=120W/(m²·K),T₀=80°C,T∞=20°C。可把几何想象为边长1cm、六面均暴露的立方体,但参数仍是合成值。

由此m=ρV=0.0025kg,C=2J/K,hAₛ=0.072W/K,Rₕ=125/9K/W,τ=250/9s≈27.778s。Lc=1/600m,Bi=0.002,很小;α=5×10⁻⁵m²/s。这里几何与面积已明确,不能把1cm边长直接替代Lc,得到的Bi会相差六倍。

时刻 剩余温差比例 预测温度°C 累计散热J
0 1 80 0
τ≈27.778s e⁻¹≈0.367879 42.073 75.854
2τ≈55.556s e⁻²≈0.135335 28.120 103.760
3τ≈83.333s e⁻³≈0.049787 22.987 114.026

目标30°C位于初温与环境温度之间,温差比例是(30−20)/(80−20)=1/6,因此t*=τln6≈49.771s。比2τ稍短,与上表2τ时已经降到28.120°C一致。到达环境20°C则只在t趋于无穷时发生;若问“20°C附近”,必须给容许温差,不能把“附近”当成精确等号。

归一化温差在零、一、二、三个时间常数处为1、0.367879、0.135335、0.049787,目标30摄氏度位于约1.792个时间常数。
原创合成模型曲线;横轴t/τ、纵轴温差比,标注点由公式C1计算,不是材料实测。查看原图
公式C3|目标时刻与累计散热
t* = τ ln[(T₀ − T∞)/(T* − T∞)]
Q出(0→t) = C[T₀ − T(t)]

第一式要求目标严格位于T₀与T∞之间,且T₀≠T∞;目标等于初温时t*=0。第二式取散热为正,加热时结果为负。t单位s,Q单位J;公式来自C1积分与能量守恒。

到30°C时,累计散热为2×(80−30)=100J;初始散热率为0.072×60=4.32W,目标时只有0.72W。用初始4.32W一直乘时间会高估散热,因为温差不断减小。全过程散热上限是C(T₀−T∞)=120J,也能用于发现单位或符号错误。

环境分两段时,温度不会跳变

前20s维持20°C环境,此后流体突然改为50°C,同时h仍为120。第一段末温T₁=20+60exp(−20/τ)≈49.205°C。环境发生阶跃,物体热容有限,物体温度仍连续;第二段应以49.205°C作为初值、以50°C作为新环境温度重新套用C1。

再经历20s,温度约为49.613°C。物体在第一段冷却,第二段反而缓慢升温,因为49.205低于新环境50。若第二段继续用80°C作初值,就等于凭空恢复了一份内能;若把T直接设成50°C,又把有限热容变成了零。

改参数时,也要重新检查模型

把h提高100倍,其他参数不变,公式给出的τ缩短到约0.278s,但Bi同时变成0.2。不能只取“100倍更快”的好处,继续忽略内部温差。公式作为一个微分方程仍能计算,并不代表它仍能描述原物体。

若几何相似地将所有长度扩大两倍,体积增八倍、面积增四倍,恒定h及物性下τ增两倍,Bi也增两倍。这个结论是在“h保持不变”的假设下成立;真实流动中的h可能随尺寸、速度和温度变化,不能自动沿用。更不能仅靠质量翻倍就推断τ翻倍而忽略面积。

基准时间常数27.778秒、Bi0.002;提高换热系数百倍后时间常数0.278秒、Bi0.2,几何相似放大两倍后时间常数55.556秒、Bi0.004。
原创合成参数对照卡;卡片面积不代表比例,数字是在其他参数按文中规定保持或缩放时的计算结果。查看原图

怎样审阅一份冷却计算

先在纸上画出物体、真正暴露的表面和环境,再写出参数单位。计算C与hAₛ后先核对τ是否为秒,随后计算Bi,才进入温度预测。对输出至少检查四件事:t=0回到初温;正时间温度保持在初温和恒定环境之间;冷却时温度单调下降;累计散热不超过有限内能差。

实验数据若明显偏离指数曲线,不应只调h把图“拟合得好看”。温度依赖物性、辐射、接触热阻、环境升温、内部温度分层和传感器滞后都可能贡献偏差。只凭一条曲线拟合出的τ,在C未知时也不能唯一识别h;只能得到hAₛ/C这一组合。

自检1:一个时间常数后温度是不是80/e?

不是。衰减的是相对环境的温差。该例应为20+60/e≈42.073°C,不能把摄氏温标的零点当作环境温度。

自检2:为什么k不出现在τ里,却仍要知道k?

τ是在内部近似等温的前提下推导的。k通过Bi=hLc/k决定这个前提是否有依据;舍去空间温差后,k才从单温度方程里消失。

自检3:两倍尺寸是否一定需要八倍冷却时间?

不是。几何相似且h不变时,热容增八倍,但换热面积也增四倍,τ只增两倍。Bi也翻倍,仍应重新核对近似。

自检4:环境从20突然变到50°C,物体温度会立即变50吗?

不会。有限热容要求温度连续,第二段继承第一段末温;变化的是新的驱动温差和温度导数。

参考文献与延伸阅读

  1. MIT Unified Engineering:18.3 Transient Heat Transfer (Convective Cooling or Heating)。官方讲义给出Biot数量级、集总能量方程与时间常数。本文另给合成参数、分段环境与独立算术核验,不把教材模型说成真实材料测试。返回公式C1/返回C2。
  2. MIT Unified Engineering:17.2 Combined Conduction and Convection。官方热阻讲义用于对流热阻与温差驱动的衔接;稳态热阻本身不提供储能随时间的解。
  3. 本站材料路线:稳态导热先区分表面与流体温度;热膨胀再把温度变化接到形变边界。温度不均时不能把单一温升直接代入均温杆模型。