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

K-eps Error

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 6, 2013, 15:20
Default K-eps Error
  #1
New Member
 
SV
Join Date: May 2009
Posts: 15
Rep Power: 8
Soheyl is on a distinguished road
When I'm trying to use the k-eps model with my MRFSimpleFOAM problem, I get the following error:

Quote:
Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 1 average: 1
bounding epsilon, min: 0 max: 20 average: 20
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#5 at kEpsilon.C:0
#6 Foam::incompressible::RASModels::kEpsilon::kEpsilo n(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::kEp silon>::New(Foam::GeometricField<Foam::Vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#8 Foam::incompressible::RASModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#9
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/MRFSimpleFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/MRFSimpleFoam"
Floating point exception (core dumped)
Any idea what is going on?
Soheyl is offline   Reply With Quote

Old   March 6, 2013, 18:37
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Soheyl,

Quote:
Code:
bounding k, min: 0 max: 1 average: 1
bounding epsilon, min: 0 max: 20 average: 20
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
About SIGFPE: http://en.wikipedia.org/wiki/SIGFPE#SIGFPE
In essence, it looks like you have incorrectly initialized the fields for "k" and "epsilon" with 0 values. Those two fields should be initiated with values larger than 0.

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   March 6, 2013, 20:04
Default
  #3
New Member
 
SV
Join Date: May 2009
Posts: 15
Rep Power: 8
Soheyl is on a distinguished road
Thank you very much Bruno. This was my 0/k file:

Quote:
internalField uniform 1;

boundaryField
{
RotorWalls
{
type kqRWallFunction;
value uniform 0;
}

Walls
{
type kqRWallFunction;
value uniform 0;
}

Inlet
{
type fixedValue;
value uniform 0;

}

Outlet
{
type zeroGradient;
}
}
and this was my 0/epsilon:

Quote:
internalField uniform 20;

boundaryField
{
RotorWalls
{
type epsilonWallFunction;
value uniform 0;
}

Walls
{
type epsilonWallFunction;
value uniform 0;
}

Inlet
{
type fixedValue;
value uniform 0;

}

Outlet
{
type zeroGradient;
}
}
My understanding was that the "initial" internal values come from internalField lines, which are nonzero. My boundary values were zero though. By changing them to nonzero values the problem is solved.

I'm curious to know why this happens though.
Soheyl is offline   Reply With Quote

Old   March 7, 2013, 09:18
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Quote:
Originally Posted by Soheyl View Post
I'm curious to know why this happens though.
All I can remember is that they cannot be zero. You can check the math behind the k-epsilon model: http://www.cfd-online.com/Wiki/K-epsilon_models
wyldckat is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh gschaider OpenFOAM 300 October 29, 2014 19:00
c++ libraries and solver compiling vaina74 OpenFOAM Installation 13 February 3, 2012 18:43
Saving ParaFoam views and case sail OpenFOAM Paraview & paraFoam 9 November 25, 2011 16:46
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51
user defined function cfduser CFX 0 April 29, 2006 10:58


All times are GMT -4. The time now is 17:13.