CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

pimpleFoam: turbulence->correct(); is not executed when using residualControl

Register Blogs Community New Posts Updated Threads Search

Like Tree4Likes
  • 4 Post By novakm

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 14, 2012, 08:46
Default pimpleFoam: turbulence->correct(); is not executed when using residualControl
  #1
hfs
Member
 
Join Date: Jul 2012
Posts: 66
Rep Power: 14
hfs is on a distinguished road
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
hfs is offline   Reply With Quote

Old   January 28, 2013, 11:40
Default
  #2
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 27
immortality is on a distinguished road
what does residualControl do when tolerances can be set at their own places?
immortality is offline   Reply With Quote

Old   January 28, 2013, 12:45
Default
  #3
hfs
Member
 
Join Date: Jul 2012
Posts: 66
Rep Power: 14
hfs is on a distinguished road
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..
hfs is offline   Reply With Quote

Old   October 29, 2013, 09:35
Default
  #4
Member
 
Martin Novák
Join Date: Dec 2012
Location: Prague
Posts: 70
Rep Power: 14
novakm is on a distinguished road
Quote:
Originally Posted by hfs View Post
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
Hi.

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
BigPapi34, kkpal, Kabra and 1 others like this.
novakm is offline   Reply With Quote

Reply

« Previous Thread | Next Thread »

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts
BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Turbulence postprocessing Mohsin FLUENT 2 October 3, 2016 15:18
Turbulence length scale and integral length scale rizhang CFX 2 April 22, 2016 08:22
Question on Turbulence Intensity Eric FLUENT 1 March 7, 2012 05:30
Discussion: Reason of Turbulence!! Wen Long Main CFD Forum 3 May 15, 2009 10:52


All times are GMT -4. The time now is 13:06.

Contact Us - CFD Online - Privacy Statement - Top