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

incorrect results with a tetrahedral mesh

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 20, 2009, 13:29
Default incorrect results with a tetrahedral mesh
  #1
New Member
 
Jonathan Masters
Join Date: Mar 2009
Posts: 4
Rep Power: 8
jcmasters is on a distinguished road
I am having some trouble getting the simple case of laminar flow over a sphere in a rectangular channel (3D) to give correct results on a tetrahedral mesh. This problem is only present on tet meshes, not hex. I've pasted an example of the results I get (in this case from potentialFoam, but all other apps give roughly the same result) plotted on a planar slice down the center:



Uniform flow should be entering from the left, top and bottom (as well as front and back not shown) are walls, and outflow has a pressureInletVelocity BC applied to it. The results in all cases seem to converge fine. It doesn't seem to be attributable to discretization schemes or linear solvers, I've tried them all.

I am running OpenFOAM-1.5 on a 32-bit dual core intel machine loaded with Ubuntu 8.2. I also reinstalled version 1.4.1 (which I had used previously with success) but I found the same thing. I have tried all the meshing apps at my disposal (Gambit, Workbench Meshing, ICEM, Salome Platform) and have found the same thing with each. Has anyone else found this? Could this simply be a post-processing problem?
jcmasters is offline   Reply With Quote

Old   May 20, 2009, 13:40
Default
  #2
Senior Member
 
Wolfgang Heydlauff
Join Date: Mar 2009
Location: Germany
Posts: 136
Rep Power: 8
wolle1982 is on a distinguished road
the same as ever:
- check Boundary conditions
- use pressure outlet
- create boundary layer
- don't use potential foam, use turbfoam. it's the most stable solver in my opinion
wolle1982 is offline   Reply With Quote

Old   May 20, 2009, 14:03
Default
  #3
New Member
 
Jonathan Masters
Join Date: Mar 2009
Posts: 4
Rep Power: 8
jcmasters is on a distinguished road
Thanks for your suggestions. I've tried several solvers (icoFoam, turbFoam, simpleFoam, ...) but the problem persists, BCs all look good (same simulation in CFX looks fine), outlet BC is pressure. I also had a more complex geometry with BL inflation and found the same thing. checkMesh gives me no errors. I've included non-ortho correctors to no avail. Using more diffusive discretization schemes (e.g. upwind) had no effect on the splotchiness of the results.
jcmasters is offline   Reply With Quote

Old   May 20, 2009, 16:14
Default
  #4
Senior Member
 
santos's Avatar
 
Jose Luis Santos
Join Date: Mar 2009
Location: Portugal
Posts: 213
Rep Power: 9
santos is on a distinguished road
Send a message via Skype™ to santos
Hi Jonathan,

I am running thousands (literally!) of simulations of creeping flow around a sphere, and did not find anything like the results you show. My setup is similar to yours, although I have symmetry on front and back.

Maybe you could double-check your BC's, as suggested by Wolfgang.

Regards,
Jose Santos
santos is offline   Reply With Quote

Old   May 20, 2009, 16:26
Default
  #5
New Member
 
Jonathan Masters
Join Date: Mar 2009
Posts: 4
Rep Power: 8
jcmasters is on a distinguished road
Hi Jose

I agree, it's definitely possible. I've been using OpenFOAM for several year now and have had no trouble. I'm wondering whether it has something to do with my OS, although I don't know why it would. It works great with hex grids (I get very stable solutions for simulations of fluids with complex rheology) and should for any grid type.

Are you're simulations on tet grids? What OS are you running on?

Thanks

Jonathan
jcmasters is offline   Reply With Quote

Old   May 20, 2009, 16:44
Default
  #6
Senior Member
 
santos's Avatar
 
Jose Luis Santos
Join Date: Mar 2009
Location: Portugal
Posts: 213
Rep Power: 9
santos is on a distinguished road
Send a message via Skype™ to santos
Hi,

I am working with meshes generated with snappyHexMesh, so I have a mix of hex, prisms and polyhedra. My OS is openSUSE 10.3 64 bits.

Regards,
Jose
santos 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
fluent add additional zones for the mesh file SSL FLUENT 2 January 26, 2008 12:55
basic of mesh refinement arya CFX 4 June 19, 2007 12:21
ICEM 10 mesh question DAK565656 CFX 6 May 8, 2007 12:16
Please Help...!need SIMPLE TETRAHEDRAL MESH GENERA necdet aslan FLUENT 0 January 2, 2006 13:11
Difference in results using different Mesh? Julie FLUENT 7 July 2, 2002 10:28


All times are GMT -4. The time now is 19:44.