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.


Reply via email to