Skip to content

定理Theorem

椭圆内部正则性

Interior elliptic regularity · 椭圆方程内部H2正则性

对平方可积右端的Laplace弱方程,用能量截断和Fourier逆乘子证明内部二阶正则性,并以凹角奇性区分内部估计与全局边界估计。

形式陈述 ​

设 Ω⊂Rd 为开集,u∈Hloc1(Ω)、f∈Lloc2(Ω) 为实值函数,且 u 是 −Δu=f 的弱解:

∫Ω∇u⋅∇ϕ=∫Ωfϕ(ϕ∈Cc∞(Ω)).

局部 Sobolev 空间的含义是:在每个紧包含于 Ω 的开集上,函数及相应弱导数属于所指定的 L2 空间。则

u∈Hloc2(Ω).

更定量地,若 U,V 是开集且 U― 紧包含于 V、V― 紧包含于 Ω,则存在只依赖 U,V,d 的常数 C,使

‖u‖H2(U)≤C(‖f‖L2(V)+‖u‖L2(V)).

这里的 H2 范数使用所有至多二阶弱导数。定理增加了弱导数的可积阶数;它不直接断言 u 是经典 C2 解。结论属于椭圆型偏微分方程的内部理论,没有对 ∂Ω 或边界数据作假设,所以也没有承诺整个 Ω 上的 H2 正则性。

直觉

方程把二阶导数的某个组合与已知右端联系起来。对 Laplace 算子,频域中的组合是 4π2|ξ|2,能同时控制各方向的二阶导数。困难在于局部弱解没有定义在整个空间,不能直接作全局 Fourier 变换。

光滑截断把这个局部问题变成一个整空间问题,但乘积求导会生成含 u 和 ∇u 的误差。先用方程自己的能量恒等式控制一阶导数,再调用整空间逆乘子,就能让这些误差保持平方可积。截断越靠近边界,导数可能越大;这解释了估计为何依赖内部区域的间隔,也解释了为什么不能把常数原样推到一个凹角上。

例子与边界

一个可直接检查的内部估计对象 ​

取 u(x)=|x|2,则 ∇u=2x、D2u=2I,并且 −Δu=−2d。它在每个有界内部区域上属于 H2,但在整空间不属于 L2。本页定理可以逐个内部区域使用;它没有要求这个局部解具有全局可积性。

若在球 B(0,R) 上把所有 Hessian 矩阵分量的平方相加,则

∫B(0,R)∑i,j|∂iju|2=4d|B(0,R)|.

二阶导数确实受到一个有限的局部估计控制。这里为方便核算使用完整 Hessian 矩阵;普通 H2 范数只按多重指标计数,非对角混合导数的计数略有不同,两者等价。

凹角怎样破坏全局二阶正则性 ​

在二维取开扇形

Ω={(r,θ):0<r<1,0<θ<3π/2},α=23.

函数 h(r,θ)=rαsin⁡(αθ) 在 Ω 内调和,在两条径向边上为零。极坐标 Laplace 公式直接验证

Δh=[α(α−1)+α−α2]rα−2sin⁡(αθ)=0.

梯度大小为 |∇h|=αrα−1,所以靠近顶点的梯度平方积分有限;例如对 0<R<1,

∫03π/2∫0R|∇h|2rdrdθ=π2R4/3.

但逐项求导可得

∑i,j=12|∂ijh|2=2α2(α−1)2r2α−4=881r−8/3.

其面积积分含 ∫0Rr−5/3dr,因而发散。它在每个远离顶点的内部区域光滑,却不在整个扇形的 H2 中。

还可使这个反例满足全边界零条件。取光滑径向函数 χ,在 r≤1/3 时等于 1、在 r≥2/3 时等于 0,令 u=χh。则 u∈H01(Ω),并有

−Δu=f=−[χ″(r)rα+(2α+1)χ′(r)rα−1]sin⁡(αθ).

f 支撑于远离顶点的环带,故属于 L2(Ω);而 u=h 在顶点附近成立,仍不属于 H2(Ω)。

这里的零边界资格也可以从闭包验证:先在 r<δ 附近截去顶点,新增梯度项的平方积分为 O(δ2α),原有梯度在小扇形中的积分也趋于零。留下的函数光滑、远离顶点,并在直边上为零;再沿直边截断,函数的线性消失抵消截断导数,得到紧支撑光滑逼近。于是 u 的确是能量空间中的零边界弱解。原方程先对紧支撑测试函数成立,再由 H01 密度扩展。失败发生在凹角处的全局二阶可积性,没有违反内部定理。

推论与应用

先用能量估计控制一阶导数 ​

取实值 ζ∈Cc∞(V),0≤ζ≤1。ζ2u 属于 H01(V);用局部光滑逼近可以把弱方程的测试范围扩展到它。展开得到

∫ζ2|∇u|2=∫fζ2u−2∫ζu∇u⋅∇ζ.

右端只在 V 内积分。Cauchy–Schwarz 与 2ab≤12a2+2b2 给出

2|∫ζu∇u⋅∇ζ|≤12∫ζ2|∇u|2+2∫u2|∇ζ|2.

另外,由 Hölder 不等式,

|∫fζ2u|≤12‖f‖L2(V)2+12‖u‖L2(V)2.

吸收左侧的一半,得到

∫ζ2|∇u|2≤‖f‖L2(V)2+(1+4‖∇ζ‖∞2)‖u‖L2(V)2.

因此,只要 ζ=1 于某个较小区域,该区域的梯度便受到 f 与 u 的外层 L2 范数控制。这种局部能量估计通常称为 Caccioppoli 估计;上式保留了本页实际需要的版本。

截断后逐项检查右端 ​

选择 U⋐W⋐V,并取 η∈Cc∞(W)、η=1 于 U 的邻域。将 v=ηu 零延拓至 Rd,它属于 H1(Rd)。因为截断支撑离边界有正距离,乘积公式可以在弱导数意义下使用,不产生边界分布。

直接展开,

(1−Δ)v=η(f+u)−2∇η⋅∇u−(Δη)u=:g.

每一项都在 L2(Rd) 中。再选择 ζ∈Cc∞(V),使它在 W― 的邻域等于 1,上一节便给出

‖g‖2≤CU,V(‖f‖L2(V)+‖u‖L2(V)).

估计中的常数包括所选截断函数的至多二阶导数界。我们只要求对固定内部区域存在这样的常数,没有声称它在 U 接近 ∂V 时仍一致。

逆乘子完成二阶提升 ​

Fourier 乘子已经证明,对任意 g∈L2(Rd),

w=F−1(g^(ξ)1+4π2|ξ|2)

是 (1−Δ)w=g 的唯一 H1 弱解,且 w∈H2,满足 ‖w‖H2≤Cd‖g‖2。上节构造的 v 也是同一方程的 H1 弱解,因此 v=w。

唯一性可以在当前接口再核对一次:差 z=v−w 满足齐次弱方程,以 z 本身作测试得到 ‖z‖22+‖∇z‖22=0,故差为零。于是 v∈H2;由于 v=u 在 U 上成立,限制回去便得到

‖u‖H2(U)≤‖v‖H2(Rd)≤C(‖f‖L2(V)+‖u‖L2(V)).

对任意内部 U 使用这一结论,就完成了 Hloc2 的证明。顺序很关键:先用已有的一阶正则性使截断右端合法,再由逆算子取得二阶正则性,没有在乘积展开时预设结论。

连续正则性给后续分析提供什么 ​

弱解的存在与弱解的正则性回答不同问题。变分形式先在能量空间求解线性边值问题;本页进一步证明,对 Laplace 算子和局部平方可积右端,所得弱解在内部具有二阶弱导数。直接法中的三维方程 −Δu+u3=f 也可接入:H01↪L6 保证 u3∈L2,故把 f−u3 当作右端,就得到 u∈Hloc2。

有限元方法的误差估计若需要整个区域的 H2 范数,必须另核对边界几何和边界数据。凹角例子表明,仅有内部估计与 L2 右端不足以支持这种全局代入。变系数的一致椭圆方程还需控制系数的正则性;本页证明的是常系数 Laplace 情形,不能直接把同一 Fourier 符号用于任意变系数。

参考资料
  • John K. Hunter, Notes on Partial Differential Equations, Chapter 4, UC Davis,2014,§4.11,定理 4.27,印刷 pp. 112–113:一致椭圆、C1 系数下的内部 H2 正则性,原文用差商证明。本文只取常系数 Laplace 情形,改用能量截断与本站已证明的 Fourier 逆乘子给出完整推导。
  • John K. Hunter, Notes on Partial Differential Equations,修订于 2014-06-18,附录 5.B.4–5.C,pp. 171–173:L2 Fourier 变换与 Sobolev 频域刻画。二次多项式和 3π/2 凹角算例均在正文独立计算。
关系图谱12 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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