diff --git a/src/TurbulenceModels/compressible/RAS/buoyantKEpsilon/buoyantKEpsilon.C b/src/TurbulenceModels/compressible/RAS/buoyantKEpsilon/buoyantKEpsilon.C
index c2bf9607a44f0b853b0f7e333670592c70cd929d..167c62e5a678a68f3f3944428bc8fae2c9ad4329 100644
--- a/src/TurbulenceModels/compressible/RAS/buoyantKEpsilon/buoyantKEpsilon.C
+++ b/src/TurbulenceModels/compressible/RAS/buoyantKEpsilon/buoyantKEpsilon.C
@@ -107,7 +107,7 @@ buoyantKEpsilon<BasicTurbulenceModel>::Gcoef() const
         lookupObject<uniformDimensionedVectorField>("g");
 
     return
-        (this->Cmu_*this->k_/this->sigmak_)
+        this->alpha_*(this->Cmu_/this->sigmak_)*this->k_
        *(g & fvc::grad(this->rho_))/(this->epsilon_ + this->epsilonMin_);
 }