I am looking to set up a simulation of spherical particles being melted by a focussed heatsource. The free-surface code seems ideal, but it lacks a solid state and the possibility to transition between solid and liquid.
Is the addition of an extra state to the free-surface code a difficult job? I have scientific coding experience, but not much experience of directly writing in any C variants. I have several months in which I can develop the underlying code, before I need to start attempting to get physically accurate results out.
Thank you in advance for any musing