Floating point exception error

lpz_michele · October 16, 2015, 07:08

I get iterations run and PIMPLE does not converge within the 20 iterations

Code:

Courant Number mean: 0.000445293 max: 6286.11
Time = 2.18166e-05

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.805531, Final residual = 0.000399167, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.259181, Final residual = 1.11701e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.259055, Final residual = 7.76482e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.172039, Final residual = 0.00111739, No Iterations 5
time step continuity errors : sum local = 6.72553e-07, global = 3.46612e-08, cumulative = -8.11236e-08
GAMG:  Solving for p, Initial residual = 0.155541, Final residual = 0.00144341, No Iterations 5
time step continuity errors : sum local = 6.73155e-07, global = 5.4537e-08, cumulative = -2.65865e-08
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 0.214015, Final residual = 0.00118249, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0435176, Final residual = 1.47477e-05, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.0435897, Final residual = 0.000215137, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.175285, Final residual = 0.00146206, No Iterations 6
time step continuity errors : sum local = 6.99109e-07, global = 5.46995e-08, cumulative = 2.8113e-08
GAMG:  Solving for p, Initial residual = 0.167401, Final residual = 0.00146482, No Iterations 7
time step continuity errors : sum local = 5.14118e-07, global = 4.024e-08, cumulative = 6.8353e-08
PIMPLE: iteration 3
DILUPBiCG:  Solving for Ux, Initial residual = 0.0397966, Final residual = 0.000267232, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0182205, Final residual = 0.000165177, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.0182424, Final residual = 3.58317e-05, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.149562, Final residual = 0.00140676, No Iterations 9
time step continuity errors : sum local = 4.84271e-07, global = 5.78833e-09, cumulative = 7.41413e-08
GAMG:  Solving for p, Initial residual = 0.144655, Final residual = 0.00136108, No Iterations 10
time step continuity errors : sum local = 3.4047e-07, global = 1.41064e-09, cumulative = 7.55519e-08
PIMPLE: iteration 4
DILUPBiCG:  Solving for Ux, Initial residual = 0.103625, Final residual = 0.000439838, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.028389, Final residual = 0.000195749, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.0284599, Final residual = 0.000233605, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.109704, Final residual = 0.00104964, No Iterations 13
time step continuity errors : sum local = 2.62813e-07, global = -6.01098e-11, cumulative = 7.54918e-08
GAMG:  Solving for p, Initial residual = 0.110735, Final residual = 0.00107394, No Iterations 16
time step continuity errors : sum local = 1.87831e-07, global = 5.53782e-11, cumulative = 7.55472e-08
PIMPLE: iteration 5
DILUPBiCG:  Solving for Ux, Initial residual = 0.0734305, Final residual = 0.000660661, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.0225332, Final residual = 0.000175618, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.0224398, Final residual = 8.654e-05, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0659701, Final residual = 0.00101109, No Iterations 20
time step continuity errors : sum local = 2.12618e-07, global = 7.06972e-11, cumulative = 7.56179e-08
GAMG:  Solving for p, Initial residual = 0.0350583, Final residual = 0.000250579, No Iterations 19
time step continuity errors : sum local = 7.24896e-08, global = -8.52339e-11, cumulative = 7.55327e-08
PIMPLE: iteration 6
DILUPBiCG:  Solving for Ux, Initial residual = 0.0128538, Final residual = 7.20424e-05, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.0128257, Final residual = 8.62813e-05, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.0205682, Final residual = 0.00014935, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0262762, Final residual = 0.00275159, No Iterations 20
time step continuity errors : sum local = 8.82026e-07, global = 3.06557e-10, cumulative = 7.58392e-08
GAMG:  Solving for p, Initial residual = 0.00603007, Final residual = 0.000797886, No Iterations 20
time step continuity errors : sum local = 8.11721e-07, global = 3.35122e-10, cumulative = 7.61743e-08
PIMPLE: iteration 7
DILUPBiCG:  Solving for Ux, Initial residual = 0.00818378, Final residual = 6.74005e-05, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.00530971, Final residual = 4.46342e-05, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.00671891, Final residual = 5.04277e-05, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.00924401, Final residual = 0.00113134, No Iterations 20
time step continuity errors : sum local = 7.99426e-07, global = 2.1622e-09, cumulative = 7.83365e-08
GAMG:  Solving for p, Initial residual = 0.0019598, Final residual = 0.000498121, No Iterations 20
time step continuity errors : sum local = 1.42565e-06, global = 1.25864e-09, cumulative = 7.95952e-08
PIMPLE: iteration 8
DILUPBiCG:  Solving for Ux, Initial residual = 0.0116867, Final residual = 5.25656e-05, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.00881643, Final residual = 9.31218e-06, No Iterations 6
DILUPBiCG:  Solving for Uz, Initial residual = 0.0125595, Final residual = 5.91165e-05, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.00276154, Final residual = 0.00141218, No Iterations 20
time step continuity errors : sum local = 4.46938e-06, global = 2.79248e-09, cumulative = 8.23876e-08
GAMG:  Solving for p, Initial residual = 0.000731615, Final residual = 0.000162034, No Iterations 20
time step continuity errors : sum local = 2.24711e-06, global = 2.33711e-09, cumulative = 8.47247e-08
PIMPLE: iteration 9
DILUPBiCG:  Solving for Ux, Initial residual = 0.0109847, Final residual = 9.57044e-05, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.00745915, Final residual = 3.29966e-05, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.0143351, Final residual = 0.00012274, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.00179665, Final residual = 0.000158175, No Iterations 20
time step continuity errors : sum local = 1.91208e-06, global = 6.72501e-09, cumulative = 9.14498e-08
GAMG:  Solving for p, Initial residual = 0.000412722, Final residual = 4.45279e-05, No Iterations 20
time step continuity errors : sum local = 2.64616e-06, global = 5.439e-09, cumulative = 9.68888e-08
PIMPLE: iteration 10
DILUPBiCG:  Solving for Ux, Initial residual = 0.0517594, Final residual = 0.000275136, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.0395561, Final residual = 0.000390823, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.0775718, Final residual = 0.000202579, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000520617, Final residual = 7.39019e-05, No Iterations 20
time step continuity errors : sum local = 5.80129e-06, global = 2.08963e-08, cumulative = 1.17785e-07
GAMG:  Solving for p, Initial residual = 8.30277e-05, Final residual = 2.53956e-05, No Iterations 20
time step continuity errors : sum local = 9.63303e-06, global = 5.84788e-09, cumulative = 1.23633e-07
PIMPLE: iteration 11
DILUPBiCG:  Solving for Ux, Initial residual = 0.0479551, Final residual = 0.000288326, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.0369187, Final residual = 3.01401e-05, No Iterations 7
DILUPBiCG:  Solving for Uz, Initial residual = 0.0405316, Final residual = 0.000322249, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000361656, Final residual = 2.35622e-05, No Iterations 20
time step continuity errors : sum local = 6.17042e-06, global = 2.85608e-08, cumulative = 1.52194e-07
GAMG:  Solving for p, Initial residual = 0.000107887, Final residual = 1.73687e-05, No Iterations 20
time step continuity errors : sum local = 2.02623e-05, global = 2.8365e-08, cumulative = 1.80559e-07
PIMPLE: iteration 12
DILUPBiCG:  Solving for Ux, Initial residual = 0.145947, Final residual = 0.000567913, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.134679, Final residual = 0.000446244, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.26212, Final residual = 0.000566838, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000139961, Final residual = 7.5937e-06, No Iterations 20
time step continuity errors : sum local = 1.07131e-05, global = 3.89594e-08, cumulative = 2.19518e-07
GAMG:  Solving for p, Initial residual = 1.89412e-05, Final residual = 2.06561e-06, No Iterations 20
time step continuity errors : sum local = 1.27439e-05, global = 3.77007e-08, cumulative = 2.57219e-07
PIMPLE: iteration 13
DILUPBiCG:  Solving for Ux, Initial residual = 0.136473, Final residual = 0.00042449, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0.100285, Final residual = 0.000653087, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.115866, Final residual = 0.00111478, No Iterations 4
GAMG:  Solving for p, Initial residual = 7.6738e-05, Final residual = 3.33884e-06, No Iterations 20
time step continuity errors : sum local = 1.67575e-05, global = 1.03064e-07, cumulative = 3.60283e-07
GAMG:  Solving for p, Initial residual = 1.75085e-05, Final residual = 1.79235e-06, No Iterations 20
time step continuity errors : sum local = 4.26782e-05, global = 1.01736e-07, cumulative = 4.62019e-07
PIMPLE: iteration 14
DILUPBiCG:  Solving for Ux, Initial residual = 0.355816, Final residual = 0.00191036, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.312617, Final residual = 0.00200449, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.426679, Final residual = 0.00270298, No Iterations 5
GAMG:  Solving for p, Initial residual = 2.80793e-05, Final residual = 3.36523e-06, No Iterations 20
time step continuity errors : sum local = 9.70931e-05, global = 3.03778e-07, cumulative = 7.65797e-07
GAMG:  Solving for p, Initial residual = 4.86248e-06, Final residual = 1.3074e-06, No Iterations 20
time step continuity errors : sum local = 0.000171439, global = 9.36991e-07, cumulative = 1.70279e-06
PIMPLE: iteration 15
DILUPBiCG:  Solving for Ux, Initial residual = 0.309695, Final residual = 0.000814464, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0.256879, Final residual = 0.00121518, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.29365, Final residual = 0.00172223, No Iterations 4
GAMG:  Solving for p, Initial residual = 7.89801e-06, Final residual = 5.69678e-07, No Iterations 20
time step continuity errors : sum local = 6.35809e-05, global = 3.37368e-07, cumulative = 2.04016e-06
GAMG:  Solving for p, Initial residual = 1.87393e-06, Final residual = 1.37319e-07, No Iterations 20
time step continuity errors : sum local = 7.52094e-05, global = 8.86364e-07, cumulative = 2.92652e-06
PIMPLE: iteration 16
DILUPBiCG:  Solving for Ux, Initial residual = 0.502919, Final residual = 0.00251847, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.468765, Final residual = 0.000525735, No Iterations 6
DILUPBiCG:  Solving for Uz, Initial residual = 0.508902, Final residual = 0.00159233, No Iterations 5
GAMG:  Solving for p, Initial residual = 4.88572e-06, Final residual = 2.37028e-07, No Iterations 20
time step continuity errors : sum local = 0.000163834, global = 3.16266e-06, cumulative = 6.08918e-06
GAMG:  Solving for p, Initial residual = 1.25873e-06, Final residual = 3.6224e-08, No Iterations 20
time step continuity errors : sum local = 0.000120213, global = 2.26185e-06, cumulative = 8.35103e-06
PIMPLE: iteration 17
DILUPBiCG:  Solving for Ux, Initial residual = 0.424699, Final residual = 0.00246043, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.420056, Final residual = 0.00195944, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.47174, Final residual = 0.00126653, No Iterations 5
GAMG:  Solving for p, Initial residual = 4.49696e-06, Final residual = 8.49181e-07, No Iterations 20
time step continuity errors : sum local = 0.00177044, global = -6.59931e-07, cumulative = 7.6911e-06
GAMG:  Solving for p, Initial residual = 9.26159e-07, Final residual = 3.15812e-07, No Iterations 20
time step continuity errors : sum local = 0.00282848, global = -1.38092e-07, cumulative = 7.55301e-06
PIMPLE: iteration 18
DILUPBiCG:  Solving for Ux, Initial residual = 0.501149, Final residual = 0.00197974, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.469467, Final residual = 0.00158737, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.521702, Final residual = 0.00233425, No Iterations 5
GAMG:  Solving for p, Initial residual = 1.63401e-06, Final residual = 4.26722e-07, No Iterations 20
time step continuity errors : sum local = 0.00438776, global = 1.68061e-05, cumulative = 2.43591e-05
GAMG:  Solving for p, Initial residual = 5.11766e-07, Final residual = 1.88297e-07, No Iterations 20
time step continuity errors : sum local = 0.00840377, global = 7.84767e-06, cumulative = 3.22068e-05
PIMPLE: iteration 19
DILUPBiCG:  Solving for Ux, Initial residual = 0.427517, Final residual = 0.00221743, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0.399281, Final residual = 0.00209523, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.475571, Final residual = 0.00229804, No Iterations 4
GAMG:  Solving for p, Initial residual = 8.19415e-07, Final residual = 4.47587e-07, No Iterations 20
time step continuity errors : sum local = 0.0170068, global = 5.09564e-05, cumulative = 8.31631e-05
GAMG:  Solving for p, Initial residual = 1.55778e-07, Final residual = 1.8941e-06, No Iterations 20
time step continuity errors : sum local = 0.347015, global = -2.01564e-06, cumulative = 8.11475e-05
PIMPLE: iteration 20
DILUPBiCG:  Solving for Ux, Initial residual = 0.851237, Final residual = 0.00741752, No Iterations 72
DILUPBiCG:  Solving for Uy, Initial residual = 0.95283, Final residual = 0.00626268, No Iterations 36
DILUPBiCG:  Solving for Uz, Initial residual = 0.902472, Final residual = 0.00562393, No Iterations 83
GAMG:  Solving for p, Initial residual = 0.0680265, Final residual = 0.000315977, No Iterations 1
time step continuity errors : sum local = 108.238, global = -0.0604695, cumulative = -0.0603884
GAMG:  Solving for p, Initial residual = 0.166241, Final residual = 7.95727e-05, No Iterations 1
time step continuity errors : sum local = 306.942, global = -0.0120188, cumulative = -0.0724072
DILUPBiCG:  Solving for omega, Initial residual = 0.998683, Final residual = 564.001, No Iterations 101
bounding omega, min: -9.5089e+20 max: 1.10262e+22 average: 5.97726e+15
DILUPBiCG:  Solving for k, Initial residual = 0.998501, Final residual = 0.00076919, No Iterations 1
PIMPLE: not converged within 20 iterations
ExecutionTime = 419.97 s  ClockTime = 420 s

fvSchemes

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(U)        Gauss linear 1;

}

divSchemes
{
    default         none;

    div(phi,U)      Gauss upwind;
    div(phi,k)      Gauss upwind;
    div(phi,omega) Gauss upwind;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear uncorrected;
}

interpolationSchemes
{
    default         linear;
    interpolate(HbyA) linear;

}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
    pcorr           ;
}


// ************************************************************************* //

lpz_michele · October 16, 2015, 07:09

Anyway I'm going to find what are the critics cell for Co

lpz_michele · October 16, 2015, 07:11

Running in parallel!

decomposeParDict

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      decomposeParDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

numberOfSubdomains     4;



method              scotch;

distributed         no;

roots             ();


// ************************************************************************* //

on bash i run it with:

mpirun -np 4 pimpleFoam -parallel > log &

RodriguezFatz · October 16, 2015, 07:21

Ok, and you run potentialFoam with the same command, right?

RodriguezFatz · October 16, 2015, 07:28

In the fvSolution, you put the "potential Flow" information into the "{" of the relavationFactors.
Since your PIMPLE does not converge, there might be a problem with the relaxation factors. I can't see any typo, but you should try if openFoam actually uses them. You could see that by running one iteration (save the log output) and running a second simulation beginning at exactly the same point but you comment out the "relaxationFactors" in your fvSolution. Now you should see a different log output (residuals). If they are the same, you got a problem...

lpz_michele · October 16, 2015, 10:34

Some fresh news! I'm writing on this to let you know and to help other foamers in their first steps.

1) I was making an error, as you mentioned. I believed that writing in the fvSolution, can automatically initialize the field.
So I followed the sequent step:

mpirun -np 4 potentialFoam -parallel > potentialFoam &

then

mpirun -np 4 pimpleFoam -parallel > pimpleFoam &

I think that is working and that pimpleFoam is reading the precedent log!

Code:

tesisti@ime042:~/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_ost_2> mpirun -np 4 potentialFoam -parallel > potentialFoam &
[1] 3668
tesisti@ime042:~/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_ost_2> mpirun -np 4 pimpleFoam -parallel > pimpleFoam &
[2] 3696
[1]   Done                    mpirun -np 4 potentialFoam -parallel > potentialFoam
tesisti@ime042:~/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_ost_2>

After running this Courant number went down, but PIMPLE converges only on first time step.

I'm going to write it as I try as you said before

Quote:

In the fvSolution, you put the "potential Flow" information into the "{" of the relavationFactors.
Since your PIMPLE does not converge, there might be a problem with the relaxation factors. I can't see any typo, but you should try if openFoam actually uses them. You could see that by running one iteration (save the log output) and running a second simulation beginning at exactly the same point but you comment out the "relaxationFactors" in your fvSolution. Now you should see a different log output (residuals). If they are the same, you got a problem...

lpz_michele · October 16, 2015, 10:50

Relaxation factors check:

Initial fvSolution

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    
      "(p|pFinal)"
     {
        solver           GAMG;
        tolerance        0;
        relTol           1e-2;

        maxIter         20;    

        smoother         DIC;
    
        nPreSweeps       1;
        nPostSweeps      1;
        nFinestSweeps    2;

        cacheAgglomeration true;

        nCellsInCoarsestLevel 1000;
        agglomerator     faceAreaPair;
        mergeLevels      1;
     }
     
     
 "(U|k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)"    
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        0;
        relTol           1.0e-2;
        maxIter         100;
        };
  


}

PIMPLE
{
    correctPhi        no;
    nOuterCorrectors    20;
    nCorrectors         2;
    nNonOrthogonalCorrectors 0;
    turbOnFinalIterOnly false;
    pRefCell         0;
    pRefValue        0;

residualControl
  {
    "(p|U|k|omega)"
     {
      tolerance 0;
      relTol 1e-2;
      }
  }
}

relaxationFactors
{
    fields
    {
       p        0.3;
    }
    equations
    {
  
    U        0.7;
    k        0.7;
    omega        0.7;
    }
    

}

potentialFlow
{
    nNonOrthogonalCorrectors 10;

}


// ************************************************************************* //

Then I commented out the relaxation factors

Code:

/*relaxationFactors
{
    fields
    {
       p        0.3;
    }
    equations
    {
  
    U        0.7;
    k        0.7;
    omega        0.7;
    }
    

}*/

Both the cases started from 0 and reached the first time step. Relaxation factors are working: in the second case there was a really fast divergence, that stopped calculating for a floating point exception

Code:

Courant Number mean: 0.000445304 max: 201.917
Time = 1.09083e-05

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 5.1893e-05, Final residual = 2.67094e-07, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.000254933, Final residual = 1.81991e-06, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.000254962, Final residual = 8.40649e-08, No Iterations 3
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00993248, No Iterations 6
time step continuity errors : sum local = 3.44445e-10, global = 1.47938e-13, cumulative = 1.47938e-13
GAMG:  Solving for p, Initial residual = 0.0475507, Final residual = 0.000425403, No Iterations 15
time step continuity errors : sum local = 3.3436e-11, global = -3.66907e-15, cumulative = 1.44269e-13
DILUPBiCG:  Solving for omega, Initial residual = 1.54384e-05, Final residual = 3.66276e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 0.000539063, No Iterations 1
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 4.38124e-08, Final residual = 3.96545e-10, No Iterations 8
DILUPBiCG:  Solving for Uy, Initial residual = 1.27014e-07, Final residual = 3.35928e-10, No Iterations 10
DILUPBiCG:  Solving for Uz, Initial residual = 8.26727e-08, Final residual = 2.56026e-10, No Iterations 8
GAMG:  Solving for p, Initial residual = 0.0883746, Final residual = 0.00119149, No Iterations 20
time step continuity errors : sum local = 9.0609e-11, global = -1.02659e-15, cumulative = 1.43242e-13
GAMG:  Solving for p, Initial residual = 0.125757, Final residual = 0.00286835, No Iterations 20
time step continuity errors : sum local = 2.37831e-10, global = 3.88479e-15, cumulative = 1.47127e-13
DILUPBiCG:  Solving for omega, Initial residual = 4.05753e-07, Final residual = 3.95868e-09, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 9.28993e-06, Final residual = 1.57289e-08, No Iterations 2
PIMPLE: iteration 3
DILUPBiCG:  Solving for Ux, Initial residual = 8.33934e-07, Final residual = 2.71442e-09, No Iterations 25
DILUPBiCG:  Solving for Uy, Initial residual = 3.5117e-06, Final residual = 2.71685e-08, No Iterations 22
DILUPBiCG:  Solving for Uz, Initial residual = 1.14694e-06, Final residual = 7.09277e-09, No Iterations 21
GAMG:  Solving for p, Initial residual = 0.195153, Final residual = 0.0111836, No Iterations 20
time step continuity errors : sum local = 9.85764e-10, global = -8.83366e-14, cumulative = 5.87901e-14
GAMG:  Solving for p, Initial residual = 0.256197, Final residual = 0.0105003, No Iterations 20
time step continuity errors : sum local = 2.50346e-09, global = 1.99406e-13, cumulative = 2.58196e-13
DILUPBiCG:  Solving for omega, Initial residual = 4.13675e-06, Final residual = 2.42137e-08, No Iterations 12
bounding omega, min: -1.14042e+06 max: 1.6918e+09 average: 64392
DILUPBiCG:  Solving for k, Initial residual = 4.71133e-06, Final residual = 3.93503e-08, No Iterations 9
bounding k, min: -1.94042 max: 201.739 average: 0.48003
PIMPLE: iteration 4
DILUPBiCG:  Solving for Ux, Initial residual = 5.05945e-05, Final residual = 4.55499e-07, No Iterations 17
DILUPBiCG:  Solving for Uy, Initial residual = 0.000132509, Final residual = 1.06735e-06, No Iterations 19
DILUPBiCG:  Solving for Uz, Initial residual = 7.25417e-05, Final residual = 7.1971e-07, No Iterations 20
GAMG:  Solving for p, Initial residual = 0.433663, Final residual = 0.0162722, No Iterations 20
time step continuity errors : sum local = 6.62259e-09, global = -1.96556e-13, cumulative = 6.16403e-14
GAMG:  Solving for p, Initial residual = 0.102424, Final residual = 0.00197756, No Iterations 20
time step continuity errors : sum local = 1.24862e-08, global = 1.71728e-12, cumulative = 1.77892e-12
DILUPBiCG:  Solving for omega, Initial residual = 0.99996, Final residual = 0.00493257, No Iterations 14
bounding omega, min: -7.12755e+07 max: 6.72841e+20 average: 3.88546e+15
DILUPBiCG:  Solving for k, Initial residual = 0.925177, Final residual = 7.76715e-10, No Iterations 1
bounding k, min: -2.96018e-09 max: 1545.46 average: 0.569534
PIMPLE: iteration 5
DILUPBiCG:  Solving for Ux, Initial residual = 0.0017509, Final residual = 1.51618e-05, No Iterations 27
DILUPBiCG:  Solving for Uy, Initial residual = 0.0054535, Final residual = 4.91513e-05, No Iterations 27
DILUPBiCG:  Solving for Uz, Initial residual = 0.00149593, Final residual = 1.15192e-05, No Iterations 29
GAMG:  Solving for p, Initial residual = 0.140851, Final residual = 0.00329212, No Iterations 20
time step continuity errors : sum local = 5.78573e-08, global = 8.80313e-11, cumulative = 8.98103e-11
GAMG:  Solving for p, Initial residual = 0.012873, Final residual = 0.000119379, No Iterations 19
time step continuity errors : sum local = 7.8355e-08, global = -2.33994e-10, cumulative = -1.44184e-10
DILUPBiCG:  Solving for omega, Initial residual = 1, Final residual = 1.80847e-14, No Iterations 1
bounding omega, min: -262144 max: 2.62256e+20 average: 2.69916e+14
DILUPBiCG:  Solving for k, Initial residual = 0.833271, Final residual = 1.69561e-06, No Iterations 1
PIMPLE: iteration 6
DILUPBiCG:  Solving for Ux, Initial residual = 0.00689325, Final residual = 6.68045e-05, No Iterations 51
DILUPBiCG:  Solving for Uy, Initial residual = 0.0365229, Final residual = 0.000300997, No Iterations 37
DILUPBiCG:  Solving for Uz, Initial residual = 0.0162189, Final residual = 0.000127679, No Iterations 36
GAMG:  Solving for p, Initial residual = 0.410032, Final residual = 0.00324479, No Iterations 7
time step continuity errors : sum local = 1.01285e-05, global = -3.63498e-08, cumulative = -3.6494e-08
GAMG:  Solving for p, Initial residual = 0.426942, Final residual = 0.00251159, No Iterations 1
time step continuity errors : sum local = 0.000357533, global = -1.10138e-06, cumulative = -1.13788e-06
DILUPBiCG:  Solving for omega, Initial residual = 1, Final residual = 1.70247e-12, No Iterations 1
bounding omega, min: -491520 max: 1.94418e+20 average: 3.79498e+14
DILUPBiCG:  Solving for k, Initial residual = 0.998041, Final residual = 6.03426e-08, No Iterations 1
bounding k, min: -3.9135e-05 max: 9.94741e+06 average: 7.32174
PIMPLE: iteration 7
DILUPBiCG:  Solving for Ux, Initial residual = 0.984971, Final residual = 0.000535075, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.984951, Final residual = 0.000303287, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.984989, Final residual = 0.000325319, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.999745, Final residual = 0.0097744, No Iterations 1
time step continuity errors : sum local = 19.1033, global = -0.00118446, cumulative = -0.0011856
GAMG:  Solving for p, Initial residual = 0.619635, Final residual = 0.00387419, No Iterations 2
time step continuity errors : sum local = 26.3709, global = 0.00229799, cumulative = 0.00111239
DILUPBiCG:  Solving for omega, Initial residual = 0.999847, Final residual = 7.54423e-08, No Iterations 1
bounding omega, min: -5.16443e+09 max: 5.22651e+22 average: 2.04419e+16
DILUPBiCG:  Solving for k, Initial residual = 0.603089, Final residual = 0.00267977, No Iterations 15
bounding k, min: -1.70743e+08 max: 3.33391e+11 average: 110687
PIMPLE: iteration 8
DILUPBiCG:  Solving for Ux, Initial residual = 0.999958, Final residual = 0.00956335, No Iterations 17
DILUPBiCG:  Solving for Uy, Initial residual = 0.999994, Final residual = 0.003945, No Iterations 20
DILUPBiCG:  Solving for Uz, Initial residual = 0.999978, Final residual = 0.00496545, No Iterations 20
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00683389, No Iterations 1
time step continuity errors : sum local = 3.18098e+11, global = 6.5654e+07, cumulative = 6.5654e+07
GAMG:  Solving for p, Initial residual = 0.0688653, Final residual = 0.00043779, No Iterations 1
time step continuity errors : sum local = 5.48273e+11, global = 1.26013e+08, cumulative = 1.91667e+08
DILUPBiCG:  Solving for omega, Initial residual = 0.999284, Final residual = 0.00870129, No Iterations 31
bounding omega, min: -3.30025e+42 max: 9.25553e+40 average: -1.29461e+36
DILUPBiCG:  Solving for k, Initial residual = 0.00188632, Final residual = 1.59517e-09, No Iterations 1
bounding k, min: -7.27495e-12 max: 3.19029e+24 average: 2.03602e+18
PIMPLE: iteration 9
DILUPBiCG:  Solving for Ux, Initial residual = 0.999994, Final residual = 0.00750712, No Iterations 22
DILUPBiCG:  Solving for Uy, Initial residual = 0.999987, Final residual = 0.00816343, No Iterations 11
DILUPBiCG:  Solving for Uz, Initial residual = 0.999996, Final residual = 0.00889624, No Iterations 19

lpz_michele · October 16, 2015, 12:10

As you stated here http://www.cfd-online.com/Forums/ope...tml#post568584

I first reduced by 0.1 the time step. This is actually my controlDict

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     simpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         0.7853976;

deltaT          0.00000109083;

writeControl    adjustableRunTime;

writeInterval   0.00000109083;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  no;

maxCo           1;

// ************************************************************************* //

Using the precedent fvSchemes

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(U)        Gauss linear 1;

}

divSchemes
{
    default         none;

    div(phi,U)      Gauss upwind;
    div(phi,k)      Gauss upwind;
    div(phi,omega) Gauss upwind;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear uncorrected;
}

interpolationSchemes
{
    default         linear;
    interpolate(HbyA) linear;

}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
    pcorr           ;
}


// ************************************************************************* //

I've set a strictly condition on pressure in fvSolution

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    
      "(p|pFinal)"
     {
        solver           GAMG;
        tolerance        1e-5;
        relTol           1e-2;

        maxIter         20;    

        smoother         DIC;
    
        nPreSweeps       1;
        nPostSweeps      1;
        nFinestSweeps    2;

        cacheAgglomeration true;

        nCellsInCoarsestLevel 1000;
        agglomerator     faceAreaPair;
        mergeLevels      1;
     }
     
     
 "(U|k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)"    
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        0;
        relTol           1.0e-2;
        maxIter         100;
        };
  


}

PIMPLE
{
    correctPhi        no;
    nOuterCorrectors    20;
    nCorrectors         2;
    nNonOrthogonalCorrectors 0;
    turbOnFinalIterOnly false;
    pRefCell         0;
    pRefValue        0;

residualControl
  {
    "(p|U|k|omega)"
     {
      tolerance 0;
      relTol 1e-2;
      }
  }
}

relaxationFactors
{
    fields
    {
       p        0.3;
    }
    equations
    {
  
    U        0.7;
    k        0.7;
    omega        0.7;
    }
    

}

potentialFlow
{
    nNonOrthogonalCorrectors 10;
}

// ************************************************************************* //

lpz_michele · October 16, 2015, 12:13

potentialFoam

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.1-262087cdf8db
Exec   : potentialFoam -parallel
Date   : Oct 16 2015
Time   : 16:36:55
Host   : "ime042"
PID    : 5875
Case   : /home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_ost_5
nProcs : 4
Slaves : 
3
(
"ime042.5876"
"ime042.5877"
"ime042.5878"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading field p

Reading field U

No finite volume options present


Calculating potential flow
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.00681261, No Iterations 12
GAMG:  Solving for p, Initial residual = 5.14105e-05, Final residual = 8.99497e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
GAMG:  Solving for p, Initial residual = 8.98793e-06, Final residual = 8.98793e-06, No Iterations 0
continuity error = 0.00059598
Interpolated U error = 4.03859e-07
ExecutionTime = 13.56 s  ClockTime = 13 s

End

Finalising parallel run

lpz_michele · October 16, 2015, 12:14

pimpleFoam log - part 1

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.1-262087cdf8db
Exec   : pimpleFoam -parallel
Date   : Oct 16 2015
Time   : 16:37:46
Host   : "ime042"
PID    : 5887
Case   : /home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_ost_5
nProcs : 4
Slaves : 
3
(
"ime042.5888"
"ime042.5889"
"ime042.5890"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing 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 kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1         0.85;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.856;
    gamma1          0.555556;
    gamma2          0.44;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

No finite volume options present


PIMPLE: max iterations = 20
    field "(p|U|k|omega)"    : relTol 0.01, tolerance 0


Starting time loop

Courant Number mean: 4.453e-05 max: 37.7229
Time = 1.09083e-06

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 3.56198e-06, Final residual = 1.33024e-08, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1.77937e-05, Final residual = 1.20302e-07, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1.78505e-05, Final residual = 3.97403e-08, No Iterations 2
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0545833, No Iterations 20
time step continuity errors : sum local = 6.65761e-11, global = -1.78856e-14, cumulative = -1.78856e-14
GAMG:  Solving for p, Initial residual = 0.244404, Final residual = 0.0139016, No Iterations 20
time step continuity errors : sum local = 4.93176e-11, global = -1.16467e-14, cumulative = -2.95323e-14
DILUPBiCG:  Solving for omega, Initial residual = 1.08174e-06, Final residual = 2.81527e-11, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 7.84715e-05, No Iterations 1
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 2.76007e-06, Final residual = 2.73864e-08, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 7.95074e-06, Final residual = 3.50348e-09, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 7.94814e-06, Final residual = 3.14482e-10, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.250384, Final residual = 0.0117232, No Iterations 20
time step continuity errors : sum local = 4.20127e-11, global = -5.571e-15, cumulative = -3.51033e-14
GAMG:  Solving for p, Initial residual = 0.126825, Final residual = 0.00643104, No Iterations 20
time step continuity errors : sum local = 3.22262e-11, global = -4.01136e-15, cumulative = -3.91146e-14
DILUPBiCG:  Solving for omega, Initial residual = 3.15242e-07, Final residual = 5.41003e-11, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00135712, Final residual = 3.01213e-07, No Iterations 1
PIMPLE: iteration 3
DILUPBiCG:  Solving for Ux, Initial residual = 2.2663e-06, Final residual = 4.26908e-10, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 5.04265e-06, Final residual = 5.00647e-08, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 5.048e-06, Final residual = 2.72827e-11, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.128728, Final residual = 0.00571153, No Iterations 20
time step continuity errors : sum local = 2.86291e-11, global = -3.29869e-15, cumulative = -4.24133e-14
GAMG:  Solving for p, Initial residual = 0.0785473, Final residual = 0.00378209, No Iterations 20
time step continuity errors : sum local = 2.20351e-11, global = -2.32816e-15, cumulative = -4.47415e-14
DILUPBiCG:  Solving for omega, Initial residual = 9.33508e-08, Final residual = 3.16818e-10, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.000190374, Final residual = 7.89671e-08, No Iterations 1
PIMPLE: iteration 4
DILUPBiCG:  Solving for Ux, Initial residual = 1.4985e-06, Final residual = 5.52762e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 3.13039e-06, Final residual = 1.76845e-09, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 3.13365e-06, Final residual = 1.93912e-08, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.079989, Final residual = 0.00324281, No Iterations 20
time step continuity errors : sum local = 1.88719e-11, global = -1.84788e-15, cumulative = -4.65894e-14
GAMG:  Solving for p, Initial residual = 0.0526402, Final residual = 0.00231035, No Iterations 20
time step continuity errors : sum local = 1.44831e-11, global = -1.24516e-15, cumulative = -4.78345e-14
DILUPBiCG:  Solving for omega, Initial residual = 2.75178e-08, Final residual = 2.16837e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 3.71172e-05, Final residual = 5.44215e-09, No Iterations 1
PIMPLE: iteration 5
DILUPBiCG:  Solving for Ux, Initial residual = 8.7526e-07, Final residual = 4.59003e-11, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 1.80013e-06, Final residual = 7.75921e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.79727e-06, Final residual = 5.99696e-09, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0535954, Final residual = 0.00202278, No Iterations 20
time step continuity errors : sum local = 1.26625e-11, global = -1.03178e-15, cumulative = -4.88663e-14
GAMG:  Solving for p, Initial residual = 0.0365161, Final residual = 0.00150673, No Iterations 20
time step continuity errors : sum local = 9.8014e-12, global = -6.75832e-16, cumulative = -4.95421e-14
DILUPBiCG:  Solving for omega, Initial residual = 8.13773e-09, Final residual = 6.50322e-11, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 8.40525e-06, Final residual = 6.76449e-10, No Iterations 1
PIMPLE: iteration 6
DILUPBiCG:  Solving for Ux, Initial residual = 4.74742e-07, Final residual = 9.03314e-10, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 9.74642e-07, Final residual = 1.87282e-11, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 9.68138e-07, Final residual = 3.85607e-09, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0368236, Final residual = 0.00132986, No Iterations 20
time step continuity errors : sum local = 8.63954e-12, global = -5.53251e-16, cumulative = -5.00954e-14
GAMG:  Solving for p, Initial residual = 0.0255278, Final residual = 0.0010117, No Iterations 20
time step continuity errors : sum local = 6.70878e-12, global = -3.63423e-16, cumulative = -5.04588e-14
DILUPBiCG:  Solving for omega, Initial residual = 2.48855e-09, Final residual = 1.75312e-11, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 2.07143e-06, Final residual = 3.2108e-10, No Iterations 1
PIMPLE: iteration 7
DILUPBiCG:  Solving for Ux, Initial residual = 2.4588e-07, Final residual = 1.29179e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 5.0857e-07, Final residual = 2.94872e-10, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 5.01872e-07, Final residual = 1.18822e-09, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0254199, Final residual = 0.000903709, No Iterations 20
time step continuity errors : sum local = 5.98628e-12, global = -3.08861e-16, cumulative = -5.07677e-14
GAMG:  Solving for p, Initial residual = 0.0177908, Final residual = 0.000691406, No Iterations 20
time step continuity errors : sum local = 4.63208e-12, global = -1.91796e-16, cumulative = -5.09595e-14
DILUPBiCG:  Solving for omega, Initial residual = 8.0241e-10, Final residual = 5.07339e-12, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 5.42468e-07, Final residual = 2.60713e-10, No Iterations 1
PIMPLE: iteration 8
DILUPBiCG:  Solving for Ux, Initial residual = 1.24042e-07, Final residual = 9.36046e-11, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 2.62234e-07, Final residual = 2.0577e-10, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 2.56914e-07, Final residual = 2.0301e-09, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0176017, Final residual = 0.00062884, No Iterations 20
time step continuity errors : sum local = 4.20965e-12, global = -1.76795e-16, cumulative = -5.11363e-14
GAMG:  Solving for p, Initial residual = 0.0123897, Final residual = 0.000480735, No Iterations 20
time step continuity errors : sum local = 3.23905e-12, global = -9.99608e-17, cumulative = -5.12362e-14
DILUPBiCG:  Solving for omega, Initial residual = 2.80956e-10, Final residual = 1.64777e-12, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.50293e-07, Final residual = 2.38917e-10, No Iterations 1
PIMPLE: iteration 9
DILUPBiCG:  Solving for Ux, Initial residual = 6.20962e-08, Final residual = 1.63488e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1.37383e-07, Final residual = 4.79795e-10, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1.33574e-07, Final residual = 1.05404e-11, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0122238, Final residual = 0.000441077, No Iterations 20
time step continuity errors : sum local = 2.97022e-12, global = -9.49174e-17, cumulative = -5.13311e-14
GAMG:  Solving for p, Initial residual = 0.00863277, Final residual = 0.000338072, No Iterations 20
time step continuity errors : sum local = 2.28529e-12, global = -5.0999e-17, cumulative = -5.13821e-14
DILUPBiCG:  Solving for omega, Initial residual = 1.12514e-10, Final residual = 6.26722e-13, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 4.48064e-08, Final residual = 2.12997e-10, No Iterations 1
PIMPLE: iteration 10
DILUPBiCG:  Solving for Ux, Initial residual = 4.50699e-08, Final residual = 2.2968e-11, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 1.06932e-07, Final residual = 9.55202e-10, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.03454e-07, Final residual = 6.06791e-10, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0935803, Final residual = 0.000910734, No Iterations 8
time step continuity errors : sum local = 8.17016e-12, global = -5.17538e-14, cumulative = -1.03136e-13
GAMG:  Solving for p, Initial residual = 0.00198551, Final residual = 0.000133573, No Iterations 20
time step continuity errors : sum local = 1.14564e-12, global = -3.87571e-17, cumulative = -1.03175e-13
DILUPBiCG:  Solving for omega, Initial residual = 8.72741e-11, Final residual = 2.14913e-13, No Iterations 3
DILUPBiCG:  Solving for k, Initial residual = 2.25658e-08, Final residual = 1.44388e-10, No Iterations 3
PIMPLE: converged in 10 iterations
ExecutionTime = 186.93 s  ClockTime = 187 s

lpz_michele · October 16, 2015, 12:15

Code:

Courant Number mean: 4.45303e-05 max: 20.623
Time = 2.18166e-06

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 2.01829e-05, Final residual = 4.48658e-08, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 4.58613e-05, Final residual = 1.18329e-07, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 4.57788e-05, Final residual = 2.81925e-08, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.172653, Final residual = 0.00722223, No Iterations 20
time step continuity errors : sum local = 4.57833e-11, global = 9.15527e-15, cumulative = -9.40195e-14
GAMG:  Solving for p, Initial residual = 0.158842, Final residual = 0.00718898, No Iterations 20
time step continuity errors : sum local = 3.51171e-11, global = 6.68425e-15, cumulative = -8.73352e-14
DILUPBiCG:  Solving for omega, Initial residual = 9.8987e-07, Final residual = 1.70521e-09, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.000476854, Final residual = 8.28099e-07, No Iterations 1
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 4.38341e-06, Final residual = 7.60603e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1.03332e-05, Final residual = 1.87359e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.02862e-05, Final residual = 1.29873e-08, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.18753, Final residual = 0.00850147, No Iterations 20
time step continuity errors : sum local = 4.25285e-11, global = 6.98942e-15, cumulative = -8.03458e-14
GAMG:  Solving for p, Initial residual = 0.180066, Final residual = 0.0087616, No Iterations 20
time step continuity errors : sum local = 3.24236e-11, global = 4.86326e-15, cumulative = -7.54825e-14
DILUPBiCG:  Solving for omega, Initial residual = 2.875e-07, Final residual = 2.47576e-09, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.000100847, Final residual = 1.73378e-07, No Iterations 1
PIMPLE: iteration 3
DILUPBiCG:  Solving for Ux, Initial residual = 7.01483e-07, Final residual = 3.94168e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 2.14442e-06, Final residual = 1.65964e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 2.16971e-06, Final residual = 7.10652e-09, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.17134, Final residual = 0.00867445, No Iterations 20
time step continuity errors : sum local = 3.17189e-11, global = 4.94761e-15, cumulative = -7.05349e-14
GAMG:  Solving for p, Initial residual = 0.159231, Final residual = 0.00865186, No Iterations 20
time step continuity errors : sum local = 2.42375e-11, global = 3.3595e-15, cumulative = -6.71754e-14
DILUPBiCG:  Solving for omega, Initial residual = 8.3793e-08, Final residual = 7.73325e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 2.54239e-05, Final residual = 5.73256e-08, No Iterations 1
PIMPLE: iteration 4
DILUPBiCG:  Solving for Ux, Initial residual = 1.36798e-06, Final residual = 7.33894e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 2.86256e-06, Final residual = 3.20603e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 2.89324e-06, Final residual = 3.45842e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.139523, Final residual = 0.00798955, No Iterations 20
time step continuity errors : sum local = 2.18564e-11, global = 2.91487e-15, cumulative = -6.42605e-14
GAMG:  Solving for p, Initial residual = 0.121528, Final residual = 0.00747932, No Iterations 20
time step continuity errors : sum local = 1.67675e-11, global = 1.93698e-15, cumulative = -6.23236e-14
DILUPBiCG:  Solving for omega, Initial residual = 2.4564e-08, Final residual = 2.32798e-10, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 6.9101e-06, Final residual = 1.26926e-08, No Iterations 1
PIMPLE: iteration 5
DILUPBiCG:  Solving for Ux, Initial residual = 1.10448e-06, Final residual = 2.34607e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 2.24041e-06, Final residual = 1.31573e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 2.25776e-06, Final residual = 2.05578e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.107127, Final residual = 0.00668879, No Iterations 20
time step continuity errors : sum local = 1.4682e-11, global = 1.70257e-15, cumulative = -6.0621e-14
GAMG:  Solving for p, Initial residual = 0.0871417, Final residual = 0.00583061, No Iterations 20
time step continuity errors : sum local = 1.12561e-11, global = 1.0663e-15, cumulative = -5.95547e-14
DILUPBiCG:  Solving for omega, Initial residual = 7.24216e-09, Final residual = 3.4381e-11, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1.94571e-06, Final residual = 2.26542e-09, No Iterations 1
PIMPLE: iteration 6
DILUPBiCG:  Solving for Ux, Initial residual = 7.02727e-07, Final residual = 1.31747e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1.41072e-06, Final residual = 1.02309e-08, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1.41593e-06, Final residual = 2.45615e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0851295, Final residual = 0.00510554, No Iterations 20
time step continuity errors : sum local = 9.71929e-12, global = 8.47027e-16, cumulative = -5.87077e-14
GAMG:  Solving for p, Initial residual = 0.0583794, Final residual = 0.00376808, No Iterations 20
time step continuity errors : sum local = 7.47691e-12, global = 5.26086e-16, cumulative = -5.81816e-14
DILUPBiCG:  Solving for omega, Initial residual = 2.17649e-09, Final residual = 1.21114e-11, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 5.57925e-07, Final residual = 1.09515e-09, No Iterations 1
PIMPLE: iteration 7
DILUPBiCG:  Solving for Ux, Initial residual = 3.98407e-07, Final residual = 1.28365e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 7.94267e-07, Final residual = 3.43745e-09, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 7.931e-07, Final residual = 2.90238e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0598663, Final residual = 0.00331439, No Iterations 20
time step continuity errors : sum local = 6.53467e-12, global = 4.34312e-16, cumulative = -5.77473e-14
GAMG:  Solving for p, Initial residual = 0.0413132, Final residual = 0.00250954, No Iterations 20
time step continuity errors : sum local = 5.11776e-12, global = 3.12394e-16, cumulative = -5.74349e-14
DILUPBiCG:  Solving for omega, Initial residual = 6.78868e-10, Final residual = 8.64119e-13, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1.60556e-07, Final residual = 5.85683e-10, No Iterations 1
PIMPLE: iteration 8
DILUPBiCG:  Solving for Ux, Initial residual = 2.10485e-07, Final residual = 7.56605e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 4.17688e-07, Final residual = 1.25136e-10, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 4.14293e-07, Final residual = 2.64208e-10, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0421072, Final residual = 0.00220901, No Iterations 20
time step continuity errors : sum local = 4.49198e-12, global = 2.46408e-16, cumulative = -5.71885e-14
GAMG:  Solving for p, Initial residual = 0.0295831, Final residual = 0.00169982, No Iterations 20
time step continuity errors : sum local = 3.5079e-12, global = 1.62539e-16, cumulative = -5.70259e-14
DILUPBiCG:  Solving for omega, Initial residual = 2.15193e-10, Final residual = 1.60704e-12, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 4.58345e-08, Final residual = 2.45942e-10, No Iterations 1
PIMPLE: iteration 9
DILUPBiCG:  Solving for Ux, Initial residual = 1.06134e-07, Final residual = 5.51975e-10, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 2.10705e-07, Final residual = 9.12304e-11, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 2.07207e-07, Final residual = 3.31069e-10, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.0296584, Final residual = 0.00149217, No Iterations 20
time step continuity errors : sum local = 3.07575e-12, global = 1.17762e-16, cumulative = -5.69082e-14
GAMG:  Solving for p, Initial residual = 0.0210013, Final residual = 0.00115708, No Iterations 20
time step continuity errors : sum local = 2.40009e-12, global = 7.40392e-17, cumulative = -5.68341e-14
DILUPBiCG:  Solving for omega, Initial residual = 7.22376e-11, Final residual = 5.88889e-13, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 1.32949e-08, Final residual = 6.86659e-11, No Iterations 1
PIMPLE: iteration 10
DILUPBiCG:  Solving for Ux, Initial residual = 5.19546e-08, Final residual = 6.23999e-12, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1.0435e-07, Final residual = 3.88881e-10, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.01569e-07, Final residual = 4.0489e-10, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0208409, Final residual = 0.00101005, No Iterations 20
time step continuity errors : sum local = 2.09421e-12, global = 4.77486e-17, cumulative = -5.67864e-14
GAMG:  Solving for p, Initial residual = 0.0148078, Final residual = 0.000786774, No Iterations 20
time step continuity errors : sum local = 1.63558e-12, global = 2.91142e-17, cumulative = -5.67573e-14
DILUPBiCG:  Solving for omega, Initial residual = 2.88359e-11, Final residual = 2.55424e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 4.32863e-09, Final residual = 1.8068e-12, No Iterations 2
PIMPLE: iteration 11
DILUPBiCG:  Solving for Ux, Initial residual = 2.51113e-08, Final residual = 8.10394e-11, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 5.21966e-08, Final residual = 7.76541e-11, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 5.0185e-08, Final residual = 6.46474e-11, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0145917, Final residual = 0.000688791, No Iterations 20
time step continuity errors : sum local = 1.43172e-12, global = 1.77394e-17, cumulative = -5.67395e-14
GAMG:  Solving for p, Initial residual = 0.0103817, Final residual = 0.000542089, No Iterations 20
time step continuity errors : sum local = 1.12803e-12, global = 1.13038e-17, cumulative = -5.67282e-14
DILUPBiCG:  Solving for omega, Initial residual = 1.43633e-11, Final residual = 2.65357e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 2.16785e-09, Final residual = 2.63317e-12, No Iterations 2
PIMPLE: iteration 12
DILUPBiCG:  Solving for Ux, Initial residual = 1.22392e-08, Final residual = 1.03571e-11, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 2.71971e-08, Final residual = 1.91492e-10, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 2.58707e-08, Final residual = 1.29402e-10, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.0101924, Final residual = 0.000485888, No Iterations 20
time step continuity errors : sum local = 1.01108e-12, global = 7.89178e-18, cumulative = -5.67203e-14
GAMG:  Solving for p, Initial residual = 0.00725867, Final residual = 0.000381229, No Iterations 20
time step continuity errors : sum local = 7.93681e-13, global = 4.36989e-18, cumulative = -5.6716e-14
DILUPBiCG:  Solving for omega, Initial residual = 9.89107e-12, Final residual = 2.77654e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.80208e-09, Final residual = 3.11185e-12, No Iterations 2
PIMPLE: iteration 13
DILUPBiCG:  Solving for Ux, Initial residual = 6.17815e-09, Final residual = 2.99393e-11, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1.50158e-08, Final residual = 3.76586e-11, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.4198e-08, Final residual = 9.78218e-12, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00711642, Final residual = 0.000341587, No Iterations 20
time step continuity errors : sum local = 7.11174e-13, global = 2.52633e-18, cumulative = -5.67134e-14
GAMG:  Solving for p, Initial residual = 0.00507087, Final residual = 0.000267721, No Iterations 20
time step continuity errors : sum local = 5.57517e-13, global = 9.78865e-19, cumulative = -5.67125e-14
DILUPBiCG:  Solving for omega, Initial residual = 9.10627e-12, Final residual = 2.58461e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.56521e-09, Final residual = 4.97523e-12, No Iterations 2
PIMPLE: iteration 14
DILUPBiCG:  Solving for Ux, Initial residual = 3.30064e-09, Final residual = 2.68999e-12, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 8.75134e-09, Final residual = 2.26527e-12, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 8.26568e-09, Final residual = 6.36412e-11, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.00497544, Final residual = 0.000239656, No Iterations 20
time step continuity errors : sum local = 4.99092e-13, global = 8.16531e-20, cumulative = -5.67124e-14
GAMG:  Solving for p, Initial residual = 0.00354679, Final residual = 0.000187809, No Iterations 20
time step continuity errors : sum local = 3.91169e-13, global = -4.44975e-19, cumulative = -5.67128e-14
DILUPBiCG:  Solving for omega, Initial residual = 7.76877e-12, Final residual = 2.72794e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.69025e-09, Final residual = 5.96954e-12, No Iterations 2
PIMPLE: iteration 15
DILUPBiCG:  Solving for Ux, Initial residual = 1.87421e-09, Final residual = 4.4777e-12, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 5.3107e-09, Final residual = 3.64066e-12, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 5.02056e-09, Final residual = 2.89241e-11, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00348813, Final residual = 0.000168407, No Iterations 20
time step continuity errors : sum local = 3.5077e-13, global = -7.99577e-19, cumulative = -5.67136e-14
GAMG:  Solving for p, Initial residual = 0.00248732, Final residual = 0.000132416, No Iterations 20
time step continuity errors : sum local = 2.75829e-13, global = -8.76065e-19, cumulative = -5.67145e-14
DILUPBiCG:  Solving for omega, Initial residual = 7.73304e-12, Final residual = 3.49189e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.59603e-09, Final residual = 6.93749e-12, No Iterations 2
PIMPLE: iteration 16
DILUPBiCG:  Solving for Ux, Initial residual = 1.11988e-09, Final residual = 2.45525e-12, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 3.31412e-09, Final residual = 1.23905e-13, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 3.13656e-09, Final residual = 1.4239e-11, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00245231, Final residual = 0.000119934, No Iterations 20
time step continuity errors : sum local = 2.49835e-13, global = -9.64969e-19, cumulative = -5.67155e-14
GAMG:  Solving for p, Initial residual = 0.00174952, Final residual = 9.45764e-05, No Iterations 20
time step continuity errors : sum local = 1.97023e-13, global = -9.06109e-19, cumulative = -5.67164e-14
DILUPBiCG:  Solving for omega, Initial residual = 6.69829e-12, Final residual = 3.62339e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.68986e-09, Final residual = 1.01524e-11, No Iterations 2
PIMPLE: iteration 17
DILUPBiCG:  Solving for Ux, Initial residual = 6.94211e-10, Final residual = 8.94132e-13, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 2.10917e-09, Final residual = 2.46274e-12, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.99698e-09, Final residual = 2.06442e-12, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00172909, Final residual = 8.59939e-05, No Iterations 20
time step continuity errors : sum local = 1.79147e-13, global = -9.2428e-19, cumulative = -5.67173e-14
GAMG:  Solving for p, Initial residual = 0.00123404, Final residual = 6.77688e-05, No Iterations 20
time step continuity errors : sum local = 1.41185e-13, global = -8.17118e-19, cumulative = -5.67181e-14
DILUPBiCG:  Solving for omega, Initial residual = 6.67461e-12, Final residual = 3.71905e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.65365e-09, Final residual = 9.26283e-12, No Iterations 2
PIMPLE: iteration 18
DILUPBiCG:  Solving for Ux, Initial residual = 4.41979e-10, Final residual = 3.17188e-13, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1.36394e-09, Final residual = 2.4338e-12, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.29087e-09, Final residual = 9.49721e-13, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00122266, Final residual = 6.16876e-05, No Iterations 20
time step continuity errors : sum local = 1.28518e-13, global = -7.97639e-19, cumulative = -5.67189e-14
GAMG:  Solving for p, Initial residual = 0.000873022, Final residual = 4.85709e-05, No Iterations 20
time step continuity errors : sum local = 1.01194e-13, global = -6.81368e-19, cumulative = -5.67196e-14
DILUPBiCG:  Solving for omega, Initial residual = 6.12744e-12, Final residual = 4.18543e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.71493e-09, Final residual = 1.22446e-11, No Iterations 2
PIMPLE: iteration 19
DILUPBiCG:  Solving for Ux, Initial residual = 2.8718e-10, Final residual = 5.39153e-14, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 8.94917e-10, Final residual = 1.5274e-12, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 8.45365e-10, Final residual = 8.00087e-12, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.000867173, Final residual = 4.43466e-05, No Iterations 20
time step continuity errors : sum local = 9.23937e-14, global = -6.4974e-19, cumulative = -5.67202e-14
GAMG:  Solving for p, Initial residual = 0.000619553, Final residual = 3.48971e-05, No Iterations 20
time step continuity errors : sum local = 7.27078e-14, global = -5.31361e-19, cumulative = -5.67208e-14
DILUPBiCG:  Solving for omega, Initial residual = 5.82169e-12, Final residual = 3.91473e-14, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 1.81952e-09, Final residual = 1.09559e-11, No Iterations 2
PIMPLE: iteration 20
DILUPBiCG:  Solving for Ux, Initial residual = 2.70768e-10, Final residual = 3.35663e-13, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 8.49185e-10, Final residual = 8.19636e-12, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 8.00204e-10, Final residual = 5.51949e-12, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.40613, Final residual = 0.00372296, No Iterations 5
time step continuity errors : sum local = 8.22605e-12, global = -6.09965e-14, cumulative = -1.17717e-13
GAMG:  Solving for p, Initial residual = 0.0108458, Final residual = 0.000605357, No Iterations 20
time step continuity errors : sum local = 7.55832e-13, global = -2.66701e-17, cumulative = -1.17744e-13
DILUPBiCG:  Solving for omega, Initial residual = 3.88552e-11, Final residual = 3.85356e-13, No Iterations 4
DILUPBiCG:  Solving for k, Initial residual = 2.0178e-09, Final residual = 8.02294e-12, No Iterations 5
PIMPLE: not converged within 20 iterations
ExecutionTime = 547.42 s  ClockTime = 547 s

lpz_michele · October 16, 2015, 12:16

Code:

Courant Number mean: 4.45303e-05 max: 18.5461
Time = 3.27249e-06

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 4.278e-06, Final residual = 1.69588e-08, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1.72667e-05, Final residual = 6.02362e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1.7335e-05, Final residual = 4.41677e-08, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.542532, Final residual = 0.00501066, No Iterations 8
time step continuity errors : sum local = 7.10968e-12, global = 1.61861e-14, cumulative = -1.01558e-13
GAMG:  Solving for p, Initial residual = 0.328845, Final residual = 0.00305906, No Iterations 8
time step continuity errors : sum local = 5.01042e-12, global = 1.27342e-14, cumulative = -8.88236e-14
DILUPBiCG:  Solving for omega, Initial residual = 9.42601e-07, Final residual = 1.61885e-09, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 3.25602e-05, Final residual = 1.11303e-07, No Iterations 2
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 1.34112e-06, Final residual = 5.07226e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 5.22317e-06, Final residual = 1.53346e-08, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 5.25385e-06, Final residual = 1.48517e-08, No Iterations 2

lpz_michele · October 16, 2015, 12:21

Finally I analized the Co.

OpenFOAM provides a native utility. After reconstructPar, in the bash type:

Co

Opened in paraFoam, and used Co as a volScalarField at first time step

Any ideas?

RodriguezFatz · October 19, 2015, 03:50

1) In your last log post, PIMPLE seems to converge. Don't you think so?
2) You Co is probably 20 because some cells are very small where the flow is fast, so why do you think this is wrong? If PIMPLE remains stable I don't see any problem.

October 16, 2015, 07:21		#44
RodriguezFatz Senior Member Philipp Join Date: Jun 2011 Location: Germany Posts: 1,297 Rep Power: 27	Ok, and you run potentialFoam with the same command, right? __________________ The skeleton ran out of shampoo in the shower.

October 16, 2015, 07:28		#45
RodriguezFatz Senior Member Philipp Join Date: Jun 2011 Location: Germany Posts: 1,297 Rep Power: 27	In the fvSolution, you put the "potential Flow" information into the "{" of the relavationFactors. Since your PIMPLE does not converge, there might be a problem with the relaxation factors. I can't see any typo, but you should try if openFoam actually uses them. You could see that by running one iteration (save the log output) and running a second simulation beginning at exactly the same point but you comment out the "relaxationFactors" in your fvSolution. Now you should see a different log output (residuals). If they are the same, you got a problem... lpz_michele likes this. __________________ The skeleton ran out of shampoo in the shower.

October 19, 2015, 03:50		#54
RodriguezFatz Senior Member Philipp Join Date: Jun 2011 Location: Germany Posts: 1,297 Rep Power: 27	1) In your last log post, PIMPLE seems to converge. Don't you think so? 2) You Co is probably 20 because some cells are very small where the flow is fast, so why do you think this is wrong? If PIMPLE remains stable I don't see any problem. __________________ The skeleton ran out of shampoo in the shower.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
error compiling modified applications	yvyan	OpenFOAM Programming & Development	21	March 1, 2016 05:53
[OpenFOAM] Native ParaView Reader Bugs	tj22	ParaView	270	January 4, 2016 12:39
[General] 2 datas on one plot	Akuji	ParaView	46	December 1, 2013 15:06
Floating point error and divergence detected	aannjj	FLUENT	0	July 2, 2013 04:44
[Gmsh] Gmsh and samplesurface	touf	OpenFOAM Meshing & Mesh Conversion	2	December 10, 2007 03:27

October 16, 2015, 07:09		#42
lpz_michele Member Mike Lob Join Date: Oct 2015 Location: It Posts: 51 Rep Power: 11	Anyway I'm going to find what are the critics cell for Co

October 16, 2015, 12:21		#53
lpz_michele Member Mike Lob Join Date: Oct 2015 Location: It Posts: 51 Rep Power: 11	Finally I analized the Co. OpenFOAM provides a native utility. After reconstructPar, in the bash type: Co Opened in paraFoam, and used Co as a volScalarField at first time step Any ideas?