Project 3- Adding a temperature dictionary to the incompressible icoFoam case to solve the energy equation

relpce DT as K

dimensionedScalar DT
(
    
    
    transportProperties.lookup("DT")
);

Info<< "Reading field T\n" << endl;
volScalarField T
(
    IOobject
    (
        "T",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

 fvScalarMatrix TEqn
        (
            fvm::ddt(T)
          + fvm::div(phi, T)
          -fvm::laplacian(DT, T)
        );
      TEqn.solve();
        }

Please find the attched latex report