Using rhok for a modified kepsilon model. How to access to the value?
 

Hi all,

Based on buoyantBoussinesqSimpleFoam, I want to used rhok variable to calculate a extra term in the k equation of kEpsilon model. The problem is that what is more or less easy in 2.1.1, in 1612+ is not so easy for me. I can not access to the rhok variable.

OF211 rhok_(U.db().lookupObject<volScalarField>("rhok")) ,
OF1612+ ??

The target is to obtain a buoyantkEpsilon model for incompressible based on Boussinesq approximation.

Thank you very much

replace U with the variable you use in ke model constructor