非稳态导热的无量纲数：傅里叶数与毕渥数
在非稳态（瞬态）导热过程中，物体的温度不仅随空间位置变化，还随时间不断演变。处理这类偏微分方程（PDEs）时，直接求解往往极其复杂。为了简化数学模型并揭示物理本质，热传导学引入了无量纲数。其中，**傅里叶数（Fourier Number, $Fo$）和毕渥数（Biot Number, $Bi$）**是最核心的两个参数。它们分别从“时间尺度”和“热阻竞争”两个维度，决定了我们应当采用何种数学模型来描述传热过程。
傅里叶数 $Fo$ 是描述热量在物体内部扩散快慢的一个无量纲时间参数。它反映了热扩散过程在特定时间尺度下的进展程度。
1. 数学定义
傅里叶数的定义公式为：
$$Fo = \frac{\alpha t}{L^2}$$
其中：
- $\alpha$ 是物体的热扩散率（Thermal Diffusivity），单位为 $\text{m}^2/\text{s}$，定义为 $\alpha = \frac{k}{\rho c_p}$（$k$ 为导热系数，$\rho$ 为密度，$c_p$ 为比热容）。
- $t$ 是经过的时间，单位为 $\text{s}$。
- $L$ 是物体的特征长度（Characteristic Length），单位为 $\text{m}$。
2. 物理意义
从物理本质上看，$Fo$ 是“热扩散时间”与“特征扩散时间”的比值。
- 热扩散率 $\alpha$ 衡量了物体传递热量与储存热量之间能力的平衡。$\alpha$ 越大，热量在物体内部传播得越快。
- $L^2$ 代表了热量扩散所需空间尺度平方的量级。
- 当 $Fo \ll 1$ 时，表示时间尚短，热量仅在物体表面附近发生变化，尚未渗透到物体中心。
- 当 $Fo$ 增大时，表示热量在物体内部的分布趋于均匀，热传导过程进入了更深层次的演变阶段。
毕渥数（Biot Number）：热阻的竞争
如果说傅里叶数关注的是“时间”，那么毕渥数 $Bi$ 关注的则是“空间上的热阻分布”。它描述了物体内部导热阻力与物体表面对流换热阻力之间的相对大小。
1. 数学定义
毕渥数的定义公式为：
$$Bi = \frac{h L_c}{k}$$
其中：
- $h$ 是物体表面的对流换热系数，单位为 $\text{W}/(\text{m}^2 \cdot \text{K})$。
- $k$ 是物体材料的导热系数，单位为 $\text{W}/(\text{m} \cdot \text{K})$。
- $L_c$ 是物体的特征长度，对于复杂形状，通常定义为体积与表面积之比：$L_c = \frac{V}{A}$。
2. 物理意义
毕渥数本质上是两种热阻的比值：
$$Bi = \frac{\text{对流热阻}}{\text{内部导热热阻}} = \frac{1/hA}{L_c/k}$$
- 当 $Bi \ll 1$ 时：意味着物体内部的导热热阻远小于表面的对流热阻。在这种情况下，热量在物体内部的传递速度极快，几乎可以忽略内部温差，物体内部各点的温度在某一时刻可以视为是均匀的。
- 当 $Bi \gg 1$ 时：意味着物体内部的导热热阻远大于表面的对流热阻。此时，热量在物体表面的传递很快，但由于内部导热缓慢，物体内部会存在显著的温度梯度。
集中参数法（Lumped Capacitance Method）
在工程实践中，毕渥数最重要的应用就是决定是否可以使用集中参数法来简化非稳态导热问题。
1. 适用条件
当满足以下判据时，可以使用集中参数法：
$$Bi = \frac{h L_c}{k} < 0.1$$
在这种情况下，我们可以将物体看作一个“集中热容”的整体，忽略物体内部的温度梯度，仅考虑物体整体温度 $T(t)$ 随时间的变化。
2. 数学模型
在集中参数法下，非稳态导热的偏微分方程简化为一阶常微分方程：
$$\rho V c_p \frac{dT}{dt} = -h A (T - T_\infty)$$
通过积分求解，得到温度随时间变化的解析解：
$$\frac{T(t) - T_\infty}{T_i - T_\infty} = \exp\left( -\frac{h A}{\rho V c_p} t \right)$$
为了与傅里叶数建立联系，我们可以将其改写为：
$$\frac{T(t) - T_\infty}{T_i - T_\infty} = \exp\left( -Bi \cdot Fo \cdot \frac{L^2}{L_c^2} \right)$$
（注：在球体或平板等规则几何体中，该指数项通常简化为 $Bi \cdot Fo$ 的某种组合形式）。
综合应用示例
案例描述：
假设一个半径为 $R = 0.05,\text{m}$ 的铜球（导热系数 $k = 385,\text{W/m}\cdot\text{K}$，密度 $\rho = 8960,\text{kg/m}^3$，比热 $c_p = 385,\text{J/kg}\cdot\text{K}$）被置于空气中进行冷却，空气的对流换热系数 $h = 50,\text{W/m}^2\cdot\text{K}$。请判断是否可以使用集中参数法，并计算冷却 100 秒后的温度变化情况（假设初始温度 $T_i = 100^\circ\text{C}$，环境温度 $T_\infty = 20^\circ\text{C}$）。
第一步：计算特征长度 $L_c$
对于球体：
$$L_c = \frac{V}{A} = \frac{\frac{4}{3}\pi R^3}{4\pi R^2} = \frac{R}{3} = \frac{0.05}{3} \approx 0.0167,\text{m}$$
第二步：计算毕渥数 $Bi$
$$Bi = \frac{h L_c}{k} = \frac{50 \times 0.0167}{385} \approx 0.00217$$
由于 $Bi = 0.00217 < 0.1$，满足集中参数法的适用条件。我们可以认为铜球内部温度是均匀的。
第三步：计算傅里叶数 $Fo$（或直接使用时间项）
首先计算热扩散率 $\alpha$：
$$\alpha = \frac{k}{\rho c_p} = \frac{385}{8960 \times 385} \approx 1.116 \times 10^{-4},\text{m}^2/\text{s}$$
计算 $t = 100,\text{s}$ 时的 $Fo$：
$$Fo = \frac{\alpha t}{R^2} = \frac{1.116 \times 10^{-4} \times 100}{(0.05)^2} \approx 4.464$$
第四步：求解温度 $T(t)$
使用集中参数公式：
$$\frac{T(100) - 20}{100 - 20} = \exp\left( -\frac{h A}{\rho V c_p} \times 100 \right)$$
注意到 $\frac{h A}{\rho V c_p} = \frac{h}{\rho c_p L_c}$，代入数值：
$$\text{指数项} = -\frac{50}{8960 \times 385 \times 0.0167} \times 100 \approx -0.0095$$
$$\frac{T(100) - 20}{80} = e^{-0.0095} \approx 0.9905$$
$$T(100) = 20 + 80 \times 0.9905 \approx 99.24^\circ\text{C}$$
结论：
由于铜的导热能力极强，毕渥数极小，冷却过程在球体内表现为温度均匀下降。经过 100 秒后，球体温度仅从 $100^\circ\text{C}$ 降至约 $99.24^\circ\text{C}$。
总结
在非稳态导热分析中，这两个无量纲数构成了决策的基石：
- 毕渥数 $Bi$ 决定了模型的复杂程度：它告诉我们是可以用简单的常微分方程（集中参数法）解决，还是必须使用复杂的偏微分方程（或 Heisler 图表/数值模拟）。
- 傅里叶数 $Fo$ 决定了过程的演进阶段：它将物理时间转化为热扩散的进度，是分析瞬态响应随时间演变的核心变量。
理解并熟练运用这两个参数，是掌握热传导学从理论推导走向工程设计的关键一步。