NACA0012 grid dependence study by Michael Alletto - Revision history

Michael Alletto at 08:53, 6 January 2022

2022-01-06T08:53:04Z

Jozsef Nagy at 05:52, 14 May 2020

2020-05-14T05:52:02Z

Jozsef Nagy at 05:36, 14 May 2020

2020-05-14T05:36:19Z

Jozsef Nagy at 05:35, 14 May 2020

2020-05-14T05:35:42Z

Jozsef Nagy at 05:35, 14 May 2020

2020-05-14T05:35:10Z

Jozsef Nagy at 05:34, 14 May 2020

2020-05-14T05:34:31Z

Jozsef Nagy: Created page with "* '''contributor''': Michael Alletto * '''affiliation''': WRD MS * '''contact''': click here for email address
2020-05-14T05:34:09Z

Created page with "* '''contributor''': Michael Alletto * '''affiliation''': WRD MS * '''contact''': <mail address='michael.alletto@gmx.de' description='author'>click here for email address</mai..."

New page
* '''contributor''': Michael Alletto
* '''affiliation''': WRD MS
* '''contact''': <mail address='michael.alletto@gmx.de' description='author'>click here for email address</mail>
* '''OpenFOAM versions''': v1912
* '''Published under''': CC BY-NC-ND license ([https://creativecommons.org/licenses creative commons licenses])

Go back to [https://wiki.openfoam.com/NACA0012_by_Michael_Alletto NAVA0012 tutorial page].

=NACA0012 grid dependence study=

==Results stationary airfoil==

As stated in the introduction, the best way to evaluate the performance of a given solver to treat a complex problem, is to have a series of test cases with increasing complexity. Each test case can be evaluated separately allowing to track down the strength and limitations of the solver. Furthermore the interpretation of the results is much easier in a simple case compared to a complex one. The gained knowledge of the simple test cases can then be applied to understand the more complex problem. In this spirit, the first evaluation performed is the one using a fixed grid.

Note that all python scrips which are used to generate the plots shown in this document are included in the provided tutorial.

==[https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/MichaelAlletto/gridConvergenceStudy%20(1).tar.gz Grid size variation]==

Files can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/MichaelAlletto/gridConvergenceStudy%20(1).tar.gz here].

This section shows the results of the grid independence study. The size and the labels of the grids evaluate are shown in the table below. The coarsest grid had 170 points in chord wise direction and 40 points in direction normal to the airfoil surface. The finest grid had 400 points in chord wise direction and 100 points in direction normal to the airfoil surface

{|
!align="center"|

!align="center"| grid0
!align="center"| grid1
!align="center"| grid2
!align="center"| grid3
|-
|align="center"| size
|align="center"| 170x40
|align="center"| 200x50
|align="center"| 299x75
|align="center"| 400x100
|-
|align="center"|
|align="center"|

|align="center"|

|align="center"|

|align="center"|

|}

The next figure shows the coarsest grid used (left) and a zoom of the grid around the airfoil.

[[File:Grid0_alettonaca.png|1000px|center| Grid used]]

The figure below shows the variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom left) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. Since we use local time stepping we have to ensure that the solution is sufficiently converged. For this purpose it is suitable to analyze the evolution of some integral quantity. Since the lift is the quantity we are most interested in, we plotted it over the time. It is evident that the lift coefficient obtained a steady state before the end of the simulation. The experiments of (labeled as '''exp''') and the simulation of (labeled as '''ref sim''') are included as reference. The reference simulation used a grid of approximately 40 000 cells where the wall was fully resolved (all first cell centers lied at y+ <math><</math> 1). They used the Spalart Allmaras model with the modifications introduced by together with a density based solver. It is evident that already for grid1 there is no noticeable difference to the finer grid regarding the pressure coefficient. The friction coefficient still changes with the grid size. Noticeable is that the present set up with equilibrium wall functions predicts a small shock induces recirculation bubble which reattaches at approximately 0.7 x/c. The extend of the recirculation bubble is similar to the reference simulation which fully resolved the wall. Y+ is in most regions of the flow between 40 and 300 for all grids where the assumption made to derive the wall functions are valid.

[[File:gridSTudy_alettonaca.png|450px|center|variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. The experiments of [3] (labeled as exp) and the simulation of [2] (labeled as ref sim) are included as reference.]

Since there is almost no difference in the pressure coefficient between grid1 and the finer grids we will continue our evaluation with this grid. This choice is motivated as follows: 1) We are interest in evaluating the performance of the solver to predict the loads on periodically pitching airfoils. Since the loads are predetermined predominantly by the pressure distribution around the airfoil the accurate calculation of the skin friction is not required. 2) We will consider only small angles of attack where the flow remains attached or only a small recirculation bubble appears. For this kind of flows wall functions and a grid which is not too fine is evidently good enough. 3) For the purpose of a tutorial we want to keep the computational time as limited a possible.

At this point we want to spend a few words about what is a grid which is good enough. In this tutorial we do only visual comparison to keep things as simple as possible. In practice one will do quantitative comparison on the quantities of interests, e.g. the lift computed and calculate the error compared to experiments. Once one has obtained an accurate estimate of the error made by the computation, it can be decided if the errors are small enough or too big for a given application.

==References==

[1] Jack R Edwards and Suresh Chandra. Comparison of eddy viscositytransport turbulence models for three-dimensional, shock-separated flowfields. AIAA journal, 34(4):756–763, 1996.

[2] Michael Iovnovich and Daniella E Raveh. Reynolds-averaged navierstokes study of the shock-buffet instability mechanism. AIAA journal, 50(4):880–890, 2012.

[3] John B McDevitt and Arthur F Okuno. Static and dynamic pressure measurements on a naca 0012 airfoil in the ames high reynolds number facility. 1985.

@@ Line 52: / Line 52: @@
 The figure below shows the variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom left) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. Since we use local time stepping we have to ensure that the solution is sufficiently converged. For this purpose it is suitable to analyze the evolution of some integral quantity. Since the lift is the quantity we are most interested in, we plotted it over the time. It is evident that the lift coefficient obtained a steady state before the end of the simulation. The experiments of  (labeled as '''exp''') and the simulation of  (labeled as '''ref sim''') are included as reference. The reference simulation used a grid of approximately 40 000 cells where the wall was fully resolved (all first cell centers lied at y+ < 1). They used the Spalart Allmaras model with the modifications introduced by  together with a density based solver. It is evident that already for grid1 there is no noticeable difference to the finer grid regarding the pressure coefficient. The friction coefficient still changes with the grid size. Noticeable is that the present set up with equilibrium wall functions predicts a small shock induces recirculation bubble which reattaches at approximately 0.7 x/c. The extend of the recirculation bubble is similar to the reference simulation which fully resolved the wall. Y+ is in most regions of the flow between 40 and 300 for all grids where the assumption made to derive the wall functions are valid.
-[[File:gridSTudy_alettonaca.png|450px|center|variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. The experiments of [3] (labeled as exp) and the simulation of [2] (labeled as ref sim) are included as reference.]]
+[[File:gridSTudy_alettonaca.png|750px|center|variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. The experiments of [3] (labeled as exp) and the simulation of [2] (labeled as ref sim) are included as reference.]]
 Since there is almost no difference in the pressure coefficient between grid1 and the finer grids we will continue our evaluation with this grid. This choice is motivated as follows: 1) We are interest in evaluating the performance of the solver to predict the loads on periodically pitching airfoils. Since the loads are predetermined predominantly by the pressure distribution around the airfoil the accurate calculation of the skin friction is not required. 2) We will consider only small angles of attack where the flow remains attached or only a small recirculation bubble appears. For this kind of flows wall functions and a grid which is not too fine is evidently good enough. 3) For the purpose of a tutorial we want to keep the computational time as limited a possible.

@@ Line 52: / Line 52: @@
 The figure below shows the variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom left) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. Since we use local time stepping we have to ensure that the solution is sufficiently converged. For this purpose it is suitable to analyze the evolution of some integral quantity. Since the lift is the quantity we are most interested in, we plotted it over the time. It is evident that the lift coefficient obtained a steady state before the end of the simulation. The experiments of  (labeled as '''exp''') and the simulation of  (labeled as '''ref sim''') are included as reference. The reference simulation used a grid of approximately 40 000 cells where the wall was fully resolved (all first cell centers lied at y+ < 1). They used the Spalart Allmaras model with the modifications introduced by  together with a density based solver. It is evident that already for grid1 there is no noticeable difference to the finer grid regarding the pressure coefficient. The friction coefficient still changes with the grid size. Noticeable is that the present set up with equilibrium wall functions predicts a small shock induces recirculation bubble which reattaches at approximately 0.7 x/c. The extend of the recirculation bubble is similar to the reference simulation which fully resolved the wall. Y+ is in most regions of the flow between 40 and 300 for all grids where the assumption made to derive the wall functions are valid.
-[[File:gridSTudy_alettonaca.png|450px|center|variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. The experiments of [3] (labeled as exp) and the simulation of [2] (labeled as ref sim) are included as reference.]
+[[File:gridSTudy_alettonaca.png|450px|center|variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. The experiments of [3] (labeled as exp) and the simulation of [2] (labeled as ref sim) are included as reference.]]
 Since there is almost no difference in the pressure coefficient between grid1 and the finer grids we will continue our evaluation with this grid. This choice is motivated as follows: 1) We are interest in evaluating the performance of the solver to predict the loads on periodically pitching airfoils. Since the loads are predetermined predominantly by the pressure distribution around the airfoil the accurate calculation of the skin friction is not required. 2) We will consider only small angles of attack where the flow remains attached or only a small recirculation bubble appears. For this kind of flows wall functions and a grid which is not too fine is evidently good enough. 3) For the purpose of a tutorial we want to keep the computational time as limited a possible.

@@ Line 2: / Line 2: @@
 * '''affiliation''': WRD MS
 * '''contact''': <mail address='michael.alletto@gmx.de' description='author'>click here for email address</mail>
-* '''OpenFOAM versions''': v1912
+* '''OpenFOAM versions''': v2112
 * '''Published under''': CC BY-NC-ND license ([https://creativecommons.org/licenses creative commons licenses])
@@ Line 9: / Line 9: @@
 =NACA0012 grid dependence study=
-Files can be downloaded [https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/MichaelAlletto/gridConvergenceStudy%20(1).tar.gz here].
+Files can be downloaded [https://gitlab.com/mAlletto/openfoamtutorials/-/tree/master/Naca0012transonic here].
 ==Results stationary airfoil==

@@ Line 5: / Line 5: @@
 * '''Published under''': CC BY-NC-ND license ([https://creativecommons.org/licenses creative commons licenses])
-Go back to [https://wiki.openfoam.com/NACA0012_by_Michael_Alletto NAVA0012 tutorial page].
+Go back to [https://wiki.openfoam.com/NACA0012_by_Michael_Alletto NACA0012 tutorial page].
 =NACA0012 grid dependence study=

@@ Line 50: / Line 50: @@
 [[File:Grid0_alettonaca.png|1000px|center| Grid used]]
-The figure below shows the variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom left) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. Since we use local time stepping we have to ensure that the solution is sufficiently converged. For this purpose it is suitable to analyze the evolution of some integral quantity. Since the lift is the quantity we are most interested in, we plotted it over the time. It is evident that the lift coefficient obtained a steady state before the end of the simulation. The experiments of  (labeled as '''exp''') and the simulation of  (labeled as '''ref sim''') are included as reference. The reference simulation used a grid of approximately 40 000 cells where the wall was fully resolved (all first cell centers lied at y+ <math><</math> 1). They used the Spalart Allmaras model with the modifications introduced by  together with a density based solver. It is evident that already for grid1 there is no noticeable difference to the finer grid regarding the pressure coefficient. The friction coefficient still changes with the grid size. Noticeable is that the present set up with equilibrium wall functions predicts a small shock induces recirculation bubble which reattaches at approximately 0.7 x/c. The extend of the recirculation bubble is similar to the reference simulation which fully resolved the wall. Y+ is in most regions of the flow between 40 and 300 for all grids where the assumption made to derive the wall functions are valid.
+The figure below shows the variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom left) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. Since we use local time stepping we have to ensure that the solution is sufficiently converged. For this purpose it is suitable to analyze the evolution of some integral quantity. Since the lift is the quantity we are most interested in, we plotted it over the time. It is evident that the lift coefficient obtained a steady state before the end of the simulation. The experiments of  (labeled as '''exp''') and the simulation of  (labeled as '''ref sim''') are included as reference. The reference simulation used a grid of approximately 40 000 cells where the wall was fully resolved (all first cell centers lied at y+ < 1). They used the Spalart Allmaras model with the modifications introduced by  together with a density based solver. It is evident that already for grid1 there is no noticeable difference to the finer grid regarding the pressure coefficient. The friction coefficient still changes with the grid size. Noticeable is that the present set up with equilibrium wall functions predicts a small shock induces recirculation bubble which reattaches at approximately 0.7 x/c. The extend of the recirculation bubble is similar to the reference simulation which fully resolved the wall. Y+ is in most regions of the flow between 40 and 300 for all grids where the assumption made to derive the wall functions are valid.
 [[File:gridSTudy_alettonaca.png|450px|center|variation of the pressure coefficient (left), the skin friction coefficient (right) and y+ (figure at the bottom) with the grid size. At the bottom right the evolution of the lift coefficient with time is shown. The experiments of [3] (labeled as exp) and the simulation of [2] (labeled as ref sim) are included as reference.]

@@ Line 8: / Line 8: @@
 =NACA0012 grid dependence study=
 ==Results stationary airfoil==
@@ Line 16: / Line 18: @@
 ==[https://github.com/jnmlujnmlu/OpenFOAMTeaching/blob/master/MichaelAlletto/gridConvergenceStudy%20(1).tar.gz Grid size variation]==
 This section shows the results of the grid independence study. The size and the labels of the grids evaluate are shown in the table below. The coarsest grid had 170 points in chord wise direction and 40 points in direction normal to the airfoil surface. The finest grid had 400 points in chord wise direction and 100 points in direction normal to the airfoil surface