# twoPhaseEulerFoam kEpsilon sigFpe error!!

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

 March 22, 2017, 05:43 twoPhaseEulerFoam kEpsilon sigFpe error!! #1 New Member   Join Date: Jan 2010 Location: Netherlands Posts: 28 Rep Power: 15 I am trying to run twoPhaseEulerFoam for a turbulent multiphase flow case with kEpsilon model, but I keep getting sigFpe error (I think) on the first time step. I have made sure that the initial k, epsilon, alphat and nut fields are non-zero. If I switch the simulation to laminar, it proceeds just fine; but if I switch either one or both phases to turbulent, I get the error below. (Eventually, I want only one phase as turbulent) I have even tried running the simulation laminar for an initial time period, and then later introducing the non-zero turbulent fields (k, eps,...) and switching the simulation to turbulent. I receive the same error. If I understand correctly, this seems to be a divide-by-zero error, but I do not understand where/how is it dividing by zero. Does anyone have suggestions where to look? Code: ```PIMPLE: iteration 1 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.260475 Min(alpha.air) = -9.79204e-15 Max(alpha.air) = 1 Constructing momentum equations min T.air 300 min T.water 300 Selecting patchDistMethod meshWave GAMG: Solving for p_rgh, Initial residual = 0.000842622, Final residual = 4.12298e-06, No Iterations 7 GAMG: Solving for p_rgh, Initial residual = 0.000746298, Final residual = 3.76317e-06, No Iterations 7 GAMG: Solving for p_rgh, Initial residual = 0.00074733, Final residual = 4.01164e-06, No Iterations 7 PIMPLE: iteration 2 MULES: Solving for alpha.air MULES: Solving for alpha.air alpha.air volume fraction = 0.260474 Min(alpha.air) = -2.1628e-12 Max(alpha.air) = 1 Constructing momentum equations min T.air 300 min T.water 300 GAMG: Solving for p_rgh, Initial residual = 0.00073935, Final residual = 7.29952e-06, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 0.00073454, Final residual = 7.21356e-06, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 0.00073508, Final residual = 9.83213e-07, No Iterations 11 #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib64/libc.so.6" #3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field&, Foam::lduMatrix const&, Foam::Field const&, Foam::FieldField const&, Foam::UPtrList const&, unsigned char, int) at ??:? #4 Foam::GaussSeidelSmoother::smooth(Foam::Field&, Foam::Field const&, unsigned char, int) const at ??:? #5 Foam::smoothSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const at ??:? #6 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) at ??:? #7 Foam::fvMatrix::solve(Foam::dictionary const&) at ??:? #8 Foam::fvMatrix::solve() at ??:? #9 Foam::SolverPerformance Foam::solve(Foam::tmp > const&) at ??:? #10 Foam::RASModels::kEpsilon > > >::correct() at ??:? #11 ? at ??:? #12 __libc_start_main in "/lib64/libc.so.6" #13 ? at ??:? Floating exception``` stamufa likes this.

 June 18, 2017, 07:35 #2 New Member   farzad khallaghi Join Date: Dec 2015 Posts: 10 Rep Power: 9 hi does anybody find the answer? i have same problem. if use default turbulence model, the solver run correctly. but when i want change the turbulence model to the other possible models i find this error. so please help me how can i solve this?

 May 2, 2019, 07:32 #3 New Member   Mustafa Join Date: Jun 2015 Location: Aachen Posts: 26 Rep Power: 10 Get another version of the same error for nuEff. Any ideas? Code: ```PIMPLE: iteration 2 [15] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [15] #1 Foam::sigFpe::sigHandler(int) at ??:? [15] #2 ? in "/lib64/libc.so.6" [15] #3 Foam::divide(Foam::Field&, Foam::UList const&, Foam::UList const&) at ??:? [15] #4 Foam::tmp > Foam::operator/(Foam::tmp > const& , Foam::GeometricField const&) at ??:? [15] #5 Foam::CompressibleTurbulenceModel::nu() const at ??:? [15] #6 Foam::RASModel > > >::nuEff() const at ??:? [15] #7 Foam::linearViscousStress > > > >::divDevRhoReff(Foam::GeometricField, Fo am::fvPatchField, Foam::volMesh>&) const at ??:? [15] #8 ? at ??:? [15] #9 __libc_start_main in "/lib64/libc.so.6" [15] #10 ? at ??:?```

 Tags kepsilon, multiphase flow, sigfpe, twophaseeulerfoam