
[Sponsors] 
April 12, 2011, 06:50 
Extrusion with OpenFoam problem No. Iterations 0

#1 
New Member

Hi, I'm trying to simulate a laminar flow in a cylinder (hole of the extrusion die) with inlet and outlet (bases) and wall (lateral). I modeled it with Salome 5 (CAE linux) and exported the mesh in UNV.
I post the files: RaspProperties /** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.7.0   \\ / A nd  Web: www.OpenFOAM.com   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object RASProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // RASModel laminar; turbulence off; printCoeffs on; // ************************************************** *********************** // transportProperties /** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.7.0   \\ / A nd  Web: www.OpenFOAM.com   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object transportProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // transportModel powerLaw; nu nu [ 0 2 1 0 0 0 0 ] 1e05; powerLawCoeffs { k k [0 2 1 0 0 0 0] 2800; n n [0 0 0 0 0 0 0] 0.4; nuMin nuMin [0 2 1 0 0 0 0] .1; nuMax nuMax [0 2 1 0 0 0 0] 1000000; } CrossPowerLawCoeffs { nu0 nu0 [ 0 2 1 0 0 0 0 ] 1e06; nuInf nuInf [ 0 2 1 0 0 0 0 ] 1e06; m m [ 0 0 1 0 0 0 0 ] 1; n n [ 0 0 0 0 0 0 0 ] 1; } BirdCarreauCoeffs { nu0 nu0 [ 0 2 1 0 0 0 0 ] 1e06; nuInf nuInf [ 0 2 1 0 0 0 0 ] 1e06; k k [ 0 0 1 0 0 0 0 ] 0; n n [ 0 0 0 0 0 0 0 ] 1; } // ************************************************** *********************** // p /** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.7.0   \\ / A nd  Web: www.OpenFOAM.com   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } wall { type zeroGradient; } } // ************************************************** *********************** // U /** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.7.0   \\ / A nd  Web: www.OpenFOAM.com   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type fixedValue; value uniform (.016 0 0); } outlet { type zeroGradient; } wall { type fixedValue; value uniform (0 0 0); } } // ************************************************** *********************** // fvSchemes /** C++ **\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 1.7.0   \\ / A nd  Web: www.OpenFOAM.com   \\/ M anipulation   \**/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default Euler; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) Gauss linear; } divSchemes { default none; div(phi,U) Gauss upwind; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div(phi,R) Gauss upwind; div(R) Gauss linear; div(phi,nuTilda) Gauss upwind; div((nuEff*dev(grad(U).T()))) Gauss linear; } laplacianSchemes { default none; laplacian(nuEff,U) Gauss linear corrected; laplacian((1A(U)),p) Gauss linear corrected; laplacian(DkEff,k) Gauss linear corrected; laplacian(DepsilonEff,epsilon) Gauss linear corrected; laplacian(DREff,R) Gauss linear corrected; laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; } interpolationSchemes { default linear; interpolate(U) linear; } snGradSchemes { default corrected; } fluxRequired { default no; p ; } // ************************************************** *********************** // then I run simpleFoam in the directory of the simulation that I created and this is the results: [...............] Time = 979 DILUPBiCG: Solving for Ux, Initial residual = 1.23786e16, Final residual = 1.23786e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.12365e16, Final residual = 1.12365e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.22862e16, Final residual = 1.22862e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30448e08 ExecutionTime = 4.93 s ClockTime = 9 s Time = 980 DILUPBiCG: Solving for Ux, Initial residual = 1.26725e16, Final residual = 1.26725e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.12963e16, Final residual = 1.12963e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.21408e16, Final residual = 1.21408e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30457e08 ExecutionTime = 4.94 s ClockTime = 9 s Time = 981 DILUPBiCG: Solving for Ux, Initial residual = 1.23068e16, Final residual = 1.23068e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.11831e16, Final residual = 1.11831e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.23645e16, Final residual = 1.23645e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30465e08 ExecutionTime = 4.94 s ClockTime = 9 s Time = 982 DILUPBiCG: Solving for Ux, Initial residual = 1.26931e16, Final residual = 1.26931e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.0935e16, Final residual = 1.0935e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.22962e16, Final residual = 1.22962e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30474e08 ExecutionTime = 4.95 s ClockTime = 9 s Time = 983 DILUPBiCG: Solving for Ux, Initial residual = 1.23579e16, Final residual = 1.23579e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.15394e16, Final residual = 1.15394e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.2132e16, Final residual = 1.2132e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30483e08 ExecutionTime = 4.96 s ClockTime = 9 s Time = 984 DILUPBiCG: Solving for Ux, Initial residual = 1.26194e16, Final residual = 1.26194e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.08941e16, Final residual = 1.08941e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.23718e16, Final residual = 1.23718e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30491e08 ExecutionTime = 4.96 s ClockTime = 9 s Time = 985 DILUPBiCG: Solving for Ux, Initial residual = 1.23803e16, Final residual = 1.23803e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.12121e16, Final residual = 1.12121e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.22867e16, Final residual = 1.22867e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.305e08 ExecutionTime = 4.96 s ClockTime = 9 s Time = 986 DILUPBiCG: Solving for Ux, Initial residual = 1.26728e16, Final residual = 1.26728e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.12798e16, Final residual = 1.12798e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.21437e16, Final residual = 1.21437e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30508e08 ExecutionTime = 4.96 s ClockTime = 9 s Time = 987 DILUPBiCG: Solving for Ux, Initial residual = 1.23051e16, Final residual = 1.23051e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.11958e16, Final residual = 1.11958e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.23644e16, Final residual = 1.23644e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30517e08 ExecutionTime = 4.97 s ClockTime = 9 s Time = 988 DILUPBiCG: Solving for Ux, Initial residual = 1.26928e16, Final residual = 1.26928e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.09385e16, Final residual = 1.09385e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.22949e16, Final residual = 1.22949e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30525e08 ExecutionTime = 4.98 s ClockTime = 9 s Time = 989 DILUPBiCG: Solving for Ux, Initial residual = 1.23595e16, Final residual = 1.23595e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.15446e16, Final residual = 1.15446e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.21323e16, Final residual = 1.21323e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30534e08 ExecutionTime = 4.98 s ClockTime = 9 s Time = 990 DILUPBiCG: Solving for Ux, Initial residual = 1.26196e16, Final residual = 1.26196e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.08924e16, Final residual = 1.08924e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.23746e16, Final residual = 1.23746e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30542e08 ExecutionTime = 4.98 s ClockTime = 9 s Time = 991 DILUPBiCG: Solving for Ux, Initial residual = 1.23788e16, Final residual = 1.23788e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.12266e16, Final residual = 1.12266e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.22859e16, Final residual = 1.22859e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30551e08 ExecutionTime = 4.98 s ClockTime = 9 s Time = 992 DILUPBiCG: Solving for Ux, Initial residual = 1.26728e16, Final residual = 1.26728e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.1293e16, Final residual = 1.1293e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.21393e16, Final residual = 1.21393e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30559e08 ExecutionTime = 4.99 s ClockTime = 9 s Time = 993 DILUPBiCG: Solving for Ux, Initial residual = 1.23065e16, Final residual = 1.23065e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.11725e16, Final residual = 1.11725e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.23649e16, Final residual = 1.23649e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30568e08 ExecutionTime = 5 s ClockTime = 9 s Time = 994 DILUPBiCG: Solving for Ux, Initial residual = 1.26933e16, Final residual = 1.26933e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.09154e16, Final residual = 1.09154e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.22979e16, Final residual = 1.22979e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30576e08 ExecutionTime = 5 s ClockTime = 9 s Time = 995 DILUPBiCG: Solving for Ux, Initial residual = 1.23577e16, Final residual = 1.23577e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.15583e16, Final residual = 1.15583e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.21322e16, Final residual = 1.21322e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30585e08 ExecutionTime = 5.01 s ClockTime = 9 s Time = 996 DILUPBiCG: Solving for Ux, Initial residual = 1.26195e16, Final residual = 1.26195e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.09024e16, Final residual = 1.09024e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.2373e16, Final residual = 1.2373e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30593e08 ExecutionTime = 5.02 s ClockTime = 9 s Time = 997 DILUPBiCG: Solving for Ux, Initial residual = 1.23799e16, Final residual = 1.23799e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.12231e16, Final residual = 1.12231e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.22866e16, Final residual = 1.22866e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30602e08 ExecutionTime = 5.02 s ClockTime = 9 s Time = 998 DILUPBiCG: Solving for Ux, Initial residual = 1.2673e16, Final residual = 1.2673e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.12967e16, Final residual = 1.12967e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.21424e16, Final residual = 1.21424e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.3061e08 ExecutionTime = 5.03 s ClockTime = 9 s Time = 999 DILUPBiCG: Solving for Ux, Initial residual = 1.23049e16, Final residual = 1.23049e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.11804e16, Final residual = 1.11804e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.2364e16, Final residual = 1.2364e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30619e08 ExecutionTime = 5.03 s ClockTime = 9 s Time = 1000 DILUPBiCG: Solving for Ux, Initial residual = 1.26929e16, Final residual = 1.26929e16, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 1.09316e16, Final residual = 1.09316e16, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1.22938e16, Final residual = 1.22938e16, No Iterations 0 DICPCG: Solving for p, Initial residual = 7.10497e07, Final residual = 7.10497e07, No Iterations 0 time step continuity errors : sum local = 3.29168e12, global = 8.50298e13, cumulative = 1.30627e08 ExecutionTime = 5.05 s ClockTime = 10 s End I can't find the problem, the issue. Someone can help me? Thank U! 

April 12, 2011, 07:14 

#2 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 253
Rep Power: 13 
Hi,
your simulation seems to be fully converged... how does it look like in Paraview? Martin 

April 12, 2011, 07:19 

#3  
New Member

Quote:
I don't understand this issue. absolutely 

April 12, 2011, 07:27 

#4 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 253
Rep Power: 13 
If you decrease the tolerances for the linear equation solver in system/fvSolution even further to let's say "tolerance 1e12;", then the initial residual is higher then desired and more iterations will be done...
But it shouldn't be necessary... your simulation results should be fine enough after a few hundred iterations. You can add "convergence 1e5;" to the SIMPLE dictionary in the fvSolution. simpleFoam will stop as soon as convergence criteria is reached: Code:
SIMPLE { nNonOrthogonalCorrectors 0; convergence 1.0e5; } 

April 12, 2011, 07:52 

#5  
New Member

Quote:
so the problem is in the solid model? 

April 12, 2011, 07:58 

#6 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 253
Rep Power: 13 
If simpleFoam makes 0 iterations for solving the linear equations system in the very first time steps, too, than there is a problem in your case. On the other hand if simpleFoam do make more than 0 iterations in the first time steps, everything is fine: simpleFoam found the correct solution and further iterations are not necessary.
You can send me your case, if you like, so I can have a look at it... my eMail comes via boardmail to prevent spamming... Martin 

April 12, 2011, 11:24 

#7  
New Member

Quote:
I want to give him REP+ 

March 27, 2016, 11:26 
what was the outcome

#8 
New Member
Join Date: Feb 2016
Posts: 5
Rep Power: 2 
Thank you for posting all of your code.
Could you please post the changes that you made to make the model run correctly. I would really appreciate it 

March 28, 2016, 11:08 

#9  
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,800
Blog Entries: 39
Rep Power: 106 
Quick answer:
Quote:
Nonetheless, all of the critical details and solutions are fully explained by MartinB. If you don't understand the explanations, then please provide information about your own case, namely the problem you're facing.
__________________


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
SimpleFoam k and epsilon bounded  nedved  OpenFOAM Running, Solving & CFD  15  December 2, 2016 03:40 
Moving mesh  Niklas Wikstrom (Wikstrom)  OpenFOAM Running, Solving & CFD  122  June 15, 2014 06:20 
icoLagrangianFoam OF1.6 myNewParticleSolver  heavy_user  OpenFOAM  16  February 11, 2012 06:15 
Low Mach number Compressible jet flow using LES  ankgupta8um  OpenFOAM Running, Solving & CFD  7  January 15, 2011 14:38 
lift and drag on ship superstructures  vaina74  OpenFOAM Running, Solving & CFD  3  June 8, 2010 12:30 