1D сферическая сетка в FiPy

Question

Я хотел бы решить уравнение диффузии в FiPy в сферических координатах на одномерной сетке. Я также хотел бы, чтобы левая граница была при r = 0,1, а не r = 0.

Я не могу найти модуль для 1D сферической симметрии; только цилиндрический. Я полагаю, что я делаю это с помощью Grid1D и просто пишу оператор del ^ 2 в сферических координатах, затем умножаем на r ^ 2 (как упоминалось здесь ). Однако я до сих пор не знаю, как указать расположение границ.

Может кто-нибудь посоветовать мне, как это сделать? Большое спасибо.

Answer 1

Все фипи-сетки могут быть смещены на вектор соответствующего размера, например,

>>> 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 , будет желательным запросом на извлечение.