Hi, How are "mesh.cellVolumes" defined in fp.CylindricalGrid1D case? The surface of a circle grid is not equal pi*r^2.
import fipy as fp L = 1.; nx = 1000 mesh = fp.CylindricalGrid1D(dx=L/nx, nx=nx) print(mesh.cellVolumes.sum()) >> 0.5 Why not pi? How should I integrate Variables on such a grid? (2*pi*Var*mesh.cellVolumes).sum()? Thank you, Pavel PS. The solution of a diffusion equation on the CylindricalGrid1D seems fine, i.e. it converges to Bessel J_0 eigenvalue fairly well. PPS. I've accidentally sent this message from an unregistered email first and the cancelation link in the autoreply does not work. So, sorry if this is the duplicate. _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]