[infrared]

Hello everyone,

I successfully ran a combustion case with OpenFOAM using C8H16 as the fuel.

I have four boundary conditions (inlet, outlet, top wall, and bottom wall) where C8H16 is injected from the inlet and impinges on the bottom wall, ensuring that no flow exits the domain. The walls are set to a temperature of 1000K.

The simulation runs well, but I cannot visualize the 'T' field in the results. I have the 0/T, fvSchemes, fvSolution, and controlDict files set up.

Please review them below and provide suggestions to improve my simulation. Your help is very valuable to me. Thank you.
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2312 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
arch "LSB;label=32;scalar=64";
class volScalarField;
location "0";
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 0 0 0 0 0];

internalField uniform 273;

boundaryField
{
bottom_wall
{
type fixedValue;
value uniform 273;
}
inlet
{
type fixedValue;
value uniform 273;
}
top_wall
{
type fixedValue;
value uniform 273;
}
outlet
{
type zeroGradient;
}
}


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

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1812 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default CrankNicolson 1;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,C1) Gauss linearUpwind grad(C1);
div(phi,C2) Gauss linearUpwind grad(C2);
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}


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

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1812 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
p
{
solver GAMG;
tolerance 1e-06;
relTol 0.1;
smoother GaussSeidel;
}

pFinal
{
$p;
smoother DICGaussSeidel;
tolerance 1e-06;
relTol 0;
}

U
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-05;
relTol 0;
}

UFinal
{
$U;
}

T
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-06;
relTol 0;
}

TFinal
{
$T;
}

"(C1|C2|...)"
{
solver smoothSolver;
smoother GaussSeidel;
tolerance 1e-05;
relTol 0;
minIter 1;
}
}

PIMPLE
{
nOuterCorrectors 1;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
}


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

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2312 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application pimpleFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 500;

deltaT 0.001;

adjustTimeStep true;

writeControl adjustable;

writeInterval 10;

writeControl adjustableRunTime;

purgeWrite 0;

writeFormat ascii;

writePrecision 12;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep yes;

maxCo 0.8;

DebugSwitches
{
SolverPerformance 0;
}

functions
{
sumMassFractions
{
type fieldExpression;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl timeStep;
outputInterval 10;
expression "C8H16 + CO2 + H2O + N2 + O2";
result Ysum;
}
}

writeObjects
(
U
p
T
C1
C2
C8H16
CO
NO
// other species fields
);

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