Fipy for 1D branched network

Question

I wish to use Fipy for a network made of 1D branches. Does Fipy allow that kind of mesh? I have not been able to find a similar example. Thank you for your help.

Answer 1

FiPy's 1D meshes don't support branched topologies, but you can get the same effective result by concatenating 2D or 3D meshes that are only one cell wide, e.g.,

mesh1 = fp.Grid2D(nx=10, ny=1)
mesh2 = fp.Grid2D(nx=1, ny=7) + [[4.], [1.]]
mesh3 = fp.Grid2D(nx=5, ny=1) + [[5.], [4.]]
mesh4 = fp.Grid2D(nx=4, ny=1) + [[0.], [6.]]
mesh = mesh1 + mesh2 + mesh3 + mesh4