Thanks that right I have not subversion how can I deal with "&" or "and" in fipy?
Regards 2008/8/20 Jonathan Guyer <[EMAIL PROTECTED]> > > > 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. > > > -- *********************************************** *********************************************** Franck Kalala Mutombo [EMAIL PROTECTED] African Institut for Mathematical Sciences, Muizenberg Cape Town, South Africa [EMAIL PROTECTED] ************************************************ ************************************************
