t

Grid Convergence Index

• 1. Grid/mesh independence   GCI

  • 1.1. Richardson extrapolation
  • 1.2. grid refinement ratio
  • 1.3. Example of Grid convergence study
  • 1.4. calculation steps
  • 1.5. Example of Grid convergence for wing:
  • 1.6. References

keywords: Richardson's extrapolation, Grid convergence index a summary of Richardson's extrapolation is here

requirement: GCI < 5%

a summary of GCI from nasa web , local downloaded file is [here](\akmkemin\Backup\ANSYS\Geo_mesh_solver_post_processing\Convergence-Verification-validation\grid convergence) ( print version is in BEM file folder)

1.1 Richardson extrapolation

1.2 grid refinement ratio

• Hexa mesh –>> grid refinement ratio

  • double nodes along each coordinates (x, y,z)

• Tetra mesh –>> effective grid refinement ratio

Definitions:

\[ r_{ij} = h_i/h_j \]

• r: grid refinement ratio,
• h: grid spacing

1.3 Example of Grid convergence study

The example is from here,
The Fortran 90 program verify.f90 was written to carry out the calculations associated with a grid convergence study involving 3 or more grids
The program is compiled on a unix system through the commands:

f90 verify.f90 -o verify

It reads in an ASCII file (prD.do) through the standard input unit (5) that contains a list of pairs of grid size and value of the observed quantity f.

input data format:
1.0 0.97050 2.0 0.96854 4.0 0.96178
verify < prD.do > prD.out

It assumes the values from the finest grid are listed first. The output is then written to the standard output unit (6) prD.out.
The output from the of {\tt verify} for the results of Appendix A are:

#+BEGINEXAMPLE

— VERIFY: Performs verification calculations —

Number of data sets read = 3

Grid Size Quantity

1.000000 0.970500
2.000000 0.968540
4.000000 0.961780

Order of convergence using first three finest grid
and assuming constant grid refinement (Eqn. 5.10.6.1)
Order of Convergence, p = 1.78618479

Richardson Extrapolation: Use above order of convergence
and first and second finest grids (Eqn. 5.4.1)
Estimate to zero grid value, fexact = 0.971300304

Grid Convergence Index on fine grids. Uses p from above.
Factor of Safety = 1.25

Grid Refinement
Step Ratio, r GCI(%)
1 2 2.000000 0.103080
2 3 2.000000 0.356244

Checking for asymptotic range using Eqn. 5.10.5.2.
A ratio of 1.0 indicates asymptotic range.

Grid Range Ratio
12 23 0.997980

— End of VERIFY —

#+END _EXAMPLE

1.4 calculation steps

1. Complete at least 3 simulations (Coarse, medium, fine) with a constant refinement ratio, r, between them (in our example we use r=2)
2. Choose a parameter indicative of grid convergence. In most cases, this should be the parameter you are studying. ie if you are studying drag, you would use drag.
3. Calculate the order of convergence, p, using:

\begin{equation}
p=ln(\frac{f_3 - f_2}{f_2- f_1}) / \ln (r)
\end{equation}

where \( f_i \) is the solution at different meshes, f1 is fine grid, \( r \) is grid refinement ratio.

1. Perform a Richardson extrapolation to predict the value at h=0

\begin{equation}
f_e = f_1 + \frac{f_1 -f_2 }{r^p - 1}
\end{equation}

fe, exact numerical value ( continuum value at zero grid spacing)

1. Calculate grid convergence index (GCI) for the medium and fine refinement levels

\begin{equation}
GCI_{fine} = \frac{F_s \vert \epsilon \vert }{r^p - 1}
\end{equation}

where \( F_s \) is a safety factor. the recommended value is 3 for two grids comparisons and 1.25 for three or more grids comparisons.

1.5 Example of Grid convergence for wing:

C:\akmkemin\Backup\academic\Notes\ANSYS\mesh\Modules\wings\wing-mesh-comparing-pointwise.doc

1.6 References

Roache, P. J. Perspective: A Method for Uniform Reporting of Grid Refinement Studies, Journal of Fluids Engineering, Vol. 116, 1994; 405-413.

Roache, P. J. Quantification of Uncertainty in Computational Fluid Dynamics, in Annual Review of Fluid Mechanics

Roache, Patrick J. Verification and validation in computational science and engineering. Vol. 895. Albuquerque, NM: Hermosa, 1998.

Author: kaiming

Created: 2019-04-16 Tue 22:15

Emacs 24.5.1 (Org mode 8.2.10)

Validate