|
[Sponsors] |
June 5, 2016, 12:16 |
Non-converging magnitudes Cylinder case
|
#1 |
New Member
Max Vorstadt
Join Date: May 2016
Posts: 28
Rep Power: 10 |
Hi! I'm running a kwSST simpleFoam simulation with a cylinder and a plate and apparently I'm having trouble with the convergence of the different magnitudes. The meshing process has been done with snappyHexMesh (all 3 steps run succesfully). The geometry of the case and a previous "lower-level" mesh can be seen in the attachments, as well as the plot of the residuals.
Does anybody have any idea why the magnitudes are not converging? Maybe a fvSchemes problem, a Boundary Conditions problem, the need to run potentialFoam before to initialise the field...? The different relevant files are attached below as well. Thanks in advance! U: Code:
dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type fixedValue; value uniform (26 0 0); } outlet { type inletOutlet; inletValue uniform (26 0 0); } lateralWall { type symmetry; } upperWall { type symmetry; } lowerWall { type symmetry; } Cylinder { type fixedValue; value uniform (0 0 0); } Plate_round { type fixedValue; value uniform (0 0 0); } } Code:
dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } upperWall { type symmetry; } lowerWall { type symmetry; } lateralWall { type symmetry; } Cylinder { type zeroGradient; } Plate_round { type zeroGradient; } } Code:
dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedValue; value uniform 0.02535; } outlet { type zeroGradient; } upperWall { type symmetry; } lowerWall { type symmetry; } lateralWall { type symmetry; } Cylinder { type kqRWallFunction; value uniform 0; } Plate_round { type kqRWallFunction; value uniform 0; } } Code:
dimensions [0 0 -1 0 0 0 0]; internalField uniform 162.5; boundaryField { inlet { type fixedValue; value uniform 162.5; } outlet { type zeroGradient; } upperWall { type symmetry; } lowerWall { type symmetry; } lateralWall { type symmetry; } Cylinder { type omegaWallFunction; value uniform 162.5; } Plate_round { type omegaWallFunction; value uniform 162.5; } } Code:
dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type calculated; value uniform 5e-4; } outlet { type zeroGradient; } upperWall { type symmetry; } lowerWall { type symmetry; } lateralWall { type symmetry; } Cylinder { type nutUSpaldingWallFunction; value uniform 0; } Plate_round { type nutUSpaldingWallFunction; value uniform 0; } } Code:
ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; grad(U) cellLimited Gauss linear 1; } divSchemes { default none; div(phi,U) bounded Gauss linearUpwindV grad(U); div(phi,k) bounded Gauss upwind; div(phi,omega) bounded Gauss upwind; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } wallDist { method meshWave; } Code:
solvers { p { solver GAMG; tolerance 1e-8; relTol 0.01; smoother GaussSeidel; nPreSweeps 0; nPostSweeps 2; cacheAgglomeration on; agglomerator faceAreaPair; nCellsInCoarsestLevel 10; mergeLevels 1; } Phi { $p; } U { solver smoothSolver; smoother GaussSeidel; tolerance 1e-8; relTol 0.1; nSweeps 1; } k { solver smoothSolver; smoother GaussSeidel; tolerance 1e-8; relTol 0.1; nSweeps 1; } omega { solver smoothSolver; smoother GaussSeidel; tolerance 1e-8; relTol 0.1; nSweeps 1; } } SIMPLE { nNonOrthogonalCorrectors 5; consistent yes; residualControl { p 1e-7; U 1e-7; "(k|epsilon|omega)" 1e-7; } } potentialFlow { nNonOrthogonalCorrectors 10; } relaxationFactors { equations { U 0.9; k 0.7; omega 0.7; } } cache { grad(U); } |
|
Tags |
convergence, cylinder, fvschemes, kwsstmodel, turbulence |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
flow over a cylinder urgent! | kevin | FLUENT | 8 | August 11, 2015 14:00 |
Problem solving a cylinder case 2d. Simple Foam | vbesteban | OpenFOAM Running, Solving & CFD | 0 | October 16, 2014 17:51 |
How to Initialise my LES case using my RANS case is there any utility for it ? | Alhasan | OpenFOAM Running, Solving & CFD | 2 | May 10, 2014 01:14 |
Compare the pressure of potential flow and a turbulent case of cylinder | ooo | OpenFOAM Running, Solving & CFD | 0 | August 2, 2013 08:25 |
Turbulent steady flow around a circular cylinder | Mirek Kabacinski | FLUENT | 0 | July 23, 2003 19:40 |