波方程能量估计 | Codeboy's Blog

显式解公式告诉我们某个点读取哪些数据；能量法不必先知道解的形状。它选一块按波速收缩的观察区域，证明能量不能从外侧赶进来。于是区域初始没有差异、内部也没有差异源时，后来的观察结果必相同。

形式陈述

收缩球中的有源估计

设c>0为常数， $u\in C^2$ 在闭锥

\mathcal C=\{(t,x):0\le t\le T,\ |x-x_0|\le R-ct\}, \qquad R\ge cT

的邻域内满足波方程 $u_{tt}-c^2\Delta u=F$ 。取 $B_t=B(x_0,R-ct)$ ，先在半径为正时定义

E(t)=\frac12\int_{B_t}(u_t^2+c^2|\nabla u|^2)\,dx, \qquad Q(t)=\sqrt{2E(t)}. \tag{1}

这里Q是 $(u_t,c\nabla u)$ 的L2范数，不是u的点值或位移范数。对0≤t<T，有

\boxed{Q(t)\le Q(0)+\int_0^t\|F(s,\cdot)\|_{L^2(B_s)}ds.} \tag{2}

在R=cT时，最后一个球缩成顶点，式(2)通过极限理解；顶点值的结论来自整个锥内的控制，而非“最后一刻球体积为零”。n=1时球就是区间，下面球面项变为左右两个端点项。

依赖域与唯一性

若两个C2解在锥内源项相同，且在底球B0上初位移和初速度都相同，则两个解在整个闭锥相同，包括R=cT时的顶点。特别地，整空间常速波方程的经典Cauchy解至多一个；这项局部结论不要求解在无穷远有界或具有有限总能量。

直觉

局部守恒先保留通量方向

将方程乘以u_t并用乘积法则：

\partial_t e-\operatorname{div}(c^2u_t\nabla u)=u_tF, \qquad e=\tfrac12(u_t^2+c^2|\nabla u|^2). \tag{3}

因此向外的能量流是 $-c^2u_t\nabla u$ 。式(3)只是一条局部恒等式，还没有声称任何固定区域中的能量守恒；区域边界可以有流入或流出。

对球半径b(t)=R−ct，在极坐标写积分后对上限求导，得到

E'(t)=\int_{B_t}\partial_t e\,dx-c\int_{\partial B_t}e\,dS.

第二项来自观察区域自身缩小。用散度定理处理式(3)，

E'(t)=\int_{B_t}u_tF\,dx +\int_{\partial B_t}(c^2u_t\partial_\nu u-ce)\,dS. \tag{4}

把梯度分成法向和切向两部分，球面被积式恰为

c^2u_t\partial_\nu u-ce =-\frac c2(u_t-c\partial_\nu u)^2 -\frac{c^3}2|\nabla_{\tan}u|^2\le0. \tag{5}

n=1没有切向部分。取相反的球面法向或遗漏移动边界项，都会破坏这个负平方。

零能量处不作除法

由Cauchy–Schwarz不等式，

E'(t)\le\sqrt{2E(t)}\,\|F(t)\|_{L^2(B_t)}.

若E可能为零，不能直接除以 $\sqrt E$ 。令 $Q_\varepsilon=\sqrt{2E+\varepsilon^2}$ ，则

Q_\varepsilon' =\frac{E'}{\sqrt{2E+\varepsilon^2}} \le\|F(t)\|_{L^2(B_t)}.

积分后令ε下降到零，得到式(2)。证明只使用有限锥内的积分，不涉及远处边界。

从导数为零走到函数相同

对两解之差w，底球初位移相同使 $\nabla w(0)=0$ ，初速度相同使 $w_t(0)=0$ ，故Q(0)=0。源差也为零，式(2)给锥内每个时刻的能量为零。导数连续，所以 $w_t=\nabla w=0$ 在锥内部处处成立。

锥内部连通，w因此为常数；底面的w=0才将这个常数固定为零。连续性再给侧面和顶点。任取一个时空点都可用以它为顶点的有限锥，因而得到整空间唯一性。

半线x>0的零Dirichlet问题也可局部处理：把收缩依赖区间与半线相交，移动端仍给非正平方；若截到墙x=0，那里w=0使w_t=0，墙端通量为零。同样得到局部差为零，无须假定整个半线上的总能量有限或限制无穷远增长。

例子与边界

什么时刻以前可以忽略远处数据

取c=2，观察点为x0=3、T=1。其初始依赖区间为[1,5]。两组初值只要在该区间相同，且源在过去锥 $|x-3|\le2(1-t)$ 内相同，顶点答案就相同；区间左侧的巨大扰动仍不能在这一时刻影响x0。

若两组初值差和源差都支撑在固定球B(0,a)内，锥不与这份支撑相交的点没有差异。于是对t≥0，解差的支撑包含于B(0,a+ct)。这里没有说球内处处都不为零，也没有说波前经过之后内部必恢复零值。后一性质要另外核空间维数和传播核。

常数位移不能从梯度能量中找出来

u=0与v=7都满足齐次波方程。它们的时间导数和空间梯度完全相同，任意局部球能量都为零，初位移却不同。这说明初位移相同不能被“初始能量差为零”替换。

有固定零Dirichlet边界时，边界条件可以固定常数；周期边界或零Neumann边界则仍可能保留它。也不能在任意维数中只由Q小推出u逐点小，那需要额外的边界、平均值或Sobolev控制。

不存在一个通用于所有PDE的速度锥

热核在任意正时间、任意远处都严格为正；它没有本页这种固定有限c的支撑界。变系数波算子的传播区域又由主部控制，不能把某一点的系数读成全域常速。数值算法的CFL条件还要比较离散模板的依赖域，连续方程的证明并未自动给任何格式稳定性。

推论与应用

固定区域上的误差预算

设Ω为有界C1区域，u、v为闭时间窗口和闭区域邻域上的C2函数，都满足零Dirichlet边界。u满足 $u_{tt}-c^2\Delta u=F$ ，v是候选近似。定义

w=v-u,\qquad \mathcal R=v_{tt}-c^2\Delta v-F.

边界上w=0对时间求导得w_t=0，所以式(3)在固定区域积分时边界通量为零。令

Q_\Omega(t)^2=\|w_t(t)\|_2^2+c^2\|\nabla w(t)\|_2^2,

同一正则化论证给

Q_\Omega(t)\le Q_\Omega(0)+\int_0^t\|\mathcal R(s)\|_2ds. \tag{6}

再由 $w(t)=w(0)+\int_0^tw_t(s)ds$ 及Minkowski积分不等式，

\|w(t)\|_2\le\|w(0)\|_2+tQ_\Omega(0) +\int_0^t(t-s)\|\mathcal R(s)\|_2ds. \tag{7}

式(6)控制速度和梯度，式(7)另保留位移误差。零Neumann边界也会消掉通量，但不能省初位移项。全空间中，若函数及所需导数在当前有限时间窗具有共同紧支撑，截断边界在支撑外，证明同样成立；若只知某种衰减，应先验证无穷远通量可消去。

一个完整的残差计算

在0<x<1、0≤t≤1上取c=2、F=0，精确问题初位移和初速度都为零，零边界。考虑候选

v(t,x)=\eta t^2x(1-x).

它满足同样的初值和边界，残差却为

\mathcal R=2\eta x(1-x)+8\eta t^2.

由 $\int_0^1x^2(1-x)^2dx=1/30$ ，

\|\mathcal R(t)\|_2\le\frac{2|\eta|}{\sqrt{30}}+8|\eta|t^2, \qquad Q_\Omega(t)\le |\eta|\left(\frac{2t}{\sqrt{30}}+\frac{8t^3}3\right). \tag{8}

本例精确解u=0，还可独立积分候选的导数，得到

Q_\Omega(t)^2=\eta^2\left(\frac{2t^2}{15}+\frac{4t^4}3\right).

这提供了对预算的交叉核验。残差不是零，但只要η小，仍有可明确报告的误差范围；把“方程没有精确满足”与“完全无法认证”混为一谈，会丢掉能量估计的用途。

参考资料

Victor Ivrii，Partial Differential Equations，§9.2，局部能量式(9.2.2)–(9.2.3)、Theorems2–4：特征锥和Cauchy/边界唯一性。本页展开移动球边界平方及有源正则化估计。
John K. Hunter，Notes on Partial Differential Equations，Chapter7，§7.1.1、印刷p.212：局部能量流与零Dirichlet边界。本文的C2论证不替代该章后续弱解存在性理论。
Julie Levandosky，Wave Equation in R，§5.4.3、印刷pp.9–11：一维收缩区间的端点平方与依赖域。

形式陈述 ​

收缩球中的有源估计 ​

依赖域与唯一性 ​

局部守恒先保留通量方向 ​

零能量处不作除法 ​

从导数为零走到函数相同 ​

什么时刻以前可以忽略远处数据 ​

常数位移不能从梯度能量中找出来 ​

不存在一个通用于所有PDE的速度锥 ​

固定区域上的误差预算 ​

一个完整的残差计算 ​

反向链接