
[Sponsors] 
pimpleFoam: turbulence>correct(); is not executed when using residualControl 

LinkBack  Thread Tools  Display Modes 
November 14, 2012, 08:46 
pimpleFoam: turbulence>correct(); is not executed when using residualControl

#1 
Member
Join Date: Jul 2012
Posts: 60
Rep Power: 5 
Hello Foamers,
I am trying to use pimpleFoam and I used residualControl in the PIMPLE settings in fvSolution file: Code:
residualControl { p { tolerance 1e4; relTol 0; absTol 0; } } Code:
/**\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 2.1.x   \\ / A nd  Web: www.OpenFOAM.org   \\/ M anipulation   \**/ Build : 2.1.xe48b747827ac Exec : pimpleFoam Date : Nov 14 2012 Time : 13:15:58 Host : "stud3" PID : 14409 Case : /home/hosam/run/hosam/Real_Geo/02_10_12/8_5ButBigDomain/maxNonOrtho_180/Simulation_Case/pisoFoam/test/Mesh26_PIMPLE_residualControl nProcs : 1 sigFpe : Floating point exception trapping  not supported on this platform fileModificationChecking : Monitoring runtime modified files using timeStampMaster allowSystemOperations : Disallowing usersupplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type RASModel Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; sigmaEps 1.3; } No field sources present PIMPLE: max iterations = 5 field p : relTol 0, tolerance 0.0001 Starting time loop Courant Number mean: 0 max: 0 deltaT = 0.012 Time = 0.012 PIMPLE: iteration 1 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 3.997576e10, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 3.648792e10, No Iterations 2 GAMG: Solving for p, Initial residual = 1, Final residual = 4.4629999e08, No Iterations 16 GAMG: Solving for p, Initial residual = 0.0038485342, Final residual = 7.4367519e08, No Iterations 7 GAMG: Solving for p, Initial residual = 0.00014119928, Final residual = 4.6734605e08, No Iterations 5 GAMG: Solving for p, Initial residual = 7.969735e06, Final residual = 6.1914261e08, No Iterations 2 time step continuity errors : sum local = 2.9519967e09, global = 8.4429551e10, cumulative = 8.4429551e10 GAMG: Solving for p, Initial residual = 0.00063887725, Final residual = 8.0013507e08, No Iterations 8 GAMG: Solving for p, Initial residual = 0.00023360957, Final residual = 4.0079707e08, No Iterations 5 GAMG: Solving for p, Initial residual = 9.8986018e06, Final residual = 4.6317862e08, No Iterations 3 GAMG: Solving for p, Initial residual = 9.9131381e07, Final residual = 4.4517288e08, No Iterations 1 time step continuity errors : sum local = 2.082721e09, global = 6.2434109e10, cumulative = 2.1995442e10 GAMG: Solving for p, Initial residual = 0.00011127464, Final residual = 6.8708458e08, No Iterations 6 GAMG: Solving for p, Initial residual = 5.2130927e05, Final residual = 4.1693241e08, No Iterations 4 GAMG: Solving for p, Initial residual = 3.0705445e06, Final residual = 4.5800285e08, No Iterations 2 GAMG: Solving for p, Initial residual = 3.3095322e07, Final residual = 2.3148007e08, No Iterations 1 time step continuity errors : sum local = 1.0827011e09, global = 4.5240398e10, cumulative = 2.3244956e10 PIMPLE: converged in 1 iterations ExecutionTime = 2.07 s ClockTime = 2 s Courant Number mean: 0.019415651 max: 1.4022251 deltaT = 0.0042789137 Time = 0.0162789 PIMPLE: iteration 1 DILUPBiCG: Solving for Ux, Initial residual = 0.4070529, Final residual = 2.8920238e07, No Iterations 3 DILUPBiCG: Solving for Uy, Initial residual = 0.14303894, Final residual = 2.9725694e07, No Iterations 4 DILUPBiCG: Solving for Uz, Initial residual = 0.14520297, Final residual = 1.5216663e07, No Iterations 4 GAMG: Solving for p, Initial residual = 0.013748865, Final residual = 8.4076433e08, No Iterations 12 GAMG: Solving for p, Initial residual = 0.1532367, Final residual = 3.190578e08, No Iterations 11 GAMG: Solving for p, Initial residual = 0.0057219265, Final residual = 6.3769806e08, No Iterations 8 GAMG: Solving for p, Initial residual = 0.00031996465, Final residual = 6.2208569e08, No Iterations 5 time step continuity errors : sum local = 1.0636136e11, global = 4.7436717e12, cumulative = 2.3719323e10 GAMG: Solving for p, Initial residual = 0.018959309, Final residual = 6.6668286e08, No Iterations 11 GAMG: Solving for p, Initial residual = 0.013452293, Final residual = 5.2039685e08, No Iterations 9 GAMG: Solving for p, Initial residual = 0.00063041257, Final residual = 5.9604576e08, No Iterations 6 GAMG: Solving for p, Initial residual = 6.4189185e05, Final residual = 4.5021156e08, No Iterations 4 time step continuity errors : sum local = 3.8996509e12, global = 1.7599916e12, cumulative = 2.3543324e10 GAMG: Solving for p, Initial residual = 0.0056431921, Final residual = 8.3164206e08, No Iterations 9 GAMG: Solving for p, Initial residual = 0.002577829, Final residual = 6.7888934e08, No Iterations 7 GAMG: Solving for p, Initial residual = 0.00015773646, Final residual = 7.1871214e08, No Iterations 5 GAMG: Solving for p, Initial residual = 1.6937761e05, Final residual = 4.2968704e08, No Iterations 3 time step continuity errors : sum local = 3.6610407e12, global = 1.3752898e12, cumulative = 2.3405795e10 PIMPLE: converged in 1 iterations ExecutionTime = 4.05 s ClockTime = 4 s Courant Number mean: 0.0069410895 max: 0.4682004 deltaT = 0.0045695323 Time = 0.0208484 Without those residualControl the turbulence is calculated as it is supposed to (Solving for epsilon... Solving for k...): Code:
/**\  =========    \\ / F ield  OpenFOAM: The Open Source CFD Toolbox   \\ / O peration  Version: 2.1.x   \\ / A nd  Web: www.OpenFOAM.org   \\/ M anipulation   \**/ Build : 2.1.xe48b747827ac Exec : pimpleFoam Date : Nov 14 2012 Time : 13:14:00 Host : "stud3" PID : 14346 Case : /home/hosam/run/hosam/Real_Geo/02_10_12/8_5ButBigDomain/maxNonOrtho_180/Simulation_Case/pisoFoam/test/Mesh26 nProcs : 1 sigFpe : Floating point exception trapping  not supported on this platform fileModificationChecking : Monitoring runtime modified files using timeStampMaster allowSystemOperations : Disallowing usersupplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type RASModel Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; sigmaEps 1.3; } No field sources present PIMPLE: no residual control data found. Calculations will employ 5 corrector loops Starting time loop Courant Number mean: 0 max: 0 deltaT = 0.012 Time = 0.012 PIMPLE: iteration 1 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 3.997576e10, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 3.648792e10, No Iterations 2 GAMG: Solving for p, Initial residual = 1, Final residual = 4.4629999e08, No Iterations 16 GAMG: Solving for p, Initial residual = 0.0038485342, Final residual = 7.4367519e08, No Iterations 7 GAMG: Solving for p, Initial residual = 0.00014119928, Final residual = 4.6734605e08, No Iterations 5 GAMG: Solving for p, Initial residual = 7.969735e06, Final residual = 6.1914261e08, No Iterations 2 time step continuity errors : sum local = 2.9519967e09, global = 8.4429551e10, cumulative = 8.4429551e10 GAMG: Solving for p, Initial residual = 0.00063887725, Final residual = 8.0013507e08, No Iterations 8 GAMG: Solving for p, Initial residual = 0.00023360957, Final residual = 4.0079707e08, No Iterations 5 GAMG: Solving for p, Initial residual = 9.8986018e06, Final residual = 4.6317862e08, No Iterations 3 GAMG: Solving for p, Initial residual = 9.9131381e07, Final residual = 4.4517288e08, No Iterations 1 time step continuity errors : sum local = 2.082721e09, global = 6.2434109e10, cumulative = 2.1995442e10 GAMG: Solving for p, Initial residual = 0.00011127464, Final residual = 6.8708458e08, No Iterations 6 GAMG: Solving for p, Initial residual = 5.2130927e05, Final residual = 4.1693241e08, No Iterations 4 GAMG: Solving for p, Initial residual = 3.0705445e06, Final residual = 4.5800285e08, No Iterations 2 GAMG: Solving for p, Initial residual = 3.3095322e07, Final residual = 2.3148007e08, No Iterations 1 time step continuity errors : sum local = 1.0827011e09, global = 4.5240398e10, cumulative = 2.3244956e10 PIMPLE: iteration 2 DILUPBiCG: Solving for Ux, Initial residual = 0.0017391822, Final residual = 1.9077619e07, No Iterations 5 DILUPBiCG: Solving for Uy, Initial residual = 0.033640316, Final residual = 4.5074921e07, No Iterations 6 DILUPBiCG: Solving for Uz, Initial residual = 0.013436332, Final residual = 2.4356566e07, No Iterations 6 GAMG: Solving for p, Initial residual = 0.00039622137, Final residual = 5.4806417e08, No Iterations 7 GAMG: Solving for p, Initial residual = 0.00011152834, Final residual = 9.7282646e08, No Iterations 4 GAMG: Solving for p, Initial residual = 7.3440028e06, Final residual = 3.857655e08, No Iterations 3 GAMG: Solving for p, Initial residual = 8.3912385e07, Final residual = 3.2069867e08, No Iterations 1 time step continuity errors : sum local = 1.4729565e09, global = 6.910362e10, cumulative = 9.2348576e10 GAMG: Solving for p, Initial residual = 8.9523276e05, Final residual = 5.6033327e08, No Iterations 6 GAMG: Solving for p, Initial residual = 3.1353052e05, Final residual = 9.7884487e08, No Iterations 3 GAMG: Solving for p, Initial residual = 2.2950887e06, Final residual = 4.32794e08, No Iterations 2 GAMG: Solving for p, Initial residual = 2.8413122e07, Final residual = 2.3968221e08, No Iterations 1 time step continuity errors : sum local = 1.1008736e09, global = 3.6977169e10, cumulative = 5.5371407e10 GAMG: Solving for p, Initial residual = 1.4179402e05, Final residual = 5.4656075e08, No Iterations 4 GAMG: Solving for p, Initial residual = 3.6482307e06, Final residual = 3.5480209e08, No Iterations 2 GAMG: Solving for p, Initial residual = 2.9214641e07, Final residual = 2.2161815e08, No Iterations 1 GAMG: Solving for p, Initial residual = 4.8976129e08, Final residual = 4.8976129e08, No Iterations 0 time step continuity errors : sum local = 2.2494978e09, global = 6.0718691e10, cumulative = 5.3472843e11 PIMPLE: iteration 3 DILUPBiCG: Solving for Ux, Initial residual = 0.00099904702, Final residual = 9.797703e07, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0.0061705212, Final residual = 3.6144404e07, No Iterations 6 DILUPBiCG: Solving for Uz, Initial residual = 0.002854862, Final residual = 8.0323395e07, No Iterations 5 GAMG: Solving for p, Initial residual = 0.00010199419, Final residual = 5.4573314e08, No Iterations 5 GAMG: Solving for p, Initial residual = 1.3747455e05, Final residual = 7.7395363e08, No Iterations 2 GAMG: Solving for p, Initial residual = 1.1527826e06, Final residual = 5.8162218e08, No Iterations 1 GAMG: Solving for p, Initial residual = 1.7115834e07, Final residual = 2.4423164e08, No Iterations 1 time step continuity errors : sum local = 1.1218673e09, global = 4.7089403e10, cumulative = 5.2436687e10 GAMG: Solving for p, Initial residual = 2.532356e05, Final residual = 6.80238e08, No Iterations 4 GAMG: Solving for p, Initial residual = 4.9178153e06, Final residual = 3.7068718e08, No Iterations 2 GAMG: Solving for p, Initial residual = 4.8898736e07, Final residual = 2.7522039e08, No Iterations 1 GAMG: Solving for p, Initial residual = 7.9070835e08, Final residual = 7.9070835e08, No Iterations 0 time step continuity errors : sum local = 3.6320854e09, global = 2.520512e10, cumulative = 2.7231567e10 GAMG: Solving for p, Initial residual = 5.5005552e06, Final residual = 7.4555709e08, No Iterations 2 GAMG: Solving for p, Initial residual = 7.3381056e07, Final residual = 4.4533923e08, No Iterations 1 GAMG: Solving for p, Initial residual = 9.9862911e08, Final residual = 9.9862911e08, No Iterations 0 GAMG: Solving for p, Initial residual = 9.9862911e08, Final residual = 9.9862911e08, No Iterations 0 time step continuity errors : sum local = 4.5871552e09, global = 3.0422843e10, cumulative = 3.1912754e11 PIMPLE: iteration 4 DILUPBiCG: Solving for Ux, Initial residual = 0.00012784799, Final residual = 6.3909398e07, No Iterations 3 DILUPBiCG: Solving for Uy, Initial residual = 0.0020538227, Final residual = 1.190295e07, No Iterations 6 DILUPBiCG: Solving for Uz, Initial residual = 0.00070249921, Final residual = 5.8821456e07, No Iterations 4 GAMG: Solving for p, Initial residual = 2.3201161e05, Final residual = 6.9299398e08, No Iterations 3 GAMG: Solving for p, Initial residual = 2.5325036e06, Final residual = 7.9765817e08, No Iterations 1 GAMG: Solving for p, Initial residual = 3.0787708e07, Final residual = 3.1822449e08, No Iterations 1 GAMG: Solving for p, Initial residual = 7.0667078e08, Final residual = 7.0667078e08, No Iterations 0 time step continuity errors : sum local = 3.2461135e09, global = 4.2469636e10, cumulative = 4.5660911e10 GAMG: Solving for p, Initial residual = 4.954581e06, Final residual = 8.4111774e08, No Iterations 2 GAMG: Solving for p, Initial residual = 8.7467545e07, Final residual = 4.6929513e08, No Iterations 1 GAMG: Solving for p, Initial residual = 1.1672592e07, Final residual = 1.8701721e08, No Iterations 1 GAMG: Solving for p, Initial residual = 3.3828785e08, Final residual = 3.3828785e08, No Iterations 0 time step continuity errors : sum local = 1.5539322e09, global = 2.4755206e10, cumulative = 7.0416117e10 GAMG: Solving for p, Initial residual = 1.1347915e06, Final residual = 5.5625608e08, No Iterations 1 GAMG: Solving for p, Initial residual = 1.4195497e07, Final residual = 1.9716658e08, No Iterations 1 GAMG: Solving for p, Initial residual = 3.4739292e08, Final residual = 3.4739292e08, No Iterations 0 GAMG: Solving for p, Initial residual = 3.4739292e08, Final residual = 3.4739292e08, No Iterations 0 time step continuity errors : sum local = 1.5957555e09, global = 1.6299208e10, cumulative = 8.6715325e10 PIMPLE: iteration 5 DILUPBiCG: Solving for Ux, Initial residual = 1.9142307e05, Final residual = 4.6805704e09, No Iterations 5 DILUPBiCG: Solving for Uy, Initial residual = 0.00057598292, Final residual = 1.8627159e08, No Iterations 6 DILUPBiCG: Solving for Uz, Initial residual = 0.00015388608, Final residual = 4.6115833e08, No Iterations 5 GAMG: Solving for p, Initial residual = 5.3084093e06, Final residual = 4.0240895e08, No Iterations 2 GAMG: Solving for p, Initial residual = 7.2502905e07, Final residual = 3.1816825e08, No Iterations 1 GAMG: Solving for p, Initial residual = 1.1078528e07, Final residual = 1.1327836e08, No Iterations 1 GAMG: Solving for p, Initial residual = 2.5294866e08, Final residual = 2.5294866e08, No Iterations 0 time step continuity errors : sum local = 1.1619251e09, global = 5.1654033e11, cumulative = 9.1880728e10 GAMG: Solving for p, Initial residual = 1.3990436e06, Final residual = 6.8777336e08, No Iterations 1 GAMG: Solving for p, Initial residual = 2.6276018e07, Final residual = 1.8556365e08, No Iterations 1 GAMG: Solving for p, Initial residual = 4.4262446e08, Final residual = 4.4262446e08, No Iterations 0 GAMG: Solving for p, Initial residual = 4.4262446e08, Final residual = 4.4262446e08, No Iterations 0 time step continuity errors : sum local = 2.033205e09, global = 2.8614673e11, cumulative = 9.4742196e10 GAMG: Solving for p, Initial residual = 3.1592008e07, Final residual = 1.7179692e08, No Iterations 1 GAMG: Solving for p, Initial residual = 4.4681906e08, Final residual = 4.4681906e08, No Iterations 0 GAMG: Solving for p, Initial residual = 4.4681906e08, Final residual = 4.4681906e08, No Iterations 0 GAMG: Solving for p, Initial residual = 4.4681906e08, Final residual = 5.8891984e09, No Iterations 1 time step continuity errors : sum local = 2.7052159e10, global = 1.8021546e11, cumulative = 9.654435e10 DILUPBiCG: Solving for epsilon, Initial residual = 0.00061820421, Final residual = 8.3080146e06, No Iterations 2 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.00082386122, No Iterations 4 ExecutionTime = 4.66 s ClockTime = 4 s Courant Number mean: 0.019482616 max: 1.8578018 deltaT = 0.0032296233 Time = 0.0152296 It is clear that the solver is not executing "turbulence>correct();" because pimple.turbCorr() is not 1 for some reason from pimpleFoam.C: Code:
if (pimple.turbCorr()) { turbulence>correct(); } You can take a look at the case with and without residualControl using the following two links: http://dl.dropbox.com/u/32336583/Mes...Control.tar.gz http://dl.dropbox.com/u/32336583/Mes...Control.tar.gz Thanks for your help 

January 28, 2013, 11:40 

#2 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,205
Rep Power: 17 
what does residualControl do when tolerances can be set at their own places?


January 28, 2013, 12:45 

#3 
Member
Join Date: Jul 2012
Posts: 60
Rep Power: 5 
I am not sure.. I am just used to it..
I am not using pimpleFoam anymore .. using pisoFoam with less timestep feels safer! Thanks for the reply.. 

October 29, 2013, 09:35 

#4  
Member
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 56
Rep Power: 4 
Quote:
In pimple control is by default set the flag turbOnFInallIterOnly to true. In order to get turb corr in every PIMPLE iter it is needed to specify turbOnFInallIterOnly false; in pimple control. Hope it helps, Martin 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Question on Turbulence Intensity  Eric  FLUENT  1  March 7, 2012 05:30 
Turbulence postprocessing  Mohsin  FLUENT  0  September 19, 2011 21:05 
Turbulence length scale and integral length scale  rizhang  CFX  1  September 10, 2009 06:38 
Discussion: Reason of Turbulence!!  Wen Long  Main CFD Forum  3  May 15, 2009 09:52 