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

floating point exception with pimpleDyMFoam

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 27, 2016, 07:58
Default floating point exception with pimpleDyMFoam
  #1
New Member
 
Join Date: Jul 2013
Posts: 3
Rep Power: 13
Elena89 is on a distinguished road
Hallo everybody,

I'm trying to simulate incompressible flow through turbomachinery with OpenFOAM using pimpleDyMFoam coupled with kOmegaSSt for turbulence modelling.
I've stuck pretty close to the propeller tutorial as to the boundary conditions adjustments.

When I run the case (in parallel as well as in serial) the first approximately 40-50 timesteps work just fine. Then all of a sudden the velocity gets very high and I get a floating point exception (please find log file attached). At least I think that the velocity is the problem...

I've been trying to figure out what causes the problem for several days now and have been trying a lot, but without any success.

I'm pretty sure that it's not a problem with my mesh since check mesh tells me "Ok". It could be something with the fvSchemes or fvSolutions, couldn't it?

Does anyone have any ideas where the problem / error might be? Any help would be greatly appreciated. Thanks in advance

Best regards

Elena
Attached Files
File Type: txt log.txt (103.1 KB, 23 views)
File Type: txt fvSolution.txt (2.6 KB, 13 views)
File Type: txt fvSchemes.txt (1.5 KB, 10 views)
Elena89 is offline   Reply With Quote

Old   March 13, 2016, 09:51
Default
  #2
Senior Member
 
Artur's Avatar
 
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 372
Rep Power: 20
Artur will become famous soon enough
Hi,

Your problem stems from the fact that continuity errors get through the roof pretty quickly:
Code:
time step continuity errors : sum local = 3.83444, global = -1.13213, cumulative = -1.60869
After which the solution goes bonkers, which results in the Co exploding, time step decreasing to crazy values and the simulation ultimately failing:
Code:
Courant Number mean: 39095.8 max: 1.52799e+11
deltaT = 1.09076e-14
--> FOAM Warning : 
    From function Time::operator++()
    in file db/Time/Time.C at line 1061
    Increased the timePrecision from 6 to 7 to distinguish between timeNames at time 0.431667
Nothing seems to be particularly alarming in your fv files, although 10 correctors seems quite a lot, given most of them don't execute at all (0 iterations in the log). My wild guess would be either the BCs or the mesh in your case causing the issue. If not you could try a smaller time step. Plus, you get some interesting warnings about "switchableNotImpleemnted", which I suspect might have something to do with the problems as well.

All the best,

A
Artur is offline   Reply With Quote

Reply

Tags
floating point exception, pimpledymfoam


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
icoFoam floating point exception (8) leizhao512 OpenFOAM Running, Solving & CFD 7 November 1, 2018 12:43
Floating point exception with pimpleDyMFoam ebah6 OpenFOAM Running, Solving & CFD 9 November 1, 2017 06:58
Inlet Velocity Profile BC - Floating Point exception during solution initialization Janshi STAR-CCM+ 4 March 14, 2012 11:21
Pipe flow in settlingFoam floating point exception jochemvandenbosch OpenFOAM Running, Solving & CFD 4 February 16, 2012 04:24
block-structured mesh for t-junction Robert@cfd ANSYS Meshing & Geometry 20 November 11, 2011 05:59


All times are GMT -4. The time now is 14:54.