如果没有物理边界上的狄利克雷条件,交替 Schwarz 不会收敛

计算科学 有限元 域分解
2021-12-17 13:50:50

注意:感谢评论,我意识到我有两个问题,每个问题都可以自己更清楚地描述。这个修订后的问题涵盖了第一个问题。


我想解决

Δu1=0,u(0,y)=0,u(1,y)=0
使用交替施瓦茨在一个矩形上。以下是我的网格,其中小圆盘表示子域的内部边界节点(具有我更新的狄利克雷条件的节点):

x=0x=1是全网格的物理边界。局部解不会收敛到全局解。问题大多出在网格的重叠部分,这里是实际解与局部解给出的解之间的差异:

数学图形

如果我施加狄利克雷条件u(x,0)=0u(x,1/ϕ)=0, 在哪里1/ϕ是网格的高度,则局部解收敛于全局解。

它不应该只适用于我之前提到的边界条件吗?

编辑:我通过以这种方式选择人工边界来解决它:

数学图形

0个回答
没有发现任何回复~