1. 课程导入与热扩散物理图像

在工程换热问题中，我们习惯先算“稳态”：给定壁面温度或热流，等温度场不再随时间变化，再求温度分布与热流量。但真实设备很少一开机就处在稳态——发动机冷启动、电子芯片上电瞬间、金属淬火入水、电池快充升温，这些过程里温度场随时间剧烈演化，热量的“扩散速度”和“边界换热强度”共同决定了局部温度峰值、达到安全温度所需时间，以及热应力的大小。本课程要解决的，正是这类瞬态热扩散问题：温度场如何随时间推进，边界上的对流、辐射、接触换热又如何作为“外部约束”参与并塑造这一过程。

一、从“热量守恒”到“热扩散方程”

瞬态导热的物理起点仍是能量守恒。取一个微元体，其内能变化率等于净导入的导热热流加上内热源生成热。对常物性、无内热源的一维情况，可写成：

ρ·c_p·∂T/∂t = ∂/∂x( k·∂T/∂x )

两边同除以 k，定义热扩散率 α = k/(ρ·c_p)（单位 m²/s），得到课程的核心控制方程：

∂T/∂t = α·∇²T        （三维直角坐标下 ∇²T = ∂²T/∂x² + ∂²T/∂y² + ∂²T/∂z²）

这个方程告诉我们两件事。第一，温度场的时间变化率正比于 α：α 越大，材料“抹平”温度差异的速度越快。铜的 α 约为 1.1×10⁻⁴ m²/s，不锈钢约为 4×10⁻⁶ m²/s，两者相差近 30 倍，同样尺寸下铜件几秒内就接近均匀，钢件却要几分钟。第二，方程右侧是空间二阶导数，说明驱动扩散的不是温度本身，而是温度的曲率——温度分布越“弯折”，局部升温或降温越快；温度呈线性分布时，内部各点温度暂时不变。

需要强调符号约定：本课程统一取外法线方向为正，即边界热流 q″ 沿外法线流出物体为正，流入为负。后面所有边界条件、集总参数能量平衡都遵循这一约定，图示方向也与此一致。

二、边界换热在瞬态过程中的角色

控制方程只描述内部扩散，真正决定“温度能升多高、降多快”的是边界条件。常见三类：第一类给定壁温 T_s(t)；第二类给定壁面热流 q″(t)；第三类对流换热 −k·∂T/∂n = h·(T_s − T_∞)，其中 h 为对流换热系数。瞬态问题中，边界条件往往本身随时间变化，例如淬火初期蒸汽膜阶段 h 很低，随后核态沸腾使 h 骤升，这种 h(t) 的阶段性变化会直接反映在表面温度曲线上。

边界换热与内部扩散的竞争，可用毕渥数 Bi = h·L_c/k 衡量（L_c 为特征长度，如平板半厚）。Bi 很小（约 Bi < 0.1）时，内部导热远快于表面换热，物体内部温度近似均匀，可用集总参数法把整个物体当作一个节点：ρ·c_p·V·dT/dt = −h·A·(T − T_∞)，解得指数衰减 T(t) − T_∞ = (T_0 − T_∞)·exp(−t/τ)，时间常数 τ = ρ·c_p·V/(h·A)。Bi 较大时，表面与内部出现明显温差，必须求解分布参数方程。这一判据贯穿全课程，是选择计算路线的第一道分岔。

三、两个无量纲数：Fo 与 Bi

把方程无量纲化，会自然出现两个数。傅里叶数 Fo = α·t/L_c² 表示“已过去的时间相对于热量扩散穿过特征长度所需时间”，Fo 小说明过程刚开始、温度扰动还集中在表面附近；Fo 大则趋于稳态。毕渥数 Bi = h·L_c/k 表示“表面换热能力相对于内部导热能力”。二者组合决定了瞬态温度场的形态：Bi 大而 Fo 小时，表面温度骤变、内部几乎未动，形成陡峭的表层梯度；Bi 小而 Fo 大时，整体温度同步缓慢趋近环境温度。

举个教学示例（数值为合理假设，非实测）：一块厚 20 mm 的钢板（k = 45 W/(m·K)，ρ = 7800 kg/m³，c_p = 470 J/(kg·K)，α ≈ 1.23×10⁻⁵ m²/s），初始 20 ℃，突然置于 200 ℃ 热风中，h = 50 W/(m²·K)。取 L_c = 0.01 m，Bi = 50×0.01/45 ≈ 0.011 < 0.1，可用集总参数：τ = 7800×470×0.01/(50) ≈ 733 s，约 12 分钟后钢板温度达到环境温度的 63%。若把 h 提高到 2000 W/(m²·K)（如水淬），Bi ≈ 0.44，集总假设失效，必须考虑厚度方向温差，表面可能在几十秒内接近水温而中心仍较热——这正是淬火开裂风险的来源。

四、课程学习路径与工程应用

本课程按“物理图像 → 数学工具 → 边界处理 → 数值实现 → 工程判读”推进：先建立扩散与边界耦合的直觉，再掌握集总参数与一维分布解的解析方法，随后处理 h(t)、辐射、接触热阻等真实边界，最后用有限差分/有限元离散求解并用 Fo、Bi 与温度曲线判读结果。工程上，这套方法直接服务于电子器件结温预估与散热设计、金属热处理冷却曲线制定、电池快充热管理与热失控预警、建筑围护结构动态热负荷计算，以及传感器保护套管的热响应时间校核。后续每一节都将围绕“给定边界，如何算出温度随时间的变化，并判断是否满足工程约束”展开。