Hello Jonathan and Daniel, Thanks for the useful suggestions. What do you mean by using Newton steps and correction variables? I don't think I have come across that before.
Best, Edwin On Fri, May 3, 2013 at 10:27 AM, Daniel Wheeler <[email protected]>wrote: > > > > 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 ] > >
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
