|
[Sponsors] |
May 24, 2010, 03:52 |
interFoam Error
|
#1 | |
Member
Join Date: Mar 2009
Location: Sydney, New South Wales, Australia
Posts: 42
Rep Power: 17 |
Hi Foamers,
I am attempting to model a device that I have meshed using snappyHexMesher and run using the interFoam solver, however I get a floating point exception error at the very first iteration. The error looks like this: Quote:
The error at #3 looks like it may be some kind of field problem. Could it be that the fields are not initialising correctly using setFields? If anybody has any ideas, I would love to hear them! R |
||
May 24, 2010, 11:35 |
|
#2 |
Senior Member
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 219
Rep Power: 21 |
Did you run checkMesh on the generated mesh to make sure nothing crazy is happening with the mesh? I would also indeed first check that you have initialized your alpha1 field correctly--make sure to patch a little into the domain on any inlet boundaries you may have--this helps the solution get started. If you are still having problems, can you post your fvSolution and fvSchemes settings?
Just a couple thoughts that may or may not help... |
|
May 25, 2010, 03:38 |
|
#3 | ||||
Member
Join Date: Mar 2009
Location: Sydney, New South Wales, Australia
Posts: 42
Rep Power: 17 |
Quote:
Thanks for replying Kent. checkMesh does give the all clear, however also identifies quite a few ployhedra. Not sure if this is suitable for use with interFoam. Quote:
Quote:
Quote:
R |
|||||
May 25, 2010, 05:04 |
|
#4 |
Senior Member
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,266
Blog Entries: 1
Rep Power: 24 |
if u use groovyBC for boundary condition add following statement in to controlDic;
libs ( "libOpenFOAM.so" "libgroovyBC.so" ) ; even though check ur boundary condition and initial condition and change them to none Zero value! maybe work |
|
May 26, 2010, 11:07 |
|
#5 |
Senior Member
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 219
Rep Power: 21 |
Polyhedral cells are not a problem--they actually should be more accurate. You haven't said what kind of problem you are trying to solve. A general description would be helpful for a diagnosis. Did you check your initialization for volume fraction and make sure it is doing what you think it is? Since you are getting a floating point error perhaps you are dividing by zero somewhere.
|
|
May 31, 2010, 02:48 |
|
#6 | |
Member
Join Date: Mar 2009
Location: Sydney, New South Wales, Australia
Posts: 42
Rep Power: 17 |
Quote:
Hi Kent, Thanks for getting back to me. I seem to have got it working now - It appeared to be some kind of scaling problem when I used the transformPoints utility - I went back to the beginning and started over, and this seems to be the sticking point, but it looks like I have it working now. Hopefully I won't have any more problems! Cheers, R |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh | gschaider | OpenFOAM Community Contributions | 300 | October 29, 2014 19:00 |
OpenFOAM on MinGW crosscompiler hosted on Linux | allenzhao | OpenFOAM Installation | 127 | January 30, 2009 20:08 |
[Netgen] Compiling Netgen on Fedora Core is driving me crazy | jango | OpenFOAM Meshing & Mesh Conversion | 3 | November 9, 2007 14:29 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 18:51 |
user defined function | cfduser | CFX | 0 | April 29, 2006 11:58 |