From 99aa00e33ca92fb060189ca2432111861c078266 Mon Sep 17 00:00:00 2001
From: sergio <sergio>
Date: Thu, 4 Oct 2012 10:55:49 +0100
Subject: [PATCH] ENH: Adding greyMeanSolidAbsorptionEmission for radiative
 properties of the solid

---
 .../radiationModels/Make/files                |   2 +
 .../radiationCoupledBase.C                    |  16 +-
 .../radiationModel/opaqueSolid/opaqueSolid.C  | 130 +++++++++
 .../radiationModel/opaqueSolid/opaqueSolid.H  | 109 ++++++++
 .../radiationModel/radiationModel.C           |   9 +-
 .../radiationModel/radiationModel.H           |   5 +-
 .../greyMeanSolidAbsorptionEmission.C         | 256 ++++++++++++++++++
 .../greyMeanSolidAbsorptionEmission.H         | 158 +++++++++++
 .../mixtures/SpecieMixture/SpecieMixture.C    |  12 +
 .../mixtures/SpecieMixture/SpecieMixture.H    |   8 +
 .../basicMultiComponentMixture.H              |   8 +
 .../constant/panelRegion/radiationProperties  |  70 +++++
 12 files changed, 776 insertions(+), 7 deletions(-)
 create mode 100644 src/thermophysicalModels/radiationModels/radiationModel/opaqueSolid/opaqueSolid.C
 create mode 100644 src/thermophysicalModels/radiationModels/radiationModel/opaqueSolid/opaqueSolid.H
 create mode 100644 src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.C
 create mode 100644 src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.H
 create mode 100644 tutorials/combustion/fireFoam/les/oppositeBurningPanels/constant/panelRegion/radiationProperties

diff --git a/src/thermophysicalModels/radiationModels/Make/files b/src/thermophysicalModels/radiationModels/Make/files
index 9822712105..9e49da1fde 100644
--- a/src/thermophysicalModels/radiationModels/Make/files
+++ b/src/thermophysicalModels/radiationModels/Make/files
@@ -8,6 +8,7 @@ radiationModel/fvDOM/radiativeIntensityRay/radiativeIntensityRay.C
 radiationModel/fvDOM/blackBodyEmission/blackBodyEmission.C
 radiationModel/fvDOM/absorptionCoeffs/absorptionCoeffs.C
 radiationModel/viewFactor/viewFactor.C
+radiationModel/opaqueSolid/opaqueSolid.C
 
 /* Scatter model */
 submodels/scatterModel/scatterModel/scatterModel.C
@@ -23,6 +24,7 @@ submodels/absorptionEmissionModel/constantAbsorptionEmission/constantAbsorptionE
 submodels/absorptionEmissionModel/binaryAbsorptionEmission/binaryAbsorptionEmission.C
 submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C
 submodels/absorptionEmissionModel/wideBandAbsorptionEmission/wideBandAbsorptionEmission.C
+submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.C
 
 
 /* Boundary conditions */
diff --git a/src/thermophysicalModels/radiationModels/derivedFvPatchFields/radiationCoupledBase/radiationCoupledBase.C b/src/thermophysicalModels/radiationModels/derivedFvPatchFields/radiationCoupledBase/radiationCoupledBase.C
index aaa301f345..2ac8a1f9fd 100644
--- a/src/thermophysicalModels/radiationModels/derivedFvPatchFields/radiationCoupledBase/radiationCoupledBase.C
+++ b/src/thermophysicalModels/radiationModels/derivedFvPatchFields/radiationCoupledBase/radiationCoupledBase.C
@@ -29,7 +29,7 @@ License
 
 #include "mappedPatchBase.H"
 #include "fvPatchFieldMapper.H"
-#include "solidThermo.H"
+#include "radiationModel.H"
 
 // * * * * * * * * * * * * * Static Member Data  * * * * * * * * * * * * * * //
 
@@ -143,10 +143,10 @@ Foam::scalarField Foam::radiationCoupledBase::emissivity() const
 
             const polyMesh& nbrMesh = mpp.sampleMesh();
 
-            const solidThermo& thermo =
-                nbrMesh.lookupObject<solidThermo>
+            const radiation::radiationModel& radiation =
+                nbrMesh.lookupObject<radiation::radiationModel>
                 (
-                    "thermophysicalProperties"
+                    "radiationProperties"
                 );
 
             // Force recalculation of mapping and schedule
@@ -157,7 +157,13 @@ Foam::scalarField Foam::radiationCoupledBase::emissivity() const
                 nbrMesh
             ).boundary()[mpp.samplePolyPatch().index()];
 
-            scalarField emissivity(thermo.emissivity(nbrPatch.index()));
+            scalarField emissivity
+            (
+                radiation.absorptionEmission().e()().boundaryField()
+                [
+                    nbrPatch.index()
+                ]
+            );
             distMap.distribute(emissivity);
 
             return emissivity;
diff --git a/src/thermophysicalModels/radiationModels/radiationModel/opaqueSolid/opaqueSolid.C b/src/thermophysicalModels/radiationModels/radiationModel/opaqueSolid/opaqueSolid.C
new file mode 100644
index 0000000000..b4396573c6
--- /dev/null
+++ b/src/thermophysicalModels/radiationModels/radiationModel/opaqueSolid/opaqueSolid.C
@@ -0,0 +1,130 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "opaqueSolid.H"
+#include "addToRunTimeSelectionTable.H"
+#include "physicoChemicalConstants.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+    namespace radiation
+    {
+        defineTypeNameAndDebug(opaqueSolid, 0);
+
+        addToRunTimeSelectionTable
+        (
+            radiationModel,
+            opaqueSolid,
+            dictionary
+        );
+    }
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::radiation::opaqueSolid::opaqueSolid
+(
+    const volScalarField& T
+)
+:
+    radiationModel(typeName, T)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::radiation::opaqueSolid::~opaqueSolid()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+bool Foam::radiation::opaqueSolid::read()
+{
+    return radiationModel::read();
+}
+
+
+void Foam::radiation::opaqueSolid::calculate()
+{
+    // Do nothing
+}
+
+
+Foam::tmp<Foam::volScalarField> Foam::radiation::opaqueSolid::Rp() const
+{
+    return tmp<volScalarField>
+    (
+        new volScalarField
+        (
+            IOobject
+            (
+                "Rp",
+                mesh_.time().timeName(),
+                mesh_,
+                IOobject::NO_READ,
+                IOobject::NO_WRITE
+            ),
+            mesh_,
+            dimensionedScalar
+            (
+                "Rp",
+                constant::physicoChemical::sigma.dimensions()/dimLength,
+                0.0
+            )
+        )
+    );
+}
+
+
+Foam::tmp<Foam::DimensionedField<Foam::scalar, Foam::volMesh> >
+Foam::radiation::opaqueSolid::Ru() const
+{
+    return tmp<DimensionedField<scalar, volMesh> >
+    (
+        new DimensionedField<scalar, volMesh>
+        (
+            IOobject
+            (
+                "Ru",
+                mesh_.time().timeName(),
+                mesh_,
+                IOobject::NO_READ,
+                IOobject::NO_WRITE
+            ),
+            mesh_,
+            dimensionedScalar
+            (
+                "Ru", dimMass/dimLength/pow3(dimTime), 0.0
+            )
+        )
+    );
+}
+
+
+// ************************************************************************* //
diff --git a/src/thermophysicalModels/radiationModels/radiationModel/opaqueSolid/opaqueSolid.H b/src/thermophysicalModels/radiationModels/radiationModel/opaqueSolid/opaqueSolid.H
new file mode 100644
index 0000000000..d50e6f6795
--- /dev/null
+++ b/src/thermophysicalModels/radiationModels/radiationModel/opaqueSolid/opaqueSolid.H
@@ -0,0 +1,109 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::radiation::opaqueSolid
+
+Description
+    Radiation for solid opaque solids - does nothing to energy equation source
+    terms (returns zeros) but creates absorptionEmissionModel and
+    scatterModel.
+
+SourceFiles
+    opaqueSolid.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef opaqueSolid_H
+#define opaqueSolid_H
+
+#include "radiationModel.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace radiation
+{
+
+/*---------------------------------------------------------------------------*\
+                             Class opaqueSolid Declaration
+\*---------------------------------------------------------------------------*/
+
+class opaqueSolid
+:
+    public radiationModel
+{
+    // Private Member Functions
+
+        //- Disallow default bitwise copy construct
+        opaqueSolid(const opaqueSolid&);
+
+        //- Disallow default bitwise assignment
+        void operator=(const opaqueSolid&);
+
+
+public:
+
+    //- Runtime type information
+    TypeName("opaqueSolid");
+
+
+    // Constructors
+
+        //- Construct from components
+        opaqueSolid(const volScalarField& T);
+
+
+    //- Destructor
+    virtual ~opaqueSolid();
+
+
+    // Member functions
+
+        // Edit
+
+            //- Solve radiation equation(s)
+            void calculate();
+
+            //- Read radiationProperties dictionary
+            bool read();
+
+            //- Source term component (for power of T^4)
+            tmp<volScalarField> Rp() const;
+
+            //- Source term component (constant)
+            tmp<DimensionedField<scalar, volMesh> > Ru() const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace radiation
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/thermophysicalModels/radiationModels/radiationModel/radiationModel/radiationModel.C b/src/thermophysicalModels/radiationModels/radiationModel/radiationModel/radiationModel.C
index 1664d8edd1..149fb19691 100644
--- a/src/thermophysicalModels/radiationModels/radiationModel/radiationModel/radiationModel.C
+++ b/src/thermophysicalModels/radiationModels/radiationModel/radiationModel/radiationModel.C
@@ -89,7 +89,7 @@ Foam::radiation::radiationModel::radiationModel
     T_(T),
     radiation_(lookup("radiation")),
     coeffs_(subDict(type + "Coeffs")),
-    solverFreq_(readLabel(lookup("solverFreq"))),
+    solverFreq_(lookupOrDefault<label>("solverFreq", 1)),
     firstIter_(true),
     absorptionEmission_(absorptionEmissionModel::New(*this, mesh_)),
     scatter_(scatterModel::New(*this, mesh_))
@@ -158,4 +158,11 @@ Foam::tmp<Foam::fvScalarMatrix> Foam::radiation::radiationModel::Sh
 }
 
 
+const Foam::radiation::absorptionEmissionModel&
+Foam::radiation::radiationModel::absorptionEmission() const
+{
+    return absorptionEmission_();
+}
+
+
 // ************************************************************************* //
diff --git a/src/thermophysicalModels/radiationModels/radiationModel/radiationModel/radiationModel.H b/src/thermophysicalModels/radiationModels/radiationModel/radiationModel/radiationModel.H
index b9963f0f70..e8aa9f5778 100644
--- a/src/thermophysicalModels/radiationModels/radiationModel/radiationModel/radiationModel.H
+++ b/src/thermophysicalModels/radiationModels/radiationModel/radiationModel/radiationModel.H
@@ -48,6 +48,7 @@ SourceFiles
 #include "fluidThermo.H"
 #include "fvMatrices.H"
 #include "blackBodyEmission.H"
+#include "absorptionEmissionModel.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -57,7 +58,6 @@ namespace radiation
 {
 
 // Forward declaration of classes
-class absorptionEmissionModel;
 class scatterModel;
 
 /*---------------------------------------------------------------------------*\
@@ -178,6 +178,9 @@ public:
 
             //- Energy source term
             virtual tmp<fvScalarMatrix> Sh(fluidThermo& thermo) const;
+
+            //- Access to absorptionEmissionModel
+            const absorptionEmissionModel& absorptionEmission() const;
 };
 
 
diff --git a/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.C b/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.C
new file mode 100644
index 0000000000..3371c57597
--- /dev/null
+++ b/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.C
@@ -0,0 +1,256 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "greyMeanSolidAbsorptionEmission.H"
+#include "addToRunTimeSelectionTable.H"
+#include "unitConversion.H"
+#include "zeroGradientFvPatchFields.H"
+#include "basicMultiComponentMixture.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+    namespace radiation
+    {
+        defineTypeNameAndDebug(greyMeanSolidAbsorptionEmission, 0);
+
+        addToRunTimeSelectionTable
+        (
+            absorptionEmissionModel,
+            greyMeanSolidAbsorptionEmission,
+            dictionary
+        );
+/*
+        template<>
+        const char* Foam::NamedEnum
+        <
+            Foam::radiation::greyMeanSolidAbsorptionEmission::
+            radiativePropertiesNumbering,
+            3
+        >::names[] =
+        {
+            "absorptivity",
+            "emissivity",
+            "emission"
+        };
+*/
+    }
+}
+
+/*
+const Foam::NamedEnum
+    <
+        Foam::radiation::greyMeanSolidAbsorptionEmission::
+        radiativePropertiesNumbering,
+        3
+    >
+    Foam::radiation::greyMeanSolidAbsorptionEmission::
+    radiativePropertiesNumberingNames_;
+*/
+
+// * * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * //
+
+Foam::tmp<Foam::scalarField> Foam::radiation::
+greyMeanSolidAbsorptionEmission::X
+(
+    const volScalarField& Yj
+) const
+{
+    const volScalarField& T = thermo_.T();
+    const volScalarField& p = thermo_.p();
+
+    const basicMultiComponentMixture& mixture =
+        dynamic_cast<const basicMultiComponentMixture&>(thermo_);
+
+    const label mySpecieI = mixture.species()[Yj.name()];
+
+    tmp<volScalarField> tXj
+    (
+        new volScalarField
+        (
+            IOobject
+            (
+                "Xj",
+                mesh().time().timeName(),
+                mesh(),
+                IOobject::NO_READ,
+                IOobject::NO_WRITE
+            ),
+            mesh(),
+            dimensionedScalar("Xj", dimless, 0.0),
+            zeroGradientFvPatchVectorField::typeName
+        )
+    );
+
+    scalarField& Xj = tXj().internalField();
+
+    tmp<scalarField> tRhoInv(Xj);
+    scalarField& rhoInv = tRhoInv();
+
+    forAll(mixture.Y(), specieI)
+    {
+        const volScalarField& Yi = mixture.Y(specieI);
+
+        forAll(Xj, iCell)
+        {
+            rhoInv[iCell] +=
+                Yi[iCell]/mixture.rho(specieI, p[iCell], T[iCell]);
+        }
+    }
+
+    forAll(Xj, iCell)
+    {
+        Xj[iCell] = Yj[iCell]/mixture.rho(mySpecieI, p[iCell], T[iCell]);
+    }
+
+    return (Xj/rhoInv);
+}
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::radiation::greyMeanSolidAbsorptionEmission::
+greyMeanSolidAbsorptionEmission
+(
+    const dictionary& dict,
+    const fvMesh& mesh
+)
+:
+    absorptionEmissionModel(dict, mesh),
+    coeffsDict_((dict.subDict(typeName + "Coeffs"))),
+    thermo_(mesh.lookupObject<fluidThermo>("thermophysicalProperties")),
+    speciesNames_(0),
+    solidData_
+    (
+        dynamic_cast
+        <
+            const basicMultiComponentMixture&
+        >(thermo_).species().size()
+    )
+{
+    if (!isA<basicMultiComponentMixture>(thermo_))
+    {
+        FatalErrorIn
+        (
+            "radiation::greyMeanSolidAbsorptionEmission::"
+            "greyMeanSolidAbsorptionEmission"
+            "("
+                "const dictionary&, "
+                "const fvMesh&"
+            ")"
+        )   << "Model requires a multi-component thermo package"
+            << abort(FatalError);
+    }
+
+    label nFunc = 0;
+    const dictionary& functionDicts = dict.subDict(typeName + "Coeffs");
+
+    forAllConstIter(dictionary, functionDicts, iter)
+    {
+        // safety:
+        if (!iter().isDict())
+        {
+            continue;
+        }
+        const word& key = iter().keyword();
+        speciesNames_.insert(key, nFunc);
+        const dictionary& dict = iter().dict();
+        dict.lookup("a") >> solidData_[nFunc][absorptivity];
+        dict.lookup("e") >> solidData_[nFunc][emissivity];
+
+        nFunc++;
+    }
+}
+
+// * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
+
+Foam::radiation::greyMeanSolidAbsorptionEmission::
+~greyMeanSolidAbsorptionEmission()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+Foam::tmp<Foam::volScalarField>
+Foam::radiation::greyMeanSolidAbsorptionEmission::
+calc(const label property) const
+{
+    const basicMultiComponentMixture& mixture =
+        dynamic_cast<const basicMultiComponentMixture&>(thermo_);
+
+    tmp<volScalarField> ta
+    (
+        new volScalarField
+        (
+            IOobject
+            (
+                "a",
+                mesh().time().timeName(),
+                mesh(),
+                IOobject::NO_READ,
+                IOobject::NO_WRITE
+            ),
+            mesh(),
+            dimensionedScalar("a", dimless/dimLength, 0.0),
+            zeroGradientFvPatchVectorField::typeName
+        )
+    );
+
+    scalarField& a = ta().internalField();
+
+    forAllConstIter(HashTable<label>, speciesNames_, iter)
+    {
+        if (mixture.contains(iter.key()))
+        {
+            const volScalarField& Y = mixture.Y(iter.key());
+            a += solidData_[iter()][property]*X(Y);
+        }
+    }
+
+    ta().correctBoundaryConditions();
+    return ta;
+}
+
+
+Foam::tmp<Foam::volScalarField>
+Foam::radiation::greyMeanSolidAbsorptionEmission::eCont
+(
+    const label bandI
+) const
+{
+   return calc(emissivity);
+}
+
+
+Foam::tmp<Foam::volScalarField>
+Foam::radiation::greyMeanSolidAbsorptionEmission::aCont
+(
+    const label bandI
+) const
+{
+   return calc(absorptivity);
+}
+
+// ************************************************************************* //
diff --git a/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.H b/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.H
new file mode 100644
index 0000000000..9e2586a88b
--- /dev/null
+++ b/src/thermophysicalModels/radiationModels/submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.H
@@ -0,0 +1,158 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011-2012 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::radiation::greyMeanSolidAbsorptionEmission
+
+Description
+    greyMeanSolidAbsorptionEmission radiation absorption and emission
+    coefficients for continuous phase
+
+    The coefficients for the species in the Look up table have to be specified
+    for use in moles x P [atm], i.e. (k[i] = species[i]*p*9.869231e-6).
+
+    The coefficients for CO and soot or any other added are multiplied by the
+    respective mass fraction being solved
+
+    All the species in the dictionary need either to be in the look-up table or
+    being solved. Conversely, all the species solved do not need to be included
+    in the calculation of the absorption coefficient
+
+    The names of the species in the absorption dictionary must match exactly the
+    name in the look-up table or the name of the field being solved
+
+SourceFiles
+    greyMeanSolidAbsorptionEmission.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef greyMeanSolidAbsorptionEmission_H
+#define greyMeanSolidAbsorptionEmission_H
+
+#include "absorptionEmissionModel.H"
+#include "fluidThermo.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+namespace radiation
+{
+
+/*---------------------------------------------------------------------------*\
+                 Class greyMeanSolidAbsorptionEmission Declaration
+\*---------------------------------------------------------------------------*/
+
+class greyMeanSolidAbsorptionEmission
+:
+    public absorptionEmissionModel
+{
+
+private:
+
+    // Private data
+
+        //- Enumering of radiative properties
+        enum radiativeProperties
+        {
+            absorptivity,
+            emissivity
+        };
+
+        //- Absorption model dictionary
+        dictionary coeffsDict_;
+
+        //- SLG thermo package
+        const fluidThermo& thermo_;
+
+        //- Hash table of species names
+        HashTable<label> speciesNames_;
+
+        //- List of solid species data
+        List<FixedList<scalar, 2> > solidData_;
+
+
+    // Private member functions
+
+       //- Calculate the volumetric fraction of Yj
+       tmp<scalarField> X(const volScalarField& Yj) const;
+
+       //- Calculate the property mixing
+       tmp<volScalarField> calc(const label) const;
+
+
+public:
+
+    //- Runtime type information
+    TypeName("greyMeanSolidAbsorptionEmission");
+
+
+    // Constructors
+
+        //- Construct from components
+        greyMeanSolidAbsorptionEmission
+        (
+            const dictionary& dict,
+            const fvMesh& mesh
+        );
+
+
+    //- Destructor
+    virtual ~greyMeanSolidAbsorptionEmission();
+
+
+    // Member Operators
+
+        // Access
+
+            // Absorption coefficient
+
+                //- Absorption coefficient for continuous phase
+                tmp<volScalarField> aCont(const label bandI = 0) const;
+
+
+            // Emission coefficient
+
+                //- Emission coefficient for continuous phase
+                tmp<volScalarField> eCont(const label bandI = 0) const;
+
+
+    // Member Functions
+
+        inline bool isGrey() const
+        {
+            return true;
+        }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace radiation
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //
diff --git a/src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.C b/src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.C
index 4cdec064bf..938244b592 100644
--- a/src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.C
+++ b/src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.C
@@ -207,4 +207,16 @@ Foam::scalar Foam::SpecieMixture<MixtureType>::alphah
 }
 
 
+template<class MixtureType>
+Foam::scalar Foam::SpecieMixture<MixtureType>::rho
+(
+    const label speciei,
+    const scalar p,
+    const scalar T
+) const
+{
+    return this->getLocalThermo(speciei).rho(p, T);
+}
+
+
 // ************************************************************************* //
diff --git a/src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.H b/src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.H
index ed237ef446..a954119dac 100644
--- a/src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.H
+++ b/src/thermophysicalModels/reactionThermo/mixtures/SpecieMixture/SpecieMixture.H
@@ -176,6 +176,14 @@ public:
                 const scalar p,
                 const scalar T
             ) const;
+
+            //- Density [kg/m3]
+            virtual scalar rho
+            (
+                const label specieI,
+                const scalar p,
+                const scalar T
+            ) const;
 };
 
 
diff --git a/src/thermophysicalModels/reactionThermo/mixtures/basicMultiComponentMixture/basicMultiComponentMixture.H b/src/thermophysicalModels/reactionThermo/mixtures/basicMultiComponentMixture/basicMultiComponentMixture.H
index 80a9e7f056..ea712fc22f 100644
--- a/src/thermophysicalModels/reactionThermo/mixtures/basicMultiComponentMixture/basicMultiComponentMixture.H
+++ b/src/thermophysicalModels/reactionThermo/mixtures/basicMultiComponentMixture/basicMultiComponentMixture.H
@@ -227,6 +227,14 @@ public:
                 const scalar p,
                 const scalar T
             ) const = 0;
+
+            //- Density [kg/m3]
+            virtual scalar rho
+            (
+                const label specieI,
+                const scalar p,
+                const scalar T
+            ) const = 0;
 };
 
 
diff --git a/tutorials/combustion/fireFoam/les/oppositeBurningPanels/constant/panelRegion/radiationProperties b/tutorials/combustion/fireFoam/les/oppositeBurningPanels/constant/panelRegion/radiationProperties
new file mode 100644
index 0000000000..c60289f7b4
--- /dev/null
+++ b/tutorials/combustion/fireFoam/les/oppositeBurningPanels/constant/panelRegion/radiationProperties
@@ -0,0 +1,70 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+| =========                 |                                                 |
+| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
+|  \\    /   O peration     | Version:  dev                                   |
+|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
+|    \\/     M anipulation  |                                                 |
+\*---------------------------------------------------------------------------*/
+
+FoamFile
+{
+    version     2.0;
+    format      ascii;
+    class       dictionary;
+    location    "constant";
+    object      radiationProperties;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+radiation       off;
+
+radiationModel  none;
+
+noRadiation
+{
+}
+
+P1Coeffs
+{
+}
+
+fvDOMCoeffs
+{
+}
+
+// Number of flow iterations per radiation iteration
+solverFreq 10;
+
+//absorptionEmissionModel constantAbsorptionEmission;
+absorptionEmissionModel greyMeanSolidAbsorptionEmission;
+
+
+greyMeanSolidAbsorptionEmissionCoeffs
+{
+
+    v
+    {
+        sigmaS      0.0;
+        a           0.0;  //opaque
+        emissivity  0.17;
+    }
+
+    char
+    {
+        sigmaS      0.0;
+        a           0.0;
+        emissivity  0.85;
+    }
+
+}
+
+scatterModel    constantScatter;
+
+constantScatterCoeffs
+{
+    sigma           sigma [ 0 -1 0 0 0 0 0 ] 0;
+    C               C [ 0 0 0 0 0 0 0 ] 0;
+}
+
+
+// ************************************************************************* //