[smog]

Hello all. I have attached a pic of the residuals. I cannot seem to get my residuals below 1e-3 for Ux and Uy, although my Uz residuals are well below 1e-3. The pressure residuals are also 1e-4<p_res<1e-3.

I am running a multiphase 3D cylindrical jet simulation. The mesh gets coarser spanwise (x and y directions), while streamwise (z direction) the mesh is much finer. I am pasting my fvSchemes and fvSolutions here as well.

fvSchemes:

ddtSchemes
{
default CrankNicolson 0.75;//Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
//div(rhoPhi,U) Gauss vanLeerV;
//div(phi,alpha) Gauss vanLeer;
//div(phi,age) Gauss linearUpwind grad(U);
//div(phirb,alpha) Gauss vanLeer;
//div(phi,k) Gauss upwind;
//div(phi,omega) Gauss upwind;
//div(phi,R) Gauss upwind;
//div(R) Gauss linear;
//div(phi,nuTilda) Gauss upwind;
//div((nuEff*dev(T(grad(U))))) Gauss linear;
//div((muEff*dev(T(grad(U))))) Gauss linear;



div(rhoPhi,U) Gauss linearUpwind grad(U); //upwind; //limitedLinearV 1;

div(phi,alpha) Gauss vanLeer;

div(phirb,alpha) Gauss linearUpwind grad(U);//interfaceCompression;

div(phi,k) Gauss limitedLinear 1;//upwind

div((nuEff*dev(T(grad(U))))) Gauss linear;
div((muEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p_rgh;
pcorr;
alpha.oil;
}


fvSolution:
solvers
{
pcorr
{
solver PCG;
preconditioner
{
preconditioner GAMG;
tolerance 1e-05;
relTol 0;
smoother DICGaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration false;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
tolerance 1e-05;
relTol 0;
maxIter 100;
}

pcorrFinal
{
solver PCG;
preconditioner
{
preconditioner GAMG;
tolerance 1e-05;
relTol 0;
smoother DICGaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration false;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

tolerance 1e-05;
relTol 0;
maxIter 100;
}

p_rgh
{
solver GAMG;
tolerance 1e-07;
relTol 0.01;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

p_rghFinal
{
solver PCG;
preconditioner
{
preconditioner GAMG;
tolerance 1e-07;
relTol 0;
nVcycles 2;
smoother DICGaussSeidel;
nPreSweeps 2;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

tolerance 1e-07;
relTol 0;
maxIter 20;
}

"alpha.oil.*"
{
nAlphaCorr 2;
nAlphaSubCycles 1;
alphaOuterCorrectors yes;
cAlpha 1;

MULESCorr yes;
nLimiterIter 3;

solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-8;
relTol 0;
}

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

}

"(U|k|omega|age)Final"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
}

nuTilda
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0;
};

"(nuTilda)Final"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-08;
relTol 0;
};
}

PIMPLE
{
momentumPredictor yes;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
nAlphaCorr 1;
nAlphaSubCycles 4;
cAlpha 1;
//correctPhi no;
//pRefCell 0;
//pRefValue 0;
}

-------------------------------------------------------------------

Please let me know what I can do to get the residuals below 1e-4. I have run 4 flow-throughs of simulations time.