Dear All, How do I set a viewer in fipy?
Regards 2008/9/12 Daniel Wheeler <[EMAIL PROTECTED]> > > 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 > > -- *********************************************** *********************************************** Franck Kalala Mutombo [EMAIL PROTECTED] African Institut for Mathematical Sciences, Muizenberg Cape Town, South Africa [EMAIL PROTECTED] ************************************************ ************************************************
