非稳态导热微分方程的建立
在热传导学的研究中，根据温度场是否随时间变化，可以将导热过程分为稳态导热（Steady-state conduction）和非稳态导热（Unsteady-state conduction，又称瞬态导热）。稳态导热研究的是系统达到热平衡后的状态，此时空间中各点的温度仅是位置的函数；而实际工程应用中，如金属淬火、电子器件加热、建筑材料受热等过程，温度场往往随时间剧烈波动。
非稳态导热问题的核心在于建立描述温度 $T$ 随空间坐标 $\mathbf{r}$ 和时间 $t$ 共同变化的偏微分方程。本文将从能量守恒定律和傅里叶导热定律出发，详细推导非稳态导热的通用微分方程。
物理基础：两大核心定律
建立非稳态导热方程需要基于两个基本的物理原理：
1. 能量守恒定律（热力学第一定律）
对于一个微小的控制体（Control Volume, $dV$），其能量的变化率必须等于通过边界进入的净热流率加上控制体内产生的内热源。用数学表达式表示为：
$$\text{能量积累率} = \text{净流入热流率} + \text{内热源产生率}$$
2. 傅里叶导热定律（Fourier's Law）
傅里叶定律描述了热量传递的驱动力与热流密度之间的关系。在各向同性介质中，热流密度矢量 $\mathbf{q}$ 与温度梯度 $\nabla T$ 成正比：
$$\mathbf{q} = -k \nabla T$$
其中，$k$ 为材料的热导率（Thermal Conductivity），$\nabla$ 为梯度算子。负号表示热量总是从高温区向低温区传递。
非稳态导热微分方程的数学推导
为了获得微分形式的方程，我们选取一个微小的控制体 $dV$。
第一步：能量积累项
假设材料的密度为 $\rho$，比热容为 $c_p$。控制体内单位体积的内能为 $\rho e$（在固体中，内能变化主要体现为显热，即 $\rho c_p T$）。那么，控制体 $dV$ 内能量随时间的变化率为：
$$\frac{\partial}{\partial t} (\rho c_p T) dV$$
第二步：热流净流入项
根据高斯散度定理（Gauss's Divergence Theorem），通过控制体表面 $A$ 的净热流通过体积积分可以转化为体积上的散度积分：
$$\oint_A \mathbf{q} \cdot \mathbf{n} , dA = \int_V (\nabla \cdot \mathbf{q}) , dV$$
由于我们要计算的是“流入”的热量，根据定义，净流入量为 $-\int_V (\nabla \cdot \mathbf{q}) , dV$。
第三步：内热源项
若控制体内存在单位体积的热源（如电加热、化学反应或核反应产生的热量），记为 $\dot{q}$，则总热产生率为：
$$\int_V \dot{q} , dV$$
第四步：合并方程
将上述三项代入能量守恒方程：
$$\int_V \rho c_p \frac{\partial T}{\partial t} , dV = -\int_V (\nabla \cdot \mathbf{q}) , dV + \int_V \dot{q} , dV$$
由于该方程对任意微小的控制体 $dV$ 都成立，我们可以去掉积分符号，得到微分形式：
$$\rho c_p \frac{\partial T}{\partial t} = -\nabla \cdot \mathbf{q} + \dot{q}$$
第五步：引入傅里叶定律
将 $\mathbf{q} = -k \nabla T$ 代入上式，得到非稳态导热的通用微分方程：
$$\rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + \dot{q}$$
方程的简化与特殊形式
在实际工程计算中，为了简化计算，通常会对参数进行假设。
1. 各向同性且热导率为常数的情况
如果材料是各向同性的，且热导率 $k$ 不随温度变化，则 $\nabla \cdot (k \nabla T) = k \nabla^2 T$（其中 $\nabla^2$ 是拉普拉斯算子）。方程变为：
$$\frac{\partial T}{\partial t} = \frac{k}{\rho c_p} \nabla^2 T + \frac{\dot{q}}{\rho c_p}$$
这里引入了一个关键的热物理参数——热扩散率（Thermal Diffusivity） $\alpha$：
$$\alpha = \frac{k}{\rho c_p}$$
热扩散率反映了材料通过热传导传递热能的能力。$\alpha$ 越大，材料温度响应越快。简化后的方程为：
$$\frac{\partial T}{\partial t} = \alpha \nabla^2 T + \frac{\dot{q}}{\rho c_p}$$
2. 一维导热情况
在许多工程问题（如薄板加热）中，可以忽略侧向热损失，仅考虑沿 $x$ 轴方向的导热。此时方程退化为一维偏微分方程：
$$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2} + \frac{\dot{q}}{\rho c_p}$$
实际应用示例：一维均匀棒的瞬态加热
问题描述：
假设有一根长度为 $L$ 的均匀金属棒，其热导率 $k$、密度 $\rho$ 和比热 $c_p$ 均为常数。初始时刻（$t=0$），整根棒的温度均匀为 $T_0$。现在将棒的两端保持在恒定温度 $T_1$ 和 $T_2$。求温度场 $T(x, t)$ 的分布。
数学模型建立：
- 控制方程：由于没有内热源（$\dot{q}=0$），方程为：
$$\frac{\partial T}{\partial t} = \alpha \frac{\partial^2 T}{\partial x^2}$$ - 初始条件 (Initial Condition)：
$$T(x, 0) = T_0 \quad (0 \le x \le L)$$ - 边界条件 (Boundary Conditions)：
$$T(0, t) = T_1 \quad (t > 0)$$
$$T(L, t) = T_2 \quad (t > 0)$$
求解思路：
这是一个典型的非齐次边界条件的非稳态问题。通常的处理方法是将其分解为两部分：一个稳态解 $T_s(x)$（满足边界条件但与时间无关）和一个瞬态解 $T_t(x, t)$（满足齐次边界条件）。通过分离变量法，可以求得该问题的解析解。
总结
非稳态导热微分方程是热传导学研究的核心。通过结合能量守恒定律和傅里叶定律，我们得到了描述温度随时间和空间演变的偏微分方程。在建立数学模型时，明确以下要素至关重要：
- 控制方程：取决于材料属性（$\alpha$）和是否存在内热源（$\dot{q}$）。
- 初始条件：描述 $t=0$ 时刻系统的状态。
- 边界条件：描述系统与外界环境的能量交换方式（如恒温、恒热流或对流换热）。
掌握这一方程的建立过程，是进行复杂热力学仿真和工程热设计的基础。