---

Konvergenz.png

Hallo Foamers,

Ich habe folgendes Problem: bei der Simulation von zwei Region (KK für Kühlkörper und bottomAir) nach Konvektion bekomme ich ständig die folgende Fehlermeldung:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.2.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.2.0
Exec  : chtMultiRegionSimpleFoam
Date  : Aug 09 2013
Time  : 13:30:00
Host  : "gpi12175.iavgroup.local"
PID    : 13403
Case  : /usr2/bottomAir_KK/case
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create fluid mesh for region bottomAir for time = 0

Create solid mesh for region KK for time = 0

*** Reading fluid mesh thermophysical properties for region bottomAir

    Adding to thermoFluid

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

    Adding to rhoFluid

    Adding to UFluid

    Adding to phiFluid

    Adding to gFluid

    Adding to turbulence

Selecting turbulence model type laminar
    Adding to ghFluid

    Adding to ghfFluid

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region KK

    Adding to thermos

Selecting thermodynamics package
{
    type            heSolidThermo;
    mixture        pureMixture;
    transport      constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to radiations

Selecting radiationModel opaqueSolid
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel none
    Adding fvOptions

No finite volume options present

Time = 1

Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.024221, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.0239137, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.0235421, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 0.0182296, No Iterations 3
Min/max T:294.282 300.428
GAMG:  Solving for p_rgh, Initial residual = 0.841536, Final residual = 0.00288963, No Iterations 6
time step continuity errors : sum local = 0.0554715, global = 2.33029e-18, cumulative = 2.33029e-18
Min/max rho:1.15697 1.18114

Solving for solid region KK
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0205543, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.988 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 1.8 s  ClockTime = 2 s

Time = 2

Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.492912, Final residual = 0.0174365, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.197266, Final residual = 0.01101, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.488943, Final residual = 0.0202762, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.426496, Final residual = 0.0100866, No Iterations 3
Min/max T:294.424 307.628
GAMG:  Solving for p_rgh, Initial residual = 0.970015, Final residual = 0.00411204, No Iterations 7
time step continuity errors : sum local = 0.0211098, global = -2.13654e-18, cumulative = 1.93753e-19
Min/max rho:1.1284 1.18805

Solving for solid region KK
DICPCG:  Solving for h, Initial residual = 0.209534, Final residual = 0.00386809, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.966 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 2.47 s  ClockTime = 3 s

Time = 3

Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.471824, Final residual = 0.0371552, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.521667, Final residual = 0.0517465, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.470758, Final residual = 0.0320942, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.479107, Final residual = 0.0243738, No Iterations 2
Min/max T:248.206 316.052
GAMG:  Solving for p_rgh, Initial residual = 0.717735, Final residual = 0.00578444, No Iterations 6
time step continuity errors : sum local = 0.0474212, global = 1.31978e-18, cumulative = 1.51353e-18
Min/max rho:0.859353 1.66003

Solving for solid region KK
DICPCG:  Solving for h, Initial residual = 0.109024, Final residual = 0.00190127, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.863 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 3.05 s  ClockTime = 3 s

Time = 4

Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.423113, Final residual = 0.0375891, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.725986, Final residual = 0.0713627, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.420529, Final residual = 0.0279138, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.710462, Final residual = 0.0287826, No Iterations 2
Min/max T:-687.938 321.57
GAMG:  Solving for p_rgh, Initial residual = 0.791067, Final residual = 0.00474987, No Iterations 5
time step continuity errors : sum local = 0.113992, global = -7.09645e-19, cumulative = 8.03888e-19
Min/max rho:0.2 2

Solving for solid region KK
DICPCG:  Solving for h, Initial residual = 0.0719346, Final residual = 0.00125857, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.145 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 3.64 s  ClockTime = 4 s

Time = 5

Solving for fluid region bottomAir
DILUPBiCG:  Solving for Ux, Initial residual = 0.348038, Final residual = 0.0172571, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.480059, Final residual = 0.0462242, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.363602, Final residual = 0.0222914, No Iterations 2
DILUPBiCG:  Solving for h, Initial residual = 0.446784, Final residual = 0.0196146, No Iterations 2

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

    From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
    in file /usr2/sw/OpenFOAM//OpenFOAM-2.2.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error: rintStack(Foam::Ostream&) in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double) const) const in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#3  Foam::heRhoThermo<Foam::rhoThermo, Foam: ureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4  Foam::heRhoThermo<Foam::rhoThermo, Foam: ureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam: erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5 
in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#6  __libc_start_main in "/lib64/libc.so.6"
#7 
in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
____________________________________________________________________________________________________________________
Dabei ist es klar zu sehen, dass der Konvergenz von h in der Luft instabil ist(Sehe Bild). Der Grund dafür habe ich nach mehrere Versuche indem ich viel BC benutzt, Relaxationfaktor geändert, div Scheme geändert, andere Netz benutzt... immer noch nicht gefunden.
____________________________________________________________________________________________________________________
Hier sind die 0-Ordner und system-Ordner Dateien:

FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField  uniform 8.24e-7;

boundaryField
{
    bottomMinX
    {
      type            compressible::epsilonWallFunction;
      value          uniform 8.24e-7;
    }
    bottomMinY
      {
type          compressible::epsilonWallFunction;
value        uniform 8.24e-7;
      }
    bottomMinZ
      {
type          compressible::epsilonWallFunction;
value        uniform 8.24e-7;
      }
    bottomMaxX
      {
type          compressible::epsilonWallFunction;
value        uniform 8.24e-7;
      }
    bottomMaxY
      {
type          compressible::epsilonWallFunction;
value        uniform 8.24e-7;
      }
    bottomMaxZ
      {
type          compressible::epsilonWallFunction;
value        uniform 8.24e-7;
      }
    bottomAir_to_KK
      {
type          compressible::epsilonWallFunction;
value        uniform 8.24e-7;
      }
}

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

FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField  uniform 3.75e-5;

boundaryField
{
    bottomMinX
    {
      type            compressible::kqRWallFunction;
      value          uniform 3.75e-5;
    }
    bottomMinY
      {
type          compressible::kqRWallFunction;
value        uniform 3.75e-5;
      }
    bottomMinZ
      {
type          compressible::kqRWallFunction;
value        uniform 3.75e-5;
      }
    bottomMaxY
      {
type          compressible::kqRWallFunction;
value        uniform 3.75e-5;
      }
    bottomMaxX
      {
type          compressible::kqRWallFunction;
value        uniform 3.75e-5;
      }
    bottomMaxZ
      {
type          compressible::kqRWallFunction;
value        uniform 3.75e-5;
      }
    bottomAir_to_KK
      {
type          compressible::kqRWallFunction;
value        uniform 3.75e-5;
      }
}

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

FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField  uniform 1e5;

boundaryField
{
    bottomMinX
      {
        type            calculated;
value          uniform 1e5;
      }
    bottomMinY
      {
        type            calculated;
value          uniform 1e5;
      }
    bottomMinZ
      {
        type            calculated;
value          uniform 1e5;
      }
    bottomMaxY
      {
        type            calculated;
value          uniform 1e5;
      }
    bottomMaxX
      {
        type            calculated;
value          uniform 1e5;
      }
    bottomMaxZ
      {
type            calculated;
value          uniform 1e5;
      }
    bottomAir_to_KK
      {
type            calculated;
value          uniform 1e5;
      }
}

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

FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField  uniform 1e5;

boundaryField
{
    bottomMinX
      {
        type            fixedFluxPressure;
value          uniform 1e5;
      }
    bottomMinY
      {
        type            fixedFluxPressure;
value          uniform 1e5;
      }
    bottomMinZ
      {
        type            fixedFluxPressure;
value          uniform 1e5;
      }
    bottomMaxX
      {
        type            fixedFluxPressure;
value          uniform 1e5;
      }
    bottomMaxY
      {
        type            fixedFluxPressure;
value          uniform 1e5;
      }
    bottomMaxZ
      {
type            fixedFluxPressure;
value          uniform 1e5;
      }
    bottomAir_to_KK
      {
type            fixedFluxPressure;
value          uniform 1e5;
      }
}

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

FoamFile
{
    version    2.0;
    format      binary;
    class      volScalarField;
    location    "1";
    object      T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField  uniform 297;

boundaryField
{
    bottomMinX
    {
      type                zeroGradient;
    }
    bottomMinY
    {
      type                zeroGradient;
    }
    bottomMinZ
    {
      type                zeroGradient;
    }
    bottomMaxX
    {
      type                zeroGradient;
    }
    bottomMaxY
    {
      type                zeroGradient;
    }
    bottomMaxZ
    {
      type                zeroGradient;
    }
    bottomAir_to_KK
    {
      type                compressible::turbulentTemperatureCoupledBaffleMixed;
      value              uniform 297;
      neighbourFieldName  T;
      kappa              fluidThermo;
      kappaName          none;
    }
}

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

FoamFile
{
    version    2.0;
    format      ascii;
    class      volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField  uniform ( 0 -0.1 0 );

boundaryField
{
    bottomMinX
      {
        type            fixedValue;
value          uniform (0 0 0);
      }
    bottomMinY
      {
        type            fixedValue;
value          uniform (0 0 0);
      }
    bottomMinZ
      {
        type            fixedValue;
value          uniform (0 0 0);
      }
    bottomMaxX
      {
        type            fixedValue;
value          uniform (0 0 0);
      }
    bottomMaxY
      {
        type            fixedValue;
value          uniform (0 0 0);
      }
    bottomMaxZ
      {
type            fixedValue;
        value          uniform (0 0 0);
      }
    bottomAir_to_KK
      {
type            fixedValue;
value          uniform (0 0 0);
      }
}

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

FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default steadyState;
}

gradSchemes
{
    default        Gauss linear;
}

divSchemes
{
    default        none;
    div(phi,U)      bounded Gauss upwind;
    div(phi,K)      bounded Gauss upwind;
    div(phi,h)      bounded Gauss upwind;
    div(phi,k)      bounded Gauss upwind;
    div(phi,K)      bounded Gauss upwind;
    div(phi,epsilon) bounded Gauss upwind;
    div(phi,R)      bounded Gauss upwind;
    div(R)          Gauss linear;
    div((muEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default        none;
    laplacian(muEff,U) Gauss linear uncorrected;
    laplacian(Dp,p_rgh) Gauss linear uncorrected;
    laplacian(alphaEff,h) Gauss linear uncorrected;
    laplacian(DkEff,k) Gauss linear uncorrected;
    laplacian(DepsilonEff,epsilon) Gauss linear uncorrected;
    laplacian(DREff,R) Gauss linear uncorrected;
}

interpolationSchemes
{
    default        linear;
}

snGradSchemes
{
    default        uncorrected;
}

fluxRequired
{
    default        no;
    p_rgh;
}

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

FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    rho
    {
        solver          PCG
        preconditioner  DIC;
        tolerance      1e-7;
        relTol          0;
    }

    p_rgh
    {
        solver          GAMG;
        tolerance        1e-7;
        relTol          0.01;

        smoother        GaussSeidel;

        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator    faceAreaPair;
        mergeLevels      1;
    }

    "(U|h|k|epsilon)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance        1e-7;
        relTol          0.1;
    }
}

SIMPLE
{
    momentumPredictor on;
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue      100000;
    rhoMin          rhoMin [1 -3 0 0 0] 0.2;
    rhoMax          rhoMax [1 -3 0 0 0] 2;
}

relaxationFactors
{
    fields
    {
        rho            1;
        p_rgh          0.3;
    }
    equations
    {
        U              0.7;
        h              0.7;
        nuTilda        0.7;
        k              0.3;
        epsilon        0.3;
        omega          0.7;
        "ILambda.*"    0.7;
    }
}

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

Falls dass jemand eine Idee hat woran könnte es liegen, dass meine Berechnung immerwieder abstürtzt. Bitte Bitte mitteilen.

Ich bedanke mich im voraus.

MfG

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Problem gelöst!! Danke für durchlesen 

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo mbay101,
ich habe ein ähnliches Problem wie du, ich bekomme die Fehlermeldung

-> FOAM FATAL ERROR:
Maximum number of iterations exceeded

    From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
    in file /usr2/sw/OpenFOAM//OpenFOAM-2.2.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

Wie hast du denn das Problem gelöst?

MfG

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo zusammen,

bitte verwendte die CODE TAGS damit man das eingefügte von euch besser lesen kann.

Dein Fehler tritt auf, wenn du nen schlechtes Netz hast, Relaxation zu hoch, falsche RB oder falsche Initialisierungswerte hast.

Weiteres fällt mir gerade nicht ein.

------------------
Grüße Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo SummerOf69,

eigentlich solle ich hallo Winter sagen aber ist ok  . Maximum number of iterations exceeded tritt auf, wenn die gerechnete variablen enorme groß werden. Was deine Simulationsergebnisse kein sinn mehr machen. (Bei mir ist die Enthalpie bzw. die Temperaturen immer nach 5 Iterationen immer ein sehr höher Zahl erreicht haben). kannst du mir die folgende Fragen beantworten:

1. kommt die Fehlermeldung nach einpaar Iterationen? oder startet OpenFoam gar nicht mit der Berechnung?
2. Wie sieht checkMesh aus? welche Netz verwendest du?
3. hast du schon mit dem Tutorial fvScheme und fvSolution die Simulation gestartet?
4. wie sieht dein 0 Ordner aus? (Boundary)

meistens sind Boundary, die das Problem verursachen, wenn checkMesh, fvScheme und fv Solution stimmen.

Gruß

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo mbay101,
danke für die Hinweise!      

1. Die Fehlermeldung kommt direkt zu Beginn, gleich nach dem ersten Zeitschritt:

Code:

---

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.1-57f3c3617a2d
Exec   : sonicFoam
Date   : Sep 02 2013
Time   : 14:20:23
Host   : "basti-MS-7693"
PID    : 4313
Case   : /OpenFOAM/user/run/2D
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0

Reading thermophysical properties

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

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

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

Creating field kinetic energy K

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Starting time loop

Time = 5e-07

Courant Number mean: 0 max: 0
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 3.16582e-06, No Iterations 8
DILUPBiCG:  Solving for Uy, Initial residual = 0.999917, Final residual = 9.72863e-06, No Iterations 8
DILUPBiCG:  Solving for e, Initial residual = 0.00113303, Final residual = 9.6655e-06, No Iterations 1
DILUPBiCG:  Solving for p, Initial residual = 1, Final residual = 5.89459e-09, No Iterations 17
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.90331e-10, global = 1.91658e-11, cumulative = 1.91658e-11
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 0.756761, Final residual = 3.48428e-06, No Iterations 10
DILUPBiCG:  Solving for Uy, Initial residual = 0.919145, Final residual = 9.66893e-06, No Iterations 9
DILUPBiCG:  Solving for e, Initial residual = 1, Final residual = 4.11551e-06, No Iterations 10
DILUPBiCG:  Solving for p, Initial residual = 0.254589, Final residual = 5.77824e-09, No Iterations 415
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.34465e-08, global = 7.04829e-09, cumulative = 7.06746e-09
DILUPBiCG:  Solving for epsilon, Initial residual = 0.973163, Final residual = 2.79042e-09, No Iterations 14
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 2.76766e-09, No Iterations 22
bounding k, min: -3.96286e+08 max: 8.97994e+07 average: -164209
ExecutionTime = 0.19 s  ClockTime = 0 s

Time = 1e-06

Courant Number mean: 0.00253351 max: 4052.93
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.871822, Final residual = 6.72156e-06, No Iterations 397
DILUPBiCG:  Solving for Uy, Initial residual = 0.90218, Final residual = 3.37918e-06, No Iterations 439
DILUPBiCG:  Solving for e, Initial residual = 0.999989, Final residual = 4.52059e-06, No Iterations 397

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

    From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error:      rintStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::species::thermo<Foam::hConstThermo<Foam:      erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam:      erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam:      erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam:      erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::*)(double) const) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#3  Foam::hePsiThermo<Foam:      siThermo, Foam:      ureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam:      erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4  Foam::hePsiThermo<Foam:      siThermo, Foam:      ureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam:      erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/sonicFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/sonicFoam"
Aborted (core dumped)

---

2. Bei dem Netz handelt es sich um ein simples 2D-Netz, erstellt mit BlockMesh, checkMesh findet keine Fehler. Hier ein Abbild:
2D_Mesh

3. Was meinst du mit "dem Tutorial fvScheme und fvSolution die Simulation gestartet" ?
Ich habe das Tutorial compressible/sonicFoam/ras/prism als Vorlage verwendet, das BlockMesh und die Randbedingungen im 0-Ordner angepasst.

4. Hier die Boundaries (0-Ordner):

Code:

---

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

internalField   uniform 266000;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 266000;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 266000;
        value           uniform 266000;
    }
    wall
    {
        type            compressible::epsilonWallFunction;
        value           uniform 266000;
    }
    defaultFaces
    {
        type            empty;
    }
}

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

internalField   uniform 1000;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 1000;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 1000;
        value           uniform 1000;
    }
    wall
    {
        type            compressible::kqRWallFunction;
        value           uniform 1000;
    }
    defaultFaces
    {
        type            empty;
    }
}

mut
boundaryField
{
    inlet
    {
        type            calculated;
        value           uniform 0;
    }
    outlet
    {
        type            calculated;
        value           uniform 0;
    }
    wall
    {
        type            mutkWallFunction;
        value           uniform 0;
    }
    defaultFaces
    {
        type            empty;
    }
}

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

internalField   uniform 100000;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 600000;
    }

    outlet
    {
        type            waveTransmissive;
        field           p;
        phi             phi;
        rho             rho;
        psi            thermo si;
        gamma           1.3;
        fieldInf        100000;
        lInf            1;
        value           uniform 100000;
    }
    wall
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

T
dimensions      [0 0 0 1 0 0 0];

internalField   uniform 300;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 300;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform 300;
        value          uniform 300;
    }
    wall
    {
        type            inletOutlet;
        inletValue      uniform 300;
        value          uniform 300;
    }

    defaultFaces
    {
        type            empty;
    }
}

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

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value           uniform (0 0 0);
    }
    wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    defaultFaces
    {
        type            empty;
    }
}

---

Meine Absicht war, die Geschwindigkeitsentwicklung in einer Düse zu Untersuchen. Dabei habe ich die Drücke vorgegeben (im inlet: P=6bar, outlet: P=1bar)

Der Fehler (Maximum number of iterations exceeded) tritt nur auf, wenn ich im BlockMesh die Maße in [mm] ändere (convertToMeters 0.001). Belasse ich sie auf [m], läuft die Simulation fehlerfrei.

Hast du eine Idee, an welchen RB es liegen könnte?

MfG

[Diese Nachricht wurde von SummerOf69 am 02. Sep. 2013 editiert.]

[Diese Nachricht wurde von SummerOf69 am 02. Sep. 2013 editiert.]

[Diese Nachricht wurde von SummerOf69 am 02. Sep. 2013 editiert.]

[Diese Nachricht wurde von SummerOf69 am 02. Sep. 2013 editiert.]

[Diese Nachricht wurde von SummerOf69 am 02. Sep. 2013 editiert.]

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo,

also einige Fehler und Anmerkungen:

T - File -> Wall: wieso inletOutlet? Versteh ich nicht. Durch die Wand kannst du keine Rückströmung erhalten weil du im U File den Geschwindigkeitsvektor U (0 0 0)^T vorgibst.

p - File -> Outlet: Hab diese BC noch nie benutzt aber musst du diese verwenden oder hast du einen Grund für diese; schon mal fixedValue auch am Outlet versucht?

U - File -> internalField: du hast hier auch eine Initialisierung von 0 m/s im gesamten Feld. Insofern du nur den stationären Zustand simulierst, empfehle ich dir das U-Profil abzuschätzen (Bernoulli) und dann einzutragen; dies löst oft viele Probleme. Oder allgemein einfach ein ganz kleines U-Profil vorzugeben (0 0 0.005) etc. Aber nach deinem Solver gehe ich auf eine transiente Betrachtung aus.

Epsilon / k - File -> Hast du die Werte berechnet? Möglich wäre am Inlet auch mit anderen BC zu arbeiten. Solltest du diese allerdings berechnet haben, dürfte es auch funktionieren.

Allgemein:

Wie du siehst ist in deiner zweiten Iteration die CO-Zahl schon über dem akzeptablen Bereich Co < 1 sollte gewährleistet sein (in den meisten Fällen). Du bist bereits bei Co > 4000.

Strömungsart? Subsonisch, Hypersonisch etc.? Sollte Ma < 0,7 würd ich den rhoSimpleFoam Löser für dein Problem verwenden; bzw. für transiente Problemstellungen den rhoPimpleFoam.

------------------
Grüße Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo,
danke für eure Anmerkungen 

Ich möchte transsonisch/hypersonisch sowie transient simulieren, daher der sonicFoam solver.
Die Werte für alphaT, k, und epsilon entprechen denen des prism-Tutorials, ich habe auch mit selbst berechneten Werten simuliert -> gleicher Fehler.
Sogar mit ausgeschalteter Turbulenz (simulationType  laminar -> kein alphaT, k und epsilon File) kommt die selbe Fehlermeldung.

Bezüglich der Randbedingungen für T, P und U:
Ich habe mich da an das prism Tutorial gehalten, weil das meinem Rechenfall am ähnlichsten ist. Ich habe es auch schon mit anderen Randbedingungen getestet (am outlet zeroGradient, fixedValue), ohne Erfolg.

Der Fehler tritt auch auf, wenn ich ein kleines U-Profil vorgebe (mit Bernoulli abzuschätzen macht wenig Sinn, da die Luft anfangs ruht).
Mit deltaT habe ich auch schon rumgespielt, bei Co-Zahlen kleiner 1 kam der Fehler ebenfalls.

Für weitere Hinweise wäre ich euch sehr dankbar, ich beschäftige mich seit einigen Wochen intensiv mit diesem Thema, finde hierfür aber keine Lösung.
Wäre denn der rhoCentralFoam solver auch gut für hypersonische Strömungen geeignet (transiente Simulation)?

MfG

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Das Beste wäre dein Case zu haben. Dann kann man kurz drüberschaun.

------------------
Grüße Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

nozzle_2D.zip

Hey,
ich wäre euch wirklich dankbar, wenn ihr mein case mal angucken könntet 
Habe es als .zip Datei angehängt.

MfG

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

sonicFoam.png

Die Antwort auf deine Frage ist nich so einfach aber erstmal:

Der Solver hat keine Implementierung für die Adjustierung der Co-Zahl. Das müsstest du schon selber vornehmen und bemerkt haben.

Du brauchst für solche Fälle ein extrem kleinen Zeitschritt (siehe Anhang);
In so einem Problem hast du eine Druchwelle die sich ausbreitet - laut Ferziger ist das glaub ein hyperbolisches Problem in dem dann Informationen zur Strömungsrichtung und entgegen diese ausgebreitet werden. Müsst ich aber nochmals nachlesen was dann genau zu beachten ist.

Ach ja dein Netz ist am Übergang auch sehr schlecht. Wieso wirst du bei der Düse vorne so extrem weitmaschig?

U: pressureInletVelocity und pressureInletOutletVelocity hab ich mal ausprobiert und dadurch hab ich mal nen Geschwindigkeitsfeld erhalten.

Bei mir bricht er immer noch ab aber ich hab zumindest mal ein paar Analysepunkte von der Strömung bekommen.

------------------
Grüße Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo,
vielen Dank für die Bemühungen! 

Die Adjustierung der Co-Zahl hatte ich versuchsweise mal implementiert, habe vergessen sie wieder zu entfernen 
Das Netz habe ich etwas homogener und feiner gestaltet, leider auch keine Verbesserung ...
Ein extrem kleines DeltaT scheint auch keine Lösung auf das Problem zu sein (es laufen deutlich mehr Iterationen durch, die Fehlermeldung erscheint jedoch weiterhin, nur "später", längere Rechenzeit). Bei einem DeltaT von 5e-13 war die maximale Co-Zahl 0.1

FerzigerŽs "Numerische Strömungsmechanik" werde ich mir mal ausleihen und das Thema nachlesen, aber was könnte man denn großartig daran ändern?

Ein weiterer Punkt der mir noch eingefallen ist:
Wäre eine Änderung der "ThermophysicalProperties" sinnvoll?
Aktuell verwende ich den Typ "hePsiThermo", ich bin mir allerdings nicht sicher, welche ThermoModel angebracht sind, und welche Attribute (mixture, transport, specie) diese verwenden müssen.

MfG

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo,

hast du versucht mit einer höhere Viskosität zu simulieren? Statt bei Luft 1,85 . 10^(-5) versuch mal mit 1,85 . 10^(-1).

Gruß

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Wie bereits erwähnt, habe ich in diesem Bereich noch keine Simulationen durchgeführt.

Die Viskosität zu ändern ist für mich nicht sonderlich sinnvoll.
Ich würde vielmehr deine RB nochmals prüfen.
Dann die Schemen auf Upwind setzen.
Dein Netz nochmals prüfen und ggf. im Netz suchen ob in solch einer Aufgabenstellung möglicherweise bestimmte Sachen beachtet werden müssen.
Vielleicht gibst du zu viele RB vor. Ich erinnere mich da an eine Gegenüberstellung der verschiedenen Strömungsformen und welche RB gesetzt werden müssen.

Wenn ein System überbestimmt ist, kann es sein das die Lösung auch falsch ist.

------------------
Grüße Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo Tobi,

bei mir war es so, wenn ich mit höhe Viskosität simulierte hat es immer problemlos funktioniert. Die Strömungslinie sahen sogar ganz gut aus. Aber das ganze war falsche, da ich nicht mit Luft sondern mehr mit Hönig simuliert habe  .

Re= V.l/mu vobei V die gesuchte ist und mu eine Temperaturabhängige Konstante ist. Da war mir klar das ich l falsche habe und das war tätsachlich so. Ich habe die Einheit mm von CATIA zu Salome nicht übernommen. Also in OpenFoam war es am Ende statt 500mm, 500m.

Das habe ich korrigiert und...funktioniert   

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Wie sich die Reynoldszahl zusammensetzt ist mir klar.

Aber das mit der Viskosität nicht!
hast du dein Netz falsch skaliert gehabt?

------------------
Grüße Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hey Leute,
ich habe eine Lösung für das Problem gefunden.

In den "fvSchemes" habe ich folgende Einträge gemacht:

Code:

---

ddtSchemes
{
    default        Euler;
    ddt(rho,h)      bounded Euler;
}

gradSchemes
{
    default        cellMDLimited Gauss linear 1.0;
}

divSchemes
{
    default        none;
    div(phi,U)      Gauss linearUpwind grad(U);
    div(phi,e)      Gauss linearUpwind grad(U);
    div(phid,p)    Gauss linearUpwind grad(U);
    div(phi,K)      Gauss linearUpwind grad(U);
    div(phiv,p)    Gauss linearUpwind grad(U);
    div(phi,k)      Gauss linearUpwind grad(U);
    div(phi,epsilon) Gauss linearUpwind grad(U);
    div((muEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default        Gauss linear limited 1;
}

interpolationSchemes
{
    default        linear;
}

snGradSchemes
{
    default        linear;
}

fluxRequired
{
    default        no;
    p;

---

Hauptsächlich habe ich in den divSchemes die Werte "bounded Gauss upwind" in "Gauss linearUpwind grad(U)" geändert, ich spiele momentan damit herum, um eine gute Balance zwischen Stabilität und Genauigkeit zu erhalten.

Ich danke euch sehr für eure guten Tipps! 

Eine Erhöhung der Viskosität halte ich übrigens unter Berücksichtigung des Maßstabs für durchaus gerechtfertigt, um Abschätzungen zu erhalten.

MfG

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Das mit der Viskosität leuchtet mir nicht ein?

Ihr skaliert die Geometrie und verwendet dann eine andere Viskosität?
Wenn das so ist - möglich ja.

Ich dachte das die original Geometrie einfach mit einer höheren Viskosität berechnet wird.

------------------
Grüße Tobias Holzmann

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

hi

i have this problem
can anyone help me?

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

From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.0.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 67.

best regards

masoud

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP