Dear all, I am solving an equation of this type
phi_t=div(phi**3*div\theta) in the right hand side of this equation phi is to the power 3. in the manual it is explain how to deal with this term when phi is to the power 1. but here I don't know how to represent it in fipy the second hand side, I try this PowerLowerConvectionTerm((phi**3).getArithmeticFaceValue()*theta.getFaceGrad())) but with tha I don't have expected result for my codes. Is there any other better way to represent such term? Regards 2008/9/3 Daniel Wheeler <[EMAIL PROTECTED]> > > On Wed, Sep 3, 2008 at 7:07 AM, franck kalala > <[EMAIL PROTECTED]> wrote: > > In short > > > > I have an equation of this kind in my sets my equations > > > > > > -\nabla\nabla p=nabla(D1(phi)\nabla D2(phi))+q1(phi)+q2(phi) > > > > D1(phi)=1-phi > > > > D2(phi)=1-phi > > > > > > the unknown is p how do the second term of the above equation? > > The right hand side is just a source term. Write it in fipy as you > would write it mathematically. > > Cheers > > -- > Daniel Wheeler > > -- *********************************************** *********************************************** Franck Kalala Mutombo [EMAIL PROTECTED] African Institut for Mathematical Sciences, Muizenberg Cape Town, South Africa [EMAIL PROTECTED] ************************************************ ************************************************
