|
[Sponsors] |
March 11, 2015, 11:32 |
Convergence with Simplefoam
|
#1 |
Member
Rubén
Join Date: Oct 2014
Location: Munich
Posts: 47
Rep Power: 11 |
Hi Foamers!
I don't know why, but I am having big problems in my case. I don't know why but simplefoam don't converges. I changed the relaxation factors, did the checkMesh... And it doesn't works! With a diferent type of refinement I obtain a converged solution, but with this refinement not. And checkMesh says that my mesh is OK as you can see. I don't knwo really what to do ! I put bounded Gauss upwind schemes for every variable, I think it's the best in this case. I attach my log file an some files more! I will be pleased for every answer... I'm a bit desesperate Thank you very much in advance! End of log fie: Code:
Time = 893 smoothSolver: Solving for Ux, Initial residual = 0.202157, Final residual = 0.0112941, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.515388, Final residual = 0.0439188, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.380651, Final residual = 0.00752332, No Iterations 1 GAMG: Solving for p, Initial residual = 1.18088e-07, Final residual = 6.93975e-08, No Iterations 1 time step continuity errors : sum local = 2.89175e+12, global = -4.73148e+10, cumulative = -4.73148e+10 smoothSolver: Solving for epsilon, Initial residual = 2.37023e-10, Final residual = 2.37023e-10, No Iterations 0 smoothSolver: Solving for k, Initial residual = 0.263584, Final residual = 0.00591683, No Iterations 1 bounding k, min: -5.67569e+30 max: 4.98836e+36 average: 5.91692e+31 ExecutionTime = 4970.94 s ClockTime = 5152 s Time = 894 smoothSolver: Solving for Ux, Initial residual = 0.202334, Final residual = 0.00481158, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.205201, Final residual = 0.0113306, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.0962262, Final residual = 0.00346977, No Iterations 2 GAMG: Solving for p, Initial residual = 1.12762e-18, Final residual = 1.12762e-18, No Iterations 0 time step continuity errors : sum local = 1.3093e+59, global = -9.35364e+42, cumulative = -9.35364e+42 smoothSolver: Solving for epsilon, Initial residual = 6.84338e-11, Final residual = 6.84338e-11, No Iterations 0 bounding epsilon, min: 1.39235e-19 max: 3.27976e+57 average: 6.52687e+51 smoothSolver: Solving for k, Initial residual = 0.87113, Final residual = 0.042468, No Iterations 1 bounding k, min: -7.32165e+31 max: 1.33669e+51 average: 3.71993e+45 ExecutionTime = 4974.96 s ClockTime = 5156 s Time = 895 smoothSolver: Solving for Ux, Initial residual = 8.65286e-09, Final residual = 8.65286e-09, No Iterations 0 smoothSolver: Solving for Uy, Initial residual = 3.91018e-10, Final residual = 3.91018e-10, No Iterations 0 smoothSolver: Solving for Uz, Initial residual = 2.95985e-08, Final residual = 2.44359e-09, No Iterations 1 GAMG: Solving for p, Initial residual = 1.27174e-20, Final residual = 1.27174e-20, No Iterations 0 time step continuity errors : sum local = 6.53034e+56, global = -3.47895e+40, cumulative = -9.38843e+42 smoothSolver: Solving for epsilon, Initial residual = 2.62477e-07, Final residual = 1.58738e-08, No Iterations 2 bounding epsilon, min: -5.52704e+76 max: 4.64178e+77 average: 1.9647e+72 smoothSolver: Solving for k, Initial residual = 0.0205769, Final residual = 0.0011484, No Iterations 3 bounding k, min: -6.71005e+39 max: 7.39668e+60 average: 1.24645e+55 ExecutionTime = 4979.05 s ClockTime = 5160 s Time = 896 smoothSolver: Solving for Ux, Initial residual = 0.98904, Final residual = 0.0714643, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.996664, Final residual = 0.0839625, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.616932, Final residual = 0.0181629, No Iterations 2 GAMG: Solving for p, Initial residual = 0.0991486, Final residual = 1.02504e-05, No Iterations 1 time step continuity errors : sum local = 1.2079e+57, global = 5.15769e+41, cumulative = -8.87266e+42 Code:
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.3.0-f5222ca19ce6 Exec : checkMesh Date : Mar 11 2015 Time : 17:27:58 Host : "usuario-SATELLITE-P50-A-14G" PID : 6239 Case : /home/usuario/Frisbees/Cilindro/DiscraftUltraStar_kepsilon_0.1vr nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 3601375 faces: 10206962 internal faces: 10062872 cells: 3308975 faces per cell: 6.12571 boundary patches: 6 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 3059355 prisms: 12075 wedges: 0 pyramids: 0 tet wedges: 58 tetrahedra: 0 polyhedra: 237487 Breakdown of polyhedra by number of faces: faces number of cells 4 11685 5 7717 6 17983 7 84921 8 25746 9 75748 10 150 11 17 12 8944 15 4552 18 24 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology inlet 76 97 ok (non-closed singly connected) outlet 64 81 ok (non-closed singly connected) top 96 117 ok (non-closed singly connected) bottom 96 117 ok (non-closed singly connected) frontandback 192 234 ok (non-closed singly connected) frisbee 143566 155426 ok (closed singly connected) Checking geometry... Overall domain bounding box (-1.5 -2 -2) (4.5 2 2) Mesh (non-empty, non-wedge) directions (1 1 1) Mesh (non-empty) directions (1 1 1) Boundary openness (-7.91179e-17 1.45473e-17 -1.85313e-18) OK. Max cell openness = 4.42999e-16 OK. Max aspect ratio = 7.52203 OK. Minimum face area = 7.81257e-08. Maximum face area = 0.250442. Face area magnitudes OK. Min volume = 3.23513e-11. Max volume = 0.125241. Total volume = 95.9998. Cell volumes OK. Mesh non-orthogonality Max: 64.9948 average: 6.55222 Non-orthogonality check OK. Face pyramids OK. Max skewness = 1.24885 OK. Coupled point location match (average 0) OK. Mesh OK. Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; tolerance 1e-07; relTol 0.01; 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-08; relTol 0.1; nSweeps 1; } } SIMPLE { nNonOrthogonalCorrectors 0; residualControl { p 1e-8; U 1e-8; "(k|epsilon|omega)" 1e-8; } } relaxationFactors { fields { p 0.3; } equations { U 0.7; k 0.7; epsilon 0.7; R 0.7; nuTilda 0.7; } } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) bounded Gauss upwind; div(phi,k) bounded Gauss upwind; div(phi,epsilon) bounded Gauss upwind; div(phi,R) bounded Gauss upwind; div(R) Gauss linear; div(phi,nuTilda) 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 ; } // ************************************************************************* // Last edited by Yuby; March 11, 2015 at 11:39. Reason: complete |
|
March 11, 2015, 11:38 |
|
#2 |
Member
Rubén
Join Date: Oct 2014
Location: Munich
Posts: 47
Rep Power: 11 |
completeAnd, of course, my 0.org with properties calculated with http://www.cfd-online.com/Tools/turbulence.php for a eddy viscosity ratio of 0.1
k: Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object k; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #include "include/initialConditions" dimensions [0 2 -2 0 0 0 0]; internalField uniform $turbulentKE; boundaryField { inlet { type fixedValue; value uniform $turbulentKE; } outlet { type zeroGradient; } frisbee { type kqRWallFunction; value $turbulentKE; } #include "include/frontBackUpperPatches" } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object epsilon; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #include "include/initialConditions" dimensions [0 2 -3 0 0 0 0]; internalField uniform $turbulentEpsilon; boundaryField { inlet { type fixedValue; value uniform $turbulentEpsilon; } outlet { type zeroGradient; } frisbee { type epsilonWallFunction; value uniform $turbulentEpsilon; } #include "include/frontBackUpperPatches" } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; location "0"; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #include "include/initialConditions" dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { #include "include/fixedInlet" outlet { type zeroGradient; } frisbee { type rotatingWallVelocity; origin (0 0 0); axis (0 0 -1); omega $omegadisk; } #include "include/frontBackUpperPatches" } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #include "include/initialConditions" dimensions [0 2 -2 0 0 0 0]; internalField uniform $pressure; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } frisbee { type zeroGradient; } #include "include/frontBackUpperPatches" } // ************************************************************************* // Code:
| ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ inlet { type fixedValue; value uniform $flowVelocity; } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ top { type slip; } bottom { type slip; } frontandback { type slip; } // ************************************************************************* // Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ flowVelocity (16 0 0); omegadisk 37.7; pressure 0; turbulentKE 0.0096; turbulentEpsilon 5.5296; #inputMode merge // ************************************************************************* // Last edited by Yuby; March 11, 2015 at 11:41. Reason: complete |
|
March 18, 2015, 05:00 |
|
#3 |
Senior Member
Jan
Join Date: Jul 2009
Location: Hamburg
Posts: 137
Rep Power: 19 |
Try to use a MRF zone around your frisbee.
Regrads, Jan |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
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 |
Convergence Problems SimpleFOAM | Kutti | OpenFOAM | 16 | June 14, 2010 08:12 |
Getting faster convergence in simpleFoam | basneb | OpenFOAM | 8 | February 9, 2010 04:20 |
Definition of convergence criterion in simpleFoam | titio | OpenFOAM Running, Solving & CFD | 1 | February 6, 2010 01:34 |