diff --git a/applications/solvers/rcfdemSolverRhoSteadyPimpleChem/createFields.H b/applications/solvers/rcfdemSolverRhoSteadyPimpleChem/createFields.H
index e7640d48..ac35d627 100644
--- a/applications/solvers/rcfdemSolverRhoSteadyPimpleChem/createFields.H
+++ b/applications/solvers/rcfdemSolverRhoSteadyPimpleChem/createFields.H
@@ -152,6 +152,22 @@ Info<< "Reading thermophysical properties\n" << endl;
         dimensionedScalar("zero", dimensionSet(1,-1,-3,-1,0,0,0), 0.0)
     );
 
+    Info<< "\nCreating fluid thermal conduction field\n" << endl;
+    volScalarField QFluidCond
+    (
+        IOobject
+        (
+            "QFluidCond",
+            runTime.timeName(),
+            mesh,
+            IOobject::NO_READ,
+            IOobject::NO_WRITE
+        ),
+        mesh,
+        dimensionedScalar("zero", dimensionSet(1,-1,-3,0,0,0,0), 0.0)
+    );
+
+
     Info<< "\nCreating thermal conductivity field\n" << endl;
     volScalarField thCond
     (
diff --git a/applications/solvers/rcfdemSolverRhoSteadyPimpleChem/monitorMass.H b/applications/solvers/rcfdemSolverRhoSteadyPimpleChem/monitorMass.H
index 5940864f..40b6d36f 100644
--- a/applications/solvers/rcfdemSolverRhoSteadyPimpleChem/monitorMass.H
+++ b/applications/solvers/rcfdemSolverRhoSteadyPimpleChem/monitorMass.H
@@ -4,4 +4,6 @@
     counter++;
     Info << "\ncurrent gas mass = " << m << "\n" << endl;
     Info << "\ncurrent added gas mass = " << m-m0 << "\n" << endl;
+
+    QFluidCond = fvc::laplacian(voidfractionRec*thCond,T);
 }