[Gbiz123]

Hi, I am new to this community, CFD, and OpenFOAM. I'm mostly self taught so far, but I'm planning on studying it further in grad school.

I am trying to model a flow through a jet nozzle using OpenFOAM on Windows Subsystem for Linux. I successfully ran simpleFoam. The solver I want to try next, rhoSimpleFoam, quits after reading the thermophysicalProperties without sending an error message. I'm not sure what could be causing this. Any help would be appreciated. I've included screenshots of the problem, my mesh, and code blocks of my dictionaries.

Screenshot of the problem:
https://imgur.com/a/5pCYKfz

Screenshot of mesh:
https://imgur.com/eFe7LJC

thermophysicalProperties:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo          hConst; 
    equationOfState perfectGas; 
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture // CO2 at room temperature (1293 K) estimated temperature for jet engine exhaust
{
    specie
    {
        molWeight   44.01;
    }
    thermodynamics
    {
        Cp          1298; //J/kg K
        Hf           0;
    }
    transport
    {
        As          1.67212e-6; //Sutherland coefficient
        Ts          170.672;      //sutherland temperature
    }
}

transportProperties:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
transportModel  Newtonian;

nu              .7424;


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

fvSolution:

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

    "(U|k|omega|e|epsilon)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       1e-6;
        relTol          0.1;
        nSweeps         1;
    }
}

SIMPLE
{
    residualControl
    {
        p               1e-4;
        U               1e-4;
        "(k|omega|e|epsilon)"   1e-4;
    }

    nNonOrthogonalCorrectors 0;
    pMinFactor      0.1;
    pMaxFactor      2;
}

relaxationFactors
{
    fields
    {
        p               0.01;
        rho             0.01;
    }
    equations
    {
        U               0.5;
        e         0.1;
        epsilon         0.1;
        k               0.1;
        h               0.1;
        omega           0.1;
    }
}

fvSchemes

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;

    limited         cellLimited Gauss linear 1;
    grad(U)         $limited;
    grad(k)         $limited;
    grad(omega)     $limited;
}

divSchemes
{
    default         none;

    div(phi,U)      bounded Gauss linearUpwind limited;

    turbulence      bounded Gauss upwind;
    energy          bounded Gauss linearUpwind limited;

    div(phi,k)      $turbulence;
    div(phi,omega)  $turbulence;

    div(phi,e)      $energy;
    div(phi,K)      $energy;
    div(phi,Ekp)    $energy;

    div(phid,p)     Gauss upwind;
    div((phi|interpolate(rho)),p)  bounded Gauss upwind;

    div(((rho*nuEff)*dev2(T(grad(U)))))    Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}
wallDist
{
    method meshWave;
}

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

fvOptions:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

limitT
{
    type       limitTemperature;
    min        0;
    max        3000;
    selectionMode all;
}

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

Cheers
Greg