Hello all, This question is related to the one posed by Thomas Barkar last week.
We are trying to solve the Cahn-Hilliard problem using thermodynamics and mobilities obtained from external routines. To illustrate the problem 3 sample codes are included. In code1.py the expression for the thermodynamics is entered in the equation definition. In code2.py (included to illustrate the difference) the expression for thermodynamics is moved to a function. And in code3.py the evaluation of the thermodynamic properties in each point has been split up in a loop. All these samples produce different results (where code1.py is the correct). In order to get the coupling with the external thermodynamics package we need to use an approach that is shown in code3.py. How can we define the problem in such a way so we can obtain the same results in code1.py and code3.py? Thanks! Joakim
code1.py
Description: code1.py
code2.py
Description: code2.py
code3.py
Description: code3.py
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
