I apologize. I apparently never hit "send" in reply to your email to me on Friday. I've replied to you now and will look at your attached code when I get the chance.
I did not write the code in that paper, but advised the lead author in formulating the code and debugging it. If I've got a copy, I can't find it. I've sent you a current address for the lead author; he may be able to provide you with a script. On Apr 29, 2014, at 4:17 PM, "Johnson, Luke A" <lukejohn...@neo.tamu.edu> wrote: > Daniel, > The code I was referring to is the code used in the paper titled "Diffusion > under temperature gradient: A phase-field model study" by Mohanty, Guyer and > Sohn. I found an archived conversation from this mailing list > "http://article.gmane.org/gmane.comp.python.fipy/3064/match=thermal+gradient" > where Jonathan Guyer mentioned digging up the FiPy code used for the paper > and was wondering if that would still be possible. > > I have attached my code as well, in response to your other email. The code is > quite temperamental as I have combined elements from quite a few sources but > it should converge to a solution as it is written right now, albeit an > unphysical one. > > Thank you for your time and interest, > Luke Johnson > > > From: "Daniel Wheeler" <daniel.wheel...@gmail.com> > To: "Multiple recipients of list" <email@example.com> > Sent: Tuesday, April 29, 2014 2:52:29 PM > Subject: Re: Question about Diffusion Under Temperature Gradient > > On Mon, Apr 28, 2014 at 7:00 PM, Johnson, Luke A > <lukejohn...@neo.tamu.edu> wrote: > Hello, > I am trying to create a relatively simple thermomigration model in FiPy for > a final project in my Kinetics class and have been using a paper titled > Diffusion Under Temperature Gradient as a reference. I can get my system to > converge and segregate if I assume constant values for the Chemical Mobility > and Mobility of Thermotransport but am having convergence problems when the > mobilities become functions of temperature and composition. I think my code > is right but the problem is that I am using unrealistic values for the > parameters not explicitly listed in the paper. Could you provide the > thermodynamic system parameters you used and maybe a workflow for the code > so I can check with my work? > > Unfortunately, I am not sure what you're referring to. Which code / > > example are you asking about? > > -- > Daniel Wheeler > _______________________________________________ > fipy mailing list > firstname.lastname@example.org > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > <finalproject_dimension.py>_______________________________________________ > fipy mailing list > email@example.com > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] _______________________________________________ fipy mailing list firstname.lastname@example.org http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]