我想求解一维网格上球坐标系中的FiPy中的扩散方程。我也希望左边界为r = 0.1,而不是r = 0。
我找不到一维球对称的模块;仅圆柱形。我想我用Grid1D做到了,只需在球坐标中编写del ^ 2运算符,然后乘以r ^ 2(如here所述)。但是,我仍然不知道如何指定边界的位置。
有人可以建议我该怎么做吗?非常感谢。
答案 0 :(得分: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上镜像的球对称网格将是一个受欢迎的请求请求。