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

potentialFoam giving strange error when initialising simpleFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 12, 2013, 07:40
Default potentialFoam giving strange error when initialising simpleFoam
  #1
New Member
 
Anonymous
Join Date: Aug 2013
Location: Europe
Posts: 24
Rep Power: 12
maero21 is on a distinguished road
Dear OpenFOAM users,

I want to run a 3D simulation of a double element wing with endplates. The mesh is around 5 million cells. However when I run simpleFoam it will immediately diverge. Therefore I want to initialize the simpleFoam run by first running potentialFoam. However, I keep getting a very strange error which does not give me any information:

Code:
Calculating potential flow
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::adjustPhi(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4  
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/potentialFoam"
#5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6  
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/potentialFoam"
Floating point exception (core dumped)
I have attached my fvSolution, fvSchemes and controlDict.

I hope someone has had this error before and knows what to do.

By the way, I am running OpenFoam 2.1.1.
Attached Files
File Type: txt controlDict .txt (3.4 KB, 31 views)
File Type: txt fvSchemes .txt (1.7 KB, 24 views)
File Type: txt fvSolution .txt (2.1 KB, 35 views)

Last edited by wyldckat; October 12, 2013 at 15:36. Reason: Added [CODE][/CODE]
maero21 is offline   Reply With Quote

Old   October 12, 2013, 15:41
Default
  #2
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Greetings maero21,

"sigFpe" means this: http://en.wikipedia.org/wiki/SIGFPE#SIGFPE
Quote:
The SIGFPE signal is sent to a process when it executes an erroneous arithmetic operation, such as division by zero (the FPE stands for floating point error).
In the line below that it says that it was running the method "adjustPhi". Therefore, something very wrong is happening with the pressure values and/or flow velocity, since it triggered a division by zero.

This usually has to do with bad boundary conditions or a damaged mesh. You can use checkMesh to ascertain the quality of the mesh:
Code:
checkMesh -allTopology -allGeometry
Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   October 15, 2013, 12:47
Default
  #3
New Member
 
Anonymous
Join Date: Aug 2013
Location: Europe
Posts: 24
Rep Power: 12
maero21 is on a distinguished road
Thank you! Indeed the mesh failed 1 mesh check. It says that it has found "cells with small determinant". I am guessing that is the reason it divides something by 0 in the solving.
How can I see those small determinant cells? I have generated the mesh in Pointwise, but apparently Pointwise is not good with 3D meshes? That would be strange, right?
maero21 is offline   Reply With Quote

Old   October 15, 2013, 17:19
Default
  #4
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Quick answer:
Code:
foamToVTK -cellSet nameOfCellSetIndicatedByCheckMesh
__________________
wyldckat is offline   Reply With Quote

Reply

Tags
initialization, potentialfoam, simplefoam stability

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How to run potentialFoam and simpleFoam together . vmsandip2011 OpenFOAM Running, Solving & CFD 11 April 2, 2021 11:56
simpleFoam: simple 1-D channel flow, yet very strange convergence behavior kishpishar OpenFOAM Running, Solving & CFD 2 June 20, 2013 14:55
simpleFoam: strange error samiam1000 OpenFOAM 7 December 11, 2012 04:52
simpleFOAM not giving steadystate solution aerospain OpenFOAM 5 August 6, 2012 08:50
BC for simpleFoam from potentialFoam results Geon-Hong OpenFOAM Running, Solving & CFD 0 April 5, 2011 23:23


All times are GMT -4. The time now is 06:39.