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

Why Differences between tetrahedra polyhedra and prism mesh pipe flow computation with sonicLiquidFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 15, 2008, 13:43
Default I want to simulate flow in a p
  #1
New Member
 
Martin/Run DU
Join Date: Mar 2009
Location: Hong Kong
Posts: 25
Rep Power: 7
chnrdu is on a distinguished road
Send a message via Skype™ to chnrdu
I want to simulate flow in a pipe system which should reflect the pressure pulsation, so I use the solver sonicLiquidFoam.

I have used three types of mesh to simulate it in order to choose a best mesh to be used further: tetrahedra, prism, polyhedra. But the polyhedra, which I am very interested in, cause abnormal pressure:

meshes: 11484
left boundary: pulsate pressure
right boundary: wall
time: 0.003125 s
pressure range: -5.96216e+6, 1.39506e+7


meshes: 44500
boundaries as above, time is the same to above
pressure range: -78, 120


meshes: 51225
boundaries as above, time is the same to above
pressure range: -66, 106

I learned, from a commercial CFD code company, the polyhedral mesh have more advantages: more accuracy, less time cost, fewer meshes, generating mesh automatically. However, the case I simulated is far from my impression. In deed, the time cost is decreased, but the error is too large.

I'd like to ask several questions:
1. why the result computed with polyhedral mesh are far away from the tetra and prism mesh?
2. How to improve the result computed with polyhedral mesh?
__________________
rdu
------------------
Martin/Run Du
chnrdu is offline   Reply With Quote

Old   June 16, 2008, 02:41
Default Hello! The first question i
  #2
New Member
 
Gregor Veble
Join Date: Mar 2009
Posts: 11
Rep Power: 7
gregorv is on a distinguished road
Hello!

The first question is of course the quality of the polyhedral mesh. Did you run checkMesh to see if there are any problems with it? It seems to me that the simulation might be diverging, which can be caused even by only a few bad cells.

Polyhedral meshes may indeed be advantageous in some respects, but getting the polyhedra of adequate shapes is a bit harder than with e.g. tetrahedra.
gregorv is offline   Reply With Quote

Old   June 16, 2008, 03:37
Default Hi, Gregor. Thanks for your
  #3
New Member
 
Martin/Run DU
Join Date: Mar 2009
Location: Hong Kong
Posts: 25
Rep Power: 7
chnrdu is on a distinguished road
Send a message via Skype™ to chnrdu
Hi, Gregor.

Thanks for your reply.

I have checked the mesh. the checkMesh result is:

Create polyMesh for time = constant

Time = constant

Mesh stats
points: 61690
edges: 123380
faces: 73175
internal faces: 67911
cells: 11484
boundary patches: 3
point zones: 0
face zones: 0
cell zones: 0

Number of cells of each type:
hexahedra: 31
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 11453

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface
F1 74 146 ok (not multiply connected)
F2 88 174 ok (not multiply connected)
Wall 5102 10204 ok (not multiply connected)

Checking geometry...
Domain bounding box: (0 -0.0159475 -0.0159343) (1 0.0159326 0.0159341)
Boundary openness (5.1346043e-18 -7.9059664e-18 5.7280395e-16) OK.
Max cell openness = 2.2111137e-16 OK.
Max aspect ratio = 7.2639034 OK.
Minimum face area = 5.2550992e-07. Maximum face area = 4.7376251e-05. Face area magnitudes OK.
Min volume = 2.0585726e-09. Max volume = 6.2445933e-07. Total volume = 0.00078484389. Cell volumes OK.
Mesh non-orthogonality Max: 43.57726 average: 12.389229
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.3541358 OK.
Min/max edge length = 0.00058204331 0.0040980817 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : average = 0.99289412 min = 0.85956038
All face flatness OK.

Mesh OK.

Maybe there are some respects to improve. Please point out.

Thank you!
__________________
rdu
------------------
Martin/Run Du
chnrdu is offline   Reply With Quote

Old   June 16, 2008, 16:41
Default Hello again, as far as I ca
  #4
New Member
 
Gregor Veble
Join Date: Mar 2009
Posts: 11
Rep Power: 7
gregorv is on a distinguished road
Hello again,

as far as I can see the mesh itself is perfectly OK.

Are the sizes of the elements for the tetrahedral, prismatic and polyhedral meshes comparable? I still think your solution is simply not converging for whatever reason. It might be something as simple as having too large a timestep, or something more intricate such as various discretisation schemes not being stable enough in the polyhedral formulation.

Unfortunately I have no experience with the actual solver you are using so I can't say much more of use. It is, however, most likely a stability issue, not an accuracy one.
gregorv is offline   Reply With Quote

Old   June 17, 2008, 07:21
Default try running checkMesh with op
  #5
Senior Member
 
mkraposhin's Avatar
 
Matvej Kraposhin
Join Date: Mar 2009
Location: Moscow, Russian Federation
Posts: 172
Rep Power: 7
mkraposhin is on a distinguished road
try running
checkMesh with option -cellDeterminant
mkraposhin 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
Problems on Tetrahedra to Polyhedra using polyDualMesh on concave surface object chnrdu OpenFOAM Pre-Processing 2 July 18, 2009 13:58
prism mesh in 2D icem-noob CFX 4 May 2, 2008 12:11
How to generate refined 3D tetrahedra etc mesh for free kar OpenFOAM Meshing Format & General Technical 8 January 29, 2008 05:21
Error in prism mesh Joao Lourenco CFX 3 October 30, 2007 07:41
Differences between solid/ fluid mesh Vidya Raja FLUENT 0 March 14, 2006 18:49


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