diff --git a/src/TurbulenceModels/phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C b/src/TurbulenceModels/phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
index 83edf7de1574d1416c9a51cf0454e71414e3667d..91c401ee54ea15340c4ed00bddd8c7f45c9810e8 100644
--- a/src/TurbulenceModels/phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
+++ b/src/TurbulenceModels/phaseIncompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
@@ -194,7 +194,6 @@ tmp<volScalarField> LaheyKEpsilon<BasicTurbulenceModel>::bubbleG() const
     tmp<volScalarField> bubbleG
     (
         Cp_
-       *liquid*liquid.rho()
        *(
             pow3(magUr)
           + pow(fluid.drag(gas).CdRe()*liquid.nu()/gas.d(), 4.0/3.0)