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

Problem with using icoFoam for flow in a straight tube

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 10, 2011, 06:01
Question Problem with using icoFoam for flow in a straight tube
  #1
New Member
 
Join Date: Feb 2011
Posts: 2
Rep Power: 0
edificelp is on a distinguished road
Hi:

I am trying to use OpenFoam to simulate incompressible, laminar flow of Newtonian fluids in a straight tube. And the 3-d mesh is converted from .msh file of Gmsh. After setting up boundary condition and running the simulation, the courant number mean increases so fast and becomes nearly infinite. Then the simulation stops with the information below:

Code:
Courant Number mean: 4.47866e+122 max: 1.17861e+124
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#2   in "/lib/libc.so.6"
#3  Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#4  
 in "/opt/openfoam171/applications/bin/linux64GccDPOpt/icoFoam"
#5  
 in "/opt/openfoam171/applications/bin/linux64GccDPOpt/icoFoam"
#6  __libc_start_main in "/lib/libc.so.6"
#7  
 in "/opt/openfoam171/applications/bin/linux64GccDPOpt/icoFoam"
[1]+  Done                    kate ~/.bashrc
Floating point exception
The file of this simulation has been attached. Could someone help me solve this problem?
Thanks a lot!
Attached Files
File Type: gz singleTubeFoam.tar.gz (62.8 KB, 16 views)
edificelp is offline   Reply With Quote

Old   February 10, 2011, 09:27
Default
  #2
Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 78
Rep Power: 7
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
A blind guess: Try to reduce your time step to suffice the CFL condition.
jhoepken is offline   Reply With Quote

Old   February 10, 2011, 14:07
Default
  #3
New Member
 
Jean El-Hajal
Join Date: Jun 2010
Location: Ulm
Posts: 16
Rep Power: 6
Jean El-Hajal is on a distinguished road
Hi,

it seems your problem came from the fvSolution and fvSchemes file. Replace them with the original ones from the cavity tutorial.

However the mesh is very bad !

Regards,

Jean
Jean El-Hajal is offline   Reply With Quote

Old   February 10, 2011, 17:23
Default
  #4
New Member
 
Join Date: Feb 2011
Posts: 2
Rep Power: 0
edificelp is on a distinguished road
Thanks a lot! After reading user guide, the incompressible steady state flow just can be simulated by simpleFoam. If icoFoam is used with steadyState, the solution will not converge.

By the way, Jean EI-Hajal, could you tell me how do you determine the mesh is bad or not? This mesh is generated by Gmsh automatically with the characteristic length. I feel headache to check this 3-d mesh.........
edificelp is offline   Reply With Quote

Old   February 11, 2011, 03:27
Default
  #5
Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 78
Rep Power: 7
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
Run checkMesh and post the output. That tool checks the quality of the mesh and prints some statistics to your console.
jhoepken 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
analytical solution of 1D shock tube problem ma Main CFD Forum 8 July 12, 2013 10:34
laplacian1%7cAUp is undefined in dictionary problem with icoFoam derath OpenFOAM Pre-Processing 3 June 14, 2013 06:24
problem when running icoFoam on a complex shape flow field wendywu OpenFOAM 1 May 20, 2009 23:40
VOF problem, liquid film thickness Pedro24 FLUENT 0 May 13, 2009 05:57
Shock tube problem ger Main CFD Forum 1 September 24, 2003 07:41


All times are GMT -4. The time now is 05:28.