On Wed, Jun 12, 2013 at 1:24 PM, John Assael <[email protected]> wrote:
> The COMSOL resutls make a 99% fit to the experimental results so they > should be right. Now, when I am trying to fit them in FiPy and although I > have generated the same mesh i dont get the same results (so I am not quite > sure about the cellsizes as they seem to make quite a difference). > Moreover, I have also tried all the available fipy.solvers with no > success, with various tolerances and iterations. > Is there a way to check the mesh quality? > You said that you are using the same mesh for both FiPy and COMSOL and the mesh is of sufficient quality for COMSOL to give what you believe is a reasonable answer. I'm not sure that the problem is with the mesh in this case since they are the same. What else shall I try? > Is there a way to switch off the correction terms in COMSOL and see what happens? Another possibility is to solve a similar problem with a Cartesian grid and then compare FiPy and COMSOL. My guess is that COMSOL does a better job with non-orthogonality than FiPy. We know that non-orthogonal correction terms is something that needs to be addressed. -- Daniel Wheeler
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
