[bluemind]

Hello,

I'm currently trying to implement the saturation pressure into the solver twoLiquidMixingFoam, so far what I have written is :

{

forAll(psat,celli){
if (ELSATl[celli] > 200.)
{
psat[celli]=cfit*exp(c1 + c2/(Tl[celli] ) +c3*log(Tl[celli]/one ) +c4*pow(Tl[celli],c5));
}
else
{
psat[celli]=0.*pfit;
}


}

}

But when I compile i get the following error message :

error: cannot convert ‘Foam::dimensioned<double>’ to ‘double’ in assignment

I've chosen the unit of the constant c1..c5 and one, such that they give a dimensionless term inside the exponential and the log function. But cfit has a dimension of pressure. Somehow when I withdraw the 'celli' term from psat I don't have any problem, but since I need to apply the temperature condition for all control volume, I feel like I have to keep (celli) this term.

Does anybody know how to solve the problem ?