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

Maximum number of iterations exceeded - rhoSimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 5, 2016, 16:01
Default Maximum number of iterations exceeded - rhoSimpleFoam
  #1
New Member
 
A Z
Join Date: Dec 2015
Posts: 11
Rep Power: 3
don6alfonso is on a distinguished road
Hi everyone,
I'm try to simulate a Pipe using rhoSimpleFoam. I get after a few teps the "Maximum number of iterations exceeded" error.

I tried some solutions I found here but without results. (changing the relaxationFactors & schemes, setting TMin, Tmax ... , checking and changing Boundary Conditions)

The Mesh looks fine.

Below the logfile:

Code:
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p	 tolerance 0.001
    field U	 tolerance 0.001
    field h	 tolerance 0.001
    field "(k|epsilon|omega)"	 tolerance 0.001

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
}

No MRF models present

No finite volume options present


Starting time loop

Time = 1e-05

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.022547925, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.055572477, No Iterations 3
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.05735183, No Iterations 3
smoothSolver:  Solving for h, Initial residual = 0.99560741, Final residual = 0.00054553968, No Iterations 1
GAMG:  Solving for p, Initial residual = 1, Final residual = 6.8005097e-09, No Iterations 35
GAMG:  Solving for p, Initial residual = 0.029825433, Final residual = 7.8671775e-09, No Iterations 17
GAMG:  Solving for p, Initial residual = 0.0016980265, Final residual = 8.9338721e-09, No Iterations 13
time step continuity errors : sum local = 7.1954987e-10, global = -9.7011997e-12, cumulative = -9.7011997e-12
rho max/min : 1.2639489 1.2014065
smoothSolver:  Solving for epsilon, Initial residual = 0.18539279, Final residual = 0.018263051, No Iterations 2
bounding epsilon, min: -1.0623939 max: 113.81697 average: 1.8930911
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 0.097817915, No Iterations 2
ExecutionTime = 251.38 s  ClockTime = 251 s

Time = 2e-05

smoothSolver:  Solving for Ux, Initial residual = 0.028023654, Final residual = 0.00072812086, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.36155382, Final residual = 0.033268639, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.33901192, Final residual = 0.0069851927, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 2.885836e-06, Final residual = 6.91766e-09, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.021754383, Final residual = 5.2425385e-09, No Iterations 31
GAMG:  Solving for p, Initial residual = 0.88325208, Final residual = 8.0467665e-09, No Iterations 26
GAMG:  Solving for p, Initial residual = 0.20759689, Final residual = 8.2151254e-09, No Iterations 22
time step continuity errors : sum local = 2.9362986e-08, global = 2.2426915e-10, cumulative = 2.1456795e-10
rho max/min : 1.4111591 1.1097508
smoothSolver:  Solving for epsilon, Initial residual = 0.99460634, Final residual = 0.087019508, No Iterations 3
bounding epsilon, min: -4188.6099 max: 58007085 average: 1010.3673
smoothSolver:  Solving for k, Initial residual = 0.99980048, Final residual = 0.072515354, No Iterations 3
bounding k, min: -56.453502 max: 8.5420722e+08 average: 820.94691
ExecutionTime = 503.97 s  ClockTime = 504 s

Time = 3e-05

smoothSolver:  Solving for Ux, Initial residual = 0.72842759, Final residual = 0.01241609, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.63153151, Final residual = 0.028186555, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.72692847, Final residual = 0.056586838, No Iterations 2
smoothSolver:  Solving for h, Initial residual = 1, Final residual = 0.0014389051, No Iterations 1
GAMG:  Solving for p, Initial residual = 0.98094699, Final residual = 8.5390336e-09, No Iterations 54
GAMG:  Solving for p, Initial residual = 6.8142427e-06, Final residual = 7.5921063e-09, No Iterations 14
GAMG:  Solving for p, Initial residual = 5.3619733e-06, Final residual = 7.6303402e-09, No Iterations 5
time step continuity errors : sum local = 0.0049250584, global = -9.5693406e-05, cumulative = -9.5693192e-05
rho max/min : 1.5289273 1.0278007
smoothSolver:  Solving for epsilon, Initial residual = 0.83663086, Final residual = 0.080689562, No Iterations 2
bounding epsilon, min: -2182.0455 max: 1.0398029e+16 average: 9.9707159e+09
smoothSolver:  Solving for k, Initial residual = 0.99736547, Final residual = 0.069018573, No Iterations 3
bounding k, min: -199.49558 max: 9.9985416e+13 average: 1.6617361e+08
ExecutionTime = 745.56 s  ClockTime = 747 s

Time = 4e-05

smoothSolver:  Solving for Ux, Initial residual = 0.43323326, Final residual = 0.011644508, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.25795325, Final residual = 0.010723061, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.29775014, Final residual = 0.0085255883, No Iterations 1
smoothSolver:  Solving for h, Initial residual = 0.089014636, Final residual = 0.00040056378, No Iterations 1


--> FOAM FATAL ERROR: 
Maximum number of iterations exceeded

    From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar)const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>]
    in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::THs(double, double, double) const at ??:?
#3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?
#5  ? at ??:?
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  ? at ??:?
Abort (core dumped)
CheckMesh

Code:
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           3151471
    faces:            8646263
    internal faces:   8378086
    cells:            2760772
    faces per cell:   6.1665176
    boundary patches: 4
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     2297184
    prisms:        73432
    wedges:        0
    pyramids:      0
    tet wedges:    800
    tetrahedra:    2
    polyhedra:     389354
    Breakdown of polyhedra by number of faces:
        faces   number of cells
            4   24349
            5   16295
            6   21202
            7   206911
            8   47896
            9   50514
           10   1042
           11   352
           12   15902
           13   32
           14   70
           15   4696
           16   1
           17   1
           18   91

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
                   Patch    Faces   Points                  Surface topology
                   inlet     4097     4165  ok (non-closed singly connected)
        outlet_beifahrer     6283     7226  ok (non-closed singly connected)
           outlet_fahrer     6262     7211  ok (non-closed singly connected)
                    wall   251535   276060  ok (non-closed singly connected)

Checking geometry...
    Overall domain bounding box (1.315 -0.38629527 -0.10150999) (3.3232987 0.38629978 0.042797109)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (1.0739519e-16 1.9630194e-16 1.0030706e-16) OK.
    Max cell openness = 1.6036898e-15 OK.
    Max aspect ratio = 43.785327 OK.
    Minimum face area = 3.457976e-10. Maximum face area = 1.8677262e-05.  Face area magnitudes OK.
    Min volume = 8.0322233e-13. Max volume = 8.0108398e-08.  Total volume = 0.012347979.  Cell volumes OK.
    Mesh non-orthogonality Max: 64.69395 average: 8.0658553
    Non-orthogonality check OK.
    Face pyramids OK.
 ***Max skewness = 10.13589, 75 highly skew faces detected which may impair the quality of the results
  <<Writing 75 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Failed 1 mesh checks.

End
Boundary

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}

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

internalField   uniform ( 0 0 0 );

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform ( 20 0 0 );
    }
    outlet_fahrer
    {
        type            inletOutlet;
        value           uniform ( 0 0 0 );
        inletValue      uniform ( 0 0 0 );
    }
    outlet_beifahrer
    {
        type            inletOutlet;
        value           uniform ( 0 0 0 );
        inletValue      uniform ( 0 0 0 );
    }
    wall
    {
        type            fixedValue;
        value           uniform ( 0 0 0 );
    }
}
Code:
FoamFile
{
    object      T;
}

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 293.15;

boundaryField
{
    wall
    {
        type            zeroGradient;
    }

   inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet_fahrer
    {
        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;    
    }
    
    outlet_beifahrer
    {
        type            inletOutlet;
        value           $internalField;
        inletValue      $internalField;
    }
}
Code:
FoamFile
{
    object      p;
}

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

internalField   uniform 101325;

boundaryField
{
    wall
    {
        type            zeroGradient;
    }

    outlet_fahrer
    {
        type            fixedValue;
        value           uniform 101325;
    }

    outlet_beifahrer
    {
        type            fixedValue;
        value           uniform 101325;
    }
    
    inlet
    {
        type            zeroGradient;
    }
}
Code:
{
    object      k;
}

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

internalField   uniform 0.01500; //82.876

boundaryField
{
   wall
    {
        type            kqRWallFunction;
        value           $internalField;
    }
    inlet
    {
        type            fixedValue; //turbulentIntensityKineticEnergyInlet;
        //intensity       0.03804;
        value           $internalField;
    }
    outlet_fahrer
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
    outlet_beifahrer
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
}
Code:
FoamFile
{
    object      epsilon;
}

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

internalField   uniform 0.008366;

boundaryField
{
        wall
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           $internalField;
    }
    inlet
    {
	type		fixedValue;
        //type            turbulentMixingLengthFrequencyInlet;
	// mixingLength    0.0019;
        value           $internalField;
    }
    outlet_fahrer
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
    outlet_beifahrer
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
}
fvSchemes

Code:
FoamFile
{
    object      fvSchemes;
}

ddtSchemes
{
    default steadyState;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(U)		Gauss linear; //new

}

divSchemes
{
   default            		none;
   div(phi,U)                    bounded Gauss linearUpwind grad(U); //bounded Gauss upwind;
   div(phi,h)                    bounded Gauss linearUpwind grad(U); //bounded Gauss upwind;
   div(phi,k)                    bounded Gauss linearUpwind grad(U); //bounded Gauss upwind;
   div(phi,K)                    bounded Gauss linearUpwind grad(U); //bounded Gauss upwind;
   div(phi,epsilon)              bounded Gauss linearUpwind grad(U); //bounded Gauss upwind;
   div(phi,omega)                bounded Gauss linearUpwind grad(U); //bounded Gauss upwind;
   div(U,p)               	 bounded Gauss linearUpwind grad(U); //bounded Gauss upwind;
   div(phi)			 bounded Gauss linearUpwind grad(U); //bounded Gauss upwind;
   //div((muEff*dev2(T(grad(U)))))  Gauss linear;
   div(((rho*nuEff)*dev2(T(grad(U)))))  Gauss linear;  //Gauss linear corrected;
}

laplacianSchemes
{
    default                          none;

    laplacian(muEff,U)               Gauss linear corrected;
    laplacian((rho|A(U)),p)          Gauss linear corrected;
    laplacian(DkEff,k)               Gauss linear corrected;
    laplacian(DepsilonEff,epsilon)   Gauss linear corrected;
    laplacian(DomegaEff,omega)       Gauss linear corrected;
    laplacian(alphaEff,h)            Gauss linear corrected;
    laplacian(mut,U)                 Gauss linear corrected;
    laplacian((rho*(1|A(U))),p)      Gauss linear corrected;
    laplacian(yPsi)                  Gauss linear corrected;
    laplacian((rho*nuEff),U)         Gauss linear corrected;
    laplacian((rho*DomegaEff),omega) Gauss linear corrected;    
    laplacian((rho*DepsilonEff),epsilon) Gauss linear corrected;
    laplacian((rho*DkEff),k)         Gauss linear corrected;   

}

interpolationSchemes
{
    default              linear;
    interpolate(U)       linear;
}

wallDist
{
    method Poisson;
// Optional entry enabling the calculation
// of the normal-to-wall field
    nRequired false;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p;
}
fvSolution

Code:
FoamFile
{
    object          fvSolution;
}

solvers
{
   p
    {
        solver           GAMG;
	preconditioner   GAMG;
        tolerance        1e-8; //war 1e-7
        relTol           0; //war 0.1
        smoother         DICGaussSeidel; //war GaussSeidel
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration on;
        agglomerator     faceAreaPair;
        nCellsInCoarsestLevel 10; //war 10
        mergeLevels      1;
	directSolveCoarsest	false; //war nicht vorhanden
    };

    U
    {
/*
//     	solver		 smoothSolver;
//         smoother         GaussSeidel;
//         nSweeps          2;
        solver		     PBiCG;
        preconditioner   DILU;
        tolerance        1e-7;
        relTol           0.01;
*/
        solver           smoothSolver;
        smoother         GaussSeidel;
        tolerance        1e-8;
        relTol           0.1;
        nSweeps          1;
    };

    h
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-06;
        relTol          0.1;
	
	/*
        solver		 PBiCG;
        preconditioner   DILU;
        tolerance        1e-5;
        relTol           0.1;
	*/
	
    };

    k
    {
       	solver		 smoothSolver;
        smoother         GaussSeidel;
        nSweeps          1; //war 2
        tolerance        1e-8; //war 5
        relTol           0.1;
    };



    epsilon
    {
       	solver		 smoothSolver;
        smoother         GaussSeidel;
        nSweeps          1; //war 2        
	tolerance        1e-8; //war 5
        relTol           0.1;
    };


    omega
    {
       	solver		 smoothSolver;
        smoother         GaussSeidel;
        nSweeps          1; //war 2        
	tolerance        1e-8; //war 5
        relTol           0.1;
    };
    
    yPsi
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-5;
        relTol          0;
    }
    
}

SIMPLE
{
    nUCorrectors                 2;
    nNonOrthogonalCorrectors 	 2; //war 5
    //pMin                         pMin   [1 -1 -2 0 0 0 0] 0;
    TMin			 TMin	[0 0 0 1 0 0 0] 283;
    Tmax			 TMax	[0 0 0 1 0 0 0] 313;
    rhoMin 			 rhoMin [1 -3 0 0 0 0 0] 0.7;
    rhoMax                       rhoMax [1 -3 0 0 0 0 0] 2;

    residualControl
    {
        p               1e-3;
        U               1e-3;
	h               1e-3;
        "(k|epsilon|omega)" 1e-3;
    }
}

relaxationFactors
{
    p               0.2; //war 0.2
    rho             0.2; //war 0.2
    U               0.6; //war 0.6
    h               0.1; //war 0.6
    k               0.6; //war 0.6
    epsilon         0.6; //war 0.6
    omega	    0.6; //war 0.6
}
Thanks in advance
don6alfonso is offline   Reply With Quote

Old   July 6, 2016, 05:24
Default
  #2
Member
 
carno
Join Date: Mar 2009
Posts: 41
Rep Power: 10
Carno is on a distinguished road
Are you using simpleFoam solution to initialize this? Or using simpleFoam files?
Carno is offline   Reply With Quote

Old   Yesterday, 06:52
Default
  #3
Senior Member
 
Attesz's Avatar
 
Attesz
Join Date: Mar 2009
Location: Munich
Posts: 361
Rep Power: 10
Attesz is an unknown quantity at this point
Hi don6alfonso, did you find a solution to your problem at the end?
__________________
https://cfd-freelancing.com
Attesz is online now   Reply With Quote

Old   Yesterday, 08:21
Default
  #4
New Member
 
A Z
Join Date: Dec 2015
Posts: 11
Rep Power: 3
don6alfonso is on a distinguished road
Quote:
Originally Posted by Attesz View Post
Hi don6alfonso, did you find a solution to your problem at the end?
Yes, I initialized the internalField in the 0/U and 0/p with the solution from potentialFoam. And I change the relaxationFactors for rho to 0.05.



don6alfonso is offline   Reply With Quote

Reply

Tags
rhosimplefoam

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Foam::error::PrintStack almir OpenFOAM Running, Solving & CFD 69 January 20, 2017 08:28
chtMultiRegionSimpleFoam: maximum number of iterations excedeed. Nkl OpenFOAM Running, Solving & CFD 18 February 16, 2016 12:42
Floating point exception error lpz_michele OpenFOAM Running, Solving & CFD 53 October 19, 2015 02:50
Unstabil Simulation with chtMultiRegionFoam mbay101 OpenFOAM Running, Solving & CFD 13 December 28, 2013 14:12
calculation stops after few time steps sivakumar OpenFOAM Running, Solving & CFD 7 March 17, 2013 07:37


All times are GMT -4. The time now is 03:07.