CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   interFoam/kOmegaSST tank filling with printStackError/Mules (https://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 05:02.