Hi everybody, i'm trying to
 

Hi everybody,

i'm trying to simulate the flow in a vertical channel due to natural convection, with each plate at a constant temperature. I use the k epsilon turbulence model.
I have choosen a cyclic boundary condition in the streamwise and spamwise directions.
For the two heated plate, i have pd-> fixedFluxBuoyantPressure, p->calculated, no slip condition, fixedValue for T and zeroGradient for k and epsilon.
I use second order discretisation schemes.

Running buoyantSimpleFoam, i have a bad convergence with a fine mesh and the computation with a coarser much crashes after some iterations.

The number of boundary conditions being quite limited, i don't see where is the problem.

Could someone help me?

Thank you

Eric