On Wed, Jan 15, 2020 at 4:21 AM Pavel Aleynikov
wrote:
>
> 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)
>
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