扩散-传输问题 FEM

计算科学 有限元 误差估计
2021-12-02 05:27:41

我正在看一本关于扩散传输问题的 FEM 书。

div(μu)+bu+γu=fin Ω
u=0in Ω (in the boundaries)

它说如果|b|μ1那么问题就是以运输为主的问题。

取上一章的一些东西,更准确地说是稳定性和收敛性的 Galerking 分析,以及误差的近似,我们有

Mμ+|b|continuity constant
α=μcoercivity constant

然后它有

Mα1+|b|μ1

因此,得出的结论是,误差的估计是

||uuh||CMαhr|u|Hr+1(Ω)

这告诉我们,如果空间步长,Galerkin 方法可能会给出不令人满意的结果h不够小。

我不明白它怎么能得出这样的结论。任何想法?请不要给出复杂的解释,我只是想理解这些 FEM 的东西,我是一名计算机科学专业的学生,​​而不是一名数学专业的学生。

1个回答

如果您按照分析直到最后一个方程,那么您可以看到问题。你写:

||uuh||CMαhr|u|Hr+1(Ω)

从上面的表达式中也可以写成

||uuh||C(1+|b|μ)hr|u|Hr+1(Ω)

现在假设您要解决固定网格大小的问题。在第一个问题|b|=μ,那么你有

||uuh||2Chr|u|Hr+1(Ω)

现在你转到下一个问题,与扩散相比,传输更强,|b|=104μ.

||uuh||104Chr|u|Hr+1(Ω)

因此,在固定尺寸的网格上,第二种解决方案要差几个数量级。改善它的唯一方法是减少h到解决运输规模的数量(抵消巨大的常数)。