
[Sponsors] 
October 28, 2013, 07:58 
The old Hex vs Tet question

#1 
Senior Member
Join Date: May 2012
Posts: 137
Rep Power: 5 
Hi,
For nonspecific flows: 1. Hex is always better than Tet when it comes to accuracy. 2. No, Tet is as good as Hex if both solutions are mesh independent, although Tet require more cells to accomplish that. Is it possible that the question is also related to cell skewness and that many researchers have just compared apples and oranges? Cheers! 

October 28, 2013, 09:25 

#2 
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,134
Rep Power: 19 
I just finished putting together some "best practice" guidelines on this and some related topics.
My opinion based on my personal experience with commercial CFD software is: You can achieve good results with tet meshes for most engineering applications if you obey the basic rules of mesh generation (overall mesh quality, volume jump, angles, boundary layer resolution...) But you have to know the drawbacks: Solutions with first order upwind schemes on tetrahedral meshes will always suffer from high numerical diffusion. The computing times with tet meshes will be higher for the same accuracy and the overall convergence behavior is worse than with hex meshes. And when it comes to LES with commercial cfd software, I would always choose hex meshes. 

October 28, 2013, 09:57 

#3  
Senior Member
Join Date: May 2012
Posts: 137
Rep Power: 5 
Quote:
Quote:
1. How is the second statement related to the first statement? 2. Do you mean that we need more cells in order to reduce numerical diffusion? 3. Numerical diffusion is affected by mesh refinement so I would guess that it is captured when doing a mesh sensitivity analysis? 4. It seems that you are leaning towards the second type of answer in my original post, right? 5. Assume that the flow is at a 45 degree angle with regards to the Hex mesh, now how will this affect the assessment of severe numerical diffusion for the Tet mesh compared to the Hex mesh? Cheers! 

October 28, 2013, 10:43 

#4  
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,134
Rep Power: 19 
Quote:
Simply refining the mesh reduces the numerical diffusion, but only with a slope of 1. If your main focus is accuracy then you will need unreasonably fine meshes just to get rid of this error source. Quote:
Nevertheless, one of the goals when creating hex meshes is to "streamline" the grid to prevent flows with a 45 degree angle with respect to the cells. 

October 28, 2013, 11:21 

#5 
Senior Member
Join Date: May 2012
Posts: 137
Rep Power: 5 
Ok so here are some 2d simulations with a commercial software.
The setup is a square domain with inflow at west and south boundary at a 45 degree angle. A passive scalar is introduced with a value of 1 at the west boundary and a value of 0 at the south boundary. Plots are from southeast corner to northwest corner. Check attachments for results. It is clear that the quads have more numerical diffusion in this case. Convergence took much longer for the triangles case. Residuals in both cases were lowered to 1e8. 

October 28, 2013, 13:28 

#6 
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,134
Rep Power: 19 
Rotating the flow direction by 90 degrees, you will get a similar diffusion error on the tet mesh.


October 28, 2013, 14:14 

#7 
Senior Member
Join Date: May 2012
Posts: 137
Rep Power: 5 
Yes. So does a tri mesh have three preferred directions (low numerical diffusion) as opposed to two directions for quad elements? If yes, then it seems to me that polyhedral elements would be superior, given that they may have even more preferred directions.


October 31, 2013, 07:21 

#8 
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 1,134
Rep Power: 19 
I would not draw the same conclusion.
In your test case with the tet mesh, the numerical diffusion is only zero because the line of separation of the passive scalar has the correct direction AND the correct position. If you could move the line of separation by half a cell size, you would still get diffusion errors. Ergo tet elements have no preferred direction at all. 

October 31, 2013, 07:36 

#9  
Senior Member
Join Date: May 2012
Posts: 137
Rep Power: 5 
Quote:
2. Textbooks talk about numerical diffusion being worse when the flow is oblique to the direction of the grid. Should they also add "the correct position" to this statement? (perhaps I don't understand the concept of "grid direction"). Cheers! 

October 31, 2013, 07:48 

#10  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,856
Rep Power: 25 
Quote:
numerical diffusion is related to the local truncation error appearing when the convective terms are discretized in nonsymmetric way. Hex and Tet cells affects the way in which the discrete operator can be computed on the computational domain. Hex cells are generally associated to structured grids that generate "line directions", Tet cells generates nonstructured grid and you cannot uniquely define a grid direction. I suggest to do some simple 2D exercise of the solution of the linear equation df/dt+udf/dx+vdf/y=0, you can prescribe u and v. Try to solve on both triangular grid and hex grid. Try to use upwind and central discretization 

October 31, 2013, 08:01 

#11  
Senior Member
Join Date: May 2012
Posts: 137
Rep Power: 5 
Quote:


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
[ICEM] hybrid mesh/Aligning tet with hex mesh  shridharsn  ANSYS Meshing & Geometry  0  October 17, 2012 11:29 
Blockmesh error  2D scramjet  please help  ishaninair  OpenFOAM Native Meshers: blockMesh  7  March 18, 2011 01:14 
HEX or TET  Which One is better choice?  MAB  CFX  16  July 25, 2008 02:40 
combination of hex and tet mesh  Raju  CFX  1  July 27, 2006 09:10 
ANSYS ICEM CFD: Hex and Tet Element Specification  Evan  CFX  1  May 24, 2006 17:09 