[thiagopl]

Hi all,

I'm running a case with a non-reacting gas mixture. The gases are considered participating media, and I'm using the P1 radiation model with the greyMean absorption model to compute the absorption coefficients. My radiationProperties file is configured as follows:

radiation 	on;

radiationModel  P1;
solverFreq 1;
absorptionEmissionModel greyMean;
greyMeanCoeffs
{
    lookUpTableFileName     "speciesTable";

    EhrrCoeff                0.0;

    CO2
    {
        Tcommon         200;   // Common Temp
        invTemp         true;   // Is the polynomial using inverse temperature.
        Tlow            200;   // Low Temp
        Thigh           2500;  // High Temp

        loTcoeffs       // coefss for T < Tcommon
        (
            0 0 0 0 0 0
        );
        hiTcoeffs        // coefss for T > Tcommon
        (
            18.741 -121.31e3  273.5e6  -194.05e9  56.31e12  -5.8169e15
        );
    }

    H2O
    {
        Tcommon         200;
        invTemp         true;
        Tlow            200;
        Thigh           2500;

        loTcoeffs
        (
            0 0 0 0 0 0
        );
        hiTcoeffs
        (
            -0.23093  -1.12390e3  9.4153e6  -2.99885e9  0.51382e12  -1.868e10
        );
    }

scatterModel	none;
sootModel	none;

The foamInfo greyMean says that "All the species in the dictionary need either to be in the look-up table or being solved", once the case is a non-reacting case, I don't have any specie neing solved, so I want to specify them in the "speciesTable".

According to the foamInfo greyMean documentation, "All the species in the dictionary need either to be in the look-up table or be solved."
Since this is a non-reacting case, I am not solving for any species fields. Therefore, I want to define them in the speciesTable instead.

My question is:
What is the required structure and format of the speciesTable file?

P.S.: I can successfully run the case by defining species properties and molar fractions in the physicalProperties file and 0/ folder, respectivley. However, my goal is to include soot absorption in the calculation without having to consider soot into the thermodynamics of the system (properties, mass fraction...).

[zd7s18533]

Hi thiagopl
You can find more information in this file:
OpenFOAM-10/src/radiationModels/absorptionEmissionModels/interpolationLookUpTable.H

I think it may have the following form:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 10 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/

FoamFile
{
format ascii;
class dictionary;
location "constant";
object radiationProperties;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

fileName "Soot";
fields
(
{
N 2;
max 1000;
min 300;
name Temperature;
}
);

output
(
{name Temperature;}
{name largeSoot;}
{name smallSoot;}
);

values
(
3(300 500 1000) // temperature

3(0.3 0.2 0.1 ) //for large soot

3(0.1 0.2 0.3 ) //for small soot
);
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


When interpolation is performed, the table will calculate the value of large soot and small soot based on the temperature