interFoam Error
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 |
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... |
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 |
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 |
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.
|
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 |
All times are GMT -4. The time now is 20:06. |