[Kummi]

Hello Foamers,
This evaporation problem is a part of work in my main problem.

The problem is simple, a rectangular 1D channel is considered where right end is symmetry, left end is with the temperature of 373K, other faces are empty (1D problem). 10% of liquid (moisture) is considered inside the channel. The liquid is evaporated based on the condition as follows:

CONDITION: @ T= 373K ==> liquid (alpha) content =0
Based on it, mass and heat balances are calculated at the boiling plane (surface) as,
Mass balance (r) = rho * alpha * (dx/dt) [kg/m2.s]
Heat balance (-K dT/dx) = (r) * latent heat [W/m2]

Based on the condition, mass and heat balances are calculated.

I'm hereby attaching the solver (OF211). In this solver, under TEqn.H, I solved a if-else loop based on the above condition.

forAll(mesh.cells(), celli)
{
if (T[celli] >= Tsat.value())
{
alpha[celli] = scalar(0.0);
}
else
{
alpha[celli] = alphaTemp[celli];
// alpha[celli] = alpha[celli];
// alpha[celli] = alpha[nearestCellIndex];
}
}

Later, the solver is compiled and temperature varies inside the channel. However, the alpha (moisture content) is not solved properly with appropriate outcome. In my case, alpha has no expression (field operation), alpha is only solved based on the above condition as the loop (check Attachment plz). So, do we need an expression to solve alpha (like alpha depends on temperature expression) basically? any ideas ?

If someone have any ideas about it, kindly do share.

It will be helpful. Thank you