You are right to use a convection term, but make the power of phi one less. phi**2 * theta.getFaceGrad() will effectively be the velocity coefficient in the convection term.
Cheers On Fri, Sep 12, 2008 at 7:34 AM, franck kalala <[EMAIL PROTECTED]> wrote: > 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] > ************************************************ > ************************************************ > -- Daniel Wheeler
