Thanks for your reply. It works well. I have one more question. I read the manual and find that for this example, phi is 0 for liquid, 1 for solid. So, does phi=0.5 represent the solid-liquid interface? What is the phase field value phi in the interface thickness?
Thanks a lot, william On Mon, Apr 19, 2010 at 2:46 PM, Jonathan Guyer <[email protected]> wrote: > > > On Apr 19, 2010, at 2:45 PM, Santan William wrote: > > > In examples.phase.anisotropy, I want to display phase field > vaiable phi at a certain node in each time step,e.g, I want to show the > value of phi at point x=y=1.How can I do this? > > I tried adding "... print "phase",phase" into the iteration as > following. It would show some values of phi, but I don't know how to show > the phi value of a certain point. > > There are a couple of ways, but if you don't know the exact cell centers, > this is probably best: > > print phi([[1.], [1.]]) > > > > > Another question is how can I print phase field variable for all > nodes at each time step into a file? For example, print phi of all nodes > into a file with a format x,y,phi for each time step. > > I appreciate your kindly help. > > See > http://www.ctcms.nist.gov/fipy/documentation/FAQ.html#how-do-i-export-data > >
