On Aug 20, 2008, at 2:57 PM, franck kalala wrote:
Did you not forgot a term in the second equation?
eqp=ImplicitDiffusionTerm(coeff=1)==_ExplicitSourceTerm(q1s+q2s)
eqphi=TransientTerm(coeff=2)==_ExplicitSourceTerm(q1s-q2s)
after subtracting the two equations of the problem I sent in the
attached file, there a term which looks like 2\nabla(phi\nabla p),
then how to represent it?
In an equation for phi?
PowerLawConvectionTerm(coeff=2 * p.getFaceGrad())
You've also got a p.getFaceGrad().getDivergence() term that isn't in
your .pdf and the sign of your ImlicitDiffusionTerm is wrong in eqp.