How to add the heat boudary in the free surface?
 

I want to solve a free surface problem with HOL, How can I track the free surface in order to add the heat transfer boudary which is related to the moving of free surface? I try to use VFOL as a condition and program code in Groud.for, but the result shows some error. Could you give some suggestion about it? Think for you kind help.