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 ]