[Kummi]

Hello everyone,

I want to set the liquid evaporation at the interface. The interface is at the fixed saturation temperature of T = 100deg. My conditions at the interface (boiling plane) are as follows:

T= 100 ==> alpha [moisture content) =0]
Based on it, mass and heat balances are calculated at interface as,
Mass balance (r) = rho * alpha * (dx/dt) [kg/m2.s]
Heat balance (-K dT/dx) = (r) * latent heat [W/m2]

There found no variation in Temperature and alpha after resolving the below code at TEqn.H

volScalarField alphaPrime = mag(gradAlphaX);
// set liquid and interfacial temperature to Tsat
forAll(mesh.cells(), celli)
{
if (alpha[celli] = scalar(0.1)) //10% of moisture
{
T[celli] = scalar (373.0); //Tsat = 100deg
alphaPrime[celli]= scalar(0); //moisture content = 0
}
}
volScalarField VX = V.component(vector::X); //VX =interface_VELOCITY
volScalarField mDot = -(rho*alpha*VX); //mass balance
{
fvScalarMatrix TEqn
(
fvm::ddt(rho*cp, T) // fvm::ddt(cp, T)
- mDot * cp * gradTx //source term
- fvm::laplacian(K, T) // - fvm::laplacian(lambda/rho, T)
);
TEqn.relax();
TEqn.solve();
}

Some assistance will be highly helpful. Thank you