Temperature BC inside a phase
 

Folks,
Do you know how to impose a constant temperature inside one of the phases (in a VOF model) throughout the simulation? I tried the following but it did not work.

Thread *lt = pt[liquidphaseindex];

begin_c_loop_all(c,t)
{
if (C_VOF(c,lt) == 0.0)
{
C_T(c,t) = 329;
}
}
end_c_loop_all(c,t)

The temperature inside that phase is decreasing to a much lower value.