热传导方程

一、直观理解：温度差会“抹平”

热传导方程描述热量在介质中因温度梯度而传递，使高温区变凉、低温区变暖，整体温度场逐渐变平滑。一维最简模型是

$$\frac{\partial T}{\partial t}=\alpha\frac{\partial^2 T}{\partial x^2},$$

其中 $T(x,t)$ 是温度，$\alpha>0$ 是热扩散率（与导热系数、密度、比热有关）。

直观图像：

• 一根细杆，中间烫、两端凉：热流向两端，中峰下降并变宽；
• 若两端维持 $0^\circ$（或固定温度），长时间后整杆趋于边界决定的稳态；
• 若两端绝热，总热量近似守恒，温度趋于空间平均值。

同一数学形式也描述物质扩散、某些期权模型的核心部分等，因此热方程是抛物型 PDE 的标准基线。在建模中，它既是物理问题本身，也是检验 pdepe、有限差分稳定性和边界条件设置的理想对象。

二、适用场景与优缺点

类型	说明
适用场景	杆/墙/土层的温度演化；可类比为污染物在静水中的扩散；需要解析解对照的数值实验；讨论边界（定温、绝热、对流换热）对稳态与时间尺度的影响。
核心优点	物理清晰；一维有丰富解析解（分离变量）；参数 $\alpha$ 含义明确；数值方法理论成熟（显式/隐式/CN）。
局限与缺点	傅里叶定律在极端尺度/极短时间可能失效；忽略对流时不能描述流体携热；变物性、相变、辐射会使方程非线性；高维复杂几何需 FEM 等工具。

使用前应确认：

1. 主导机制是传导而非强制对流；
2. 初边值完整且与物理装置一致；
3. $\alpha$ 的单位与 $x,t$ 单位匹配；
4. 有稳态、能量或解析解用于验证。

三、核心原理：傅里叶定律与能量守恒

1. 推导提纲

一维杆横截面积 $A$，内能与温度成正比。热流密度由傅里叶定律

$$q=-k\frac{\partial T}{\partial x}$$

给出。对小区段写能量守恒并取极限，得

$$\rho c_p T_t = (k T_x)_x + Q,$$

$Q$ 为源。若 $k,\rho,c_p$ 为常数，则 $\alpha=k/(\rho c_p)$，$T_t=\alpha T_{xx}+Q/(\rho c_p)$。

2. 定解条件

• 初值：$T(x,0)=T_0(x)$；
• Dirichlet：端点温度给定；
• Neumann：端点热流给定，绝热为 $T_x=0$；
• Robin：对流换热 $-kT_x=h(T-T_{\infty})$。

3. 解析结构（齐次 Dirichlet 例）

$T(0,t)=T(1,t)=0$，$T(x,0)=\sin(n\pi x)$ 时，

$$T(x,t)=e^{-\alpha n^2\pi^2 t}\sin(n\pi x).$$

衰减率与 $n^2$ 成正比：空间变化越剧烈，抹平越快。一般初值可按正弦级数展开。

4. 时间尺度

特征扩散时间

$$t_c\sim \frac{L^2}{\alpha},$$

$L$ 为长度尺度。这比盲目调参更能指导实验与无量纲化。

四、关键决策

1. 要温度场还是只要平均温度
   空间均匀时可降阶为 ODE；梯度重要时才上 PDE。

2. 边界类型
   与实验装置一致：恒温槽、保温层、空气自然对流分别对应不同 BC。

3. 数值路线
   一维：pdepe 或自写差分；二维复杂域：PDE Toolbox / FEM。

4. 参数
   $\alpha$ 来自材料手册或由升温曲线反演；反演时注意与边界耦合。

5. 验证
   解析模态衰减、绝热时平均温度守恒、网格/时间收敛。

五、MATLAB 示例

%% pdepe 解一维热方程
% T_t = alpha * T_xx, x in [0,1]
% T(0,t)=T(1,t)=0, T(x,0)=sin(pi x)

alpha = 0.1;
m = 0;  % 笛卡尔一维
x = linspace(0, 1, 80);
t = linspace(0, 0.5, 60);

pdefun = @(x, t, u, dudx) deal( ...
    1, ...                 % c
    alpha * dudx, ...      % f（通量）
    0);                    % s 源

icfun = @(x) sin(pi * x);
bcfun = @(xl, ul, xr, ur, t) deal( ...
    ul, 0, ...             % 左：u=0  <=> p=ul, q=0
    ur, 0);                % 右：u=0

sol = pdepe(m, pdefun, icfun, bcfun, x, t);
% sol(j,i) ≈ u(x_i, t_j)

figure;
surf(x, t, sol, "EdgeColor", "none");
xlabel("x"); ylabel("t"); zlabel("T");
title("热传导：pdepe 数值解"); colorbar; view(45, 30);

% 与解析解比较
Tnum = sol(end, :)';
Ttrue = exp(-alpha * pi^2 * t(end)) * sin(pi * x);
fprintf("t=tEnd 最大误差 = %.3e\n", max(abs(Tnum - Ttrue)));

figure;
plot(x, Tnum, "o", x, Ttrue, "-", "LineWidth", 1.2); grid on;
legend("pdepe", "解析解"); title("末端时刻截面");

%% 绝热边界：平均温度应近似守恒
bcInsul = @(xl, ul, xr, ur, t) deal( ...
    0, 1, ...              % 左：p+q*f=0 且 f=alpha*ux，取 q=1,p=0 => ux=0
    0, 1);
ic2 = @(x) exp(-30*(x-0.4).^2);
sol2 = pdepe(m, pdefun, ic2, bcInsul, x, t);
meanT = mean(sol2, 2);     % 均匀网格近似
figure;
plot(t, meanT, "LineWidth", 1.2); grid on;
xlabel("t"); ylabel("空间平均 T");
title("绝热杆：平均温度近似守恒");

%% 有限差分显式对照
nx = 81; xfd = linspace(0,1,nx); dx = xfd(2)-xfd(1);
r = 0.4; dt = r*dx^2/alpha;
tEnd = 0.5; nt = ceil(tEnd/dt); dt = tEnd/nt; r = alpha*dt/dx^2;
u = sin(pi*xfd);
for n = 1:nt
    un = u;
    u(2:end-1) = un(2:end-1) + r*(un(3:end)-2*un(2:end-1)+un(1:end-2));
    u(1)=0; u(end)=0;
end
figure;
plot(xfd, u, "o-", xfd, exp(-alpha*pi^2*tEnd)*sin(pi*xfd), "-");
grid on; legend("FDM", "解析"); title(sprintf("显式差分 r=%.2f", r));

%% 时间尺度估计
L = 1; tc = L^2/alpha;
fprintf("特征扩散时间 L^2/alpha = %.3f\n", tc);

pdepe 边界格式提醒：p + q·f = 0 的形式定义 Dirichlet/Neumann，写 BC 时对照文档，避免把 $q$ 设错导致“看起来像绝热其实不是”。

六、如何解释结果

热传导结果应翻译为工程/物理语言：

1. 方程、$\alpha$、几何尺度与边界的实际含义；
2. 温度场如何随时间抹平，热点寿命是否符合 $L^2/\alpha$；
3. 稳态由谁决定（边界温度还是热源平衡）；
4. 数值误差：与解析模态或加密解的对比；
5. 能量/平均温度在绝热条件下是否守恒；
6. 参数不确定性：$\alpha$ 误差如何影响达到某阈值的时间。

若模拟用于“何时可安全触摸/出炉”，应报告阈值定义、最热点位置轨迹及其对 $\alpha$ 的敏感性。

七、常见误区

1. 单位混乱

$\alpha$ 若为 $\mathrm{m}^2/\mathrm{s}$，却把 $x$ 写成厘米、$t$ 写成分钟，数值会错数量级。

2. 边界与实验不符

实际存在换热却设成绝热，稳态与时间都会偏。

3. 显式差分不查 $r\le 1/2$

出现网格振荡时先查稳定性，再谈物理失稳。

4. 把数值黏性或绘图插值当成额外物理现象

5. 需要对流换热却坚持纯热方程

流体运动主导时应耦合对流项或改用集总参数/CFD 简化。

八、小结

热传导方程的核心思想可以概括为：

在傅里叶定律与能量守恒下，温度场以扩散方式演化，空间曲率（二阶导数）驱动随时间的抹平，速率由热扩散率 $\alpha$ 控制。

它是抛物型问题的第一范例，也是连接物理直觉、解析解与数值方法的桥梁。建模时把边界、量纲和验证放在求解器名称之前：材料与边界正确了，pdepe 或差分才能给出可信的温度史。