|
[Sponsors] |
January 19, 2016, 10:33 |
Convergence with Interfoam
|
#1 |
Member
shashank moghe
Join Date: Feb 2015
Posts: 32
Rep Power: 11 |
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. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Convergence | Centurion2011 | FLUENT | 48 | June 14, 2022 23:29 |
k-e & GAMG interFoam Schemitisation Stability Issue | JFM | OpenFOAM Running, Solving & CFD | 3 | December 1, 2015 05:58 |
Convergence Test on InterFoam | rcarmi | Main CFD Forum | 0 | December 18, 2014 15:20 |
Force can not converge | colopolo | CFX | 13 | October 4, 2011 22:03 |
Convergence of CFX field in FSI analysis | nasdak | CFX | 2 | June 29, 2009 01:17 |