On Mon, Jul 7, 2014 at 5:55 AM, Ivas Toni <toni.i...@mat.ethz.ch> wrote: > Dear All, > > I am testing my C++ code for Kobayashi-Warren-Carter model by Fipy. My > question is concerning sourceCoeff in mesh20x20.py example from Fipy > distribution: > sourceCoeff = (diffusionCoeff * thetaGradDiff).divergence > This term is explained I think in J.A. Warren et al. / Acta Materialia 51 > (2003) 6035–6058 Eqs.(A28-A31). > The subtraction operator is defined as Subt(\theta_0-\theta_1): > \theta_1 - \theta_0 + 2\pi if \theta_1 - \theta_0< -\pi > \theta_1 - \theta_0 if -\pi < \theta_1 - \theta_0< \pi > > \theta_1 - \theta_0 - 2\pi if \pi <= \theta_1 - \theta_0
Hi Toni, Thanks for the question. Unfortunately, I need to sit down and work through this again. I will hopefully get back to you in a week or two when I have some time. Regarding comparison with your C++ model, you may want to compare the output numbers directly for each field with "sourceCoeff" on and off and see. Cheers, Daniel -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
