具有纯 Neumann 边界条件的泊松方程(使用 FEM)

计算科学 有限元 芬尼克斯
2021-12-04 08:16:55

我试图推导出具有纯诺依曼边界条件的泊松方程的正确变分形式,以及一个额外的约束Ωudx=0如此链接中所述。

据我了解,问题可以定义为:

寻找(u,c)V×R这样

F(u,c)=12uufuΩgu+cu

被最小化了uc

为了最小化相对于u,我开始写

ϵF(u+ϵv,c)|ϵ=0=0

这导致

uv+cv=fv+Ωgv.

关于c这给了u=0.

只需添加这两个等式就会导致链接中提到的变分形式(其中du=0对所有人dR确保两个方程分别保持正确)。到目前为止这是正确的吗?

现在,在有限元描述中,我们还有一个未知数,c,即拉格朗日乘数,不是吗?要进入一组方程,方程的变分形式必须对所有方程都成立N测试函数(形状函数)v. 我是否正确,这将导致N×(N+1)矩阵方程?我是否也正确,因为表达式应该对所有人都有效dR我们可以在这些表格的实际评估中选择一个随机值吗?

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