Skip to content

变分形式与弱问题

Variational formulation · Weak formulation

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

从强形式到弱形式

以有界 Lipschitz 区域 ΩRd 上的 Poisson 问题为主线:

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

强形式要求微分方程与边界条件具有逐点或经典迹意义。若先假设 u 足够光滑,将方程乘以检验函数 v 并分部积分,可得

Ω(Δu)vdx=ΩuvdxΩunvds.

V=H01(Ω) 后,v 在边界上的迹为零,边界项消失。弱问题写成:求 uV,使

a(u,v)=(v)(vV),

其中

a(w,v)=Ωwvdx,(v)=Ωfvdx.

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

边界条件进入哪里

(κu)+cu=f

若边界分成 ΓDΓN,并给定

u=gD on ΓD,κun=gN on ΓN,

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

UgD={wH1(Ω):trw=gD on ΓD},V={vH1(Ω):trv=0 on ΓD}.

双线性型与右端为

a(w,v)=Ωκwvdx+Ωcwvdx,(v)=Ωfvdx+ΓNgNvds.

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

要放入标准适定性框架,还需说明数据和系数条件。例如可要求 ΓD 具有正表面测度,gD 是某个 H1 函数的可容许迹,并且

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

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

抽象适定性

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

|a(w,v)|MwVvV,a(v,v)αvV2,α>0.

若线性泛函 连续,Lax–Milgram 定理保证弱问题存在唯一解,并有

uV1αV.

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

直觉

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

能量最小化只适用于对称强制问题。此时

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

的唯一极小点是弱解。对流占优方程、鞍点系统或一般非对称问题未必存在这种图像;此时真正控制适定性的可能是 inf–sup 条件。

例子与边界

纯 Neumann 问题中,常数函数使 a(v,v)=0。这不是离散矩阵装配错误,而是连续问题本来只在模去常数后唯一。必须先检查总源项与边界通量的兼容条件,再决定固定均值还是使用商空间。

弱形式也不是唯一的。不同分部积分次数、边界数据处理或试探—检验空间会产生不同 formulation;声称两者等价时必须核对正则性与边界项。连续问题适定也不保证任意离散空间稳定,后续离散由Galerkin 方法单独处理。

参考资料
  • Lawrence C. Evans, Partial Differential Equations, 2nd ed., AMS, 2010,weak formulations and Sobolev spaces。
  • Susanne C. Brenner and L. Ridgway Scott, The Mathematical Theory of Finite Element Methods, 3rd ed., Springer, 2008,variational problems。