[andrewmichael]

Hello everyone,

I am trying to simulate the flow through a diffuser using OpenFOAM. After the simulation is complete, paraview shows the values of pressure being 1e-38. I do not have a clue why this is happening because with the exact fvSchemes and fvSolutions it did not show this.

The fvSolutions and Schemes are displayed below:

fvSolutions:
solvers
{

Phi
{
solver GAMG;
smoother DIC;
tolerance 1e-06;
relTol 0.1;
}

p
{
solver GAMG;
tolerance 1e-06;
relTol 0.1;
smoother GaussSeidel;
$Phi

}


"(U|k|epsilon|omega|f|v2)"
{
solver smoothSolver;
smoother symGaussSeidel;
// solver GAMG;
// smoother GaussSeidel;
tolerance 1e-06;
relTol 0.1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
consistent yes;

residualControl
{
p 1e-6;
U 1e-6;
"(k|epsilon|omega|f|v2)" 1e-6;
}
}

relaxationFactors
{
equations
{
"U.*" 0.95; // 0.9 is more stable but 0.95 more convergent
"p.*" 0.95; // 0.9 is more stable but 0.95 more convergent
"k.*" 0.95; // 0.9 is more stable but 0.95 more convergent
"epsilon.*" 0.95; // 0.9 is more stable but 0.95 more convergent
".*" 0.95; // 0.9 is more stable but 0.95 more convergent
}
}



potentialFlow
{
nNonOrthogonalCorrectors 3;
}


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


fvSchemes

ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linearUpwind grad(U);
div(phi,k) Gauss linearUpwind grad(k);
div(phi,epsilon) Gauss linearUpwind grad(epsilon);
div(phi,omega) Gauss linearUpwind grad(omega);
div(phi,v2) Gauss upwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
div(nonlinearStress) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

wallDist
{
method meshWave;
}

Kind regards,

Andrei