I added the .py file as attaching the notebook is somehow not supported On Monday, February 1, 2021 at 12:20:45 PM UTC+1 vitesh shah wrote:
> Hello, > > I am currently using Fipy, to solve a set of phase field equations. These > equations are based on the KWC model but with an additional extension made > to it, which enables it to account for stored energy of a material to be a > driving force. > > > > Currently, I have three equations which I attempt to solve using Fipy: > > 1. Phase equation > 2. Theta (misorientation) equation > 3. Equation for the stored energy. > > > > If I solve without the stored energy, the Fipy solver works fine and I get > expected behavior (in 1D at least). But when I also attempt to solve > equation for stored energy, I get a runtime warning: > > > “/nethome/v.shah/.local/lib/python3.8/site-packages/fipy/solvers/scipy/linearLUSolver.py:41: > > RuntimeWarning: invalid value encountered in double_scalars if > (numerix.sqrt(numerix.sum(errorVector**2)) / error0) <= self.tolerance:” > > So I suspect that solver is not handling these equations properly. And > this affects the solution that we get. > > > > We expect that for 1D case, the phase boundary should move when stored > energy is introduced (in normal phase field it should remain stationary). > But this behavior is not observed. > > Also when I compare the analytical rates and the rates that I observe from > Fipy and they have quite a different behavior. > > > > I am attaching a jupyter notebook here, so that the whole code (with > references to the papers) can be looked at. (link: if the attachment doesnt > work) > > > > It would be great, if you can give me some suggestions or point out some > mistakes that I might have made. > > > > With regards, > > Vitesh Shah > -- To unsubscribe from this group, send email to fipy+unsubscr...@list.nist.gov View this message at https://list.nist.gov/fipy --- To unsubscribe from this group and stop receiving emails from it, send an email to fipy+unsubscr...@list.nist.gov.
