4.2 偏微分方程与热扩散 · 从直觉到证明

一步一步建立直观

先预测，再点“演示本步”核对右图。每一步都重设到指定情形，随后可以自由修改参数。

常微分方程追踪有限个状态，热方程则追踪每个空间位置的温度。右图的横轴是位置，纵轴可看作相对于平均温度的偏差，负值表示低于平均。时间推进时，整条温度曲线改变，空间形状本身成为系统状态。

一个温度尖峰比两侧更热，热流向外，峰顶就会下降；谷底则因邻域更热而上升。二阶空间导数把这种邻域不均衡写成局部量。本例的方程是 $u_t=u_{xx}$ ，且区间首尾相接；边界条件不同，适用的基本模式也会改变。

将高频模式的频率加大，而保持观察时间相同。高频部分的曲率更大，所以衰减更快，速度含频率的平方。图形显示的是两个正弦模式的精确演化，不是任意材料与边界的完整仿真；物理单位和扩散系数在这里已经归一化。

一根杆在各位置有不同温度，所以状态是温度分布 $u(x,t)$ 。时间改变时，整条曲线都变。偏导数让我们分别固定位置观察时间变化，或固定时间观察空间变化，两者承担不同物理角色。

为什么热方程出现空间二阶导数，而不是温度本身？一个位置比两侧邻居都热，就会向两侧散热；比邻居冷则吸热。决定变化的是中心相对于邻域平均的偏离。三点模型给出 $u(x-h)-2u(x)+u(x+h)$ ，除以 $h^2$ 后趋向曲率。

在常系数模型中，热通量与负温度梯度成正比；某点温度变化来自流入与流出的差，所以对通量再取一次空间变化，便得到二阶导数。线性温度分布虽然有热流，但各小段流入流出平衡，内部温度可以暂时不变。

端点怎样处理会改变整个问题。保持端点温度固定，与隔热端点无热流，是不同边界条件；同一个初始曲线在不同边界下会演化不同。只给方程而不给初始和边界信息，通常不足以选出具体解。

若做离散模拟，还要检查时间步稳定性。把空间分得更细却不相应调整时间步，可能让数值解振荡甚至爆炸，那不一定是物理现象。接下来我们寻找一些形状简单、能精确追踪的模式，用来解释也用来校验模拟。

与本图相连的定义与证明：4.3 · 正弦模式满足热方程。本节的具体条件与形式化范围见下文。

温度同时依赖位置和时间。热从较热处流向较冷处，局部的“比邻居更高还是更低”决定温度怎样变化。

一维均匀介质模型为 $u_t=\kappa u_{xx}$ ， $\kappa>0$ 。离散空间中， $u_{xx}(x)\approx\frac{u(x-h)-2u(x)+u(x+h)}{h^2}.$ 这个二阶差分是局部平均恢复趋势的代数版本。

例 4.2.1 把概念用于具体计算

邻点温度为 $4,10,6$ ，中心二阶差分分子为 $-10$ ，所以中心倾向降温。若为 $6,2,8$ ，分子为 $10$ ，中心倾向升温。

为什么两个解相加仍满足齐次热方程？

展开解答与检查

时间导数与空间二阶导数都线性。

Lean 对照 4.2.1 命题、证明与适用边界

验证三点二阶差分算子对相加的线性性；连续 PDE 的分析及离散极限未由此命题覆盖。

DifferentialEquations.heat_stencil

theorem heat_stencil (a b c d e f κ : ℝ) :
    κ*((a+d)-2*(b+e)+(c+f)) = κ*(a-2*b+c)+κ*(d-2*e+f) := by ring

原课对照观课问题、视频与图示

观察中心相对邻居的高低，解释为什么温度变化由空间二阶差分决定。

点击后连接 YouTube；播放器不可用时可通过上方原链接观看。