1. 热方程物理背景与数学表述

在进入有限差分数值求解之前，必须先弄清楚我们究竟在解什么。热方程描述的是温度在空间和时间上的演化规律，它并非凭空写出的偏微分方程，而是能量守恒与傅里叶导热定律结合的必然结果。本节从物理直觉出发，逐步推导一维热方程，并明确温度场、热扩散系数、边界条件与初始条件各自的物理含义，为后续离散化建立清晰的问题原型。

从能量守恒到控制方程

考虑一根截面积为 A、密度为 ρ、比热容为 c 的均匀细杆，沿轴向取一个厚度为 Δx 的微元。设温度场为 T(x, t)，即位置 x 处、时刻 t 的温度。微元内热能的增加率等于净流入的热量（无内热源时）：

微元热能为 ρcAΔx·∂T/∂t；左侧流入热流为 q(x)A，右侧流出为 q(x+Δx)A，净流入为 −[q(x+Δx)−q(x)]A。令 Δx→0，得到能量守恒的微分形式：

ρc ∂T/∂t = −∂q/∂x

这里 q 是热流密度（单位 W/m²），负号表示热流方向与温度梯度方向相反——热量总是从高温流向低温。

引入傅里叶定律

傅里叶导热定律给出热流密度与温度梯度的关系：

q = −k ∂T/∂x

其中 k 为导热系数（W/(m·K)）。将其代入能量守恒式，得到：

ρc ∂T/∂t = ∂/∂x(k ∂T/∂x)

若材料均匀，k 为常数，两边同除以 ρc，定义热扩散系数 α = k/(ρc)，单位 m²/s，便得到标准的一维热方程：

∂T/∂t = α ∂²T/∂x²

这个方程说明：某点温度随时间的变化率，正比于该点温度分布的二阶空间导数（即温度曲线的“弯曲程度”）。温度分布越“尖”，局部变化越快；温度分布越平缓，变化越慢。这正是热传导“抹平”温度差异的数学表达。

温度场与热扩散系数的物理含义

温度场 T(x, t) 是本节的核心未知量，它是一个二元函数：在每一时刻 t，温度沿杆长 x 形成一条曲线；随着时间推进，这条曲线不断演化。数值求解的目标，就是在离散的网格点 (x_i, t_n) 上求出 T_iⁿ 的近似值。

热扩散系数 α 是唯一出现在方程中的材料参数，它把导热能力（k）与储热能力（ρc）综合在一起。直观地说，α 越大，热量在材料中“跑”得越快，温度分布被抹平得越迅速。以常见材料为例：纯铜的 α 约为 1.1×10⁻⁴ m²/s，而混凝土约为 0.5×10⁻⁶ m²/s，相差两个数量级以上。这意味着同样尺寸下，铜杆中的温度扰动几乎瞬间扩散，而混凝土墙体中的温度变化则缓慢得多。这个参数直接决定了数值求解时时间步长的选取——α 越大，显式格式对 Δt 的限制越严格。

初始条件与边界条件

热方程本身只描述“如何演化”，并不指定“从哪里出发”和“端点如何”。要构成一个适定的初边值问题，必须补充两类条件。

初始条件给出 t=0 时刻整根杆的温度分布：T(x, 0) = T₀(x)。例如一根初始处于均匀室温 20 °C 的金属杆，可写为 T₀(x) = 20。若左半段被加热到 80 °C、右半段保持 20 °C，则初始条件是一条阶跃型分布。初始条件的光滑程度会影响数值解的早期行为，阶跃处往往出现较陡的梯度。

边界条件描述两端 x=0 与 x=L 处的换热情况，常见三类：

• 第一类（Dirichlet，定温）：直接给定端点温度，如 T(0, t) = 100 °C。适用于端点与恒温源良好接触的情形。
• 第二类（Neumann，定热流）：给定端点热流密度，如 −k ∂T/∂x|_x=0 = q₀。绝热边界即 q₀ = 0，对应 ∂T/∂x = 0。
• 第三类（Robin，对流）：端点与流体对流换热，如 −k ∂T/∂x|_x=0 = h[T∞ − T(0, t)]，其中 h 为对流换热系数，T∞ 为环境温度。

边界条件的类型直接改变离散后代数方程组的构造方式：定温边界把端点值作为已知量移入右端项，而定热流与对流边界则需要用虚拟节点或单侧差分来近似导数。这是后续编程实现中最容易出错的地方之一，必须在问题原型阶段就明确标注。

一个可复现的教学示例

为便于后续章节对照，这里设定一个具体的问题原型（参数为教学假设，非实测）：一根长 L = 1 m 的均匀铝杆，热扩散系数取 α = 9.7×10⁻⁵ m²/s。初始温度全场 20 °C；t > 0 后左端 x=0 突然接触 100 °C 恒温源（第一类边界），右端 x=L 绝热（第二类边界，热流为零）。

这个组合的物理图像很清晰：热量从左端注入，向右扩散，由于右端绝热，热量不断累积，最终整根杆趋于 100 °C 的稳态。数值求解时，我们关心的是不同时刻的温度曲线形状、左端附近梯度的大小，以及达到近似稳态所需的时间量级。粗略估计，扩散时间尺度约为 L²/α ≈ 1/(9.7×10⁻⁵) ≈ 1.03×10⁴ s，即约 2.9 小时，这为选择模拟总时长提供了参考。

为离散化做好准备

至此，我们已把物理问题转化为一个明确的数学原型：在区间 [0, L] 上求解 ∂T/∂t = α ∂²T/∂x²，附带初始分布与两端边界条件。后续的有限差分工作，本质上就是用网格上的差商代替 ∂T/∂t 与 ∂²T/∂x²，把连续方程变成一组可逐点推进的代数递推式。理解本节建立的物理与数学原型，是判断差分格式是否合理、边界处理是否正确、数值结果是否可信的根本依据。