有内热源稳态导热的能量守恒表达

在热传导学研究中，理解能量如何在物体内部流动以及如何产生是建立数学模型的基础。在许多实际工程场景中，热量不仅通过物体表面进行交换，还会由于物体内部的物理或化学过程而产生，例如电加热元件中的焦耳热、核反应堆中的核能释放，以及化学反应器中的化学能转化。

本文将深入探讨在稳态条件下，含有内热源的物体如何通过能量守恒定律建立其导热的数学表达形式。
根据热力学第一定律，对于一个封闭的控制体（Control Volume, CV），能量的变化率等于净流入能量与内部能量产生之和。在热传导问题中，能量主要以热流的形式进行传递。

对于一个处于稳态（Steady-state）的系统，意味着系统内任意一点的温度不随时间变化，即 $\frac{\partial T}{\partial t} = 0$。此时，控制体内的能量储存项为零。能量守恒的物理描述可以简化为：

$$\text{净流入的热流} + \text{内部产生的热量} = 0$$

或者更严谨地表达为：
$$\dot{E}{in} - \dot{E}{out} + \dot{E}_{gen} = 0$$

其中：

• $\dot{E}_{in}$ 为通过控制体边界流入的热流率。
• $\dot{E}_{out}$ 为通过控制体边界流出的热流率。
• $\dot{E}_{gen}$ 为控制体内单位时间内产生的热量。

导热方程的数学推导

为了将上述物理描述转化为可求解的偏微分方程，我们需要引入微元分析法和傅里叶导热定律（Fourier's Law）。

1. 微元控制体的能量平衡

考虑一个微小的控制体 $dV$。在稳态条件下，该微元内的能量平衡方程为：
$$-\nabla \cdot \mathbf{q} + \dot{q} = 0$$

这里：

• $\mathbf{q}$ 是热流密度向量（Heat flux vector），单位为 $\text{W/m}^2$。
• $\nabla \cdot \mathbf{q}$ 是热流密度的散度，表示单位时间内单位体积内流出的净热量。
• $\dot{q}$ 是单位体积内的热源强度（Volumetric heat generation rate），单位为 $\text{W/m}^3$。

2. 引入傅里叶定律

根据傅里叶导热定律，热流密度与温度梯度成正比：
$$\mathbf{q} = -k \nabla T$$

其中 $k$ 是材料的热导率（Thermal conductivity）。将此式代入能量平衡方程中，得到通用的导热微分方程：
$$\nabla \cdot (k \nabla T) + \dot{q} = 0$$

如果假设材料是各向同性的且热导率 $k$ 为常数，方程可以进一步简化为：
$$\nabla^2 T + \frac{\dot{q}}{k} = 0$$
这在数学上被称为泊松方程（Poisson's Equation），而如果没有内热源（即 $\dot{q}=0$），则退化为拉普拉斯方程（Laplace's Equation）。

不同坐标系下的表达形式

由于工程问题的几何形状各异，我们需要在不同的坐标系下展开上述算子 $\nabla^2$。

一维直角坐标系 (Cartesian Coordinates)

对于仅在 $x$ 方向发生导热的情况：
$$\frac{d}{dx} \left( k \frac{dT}{dx} \right) + \dot{q} = 0$$
若 $k$ 为常数：
$$\frac{d^2T}{dx^2} + \frac{\dot{q}}{k} = 0$$

圆柱坐标系 (Cylindrical Coordinates)

对于圆柱形物体（如电缆、加热棒），假设热量仅沿半径 $r$ 方向传递：
$$\frac{1}{r} \frac{d}{dr} \left( r k \frac{dT}{dr} \right) + \dot{q} = 0$$
若 $k$ 为常数：
$$\frac{d^2T}{dr^2} + \frac{1}{r} \frac{dT}{dr} + \frac{\dot{q}}{k} = 0$$

球坐标系 (Spherical Coordinates)

对于球形物体（如核燃料球），假设热量仅沿半径 $r$ 方向传递：
$$\frac{1}{r^2} \frac{d}{dr} \left( r^2 k \frac{dT}{dr} \right) + \dot{q} = 0$$
若 $k$ 为常数：
$$\frac{d^2T}{dr^2} + \frac{2}{r} \frac{dT}{dr} + \frac{\dot{q}}{k} = 0$$

典型实例分析：圆柱形均匀内热源

为了更好地理解这些方程的应用，我们来看一个经典的工程案例：一个半径为 $R$、热导率为 $k$ 的实心圆柱体，其内部均匀产生热量 $\dot{q}$，且圆柱表面温度维持在恒定的 $T_s$。

1. 建立方程：
使用圆柱坐标系的稳态方程：
$$\frac{d}{dr} \left( r \frac{dT}{dr} \right) = -\frac{\dot{q}r}{k}$$

2. 积分求解：
对 $r$ 进行第一次积分：
$$r \frac{dT}{dr} = -\frac{\dot{q}r^2}{2k} + C_1$$
整理得：
$$\frac{dT}{dr} = -\frac{\dot{q}r}{2k} + \frac{C_1}{r}$$

进行第二次积分：
$$T(r) = -\frac{\dot{q}r^2}{4k} + C_1 \ln(r) + C_2$$

3. 应用边界条件：

• 边界条件一（对称性）： 在圆柱中心 $r=0$ 处，由于对称性，温度梯度必须为零，即 $\left. \frac{dT}{dr} \right|_{r=0} = 0$。这要求 $C_1 = 0$（否则 $\ln(0)$ 无意义且梯度趋于无穷）。
• 边界条件二（表面温度）： 在表面 $r=R$ 处，$T(R) = T_s$。
  $$T_s = -\frac{\dot{q}R^2}{4k} + C_2 \implies C_2 = T_s + \frac{\dot{q}R^2}{4k}$$

4. 最终温度分布公式：
$$T(r) = T_s + \frac{\dot{q}}{4k}(R^2 - r^2)$$

通过该公式可以清晰地看到，在均匀内热源的作用下，圆柱体内部的温度分布呈抛物线型，最高温度出现在中心点 $r=0$。

总结

有内热源稳态导热的能量守恒表达是热传导分析的核心。通过将物理上的能量平衡转化为数学上的泊松方程，工程师可以根据物体的几何形状选择合适的坐标系，并结合边界条件求解出精确的温度分布。掌握这一模型对于设计散热系统、核反应堆安全评估以及电子元件的热管理具有至关重要的意义。