On May 7, 2010, at 1:19 PM, Mike Waters wrote:
>
> I am trying to write a kinetic equation roughly in the form:
>
> \frac{\partial\phi}{\partial t}= \alpha\nabla^{2}\phi + \beta\nabla^{2}T ,
>
> Where both T and phi are cell variables. I'm not sure how to describe
> the Laplacian of T in the kinetic equation.
>
>
> What happens if I use this:
>
> TransientTerm() ==DiffusionTerm(alpha) + (beta *
> T.getFaceGrad()).getDivergence() ?
Yep, that's exactly what we suggest in the FAQ:
http://www.ctcms.nist.gov/fipy/documentation/FAQ.html#how-do-i-represent-a-diffusive-source
The only caveat is that this won't see any boundary conditions on T.
