how to use INC_DENSITY_MODEL= BOUSSINESQ to solve the laminar convection in a cavity probelm.

[ztdepztdep]

INC_DENSITY_MODEL= BOUSSINESQ
INC_ENERGY_EQUATION= YES
FLUID_MODEL= CONSTANT_DENSITY
BODY_FORCE= NO
THERMAL_EXPANSION_COEFF= 0.003333333333
INC_TEMPERATURE_INIT= 300.0
INC_TEMPERATURE_REF= 300.0

I have set the above value, but i still cann't got a converged solution.
the rms(H) still ossicilating about -0.66.

[bigfooted]

I see that we do not have a testcase for it in the testcases repository, but we do have a buoyancy problem defined using a user defined source term in py_wrapper/custom_source_buoyancy/
you can check the numerical settings there and see if it is different from yours.