CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

twoPhaseEulerFoam kEpsilon sigFpe error!!

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

Like Tree2Likes
  • 2 Post By akm

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 22, 2017, 04:43
Exclamation twoPhaseEulerFoam kEpsilon sigFpe error!!
  #1
akm
New Member
 
Join Date: Jan 2010
Location: Netherlands
Posts: 28
Rep Power: 16
akm is on a distinguished road
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<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
luccy and stamufa like this.
akm is offline   Reply With Quote

Old   June 18, 2017, 06:35
Default
  #2
New Member
 
farzad khallaghi
Join Date: Dec 2015
Posts: 10
Rep Power: 10
farzadpolytechnic is on a distinguished road
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?
farzadpolytechnic is offline   Reply With Quote

Old   May 2, 2019, 06:32
Default
  #3
New Member
 
Mustafa
Join Date: Jun 2015
Location: Aachen
Posts: 26
Rep Power: 11
stamufa is on a distinguished road
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<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 ??:?
stamufa is offline   Reply With Quote

Reply

Tags
kepsilon, multiphase flow, sigfpe, twophaseeulerfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[OpenFOAM] Native ParaView Reader Bugs tj22 ParaView 270 January 4, 2016 11:39
OpenFOAM without MPI kokizzu OpenFOAM Installation 4 May 26, 2014 09:17
Compile problem ivanyao OpenFOAM Running, Solving & CFD 1 October 12, 2012 09:31
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50
Compiling problems with hello worldC fw407 OpenFOAM Installation 21 January 6, 2008 17:38


All times are GMT -4. The time now is 22:50.