Hello,

I am currently working through the Cahn-Hillard-example and my interested peaked!

As I understand it, the equations models phase separation, so it can be used to procedural generate configurations with two phases.

So would like to now if we could:

- take the solution at a certain time step,
- declare i.e. a solid and liquid phase based on the solution function (something like )
- solve a fluid-flow on the whole domain?

For the record I am not so much interested in a complete solution, because I like to challenge my self, but rather interested how the software specific concepts of FENICS could be used to do this. If nobody dislikes this idea I would also like to use this post to link informations, code and resources to this challenge, and hopefully at one point post the complete solution.