Re: Adaptive mesh refinement in FiPy?
Hi Anders, No, sorry, there is no adaptive mesh refinement in FiPy. Cheers, Daniel On Mon, Nov 21, 2016 at 2:29 AM, Anders Ericsson wrote: > Hello, > > > I am wondering if there is any adaptive mesh refinement algorithm > implemented in FiPy? > > If so, is there some sort of manual on how to use it? > > > Thanks and best regards, > > Anders > > > ___ > fipy mailing list > fipy@nist.gov > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > -- Daniel Wheeler ___ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
RE: source term with
Thx for your reply. I was afraid it was not possible to force the solver. Do I understand you correctly that you propose to tackle this problem with something like: c_f[c_f < 0.] = 0. I tried that already, but that does not seem to help. Martin -Original Message- From: fipy-boun...@nist.gov [mailto:fipy-boun...@nist.gov] On Behalf Of Guyer, Jonathan E. Dr. (Fed) Sent: donderdag 17 november 2016 22:05 To: FIPY Subject: Re: source term with I don't believe there's anything you can do with the solver to enforce this constraint, but you can just manually "fix" the concentration in your sweep/solve loop. This is how the "double obstacle potential" (infinite energy barrier at zero concentration) is algorithmically implemented in my field. > On Nov 17, 2016, at 11:23 AM, Korevaar, Martin > wrote: > > However, due to some numerical inaccuracies, c_f sometimes becomes smaller > than zero and solutions fail. Is it possible to put a constraint on the > values of c_f so that the solver will only return solutions that comply to > this constraint? ___ 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 ]
