[Hillie]

Dear Foamers,

I have been working on a nozzle simulation with rhoCentralFoam. I have a CFX solution which I am currently using as a baseline.

The core-flow matches up beautifully with the CFX results. The wall pressure is also looks very nice, and compares very well. My current issue is with the wall temperature.
I have a large drop in wall temperature in my throat, which recovers after my throat to about 3 degrees lower then my CFX baseline.
My temperature B.C at the well is zero gradient, i.e adiabatic wall. If I compute the adiabatic wall temperature for my case it look like the CFX solution is right.

The wall shear is also over predicted w.r.t the CFX solution. Which would suggest an over prediction of the velocity wall gradient. Yplus is of order 1, and checkmesh also checks out.

My fvSchemes file is currently set to the following:

fluxScheme      Tadmor;

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
     default               leastSquares;
    grad(U)                cellLimited leastSquares 1;
    grad(rho)              cellLimited leastSquares 1;  
    grad(omega)            cellLimited leastSquares 1;  
    grad(k)                cellLimited leastSquares 1;
    
}

divSchemes
{

    default          none;
    div(tauMC)       Gauss linear;
    div(phi,U)       Gauss linear;
    div(phid,p)      Gauss linear;
    div(phi,K)       Gauss linear;
    div(phi,h)       Gauss Gamma 1;
    div(phi,k)       Gauss Gamma 1;
    div(phi,epsilon) Gauss linear;
    div(phi,R)       Gauss Gamma 1;
    div(phi,omega)   Gauss Gamma 1;
    div((rho*R))     Gauss Gamma 1;
    div(R)           Gauss Gamma 1;
    div(U)           Gauss linear;
    div((muEff*dev2(T(grad(U))))) Gauss linear;
    div(yPhi,yWall)  Gauss linear ;


}

laplacianSchemes
{
    default         none;
   laplacian(muEff,U) Gauss Gamma phi 1 corrected;
    laplacian(alphaEff,e) Gauss Gamma phi 1 corrected;
    laplacian((rho*DomegaEff),omega)  Gauss Gamma phi 1 corrected;
    laplacian((rho*DkEff),k) Gauss Gamma phi 1 corrected;
    laplacian(yPsi)   Gauss Gamma phi 1 corrected;
    laplacian(yWall)   Gauss Gamma phi 1 corrected;

 
}

interpolationSchemes
{
    default            linear;
    reconstruct(rho)   Gamma 1;
    reconstruct(U)     GammaV 1;
    reconstruct(T)     Gamma 1;
    interpolate(rhoU)  Gamma pos 1;
    interpolate(rho)   Gamma pos 1;

}

snGradSchemes
{
    default         corrected;
}


wallDist
{
      method advectionDiffusion;

        // Optional entry enabling the calculation
        // of the normal-to-wall field
        nRequired false;

         advectionDiffusionCoeffs
         {
             method    Poisson;
             epsilon   0.4;
             tolerance 1e-3;
             maxIter   20;
         }

}

If anybody has any idea's on what I could use to get the wall temperature up those last couple of degrees that would be grand.