CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   pimpleFoam: turbulence->correct(); is not executed when using residualControl (http://www.cfd-online.com/Forums/openfoam-solving/109288-pimplefoam-turbulence-correct-not-executed-when-using-residualcontrol.html)

hfs November 14, 2012 08: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

immortality January 28, 2013 11:40

what does residualControl do when tolerances can be set at their own places?

hfs January 28, 2013 12:45

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..

novakm October 29, 2013 09:35

Quote:

Originally Posted by hfs (Post 392044)
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


All times are GMT -4. The time now is 20:30.