[akm]

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?

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<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8  Foam::fvMatrix<double>::solve() at ??:?
#9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
#10  Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::PhaseCompressibleTurbulenceModel<Foam::phaseModel> > > >::correct() at ??:?
#11  ? at ??:?
#12  __libc_start_main in "/lib64/libc.so.6"
#13  ? at ??:?
Floating exception

[farzadpolytechnic]

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?

[stamufa]

Get another version of the same error for nuEff. Any ideas?

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<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[15] #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fv
PatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&
, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[15] #5  Foam::CompressibleTurbulenceModel<Foam::fluidThermo>::nu() const at ??:?
[15] #6  Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam
::fluidThermo> > > >::nuEff() const at ??:?
[15] #7  Foam::linearViscousStress<Foam::RASModel<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::Compre
ssibleTurbulenceModel<Foam::fluidThermo> > > > >::divDevRhoReff(Foam::GeometricField<Foam::Vector<double>, Fo
am::fvPatchField, Foam::volMesh>&) const at ??:?
[15] #8  ? at ??:?
[15] #9  __libc_start_main in "/lib64/libc.so.6"
[15] #10  ? at ??:?