如何推导有限元法偏微分方程的弱公式?

计算科学 pde 有限元
2021-12-01 23:13:31

我已经对有限元方法进行了基本介绍,它并没有强调对“弱公式”的复杂理解。我知道使用 galerkin 方法,我们将(椭圆)PDE 的两侧乘以一个测试函数,然后积分(通过部分或发散定理)。有时,在得出适当的弱公式之前,我需要按部分积分两次(根据书后的答案)。但是当我尝试将相同的概念应用于其他 PDE 时(假设它们仍然与时间无关),我似乎无法识别该公式何时适合离散化。是否有任何“危险信号”可以告诉我这种形式可以离散化为线性方程组?

此外,如何选择一组合适的基函数?

1个回答

问自己以下问题:

首先,分部积分如何影响问题的可解性,以及解的空间?

其次,您可以为哪个函数空间构建一系列可以实现的子空间(ansatz 函数)?

让我们考虑泊松问题u=f对于,比如说,在上,具有齐次狄利克雷边界条件。通过积分,方程的左边和右边可以被认为是上的有界泛函,比如对于我们有fL2[0,1]L2ϕL2

ϕuϕdxϕfϕdx

由于中的任何函数都可以,因此如果您只知道所有测试函数的值,则两个积分泛函都是完全已知的。但是使用测试功能,您可以进行分部集成,并将左侧转换为功能L2L2

ϕuϕdx

将其读作:“我采用一个测试函数,计算它的微分,并将它与 -u' 在 [0,1] 上进行积分,然后返回结果。” 但是这个函数没有定义和限制在上,因为你不能取任意函数的微分。一般来说,它们可能看起来非常奇怪。ϕL2L2

我们仍然观察到,这个泛函可以扩展到 Sobolev 空间,它甚至是上的有界泛函。这意味着,给定 -norm 的倍数粗略估计值。此外,函数当然不仅在上定义和有界,而且在上定义和有界。H1H01ϕH01uϕdxH01ϕϕfϕdxL2H01

现在您可以,例如,应用 Lax-Milgram 引理,因为它在任何 PDE 书中都有介绍。也有描述它的有限元书,只有功能分析,例如 Ciarlet 的经典著作,或 Braess 的相当新的书。

Lax-Milgram 引理为 PDE 人提供了一个很好的纯分析工具,但他们也使用了很多陌生的工具来实现他们的目的。尽管如此,这些工具也与数值分析相关,因为您实际上可以为这些空间构建离散化。

例如,为了具有的离散子空间,只需使用帽子函数。它们没有跳跃并且是分段可微的。它们的微分是分段常数向量场。这种构造适用于,这很好,但是你能想出一个 ansatz 空间,它的函数不仅有梯度(很好,即平方可积),而且谁的梯度又发散了?(再次,平方可积)。一般来说,这很难。H01d=1,2,3,...

因此,一般而言,您如何构建弱公式的原因是您想要应用 Lax-Milgram 引理,并拥有一个使得函数实际上可以实现的公式。(作为记录,Lax-Milgram 不是该上下文中的最后一个词, ansatz 空间也不是离散化中的最后一个词,例如,参见不连续 Galerkin 方法。)H01

对于混合边界条件的情况,自然测试空间可能与您的搜索空间不同(在分析设置中),但如果不参考分布理论,我不知道如何描述,所以我停在这里。我希望这是有帮助的。