我想解决两种简单的二维问题，即像这样的 L 形几何上的固定热方程：

并计算以下几何形状的气隙中的趋磁场：

由于重新发明轮子并不聪明，因此我搜索了一个可以处理此类问题的开源库。据我所知，在 Python、Octave、Scilab 或 R 中都没有可以处理这种类型的几何图形的有限差分 (FD) 库，它不是矩形的，但仍然应该适合矩形网格（如果我记得从 Incropera 的“热传递”中正确读取，内部和外部节点可以由 FD 处理）

这种缺乏可用性让我感到奇怪。Python、Octave、Scilab 和 R 都具有用于求解 ODE 的内置库，但为什么没有用于 PDE 的 FD？是否得出结论认为专注于有限元 (FE) 就足够/更好了？我想继续使用 FD 有两个原因：简单，并且我希望将来解决其他类型的 PDE，这可能在 FE 求解器上不可用。

考虑开发这样一个库是否有意义？FD（对简单几何形状）的限制不值得努力吗？

提前谢谢了！