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? What else shall I try? Thank you very much!!!
On Wed, Jun 12, 2013 at 5:37 PM, Daniel Wheeler <[email protected]>wrote: > On Wed, Jun 12, 2013 at 6:37 AM, John Assael <[email protected]> wrote: > >> Thank you very much its working! >> > > Great! > > >> What I m trying to do is simulate results from COMSOL with FiPy, however >> although I have a small time step 2E-04, I cannot make a perfect fit. So, I >> suppose its the mesh... Is there anyway to optimize it without trial and >> error as there are many faces? >> > > Ideally, refining the meshes for both codes would bring the results into > agreement, but there are reasons why this might not necessarily be the > case. For example, FiPy will have error associated with non-orthogonality > in the mesh. COMSOL may have ways to handle this, which reduces its error > (maybe these corrections can be switched on and off in COMSOL). There could > also be different solution tolerances, solvers and other control > parameters. As a starting point, using exactly the same mesh for both FiPy > and COMSOL would be a good place to start if that's possible. > > -- > Daniel Wheeler > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > > -- ..:: ic3man.gr ::.. software | design | development _/_/_/ email: [email protected] _/_/_/ www : http://www.ic3man.gr ------------------------------------------------------------------ This e-mail and any attachments are confidential. You may not copy or disseminate any information contained in them to anyone other than the intended recipient. If you are not the intended recipient please contact the sender by reply e-mail and destroy all copies of the original message immediately. ------------------------------------------------------------------
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
