1 次多项式的 FEM 振荡

计算科学 有限元 有限差分
2021-12-16 13:20:01

我有以下椭圆一维问题

μu+βu=1
u(0)=u(1)=1
其中对于这个特定问题,我使用以下空间步骤μ=10e5β=1h=[0.1,0.01]

令我惊讶的是,这个问题并没有出现振荡。我认为这与 Péclet 数有关。当问题的源项(右侧函数)等于 0 时,我知道 Péclet 数应该小于 1 以避免这种振荡。Pe

这种缺乏振荡是否与我的问题的本地 Péclet 数有关?

我尝试使用有限差分来查看是否可以找到为什么这不会出现振荡,但我不知道如何从这里开始

(Pe1)ui+1+2ui(Pe+1)ui1=h2μ

关于如何进行或我是否真的朝着正确的方向有任何想法?

1个回答

例如,Wesseling,计算流体动力学原理,Springer 2001 中介绍了有限差分/体积情况。主要思想是采用形式的解,填写并查看它是否允许解以使交替出现符号(摆动)。包括边界条件也很重要。ui=A+Bziz<0zi