mesh - FiPy 中的一维球形网格

Question

我想在一维网格上的球坐标中求解 FiPy 中的扩散方程。我还希望左边界位于 r=0.1，而不是 r=0。

我找不到一维球面对称的模块；只有圆柱形。我想我是用 Grid1D 来做的，只需在球坐标中编写 del^2 运算符，然后乘以 r^2（如此处所述）。但是，我仍然不知道如何指定边界的位置。

有人可以建议我怎么做吗？非常感谢。

Answer 1

所有 fipy 网格都可以通过适当维度的向量进行偏移，例如，

>>> m = fp.Grid1D(nx=10, dx=.1) + [[1.5]]
>>> print m.x
[1.55 1.65 1.75 1.85 1.95 2.05 2.15 2.25 2.35 2.45]

在CylindricalUniformGrid1D上镜像的球对称网格将是一个受欢迎的拉取请求。