椭圆方程的稳定参数

计算科学 有限元 椭圆pde
2021-12-18 17:16:48

我只是想求解椭圆方程:

κ2u+u=f
在哪里f[0,1]. 当使用带有拉格朗日元素的连续 Galerkin 时,我注意到κ必须大于与网格分辨率相关的某个值以避免u小于 0 或大于 1。我想避免这些过冲/下冲,而不必增加网格分辨率。有人建议我使用有限体积来求解这个方程,从而解决了这个问题。然而,它使用非笛卡尔二阶网格的实现是繁重的。因此我决定改用稳定的方法:
(ku,v)+(u,v)+(αh2u,v)=(f,v)
考虑到我人为地增加了扩散系数,这似乎可行,这将随着细化而消失。但是,我还没有找到关于适当值的文献α. 有人可以给我一些指示吗?谢谢。

1个回答

你想解决的问题,用一个小κ是“奇异扰动的”,即它通常具有边界层和内部层。你能想到的方式是,如果κ为零,那么你会得到u=f,并且对于κ>0, 你得到一个平滑的版本f带有平滑长度刻度L这与κ在某种程度上——我相信它是Lκ.

如果你有那个网格尺寸,这很好用h大约等于或小于L, 但如果h>L,那么您实际上是在尝试解决不连续的函数。对于这些类型的函数,除非您使用(i)是非线性的或(ii)具有足够的人工扩散以使物理扩散产生的长度尺度的方案,否则您将得到吉布斯现象κ加上人工扩散至少与网格尺寸相同。

有限体积方法(基本上是分段常数“不连续 Galerkin”方案)和添加αh2κ落入第二类。另一种方法是使用非线性冲击捕获方法,尽管这会使问题的解决变得更加复杂。