[Mondal131211]

Hi Foamers,

I am trying to reproduce a simple result by buoyantBoussinesqSimpleFoam solver (OpenFoam 3.1.1 version). I have already set the simulation but not sure whether the set up is 100% correct or not. The description of the problem is attached here. This is basically a laminar case. I took geometry according to this paper. I think the difficulty is in choosing boundary conditions. I have chosen my boundary conditions according to the attached problem is as follows,

Velocity:
Hot wall – fixedValue, uniform (0 0 0)
Top and Bottom wall - fixedValue, uniform (0 0 0)
Open boundary – zeroGradient
Default faces-empty

Temperature:
Hot wall – fixedValue, uniform 333
Top and Bottom wall - zeroGradient
Open boundary – inletOutlet, Inletvalue - uniform 289, value- 289
Default faces-empty

P_rgh:
Hot wall – fixedFluxPressure, rho-rhok, uniform 0
Top and Bottom wall- fixedFluxPressure, rho-rhok, uniform 0
Open boundary – fixedValue, uniform 0
Default faces-empty

P:
Hot wall – Calculated (0)
Top and Bottom wall - Calculated (0)
Open boundary – Calculated (0)
Default faces-empty

I have got some result but it is not similar to that related problem.

Can anyone tell what is the mistake that I choose in my boundary conditions based on the OpenFOAM solvers?