形式陈述
设 Ω ⊂ R d 为开集,u ∈ H loc 1 ( Ω ) 、f ∈ L loc 2 ( Ω ) 为实值函数,且 u 是 − Δ u = f 的弱解:
∫ Ω ∇ u ⋅ ∇ ϕ = ∫ Ω f ϕ ( ϕ ∈ C c ∞ ( Ω ) ) . 局部 Sobolev 空间 公理库 Sobolev 空间 Sobolev space 同时控制函数与弱导数的可积性,在一维中建立端点迹与能量空间,证明二分之一 Hölder 估计和点值泛函的精确范数,并说明高维边界。 的含义是:在每个紧包含于 Ω 的开集上,函数及相应弱导数属于所指定的 L 2 空间。则
u ∈ H loc 2 ( Ω ) . 更定量地,若 U , V 是开集且 U ― 紧包含于 V 、V ― 紧包含于 Ω ,则存在只依赖 U , V , d 的常数 C ,使
‖ u ‖ H 2 ( U ) ≤ C ( ‖ f ‖ L 2 ( V ) + ‖ u ‖ L 2 ( V ) ) . 这里的 H 2 范数使用所有至多二阶弱导数。定理增加了弱导数的可积阶数;它不直接断言 u 是经典 C 2 解。结论属于椭圆型偏微分方程 公理库 偏微分方程 Partial differential equation · PDE 含未知多元函数及其偏导数的方程。 的内部理论,没有对 ∂ Ω 或边界数据作假设,所以也没有承诺整个 Ω 上的 H 2 正则性。
直觉
方程把二阶导数的某个组合与已知右端联系起来。对 Laplace 算子,频域中的组合是 4 π 2 | ξ | 2 ,能同时控制各方向的二阶导数。困难在于局部弱解没有定义在整个空间,不能直接作全局 Fourier 变换。
光滑截断把这个局部问题变成一个整空间问题,但乘积求导会生成含 u 和 ∇ u 的误差。先用方程自己的能量恒等式控制一阶导数,再调用整空间逆乘子,就能让这些误差保持平方可积。截断越靠近边界,导数可能越大;这解释了估计为何依赖内部区域的间隔,也解释了为什么不能把常数原样推到一个凹角上。
例子与边界
一个可直接检查的内部估计对象
取 u ( x ) = | x | 2 ,则 ∇ u = 2 x 、D 2 u = 2 I ,并且 − Δ u = − 2 d 。它在每个有界内部区域上属于 H 2 ,但在整空间不属于 L 2 。本页定理可以逐个内部区域使用;它没有要求这个局部解具有全局可积性。
若在球 B ( 0 , R ) 上把所有 Hessian 矩阵分量的平方相加,则
∫ B ( 0 , R ) ∑ i , j | ∂ i j u | 2 = 4 d | B ( 0 , R ) | . 二阶导数确实受到一个有限的局部估计控制。这里为方便核算使用完整 Hessian 矩阵;普通 H 2 范数只按多重指标计数,非对角混合导数的计数略有不同,两者等价。
凹角怎样破坏全局二阶正则性
在二维取开扇形
Ω = { ( r , θ ) : 0 < r < 1 , 0 < θ < 3 π / 2 } , α = 2 3 . 函数 h ( r , θ ) = r α sin ( α θ ) 在 Ω 内调和,在两条径向边上为零。极坐标 Laplace 公式直接验证
Δ h = [ α ( α − 1 ) + α − α 2 ] r α − 2 sin ( α θ ) = 0. 梯度大小为 | ∇ h | = α r α − 1 ,所以靠近顶点的梯度平方积分有限;例如对 0 < R < 1 ,
∫ 0 3 π / 2 ∫ 0 R | ∇ h | 2 r d r d θ = π 2 R 4 / 3 . 但逐项求导可得
∑ i , j = 1 2 | ∂ i j h | 2 = 2 α 2 ( α − 1 ) 2 r 2 α − 4 = 8 81 r − 8 / 3 . 其面积积分含 ∫ 0 R r − 5 / 3 d r ,因而发散。它在每个远离顶点的内部区域光滑,却不在整个扇形的 H 2 中。
还可使这个反例满足全边界零条件。取光滑径向函数 χ ,在 r ≤ 1 / 3 时等于 1 、在 r ≥ 2 / 3 时等于 0 ,令 u = χ h 。则 u ∈ H 0 1 ( Ω ) ,并有
− Δ u = f = − [ χ ″ ( r ) r α + ( 2 α + 1 ) χ ′ ( r ) r α − 1 ] sin ( α θ ) . f 支撑于远离顶点的环带,故属于 L 2 ( Ω ) ;而 u = h 在顶点附近成立,仍不属于 H 2 ( Ω ) 。
这里的零边界资格也可以从闭包验证:先在 r < δ 附近截去顶点,新增梯度项的平方积分为 O ( δ 2 α ) ,原有梯度在小扇形中的积分也趋于零。留下的函数光滑、远离顶点,并在直边上为零;再沿直边截断,函数的线性消失抵消截断导数,得到紧支撑光滑逼近。于是 u 的确是能量空间中的零边界弱解。原方程先对紧支撑测试函数成立,再由 H 0 1 密度扩展。失败发生在凹角处的全局二阶可积性,没有违反内部定理。
推论与应用
先用能量估计控制一阶导数
取实值 ζ ∈ C c ∞ ( V ) ,0 ≤ ζ ≤ 1 。ζ 2 u 属于 H 0 1 ( V ) ;用局部光滑逼近可以把弱方程的测试范围扩展到它。展开得到
∫ ζ 2 | ∇ u | 2 = ∫ f ζ 2 u − 2 ∫ ζ u ∇ u ⋅ ∇ ζ . 右端只在 V 内积分。Cauchy–Schwarz 与 2 a b ≤ 1 2 a 2 + 2 b 2 给出
2 | ∫ ζ u ∇ u ⋅ ∇ ζ | ≤ 1 2 ∫ ζ 2 | ∇ u | 2 + 2 ∫ u 2 | ∇ ζ | 2 . 另外,由 Hölder 不等式 公理库 Hölder 不等式 Hölder's inequality 共轭指数下函数乘积的 L¹ 范数由各自的 Lᵖ 范数乘积控制。 ,
| ∫ f ζ 2 u | ≤ 1 2 ‖ f ‖ L 2 ( V ) 2 + 1 2 ‖ u ‖ L 2 ( V ) 2 . 吸收左侧的一半,得到
∫ ζ 2 | ∇ u | 2 ≤ ‖ f ‖ L 2 ( V ) 2 + ( 1 + 4 ‖ ∇ ζ ‖ ∞ 2 ) ‖ u ‖ L 2 ( V ) 2 . 因此,只要 ζ = 1 于某个较小区域,该区域的梯度便受到 f 与 u 的外层 L 2 范数控制。这种局部能量估计通常称为 Caccioppoli 估计;上式保留了本页实际需要的版本。
截断后逐项检查右端
选择 U ⋐ W ⋐ V ,并取 η ∈ C c ∞ ( W ) 、η = 1 于 U 的邻域。将 v = η u 零延拓至 R d ,它属于 H 1 ( R d ) 。因为截断支撑离边界有正距离,乘积公式可以在弱导数意义下使用,不产生边界分布。
直接展开,
( 1 − Δ ) v = η ( f + u ) − 2 ∇ η ⋅ ∇ u − ( Δ η ) u =: g . 每一项都在 L 2 ( R d ) 中。再选择 ζ ∈ C c ∞ ( V ) ,使它在 W ― 的邻域等于 1 ,上一节便给出
‖ g ‖ 2 ≤ C U , V ( ‖ f ‖ L 2 ( V ) + ‖ u ‖ L 2 ( V ) ) . 估计中的常数包括所选截断函数的至多二阶导数界。我们只要求对固定内部区域存在这样的常数,没有声称它在 U 接近 ∂ V 时仍一致。
逆乘子完成二阶提升
Fourier 乘子 公理库 Fourier 乘子与平方可积估计 Fourier multiplier · L2 Fourier multiplier 把频域逐点乘法搬回函数空间,证明L2乘子范数恰为符号的本质上确界,并用一个有界逆乘子构造二阶弱解。 已经证明,对任意 g ∈ L 2 ( R d ) ,
w = F − 1 ( g ^ ( ξ ) 1 + 4 π 2 | ξ | 2 ) 是 ( 1 − Δ ) w = g 的唯一 H 1 弱解,且 w ∈ H 2 ,满足 ‖ w ‖ H 2 ≤ C d ‖ g ‖ 2 。上节构造的 v 也是同一方程的 H 1 弱解,因此 v = w 。
唯一性可以在当前接口再核对一次:差 z = v − w 满足齐次弱方程,以 z 本身作测试得到
‖ z ‖ 2 2 + ‖ ∇ z ‖ 2 2 = 0 ,故差为零。于是 v ∈ H 2 ;由于 v = u 在 U 上成立,限制回去便得到
‖ u ‖ H 2 ( U ) ≤ ‖ v ‖ H 2 ( R d ) ≤ C ( ‖ f ‖ L 2 ( V ) + ‖ u ‖ L 2 ( V ) ) . 对任意内部 U 使用这一结论,就完成了 H loc 2 的证明。顺序很关键:先用已有的一阶正则性使截断右端合法,再由逆算子取得二阶正则性,没有在乘积展开时预设结论。
连续正则性给后续分析提供什么
弱解的存在与弱解的正则性回答不同问题。变分形式 公理库 变分形式与弱问题 Variational formulation · Weak formulation 通过检验函数与分部积分把强微分方程改写为低正则性空间中的连续弱问题,并明确边界条件与适定性。 先在能量空间求解线性边值问题;本页进一步证明,对 Laplace 算子和局部平方可积右端,所得弱解在内部具有二阶弱导数。直接法 公理库 变分法的直接法 Direct method in the calculus of variations 从有界极小化序列、弱极限和下半连续性证明极小点存在,以一维点值与三维紧嵌入分别完成四次势能的唯一弱解及精确算例。 中的三维方程 − Δ u + u 3 = f 也可接入:H 0 1 ↪ L 6 保证 u 3 ∈ L 2 ,故把 f − u 3 当作右端,就得到 u ∈ H loc 2 。
有限元方法 公理库 有限元方法 Finite element method · FEM · Finite element assembly 从弱形式与帽函数出发,完整构造局部单元、组装并求解小系统,再用能量与逼近解释误差阶和边界处理。 的误差估计若需要整个区域的 H 2 范数,必须另核对边界几何和边界数据。凹角例子表明,仅有内部估计与 L 2 右端不足以支持这种全局代入。变系数的一致椭圆方程还需控制系数的正则性;本页证明的是常系数 Laplace 情形,不能直接把同一 Fourier 符号用于任意变系数。
参考资料