buoyantBoussinesqPimpleFoam solution diverge at reference pressure cell
 

Hello,
I am using userbuoyantBoussinesqPimpleFoam in OF 2.3.1. to solve natural convection on heated cylinder.

BCS used are
Inlet:
T fixedValue
U zeroGradient
P fixedFluxPressure
P_rgh calculated
Cylinder
T fixedValue
U fixedValue
P fixedFluxPressure
P_rgh calculated
side Walls:
T zerogradient
U slip
P fixedFluxPressure
P_rgh calculated
outlet:
T zeroGradient
U zerogradient
P fixedFluxPressure
P_rgh calculated

When running the simulation with these BSCs, solution get diverge at reference pressure cell.

I also tried by changing outlet P_rgh BSC to fixedValue, then solution at outlet boundary is weird.

Hi,
did you manage to find a solution? I am facing the same problem for days and I fear that I am going crazy.