diff --git a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimpleChem/EEqn.H b/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimpleChem/EEqn.H
deleted file mode 100644
index ae52bae5..00000000
--- a/applications/solvers/cfdemSolverRhoPimple/cfdemSolverRhoPimpleChem/EEqn.H
+++ /dev/null
@@ -1,44 +0,0 @@
-// contributions to internal energy equation can be found in
-// Crowe et al.: "Multiphase flows with droplets and particles", CRC Press 1998
-
-{
-    // dim he = J / kg
-    volScalarField& he = thermo.he();
-    particleCloud.energyContributions(Qsource);
-
-
-    fvScalarMatrix EEqn
-    (
-        fvm::ddt(rhoeps, he) + fvm::div(phi, he)
-      + fvc::ddt(rhoeps, K) + fvc::div(phi, K)
-      + (
-            he.name() == "e"
-          ? fvc::div
-            (
-                fvc::absolute(phi/fvc::interpolate(rho), voidfraction*U),
-                p,
-                "div(phiv,p)"
-            )
-          : -dpdt
-        )
-     // net heat transfer from particles to fluid
-      - Qsource
-     // thermal conduction of the fluid with effective conductivity
-      - fvc::laplacian(rhoeps*particleCloud.thermCondM().thermDiff(),he) 
-          
-      // + particle-fluid energy transfer due to work
-      // + fluid energy dissipation due to shearing
-     ==
-        fvOptions(rho, he)
-    );
-
-    EEqn.relax();
-
-    fvOptions.constrain(EEqn);
-
-    EEqn.solve();
-
-    fvOptions.correct(he);
-
-    thermo.correct();
-}
diff --git a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/transportProperties b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/transportProperties
index 63cdca8c..47c0b1a4 100644
--- a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/transportProperties
+++ b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/transportProperties
@@ -19,6 +19,10 @@ transportModel  Newtonian;
 
 nu              nu [ 0 2 -1 0 0 0 0 ] 1.5e-04;
 
+kf		kf [ 1 1 -3 -1 0 0 0 ] 0;
+
+Cp		Cp [ 0 2 -2 -1 0 0 0 ] 1;
+
 CrossPowerLawCoeffs
 {
     nu0             nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
@@ -35,6 +39,4 @@ BirdCarreauCoeffs
     n               n [ 0 0 0 0 0 0 0 ] 1;
 }
 
-kf		kf [ 0 0 0 0 0 0 0 ] 0;
-Cp		Cp [ 0 0 0 0 0 0 0 ] 0;
 // ************************************************************************* //