非重叠域分解

计算科学 域分解
2021-12-09 00:30:19

我通过重叠域分解解决了一个简单的测试示例。问题域是一个分解为两个域的矩形。

所以交集边界上的值在第一时间被猜测出来,并在每一步更新。

现在我想通过非重叠的方法来解决这个问题。但是当我使用这个程序时,所有步骤的答案都是一样的。我错了吗?或者我应该以不同的方式解决这个问题?

请提示我。

1个回答

简短回答:是的,您必须使用不同的方法,例如Neumann-Neumann 方法一个很好的参考是Widlund 的书非重叠方法基于以下原则,如果u求解域中的泊松方程Ω已划分为两个域Ω1,Ω2, 在共享边界Γ两个域中,u它的梯度必须匹配Γ.

相比之下,对于重叠 DD 方法,域的交集Ω1Ω2有一个非空的内部。您可以使用的内部值uΩ1作为狄利克雷边界条件来求解泊松方程Ω2,然后使用内部值Ω2更新解决方案Ω1,等等。如果您尝试将重叠方法应用于实际上不重叠的域,则每个子域上的解决方案都没有分歧的余地,因此您都不会更新。

更重要的是,重叠 DD 方法的收敛速度如下H/h, 在哪里H是重叠的宽度和h是网格宽度*。有一个权衡:你可以有一个非常小的重叠和廉价的子域解决方案,但该方法不会快速收敛;或者您可以与更昂贵的子域求解有很大的重叠,但整个方法收敛得更快。

*除非您使用全局粗略求解,但现在我们离题了。