You might want to multiply by the cell volumes, "mesh.cellVolumes".
On Thu, Jan 4, 2018 at 10:22 AM, Anders Ericsson
<anders.erics...@solid.lth.se> wrote:
> Hi,
>
>
> I wonder if there is a simple way to compute the total free energy of a
> system (Phase-field modeling) in FiPy?
>
>
> That is e.g.:
>
>
> \begin{equation}
> F = \int_V f(\phi, c, T) + \frac{\epsilon_{\phi}^2}{2}\nabla \phi^2 dV
> \end{equation}
>
> I figured that it would be something in accordance with:
>
> def freeEnergyVolume(phi_, c_, T, epsSq):
> return (0.5 * epsSq * (phi_.grad.mag)**2 +
> f(phi_,c_,T)).cellVolumeAverage
>
> But cellVolumeAverage wouldn't give me the full energy of the domain if I'm
> not mistaken?
--
Daniel Wheeler
_______________________________________________
fipy mailing list
fipy@nist.gov
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
