Dear Raymond, Jonathan and Dan, Thanks a lot for all your replies.
Indeed this works ! Setting the diffusion coefficient to all interior faces helps a lot. Thanks for your help once again. Krishna -----Original Message----- From: fipy-boun...@nist.gov [mailto:fipy-boun...@nist.gov] On Behalf Of Daniel Wheeler Sent: Wednesday, July 13, 2016 3:17 PM To: Multiple recipients of list <fipy@nist.gov> Subject: Re: Possibility of setting NaN constraints for CellVariable ? On Wed, Jul 13, 2016 at 10:06 AM, Guyer, Jonathan E. Dr. (Fed) <jonathan.gu...@nist.gov> wrote: > You can't do this. NaNs immediately propagate and will break your whole > solution. > > If I had to do this, I would decouple this region from the other equations by > zeroing out the coefficients of those other equations in the appropriate > region. Seconded, remember that the coefficients at the faces need to be zero to not have propagation. This will require masking the faces based on position. -- Daniel Wheeler _______________________________________________ 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 ]