chtMultiRegionFoam: Removed Qdot from top-level solver

Qdot is only relevant for reacting cases, and even then it is only
written out for post-processing purposes. It has been removed from
chtMultiRegionFoam as this solvers' typical usage is for non-reacting
simulations.

If it is desired to obtain Qdot from a chtMultiRegionFoam simulation,
then a better way would be to implement a function object to look up the
reaction model and write it out.

applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/EEqn.H

@@ -19,7 +19,7 @@
      ==
         rho*(U&g)
       + rad.Sh(thermo, he)
-      + Qdot
+      + reaction.Qdot()
       + fvOptions(rho, he)
     );
 

applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/YEqn.H

@@ -16,13 +16,21 @@ if (Y.size())
 
 {
     reaction.correct();
-    Qdot = reaction.Qdot();
-    volScalarField Yt
-    (
-        IOobject("Yt", runTime.timeName(), mesh),
-        mesh,
-        dimensionedScalar("Yt", dimless, 0)
-    );
+
+    tmp<volScalarField> Yt(nullptr);
+
+    if (Y.size())
+    {
+        Yt = tmp<volScalarField>
+        (
+            new volScalarField
+            (
+                IOobject("Yt", runTime.timeName(), mesh),
+                mesh,
+                dimensionedScalar("Yt", dimless, 0)
+            )
+        );
+    }
 
     forAll(Y, i)
     {
@@ -49,7 +57,7 @@ if (Y.size())
             fvOptions.correct(Yi);
 
             Yi.max(0.0);
-            Yt += Yi;
+            Yt.ref() += Yi;
         }
     }
 

applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/createFluidFields.H

@@ -15,7 +15,6 @@ PtrList<volScalarField> KFluid(fluidRegions.size());
 PtrList<volScalarField> dpdtFluid(fluidRegions.size());
 PtrList<multivariateSurfaceInterpolationScheme<scalar>::fieldTable>
     fieldsFluid(fluidRegions.size());
-PtrList<volScalarField> QdotFluid(fluidRegions.size());
 
 List<scalar> initialMassFluid(fluidRegions.size());
 
@@ -245,25 +244,6 @@ forAll(fluidRegions, i)
     }
     fieldsFluid[i].add(thermoFluid[i].he());
 
-    Info<< "    Adding to QdotFluid\n" << endl;
-    QdotFluid.set
-    (
-        i,
-        new volScalarField
-        (
-            IOobject
-            (
-                "Qdot",
-                runTime.timeName(),
-                fluidRegions[i],
-                IOobject::READ_IF_PRESENT,
-                IOobject::AUTO_WRITE
-            ),
-            fluidRegions[i],
-            dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
-        )
-    );
-
     Info<< "    Adding MRF\n" << endl;
     MRFfluid.set
     (

applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/setRegionFluidFields.H

@@ -43,8 +43,6 @@
     multivariateSurfaceInterpolationScheme<scalar>::fieldTable& fields =
         fieldsFluid[i];
 
-    volScalarField& Qdot = QdotFluid[i];
-
     radiation::radiationModel& rad = radiation[i];
 
     IOMRFZoneList& MRF = MRFfluid[i];