形式陈述
从强形式到弱形式
设整数 d ≥ 1 ,以有界 Lipschitz 区域 Ω ⊂ R d 上的偏微分方程 理路 偏微分方程 Partial differential equation · PDE 含未知多元函数及其偏导数的方程。 Poisson 问题为主线:
− Δ u = f in Ω , u = 0 on ∂ Ω . 强形式要求微分方程与边界条件具有逐点或经典迹意义。先取在 Ω ― 的邻域上定义的 u ∈ C 2 、v ∈ C 1 。对 C 1 向量场 v ∇ u 使用散度定理的 Lipschitz 区域版本 理路 散度定理 Divergence theorem · Gauss theorem · Gauss–Ostrogradsky theorem 区域内部散度的体积分等于其边界上向外法向通量的积分。 ,并展开 ∇ ⋅ ( v ∇ u ) = ∇ v ⋅ ∇ u + v Δ u ,可得
∫ Ω ( − Δ u ) v d x = ∫ Ω ∇ u ⋅ ∇ v d x − ∫ ∂ Ω ∂ u ∂ n v d s . 在Sobolev 空间 理路 Sobolev 空间 Sobolev space 同时控制函数与弱导数的可积性,在一维中建立端点迹与能量空间,证明二分之一 Hölder 估计和点值泛函的精确范数,并说明高维边界。 中取 V = H 0 1 ( Ω ) 后,v 在边界上的迹为零,边界项消失。弱问题写成:求 u ∈ V ,使
a ( u , v ) = ℓ ( v ) ( ∀ v ∈ V ) , 其中
a ( w , v ) = ∫ Ω ∇ w ⋅ ∇ v d x , ℓ ( v ) = ∫ Ω f v d x . 这个改写把二阶导数转移到一次弱导数 理路 弱导数 Weak derivative 通过分部积分恒等式相对于测试函数定义的广义导数。 上,因此允许 u 只属于 H 0 1 ( Ω ) 。“弱”描述正则性和检验方式,不表示数值精度较低。足够光滑的弱解仍满足强方程;具有相应边界条件的强解也满足弱形式。
边界条件进入哪里
对
− ∇ ⋅ ( κ ∇ u ) + c u = f 若边界分成两个可测部分 Γ D 与 Γ N ,并给定
u = g D on Γ D , κ ∇ u ⋅ n = g N on Γ N , 则试探集合和检验空间分别为
U g D = { w ∈ H 1 ( Ω ) : tr w = g D on Γ D } , V = { v ∈ H 1 ( Ω ) : tr v = 0 on Γ D } . 双线性型与右端为
a ( w , v ) = ∫ Ω κ ∇ w ⋅ ∇ v d x + ∫ Ω c w v d x , ℓ ( v ) = ∫ Ω f v d x + ∫ Γ N g N v d s . Dirichlet 条件限制可容许试探函数,通常称 essential boundary condition;Neumann 数据经分部积分自然进入右端,因此称 natural boundary condition。非齐次 Dirichlet 情形中 U g D 是仿射集合而 V 是线性空间,trial 与 test 不能不加条件地视为同一对象。
若经典超调和解满足 Δ u ≤ 0 (例如 − Δ u = f ≥ 0 ),在内部严格为正、指定边界点取值为零,且该点具有 C 1 边界和内切球,解在该点具有连续的一阶导数,则对 − u 应用Hopf 引理 理路 Hopf 边界点引理 Hopf boundary point lemma · Boundary point lemma 在内切球接触的边界最大点,严格低于该最大值的次调和函数具有严格正的外法向导数。 可得严格负的外法向导数。这项点态通量结论同时需要方程符号、正则性与几何,不能仅由弱形式自动推出。
要放入标准适定性框架,还需说明数据和系数条件。例如在有界连通 Lipschitz 区域上,可要求 Γ D 具有正表面测度,g D 是某个 H 1 函数的可容许迹,实值系数 κ , c 可测,并且几乎处处满足
0 < κ 0 ≤ κ ( x ) ≤ κ 1 < ∞ , 0 ≤ c ( x ) ≤ c 1 < ∞ . 这些条件使齐次空间满足相应 Poincaré 不等式;若区域不连通,则每个没有反应项控制的连通分量都需要自己的 Dirichlet 约束。右端还需是 V 上的连续泛函,例如 f ∈ L 2 ( Ω ) 、g N ∈ L 2 ( Γ N ) 配合迹不等式即可。
先选 lifting G ∈ H 1 ( Ω ) 满足 tr G = g D ,再写 u = G + w ,就把未知量 w 放回齐次空间 V ,求解 a ( w , v ) = ℓ ( v ) − a ( G , v ) 。若 Γ D 为空、系数变号或数据没有所需迹,这条论证不能原样使用。
抽象适定性
设 V 是实Hilbert 空间 理路 Hilbert 空间 Hilbert space 关于内积诱导范数完备的实或复内积空间。 ,双线性型 理路 双线性映射与形式 Bilinear map · Bilinear form 对两个向量变量分别线性的映射;值域为标量域时称双线性形式。 a : V × V → R 满足连续性与强制性:
| a ( w , v ) | ≤ M ‖ w ‖ V ‖ v ‖ V , a ( v , v ) ≥ α ‖ v ‖ V 2 , α > 0. 连续性还意味着固定 w ∈ V 后,映射 v ↦ a ( w , v ) 是 V 上的有界线性算子 理路 有界线性算子 Bounded linear operator 把有界集映为有界集,等价地连续的线性映射。 。若线性泛函 ℓ 连续,Lax–Milgram 定理保证弱问题存在唯一解,并有
‖ u ‖ V ≤ 1 α ‖ ℓ ‖ V ′ . 齐次 Dirichlet Poisson 问题中,Cauchy–Schwarz 不等式 理路 Cauchy–Schwarz 不等式 Cauchy–Schwarz inequality · 柯西–施瓦茨不等式 内积的绝对值不超过两向量范数之积,且等号精确刻画线性相关。 给出连续性,Poincaré 不等式使 ‖ ∇ v ‖ L 2 控制 H 1 范数,从而得到强制性。纯 Neumann Poisson 问题保留常数核;只有右端满足兼容条件后,才能在零均值子空间或模去常数的商空间中恢复唯一性。
Lax–Milgram 为什么成立
将“对每个检验函数成立”的无限多个等式转换为一个 Hilbert 空间中的算子方程。由Riesz 表示定理 理路 Hilbert 空间 Riesz 表示定理 Riesz representation theorem for Hilbert spaces Hilbert 空间上的每个连续线性泛函都唯一由与某向量取内积表示。 ,存在唯一的 T w ∈ V 与 b ∈ V ,使
( T w , v ) V = a ( w , v ) , ( b , v ) V = ℓ ( v ) . 表示的唯一性使 T 为线性算子,连续性给出 ‖ T w ‖ V ≤ M ‖ w ‖ V 。另一方面,
α ‖ w ‖ V 2 ≤ ( T w , w ) V ≤ ‖ T w ‖ V ‖ w ‖ V , ‖ T w ‖ V ≥ α ‖ w ‖ V . 因此 T 单射,而且值域闭:若 T w n 收敛,则上述下界使 w n 为 Cauchy 序列;完备性给出 w n → w ,连续性给出 T w n → T w 。
还需排除值域遗漏某个方向。若 y 与 ran T 正交,特别取值域中的 T y ,便有
0 = ( T y , y ) V = a ( y , y ) ≥ α ‖ y ‖ V 2 , 所以 y = 0 。对已经证明闭的线性子空间 ran T ,Hilbert 投影定理 理路 Hilbert 空间投影定理 Hilbert projection theorem · Projection theorem Hilbert 空间中每个闭线性子空间都给出唯一的正交分解与最近点投影。 给出 V = ran T ⊕ ( ran T ) ⊥ ;正交补为零便得到 ran T = V 。于是 T u = b 有唯一解,这正是弱问题。最后取 v = u ,
α ‖ u ‖ V 2 ≤ ℓ ( u ) ≤ ‖ ℓ ‖ V ′ ‖ u ‖ V , 得到所述解估计。整个证明没有用到 a 的对称性;对称性是在下一节解释能量最小化时才需要的条件。
例子与边界
在区间上逐项验证适定性
取 V = H 0 1 ( 0 , 1 ) 、‖ v ‖ V = ‖ v ′ ‖ 2 ,并令 a ( w , v ) = ∫ 0 1 w ′ v ′ 。Sobolev 空间中的 Poincaré 证明已保证这确实是完备的范数。此时连续常数和强制常数均可取 M = α = 1 。对于 f ∈ L 2 ( 0 , 1 ) ,
| ℓ ( v ) | = | ∫ 0 1 f v | ≤ ‖ f ‖ 2 ‖ v ‖ 2 ≤ 1 2 ‖ f ‖ 2 ‖ v ‖ V . 故唯一弱解满足 ‖ u ′ ‖ 2 ≤ ‖ f ‖ 2 / 2 。当 f = 1 时,直接积分 − u ″ = 1 并施加两个端点条件,得到
u ( x ) = x ( 1 − x ) 2 , ‖ u ′ ‖ 2 = 1 12 . 这个显式解既核对弱方程,也说明适定性上界无须等于精确范数。
纯 Neumann 条件的符号与兼容性
考虑没有反应项的 − ∇ ⋅ ( κ ∇ u ) = f ,约定给定量是外向 导数通量 g = κ ∇ u ⋅ n 。此时检验空间包含常数 1 ,弱方程取 v = 1 得到
0 = ∫ Ω f d x + ∫ ∂ Ω g d s . 兼容条件中的加号来自弱形式右端的边界项;如果将物理通量定义为 − κ ∇ u ,所记边界数据的符号也随之改变。
在一维、κ = 1 时,若写成 u ′ ( 0 ) = q 0 、u ′ ( 1 ) = q 1 ,左右外法向分别为 − 1 和 1 ,所以外向数据是 g ( 0 ) = − q 0 、g ( 1 ) = q 1 ,条件成为
∫ 0 1 f ( x ) d x = q 0 − q 1 . 零源项、零通量允许任意常数解;f = 1 、两端零通量则根本没有解。兼容性解决的是存在性,固定均值解决的是常数方向造成的不唯一性,二者是不同的问题。在有界连通 Lipschitz 区域、均匀正扩散系数和适当连续右端下,零均值 Poincaré 不等式使兼容问题在零均值子空间上恢复强制性。
弱解与正则性
Lax–Milgram 给出的是所选空间中的解,并不自动把它升级为 H 2 解。区间上的常系数 Poisson 方程可以从 u ″ = − f ∈ L 2 直接得到 H 2 正则性;高维中这一步还依赖边界几何和系数正则性。凹角或系数界面可能让弱解存在唯一,而二阶导数仍不平方可积。椭圆内部正则性 理路 椭圆内部正则性 Interior elliptic regularity · 椭圆方程内部H2正则性 对平方可积右端的Laplace弱方程,用能量截断和Fourier逆乘子证明内部二阶正则性,并以凹角奇性区分内部估计与全局边界估计。 给出常系数 Laplace 情形的完整局部证明:能量截断控制一阶项,Fourier 逆乘子提升到内部 H 2 ;这项内部估计不包含全局边界正则性。
弱形式也不是唯一的。不同分部积分次数、边界数据处理或试探—检验空间会产生不同 formulation;声称两者等价时必须核对正则性与边界项。