Skip to content

原则Principle

变分形式与弱问题

Variational formulation · Weak formulation

通过检验函数与分部积分把强微分方程改写为低正则性空间中的连续弱问题,并明确边界条件与适定性。

形式陈述 ​

从强形式到弱形式 ​

设整数 d≥1,以有界 Lipschitz 区域 Ω⊂Rd 上的偏微分方程 Poisson 问题为主线:

−Δu=fin Ω,u=0on ∂Ω.

强形式要求微分方程与边界条件具有逐点或经典迹意义。先取在 Ω― 的邻域上定义的 u∈C2、v∈C1。对 C1 向量场 v∇u 使用散度定理的 Lipschitz 区域版本,并展开 ∇⋅(v∇u)=∇v⋅∇u+vΔu,可得

∫Ω(−Δu)vdx=∫Ω∇u⋅∇vdx−∫∂Ω∂u∂nvds.

在Sobolev 空间中取 V=H01(Ω) 后,v 在边界上的迹为零,边界项消失。弱问题写成:求 u∈V,使

a(u,v)=ℓ(v)(∀v∈V),

其中

a(w,v)=∫Ω∇w⋅∇vdx,ℓ(v)=∫Ωfvdx.

这个改写把二阶导数转移到一次弱导数上,因此允许 u 只属于 H01(Ω)。“弱”描述正则性和检验方式,不表示数值精度较低。足够光滑的弱解仍满足强方程;具有相应边界条件的强解也满足弱形式。

边界条件进入哪里 ​

对

−∇⋅(κ∇u)+cu=f

若边界分成两个可测部分 ΓD 与 ΓN,并给定

u=gD on ΓD,κ∇u⋅n=gN on ΓN,

则试探集合和检验空间分别为

UgD={w∈H1(Ω):trw=gD on ΓD},V={v∈H1(Ω):trv=0 on ΓD}.

双线性型与右端为

a(w,v)=∫Ωκ∇w⋅∇vdx+∫Ωcwvdx,ℓ(v)=∫Ωfvdx+∫ΓNgNvds.

Dirichlet 条件限制可容许试探函数,通常称 essential boundary condition;Neumann 数据经分部积分自然进入右端,因此称 natural boundary condition。非齐次 Dirichlet 情形中 UgD 是仿射集合而 V 是线性空间,trial 与 test 不能不加条件地视为同一对象。

若经典超调和解满足 Δu≤0(例如 −Δu=f≥0),在内部严格为正、指定边界点取值为零,且该点具有 C1 边界和内切球,解在该点具有连续的一阶导数,则对 −u 应用Hopf 引理可得严格负的外法向导数。这项点态通量结论同时需要方程符号、正则性与几何,不能仅由弱形式自动推出。

要放入标准适定性框架,还需说明数据和系数条件。例如在有界连通 Lipschitz 区域上,可要求 ΓD 具有正表面测度,gD 是某个 H1 函数的可容许迹,实值系数 κ,c 可测,并且几乎处处满足

0<κ0≤κ(x)≤κ1<∞,0≤c(x)≤c1<∞.

这些条件使齐次空间满足相应 Poincaré 不等式;若区域不连通,则每个没有反应项控制的连通分量都需要自己的 Dirichlet 约束。右端还需是 V 上的连续泛函,例如 f∈L2(Ω)、gN∈L2(ΓN) 配合迹不等式即可。

先选 lifting G∈H1(Ω) 满足 trG=gD,再写 u=G+w,就把未知量 w 放回齐次空间 V,求解 a(w,v)=ℓ(v)−a(G,v)。若 ΓD 为空、系数变号或数据没有所需迹,这条论证不能原样使用。

抽象适定性 ​

设 V 是实Hilbert 空间,双线性型 a:V×V→R 满足连续性与强制性:

|a(w,v)|≤M‖w‖V‖v‖V,a(v,v)≥α‖v‖V2,α>0.

连续性还意味着固定 w∈V 后,映射 v↦a(w,v) 是 V 上的有界线性算子。若线性泛函 ℓ 连续,Lax–Milgram 定理保证弱问题存在唯一解,并有

‖u‖V≤1α‖ℓ‖V′.

齐次 Dirichlet Poisson 问题中,Cauchy–Schwarz 不等式给出连续性,Poincaré 不等式使 ‖∇v‖L2 控制 H1 范数,从而得到强制性。纯 Neumann Poisson 问题保留常数核;只有右端满足兼容条件后,才能在零均值子空间或模去常数的商空间中恢复唯一性。

Lax–Milgram 为什么成立 ​

将“对每个检验函数成立”的无限多个等式转换为一个 Hilbert 空间中的算子方程。由Riesz 表示定理,存在唯一的 Tw∈V 与 b∈V,使

(Tw,v)V=a(w,v),(b,v)V=ℓ(v).

表示的唯一性使 T 为线性算子,连续性给出 ‖Tw‖V≤M‖w‖V。另一方面,

α‖w‖V2≤(Tw,w)V≤‖Tw‖V‖w‖V,‖Tw‖V≥α‖w‖V.

因此 T 单射,而且值域闭:若 Twn 收敛,则上述下界使 wn 为 Cauchy 序列;完备性给出 wn→w,连续性给出 Twn→Tw。

还需排除值域遗漏某个方向。若 y 与 ranT 正交,特别取值域中的 Ty,便有

0=(Ty,y)V=a(y,y)≥α‖y‖V2,

所以 y=0。对已经证明闭的线性子空间 ranT,Hilbert 投影定理给出 V=ranT⊕(ranT)⊥;正交补为零便得到 ranT=V。于是 Tu=b 有唯一解,这正是弱问题。最后取 v=u,

α‖u‖V2≤ℓ(u)≤‖ℓ‖V′‖u‖V,

得到所述解估计。整个证明没有用到 a 的对称性;对称性是在下一节解释能量最小化时才需要的条件。

直觉

强形式在每个位置读取微分方程,弱形式则询问方程对所有允许检验方向是否平衡。分部积分把一个导数转移给检验函数,也把边界通量显式暴露出来。这个过程不是把方程含糊化,而是把内部平衡、边界条件和允许正则性放进同一个可检验陈述。

对上述线性双线性型框架,若 a 对称且强制,则可用二次能量表述问题。此时

J(v)=12a(v,v)−ℓ(v)

的唯一极小点是弱解。能量最小化也可处理非线性问题:直接法从含 v4/4 的凸能量得到 −u″+u3=f 的唯一弱解,其存在性依靠弱紧性与下半连续性。对流占优方程、鞍点系统或一般非对称问题未必具有上述二次能量表述;此时真正控制适定性的可能是 inf–sup 条件。

例子与边界

在区间上逐项验证适定性 ​

取 V=H01(0,1)、‖v‖V=‖v′‖2,并令 a(w,v)=∫01w′v′。Sobolev 空间中的 Poincaré 证明已保证这确实是完备的范数。此时连续常数和强制常数均可取 M=α=1。对于 f∈L2(0,1),

|ℓ(v)|=|∫01fv|≤‖f‖2‖v‖2≤12‖f‖2‖v‖V.

故唯一弱解满足 ‖u′‖2≤‖f‖2/2。当 f=1 时,直接积分 −u″=1 并施加两个端点条件,得到

u(x)=x(1−x)2,‖u′‖2=112.

这个显式解既核对弱方程,也说明适定性上界无须等于精确范数。

纯 Neumann 条件的符号与兼容性 ​

考虑没有反应项的 −∇⋅(κ∇u)=f,约定给定量是外向导数通量 g=κ∇u⋅n。此时检验空间包含常数 1,弱方程取 v=1 得到

0=∫Ωfdx+∫∂Ωgds.

兼容条件中的加号来自弱形式右端的边界项;如果将物理通量定义为 −κ∇u,所记边界数据的符号也随之改变。

在一维、κ=1 时,若写成 u′(0)=q0、u′(1)=q1,左右外法向分别为 −1 和 1,所以外向数据是 g(0)=−q0、g(1)=q1,条件成为

∫01f(x)dx=q0−q1.

零源项、零通量允许任意常数解;f=1、两端零通量则根本没有解。兼容性解决的是存在性,固定均值解决的是常数方向造成的不唯一性,二者是不同的问题。在有界连通 Lipschitz 区域、均匀正扩散系数和适当连续右端下,零均值 Poincaré 不等式使兼容问题在零均值子空间上恢复强制性。

弱解与正则性 ​

Lax–Milgram 给出的是所选空间中的解,并不自动把它升级为 H2 解。区间上的常系数 Poisson 方程可以从 u″=−f∈L2 直接得到 H2 正则性;高维中这一步还依赖边界几何和系数正则性。凹角或系数界面可能让弱解存在唯一,而二阶导数仍不平方可积。椭圆内部正则性给出常系数 Laplace 情形的完整局部证明:能量截断控制一阶项,Fourier 逆乘子提升到内部 H2;这项内部估计不包含全局边界正则性。

弱形式也不是唯一的。不同分部积分次数、边界数据处理或试探—检验空间会产生不同 formulation;声称两者等价时必须核对正则性与边界项。

推论与应用

连续弱问题适定并不保证任意离散空间稳定。Galerkin 方法还要选择有限维 trial/test spaces 并验证离散稳定性;有限元方法进一步把这些空间落实为单元、局部自由度、数值积分和稀疏装配。这个次序把连续建模、离散逼近和代数求解分开,避免用一张收敛曲线替代任何一层的条件。

参考资料
  • Lawrence C. Evans, Partial Differential Equations, 2nd ed., AMS, 2010,§6.2.1,Lax–Milgram 与二阶椭圆方程弱解。
  • Alexandre Ern and Jean-Luc Guermond, Finite Elements, Texas A&M MATH 661 Fall 2025 课程讲义,第 25 章 Main results on well-posedness,§25.2 引理 25.2、例 25.4、定义 25.7 及其后的讨论、命题 25.8,分别对应适定性、Poisson 模型、非对称情形与能量极小化。
关系图谱28 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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