数值查找格林函数

计算科学 有限元 有限差分
2021-12-02 03:47:40

我正在尝试对这个方程的格林函数进行数值评估:

[2x2+f(x)]G(x)=δ(xx0)

具有 Neumann 边界条件。这里,函数f(x)是已知的。解决这个问题的难点在于数值近似δ功能。而不是近似δ1x我试图解决δ通过在周围的一个小区间内整合双方x0,但这给我留下了关于导数的条件Gx0我不知道如何在有限元方法中实现。

我真的很感激求解此类方程的参考(最好使用有限元方法),或者是否有人可以为求解此类方程提供指导。

1个回答

为什么要近似 delta 函数?在有限元法中,您需要评估方程的弱形式,其右侧将如下所示:

φi(x)δ(xx0),
在哪里φi是个i形状函数。但是你知道这个积分是什么:
φi(x)δ(xx0)=φi(x0).
当然,大多数形状函数在x0,但在包含的单元格上定义的形状函数x0不会,这些将是您右侧向量中唯一的非零条目。