|
[Sponsors] |
November 23, 2014, 00:07 |
No convergence - SimpleFOAM
|
#1 |
Senior Member
SinaJ
Join Date: Nov 2009
Posts: 136
Rep Power: 16 |
Hi FOAMers! :-)
I'm trying to simulate a flow with 2 inlets and 2 outlets. I can't reach a better convergence! Please see the attached for the residuals. fvSolution Code:
solvers { p { solver GAMG; tolerance 1e-07; // was 1e-6 relTol 0.05; // was 0.1 smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration on; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; mergeLevels 1; } "(U|k|epsilon|R|nuTilda)" { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-06; // was 1e-5 relTol 0.1; // was 0.1 } } SIMPLE { nNonOrthogonalCorrectors 0; // was 0 ! residualControl { p 1e-5; U 1e-5; "(k|epsilon|omega)" 1e-5; } } relaxationFactors { fields { p 0.3; } equations { U 0.5; k 0.7; epsilon 0.7; R 0.7; nuTilda 0.7; } } fvSchemes Code:
ddtSchemes { default steadyState; } gradSchemes { default edgeCellsLeastSquares; } divSchemes { default none; div(phi,U) bounded Gauss QUICK; div(phi,k) bounded Gauss upwind; div(phi,epsilon) bounded Gauss upwind; div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } fluxRequired { default no; p ; } |
|
November 23, 2014, 09:49 |
|
#2 |
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 21 |
What are your boundary conditions?
|
|
November 23, 2014, 14:01 |
|
#3 |
Senior Member
SinaJ
Join Date: Nov 2009
Posts: 136
Rep Power: 16 |
Two velocity inlets (fixed velocity, zero gradient pressure) . Two averaged pressure outlets (0kpa and 10kpa, zero-gradient velocity). The walls are no-slip (zero velocity, zero gradient pressure).
|
|
November 23, 2014, 17:42 |
|
#4 |
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 650
Rep Power: 21 |
Does the two different pressure outlets make physically sense? Meaning: Is the pressure drop between the two outlets so high?
|
|
November 23, 2014, 18:04 |
|
#5 |
Senior Member
SinaJ
Join Date: Nov 2009
Posts: 136
Rep Power: 16 |
Yes It has a physical meaning. and I could get converged solution with Fluent but I want to get the converged solution with OpenFOAM.
|
|
November 27, 2014, 07:01 |
|
#6 |
Senior Member
|
Hmmm, I would rather check these conditions first. Try with these and check your convergence...
gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) Gauss upwind; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div((nuEff*dev(T(grad(U))))) Gauss linear; } - Best Luck! |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Convergence | Centurion2011 | FLUENT | 48 | June 14, 2022 23:29 |
convergence of QUICK scheme - simpleFoam | Luis Batista | OpenFOAM Running, Solving & CFD | 10 | May 11, 2013 17:35 |
Convergence and steady state using simpleFoam | sfigato | OpenFOAM Running, Solving & CFD | 0 | February 8, 2013 04:14 |
Force can not converge | colopolo | CFX | 13 | October 4, 2011 22:03 |
Definition of convergence criterion in simpleFoam | titio | OpenFOAM Running, Solving & CFD | 1 | February 6, 2010 01:34 |