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

simpleFoam error - calculation stopped

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 9, 2010, 09:56
Default simpleFoam error - calculation stopped
  #1
Member
 
Martin
Join Date: Aug 2010
Location: Germany
Posts: 55
Rep Power: 6
RugbyGandalf is on a distinguished road
hello,

i tried to calculate the flow through two combined cylinders...
meshing tool was netgen.
it is one of the example geometries of netgen...

but simpleFoam solver always breaks the operation after 23 steps of calculation.
it points out the following errors in the terminal windows:

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/opt/openfoam170/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#7
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/simpleFoam"
#8 __libc_start_main in "/lib/libc.so.6"
#9
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/simpleFoam"
Floating point exception


could someone help me, how to get this case run?

i also postet the "0" time directory and system folder
constant folder was a bit to large

thank you very much in advance

Martin
Attached Files
File Type: zip zylinder.zip (6.2 KB, 4 views)
RugbyGandalf is offline   Reply With Quote

Old   September 9, 2010, 10:29
Default
  #2
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 6
keul is on a distinguished road
Hi,

I can have a look at your problem but I need your constant folder. You could reduce the folder size if you delete all the files in polymesh exept the blockMeshDict file.
keul is offline   Reply With Quote

Old   September 9, 2010, 10:46
Default
  #3
Member
 
Martin
Join Date: Aug 2010
Location: Germany
Posts: 55
Rep Power: 6
RugbyGandalf is on a distinguished road
thank you very much...

but i do not have a blockmesh file...

the mesh was created by netgen...

but of course i could send it via email, if that would be okay..!?
RugbyGandalf is offline   Reply With Quote

Old   September 9, 2010, 10:50
Default
  #4
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 6
keul is on a distinguished road
No problem!
keul is offline   Reply With Quote

Old   September 9, 2010, 10:59
Default
  #5
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 6
keul is on a distinguished road
I sended you my mail by private message.
keul is offline   Reply With Quote

Old   September 10, 2010, 03:43
Default
  #6
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 6
keul is on a distinguished road
Hi RugbyGandalf!

I launched your case with epsilon = 1 and it works fine. When you work with a turbulent model the solution convergence depends a lot of your initial conditions on k and epsilon. So maybe try other values until you find good ones.

Good luck!
keul is offline   Reply With Quote

Old   September 10, 2010, 03:59
Default
  #7
Member
 
Martin
Join Date: Aug 2010
Location: Germany
Posts: 55
Rep Power: 6
RugbyGandalf is on a distinguished road
thank you very much

it is working great

is there a "rule" for setting epsilon and k?

what are these values saying exactly?
i am very new to cfd and openFoam - so i do not already know about all the values
RugbyGandalf is offline   Reply With Quote

Old   September 10, 2010, 04:30
Default
  #8
New Member
 
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 6
keul is on a distinguished road
I'm not an expert in turbulent model but k and epsilon depend of initial velocity, characteristic length and Reynolds number. Maybe you should have a look at some explanations you could find on Wikipedia to fix your initial conditions.
keul 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
Laminar simpleFoam and inviscid simpleFoam herenger OpenFOAM Running, Solving & CFD 7 July 11, 2013 06:27
MPI Error - simpleFoam - Floating Point Exception scott OpenFOAM Running, Solving & CFD 3 April 13, 2012 16:34
simpleFoam ddt Euler ? Mo-ITB OpenFOAM Running, Solving & CFD 2 June 12, 2010 13:36
Naca0012 k-e mpirun gives fpe whereas simpleFoam not Pierpaolo OpenFOAM 1 May 8, 2010 03:08
Error running simpleFoam in parallel skabilan OpenFOAM Running, Solving & CFD 2 August 29, 2008 09:42


All times are GMT -4. The time now is 10:32.