On Fri, Jan 17, 2020 at 5:31 AM Pavel Aleynikov
wrote:
> > sort of ratio so the absolute value of the volume won't matter.
> Unless you have a source with a defined integral absolute value. But I agree
> that in general "theta = 1" is not a problem, as long as it's clear from the
>
Ok. I understand. I see the logic of saving on some operations.
> Argulably, these quantities should be calculated in some
> sort of ratio so the absolute value of the volume won't matter.
Unless you have a source with a defined integral absolute value. But I agree
that in general "theta = 1" is
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