[Thamali]

Hi foamers,

I am trying to run icoIbFoam including energy equation to see flow around a 1mm,spherical particle with 350K constant temperature.I refined the mesh and it is attached.
Since the results were showing very high Courant numbers and very high velocitites,etc with very high mesh refinement(attached one is after reducing mesh refinement and without rapid mesh rerinement),I changed the solver to vary deltaT to adjust to make Courant number less than 1.But, still the velocities go around 100s.My inlet velocity is (0.1569 0 0).
velocity boundary condition;

boundaryField
{
ibSphere
{
type immersedBoundary;
refValue uniform (0 0 0);
refGradient uniform (0 0 0);
fixesValue yes;

setDeadCellValue yes;
deadCellValue (0 0 0);

value uniform (0 0 0);
}

in
{
type fixedValue;
value uniform (0.15690 0 0);
}

out
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0.15690 0 0);
}

top
{
type fixedValue;
value uniform (0 0 0);
}

bottom
{
type fixedValue;
value uniform (0 0 0);
}

back
{
type fixedValue;
value uniform (0 0 0);
}

front
{
type fixedValue;
value uniform (0 0 0);
}
}

I am attaching the p and U profiles after few seconds.
My domain size is 40mm*20mm*20mm.
The results seems making an error due to mesh refinement problem.I used the method like in sphereInChannel for refinement.

Can someone help me to figure out my proble?

Thanks in advance.
Thamali