# interFoam/kOmegaSST tank filling with printStackError/Mules

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

 December 12, 2012, 07:39 interFoam/kOmegaSST tank filling with printStackError/Mules #1 Member   Simon Arne Join Date: May 2012 Posts: 42 Rep Power: 14 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::Field&, Foam::geometricOneField const&, Foam::GeometricField const&, Foam::GeometricField const&, Foam::GeometricField 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&, Foam::GeometricField const&, Foam::GeometricField&, 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

 December 12, 2012, 07:40 #2 Member   Simon Arne Join Date: May 2012 Posts: 42 Rep Power: 14 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; } } // ************************************************************************* //```

 December 12, 2012, 12:11 #3 Member   Simon Arne Join Date: May 2012 Posts: 42 Rep Power: 14 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::Field&, Foam::geometricOneField const&, Foam::GeometricField const&, Foam::GeometricField const&, Foam::GeometricField 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&, Foam::GeometricField const&, Foam::GeometricField&, 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 } Last edited by simpomann; December 12, 2012 at 12:43.

 February 17, 2014, 17:06 #4 New Member   Matej Muller Join Date: Oct 2011 Location: Slovenia Posts: 25 Rep Power: 14 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