[gomeslorenzo]

Good morning foamers,

I'm using chtmultiregionfoam to simulate a heatsink under natural convection. My problem is that the solution never converges. I ploted the temperature x time function of the tip of a fin to see what is happening and the result was a straight line. I think there is a problem with the BC but I can't find where. Could anyone help me?

Here is the BC file for the heatsink ( named here heater ):

dictionaryReplacement
{
boundary
{
minZ
{
type wall;
}

}

T
{
internalField uniform 300;

boundaryField
{
".*"
{
type zeroGradient;
value uniform 300;
}

"heater_to_.*"
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
Tnbr T;
kappa solidThermo;
kappaName none;
value uniform 300;
}

minZ
{
type compressible::turbulentHeatFluxTemperature;
heatSource flux; // power [W]; flux [W/m2]
q uniform 200; // heat power or flux
kappaName none;
kappa solidThermo; // calculate K by alphaEff*thermo.Cp
value uniform 300; // initial temperature value
}

}
}
}

Thanks for any help.