Skip to content

模型Model

Hamilton–Jacobi方程

Hamilton–Jacobi equation · Hamilton-Jacobi equation

由Hamilton相流和作用量构造一阶PDE的经典解,证明梯度恢复恒等式,并用谐振子焦散区分相空间唯一性与位置投影退化。

同一群粒子的轨道可以在位置空间重叠,而在位置与动量组成的相空间中仍各走各的。Hamilton–Jacobi方程问的是:能否用一个标量函数 u(t,x) 的空间梯度,给这些轨道统一记录动量?答案在短时间内常常是肯定的;当位置投影折叠时,标量函数的经典梯度描述会先失效。

形式陈述 ​

一个势函数描述位置与动量 ​

本页的 Hamilton–Jacobi方程 是一阶标量PDE

(1)ut+H(x,Dxu)=0,u(0,x)=g(x).

H(x,p) 称Hamilton量,p 与 x 具有相同维数。经典解 u 至少为 C1,使等式逐点有意义。下面的经典构造假设 H∈C2、g∈C2,且所有轨道暂时留在系数定义的开域内。

给每个初始标签 ξ 解Hamilton系统

(2)X˙=Hp(X,P),P˙=−Hx(X,P),X(0,ξ)=ξ,P(0,ξ)=Dg(ξ).

再记录作用量

(3)U(t,ξ)=g(ξ)+∫0t[P(s,ξ)⋅Hp(X(s,ξ),P(s,ξ))−H(X(s,ξ),P(s,ξ))]ds.

若在考察点 DξX(t,ξ) 可逆,就能在该点附近反解标签 ξ=Ξ(t,x),并定义

(4)u(t,x)=U(t,Ξ(t,x)).

此函数满足(1),且 Dxu(t,X(t,ξ))=P(t,ξ)。在 t=0,位置导数为单位矩阵,所以每个初始点附近都有一个足够短的有效时间窗。

直觉

为什么动量必须沿Hamilton方程变化 ​

先假设 u 足够光滑,令 P(t)=Dxu(t,X(t))。对(1)按空间微分,再选 X˙=Hp,得到

P˙=Dxut+Dx2uX˙=−Hx−Dx2uHp+Dx2uHp=−Hx.

沿同一条曲线,

ddtu(t,X(t))=ut+Dxu⋅X˙=−H+P⋅Hp,

正是(3)的积分量。这是从已有解推出特征方程的方向;反向构造还要排除“积分出来的 U 与 P 没有梯度关系”的可能。

作用量的标签导数补上反向证明 ​

由ODE参数可微性,(2)的流对标签为 C1。固定任意标签方向,记沿它的变化为 δX,δP,δU。对(3)的被积函数微分:

ddtδU=δP⋅Hp+P⋅δ(Hp)−Hx⋅δX−Hp⋅δP=P⋅δX˙+P˙⋅δX=ddt(P⋅δX).

初始时 δU=Dg(ξ)⋅δξ=P(0)⋅δX(0),故

(5)DξU=PTDξX.

矩阵里的 DξU 取行向量。它说明标签方向的作用量差,恰由动量作用在位置差上给出。

对映射 (t,ξ)↦(t,X(t,ξ)) 用反函数定理;其可逆性等价于 DξX 可逆。将(5)与(4)的链式求导比较,便得到 Dxu=P。再由 Ut=ut+Dxu⋅Xt 得 ut=−H。存在性证明因此没有要求事先拥有一个待求PDE解。

例子与边界

自由传播与初值曲率 ​

取 H(p)=|p|2/2。此时 P(t,ξ)=Dg(ξ),

X(t,ξ)=ξ+tDg(ξ),DξX=I+tD2g(ξ),U=g(ξ)+t2|Dg(ξ)|2.

若一维 g(ξ)=kξ2/2,在 1+kt≠0 的局部区域有

u(t,x)=kx22(1+kt).

k≥0 时正时间没有这类投影退化;k<0 时 t=−1/k 出现坍缩。局部行列式不为零只保证局部逆,不能单凭它宣称任意大范围内的标签映射全局一一对应。

谐振子在相空间里没有碰撞 ​

取一维

H(x,p)=x2+p22,g(x)=0.

Hamilton方程给出

X(t,ξ)=ξcos⁡t,P(t,ξ)=−ξsin⁡t,U(t,ξ)=−ξ24sin⁡(2t).

对 |t|<π/2,

u(t,x)=−x22tan⁡t.

直接回代:ut=−x2sec2⁡t/2,Dxu=−xtan⁡t,两项在(1)中抵消。

当 t=π/2,X=0 对所有标签都成立,P=−ξ 却彼此不同。相空间中的点 (X,P)=(0,−ξ) 没有合并;退化发生在只看 X 的投影。这里称这种位置投影奇异为焦散。短时间经典解的失败与ODE初值唯一性不矛盾。

左图只观察位置,右图同时保留动量;两图中的终时相同,右图没有把不同标签误画成同一个状态。

推论与应用

光滑梯度能转成守恒律 ​

在一维且 H=H(p) 时,若 u 为 C2,对方程按 x 求导,令 w=ux,便有

wt+∂xH(w)=0.

对 H(w)=w2/2,这是Burgers守恒律。这一步有明确的光滑性前提;在尖点之后,需要分别指定Hamilton–Jacobi方程的黏性解与守恒律的熵解,不能把逐点微分当作一般弱解证明。

一阶黏性解用上下接触函数解释不可微点;对二次自主Hamilton量,Hopf–Lax公式直接取一族代价的最小值,并给出经典特征描述失效后仍可认证的解。其唯一性依赖明确的解类与比较定理,不来自随意挑一条穿过焦散的射线。

参考资料
关系图谱7 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

下位 / 直接特例

暂未标注直接特例。

类型化关系