
[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: 252
Rep Power: 12 
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: 252
Rep Power: 12 
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: 252
Rep Power: 12 
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+ 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Moving mesh  Niklas Wikstrom (Wikstrom)  OpenFOAM Running, Solving & CFD  122  June 15, 2014 06:20 
SimpleFoam k and epsilon bounded  nedved  OpenFOAM Running, Solving & CFD  13  November 4, 2013 15:13 
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 