Skip to content

定理Theorem

Eikonal方程的距离解

Eikonal distance solution · Zero-boundary eikonal equation · 到边界距离的黏性解

在任意非空有界开域上,以到边界距离构造零Dirichlet Eikonal黏性解,用最短线段证明接触条件,再用严格缩放和双变量认证连续闭包解的唯一性。

如果从边界出发,以单位速度向区域内部传播,到一个点最早需要多长时间?在没有方向或位置速度差异时,答案就是到最近边界的距离。多个最近边界点可能同时到达,形成不可微的交汇处;这些点仍必须满足Eikonal方程的黏性接触条件。

形式陈述 ​

区域可以不凸,边界不要求光滑 ​

设 Ω⊂Rn 为非空有界开集。定义

d(x)=dist(x,∂Ω)=mina∈∂Ω|x−a|(x∈Ω―).

则 d 连续,边界上为零,内部严格为正,并且是

(1)|Du|=1 于 Ω,u=0 于 ∂Ω

的唯一连续闭包黏性解。这里解类是 u∈C(Ω―);边界条件指每个边界点的实际值为零,内部使用一阶接触测试。

无需假设 Ω 连通、凸或边界可微。唯一性针对整个有界开域及其全边界;把某些边界分支遗漏,或仅对开域内部取连续函数而不规定边界极限,已是不同问题。

直觉

最近边界给出一条有效的最短线段 ​

∂Ω 非空且紧,由极值定理取得最近点 a。若 x∈Ω,令 r=d(x)>0,则开球 B(x,r) 包含在 Ω 中。否则从 x 向一个球内的域外点走,必先碰到距离小于 r 的边界点,与定义矛盾。

令 e=(x−a)/r。当 0<h<r 时,x−he∈Ω,且

(2)d(x−he)=r−h.

上界来自同一个边界点 a,下界来自距离函数的Lipschitz性质。即使整个区域不凸,这一条走向最近边界的短线段仍全部有效。

距离函数的斜率不会超过一 ​

三角不等式给 |x−a|≤|x−y|+|y−a|,对 a 取最小值并交换 x,y,得到

|d(x)−d(y)|≤|x−y|.

因此 d 为 1-Lipschitz函数。这限制了所有接触斜率;(2)又提供一个确实按单位速度下降的方向。两部分分别负责次解和超解条件。

例子与边界

区间中央的尖点不能略过 ​

在 Ω=(−1,1),d(x)=1−|x|。零点以外斜率为 ±1;在零点,上接触斜率允许全部 p∈[−1,1],满足 |p|−1≤0,下接触不存在。于是尖点通过黏性检查。

v(x)=|x|−1 的边界值也为零,并在零点以外满足 |v′|=1,但常数 −1 从下接触它,违反 |ϕ′|−1≥0。到达时间不能用这个负帐篷代替。

圆环同时接收两侧边界的传播 ​

若 Ω={x:1<|x|<3},则

d(x)=min{|x|−1, 3−|x|}.

在 |x|=2 的圆周上,两侧边界同时到达,距离不可微。若只写成 3−|x|,它在外圆边界为零,却在内圆边界等于二,未满足本题边界条件。边界的所有分支都会参与最早到达值。

非零边界不能任意指定 ​

一维 |u′|=1 的连续黏性解具有相应的斜率限制。若把区间两端相距二的边界值差指定为十,就不再是本页零边界结论的直接推广。一般边界数据与传播代价之间需要兼容条件;本页不宣称任意连续边界值都可用同一距离函数处理。

推论与应用

上接触来自Lipschitz界,下接触来自最短方向 ​

若 ϕ 在内部点 x 从上接触 d,Lipschitz估计和任意方向小增量给 |Dϕ(x)|≤1,所以次解条件成立。

若从下接触,取前述最近边界方向 e。由(2),

ϕ(x−he)−ϕ(x)≤d(x−he)−d(x)=−h.

除以 h>0 并令其趋零,得到 Dϕ(x)⋅e≥1。因为 |e|=1,所以 |Dϕ(x)|≥1,超解条件成立。这也解释了某些点为什么没有下接触:若有多个相反的最短方向,它们可能要求同一个梯度同时朝相反方向投影至少一。

严格缩放把比较变成不可能的双斜率 ​

设 u,v∈C(Ω―),分别为(1)的次解和超解,边界值都为零。固定 0<μ<1。若 ϕ 从上接触 μu,则 ϕ/μ 从上接触 u,故

(3)|Dϕ|≤μ.

它成为一个严格次解。我们证明 μu≤v。

反设某个内部点 x0 满足 μu(x0)−v(x0)=δ>0。在紧集 Ω―2 上最大化

(4)Φα(x,y)=μu(x)−v(y)−α2|x−y|2.

最大值至少为 δ,由连续函数有界性可得 |xα−yα|→0。取任意收敛子列,两点趋于同一个 z;若 z∈∂Ω,连续性和零边界使差 μu(xα)−v(yα)→0,与最大值至少 δ 矛盾。因此对充分大的 α,两个最大点都在域内。

固定另一变量,(4)分别给 μu 的上接触与 v 的下接触,二者空间梯度恰为同一个

p=α(xα−yα).

由(3)得 |p|≤μ,由超解条件得 |p|≥1,矛盾。所以 μu≤v;令 μ↑1 得 u≤v。若两者都是解,交换角色便得相等。

证明只使用闭包连续、域有界和零边界,没有预先要求一般候选解Lipschitz;因而唯一性解类确实如形式陈述所写。存在性已由距离函数的两个接触检查完成。

单位速度怎样换成固定速度 ​

若传播速度为常数 c>0,到达时间为 τ=d/c,满足 c|Dτ|=1 和零边界。将任意候选乘以 c 即回到(1),同样给唯一性。若速度随位置变化,直线距离通常不再是时间,需要沿路径累计代价;不能仅把 d(x) 除以终点的速度就结束计算。

特征与前沿证书练习把最短方向、尖点测试和Hamilton–Jacobi包络放在一起,要求分别交付局部经典构造与全域黏性唯一性的证据。

参考资料
关系图谱5 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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