hfs |
November 14, 2012 07:46 |
pimpleFoam: turbulence->correct(); is not executed when using residualControl
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;
}
}
When I added that, pimpleFoam did not solve turbulance:
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();
}
[Q]: Is this a bug or am I misusing pimpleFoam ??
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
|