ThermophysicalTransportModels: Corrected the documentation of the units of the q() function

heSolidThermo: Changed the q() function to return an intensive flux consistent
with ThermophysicalTransportModels.

src/ThermophysicalTransportModels/laminar/Fourier/Fourier.H

@@ -112,7 +112,7 @@ public:
             const label patchi
         ) const;
 
-        //- Return the heat flux [W]
+        //- Return the heat flux [W/m^2]
         virtual tmp<surfaceScalarField> q() const;
 
         //- Return the source term for the energy equation

src/ThermophysicalTransportModels/laminar/unityLewisFourier/unityLewisFourier.H

@@ -118,7 +118,7 @@ public:
             return this->alphaEff(patchi);
         }
 
-        //- Return the heat flux [W]
+        //- Return the heat flux [W/m^2]
         virtual tmp<surfaceScalarField> q() const;
 
         //- Return the source term for the energy equation

src/ThermophysicalTransportModels/thermophysicalTransportModel/thermophysicalTransportModel.H

@@ -127,7 +127,7 @@ public:
             const label patchi
         ) const = 0;
 
-        //- Return the heat flux [W]
+        //- Return the heat flux [W/m^2]
         virtual tmp<surfaceScalarField> q() const = 0;
 
         //- Return the source term for the energy equation

src/ThermophysicalTransportModels/turbulence/eddyDiffusivity/eddyDiffusivity.H

@@ -181,7 +181,7 @@ public:
             const label patchi
         ) const;
 
-        //- Return the heat flux [W]
+        //- Return the heat flux [W/m^2]
         virtual tmp<surfaceScalarField> q() const;
 
         //- Return the source term for the energy equation

src/ThermophysicalTransportModels/turbulence/nonUnityLewisEddyDiffusivity/nonUnityLewisEddyDiffusivity.H

@@ -143,7 +143,7 @@ public:
             );
         }
 
-        //- Return the heat flux [W]
+        //- Return the heat flux [W/m^2]
         virtual tmp<surfaceScalarField> q() const;
 
         //- Return the source term for the energy equation

src/ThermophysicalTransportModels/turbulence/unityLewisEddyDiffusivity/unityLewisEddyDiffusivity.H

@@ -199,7 +199,7 @@ public:
             return alphaEff(patchi);
         }
 
-        //- Return the heat flux [W]
+        //- Return the heat flux [W/m^2]
         virtual tmp<surfaceScalarField> q() const;
 
         //- Return the source term for the energy equation

src/functionObjects/field/wallHeatFlux/wallHeatFlux.C

@@ -79,13 +79,12 @@ Foam::functionObjects::wallHeatFlux::calcWallHeatFlux
         twallHeatFlux.ref().boundaryFieldRef();
 
     const surfaceScalarField::Boundary& qBf = q.boundaryField();
-    const surfaceScalarField::Boundary& magSf = mesh_.magSf().boundaryField();
 
     forAllConstIter(labelHashSet, patchSet_, iter)
     {
         const label patchi = iter.key();
 
-        wallHeatFluxBf[patchi] = -qBf[patchi]/magSf[patchi];
+        wallHeatFluxBf[patchi] = -qBf[patchi];
     }
 
     if (foundObject<volScalarField>("qr"))

src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.C

@@ -340,7 +340,7 @@ Foam::heSolidThermo<BasicSolidThermo, MixtureType>::q() const
             isotropic()
           ? fvm::laplacian(this->kappa(), this->T_)().flux()
           : fvm::laplacian(KappaLocal(), this->T_)().flux()
-        );
+        )/mesh.magSf();
 }
 
 

src/thermophysicalModels/solidThermo/solidThermo/heSolidThermo.H

@@ -115,7 +115,7 @@ public:
         //  in the local coordinate system [W/m/K]
         virtual tmp<symmTensorField> KappaLocal(const label patchi) const;
 
-        //- Return the heat flux [W]
+        //- Return the heat flux [W/m^2]
         virtual tmp<surfaceScalarField> q() const;
 
         //- Return the source term for the energy equation