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

URGENT! Error with icoFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 7, 2012, 14:52
Default URGENT! Error with icoFoam
  #1
New Member
 
Mattia
Join Date: Oct 2012
Location: Milan
Posts: 28
Rep Power: 4
gara1988 is on a distinguished road
Hello, I did a mesh around an airfoil. BlockMesh has gone well.

When I run icoFoam I have this error:


Starting time loop

Time = 0.001

Courant Number mean: 0.0568875 max: 7.3117e+297
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 void Foam::fvc::surfaceIntegrate<Foam::Vector<double> >(Foam::Field<Foam::Vector<double> >&, Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::tmp<Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::tmp<Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::tmp<Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<Foam::Vector<double> >(Foam::tmp<Foam::GeometricField<Foam::Vector<doub le>, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::gaussLaplacianScheme<Foam::Vector<double >, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/icoFoam"
#9
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/icoFoam"
#10
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/icoFoam"
#11
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/icoFoam"
#12 __libc_start_main in "/lib/libc.so.6"
#13
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/icoFoam"
Floating point exception

Is Courant number the problem? But max courant number is terrible high!!

Thanks for the reply!

Bye
Mattia
gara1988 is offline   Reply With Quote

Old   October 7, 2012, 15:06
Default
  #2
New Member
 
Mattia
Join Date: Oct 2012
Location: Milan
Posts: 28
Rep Power: 4
gara1988 is on a distinguished road
I have seen tha,t in other post, Courant number is indicated as the problem. Here I have C number equal to 7.3117e+297, it is terrible high. I don't know why it is so high.
gara1988 is offline   Reply With Quote

Old   October 7, 2012, 15:34
Default
  #3
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Mattia,

Since it's icoFoam, there are at least two reasons this error occurs:
  1. The mesh might be damaged. To check this, run:
    Code:
    checkMesh
  2. You didn't indicate anything about how you configured parameters in "controlDict". Which probably means that you changed nothing or very little in it. Therefore, I suggest you go read the User Guide once again, namely this section: 2.1.5 Increasing the mesh resolution
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   October 7, 2012, 15:42
Default
  #4
New Member
 
Mattia
Join Date: Oct 2012
Location: Milan
Posts: 28
Rep Power: 4
gara1988 is on a distinguished road
Thanks Bruno for the reply.

When I run checkMesh I have "failed 3 mesh checks".
What means that the mesh is damaged?

Best regards
Mattia.
gara1988 is offline   Reply With Quote

Old   October 7, 2012, 15:48
Default
  #5
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
I meant "damaged" in a generic way.
As you've seen, it failed 3 mesh checks, which means that there seems to be something wrong with the mesh... which could be considered as "damaged".

You'll have to diagnose yourself what checkMesh tells you is wrong and fix "blockMeshDict" so that checkMesh no longer complains.
wyldckat is offline   Reply With Quote

Old   October 7, 2012, 17:07
Default
  #6
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,077
Blog Entries: 4
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Mattia,

which error messages do you get? Some messages are not as critical as other.
And set your start time to 1e-6 or 1e-7 for starting your solution.

But first check your mesh again. Like bruno said, you have to improve your mesh quality if the errors are critical errors
Tobi is offline   Reply With Quote

Old   October 8, 2012, 02:33
Default
  #7
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hints:

- the solution fails when calculating a Laplacian, and the mesh fails 3 tests, so you should report what tests it fails (or the whole output of checkMesh), and most likely re-mesh until those errors disappear. Tweaking the solver is unlikely to help.

- Please, next time don't write your question is "urgent". All questions are "urgent" in the same way, and trying to get more attention usually results in the opposite.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   October 8, 2012, 03:36
Default
  #8
New Member
 
Mattia
Join Date: Oct 2012
Location: Milan
Posts: 28
Rep Power: 4
gara1988 is on a distinguished road
Ok. I'm sorry Alberto.

Thanks all for the reply!

Best regards.
Mattia
gara1988 is offline   Reply With Quote

Old   October 8, 2012, 03:59
Default
  #9
New Member
 
Mattia
Join Date: Oct 2012
Location: Milan
Posts: 28
Rep Power: 4
gara1988 is on a distinguished road
How can I undestand where it's the problem of mesh? For example what is the face that has zero area.

Sorry but I'm studying openfoam from a week.


Checking geometry...
Overall domain bounding box (-5 -5 0) (15 5 0)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (0 0 2.45443e-16) OK.
Max cell openness = 9.72873e-16 OK.
Max aspect ratio = 0 OK.
***Zero or negative face area detected. Minimum area: 0
<<Writing 28416 zero area faces to set zeroAreaFaces
Min volume = 1.33333e-300. Max volume = 2e-300. Total volume = 3.95333e-296. Cell volumes OK.
Mesh non-orthogonality Max: 90 average: 90
***Number of non-orthogonality errors: 39300.
<<Writing 39300 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
***Max skewness = 706.342, 42336 highly skew faces detected which may impair the quality of the results
<<Writing 42336 skew faces to set skewFaces
Coupled point location match (average 0) OK.

Failed 3 mesh checks.

End


Thanks.
Best regards.
Mattia
gara1988 is offline   Reply With Quote

Old   October 8, 2012, 08:27
Default
  #10
New Member
 
Mattia
Join Date: Oct 2012
Location: Milan
Posts: 28
Rep Power: 4
gara1988 is on a distinguished road
With checkMesh and doing the mesh better I solve the problem. Thanks a lot.

Mattia
gara1988 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
icoFoam crash with unreasonable velocity. Bylund OpenFOAM Running, Solving & CFD 2 November 20, 2011 21:48
Urgent - Polyflow - particle tracking visualization -Urgent shafaatht ANSYS 0 October 13, 2010 04:56
Density in icoFoam Densidad en icoFoam manuel OpenFOAM Running, Solving & CFD 8 September 22, 2010 04:10
Kubuntu uses dash breaks All scripts in tutorials platopus OpenFOAM Bugs 8 April 15, 2008 07:52
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58


All times are GMT -4. The time now is 12:03.