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

Floating point exception error

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree2Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 15, 2015, 05:25
Default Floating point exception error
  #1
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
Hi everyone! I'm tring to simulate the flowfield on a stationary VAWT for some debugging. In the complete simulation (VAWT rotating) iterations start, but do not make sense. So I'm trying the same dictionaries with no rotations. Using pimpleFoam unstead of pimpleDyMFoam.

After running the first iteration for U, this appears:

Code:
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib64/libc.so.6"
#3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
#6  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#7  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#9  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10  
 at ??:?
#11  __libc_start_main in "/lib64/libc.so.6"
#12  
 at /home/abuild/rpmbuild/BUILD/glibc-2.19/csu/../sysdeps/x86_64/start.S:125
Floating point exception
How can I fix it?
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 05:40
Default
  #2
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
We need more information.
1) How does the rest of the log file look like?
2) Boundary conditions
3) fvSchemes
4) fvSolution
Cheers,
Philipp.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   October 15, 2015, 06:19
Default
  #3
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
rest of the log file
Code:
tesisti@ime042:~/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./flussoOltreOstacolo_tent1_laminare> of231
tesisti@ime042:~/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./flussoOltreOstacolo_tent1_laminare> pimpleFoam 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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
Date   : Oct 15 2015
Time   : 11:26:00
Host   : "ime042"
PID    : 27789
Case   : /home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./flussoOltreOstacolo_tent1_laminare
nProcs : 1
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 = 100                                                                                                                                                                                                                                                   
    field "(p|U|k|omega)"       : relTol 0, tolerance 0.01                                                                                                                                                                                                                     
                                                                                                                                                                                                                                                                               
                                                                                                                                                                                                                                                                               
Starting time loop                                                                                                                                                                                                                                                             
                                                                                                                                                                                                                                                                               
forces forces:                                                                                                                                                                                                                                                                 
    Not including porosity effects                                                                                                                                                                                                                                             
                                                                                                                                                                                                                                                                               
Courant Number mean: 2.72708e-09 max: 2.97915e-07                                                                                                                                                                                                                              
Time = 1.09083e-05                                                                                                                                                                                                                                                             
                                                                                                                                                                                                                                                                               
PIMPLE: iteration 1                                                                                                                                                                                                                                                            
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 8.13576e-14, No Iterations 1                                                                                                                                                                             
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 6.4221e-14, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 6.59122e-14, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
...ETC...
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 linearUpwind grad(U);
    div(phi,k)      Gauss upwind;
    div(phi,omega) Gauss upwind;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
    interpolate(HbyA) linear;

}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p               ;
    pcorr           ;
}


// ************************************************************************* //
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
{
    pcorr
    {
        solver GAMG;
    smoother GaussSeidel;
    nPreSweeps 0;
    nPostSweeps 2;
    cacheAgglomeration off;
    agglomerator faceAreaPair;
    nCellsInCoarsestLevel 10;
    mergeLevels 1;

    tolerance 0.02;
    relTol 0;
    }

    p
    {
    $pcorr;
    tolerance 1e-07;
    relTol 0.01;
    }  
    
    pFinal
    {
        $p;
    tolerance     1e-07;
        relTol          0;
    }


    U
    {
    solver        smoothSolver;
    smoother    GaussSeidel;
    tolerance    1e-8;
    relTol        0.1;

    }
    
    UFinal
    {
    solver        smoothSolver;
    smoother    GaussSeidel;
    tolerance    1e-8;
    relTol        0.1;

    }
    
    "(k|omega)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       1e-7;
        relTol          0.1;
    }

    "(k|omega)Final"
    {
        $U;
    tolerance     1e-7;
        relTol          0;
    }


}

PIMPLE
{
    nNonOrthogonalCorrectors 0;
    correctPhi        no;
    nOuterCorrectors    100;
    nCorrectors         2;
    pRefCell         0;
    pRefValue        0;

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

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


// ************************************************************************* //
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 06:22
Default
  #4
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
boundary
Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}

15        
(
    cubodx
    {
        type patch;
        nFaces 592;
        startFace 7406234;
    }
    cubosopra
    {
        type patch;
        nFaces 588;
        startFace 7406826;
    }
    cubosotto
    {
        type patch;
        nFaces 588;
        startFace 7407414;
    }
    cubosx
    {
        type patch;
        nFaces 590;
        startFace 7408002;
    }
    inlet
    {
        type patch;
        nFaces 478;
        startFace 7408592;
    }
    outlet
    {
        type patch;
        nFaces 476;
        startFace 7409070;
    }
    pala1
    {
        type wall;
        nFaces 9506;
        startFace 7409546;
    }
    pala1sopra
    {
        type wall;
        nFaces 301;
        startFace 7419052;
    }
    pala1sotto
    {
        type wall;
        nFaces 301;
        startFace 7419353;
    }
    pala2
    {
        type wall;
        nFaces 9506;
        startFace 7419654;
    }
    pala2sopra
    {
        type wall;
        nFaces 301;
        startFace 7429160;
    }
    pala2sotto
    {
        type wall;
        nFaces 301;
        startFace 7429461;
    }
    pala3
    {
        type wall;
        nFaces 9506;
        startFace 7429762;
    }
    pala3sopra
    {
        type wall;
        nFaces 301;
        startFace 7439268;
    }
    pala3sotto
    {
        type wall;
        nFaces 301;
        startFace 7439569;
    }
)
k
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       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0.375;

boundaryField
{
    cubodx
    {
        type            zeroGradient;
    }
    cubosopra
    {
        type            zeroGradient;
    }
    cubosotto
    {
        type            zeroGradient;
    }
    cubosx
    {
        type            zeroGradient;
    }
    inlet
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity       0.05;
        value           uniform 0.375;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0.375;
        value           uniform 0.375;
    }
    pala1
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
    pala1sopra
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
    pala1sotto
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
    pala2
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
    pala2sopra
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
    pala2sotto
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
    pala3
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
    pala3sopra
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
    pala3sotto
    {
        type            kqRWallFunction;
        value           uniform 0.375;
    }
}


// ************************************************************************* //
omega
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       volScalarField;
    location    "0";
    object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 -1 0 0 0 0];

internalField   uniform 1118.03;

boundaryField
{
    cubodx
    {
        type            zeroGradient;
    }
    cubosopra
    {
        type            zeroGradient;
    }
    cubosotto
    {
        type            zeroGradient;
    }
    cubosx
    {
        type            zeroGradient;
    }
    inlet
    {
        type            fixedValue;
        value           uniform 1118.03;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 1118.03;
        value           uniform 1118.03;
    }
    pala1
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
    pala1sopra
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
    pala1sotto
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
    pala2
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
    pala2sopra
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
    pala2sotto
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
    pala3
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
    pala3sopra
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
    pala3sotto
    {
        type            omegaWallFunction;
        value           uniform 1118.03;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        beta1           0.075;
    }
}


// ************************************************************************* //
p
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       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    cubodx
    {
        type zeroGradient;
        
    }
    
    
    cubosopra
    {
        type zeroGradient;
        
    }
    
    
    cubosotto
    {
        type zeroGradient;
       
    }
    
    
    cubosx
    {
        type zeroGradient;
        
    }

    inlet
    {
        type zeroGradient;       
    }
    
    
    outlet
    {
        type fixedValue;
    value $internalField;        
    }
  
  
    pala1
    {
        type zeroGradient;
    }
    
    
    pala1sopra
    {
        type zeroGradient;        
    }
    
    
    pala1sotto
    {
    type zeroGradient;
    }
    
    
    pala2
    {
        type zeroGradient;      
    }
    
    
    pala2sopra
    {
        
       type    zeroGradient;
    }
    
    
    pala2sotto
    {
       type  zeroGradient; 
    }
    
    
    pala3
    {
        type    zeroGradient;
    }
    
    
    pala3sopra
    {
        type zeroGradient;      
    }
    
    
    pala3sotto
    {
     type zeroGradient;
    }
  
}

// ************************************************************************* //
U
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       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    cubodx
    {
        type slip;        
    }
    
    
    cubosopra
    {
        type slip;        
    }
    
    
    
    cubosotto
    {
        type slip;       
    }
    
    
    
    cubosx
    {
        type slip;        
    }
 
    inlet
    {
        type fixedValue;
    value uniform (10.0 0 0);       
    }
    
    
    
    outlet
    {
        type zeroGradient;
        
    }
  
    pala1
    {
        type fixedValue;
    value $internalField;
    }
    
    
    pala1sopra
    {
        type fixedValue;
    value $internalField;
        
    }
    
    
    pala1sotto
    {
    type fixedValue;
    value $internalField; 
    }
    
    
    
    pala2
    {
        type fixedValue;
    value $internalField;
      
    }
    
    
    
    pala2sopra
    {
        
       type    fixedValue;
       value     $internalField;
    }
    
    
    
    pala2sotto
    {
       type  fixedValue;
       value $internalField; 
    }
    
    
    
    pala3
    {
        type    fixedValue ;
    value $internalField;
    }
    
    
    
    pala3sopra
    {
        type fixedValue;
    value $internalField;
      
    }
    
    
    
    pala3sotto
    {
     type fixedValue;
     value $internalField;
    }

}

// ************************************************************************* //
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 06:30
Default
  #5
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Michele, you use a wrong syntax for the relaxation factors in your fvSolution file.
Firstly, "p" must be written under "fields" instead of "equations" as far as I know. Thus, I guess your pressure field is not under-relaxated at all.
Secondly, using the ".*" notation relaxes all iterations, also the last one, as someone stated to one of my posts here:
http://www.cfd-online.com/Forums/ope...tml#post481037
So you should not use the ".*".
I will write some more things if you have posted the rest of the log file.
lpz_michele likes this.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   October 15, 2015, 06:32
Default
  #6
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
I have already posted the rest of the log file.

Thank you for the help.
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 07:11
Default
  #7
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Great, didn't see it.
The pressure solver doesn't work. Could you try to run it with some other setting? If the error remains we need to search somewhere else...
Code:
    "(p|pFinal)"
     {
        solver           GAMG;
        tolerance        1e-12;
        relTol           0.1;

        maxIter         20;    

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

        cacheAgglomeration true;

        nCellsInCoarsestLevel 1000;
        agglomerator     faceAreaPair;
        mergeLevels      1;
     }
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   October 15, 2015, 07:38
Default
  #8
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
I've replaced the code you pasted in the precedent post with the precedent one. And changed all the settings in relaxation factors. My fvSolution file now looks like:
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-12;
        relTol           0.1;

        maxIter         20;    

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

        cacheAgglomeration true;

        nCellsInCoarsestLevel 1000;
        agglomerator     faceAreaPair;
        mergeLevels      1;
     }
     
     
    p
    {
    $pcorr;
    tolerance 1e-07;
    relTol 0.01;
    }  
    
    pFinal
    {
        $p;
    tolerance     1e-07;
        relTol          0;
    }


    U
    {
    solver        smoothSolver;
    smoother    GaussSeidel;
    tolerance    1e-8;
    relTol        0.1;

    }
    
    UFinal
    {
    solver        smoothSolver;
    smoother    GaussSeidel;
    tolerance    1e-8;
    relTol        0.1;

    }
    
    "(k|omega)"
    {
        solver          smoothSolver;
        smoother        GaussSeidel;
        tolerance       1e-7;
        relTol          0.1;
    }

    "(k|omega)Final"
    {
        $U;
    tolerance     1e-7;
        relTol          0;
    }


}

PIMPLE
{
    correctPhi        no;
    nOuterCorrectors    100;
    nCorrectors         2;
    nNonOrthogonalCorrectors 0;
    
    pRefCell         0;
    pRefValue        0;

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

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


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

Trying to launch but it gives me
Code:
 PIMPLE: iteration 1                                                                                                                                                                                                                                                            
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 8.13576e-14, No Iterations 1                                                                                                                                                                             
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 6.4221e-14, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 6.59122e-14, No Iterations 1


--> FOAM FATAL IO ERROR: 
keyword solver is undefined in dictionary "/home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_1_ostacolo_simple_forum/system/fvSolution.solvers.p"

file: /home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_1_ostacolo_simple_forum/system/fvSolution.solvers.p from line 46 to line 47.

    From function dictionary::lookupEntry(const word&, bool, bool) const
    in file db/dictionary/dictionary.C at line 437.

FOAM exiting
Maybe a sytax error? Or just replaced the wrong part?
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 07:43
Default
  #9
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Yes. You need to delete the two old entries "p{ ... }" and "pFinal{ ... }" directly below the new one. Otherwise "my" code is overwritten by your old entries and the solver keyword is missing.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   October 15, 2015, 09:53
Default
  #10
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
Well, this time it gets 7 iterations to create problems

error
Code:
 #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[3] #1  Foam::sigFpe::sigHandler(int) at ??:?
[3] #2   in "/lib64/libc.so.6"
[3] #3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[3] #4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[3] #5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[3] #6  
[3]  at ??:?
[3] #7  
[3]  at ??:?
[3] #8  
[3]  at ??:?
[3] #9  
[3]  at ??:?
[3] #10  __libc_start_main in "/lib64/libc.so.6"
[3] #11  
[3]  at /home/abuild/rpmbuild/BUILD/glibc-2.19/csu/../sysdeps/x86_64/start.S:125
[ime042:32119] *** Process received signal ***
[ime042:32119] Signal: Floating point exception (8)
[ime042:32119] Signal code:  (-6)
[ime042:32119] Failing at address: 0x3e800007d77
[ime042:32119] [ 0] /lib64/libc.so.6(+0x35200) [0x7fe8b18e0200]
[ime042:32119] [ 1] /lib64/libc.so.6(gsignal+0x37) [0x7fe8b18e0187]
[ime042:32119] [ 2] /lib64/libc.so.6(+0x35200) [0x7fe8b18e0200]
[ime042:32119] [ 3] /home/tesisti/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam19GaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x369) [0x7fe8b2a17719]
[ime042:32119] [ 4] /home/tesisti/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam19GaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x2d) [0x7fe8b2a1795d]
[ime042:32119] [ 5] /home/tesisti/OpenFOAM/OpenFOAM-2.3.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x4e7) [0x7fe8b2a11397]
[ime042:32119] [ 6] pimpleFoam() [0x439c5a]
[ime042:32119] [ 7] pimpleFoam() [0x44e78f]
[ime042:32119] [ 8] pimpleFoam() [0x44ea7e]
[ime042:32119] [ 9] pimpleFoam() [0x42020d]
[ime042:32119] [10] /lib64/libc.so.6(__libc_start_main+0xf5) [0x7fe8b18ccb05]
[ime042:32119] [11] pimpleFoam() [0x421071]
[ime042:32119] *** End of error message ***                                                                                                                                                                                                                                    
--------------------------------------------------------------------------                                                                                                                                                                                                     
mpirun noticed that process rank 3 with PID 32119 on node ime042 exited on signal 8 (Floating point exception).                                                                                                                                                                
--------------------------------------------------------------------------
log file
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 15 2015
Time   : 14:50:22
Host   : "ime042"
PID    : 32116
Case   : /home/tesisti/OpenFOAM/tesisti-2.3.1/run/LOPEZ_tent./tent_1_ostacolo_simple_forum
nProcs : 4
Slaves : 
3
(
"ime042.32117"
"ime042.32118"
"ime042.32119"
)

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 = 100
    field "(p|U|k|omega)"    : relTol 0, tolerance 0.01


Starting time loop

forces forces:
    Not including porosity effects

Courant Number mean: 2.72708e-09 max: 2.97915e-07
Time = 1.09083e-05

PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 8.10494e-14, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 6.42711e-14, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 6.55968e-14, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0941324, No Iterations 6
time step continuity errors : sum local = 5.13413e-10, global = -6.965e-12, cumulative = -6.965e-12
GAMG:  Solving for p, Initial residual = 0.18443, Final residual = 0.0162853, No Iterations 2
time step continuity errors : sum local = 4.27511e-09, global = -4.4206e-10, cumulative = -4.49025e-10
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.551614, Final residual = 0.000245009, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.117547, Final residual = 6.71312e-05, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.130338, Final residual = 3.37683e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.155643, Final residual = 0.0111007, No Iterations 3
time step continuity errors : sum local = 2.67407e-09, global = -1.67831e-10, cumulative = -6.16856e-10
GAMG:  Solving for p, Initial residual = 0.131887, Final residual = 0.00901081, No Iterations 2
time step continuity errors : sum local = 3.61971e-09, global = -3.24685e-10, cumulative = -9.41541e-10
PIMPLE: iteration 3
smoothSolver:  Solving for Ux, Initial residual = 0.293694, Final residual = 0.000130662, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0586409, Final residual = 3.56826e-05, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.0660638, Final residual = 1.2281e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.119831, Final residual = 0.00790474, No Iterations 2
time step continuity errors : sum local = 3.18576e-09, global = -2.72153e-10, cumulative = -1.21369e-09
GAMG:  Solving for p, Initial residual = 0.0971771, Final residual = 0.00894128, No Iterations 1
time step continuity errors : sum local = 4.45609e-09, global = -5.33984e-10, cumulative = -1.74768e-09
PIMPLE: iteration 4
smoothSolver:  Solving for Ux, Initial residual = 0.0963746, Final residual = 3.0991e-05, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0294886, Final residual = 1.47947e-05, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.0324897, Final residual = 4.06102e-06, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0889663, Final residual = 0.00849765, No Iterations 1
time step continuity errors : sum local = 4.15247e-09, global = -4.6771e-10, cumulative = -2.21539e-09
GAMG:  Solving for p, Initial residual = 0.0697397, Final residual = 0.00692439, No Iterations 1
time step continuity errors : sum local = 3.73189e-09, global = -4.04551e-10, cumulative = -2.61994e-09
PIMPLE: iteration 5
smoothSolver:  Solving for Ux, Initial residual = 0.0491208, Final residual = 3.19173e-05, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.0162711, Final residual = 1.53472e-05, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.0181073, Final residual = 4.73704e-06, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0667263, Final residual = 0.00600909, No Iterations 1
time step continuity errors : sum local = 3.28066e-09, global = -3.66243e-10, cumulative = -2.98618e-09
GAMG:  Solving for p, Initial residual = 0.0517721, Final residual = 0.00505175, No Iterations 1
time step continuity errors : sum local = 2.93846e-09, global = -3.62277e-10, cumulative = -3.34846e-09
PIMPLE: iteration 6
smoothSolver:  Solving for Ux, Initial residual = 0.015568, Final residual = 0.000222185, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00836338, Final residual = 4.93298e-05, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.0108277, Final residual = 7.2597e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0534809, Final residual = 0.00328122, No Iterations 3
time step continuity errors : sum local = 1.92493e-09, global = 6.20145e-11, cumulative = -3.28644e-09
GAMG:  Solving for p, Initial residual = 0.0597616, Final residual = 0.00585955, No Iterations 1
time step continuity errors : sum local = 3.77776e-09, global = -2.90223e-11, cumulative = -3.31547e-09
PIMPLE: iteration 7
Error must be somewhere else.
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 09:57
Default
  #11
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Looks like it's the velocity solver this time, right?
I am not a solver expert, so I can't tell you why the smooth solver creates these problems. Usually, I use this setting:


Code:
solvers
{
        "(p|pFinal)"
     {
        solver           GAMG;
        tolerance        1e-12;
        relTol           0.1;

        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        1e-12;
        relTol           1.0e-2;
        maxIter         100;
        };

}
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   October 15, 2015, 09:59
Default
  #12
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
Ok, so actually my fvSolution looks like

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-12;
        relTol           0.1;

        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        1e-12;
        relTol           1.0e-2;
        maxIter         100;
        };
  


}

PIMPLE
{
    correctPhi        no;
    nOuterCorrectors    100;
    nCorrectors         2;
    nNonOrthogonalCorrectors 0;
    
    pRefCell         0;
    pRefValue        0;

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

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


// ************************************************************************* //
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 10:02
Default
  #13
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
I would recommend a different residualControl setting. I usually use the relTol rather than the tolerance, because you don't have to concern yourself with the absolute values of the residuals. But maybe others think differently about that !?
Then, I would also lower the nOuterCorrectors to something like 20. If you can't get residuals down in 20 iterations, 100 won't work, too.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   October 15, 2015, 10:10
Default
  #14
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
Well, now it keeps going on.

First time step converged in 12 iterations. But numbers don't match with physics.

Second time step starts with a max Courant number of 76, further more time step continuity errors keep growing.

Code:
 
forces forces output:
    sum of forces:
        pressure : (3.20219e+13 6.22034e+11 -4.0116e+11)
        viscous  : (226981 3372.33 7927.69)
        porous   : (0 0 0)
    sum of moments:
        pressure : (-1.15048e+14 4.54921e+13 -3.75908e+14)
        viscous  : (-1.4342e+06 -2.72523e+06 -1.60083e+07)
        porous   : (0 0 0)

Courant Number mean: 4.59332e-07 max: 76.5398
Time = 2.18166e-05

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.866664, Final residual = 3.15156e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.257985, Final residual = 6.31871e-06, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.25924, Final residual = 4.27647e-07, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0638671, Final residual = 0.00387404, No Iterations 3
time step continuity errors : sum local = 2.82518e-09, global = 2.29057e-10, cumulative = -4.71539e-09
GAMG:  Solving for p, Initial residual = 0.0570386, Final residual = 0.00406766, No Iterations 3
time step continuity errors : sum local = 2.34492e-09, global = 2.2376e-10, cumulative = -4.49163e-09
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 0.260287, Final residual = 2.3509e-05, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0472878, Final residual = 6.97051e-06, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.0472358, Final residual = 6.89489e-06, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0639765, Final residual = 0.00506699, No Iterations 3
time step continuity errors : sum local = 2.88782e-09, global = 2.81694e-10, cumulative = -4.20994e-09
GAMG:  Solving for p, Initial residual = 0.0966392, Final residual = 0.00895581, No Iterations 2
time step continuity errors : sum local = 3.7894e-09, global = 4.5091e-10, cumulative = -3.75903e-09
PIMPLE: iteration 3
DILUPBiCG:  Solving for Ux, Initial residual = 0.0180231, Final residual = 4.31519e-05, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.00996661, Final residual = 1.28832e-05, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.00989444, Final residual = 3.77719e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0819241, Final residual = 0.00574263, No Iterations 3
time step continuity errors : sum local = 2.42303e-09, global = 2.66958e-10, cumulative = -3.49207e-09
GAMG:  Solving for p, Initial residual = 0.132618, Final residual = 0.0101418, No Iterations 2
time step continuity errors : sum local = 3.15306e-09, global = 4.19991e-10, cumulative = -3.07208e-09
PIMPLE: iteration 4
DILUPBiCG:  Solving for Ux, Initial residual = 0.0849585, Final residual = 0.000508934, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0167581, Final residual = 1.4555e-06, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.0169335, Final residual = 8.71885e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.115883, Final residual = 0.0115837, No Iterations 1
time step continuity errors : sum local = 3.59584e-09, global = 1.59479e-09, cumulative = -1.47728e-09
GAMG:  Solving for p, Initial residual = 0.0900661, Final residual = 0.00743689, No Iterations 3
time step continuity errors : sum local = 2.21426e-09, global = 2.69514e-10, cumulative = -1.20777e-09
PIMPLE: iteration 5
DILUPBiCG:  Solving for Ux, Initial residual = 0.078644, Final residual = 0.000613499, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.0161631, Final residual = 0.00011581, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.016418, Final residual = 2.26945e-05, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.163533, Final residual = 0.015967, No Iterations 1
time step continuity errors : sum local = 3.72005e-09, global = 1.51401e-09, cumulative = 3.06234e-10
GAMG:  Solving for p, Initial residual = 0.0970837, Final residual = 0.00890793, No Iterations 3
time step continuity errors : sum local = 2.17655e-09, global = 3.06431e-10, cumulative = 6.12665e-10
PIMPLE: iteration 6
DILUPBiCG:  Solving for Ux, Initial residual = 0.0531039, Final residual = 2.68521e-05, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.00897067, Final residual = 7.4357e-05, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.0113846, Final residual = 2.54919e-05, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.18581, Final residual = 0.0142419, No Iterations 1
time step continuity errors : sum local = 2.69513e-09, global = 9.0678e-10, cumulative = 1.51945e-09
GAMG:  Solving for p, Initial residual = 0.122813, Final residual = 0.0107259, No Iterations 2
time step continuity errors : sum local = 2.05469e-09, global = 5.09273e-10, cumulative = 2.02872e-09
PIMPLE: iteration 7
DILUPBiCG:  Solving for Ux, Initial residual = 0.0142203, Final residual = 1.11577e-05, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.00513482, Final residual = 7.10856e-06, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.00550852, Final residual = 5.4645e-05, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.155244, Final residual = 0.0127537, No Iterations 1
time step continuity errors : sum local = 2.21653e-09, global = 7.59961e-10, cumulative = 2.78868e-09
GAMG:  Solving for p, Initial residual = 0.105718, Final residual = 0.00990871, No Iterations 2
time step continuity errors : sum local = 1.74794e-09, global = 3.39822e-10, cumulative = 3.1285e-09
PIMPLE: iteration 8
DILUPBiCG:  Solving for Ux, Initial residual = 0.00490551, Final residual = 1.22639e-05, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.00136004, Final residual = 7.06376e-07, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.00267391, Final residual = 5.58931e-06, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.142433, Final residual = 0.00811321, No Iterations 3
time step continuity errors : sum local = 1.27769e-09, global = 4.30803e-12, cumulative = 3.13281e-09
GAMG:  Solving for p, Initial residual = 0.169433, Final residual = 0.016526, No Iterations 2
time step continuity errors : sum local = 2.1341e-09, global = 2.21144e-10, cumulative = 3.35395e-09
PIMPLE: iteration 9
DILUPBiCG:  Solving for Ux, Initial residual = 0.00287841, Final residual = 1.57161e-05, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.000483634, Final residual = 4.30534e-06, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.00230007, Final residual = 1.06611e-06, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.157906, Final residual = 0.0155495, No Iterations 3
time step continuity errors : sum local = 2.3094e-09, global = 5.37005e-10, cumulative = 3.89096e-09
GAMG:  Solving for p, Initial residual = 0.141294, Final residual = 0.0131548, No Iterations 4
time step continuity errors : sum local = 2.03133e-09, global = 4.90192e-10, cumulative = 4.38115e-09
PIMPLE: iteration 10
DILUPBiCG:  Solving for Ux, Initial residual = 0.000192267, Final residual = 2.26198e-07, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.000125139, Final residual = 1.72955e-07, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.000595267, Final residual = 2.00749e-06, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.165359, Final residual = 0.0152966, No Iterations 4
time step continuity errors : sum local = 2.71149e-09, global = 3.18975e-11, cumulative = 4.41305e-09
GAMG:  Solving for p, Initial residual = 0.300512, Final residual = 0.0231589, No Iterations 3
time step continuity errors : sum local = 4.86893e-09, global = -7.42521e-11, cumulative = 4.3388e-09
PIMPLE: iteration 11
DILUPBiCG:  Solving for Ux, Initial residual = 0.000350208, Final residual = 1.94364e-06, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.00032242, Final residual = 2.32274e-06, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.00115122, Final residual = 2.25507e-06, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.348331, Final residual = 0.0293308, No Iterations 4
time step continuity errors : sum local = 7.32133e-09, global = -2.00533e-10, cumulative = 4.13826e-09
GAMG:  Solving for p, Initial residual = 0.17185, Final residual = 0.0124721, No Iterations 3
time step continuity errors : sum local = 8.7529e-09, global = -1.76329e-10, cumulative = 3.96193e-09
PIMPLE: iteration 12
DILUPBiCG:  Solving for Ux, Initial residual = 0.000558672, Final residual = 4.85112e-06, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.000562371, Final residual = 5.38573e-06, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.00113069, Final residual = 6.76897e-06, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.386677, Final residual = 0.034756, No Iterations 5
time step continuity errors : sum local = 1.65232e-08, global = 1.55496e-10, cumulative = 4.11743e-09
GAMG:  Solving for p, Initial residual = 0.0460313, Final residual = 0.00430096, No Iterations 2
time step continuity errors : sum local = 5.36873e-08, global = -5.97013e-10, cumulative = 3.52042e-09

...

PIMPLE: iteration 25
DILUPBiCG:  Solving for Ux, Initial residual = 0.45703, Final residual = 0.00364494, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.466817, Final residual = 0.0032379, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.448066, Final residual = 0.00264534, No Iterations 4
GAMG:  Solving for p, Initial residual = 4.89754e-07, Final residual = 4.45106e-08, No Iterations 17
time step continuity errors : sum local = 0.0126292, global = -0.000375018, cumulative = -0.000298839
GAMG:  Solving for p, Initial residual = 1.2146e-07, Final residual = 7.14817e-08, No Iterations 20
time step continuity errors : sum local = 0.128873, global = 0.00011411, cumulative = -0.00018473
PIMPLE: iteration 26
DILUPBiCG:  Solving for Ux, Initial residual = 0.488462, Final residual = 0.00410571, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.489245, Final residual = 0.00398358, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.468887, Final residual = 0.00105651, No Iterations 4
GAMG:  Solving for p, Initial residual = 1.55761e-07, Final residual = 1.28471e-08, No Iterations 6
time step continuity errors : sum local = 0.0261039, global = 0.000406571, cumulative = 0.000221841
GAMG:  Solving for p, Initial residual = 4.17219e-08, Final residual = 3.32099e-09, No Iterations 12
time step continuity errors : sum local = 0.0568839, global = -0.00190224, cumulative = -0.0016804
PIMPLE: iteration 27
DILUPBiCG:  Solving for Ux, Initial residual = 0.485182, Final residual = 0.00444568, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.477219, Final residual = 0.00236323, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.452158, Final residual = 0.000454943, No Iterations 5
GAMG:  Solving for p, Initial residual = 9.29377e-08, Final residual = 7.82375e-09, No Iterations 7
time step continuity errors : sum local = 0.08866, global = 0.00500549, cumulative = 0.00332509
GAMG:  Solving for p, Initial residual = 6.04943e-08, Final residual = 5.50121e-09, No Iterations 2
time step continuity errors : sum local = 0.142202, global = 0.00511167, cumulative = 0.00843676
PIMPLE: iteration 28
DILUPBiCG:  Solving for Ux, Initial residual = 0.5799, Final residual = 0.0032216, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.56742, Final residual = 0.00354331, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.464652, Final residual = 0.00279237, No Iterations 3
GAMG:  Solving for p, Initial residual = 1.15256e-07, Final residual = 1.73304e-08, No Iterations 20
time step continuity errors : sum local = 0.296138, global = -0.000460135, cumulative = 0.00797662
GAMG:  Solving for p, Initial residual = 1.43725e-08, Final residual = 1.9104e-08, No Iterations 20
time step continuity errors : sum local = 3.256, global = -4.9317e-05, cumulative = 0.0079273
PIMPLE: iteration 29
DILUPBiCG:  Solving for Ux, Initial residual = 0.475133, Final residual = 0.00420981, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.514787, Final residual = 0.00477123, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.387523, Final residual = 0.00241095, No Iterations 4
GAMG:  Solving for p, Initial residual = 7.66632e-08, Final residual = 4.95751e-09, No Iterations 5
time step continuity errors : sum local = 0.351647, global = 0.00446064, cumulative = 0.0123879
GAMG:  Solving for p, Initial residual = 1.11533e-08, Final residual = 8.46322e-10, No Iterations 4
time step continuity errors : sum local = 0.669685, global = -0.0141515, cumulative = -0.00176359
PIMPLE: iteration 30
DILUPBiCG:  Solving for Ux, Initial residual = 0.482939, Final residual = 0.00169732, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.492853, Final residual = 0.00167008, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.423742, Final residual = 0.0033366, No Iterations 3
GAMG:  Solving for p, Initial residual = 7.33705e-08, Final residual = 7.02836e-09, No Iterations 4
time step continuity errors : sum local = 1.208, global = 0.0817156, cumulative = 0.079952
GAMG:  Solving for p, Initial residual = 9.10694e-09, Final residual = 7.2643e-10, No Iterations 8
time step continuity errors : sum local = 1.11649, global = 0.0696683, cumulative = 0.14962
PIMPLE: iteration 31
DILUPBiCG:  Solving for Ux, Initial residual = 0.434401, Final residual = 0.00276104, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.438609, Final residual = 0.0035577, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.428942, Final residual = 0.00169325, No Iterations 5
GAMG:  Solving for p, Initial residual = 6.13821e-06, Final residual = 8.82587e-07, No Iterations 20
time step continuity errors : sum local = 3.98246, global = 0.277155, cumulative = 0.426775
GAMG:  Solving for p, Initial residual = 2.69328e-08, Final residual = 2.1687e-09, No Iterations 1
time step continuity errors : sum local = 4.15637, global = 0.0797123, cumulative = 0.506488
PIMPLE: iteration 32
DILUPBiCG:  Solving for Ux, Initial residual = 0.728414, Final residual = 0.00296967, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.541076, Final residual = 0.00394594, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.79722, Final residual = 0.00145519, No Iterations 3
GAMG:  Solving for p, Initial residual = 7.0245e-08, Final residual = 6.54451e-09, No Iterations 9
time step continuity errors : sum local = 5.78199, global = 0.245586, cumulative = 0.752073
GAMG:  Solving for p, Initial residual = 5.14576e-09, Final residual = 3.40823e-10, No Iterations 8
time step continuity errors : sum local = 8.29289, global = -0.468393, cumulative = 0.283681
PIMPLE: iteration 33
DILUPBiCG:  Solving for Ux, Initial residual = 0.431882, Final residual = 0.000428693, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0.451638, Final residual = 0.00125437, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.500052, Final residual = 0.00311681, No Iterations 4
GAMG:  Solving for p, Initial residual = 9.34601e-08, Final residual = 1.17991e-08, No Iterations 20
time step continuity errors : sum local = 23.4012, global = -15.7156, cumulative = -15.4319
GAMG:  Solving for p, Initial residual = 4.27502e-08, Final residual = 1.25665e-08, No Iterations 20
time step continuity errors : sum local = 85.7935, global = -34.1987, cumulative = -49.6306
PIMPLE: iteration 34
DILUPBiCG:  Solving for Ux, Initial residual = 0.40978, Final residual = 0.00235957, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.451869, Final residual = 0.001927, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.419105, Final residual = 0.00210072, No Iterations 5
GAMG:  Solving for p, Initial residual = 1.92655e-06, Final residual = 8.91561e-07, No Iterations 20
time step continuity errors : sum local = 208.798, global = -51.2509, cumulative = -100.881
GAMG:  Solving for p, Initial residual = 7.43198e-07, Final residual = 2.23723e-06, No Iterations 20
time step continuity errors : sum local = 1943.26, global = 18.5426, cumulative = -82.3388
PIMPLE: iteration 35
DILUPBiCG:  Solving for Ux, Initial residual = 0.330767, Final residual = 0.0021958, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.343671, Final residual = 0.00117926, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.322184, Final residual = 0.00102286, No Iterations 4
GAMG:  Solving for p, Initial residual = 2.22497e-07, Final residual = 2.12035e-08, No Iterations 6
time step continuity errors : sum local = 118.654, global = -61.5736, cumulative = -143.912
GAMG:  Solving for p, Initial residual = 7.57818e-08, Final residual = 6.6613e-09, No Iterations 10
time step continuity errors : sum local = 186.871, global = 43.3242, cumulative = -100.588
PIMPLE: iteration 36
DILUPBiCG:  Solving for Ux, Initial residual = 0.557896, Final residual = 0.00349763, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.566518, Final residual = 0.000803787, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.484968, Final residual = 0.00214561, No Iterations 4
GAMG:  Solving for p, Initial residual = 3.62432e-07, Final residual = 3.45964e-08, No Iterations 13
time step continuity errors : sum local = 388.037, global = -4.52998, cumulative = -105.118
GAMG:  Solving for p, Initial residual = 5.9511e-08, Final residual = 1.23765e-08, No Iterations 20
time step continuity errors : sum local = 1321.47, global = -0.315589, cumulative = -105.434
PIMPLE: iteration 37
DILUPBiCG:  Solving for Ux, Initial residual = 0.486504, Final residual = 0.00191651, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.48672, Final residual = 0.00204066, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.427991, Final residual = 0.000811272, No Iterations 5
GAMG:  Solving for p, Initial residual = 4.65018e-08, Final residual = 4.45178e-09, No Iterations 20
time step continuity errors : sum local = 731.468, global = 198.436, cumulative = 93.0026
GAMG:  Solving for p, Initial residual = 9.32307e-09, Final residual = 4.49918e-09, No Iterations 20
time step continuity errors : sum local = 4870.59, global = -1325.02, cumulative = -1232.01
PIMPLE: iteration 38
DILUPBiCG:  Solving for Ux, Initial residual = 0.479185, Final residual = 0.00165469, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.476365, Final residual = 0.000516798, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.489416, Final residual = 0.00472329, No Iterations 3
GAMG:  Solving for p, Initial residual = 1.9266e-06, Final residual = 5.08132e-07, No Iterations 20
time step continuity errors : sum local = 7383.56, global = 3325.75, cumulative = 2093.74
GAMG:  Solving for p, Initial residual = 4.86745e-07, Final residual = 4.49227e-07, No Iterations 20
time step continuity errors : sum local = 57164.3, global = -1297.7, cumulative = 796.039
PIMPLE: iteration 39
DILUPBiCG:  Solving for Ux, Initial residual = 0.469658, Final residual = 0.00286888, No Iterations 4
DILUPBiCG:  Solving for Uy, Initial residual = 0.486166, Final residual = 0.000287571, No Iterations 6
DILUPBiCG:  Solving for Uz, Initial residual = 0.420831, Final residual = 0.00308034, No Iterations 3
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 10:11
Default
  #15
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
Trying different settings like you suggested. I'll write on it.
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 10:14
Default
  #16
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Don't cut the log file, I miss the most intersting parts
What about the k and omaga solvers?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   October 15, 2015, 10:21
Default
  #17
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
this one, as you suggested, right?

Code:
 "(U|k|epsilon|omega|UFinal|kFinal|epsilonFinal|omegaFinal)"    
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        1e-12;
        relTol           1.0e-2;
        maxIter         100;
        };
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 10:23
Default
  #18
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
No, I mean the log file is missing the information about the turbulence.
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Old   October 15, 2015, 10:27
Default
  #19
Member
 
Mike Lob
Join Date: Oct 2015
Location: It
Posts: 51
Rep Power: 10
lpz_michele is on a distinguished road
As far as I remember, all the log files I've seen did not write nothing about k-omega. Did I forgot to write something?


Anyway, I tried with the following 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-5;

        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-5;
        maxIter         100;
        };
  


}

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

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

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


// ************************************************************************* //
log file
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 1e-05, tolerance 0


Starting time loop

forces forces:
    Not including porosity effects

Courant Number mean: 2.72708e-09 max: 2.97915e-07
Time = 1.09083e-05

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 5.38097e-17, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 2.19952e-16, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 8.98804e-17, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.000850538, No Iterations 20
time step continuity errors : sum local = 4.63897e-12, global = -2.99077e-14, cumulative = -2.99077e-14
GAMG:  Solving for p, Initial residual = 0.182951, Final residual = 1.48404e-05, No Iterations 20
time step continuity errors : sum local = 3.92995e-12, global = -4.25516e-14, cumulative = -7.24593e-14
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 0.575869, Final residual = 5.20509e-09, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.130436, Final residual = 5.97277e-10, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.130436, Final residual = 6.17317e-10, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.14681, Final residual = 1.26348e-05, No Iterations 20
time step continuity errors : sum local = 3.58e-12, global = -4.61894e-14, cumulative = -1.18649e-13
GAMG:  Solving for p, Initial residual = 0.123814, Final residual = 6.48693e-06, No Iterations 20
time step continuity errors : sum local = 2.91471e-12, global = -4.44487e-14, cumulative = -1.63097e-13
PIMPLE: iteration 3
DILUPBiCG:  Solving for Ux, Initial residual = 0.328115, Final residual = 2.9441e-09, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0652315, Final residual = 2.92304e-10, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0652764, Final residual = 3.10868e-10, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.113768, Final residual = 5.22568e-06, No Iterations 20
time step continuity errors : sum local = 2.41215e-12, global = -3.75052e-14, cumulative = -2.00603e-13
GAMG:  Solving for p, Initial residual = 0.091437, Final residual = 3.3779e-06, No Iterations 20
time step continuity errors : sum local = 1.89695e-12, global = -3.15321e-14, cumulative = -2.32135e-13
PIMPLE: iteration 4
DILUPBiCG:  Solving for Ux, Initial residual = 0.17629, Final residual = 1.56232e-09, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0338782, Final residual = 1.51358e-10, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0339012, Final residual = 1.6177e-10, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0863256, Final residual = 2.64111e-06, No Iterations 20
time step continuity errors : sum local = 1.50007e-12, global = -2.45599e-14, cumulative = -2.56695e-13
GAMG:  Solving for p, Initial residual = 0.0670754, Final residual = 1.84298e-06, No Iterations 20
time step continuity errors : sum local = 1.15507e-12, global = -1.90806e-14, cumulative = -2.75775e-13
PIMPLE: iteration 5
DILUPBiCG:  Solving for Ux, Initial residual = 0.0931013, Final residual = 8.08966e-10, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0182765, Final residual = 7.92359e-11, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0182847, Final residual = 8.65505e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0635045, Final residual = 1.43359e-06, No Iterations 20
time step continuity errors : sum local = 9.02351e-13, global = -1.44497e-14, cumulative = -2.90225e-13
GAMG:  Solving for p, Initial residual = 0.048217, Final residual = 1.04543e-06, No Iterations 20
time step continuity errors : sum local = 6.94057e-13, global = -1.09189e-14, cumulative = -3.01144e-13
PIMPLE: iteration 6
DILUPBiCG:  Solving for Ux, Initial residual = 0.048723, Final residual = 4.10896e-10, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0105615, Final residual = 4.24163e-11, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0105538, Final residual = 4.84504e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0454965, Final residual = 8.2809e-07, No Iterations 20
time step continuity errors : sum local = 5.50428e-13, global = -8.40717e-15, cumulative = -3.09551e-13
GAMG:  Solving for p, Initial residual = 0.0339978, Final residual = 6.22983e-07, No Iterations 20
time step continuity errors : sum local = 4.26835e-13, global = -6.33486e-15, cumulative = -3.15886e-13
PIMPLE: iteration 7
DILUPBiCG:  Solving for Ux, Initial residual = 0.025365, Final residual = 2.14062e-10, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00664658, Final residual = 2.96478e-11, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00662996, Final residual = 2.91403e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0319445, Final residual = 5.04916e-07, No Iterations 20
time step continuity errors : sum local = 3.45932e-13, global = -4.91469e-15, cumulative = -3.208e-13
GAMG:  Solving for p, Initial residual = 0.0236159, Final residual = 3.87828e-07, No Iterations 20
time step continuity errors : sum local = 2.70513e-13, global = -3.69509e-15, cumulative = -3.24496e-13
PIMPLE: iteration 8
DILUPBiCG:  Solving for Ux, Initial residual = 0.0132529, Final residual = 3.12824e-10, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00443574, Final residual = 1.94406e-11, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.00442741, Final residual = 2.19916e-11, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0221197, Final residual = 3.17427e-07, No Iterations 20
time step continuity errors : sum local = 2.21317e-13, global = -2.41478e-15, cumulative = -3.2691e-13
GAMG:  Solving for p, Initial residual = 0.016238, Final residual = 2.47359e-07, No Iterations 20
time step continuity errors : sum local = 1.74362e-13, global = -1.75134e-15, cumulative = -3.28662e-13
PIMPLE: iteration 9
DILUPBiCG:  Solving for Ux, Initial residual = 0.00708739, Final residual = 1.48741e-10, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.00301968, Final residual = 8.01464e-09, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0030182, Final residual = 1.53244e-10, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0151832, Final residual = 2.32595e-07, No Iterations 20
time step continuity errors : sum local = 1.6388e-13, global = -6.4121e-16, cumulative = -3.29303e-13
GAMG:  Solving for p, Initial residual = 0.0110994, Final residual = 1.86506e-07, No Iterations 20
time step continuity errors : sum local = 1.32314e-13, global = -3.65334e-16, cumulative = -3.29668e-13
PIMPLE: iteration 10
DILUPBiCG:  Solving for Ux, Initial residual = 0.00396473, Final residual = 9.72602e-12, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.00205423, Final residual = 3.17736e-11, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.00205473, Final residual = 2.35815e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.0103786, Final residual = 2.70309e-07, No Iterations 20
time step continuity errors : sum local = 1.91695e-13, global = 6.29131e-17, cumulative = -3.29605e-13
GAMG:  Solving for p, Initial residual = 0.00757311, Final residual = 2.43872e-07, No Iterations 20
time step continuity errors : sum local = 1.73711e-13, global = 1.07779e-16, cumulative = -3.29498e-13
PIMPLE: iteration 11
DILUPBiCG:  Solving for Ux, Initial residual = 0.00234174, Final residual = 2.28446e-09, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 0.00138961, Final residual = 1.31493e-09, No Iterations 4
DILUPBiCG:  Solving for Uz, Initial residual = 0.00139063, Final residual = 2.06147e-09, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.00709553, Final residual = 5.57821e-07, No Iterations 20
time step continuity errors : sum local = 3.97233e-13, global = 2.16375e-16, cumulative = -3.29281e-13
GAMG:  Solving for p, Initial residual = 0.00518617, Final residual = 5.17207e-07, No Iterations 20
time step continuity errors : sum local = 3.69376e-13, global = 1.81566e-16, cumulative = -3.291e-13
PIMPLE: iteration 12
DILUPBiCG:  Solving for Ux, Initial residual = 0.001446, Final residual = 1.07055e-08, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.000934765, Final residual = 2.94952e-09, No Iterations 6
DILUPBiCG:  Solving for Uz, Initial residual = 0.000936384, Final residual = 1.07785e-09, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00488315, Final residual = 1.22038e-06, No Iterations 20
time step continuity errors : sum local = 8.71415e-13, global = 3.63095e-17, cumulative = -3.29063e-13
GAMG:  Solving for p, Initial residual = 0.00361049, Final residual = 1.10981e-06, No Iterations 20
time step continuity errors : sum local = 7.93999e-13, global = 2.2947e-17, cumulative = -3.2904e-13
PIMPLE: iteration 13
DILUPBiCG:  Solving for Ux, Initial residual = 0.000913855, Final residual = 8.92874e-10, No Iterations 7
DILUPBiCG:  Solving for Uy, Initial residual = 0.000624009, Final residual = 6.93498e-10, No Iterations 6
DILUPBiCG:  Solving for Uz, Initial residual = 0.000628673, Final residual = 2.40093e-09, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.00344592, Final residual = 2.92139e-06, No Iterations 20
time step continuity errors : sum local = 2.08991e-12, global = 1.37988e-17, cumulative = -3.29027e-13
GAMG:  Solving for p, Initial residual = 0.00265492, Final residual = 2.466e-06, No Iterations 20
time step continuity errors : sum local = 1.76652e-12, global = -4.18087e-18, cumulative = -3.29031e-13
PIMPLE: iteration 14
DILUPBiCG:  Solving for Ux, Initial residual = 0.000574472, Final residual = 2.83792e-09, No Iterations 7
DILUPBiCG:  Solving for Uy, Initial residual = 0.000391953, Final residual = 8.47879e-10, No Iterations 6
DILUPBiCG:  Solving for Uz, Initial residual = 0.000421002, Final residual = 1.50038e-09, No Iterations 8
GAMG:  Solving for p, Initial residual = 0.00263674, Final residual = 5.59006e-06, No Iterations 20
time step continuity errors : sum local = 4.00474e-12, global = 2.33869e-17, cumulative = -3.29007e-13
GAMG:  Solving for p, Initial residual = 0.00233292, Final residual = 3.80849e-06, No Iterations 20
time step continuity errors : sum local = 2.73149e-12, global = -1.58426e-17, cumulative = -3.29023e-13
PIMPLE: iteration 15
DILUPBiCG:  Solving for Ux, Initial residual = 0.000263722, Final residual = 4.01094e-10, No Iterations 8
DILUPBiCG:  Solving for Uy, Initial residual = 0.000157811, Final residual = 1.16427e-09, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 0.000278886, Final residual = 2.53124e-09, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.00328214, Final residual = 1.48079e-05, No Iterations 20
time step continuity errors : sum local = 1.06313e-11, global = 1.10572e-16, cumulative = -3.28913e-13
GAMG:  Solving for p, Initial residual = 0.00393972, Final residual = 1.76433e-05, No Iterations 20
time step continuity errors : sum local = 1.26889e-11, global = -9.98943e-17, cumulative = -3.29012e-13
PIMPLE: iteration 16
DILUPBiCG:  Solving for Ux, Initial residual = 3.70449e-05, Final residual = 3.67244e-10, No Iterations 8
DILUPBiCG:  Solving for Uy, Initial residual = 2.11178e-05, Final residual = 1.22973e-10, No Iterations 9
DILUPBiCG:  Solving for Uz, Initial residual = 0.0001729, Final residual = 8.05336e-10, No Iterations 8
GAMG:  Solving for p, Initial residual = 0.00536874, Final residual = 9.37444e-05, No Iterations 20
time step continuity errors : sum local = 6.75098e-11, global = -7.25852e-16, cumulative = -3.29738e-13
GAMG:  Solving for p, Initial residual = 0.0132991, Final residual = 0.000106057, No Iterations 20
time step continuity errors : sum local = 7.74716e-11, global = 9.05626e-16, cumulative = -3.28833e-13
PIMPLE: iteration 17
DILUPBiCG:  Solving for Ux, Initial residual = 0.000318289, Final residual = 3.27937e-10, No Iterations 12
DILUPBiCG:  Solving for Uy, Initial residual = 0.000168046, Final residual = 1.47304e-09, No Iterations 8
DILUPBiCG:  Solving for Uz, Initial residual = 2.59715e-05, Final residual = 8.5052e-11, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.0187933, Final residual = 0.00361403, No Iterations 20
time step continuity errors : sum local = 2.62536e-09, global = -2.1167e-13, cumulative = -5.40503e-13
GAMG:  Solving for p, Initial residual = 0.0783631, Final residual = 0.00585178, No Iterations 20
time step continuity errors : sum local = 5.31193e-09, global = -4.36381e-13, cumulative = -9.76884e-13
PIMPLE: iteration 18
DILUPBiCG:  Solving for Ux, Initial residual = 0.000485379, Final residual = 1.97349e-09, No Iterations 10
DILUPBiCG:  Solving for Uy, Initial residual = 0.000632127, Final residual = 2.95446e-09, No Iterations 10
DILUPBiCG:  Solving for Uz, Initial residual = 0.000310312, Final residual = 1.76645e-09, No Iterations 9
GAMG:  Solving for p, Initial residual = 0.0984537, Final residual = 0.00645455, No Iterations 20
time step continuity errors : sum local = 4.64713e-09, global = 3.52369e-14, cumulative = -9.41647e-13
GAMG:  Solving for p, Initial residual = 0.0469483, Final residual = 0.00256194, No Iterations 20
time step continuity errors : sum local = 6.73626e-09, global = -8.1758e-14, cumulative = -1.02341e-12
PIMPLE: iteration 19
DILUPBiCG:  Solving for Ux, Initial residual = 0.000620827, Final residual = 1.78317e-09, No Iterations 11
DILUPBiCG:  Solving for Uy, Initial residual = 0.000546878, Final residual = 1.73103e-09, No Iterations 10
DILUPBiCG:  Solving for Uz, Initial residual = 0.000262496, Final residual = 6.11569e-10, No Iterations 11
GAMG:  Solving for p, Initial residual = 0.158795, Final residual = 0.0197291, No Iterations 20
time step continuity errors : sum local = 3.89702e-08, global = -5.62443e-13, cumulative = -1.58585e-12
GAMG:  Solving for p, Initial residual = 0.0769585, Final residual = 0.00629057, No Iterations 20
time step continuity errors : sum local = 4.59415e-08, global = -5.9324e-13, cumulative = -2.17909e-12
PIMPLE: iteration 20
DILUPBiCG:  Solving for Ux, Initial residual = 0.0033839, Final residual = 2.3056e-08, No Iterations 65
DILUPBiCG:  Solving for Uy, Initial residual = 0.00399156, Final residual = 2.71118e-08, No Iterations 56
DILUPBiCG:  Solving for Uz, Initial residual = 0.00267406, Final residual = 1.59844e-08, No Iterations 65
GAMG:  Solving for p, Initial residual = 0.191895, Final residual = 0.00318539, No Iterations 20
time step continuity errors : sum local = 3.28746e-08, global = -6.25726e-12, cumulative = -8.43635e-12
GAMG:  Solving for p, Initial residual = 0.0450605, Final residual = 0.00172445, No Iterations 20
time step continuity errors : sum local = 3.66081e-07, global = -5.65159e-11, cumulative = -6.49523e-11
DILUPBiCG:  Solving for omega, Initial residual = 0.0624003, Final residual = 4.45631e-07, No Iterations 49
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 8.95703e-06, No Iterations 34
PIMPLE: not converged within 20 iterations
ExecutionTime = 374.04 s  ClockTime = 374 s

forces forces output:
    sum of forces:
        pressure : (-2.28408e+17 5.32264e+17 6.93677e+17)
        viscous  : (-2.95957e+09 9.28817e+07 -1.33636e+10)
        porous   : (0 0 0)
    sum of moments:
        pressure : (4.85726e+20 2.13491e+20 1.42405e+19)
        viscous  : (2.56719e+11 2.87508e+11 1.30512e+11)
        porous   : (0 0 0)

Courant Number mean: 5.74294e-06 max: 1.27126e+06
Time = 2.18166e-05

PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.0575416, Final residual = 2.68272e-07, No Iterations 10
DILUPBiCG:  Solving for Uy, Initial residual = 0.0395307, Final residual = 2.18704e-07, No Iterations 10
DILUPBiCG:  Solving for Uz, Initial residual = 0.0887286, Final residual = 6.59463e-07, No Iterations 11
GAMG:  Solving for p, Initial residual = 0.0210242, Final residual = 0.000440192, No Iterations 20
time step continuity errors : sum local = 6.20467e-07, global = -1.58818e-09, cumulative = -1.65313e-09
GAMG:  Solving for p, Initial residual = 0.00116681, Final residual = 1.03143e-05, No Iterations 20
time step continuity errors : sum local = 4.39509e-07, global = -2.29711e-10, cumulative = -1.88284e-09
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 0.104692, Final residual = 7.13485e-07, No Iterations 10
DILUPBiCG:  Solving for Uy, Initial residual = 0.293159, Final residual = 1.94508e-06, No Iterations 9
DILUPBiCG:  Solving for Uz, Initial residual = 0.0575179, Final residual = 1.97897e-07, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.000959239, Final residual = 0.000823077, No Iterations 20
time step continuity errors : sum local = 3.69348e-05, global = -2.23494e-08, cumulative = -2.42322e-08
lpz_michele is offline   Reply With Quote

Old   October 15, 2015, 10:43
Default
  #20
Senior Member
 
RodriguezFatz's Avatar
 
Philipp
Join Date: Jun 2011
Location: Germany
Posts: 1,297
Rep Power: 26
RodriguezFatz will become famous soon enough
Firstly, this is not the reason for your problems but: In my view the relTol of 1e-5 is insane. You see it in the later iterations that the linear solvers take 10 or even 65 iterations to converge. Also pressure takes all 20 iterations. In some ordinary simulations all solvers should take something like less than 5 iterations, often just 1 or two. And the outer algorithm still remains stable. There is no increase in overall accuracy by setting relTol to such a low value.

Now, to the actual problem. PIMPLE doesn't converge. You can see the residuals increasing after some initial decrease. You can try to use more stable numerics in the fvSchemes file. Use (1st order) upwind also for the velocity and uncorrected scheme for laplacian instead of corrected. Did you try to initialize the flow field with potentialFoam?
__________________
The skeleton ran out of shampoo in the shower.
RodriguezFatz is offline   Reply With Quote

Reply

Tags
debug information, floating point exception

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
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


All times are GMT -4. The time now is 22:43.