[don6alfonso]

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:

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

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

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 );
    }
}

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;
    }
}

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;
    }
}

{
    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;
    }
}

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

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

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

[Carno]

Are you using simpleFoam solution to initialize this? Or using simpleFoam files?

[Attesz]

Hi don6alfonso, did you find a solution to your problem at the end?

[don6alfonso]

Originally Posted by Attesz

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.