非稳态导热的物理特征与温度场变化
在实际的工程热物理过程中，绝大多数热传导现象都是随时间变化的。非稳态导热（Transient Heat Conduction），又称瞬态导热，是指物体内部的温度分布不仅是空间坐标的函数，而且是时间的函数，即 $T = T(x, y, z, t)$。
与稳态导热（温度场不随时间变化）不同，非稳态导热涉及能量的累积与释放。当物体受到外部热扰动（如环境温度改变或内部热源波动）时，物体内部会产生温度梯度，热量在物体内部重新分布，直到达到新的热平衡状态。
非稳态导热的物理特征
非稳态导热的核心物理特征在于“热惯性”与“扩散速度”。其关键物理量是热扩散率（Thermal Diffusivity），通常用 $\alpha$ 表示，定义为：
$$\alpha = \frac{k}{\rho c_p}$$
其中，$k$ 为导热系数，$\rho$ 为密度，$c_p$ 为比热容。
1. 热扩散率的物理意义
热扩散率衡量了材料在非稳态过程中传递热量的快慢。
- 高热扩散率：材料能够迅速响应外部温度变化，内部温度场能快速趋于均匀或达到稳态。
- 低热扩散率：材料具有较强的热惯性，温度波在内部传播缓慢，导致物体内部在较长时间内存在显著的温度梯度。
2. 温度场的时空演变
在非稳态导热过程中，温度场的演变通常经历以下阶段：
- 初始响应阶段：热扰动首先作用于边界，边界处产生极大的温度梯度，热量迅速向内部渗透。
- 渗透阶段：热量逐渐向物体深处传递，温度分布曲线随时间平缓化，受影响的区域（热渗透深度）不断扩大。
- 趋稳阶段：随着时间增加，物体内部的温度变化率 $\partial T/\partial t$ 逐渐减小，最终温度场趋向于一个不随时间变化的稳态分布。
非稳态导热的数学模型
非稳态导热的数学描述基于能量守恒定律和傅里叶导热定律。对于一个各向同性的均匀介质，其通用控制方程（热传导方程）为：
$$\rho c_p \frac{\partial T}{\partial t} = k \left( \frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} + \frac{\partial^2 T}{\partial z^2} \right) + q$$
其中 $q$ 为单位体积的内热源强度。若考虑一维简化模型且无内热源，方程可简化为：
$$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2}$$
这是一个典型的抛物线型偏微分方程。为了求解该方程，必须定义两个关键条件：
- 初始条件：$t=0$ 时，物体内部全空间的温度分布 $T(x, 0) = T_i(x)$。
- 边界条件：在 $t > 0$ 时，物体表面在不同时间点的热状态（如给定温度、给定热流密度或对流换热条件）。
典型分析方法：集中参数法
在某些特定条件下，非稳态导热的分析可以极大简化。当物体的内部导热能力远高于其表面的对流换热能力时，可以认为物体内部温度分布在任何时刻都是均匀的。这种方法被称为集中参数法（Lumped Capacitance Method）。
1. 判定准则：毕奥数 (Biot Number)
判断是否可以使用集中参数法的标准是毕奥数 $Bi$：
$$Bi = \frac{h L_c}{k}$$
其中 $h$ 为对流换热系数，$L_c$ 为特征长度（通常为 $\text{体积}/\text{表面积}$），$k$ 为材料导热系数。
- 当 $Bi < 0.1$ 时，内部温度梯度可忽略，适用集中参数法。
- 当 $Bi \ge 0.1$ 时，必须考虑空间温度分布，采用解析法（如分离变量法）或数值法（如有限差分法）。
2. 温度变化示例
假设一个球形金属小球（$Bi < 0.1$）初始温度为 $T_i$，将其放入温度为 $T_\infty$ 的流体中。其温度随时间的变化遵循指数衰减规律：
$$\frac{T(t) - T_\infty}{T_i - T_\infty} = \exp\left( -\frac{h A_s}{\rho V c_p} t \right)$$
该公式表明，物体的温度是以指数形式趋近于环境温度的。时间常数 $\tau = \frac{\rho V c_p}{h A_s}$ 决定了冷却或加热的速度。
总结
非稳态导热研究的是温度在时间和空间上的双重演变。其物理核心在于热扩散率 $\alpha$ 对热量传递速度的控制。在实际应用中，通过分析毕奥数 $Bi$，工程师可以决定是采用简单的集中参数模型进行快速估算，还是采用复杂的偏微分方程模型进行精确分析。理解非稳态导热对于设计散热器、分析淬火工艺以及研究生物组织的热响应具有至关重要的意义。