CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   interFoam/kOmegaSST tank filling with printStackError/Mules (http://www.cfd-online.com/Forums/openfoam-solving/110489-interfoam-komegasst-tank-filling-printstackerror-mules.html)

simpomann December 12, 2012 08:39

interFoam/kOmegaSST tank filling with printStackError/Mules
 
Hey there,

I am simulating the filling of a tank with interfoam and try to use k-Omega-SST model for turbulence. The mesh has around 400.000 cells, checkMesh says okay. Unluckily I encounter these errors a lot (sometimes after 0.25s, with different transport properties after 0.4s):
Code:

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.524686, Final residual = 0.000681733, No Iterations 2
time step continuity errors : sum local = 1.32975e-10, global = 2.51813e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0920533, Final residual = 7.70789e-08, No Iterations 5
time step continuity errors : sum local = 2.09622e-14, global = -1.45744e-18, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.000158889, Final residual = 7.80272e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00244614, Final residual = 5.717e-10, No Iterations 2
bounding k, min: -0.000200721 max: 3.00901e+20 average: 1.31669e+15
ExecutionTime = 24866.4 s  ClockTime = 24916 s

Courant Number mean: 2.01358e-08 max: 0.712635
Interface Courant Number mean: 8.29938e-09 max: 0.282305
deltaT = 2.2569e-14
Time = 0.2550889169617

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.408204, Final residual = 0.000561119, No Iterations 2
time step continuity errors : sum local = 2.68439e-10, global = -5.36369e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0698315, Final residual = 5.11837e-08, No Iterations 5
time step continuity errors : sum local = 3.27413e-14, global = -7.29092e-18, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.00329008, Final residual = 2.10387e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.0124975, Final residual = 7.00019e-09, No Iterations 2
ExecutionTime = 24868.6 s  ClockTime = 24919 s

Courant Number mean: 3.59009e-08 max: 3.49031
Interface Courant Number mean: 1.43023e-08 max: 0.550178
deltaT = 5.17294e-15
Time = 0.255088916961705

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28037e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.587722, Final residual = 0.000830391, No Iterations 2
time step continuity errors : sum local = 7.64239e-11, global = 4.21113e-16, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0960292, Final residual = 3.83386e-08, No Iterations 6
time step continuity errors : sum local = 5.13861e-15, global = -9.16091e-19, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.000484156, Final residual = 1.42813e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00071846, Final residual = 1.26128e-09, No Iterations 2
ExecutionTime = 24870.9 s  ClockTime = 24921 s

Courant Number mean: 1.32675e-08 max: 0.702258
Interface Courant Number mean: 5.5049e-09 max: 0.202586
deltaT = 5.76223e-15
Time = 0.25508891696171

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.143881, Final residual = 0.000365565, No Iterations 2
time step continuity errors : sum local = 4.46931e-11, global = -4.75218e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0402696, Final residual = 3.67927e-08, No Iterations 6
time step continuity errors : sum local = 4.85566e-15, global = -1.63317e-18, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 8.34399e-05, Final residual = 9.95434e-12, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00146217, Final residual = 2.85304e-10, No Iterations 2
ExecutionTime = 24873.2 s  ClockTime = 24923 s

Courant Number mean: 1.83689e-08 max: 0.546694
Interface Courant Number mean: 7.7516e-09 max: 0.262968
deltaT = 6.60544e-15
Time = 0.255088916961717

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0613232, Final residual = 0.000199917, No Iterations 2
time step continuity errors : sum local = 3.64953e-11, global = 1.40235e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00862543, Final residual = 7.72936e-08, No Iterations 5
time step continuity errors : sum local = 1.45349e-14, global = 4.97913e-17, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.000305148, Final residual = 9.20208e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00171966, Final residual = 6.46387e-10, No Iterations 2
bounding k, min: -2.15624e+09 max: 2.68344e+21 average: 1.08957e+16
ExecutionTime = 24875.4 s  ClockTime = 24925 s

Courant Number mean: 2.361e-08 max: 0.507769
Interface Courant Number mean: 1.01369e-08 max: 0.318716
deltaT = 7.64614e-15
Time = 0.255088916961725

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00200676  Min(alpha1) = -6.28036e-20  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.24767, Final residual = 0.000361503, No Iterations 2
time step continuity errors : sum local = 1.34031e-10, global = 3.46138e-15, cumulative = -1.01128e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0585152, Final residual = 3.09369e-08, No Iterations 5
time step continuity errors : sum local = 1.38894e-14, global = 8.65978e-17, cumulative = -1.01128e-05
DILUPBiCG:  Solving for omega, Initial residual = 0.00217707, Final residual = 1.56868e-10, No Iterations 2
DILUPBiCG:  Solving for k, Initial residual = 0.00609363, Final residual = 2.66106e-09, No Iterations 2
ExecutionTime = 24877.6 s  ClockTime = 24928 s

Courant Number mean: 3.25704e-08 max: 1.92022
Interface Courant Number mean: 1.36638e-08 max: 0.455131
deltaT = 3.18553e-15
Time = 0.255088916961728

[3] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2  in "/lib64/libc.so.6"
[3] #3  void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #4  at MULES.C:0
[3] #5  Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #6 
[3]  in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/interFoam"
[3] #7  __libc_start_main in "/lib64/libc.so.6"
[3] #8 
[3]  at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
[c103-06:14726] *** Process received signal ***
[c103-06:14726] Signal: Floating point exception (8)
[c103-06:14726] Signal code:  (-6)
[c103-06:14726] Failing at address: 0x35c600003986
[c103-06:14726] [ 0] /lib64/libc.so.6(+0x34e10) [0x7fd91ae43e10]
[c103-06:14726] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7fd91ae43d95]
[c103-06:14726] [ 2] /lib64/libc.so.6(+0x34e10) [0x7fd91ae43e10]
[c103-06:14726] [ 3] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geometricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMeshEEESK_RKT0_RKT1_ddi+0x1068) [0x7fd91ca71fa8]
[c103-06:14726] [ 4] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(+0x7c007e) [0x7fd91ca6507e]
[c103-06:14726] [ 5] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERS8_dd+0x80) [0x7fd91ca65310]
[c103-06:14726] [ 6] interFoam() [0x42d35e]
[c103-06:14726] [ 7] /lib64/libc.so.6(__libc_start_main+0xed) [0x7fd91ae3023d]
[c103-06:14726] [ 8] interFoam() [0x425119]
[c103-06:14726] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 14726 on node c103-06 exited on signal 8 (Floating point exception).

Timestep decreases and this happens. I checked velocities and pressure at the last working timestep, nothing out of hand here.
I guess I have to supply my fvSchemes, fvSolutions and my initial conditions to trace this. My guess is a connection to the k-Omega-SST model. Probably my initial settings are bad, but I calculated them several times following the wiki here and can't manage to get it work.
Case is currently starting laminar, if this works I will update you asap.

U
Code:

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

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

internalField  uniform (0 0 0);

boundaryField

{
    walls
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
    //Durchströmt
   
    inlet
    {
      type            flowRateInletVelocity;
      flowRate          constant 0.000666667;//incompressible solver: [m³/s]
      value          uniform (0 0 0); //placeholder
   
    //type            fixedValue;            //Alternative, Sascha
        //value          (1.36 -2.56 0);
    }
       
    atmosphaere_CATIASTL
    {
        type            zeroGradient;
    }
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
   
    efr_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }

    entl_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
   
    pistole_vorn_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
   
   
    pistole_hinten_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
   
    innenteile_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
     
    flv_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
   
  einsatz_efr_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
   
    verschlusskopf_CATIASTL
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
   
}

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

p_rgh

Code:

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

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

internalField  uniform 0;

boundaryField

{
    walls
    {
        type            buoyantPressure;
        value          uniform 0;
    }
    //Durchströmt
   
    inlet
    {
        type            buoyantPressure;
        value          uniform 0;
    }
   
    atmosphaere_CATIASTL
  {
        type            totalPressure;
        p0              uniform 0;
        U              U;
        phi            phi;
        rho            rho;
        psi            none;
        gamma          1;
        value          uniform 0;   
    }
   
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            buoyantPressure;
        value          uniform 0;
    }
   
    efr_CATIASTL
    {
        type            buoyantPressure;
        value          uniform 0;
    }

    entl_CATIASTL
    {
        type            buoyantPressure;
        value          uniform 0;
    }
   
    pistole_vorn_CATIASTL
      {
        type            buoyantPressure;
        value          uniform 0;
    }
   
  // sd_CATIASTL
    //{
    //    type            buoyantPressure;
    //    value          uniform 0;
    //}
   
    pistole_hinten_CATIASTL
    {
        type            buoyantPressure;
        value          uniform 0;
    }
   
    innenteile_CATIASTL
    {
        type            buoyantPressure;
        value          uniform 0;
    }
   
   
    flv_CATIASTL
    {
        type            buoyantPressure;
        value          uniform 0;
    }
   
      einsatz_efr_CATIASTL
    {
        type            buoyantPressure;
        value          uniform 0;
    }
   
    verschlusskopf_CATIASTL
    {
        type            buoyantPressure;
        value          uniform 0;
    }
   
}

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

k
Code:

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

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

internalField  uniform 0.0286;

boundaryField

{
    walls
    {
        type            kqRWallFunction;
        value          $internalField;
    }
    //Durchströmt
   
    inlet
    {
        type            turbulentIntensityKineticEnergyInlet;
        intensity      0.041;      // 5% turbulent intensity
        value          $internalField;
    }
   
        atmosphaere_CATIASTL
    {
        type            inletOutlet;
        inletValue      $internalField;
    }
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
   
    efr_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }

    entl_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
   
    pistole_vorn_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
   
    sd_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
   
    pistole_hinten_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
   
    innenteile_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
   
   
    flv_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
   
      einsatz_efr_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
   
    verschlusskopf_CATIASTL
    {
        type            kqRWallFunction;
        value          $internalField;
    }
}


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

Omega
Code:

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


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

internalField  uniform 1643.71;

boundaryField

{
    walls
    {
        type            omegaWallFunction;
        value          $internalField;
    }
    //Durchströmt
   
    inlet
    {
        type            fixedValue;
        value          uniform 280.78;
    }
     
        atmosphaere_CATIASTL
    {
        type            inletOutlet;
        inletValue      $internalField;
        value          $internalField;
    }
   
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
   
    efr_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }

    entl_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
   
    pistole_vorn_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
   
    sd_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
   
    pistole_hinten_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
   
    innenteile_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
   
   
    flv_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
   
    einsatz_efr_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
   
    verschlusskopf_CATIASTL
    {
        type            omegaWallFunction;
        value          $internalField;
    }
}


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

Nut
Code:

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

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

internalField  uniform 0.00010198;

boundaryField

{
    walls
    {
        type            nutkWallFunction;
        value          $internalField;
    }
    //Durchströmt
   
    inlet
    {
        type            calculated;
        value          uniform 0.00010198;
    }
       
        atmosphaere_CATIASTL
    {
        type            calculated;
        value          $internalField;
    }
   
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   
    efr_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }

    entl_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   
    pistole_vorn_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   
    sd_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   
    pistole_hinten_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   
    innenteile_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   

   
    flv_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   
    einsatz_efr_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   
    verschlusskopf_CATIASTL
    {
        type            nutkWallFunction;
        value          $internalField;
    }
   
}

I guess inlet is self explaining. The only outlet is "atmosphaere". All the other stuff are pipes etc, so basically walls.
If somebody spots a mistake in one of the files please tell me! I am super lost with this.

All help is highly appreciated! Big thanks and best regards!
I hope somebody with more experience can help me to fix this, unluckily it's a bit urgent for my thesis.

Greetings!
Simon

simpomann December 12, 2012 08:40

fvSchemes
Code:

ddtSchemes
{
    default        Euler;
}

gradSchemes
{
    default        Gauss linear;            //Alternative:default cellMDLimited Gauss linear 1
}

divSchemes
{
    div(rho*phi,U)  Gauss limitedLinearV 1;
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
    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;
}

snGradSchemes
{
    default        corrected;
}

fluxRequired
{
    default        no;
    p;
    pcorr;
    alpha1;
    p_rgh;
}

fvSolutions
Code:

solvers
{
    p_rgh                //Lösung der Druckgleichung
    {
        solver          GAMG;
        tolerance      1e-07;
        relTol          0.01;
        smoother        DIC;
        nPreSweeps      0;
        nPostSweeps    2;
        nFinestSweeps  2;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator    faceAreaPair;
        mergeLevels    1;
    }

    p_rghFinal                //finaler Loop
    {
        solver          PCG;
        preconditioner
        {
            preconditioner  GAMG;
            tolerance      1e-07;
            relTol          0;
            nVcycles        2;
            smoother        DICGaussSeidel;
            nPreSweeps      2;
            nPostSweeps    2;
            nFinestSweeps  2;
            cacheAgglomeration true;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels    1;
        }
        tolerance      1e-07;
        relTol          0;
        maxIter        10;
    }

    pcorr                //vor dem Loop gelöste Druckvergleichung für Startwerte des Feldes innerhalb des Iterationsschritts
    {
        solver          PCG;
        preconditioner
        {
            preconditioner  GAMG;
            tolerance      1e-05;
            relTol          0;
            smoother        DICGaussSeidel;
            nPreSweeps      0;
            nPostSweeps    2;
            nFinestSweeps  2;
            cacheAgglomeration false;
            nCellsInCoarsestLevel 10;
            agglomerator    faceAreaPair;
            mergeLevels    1;
        }
        tolerance      1e-05;
        relTol          0;
        maxIter        10;
    }


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

    "(U|k|omega)Final"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance      1e-07;
        relTol          0;
    }
}

PIMPLE
{
    momentumPredictor no;
    nCorrectors    2;
    nNonOrthogonalCorrectors 0;
    nAlphaCorr      1;
    nAlphaSubCycles 2;
    cAlpha          1;  // Alternative 0
}


And alpha (dont expect problems here)
Code:

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

dimensions      [0 0 0 0 0 0 0];

internalField  uniform 0;

boundaryField
{
    walls
    {
        type            zeroGradient;
    }

    //Durchströmt
   
    inlet
    {
        type            fixedValue;
        value          uniform 1;
    }
       
        atmosphaere_CATIASTL
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value          uniform 0;
    }
   
    //Alles was umströmt wird

    tank_CATIASTL
    {
        type            zeroGradient;
    }
   
    efr_CATIASTL
    {
        type            zeroGradient;
    }

    entl_CATIASTL
    {
        type            zeroGradient;
    }
   
    pistole_vorn_CATIASTL
    {
        type            zeroGradient;
    }
   
    sd_CATIASTL
    {
        type            zeroGradient;
    }
   
    pistole_hinten_CATIASTL
    {
        type            zeroGradient;
    }
   
    innenteile_CATIASTL
    {
        type            zeroGradient;
    }
   

   
    flv_CATIASTL
    {
        type            zeroGradient;
    }
   
      einsatz_efr_CATIASTL
    {
        type            zeroGradient;
    }
   
    verschlusskopf_CATIASTL
    {
        type            zeroGradient;
    }

}

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


simpomann December 12, 2012 13:11

Hey,

I checked a lamniar case. Same error but at a different timestep
Code:

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62932e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62932e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.065834, Final residual = 0.000307154, No Iterations 2
time step continuity errors : sum local = 2.85125e-12, global = -6.3837e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.02062, Final residual = 3.97336e-08, No Iterations 7
time step continuity errors : sum local = 3.9867e-16, global = -5.09275e-19, cumulative = -3.90402e-05
ExecutionTime = 15548.4 s  ClockTime = 15612 s

Courant Number mean: 2.84517e-08 max: 0.867998
Interface Courant Number mean: 5.5157e-09 max: 0.0944208
deltaT = 4.14157e-15
Time = 0.202400597327766

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0662251, Final residual = 0.000298758, No Iterations 2
time step continuity errors : sum local = 2.77135e-12, global = -4.86472e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0204396, Final residual = 3.2508e-08, No Iterations 7
time step continuity errors : sum local = 3.2613e-16, global = -4.25822e-19, cumulative = -3.90402e-05
ExecutionTime = 15550.6 s  ClockTime = 15614 s

Courant Number mean: 2.84501e-08 max: 0.867991
Interface Courant Number mean: 5.51484e-09 max: 0.0943805
deltaT = 3.81716e-15
Time = 0.202400597327769

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0657749, Final residual = 0.00031065, No Iterations 2
time step continuity errors : sum local = 2.88026e-12, global = -6.24636e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0206025, Final residual = 3.60495e-08, No Iterations 7
time step continuity errors : sum local = 3.61229e-16, global = -3.50359e-19, cumulative = -3.90402e-05
ExecutionTime = 15552.7 s  ClockTime = 15616 s

Courant Number mean: 2.84444e-08 max: 0.867923
Interface Courant Number mean: 5.6614e-09 max: 0.094348
deltaT = 3.51843e-15
Time = 0.202400597327773

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0661369, Final residual = 0.000300055, No Iterations 2
time step continuity errors : sum local = 2.77963e-12, global = -4.82001e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0204266, Final residual = 3.17341e-08, No Iterations 7
time step continuity errors : sum local = 3.17909e-16, global = -3.68428e-19, cumulative = -3.90402e-05
ExecutionTime = 15554.9 s  ClockTime = 15619 s

Courant Number mean: 2.8439e-08 max: 0.867925
Interface Courant Number mean: 5.317e-09 max: 0.094304
deltaT = 3.24307e-15
Time = 0.202400597327776

MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
MULES: Solving for alpha1
Liquid phase volume fraction = 0.00157538  Min(alpha1) = -1.62931e-21  Max(alpha1) = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0657622, Final residual = 0.00031015, No Iterations 2
time step continuity errors : sum local = 2.87193e-12, global = -6.3359e-15, cumulative = -3.90402e-05
GAMGPCG:  Solving for p_rgh, Initial residual = 0.0205844, Final residual = 3.79572e-08, No Iterations 7
time step continuity errors : sum local = 3.79847e-16, global = -3.64393e-19, cumulative = -3.90402e-05
ExecutionTime = 15557.1 s  ClockTime = 15621 s

Courant Number mean: 2.84311e-08 max: 0.867857
Interface Courant Number mean: 5.31414e-09 max: 0.0942674
deltaT = 2.9895e-15
Time = 0.202400597327779

[3] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2  in "/lib64/libc.so.6"
[3] #3  void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #4  at MULES.C:0
[3] #5  Foam::MULES::explicitSolve(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #6 
[3]  in "/opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/bin/interFoam"
[3] #7  __libc_start_main in "/lib64/libc.so.6"
[3] #8 
[3]  at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
[c103-06:28794] *** Process received signal ***
[c103-06:28794] Signal: Floating point exception (8)
[c103-06:28794] Signal code:  (-6)
[c103-06:28794] Failing at address: 0x35c60000707a
[c103-06:28794] [ 0] /lib64/libc.so.6(+0x34e10) [0x7f3d0d63ae10]
[c103-06:28794] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7f3d0d63ad95]
[c103-06:28794] [ 2] /lib64/libc.so.6(+0x34e10) [0x7f3d0d63ae10]
[c103-06:28794] [ 3] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES7limiterINS_17geometricOneFieldENS_9zeroFieldES3_EEvRNS_5FieldIdEERKT_RKNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNSA_IdNS_13fvsPatchFieldENS_11surfaceMeshEEESK_RKT0_RKT1_ddi+0x1068) [0x7f3d0f268fa8]
[c103-06:28794] [ 4] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(+0x7c007e) [0x7f3d0f25c07e]
[c103-06:28794] [ 5] /opt/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam5MULES13explicitSolveERNS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERS8_dd+0x80) [0x7f3d0f25c310]
[c103-06:28794] [ 6] interFoam() [0x42d35e]
[c103-06:28794] [ 7] /lib64/libc.so.6(__libc_start_main+0xed) [0x7f3d0d62723d]
[c103-06:28794] [ 8] interFoam() [0x425119]
[c103-06:28794] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 28794 on node c103-06 exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

As far as I understand, there is an error during the alpha calculation. But what can cause it? If the turbulence model is allright, checkMesh says okay.... where can I start searching my answers?
Could it be this part from my fvSolution?
PIMPLE { momentumPredictor no; nCorrectors 2; nNonOrthogonalCorrectors 0; nAlphaCorr 1; nAlphaSubCycles 2; cAlpha 1; // Alternative 0 }

matejmuller February 17, 2014 18:06

Hey Simon!

Have you perhaps solved this issue? I'm also trying to run interFoam with the kOmega model and also don't know which solvers and schemes to use, since the tutorial of kOmegaSST exists only for the LTSInterFoam. Have you figured out the right solvers and schemes?

thanks for the help,

Matej


All times are GMT -4. The time now is 16:26.