Used Heat Flux instead of tempearture
 

Hi
I am working on multi phase change solver called, reactingTwoPhaseEulerFoam.I have boiling and condensation at boundary. instead of constant temperature at wall, I want to apply heat flux and openFoam automatically calculate wall temperature. The library that I used is called "fixedMultiPhaseHeatFlux, for example when I have boiling, I set this boundary for T.water at wall and for another phase T.steam I applied, "copyFixedValue"(see end of message). My problem is that, when I set heat flux of 2000 W/m2 for wall, when my running is complete, the value of heat flux that I gather from all control volume at wall is not exactly the value that I set as input, it has error of 4-15 %. I want to know anyone has ever worked with library before, what is your experience about this?
thanks a lot for your attention

************************************************
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object T.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 286.17;

boundaryField
{
evaporatorLiq
{
type fixedMultiPhaseHeatFlux;
relax 0.5;
q uniform 2000;
phase "liquid";
value uniform 289.6;

}
***********************************************
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object T.steam;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 286.17;

boundaryField
{
evaporatorLiq
{
type copiedFixedValue;
sourceFieldName T.water;
value uniform 289.6;


}

Hey! Have you solve the problem? I met the same question.