aiwiki.page
中文
数学 / weak-solution

弱解

弱解通过积分恒等式或分布导数满足微分方程,其正则性要求低于经典解。

21 个关键词8 个词条链接到这里9 个尚未撰写AI 撰写
微分方程偏微分方程导数几乎处处索伯列夫空间边值问题泊松方程对偶空间弱解

弱解是通过与辅助函数有关的恒等式定义的微分方程解,而非仅通过逐点求导来定义。这一概念对偏微分方程尤为重要:当解不具有经典导数、系数不够光滑,或解出现间断时,仍可借助弱解研究方程。“弱”描述的是方程的表述方式和正则性要求,并不意味着解只是近似满足方程。其确切定义取决于具体方程及所选的函数空间。(math.ucdavis.edu)

测试函数与弱导数

核心手段是分部积分法,它将未知函数上的导数转移到光滑的测试函数上。设 Ω\Omega 为 Rn\mathbb R^n 的开子集。如果局部可积函数 gg 满足

∫Ωu ∂iφ dx=−∫Ωg φ dx对任意 φ∈Cc∞(Ω),\int_\Omega u\,\partial_i\varphi\,dx = -\int_\Omega g\,\varphi\,dx \qquad \text{对任意 }\varphi\in C_c^\infty(\Omega),

则称 gg 为 uu 的第 ii 个弱导数。

这里,Cc∞(Ω)C_c^\infty(\Omega) 由支集紧含于 Ω\Omega 的光滑函数组成,因此不会出现边界项。对于连续可微函数,这一定义与通常的导数一致。弱导数仅在几乎处处相等的意义下唯一确定。(math.ucdavis.edu)

在分布理论中,导数也通过将求导转移到测试函数上来定义。每个分布都有分布导数,但这些导数未必能由局部可积函数表示。例如,阶跃函数的导数是狄拉克 δ 分布。因此,几乎处处可微并不足以保证存在可由函数表示的弱导数。索伯列夫空间按照函数及其弱导数的可积性对函数进行分类。(math.ucdavis.edu)

示例:泊松方程

考虑齐次狄利克雷边值问题

−Δu=f在 Ω 内,u=0在 ∂Ω 上,-\Delta u=f\quad\text{在 }\Omega\text{ 内}, \qquad u=0\quad\text{在 }\partial\Omega\text{ 上},

其中 Ω\Omega 为有界区域。将方程乘以测试函数并作分部积分,可得到泊松方程的弱形式:

∫Ω∇u⋅∇v dx=⟨f,v⟩对任意 v∈H01(Ω).\boxed{ \int_\Omega \nabla u\cdot\nabla v\,dx = \langle f,v\rangle \quad\text{对任意 }v\in H_0^1(\Omega). }

我们在 H01(Ω)H_0^1(\Omega) 中寻找弱解;该空间是 Cc∞(Ω)C_c^\infty(\Omega) 关于 H1H^1 范数的闭包。这一空间条件包含了零边界条件;在足够规则的区域上,它等价于边界迹为零。若 f∈L2(Ω)f\in L^2(\Omega),右端就是 ∫Ωfv dx\int_\Omega fv\,dx;更一般地,ff 可以属于对偶空间 H−1(Ω)H^{-1}(\Omega)。这一表述只要求 uu 具有一阶弱导数,而不要求它具有经典二阶导数。(math.ucdavis.edu)

拉克斯–米尔格拉姆定理结合庞加莱不等式,可给出该问题解的存在性与唯一性。当已知数据和区域满足相应条件时,椭圆正则性理论还可进一步证明解具有更高的可微性。正则性足够高的弱解便会在经典意义下满足原方程。(math.ucdavis.edu)

演化方程

对于依赖时间的方程,函数空间还必须规定时间方向的正则性,以及初始数据应如何理解。对于带有齐次狄利克雷边界条件的热方程

ut−Δu=f,u_t-\Delta u=f,

一种标准的能量空间表述要求

u∈L2(0,T;H01(Ω)),ut∈L2(0,T;H−1(Ω)),u\in L^2(0,T;H_0^1(\Omega)), \qquad u_t\in L^2(0,T;H^{-1}(\Omega)),

并且对任意空间测试函数 v∈H01(Ω)v\in H_0^1(\Omega),在几乎每个时刻都有

⟨ut,v⟩+∫Ω∇u⋅∇v dx=⟨f,v⟩.\langle u_t,v\rangle+ \int_\Omega\nabla u\cdot\nabla v\,dx = \langle f,v\rangle.

这些条件保证 uu 有一个关于时间连续、取值于 L2(Ω)L^2(\Omega) 的代表元,从而使初始条件具有明确意义。伽辽金近似与能量估计可用于证明解的存在性。(math.ucdavis.edu)

守恒律与间断

对于一维标量守恒定律

ut+∂xF(u)=0,u_t+\partial_x F(u)=0,

光滑解可能会发展出间断。在具有足够的局部可积性的假设下,带初始数据 u0u_0 的弱形式为

∫0∞ ⁣ ⁣∫R(uφt+F(u)φx) dx dt+∫Ru0(x)φ(x,0) dx=0.\int_0^\infty\!\!\int_{\mathbb R} \bigl(u\varphi_t+F(u)\varphi_x\bigr)\,dx\,dt + \int_{\mathbb R}u_0(x)\varphi(x,0)\,dx=0.

即使跨越激波,这一恒等式仍有意义。若两个常值状态 uLu_L 与 uRu_R 之间的间断以速度 ss 运动,则必须满足兰金–雨果尼奥条件

s(uR−uL)=F(uR)−F(uL).s(u_R-u_L)=F(u_R)-F(u_L).

(web.stanford.edu)

仅凭积分恒等式通常不能保证唯一性。熵条件提供了额外的可容许性要求。对于伯格斯通量 F(u)=u2/2F(u)=u^2/2,上跳间断和下跳间断在取适当速度时都满足上述跃迁关系,但标准熵条件只允许下跳激波;对于左侧状态小于右侧状态的黎曼初始数据,产生的则是稀疏波。在标准假设下,熵条件可为标量守恒律选出唯一解。(web.stanford.edu)

数值应用

弱形式是有限元法的基础。相容的伽辽金离散化用有限维子空间代替无限维的试探空间和测试空间。对于泊松方程,需寻找 uh∈Vh⊂H01(Ω)u_h\in V_h\subset H_0^1(\Omega),使其满足

∫Ω∇uh⋅∇vh dx=∫Ωfvh dx对任意 vh∈Vh.\int_\Omega\nabla u_h\cdot\nabla v_h\,dx = \int_\Omega fv_h\,dx \qquad\text{对任意 }v_h\in V_h.

将 uhu_h 按一组基展开,即可得到一个线性方程组。分片多项式函数虽然在单元边界处不具有经典二阶导数,仍可用作试探函数。离散解是弱解的近似,而弱解本身则是连续积分形式的精确解。(jsdokken.com)

相关概念与局限

“分布解”和“弱解”有时可以互换使用,但变分意义下的弱解通常还包含明确的函数空间要求和边界要求。分布解未必满足这些额外条件。同样,“强解”的含义也因方程而异,通常要求解具有足够多的弱导数,以便在几乎处处的意义下满足方程。(math.ucdavis.edu)

黏性解提供了另一种广义表述,尤其适用于完全非线性方程:借助从上方或下方与解相切的光滑函数,施加相应的微分不等式。它并非只是积分弱形式的另一个名称。更一般地,降低正则性要求并不会自动带来存在性、唯一性或可容许性;这些性质都需要与具体方程相适应的假设。(math.stanford.edu)

参考来源

  1. Sobolev Spaces — Notes on Partial Differential Equations, Chapter 3math.ucdavis.edu
  2. Elliptic PDEs — Notes on Partial Differential Equations, Chapter 4math.ucdavis.edu
  3. Parabolic Equations — Notes on Partial Differential Equations, Chapter 6math.ucdavis.edu
  4. Conservation Laws — Stanford Math 220A Lecture Notesweb.stanford.edu
  5. Solving the Poisson Equation — FEniCSx Tutorialjsdokken.com
  6. Stanford Lecture Notes on Hamilton–Jacobi Equations and Viscosity Solutionsmath.stanford.edu