是否有解决 Neumann 问题且收敛速度与层数无关的多重网格算法?

计算科学 pde 多重网格
2021-11-30 01:00:46

多重网格方法通常解决水平上的狄利克雷问题(例如点 Jacobi 或 Gauss-Seidel)。当使用连续有限元方法时,组装小的 Neumann 问题比组装小的 Dirichlet 问题要便宜得多。诸如 BDDC(如 FETI-DP)之类的非重叠域分解方法可以解释为解决层次上的“固定”诺依曼问题的多重网格方法。不幸的是,多级 BDDC 的条件数缩放为

C(1+log(Hh))2L

在哪里L是层数和H/h是粗化率。相比之下,基于 Dirichlet 问题的具有平滑器的多重网格方法的条件数具有独立于级别数的条件数。

有没有办法在不失去级别独立性的情况下解决“固定”诺伊曼问题?

1个回答

我不确定这与 BDDC 有什么不同,也没有进行过彻底的分析,但是当我之前阅读它时,这似乎很有趣:

用于大型网格流体模拟的并行多网格 Poisson 求解器