rhoCentralFoam: Updated the handling of thermal transport

Note that smoluchowskiJumpTFvPatchScalarField still uses the specified Pr from the
thermophysicalProperties whereas it should derive appropriate information from
the thermo package

applications/solvers/compressible/rhoCentralFoam/BCs/T/smoluchowskiJumpTFvPatchScalarField.C

@@ -185,12 +185,10 @@ void Foam::smoluchowskiJumpTFvPatchScalarField::updateCoeffs()
 
     dimensionedScalar Pr
     (
-        dimensionedScalar::lookupOrDefault
-        (
-            "Pr",
-            thermophysicalProperties,
-            1.0
-        )
+        "Pr",
+        dimless,
+        thermophysicalProperties.subDict("mixture").subDict("transport")
+        .lookup("Pr")
     );
 
     Field<scalar> C2

applications/solvers/compressible/rhoCentralFoam/readThermophysicalProperties.H

@@ -1,23 +0,0 @@
-Info<< "Reading thermophysicalProperties\n" << endl;
-
-// Pr defined as a separate constant to enable calculation of k, currently
-// inaccessible through thermo
-IOdictionary thermophysicalProperties
-(
-    IOobject
-    (
-        "thermophysicalProperties",
-        runTime.constant(),
-        mesh,
-        IOobject::MUST_READ_IF_MODIFIED,
-        IOobject::NO_WRITE
-    )
-);
-
-dimensionedScalar Pr
-(
-    "Pr",
-    dimless,
-    thermophysicalProperties.subDict("mixture").subDict("transport")
-   .lookup("Pr")
-);

applications/solvers/compressible/rhoCentralFoam/rhoCentralDyMFoam/rhoCentralDyMFoam.C

@@ -46,7 +46,6 @@ int main(int argc, char *argv[])
     #include "createTime.H"
     #include "createMesh.H"
     #include "createFields.H"
-    #include "readThermophysicalProperties.H"
     #include "readTimeControls.H"
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@@ -236,13 +235,10 @@ int main(int argc, char *argv[])
 
         if (!inviscid)
         {
-            volScalarField k("k", thermo.Cp()*muEff/Pr);
             solve
             (
                 fvm::ddt(rho, e) - fvc::ddt(rho, e)
               - fvm::laplacian(turbulence->alphaEff(), e)
-              + fvc::laplacian(turbulence->alpha(), e)
-              - fvc::laplacian(k, T)
             );
             thermo.correct();
             rhoE = rho*(e + 0.5*magSqr(U));

applications/solvers/compressible/rhoCentralFoam/rhoCentralFoam.C

@@ -45,7 +45,6 @@ int main(int argc, char *argv[])
     #include "createTime.H"
     #include "createMesh.H"
     #include "createFields.H"
-    #include "readThermophysicalProperties.H"
     #include "readTimeControls.H"
 
     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //