On Apr 27, 2006, at 3:51 PM, Buddy Damm wrote:
Your convection *coefficient* (not equation!!!) is then:
Mc * yva * d2Gmducdphi * phi.getFaceGrad()
should this be (in my case of factoring out uc)
convTerm = (1/uc)*Mc * yva * d2Gmducdphi # added the (1/uc)
Why are you adding the (1/uc)?
Neither Daniel nor I remember telling you to divide by uc, but if we
did, it's not what we meant. Factor it out.
Your convection term is \nabla\cdot(uc * something)
The convection coeff is 'something'. Just 'something'. Not
'something'/uc or 'something' * uc. Just 'something'.
Other than that, this looks fine, although I don't see any advantage
to separately writing convTerm and phaseVelocity, especially when
convTerm is not a `Term` in FiPy's sense of the word.
