|
[Sponsors] |
pimpleFoam: turbulence->correct(); is not executed when using residualControl |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 14, 2012, 07:46 |
pimpleFoam: turbulence->correct(); is not executed when using residualControl
|
#1 |
Member
Join Date: Jul 2012
Posts: 66
Rep Power: 14 |
Hello Foamers,
I am trying to use pimpleFoam and I used residualControl in the PIMPLE settings in fvSolution file: Code:
residualControl { p { tolerance 1e-4; 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.x-e48b747827ac 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 run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied 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.997576e-10, 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.648792e-10, No Iterations 2 GAMG: Solving for p, Initial residual = 1, Final residual = 4.4629999e-08, No Iterations 16 GAMG: Solving for p, Initial residual = 0.0038485342, Final residual = 7.4367519e-08, No Iterations 7 GAMG: Solving for p, Initial residual = 0.00014119928, Final residual = 4.6734605e-08, No Iterations 5 GAMG: Solving for p, Initial residual = 7.969735e-06, Final residual = 6.1914261e-08, No Iterations 2 time step continuity errors : sum local = 2.9519967e-09, global = -8.4429551e-10, cumulative = -8.4429551e-10 GAMG: Solving for p, Initial residual = 0.00063887725, Final residual = 8.0013507e-08, No Iterations 8 GAMG: Solving for p, Initial residual = 0.00023360957, Final residual = 4.0079707e-08, No Iterations 5 GAMG: Solving for p, Initial residual = 9.8986018e-06, Final residual = 4.6317862e-08, No Iterations 3 GAMG: Solving for p, Initial residual = 9.9131381e-07, Final residual = 4.4517288e-08, No Iterations 1 time step continuity errors : sum local = 2.082721e-09, global = 6.2434109e-10, cumulative = -2.1995442e-10 GAMG: Solving for p, Initial residual = 0.00011127464, Final residual = 6.8708458e-08, No Iterations 6 GAMG: Solving for p, Initial residual = 5.2130927e-05, Final residual = 4.1693241e-08, No Iterations 4 GAMG: Solving for p, Initial residual = 3.0705445e-06, Final residual = 4.5800285e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 3.3095322e-07, Final residual = 2.3148007e-08, No Iterations 1 time step continuity errors : sum local = 1.0827011e-09, global = 4.5240398e-10, cumulative = 2.3244956e-10 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.8920238e-07, No Iterations 3 DILUPBiCG: Solving for Uy, Initial residual = 0.14303894, Final residual = 2.9725694e-07, No Iterations 4 DILUPBiCG: Solving for Uz, Initial residual = 0.14520297, Final residual = 1.5216663e-07, No Iterations 4 GAMG: Solving for p, Initial residual = 0.013748865, Final residual = 8.4076433e-08, No Iterations 12 GAMG: Solving for p, Initial residual = 0.1532367, Final residual = 3.190578e-08, No Iterations 11 GAMG: Solving for p, Initial residual = 0.0057219265, Final residual = 6.3769806e-08, No Iterations 8 GAMG: Solving for p, Initial residual = 0.00031996465, Final residual = 6.2208569e-08, No Iterations 5 time step continuity errors : sum local = 1.0636136e-11, global = 4.7436717e-12, cumulative = 2.3719323e-10 GAMG: Solving for p, Initial residual = 0.018959309, Final residual = 6.6668286e-08, No Iterations 11 GAMG: Solving for p, Initial residual = 0.013452293, Final residual = 5.2039685e-08, No Iterations 9 GAMG: Solving for p, Initial residual = 0.00063041257, Final residual = 5.9604576e-08, No Iterations 6 GAMG: Solving for p, Initial residual = 6.4189185e-05, Final residual = 4.5021156e-08, No Iterations 4 time step continuity errors : sum local = 3.8996509e-12, global = -1.7599916e-12, cumulative = 2.3543324e-10 GAMG: Solving for p, Initial residual = 0.0056431921, Final residual = 8.3164206e-08, No Iterations 9 GAMG: Solving for p, Initial residual = 0.002577829, Final residual = 6.7888934e-08, No Iterations 7 GAMG: Solving for p, Initial residual = 0.00015773646, Final residual = 7.1871214e-08, No Iterations 5 GAMG: Solving for p, Initial residual = 1.6937761e-05, Final residual = 4.2968704e-08, No Iterations 3 time step continuity errors : sum local = 3.6610407e-12, global = -1.3752898e-12, cumulative = 2.3405795e-10 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.x-e48b747827ac 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 run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied 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.997576e-10, 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.648792e-10, No Iterations 2 GAMG: Solving for p, Initial residual = 1, Final residual = 4.4629999e-08, No Iterations 16 GAMG: Solving for p, Initial residual = 0.0038485342, Final residual = 7.4367519e-08, No Iterations 7 GAMG: Solving for p, Initial residual = 0.00014119928, Final residual = 4.6734605e-08, No Iterations 5 GAMG: Solving for p, Initial residual = 7.969735e-06, Final residual = 6.1914261e-08, No Iterations 2 time step continuity errors : sum local = 2.9519967e-09, global = -8.4429551e-10, cumulative = -8.4429551e-10 GAMG: Solving for p, Initial residual = 0.00063887725, Final residual = 8.0013507e-08, No Iterations 8 GAMG: Solving for p, Initial residual = 0.00023360957, Final residual = 4.0079707e-08, No Iterations 5 GAMG: Solving for p, Initial residual = 9.8986018e-06, Final residual = 4.6317862e-08, No Iterations 3 GAMG: Solving for p, Initial residual = 9.9131381e-07, Final residual = 4.4517288e-08, No Iterations 1 time step continuity errors : sum local = 2.082721e-09, global = 6.2434109e-10, cumulative = -2.1995442e-10 GAMG: Solving for p, Initial residual = 0.00011127464, Final residual = 6.8708458e-08, No Iterations 6 GAMG: Solving for p, Initial residual = 5.2130927e-05, Final residual = 4.1693241e-08, No Iterations 4 GAMG: Solving for p, Initial residual = 3.0705445e-06, Final residual = 4.5800285e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 3.3095322e-07, Final residual = 2.3148007e-08, No Iterations 1 time step continuity errors : sum local = 1.0827011e-09, global = 4.5240398e-10, cumulative = 2.3244956e-10 PIMPLE: iteration 2 DILUPBiCG: Solving for Ux, Initial residual = 0.0017391822, Final residual = 1.9077619e-07, No Iterations 5 DILUPBiCG: Solving for Uy, Initial residual = 0.033640316, Final residual = 4.5074921e-07, No Iterations 6 DILUPBiCG: Solving for Uz, Initial residual = 0.013436332, Final residual = 2.4356566e-07, No Iterations 6 GAMG: Solving for p, Initial residual = 0.00039622137, Final residual = 5.4806417e-08, No Iterations 7 GAMG: Solving for p, Initial residual = 0.00011152834, Final residual = 9.7282646e-08, No Iterations 4 GAMG: Solving for p, Initial residual = 7.3440028e-06, Final residual = 3.857655e-08, No Iterations 3 GAMG: Solving for p, Initial residual = 8.3912385e-07, Final residual = 3.2069867e-08, No Iterations 1 time step continuity errors : sum local = 1.4729565e-09, global = 6.910362e-10, cumulative = 9.2348576e-10 GAMG: Solving for p, Initial residual = 8.9523276e-05, Final residual = 5.6033327e-08, No Iterations 6 GAMG: Solving for p, Initial residual = 3.1353052e-05, Final residual = 9.7884487e-08, No Iterations 3 GAMG: Solving for p, Initial residual = 2.2950887e-06, Final residual = 4.32794e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 2.8413122e-07, Final residual = 2.3968221e-08, No Iterations 1 time step continuity errors : sum local = 1.1008736e-09, global = -3.6977169e-10, cumulative = 5.5371407e-10 GAMG: Solving for p, Initial residual = 1.4179402e-05, Final residual = 5.4656075e-08, No Iterations 4 GAMG: Solving for p, Initial residual = 3.6482307e-06, Final residual = 3.5480209e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 2.9214641e-07, Final residual = 2.2161815e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 4.8976129e-08, Final residual = 4.8976129e-08, No Iterations 0 time step continuity errors : sum local = 2.2494978e-09, global = -6.0718691e-10, cumulative = -5.3472843e-11 PIMPLE: iteration 3 DILUPBiCG: Solving for Ux, Initial residual = 0.00099904702, Final residual = 9.797703e-07, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0.0061705212, Final residual = 3.6144404e-07, No Iterations 6 DILUPBiCG: Solving for Uz, Initial residual = 0.002854862, Final residual = 8.0323395e-07, No Iterations 5 GAMG: Solving for p, Initial residual = 0.00010199419, Final residual = 5.4573314e-08, No Iterations 5 GAMG: Solving for p, Initial residual = 1.3747455e-05, Final residual = 7.7395363e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 1.1527826e-06, Final residual = 5.8162218e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 1.7115834e-07, Final residual = 2.4423164e-08, No Iterations 1 time step continuity errors : sum local = 1.1218673e-09, global = -4.7089403e-10, cumulative = -5.2436687e-10 GAMG: Solving for p, Initial residual = 2.532356e-05, Final residual = 6.80238e-08, No Iterations 4 GAMG: Solving for p, Initial residual = 4.9178153e-06, Final residual = 3.7068718e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 4.8898736e-07, Final residual = 2.7522039e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 7.9070835e-08, Final residual = 7.9070835e-08, No Iterations 0 time step continuity errors : sum local = 3.6320854e-09, global = 2.520512e-10, cumulative = -2.7231567e-10 GAMG: Solving for p, Initial residual = 5.5005552e-06, Final residual = 7.4555709e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 7.3381056e-07, Final residual = 4.4533923e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 9.9862911e-08, Final residual = 9.9862911e-08, No Iterations 0 GAMG: Solving for p, Initial residual = 9.9862911e-08, Final residual = 9.9862911e-08, No Iterations 0 time step continuity errors : sum local = 4.5871552e-09, global = 3.0422843e-10, cumulative = 3.1912754e-11 PIMPLE: iteration 4 DILUPBiCG: Solving for Ux, Initial residual = 0.00012784799, Final residual = 6.3909398e-07, No Iterations 3 DILUPBiCG: Solving for Uy, Initial residual = 0.0020538227, Final residual = 1.190295e-07, No Iterations 6 DILUPBiCG: Solving for Uz, Initial residual = 0.00070249921, Final residual = 5.8821456e-07, No Iterations 4 GAMG: Solving for p, Initial residual = 2.3201161e-05, Final residual = 6.9299398e-08, No Iterations 3 GAMG: Solving for p, Initial residual = 2.5325036e-06, Final residual = 7.9765817e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 3.0787708e-07, Final residual = 3.1822449e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 7.0667078e-08, Final residual = 7.0667078e-08, No Iterations 0 time step continuity errors : sum local = 3.2461135e-09, global = 4.2469636e-10, cumulative = 4.5660911e-10 GAMG: Solving for p, Initial residual = 4.954581e-06, Final residual = 8.4111774e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 8.7467545e-07, Final residual = 4.6929513e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 1.1672592e-07, Final residual = 1.8701721e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 3.3828785e-08, Final residual = 3.3828785e-08, No Iterations 0 time step continuity errors : sum local = 1.5539322e-09, global = 2.4755206e-10, cumulative = 7.0416117e-10 GAMG: Solving for p, Initial residual = 1.1347915e-06, Final residual = 5.5625608e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 1.4195497e-07, Final residual = 1.9716658e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 3.4739292e-08, Final residual = 3.4739292e-08, No Iterations 0 GAMG: Solving for p, Initial residual = 3.4739292e-08, Final residual = 3.4739292e-08, No Iterations 0 time step continuity errors : sum local = 1.5957555e-09, global = 1.6299208e-10, cumulative = 8.6715325e-10 PIMPLE: iteration 5 DILUPBiCG: Solving for Ux, Initial residual = 1.9142307e-05, Final residual = 4.6805704e-09, No Iterations 5 DILUPBiCG: Solving for Uy, Initial residual = 0.00057598292, Final residual = 1.8627159e-08, No Iterations 6 DILUPBiCG: Solving for Uz, Initial residual = 0.00015388608, Final residual = 4.6115833e-08, No Iterations 5 GAMG: Solving for p, Initial residual = 5.3084093e-06, Final residual = 4.0240895e-08, No Iterations 2 GAMG: Solving for p, Initial residual = 7.2502905e-07, Final residual = 3.1816825e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 1.1078528e-07, Final residual = 1.1327836e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 2.5294866e-08, Final residual = 2.5294866e-08, No Iterations 0 time step continuity errors : sum local = 1.1619251e-09, global = 5.1654033e-11, cumulative = 9.1880728e-10 GAMG: Solving for p, Initial residual = 1.3990436e-06, Final residual = 6.8777336e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 2.6276018e-07, Final residual = 1.8556365e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 4.4262446e-08, Final residual = 4.4262446e-08, No Iterations 0 GAMG: Solving for p, Initial residual = 4.4262446e-08, Final residual = 4.4262446e-08, No Iterations 0 time step continuity errors : sum local = 2.033205e-09, global = 2.8614673e-11, cumulative = 9.4742196e-10 GAMG: Solving for p, Initial residual = 3.1592008e-07, Final residual = 1.7179692e-08, No Iterations 1 GAMG: Solving for p, Initial residual = 4.4681906e-08, Final residual = 4.4681906e-08, No Iterations 0 GAMG: Solving for p, Initial residual = 4.4681906e-08, Final residual = 4.4681906e-08, No Iterations 0 GAMG: Solving for p, Initial residual = 4.4681906e-08, Final residual = 5.8891984e-09, No Iterations 1 time step continuity errors : sum local = 2.7052159e-10, global = 1.8021546e-11, cumulative = 9.654435e-10 DILUPBiCG: Solving for epsilon, Initial residual = 0.00061820421, Final residual = 8.3080146e-06, 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, 10:40 |
|
#2 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27 |
what does residualControl do when tolerances can be set at their own places?
|
|
January 28, 2013, 11:45 |
|
#3 |
Member
Join Date: Jul 2012
Posts: 66
Rep Power: 14 |
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, 08:35 |
|
#4 | |
Member
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 70
Rep Power: 13 |
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 | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Turbulence postprocessing | Mohsin | FLUENT | 2 | October 3, 2016 14:18 |
Turbulence length scale and integral length scale | rizhang | CFX | 2 | April 22, 2016 07:22 |
Question on Turbulence Intensity | Eric | FLUENT | 1 | March 7, 2012 04:30 |
Discussion: Reason of Turbulence!! | Wen Long | Main CFD Forum | 3 | May 15, 2009 09:52 |