我正在使用 Galerkin 的方法为我手头的问题(1D/2D)设置一些 pde 求解器。
确实存在一些突然的材质界面,到目前为止,我对此所做的只是:
1) 界面上没有元素,即节点正好位于界面上。
2)没有对接口施加特殊的BC,只是自然的BC(原样)。
我想知道:
1)我所做的足以处理接口吗?
2)有人建议沿界面使用零厚度元素,这样界面上的每个节点都关联一个重复节点,有必要吗?理由是什么?
我正在使用 Galerkin 的方法为我手头的问题(1D/2D)设置一些 pde 求解器。
确实存在一些突然的材质界面,到目前为止,我对此所做的只是:
1) 界面上没有元素,即节点正好位于界面上。
2)没有对接口施加特殊的BC,只是自然的BC(原样)。
我想知道:
1)我所做的足以处理接口吗?
2)有人建议沿界面使用零厚度元素,这样界面上的每个节点都关联一个重复节点,有必要吗?理由是什么?
这取决于你的问题的物理性质。让我们以二阶椭圆方程为例,其中和在界面上是不连续的。
在最常见的情况下,在整个界面中,您和,因此您的无所事事和组装(界面上的唯一自由度)策略符合到物理学。
但是在某些情况下,物理是不同的。例如,在电磁中,可能存在表面(界面)电荷,因此对于横向场,您会得到一个跳跃其中和是来自的轨迹,用于和之间的接口,而在接口之间保持连续。
更奇怪的是,对于接口正常的字段组件,您会得到不同的东西:和这些是非常不寻常的。因此,人们应该要求的不是数学,而是物理。
对于椭圆 PDE 的材料属性更改,您不必做任何特别的事情。最好将节点放在接口上,但您甚至不需要强加 BC。这不是一个边界。只是系数的变化。