我只是想求解椭圆方程:
在哪里. 当使用带有拉格朗日元素的连续 Galerkin 时,我注意到必须大于与网格分辨率相关的某个值以避免小于 0 或大于 1。我想避免这些过冲/下冲,而不必增加网格分辨率。有人建议我使用有限体积来求解这个方程,从而解决了这个问题。然而,它使用非笛卡尔二阶网格的实现是繁重的。因此我决定改用稳定的方法:
考虑到我人为地增加了扩散系数,这似乎可行,这将随着细化而消失。但是,我还没有找到关于适当值的文献. 有人可以给我一些指示吗?谢谢。
我只是想求解椭圆方程:
你想解决的问题,用一个小是“奇异扰动的”,即它通常具有边界层和内部层。你能想到的方式是,如果为零,那么你会得到,并且对于, 你得到一个平滑的版本带有平滑长度刻度这与在某种程度上——我相信它是.
如果你有那个网格尺寸,这很好用大约等于或小于, 但如果,那么您实际上是在尝试解决不连续的函数。对于这些类型的函数,除非您使用(i)是非线性的或(ii)具有足够的人工扩散以使物理扩散产生的长度尺度的方案,否则您将得到吉布斯现象加上人工扩散至少与网格尺寸相同。
有限体积方法(基本上是分段常数“不连续 Galerkin”方案)和添加到落入第二类。另一种方法是使用非线性冲击捕获方法,尽管这会使问题的解决变得更加复杂。