CFD Online Discussion Forums

CFD Online Discussion Forums (
-   CFX (
-   -   Grid refinement ratio for unstructured mesh (

oj.bulmer February 28, 2013 12:17

Grid refinement ratio for unstructured mesh

I am trying to streamline the mesh independence practices using Grid Convergence Index and Richardson Extrapolation as they are considered the best in practical CFD. (Ref: Roache, P. J. "Perspective: a method for uniform reporting of grid refinement studies." TRANSACTIONS-AMERICAN SOCIETY OF MECHANICAL ENGINEERS JOURNAL OF FLUIDS ENGINEERING 116 (1994): 405-405.)

Now, element sizes for the consecutive meshes should be known to get the grid refinement ratio used in the calculations in this procedure. But since we use unstructured mesh along with the extruded mesh; I am not sure if the value of element local element size (in the important regions of high gradients) is representative of the size for that particular mesh.

Roache suggests the following formula for grid refinement ratio r for unstructured meshes, using their element-counts N_1 and N_2:

r = \left(\frac{N_1}{N_2}\right)^D

Now, he defines D as Dimensionality. What does it mean and how to calculate it from the fluid domain?? :confused:


ghorrocks February 28, 2013 17:47

You would have to read the paper for the definition of D. But I suspect it will be 2 for 2D models and 3 for 3D models.

oj.bulmer March 1, 2013 05:22

Thanks Glenn.

The only definition Roache gives of D is Dimensionality! But a little more research online after posting my query here landed me to some references that reciprocate the similar definition that you mentioned. Dimensionality is the no. of principle directions in which there is a significant variation of the results. Essentially, a 2D case can have single dimensionality if there are gradients only in single directions.

If I may ask, do you personally recommend this method for mesh independence when it comes to unstructured meshes? Or you rely on more general methods of assessing the parameters of interest from results of successive meshes and choose the mesh when you see no significant variation in them? The only problem in this approach is there seems to be no agreed quantification of word "significant".

In this paper, Richardson's extrapolation takes into account the results of meshes that are successively refined and then extrapolates this to predict the value of result at zero grid spacing (i.e. continuum). At the same time, the Grid Convergence Index helps understand if the results generated in the current successive meshes are sufficiently close to the asymptotic range, where any refinement in the mesh won't result in significant change in results. This way we not only ensure the mesh independence but also predict the numerical solution at continuum. Obviously, there will be a difference between actual result and numerical result at continuum, but then, CFD is about getting the best results within its caveats.


ghorrocks March 3, 2013 06:11

The definition of significant is what you define as significant. If you are doing initial design calcs then 20% error might be fine - you just want a general ida of what is going on. But for final design work you might want 5% error. Or if you are doing detailed turbomachinery work then 1% error is suiable.

oj.bulmer March 4, 2013 05:28

Thanks. When you say 5% change in solution by refinement, what is the refinement ratio for the mesh size?

I go a bit detailed. If f_1 and f_2 are the solutions of two successive meshes ( f_1 smaller), then is the error you mentioned, say epsilon_{12}, must be following:

\epsilon_{12} = \frac{f_2-f_1}{f_1} * 100 \%

The fundamental problem with this, accordign to Roahce, is that this error doesn't take into account the refinement ratio r and the global order of accuracy p, which is typically smaller than the local order of accuracy of stencil used for discretization, due to grid stretching and quality, non-linearities and assumptions in turbulence modelling etc.

The Roache then defines Grid Convergence Index which takes into account these uncertainties as:

GCI_{12} = \frac{F_s |\epsilon_{12}|}{r^p - 1} * 100 \%

F_s being the factor of safety, 3 for two meshes, 1.25 for three meshes used. In case of three meshes, the GCI for second and third mesh can be given similarly as:

GCI_{23} = \frac{F_s |\epsilon_{23}|}{r^p - 1} * 100 \%

Finally, with the relation GCI_{23} = r^p GCI_{12}, if the both sides are close to each other, it can be concluded that the meshes are in asymptotic range! The refinement ratio r is easy to calculate from the formula I mentioned in earlier post, and Roache suggests a minimum value of 1.1. So method is not too much demanding either.



ghorrocks March 4, 2013 05:36

I did not make myself clear enough in my previous post. I was commenting on your quote "The only problem in this approach is there seems to be no agreed quantification of word "significant".

For some people or cases the detailed analysis of Roache is not possible or practical so they fall back to the simple sensitivity analysis approach of two meshes with significnatly different meshes, and keep refining until the difference between two successive refinements in the parameters of interest (drag, lift, pressure drop, whatever is important) is less than an accuracy tolerance you are happy with.

Roache's approach does not require this as it checks you are in the convergent region and if so extrapolates to zero mesh. This is obviously a much better approach.

So my 5% was referring to changes in parameters of interest, not of changes in mesh size.

oj.bulmer March 4, 2013 06:45

Thanks, I understood that 5% was for the change in parameters of interest. I was asking about what should be percentage change in mesh sizes (% change in no. of elements for unstructured mesh) that you suggest for the two successive meshes in consideration :)

ghorrocks March 4, 2013 17:49

Oh - Roache gives guidelines on this as well.

But I generally double or halve the length dimension of the element. So a single hex element would become 8 for a 3D model, then 64. Often you are limited in how far you can take this before you run out of capacity, so one model is just "as big as I can fit in".

All times are GMT -4. The time now is 11:33.