非线性偏微分方程多重网格法中的平滑算子问题

计算科学 线性代数 pde 多重网格 迭代法 参考请求
2021-12-21 20:02:29

假设我们正在处理一个非线性问题,比如说

Au:=Lu+G(u)=f
的非线性是多项式类型,即是线性算子,并且 ,或者更一般地说,我们有一个 Fréchet 导数界限,如: 其中D是算子A的 Fréchet 导数。ALG(u)=uk
A(u)A(v)DA(u)(uv)CA(u)A(v)
DA

我们想使用一些数值方法和迭代求解器来获得这个问题的有限维近似。对于两个网格循环的多重网格,根据Saad 第 13 章所写的书,在我们计算更精细级别的网格尺寸h(与粗级别网格尺寸2h )的残差之前,我们想应用一个预平滑算子\mathcal {S}到更精细级别的初始猜测u^h_0 :2hSu0h

  1. Presmooth: uh:=S(Ah,u0h,fh)
  2. 计算残差:rh=fhAhuh

然后在接下来的几个步骤中,我们可以返回粗网格2h以在更精细的级别求解和校正uh这里Ah是算子A的有限维线性化矩阵。

现在我的问题是:为非线性多重网格方法选择“好的”预平滑算子的标准是什么?以及如何构造这个 presmooth 算子?有规律可循吗?感谢您是否可以与我分享一些具有与此类问题类似的示例的参考资料。

1个回答

Brandt 的书Multigrid technologies对 Full Approximation Scheme 有很好的讨论,1982 Guide to Multigrid Development也有一些讨论。与线性问题一样,要求平滑器去除残差的高频分量。形式谱分析只适用于线性化问题,而 FAS 理论只需要局部线性化。Trottenberg、Oosterlee 和 Schueller 于 2001 年出版的关于多重网格的书也是一本很好的参考书,其中包含非线性问题的局部傅里叶分析。