diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/calculateStress.H b/applications/solvers/stressAnalysis/solidDisplacementFoam/calculateStress.H
index 0703b841e2..dd0640aed4 100644
--- a/applications/solvers/stressAnalysis/solidDisplacementFoam/calculateStress.H
+++ b/applications/solvers/stressAnalysis/solidDisplacementFoam/calculateStress.H
@@ -15,8 +15,7 @@ if (runTime.writeTime())
 
     if (thermo.thermalStress())
     {
-        const volScalarField& T = Tptr();
-        sigma = sigma - I*(rho*threeKalpha*T);
+        sigma = sigma - I*(rho*threeKalpha*thermo.T());
     }
 
     volScalarField sigmaEq
diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/createFields.H b/applications/solvers/stressAnalysis/solidDisplacementFoam/createFields.H
index 40b16ad481..5b800e1780 100644
--- a/applications/solvers/stressAnalysis/solidDisplacementFoam/createFields.H
+++ b/applications/solvers/stressAnalysis/solidDisplacementFoam/createFields.H
@@ -14,30 +14,6 @@ volVectorField D
     mesh
 );
 
-
-autoPtr<volScalarField> Tptr(nullptr);
-
-if (thermo.thermalStress())
-{
-    Info<< "Reading field T\n" << endl;
-    Tptr.reset
-    (
-        new volScalarField
-        (
-            IOobject
-            (
-                "T",
-                runTime.timeName(),
-                mesh,
-                IOobject::MUST_READ,
-                IOobject::AUTO_WRITE
-            ),
-            mesh
-        )
-    );
-}
-
-
 Info<< "Calculating stress field sigmaD\n" << endl;
 volSymmTensorField sigmaD
 (
diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/solidDisplacementFoam.C b/applications/solvers/stressAnalysis/solidDisplacementFoam/solidDisplacementFoam.C
index 31df2258f1..fec2a0941f 100644
--- a/applications/solvers/stressAnalysis/solidDisplacementFoam/solidDisplacementFoam.C
+++ b/applications/solvers/stressAnalysis/solidDisplacementFoam/solidDisplacementFoam.C
@@ -69,7 +69,7 @@ int main(int argc, char *argv[])
         {
             if (thermo.thermalStress())
             {
-                volScalarField& T = Tptr();
+                volScalarField& T = thermo.T();
                 fvScalarMatrix TEqn
                 (
                     fvm::ddt(rho, Cp, T)
@@ -96,8 +96,7 @@ int main(int argc, char *argv[])
 
                 if (thermo.thermalStress())
                 {
-                    const volScalarField& T = Tptr();
-                    DEqn += fvc::grad(threeKalpha*T);
+                    DEqn += fvc::grad(threeKalpha*thermo.T());
                 }
 
                 fvOptions.constrain(DEqn);