用 Dirichlet-Neumann 的混合 BC 模拟 3-D 拉普拉斯方程

计算科学 边界条件 泊松 椭圆pde
2021-12-12 19:18:36

我在单位 3-D 立方体上使用有限差分法模拟混合边界值 (Dirichlet-Neumann) 问题,使得左、下和前平面具有u=u(x,y,z)=1(Dirichlet) 和右、上、背平面有un=0. 我知道真正的解决方案是1无处不在,因此我使用了停止条件作为errmax<0.00001我发现像这样的错误。当输入为时问题收敛,但当输入为时,最大误差不会降低到以下。同样,当问题大小为时,误差不会降低到以下。我的问题是:这是预期的行为吗? err=1.0ucurrent16×16×1632×32×3264×64×640.000019128×128×1280.000040

1个回答

不,这不是预期的。您的迭代需要收敛。由于您没有说您的迭代是什么,因此很难说我错了什么。但是很容易说你一定有一个bug:只要你的迭代是任何合理的方法,它就需要收敛!