I'm back again :P
I have two questions:
(1) If I define a mesh Grid1D(nx=30, dx=5) and solve a CellVariable
using this mesh, the CellVariable.phi value would be an array of 30.
Does phi.value[0] indicate the value at r = 0? or r = 5? If it is r =
0, doesn't that make the last value in the array r = (nx-1)*dx?
(2) If I set the boundary condition of a CellVariable phi to
BCs = (FixedValue(faces=mesh.getFacesRight(), value = 100.),
FixedFlux(faces=mesh.getFacesLeft(),value = 0.))
I'm modeling this for a sphere but assuming uniform diffusion (so it is
a 1D-diffusion problem). Hopefully by declaring this I mean that at r =
0 the core, the flux is 0. But at r = R which is the surface of the
sphere, the value of phi should be 100.
I noticed this problem because when I set my right boundary to some
value C, and then look at the solved variable, the last value in the
array is not C, but slightly lower than it.
Elizabeth
