跳到正文
格致开物MATHWIKI

热核

热核(heat kernel)描述单位点源经过一定时间后的扩散分布。在整个欧氏空间、扩散系数为正常数的热方程中,它是一个随时间变宽、峰值降低、总质量保持一的高斯核。核本身与某个具体初始温度的解不同;后者由核与初始数据卷积得到。

设空间维数为正整数 n,扩散系数 κ>0。对 t>0,全空间热核为 Kt(x)=1(4πκt)n/2exp⁡(−|x|24κt). 它满足 ∂tKt=κΔKt。在概率意义下,这是均值零、协方差矩阵 2κtIn 的正态密度;方差中的二不能遗漏。

方程和总质量的逐项检查

对数求导给出 ∂tKt=(−n2t+|x|24κt2)Kt. 对每个空间坐标分别求两次导数再相加,得到 ΔKt=(−n2κt+|x|24κ2t2)Kt. 第二式乘以 κ 恰好等于第一式,验证方程。再作换元 x=4κty,利用高斯积分 ∫ℝne−|y|2dy=πn/2,可得 ∫ℝnKt(x)dx=1. 峰值随时间下降,并不意味着质量消失;下降与宽度增加相互补偿。

一维扩散系数一的热核在时间四分之一、一和四的高斯曲线,越晚越宽且峰越低
三条曲线在整个实轴下的面积均为一,图窗外的尾部仍计入总质量。纵轴是单位点源的密度,不是固定物体的实际温度。

初始点源的正确极限

对任意光滑紧支集测试函数 φ,同样换元得 ∫Kt(x)φ(x)dx=π−n/2∫e−|y|2φ(4κty)dy. 测试函数有界,乘上可积的高斯函数可用控制收敛定理,因此 t↓0 时结果趋于 φ(0)。这就是 Kt→δ0 的分布意义。零以外每个固定点处函数趋零,零点处却发散;点态极限不能替代这项单位质量极限。

若把核在负时间定义为零,得到因果时空分布 Γ(x,t)=𝟏t>0Kt(x),则 (∂t−κΔ)Γ=δ(0,0). 对测试函数从时间 ε 起分部积分,内部项由热方程消去,留下初始边界的核平均;再令 ε↓0,正好得到测试函数在时空原点的值。这把正时间的核与基本解定义联系起来。

高斯初值的完整卷积

在一维取 κ=1、初值 f(x)=e−x2。解为 u(x,t)=14πt∫ℝexp⁡(−(x−y)24t−y2)dy. 指数配方: −(x−y)24t−y2=−1+4t4t(y−x1+4t)2−x21+4t. 第一项的积分是 4πt/(1+4t),所以 u(x,t)=11+4texp⁡(−x21+4t). 当 t=3/4,得到 u(x,t)=e−x2/4/2,故中心值为 1/2,x=2 时为 1/(2e)。整条曲线的积分仍为初值的 π,不是一;卷积保留输入质量,只有单位点源核本身被归一化为一。

初始曲线e的负x平方与时间四分之三后的二分之一乘e的负x平方除四,后者更宽且中心降为二分之一
本图的初始总质量为根号π,两条曲线都保持这项质量。它与上一幅单位质量核图的对象不同。

变式。 改为 κ=2、t=3/8,扩散时间尺度 κt 仍为 3/4,同一初值产生完全相同的剖面。一般式只需将 1+4t 替换成 1+4κt。

连续演化与边界限制

高斯配方还给出半群关系 Ks∗Kt=Ks+t:分两段时间扩散,与一次扩散相同。这对应方差 2κs+2κt 相加,不是把峰值简单相乘。对有界连续初值,核的近似恒等性质恢复初始值,正时间的高斯衰减保证可在卷积中求导。

上述核适用于整个空间。有限区间、吸收或反射边界、曲面几何及变系数方程会产生不同热核;不能把这个高斯表达式不加调整地当作所有边值问题的解。反向热演化还会放大高频误差,正时间的平滑性质不代表能够稳定地从模糊数据恢复过去。

参考文献

Stanford University,Math220B课程讲义,Heat Equation — Math220B course notes,§2.4.2,式(2.21);§2.5.4,式(2.22)—(2.23)。