|
[Sponsors] |
September 9, 2010, 09:56 |
simpleFoam error - calculation stopped
|
#1 |
Member
Martin
Join Date: Aug 2010
Location: Germany
Posts: 55
Rep Power: 15 |
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 |
|
September 9, 2010, 10:29 |
|
#2 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 15 |
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. |
|
September 9, 2010, 10:46 |
|
#3 |
Member
Martin
Join Date: Aug 2010
Location: Germany
Posts: 55
Rep Power: 15 |
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..!? |
|
September 9, 2010, 10:50 |
|
#4 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 15 |
No problem!
|
|
September 9, 2010, 10:59 |
|
#5 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 15 |
I sended you my mail by private message.
|
|
September 10, 2010, 03:43 |
|
#6 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 15 |
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! |
|
September 10, 2010, 03:59 |
|
#7 |
Member
Martin
Join Date: Aug 2010
Location: Germany
Posts: 55
Rep Power: 15 |
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 |
|
September 10, 2010, 04:30 |
|
#8 |
New Member
mikael
Join Date: Jun 2010
Posts: 18
Rep Power: 15 |
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.
|
|
|
|
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 |