mesh - FiPy の 1D 球面グリッド

Question

1D グリッド上の球座標で FiPy の拡散方程式を解きたいと思います。また、左の境界を r=0 ではなく r=0.1 にしたいと考えています。

1D 球対称のモジュールが見つかりません。円筒形のみ。私は Grid1D でそれを行い、単純に del^2 演算子を球面座標で記述し、次に r^2 を乗算します (ここで説明したように)。ただし、境界の位置を指定する方法はまだわかりません。

誰かがこれを行う方法を教えてもらえますか? どうもありがとう。