This is where I don't understand because doing ImplicitSourceTerm(coeff=(phi+1))
is like doing ImplicitSourceTerm(coeff=2) because I initially defined phi as phi = CellVariable(mesh=mesh, value=1) How should I define phi so that the source term is not a constant? ________________________________________ From: fipy-boun...@nist.gov <fipy-boun...@nist.gov> on behalf of Guyer, Jonathan E. Dr. (Fed) <jonathan.gu...@nist.gov> Sent: Friday, April 22, 2016 1:22:03 PM To: FIPY Subject: Re: Diffusion equation with additional terms -ImplicitSourceTerm(coeff=(phi+1)) will do it, although you may want to look at http://www.ctcms.nist.gov/fipy/examples/phase/generated/examples.phase.simple.html for a discussion of the most efficient way to linearize a source. > On Apr 22, 2016, at 12:22 PM, Gaury, Benoit H. (IntlAssoc) > <benoit.ga...@nist.gov> wrote: > > Hi everyone, > > I want to solve an equation of the form > > \nabla^2 phi - phi^2 - phi = 0, > > where phi is the unknown we solve for. I found out that an additional linear > term can be added do the stationary diffusion equation using > ImplicitSourceTerm(). However I didn't find anything to add phi^2 + phi. Did > I miss something? > > Thanks, > Benoit > > _______________________________________________ > fipy mailing list > fipy@nist.gov > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]