On Thu, May 2, 2013 at 3:18 PM, Jonathan Guyer <[email protected]> wrote:
> > On May 2, 2013, at 12:19 PM, Edwin Sze Lun Khoo <[email protected]> wrote: > > > Assigning the coefficient of the convection term as -alpha/C**2 seems > incorrect because the coefficient of the convection term should be that of > whatever multiplies the dependent variable inside the divergence operator; > -alpha/C**2 here, however, is outside and not inside the d/dx operator. > > You're right, of course. I don't see a way to make this term implicit. > Wheeler might know a way. > What about using a Newton steps for this? Then the transient and diffusion stay the same but the convection becomes - (2 C \delta C)_x where \delta C is the correction variable. -- Daniel Wheeler
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
