Kevin -
You are correct that there's no way to project a FaceVariable back to a
CellVariable. Since this is only for the initial condition, I would just
explicitly set it, e.g.,
phi.value = (1 + x) / (gamma * psi.grad[0])
As for the boundary condition, I think FiPy's default no-flux boundar
Hi,
I have two questions.
I have a complicated set of 1D equations. The diffusion coefficient of the
1st is dependent on the derivative of the 2nd's variable, and the 2nd has a
source term of the first. In addition, the initial condition of one of the
variables is dependent on this diffusion coef