Skip to content

定理Theorem

热方程的 Duhamel 原理

Duhamel principle for the heat equation

把各时刻注入的空间热源先演化到共同终时再积分,构造受迫热方程的解,并区分经典解与积分意义的 mild 解。

形式陈述 ​

一个热源持续工作,每一小份热量从不同时间开始扩散。如何把它们叠到同一个观察时刻?

令 S(t)f=Γt∗f 是整空间的热核演化。考虑

ut−Δu=F(t,x),u(0,x)=f(x).

Duhamel 公式为

u(t)=S(t)f+∫0tS(t−s)F(s)ds.

为保证下面逐项经典求导,一组方便的充分条件是:f∈Cb2(Rd),F:[0,T]→Cb2(Rd) 连续;这里 Cb2 表示函数和至多二阶导数连续且有界,并用它们的上确界之和度量。公式此时产生连续到初始时刻的经典解,正时间满足方程。

若只要求 f∈Lp、1≤p<∞,且 F 作为 Lp 值函数强可测,并满足 ∫0T‖F(s)‖pds<∞,可先把右侧定义为函数空间中的积分解,通常称 mild 解。它不凭公式外观就自动具有逐点的 ut 与 Δu;升级到经典解需要相应正则性。

直觉

时刻 s 注入的热量是 F(s)ds。到时刻 t,这份热量已经扩散了 t−s,所以贡献为 S(t−s)F(s)ds。积分将所有注入时刻的贡献加起来,再叠上最初已有热量的自由演化。

因此演化算子里的时间必须是 t−s。若写成 S(s),就让越晚注入的热量扩散得越久,颠倒了过程的时间方向。

例子与边界

共振式累积不是指数发散 ​

一维中取零初值,热源

F(s,x)=e−scos⁡x.

因为 ∂xxcos⁡x=−cos⁡x,热半群使 S(t−s)cos⁡x=e−(t−s)cos⁡x。于是

u(t,x)=∫0te−(t−s)e−scos⁡xds=te−tcos⁡x.

直接求导:

ut=(1−t)e−tcos⁡x,uxx=−te−tcos⁡x,

两者相减恰为 e−tcos⁡x。因源项衰减率与该模态的自然衰减率相同,各时刻到达终时的贡献相同,积分留下一个 t 因子;整体仍随长时间趋零。

如果源改成恒定的 cos⁡x,结果是 (1−e−t)cos⁡x,而不是 tcos⁡x。忽略注入之后的扩散,会得到错误的长期行为。

局部热源可以从两次积分读出传播 ​

若 F(s,y)=q(s)ψ(y),则

u(x,t)=∫0tq(s)[∫RdΓt−s(x−y)ψ(y)dy]ds.

方括号先在空间传播,外层再按时间累加。若 q,ψ≥0,每份贡献都非负;若 q 在某段时间为零,只是没有新注入,先前的热量仍继续扩散。

推论与应用

方程右端来自积分的上限 ​

记 v(t)=∫0tS(t−s)F(s)ds。在已声明的光滑有界条件下,可以把空间二阶导数移到 F(s) 上;‖ΔF(s)‖∞ 在有限时间段有界,因而控制 s↑t 附近的求导。

用微积分基本定理对移动上限求导,

vt(t)=F(t)+∫0tΔS(t−s)F(s)ds=F(t)+Δv(t).

上限项是 S(0)F(t)=F(t)。若忘了它,就会错误地把受迫解算成齐次解。初值也需要检查:

‖v(t)‖∞≤∫0t‖F(s)‖∞ds⟶0.

第一项 S(t)f 恢复原初值,所以二者合成完整解。

一个有用的稳定性估计 ​

热半群是 Lp 收缩,再用Minkowski 积分不等式控制时间叠加。在前述有限 p 的积分解条件下,或对有界经典解取 p=∞,两组初值与源项 (f,F)、(g,G) 的解满足

‖u(t)−v(t)‖p≤‖f−g‖p+∫0t‖F(s)−G(s)‖pds.

这份误差预算将初值误差和累计注入误差分开。它不承诺误差随时间指数增长,也不要求把源项冻结在某个时刻。

有边界时必须换整套演化算子 ​

在有限区域中,若 SD(t) 是零 Dirichlet 热演化,那么同一公式

u(t)=SD(t)f+∫0tSD(t−s)F(s)ds

自然保持零边界。把自由空间 S(t) 直接代入,一般不能满足墙上的条件。非零且随时间变化的边界数据,则可先选一个满足边界的光滑延拓,减去它,把其时间导数和 Laplacian 合并到新源项中,再用零边界演化。

线性常微分方程组的常数变易公式是同一个叠加思想的有限维版本。PDE 中演化的状态是一整个函数,Laplacian 通常是无界算子;这就是为什么写出积分公式后,还要区分 mild 解与具有经典导数的解。

参考资料
  • Stanford Math 220B,The Heat Equation,§2.7,尤其式 (2.25)–(2.28):受迫热方程的 Duhamel 构造、上限项与初值核验。
  • John K. Hunter,Notes on Partial Differential Equations, Chapter 5,§5.4.5,Theorem 5.45(有界生成元的经典公式)、Definition 5.47(强连续半群的 strong、mild 解):无界生成元时不能仅由公式断言经典可微性。
关系图谱7 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系