|
[Sponsors] |
March 22, 2017, 04:43 |
twoPhaseEulerFoam kEpsilon sigFpe error!!
|
#1 |
New Member
Join Date: Jan 2010
Location: Netherlands
Posts: 28
Rep Power: 16 |
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 |
|
June 18, 2017, 06:35 |
|
#2 |
New Member
farzad khallaghi
Join Date: Dec 2015
Posts: 10
Rep Power: 10 |
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, 06:32 |
|
#3 |
New Member
Mustafa
Join Date: Jun 2015
Location: Aachen
Posts: 26
Rep Power: 11 |
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 ??:? |
|
Tags |
kepsilon, multiphase flow, sigfpe, twophaseeulerfoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |