Re: Taking Dot Product of faceGrad

2013-12-19 Thread Daniel Wheeler
On Wed, Dec 18, 2013 at 1:37 PM, Jane Hung jyh...@mit.edu wrote: OK, I can write the equations with grad instead of faceGrad, but it doesn't seem like the equations can be solved with what I've tried. I tried to use eq.solve, which gives a TypeError, This could be a FiPy issue if it works with

Re: Taking Dot Product of faceGrad

2013-12-19 Thread Guyer, Jonathan E. Dr.
On Dec 19, 2013, at 12:59 PM, Daniel Wheeler daniel.wheel...@gmail.com wrote: I tried this and the residuals decrease quite a bit. Would you expect the residuals to approach zero at every time step? This seldom happens with non-linear equations. The absolute magnitude of the residuals isn't

Re: Taking Dot Product of faceGrad

2013-12-19 Thread Daniel Wheeler
On Thu, Dec 19, 2013 at 1:33 PM, Guyer, Jonathan E. Dr. jonathan.gu...@nist.gov wrote: Do 4th order diffusion terms even work with coupled? I thought they didn't. You're right, I forgot. So Jane, you'll need to split the equations as I suggested. Your current implementation combining fourth