ENH: radiation - temperature dependent properties

ENH: radiationModel - expose T()

src/thermophysicalModels/radiation/derivedFvPatchFields/MarshakRadiation/MarshakRadiationFvPatchScalarField.C

@@ -139,7 +139,7 @@ void Foam::radiation::MarshakRadiationFvPatchScalarField::updateCoeffs()
 
     const tmp<scalarField> temissivity
     (
-        boundaryRadiation.emissivity(patch().index())
+        boundaryRadiation.emissivity(patch().index(), 0, nullptr, &Tp)
     );
 
     const scalarField& emissivity = temissivity();

src/thermophysicalModels/radiation/derivedFvPatchFields/greyDiffusiveRadiation/greyDiffusiveRadiationMixedFvPatchScalarField.C

@@ -178,14 +178,14 @@ updateCoeffs()
 
     const tmp<scalarField> temissivity
     (
-        boundaryRadiation.emissivity(patch().index())
+        boundaryRadiation.emissivity(patch().index(), 0, nullptr, &Tp)
     );
 
     const scalarField& emissivity = temissivity();
 
     const tmp<scalarField> ttransmissivity
     (
-        boundaryRadiation.transmissivity(patch().index())
+        boundaryRadiation.transmissivity(patch().index(), 0, nullptr, &Tp)
     );
 
     const scalarField& transmissivity = ttransmissivity();

src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C

@@ -167,15 +167,17 @@ updateCoeffs()
         boundaryRadiationProperties::New(internalField().mesh());
 
 
+    const auto& Tp = radiation.T().boundaryField()[patchi];
+
     const tmp<scalarField> temissivity
     (
-        boundaryRadiation.emissivity(patch().index(), lambdaId)
+        boundaryRadiation.emissivity(patchi, lambdaId, nullptr, &Tp)
     );
     const scalarField& emissivity = temissivity();
 
     const tmp<scalarField> ttransmissivity
     (
-        boundaryRadiation.transmissivity(patch().index(), lambdaId)
+        boundaryRadiation.transmissivity(patchi, lambdaId, nullptr, &Tp)
     );
     const scalarField& transmissivity = ttransmissivity();
 
@@ -217,7 +219,7 @@ updateCoeffs()
              const volScalarField& qSec =
                 this->db().lookupObject<volScalarField>(qSecName);
 
-            Ir += qSec.boundaryField()[patch().index()];
+            Ir += qSec.boundaryField()[patchi];
         }
     }
 

src/thermophysicalModels/radiation/radiationModels/radiationModel/radiationModel.H

@@ -240,6 +240,9 @@ public:
                 return radiation_;
             }
 
+            //- Return access to the temperature field
+            const volScalarField& T() const noexcept { return T_; }
+
             //- Source term component (for power of T^4)
             virtual tmp<volScalarField> Rp() const = 0;
 

src/thermophysicalModels/radiation/radiationModels/viewFactor/viewFactor.C

@@ -738,14 +738,14 @@ void Foam::radiation::viewFactor::calculate()
 
             fvPatchScalarField& qrPatch = qrBf[patchID];
 
-            greyDiffusiveViewFactorFixedValueFvPatchScalarField& qrp =
+            auto& qrp =
                 refCast
                 <
                     greyDiffusiveViewFactorFixedValueFvPatchScalarField
                 >(qrPatch);
 
             const tmp<scalarField> teb =
-                boundaryRadiation.emissivity(patchID, bandI);
+                boundaryRadiation.emissivity(patchID, bandI, nullptr, &Tp);
 
             const scalarField& eb = teb();