Simulation of Turbulent Flow around an Airfoil

ENGR:2510 Mechanics of Fluids and Transfer Processes

CFD Lab 2
(ANSYS 17.1; Last Updated: Nov. 7, 2016)

By Timur Dogan, Michael Conger, Andrew Opyd, Dong-Hwan Kim

Maysam Mousaviraad, Tao Xing and Fred Stern

IIHR-Hydroscience & Engineering

The University of Iowa
C. Maxwell Stanley Hydraulics Laboratory
Iowa City, IA 52242-1585

1. Purpose
The Purpose of CFD Lab 2 is to conduct parametric studies for turbulent flow around Clark-Y airfoil
following the “CFD process” by an interactive step-by-step approach. Students will have “hands-on”
experiences using ANSYS to investigate the effect of angle of attack and effect of different turbulence
models on the simulations results. These effects will be studied by comparing simulation results with EFD
data. Students will analyze the differences and possible numerical errors, and present results in Lab report.

Geometry Physics Mesh Solution Results

Airfoil (ANSYS General (ANSYS Structured Solution Plots (ANSYS

Design Modeler) Fluent - Setup) (ANSYS Mesh) Methods Fluent- Results)
(ANSYS Fluent -
O-Domain Model (ANSYS
Solution)
(ANSYS Design Fluent - Setup) Graphics and
Modeler) Animations
Boundary Solution Controls (ANSYS Fluent-
Conditions (ANSYS Fluent - Results)
(ANSYS Fluent - Solution)
Setup)
Monitors
Reference Values (ANSYS Fluent -
(ANSYS Fluent - Solution)
Setup)
Turbulent Solution
Solution Initialization
Initialization (ANSYS Fluent -
(ANSYS Fluent - Solution)
Solution)
Run Calculation
(ANSYS Fluent -
Solution)

Flow chart for “CFD Process” for airfoil flow

1
2. Simulation Design
In EFD Lab 3, you have conducted experimental study for turbulent airfoil flow around a Clark-Y airfoil
(Re300,000). The data you have measured were used for CFD PreLab 2. In CFD Lab 2, simulation will be
conducted under the same conditions of EFD Lab 3, except angle of attack and turbulent models that will be
changed in this lab.

The problem to be solved is turbulent flow around the Clark-Y airfoil with angle of attack ()

Table 1 - Geometry dimensions

Parameters Symbol Unit Value
Chord Length C m 0.3048
Domain radius Rc m 12
Angle of attack α m 16

Figure 1 – Geometry

2
3. Opening ANSYS Workbench
3.1. Start > All Programs > ANSYS 17.1 > Workbench 17.1

3.2. From the ANSYS Workbench home screen (Project Schematic), drag and drop a Geometry, Mesh, and
two Fluent component from the Component Systems drop down menu onto the Project Schematic.
Project Schematic should resemble the schematic below. Rename the components as per below.

3.3. Create a Folder on the H: Drive called CFD Lab 2.

3.4. Save the project file by clicking File > Save As…
3.5. Save the project onto the H: Drive in the folder you just created and name it CFD Lab 2 Turbulent Flow.

3
4. Geometry
4.1. Right click Geometry then select Import Geometry > Browse… Select airfoil.igs and click OK.

4.2. Right click Geometry and select Edit Geometry in DesignModeler…

4.3. Click Generate.

4
4.4. Select Sketching > Constraints > Auto Constraints. Enable the auto constraints option to pick the
exact point as below.

5
4.5. Create > Body Transformation > Rotate. Select the airfoil and click Apply. Click the yellow box
labeled Axis Selection then click the XYPlane in the Tree Outline, then click Apply. Change the
Angle to 16° and click Generate.

1 2

3 4

6
4.6. Select XYPlane and click New Sketch button.

4.7. Select the sketch you created and click sketching button.

7
4.8. Sketching Toolbox > Draw > Circle. Click on the xy-plane origin and click behind the airfoil. (Click z-
arrow at right bottom to set the view as perpendicular to xy-plane)

4.9. Sketching Toolboxes > Dimensions > General. Click on the circle and change the diameter to 12m.

8
4.10. Concept > Surface From Sketches. Select Sketch1 from the Tree Outline, click Apply, then click
the Generate button.

4.11. Create > Boolean. Change operation to Subtract then select the circle for Target Bodies and click
Apply then select airfoil (click on the airfoil until only the airfoil’s color turns yellow) for Tool Bodies
and click Apply then click Generate. This will subtract the airfoil surface from the circle.

9
4.12. Concept > Split Edges. Select the perimeter of the circle and click Apply. Select Generate. This
should split the circle into two semicircles. You can see the semicircles by selecting the perimeter above
and below the x-axis.

4.13. Repeat the process from 4.12 on the two semicircles. This should yield four circular quadrants.

4.14. Repeat this process for the arc in quadrant IV (lower right). Change the Fraction to 0.822222. This
splits the arc into a 16° and a 74° arc.

10
4.15. Concept > Lines From Points. Draw a line from the point on the circle to the point on the airfoil
making sure to hold Ctrl while doing so. (Note: The point on the front part of the airfoil is not exactly on
the origin, you need to make sure to select the point on the airfoil and not the origin. Zoom in and find the
point just below the origin and select that point. The images below show the locations of the points circled.
When selecting points to generate lines with, be sure to select the point on circle and then the point on the
airfoil as to avoid complications when sizing mesh).

fluid domain

airfoil
origin

Line should be connected to this point

11
4.16. Once you select both points click Apply, then click Generate.

4.17. Repeat this process creating lines from the edge of the circle to the airfoil starting from the circle and
ending at the airfoil (A to A’, B to B’ and C to C’). The images below show the locations of the points on
the airfoil and the points on the circle.
A

A’

B’
C’
C

4.18. Tools > Projections. Select the four lines you created for Edges and select the circle for Target then
click Generate. This will split your geometry into four sections

12
4.19. Tools > Merge. Change the Merge Type to Edges. Select the 16° arc and the arc in quadrant I and
select Apply. Click Generate. This merges the lines into one line which can be sized for meshing easier.

4.20. File > Save project and exit.

5. Mesh Generation
5.1. From the Project Schematic right click on Mesh and select Edit… from the dropdown menu.

13
5.2. Right click Mesh then Insert > Face Meshing. Select the four surfaces while then click Apply.

14
5.3. Right click Mesh and Insert > Sizing. Select two edges as per below and change the parameters as per
below. You might need to change the cursor to “Edge selector” at this moment.

15
5.4. Right click Mesh and Insert > Sizing. Select two edges as per below and change the parameters as per
below.

5.5. Right click Mesh and Insert > Sizing. Select all for lines leading from the circle to the airfoil surface, and
click Apply. Change parameters as per below. Note: If you did not create the lines starting from the outer
circle and ending on the airfoil surface, you may have issues with biasing. If this is your case, size the
lines individually making sure that the sizing is finest at the surface of the airfoil.

16
5.6. Right click Mesh and Insert > Sizing. Select the surface at the top of leading edge of the airfoil and
click Apply. Change Parameters as per below.

5.7. Right click Mesh and Insert > Sizing. Select the surface at the bottom of leading edge of the airfoil and
click Apply. Change Parameters as per below.

5.8. Right click Mesh and Insert > Sizing. Select the surface at the top of trailing edge of the airfoil and
click Apply. Change Parameters as per below.

17
5.9. Right click Mesh and Insert > Sizing. Select the surface at the bottom of trailing edge of the airfoil and
click Apply. Change Parameters as per below.

5.10. Click on Mesh under the Outline, under the Details of “Mesh”, change the Physics Preference from
Mechanical to CFD. This changes the grid solver to a fluids style solver rather than a FEA style solver.

5.11. Click Generate Mesh. Click on the Mesh button under the Outline and make sure it resembles the
mesh below.

18
5.12. Hold Ctrl and select the two left most semicircle arcs, right click on them and select Create Named
Selection, name the selection inlet. Use the edge select button from the toolbar.

5.13. Do the same for the two right most semicircle arcs and name them outlet.

19
5.14. Select the four arcs that make the airfoil, right click and Create Named Selection and name it airfoil.

5.15. Now use the face button to select the four semicircle quadrants and Create Named Selection and
name them fluid.

5.16. File > Save Project. Exit the window.

5.17. Right click Mesh and select Update from the dropdown menu.

20
6. Setup
6.1. From the Project Schematic right click on Setup and select Edit…

6.2. Select Double Precision and click OK.

21
6.3. Solution Setup > Models > Viscous –Laminar > Edit… Change the parameters as per below and click
OK. (For the k-ω case, you will select k-omega (2 eqn)).

k-ε model

k-ω model

22
6.4. Solution Setup > Materials > air > Create/Edit… Change parameters as per experimental data and
click Change/Create.

The values in the figure above are for 24° temperature.

NOTE: viscosity used in ANSYS is the dynamic viscosity ( kg m  s ), NOT kinematic viscosity ( m 2 s )

23
6.5. Solution Setup > Boundary Conditions > inlet > Edit… Change parameters as per below.

k-ε model

k-ω model

24
6.6. Solution Setup > Boundary Conditions > Outlet > Edit… Change the parameters as per below and
click OK. (Use outlet B.C. for both k-ε and k-ω models)

6.7. Solution Setup > Reference Values. Change parameters as per below. The Velocity, Temperature,
Density, and Viscosity should be entered from EFD data.

25
6.8. Solution > Solution Methods. Change parameters as per below for both k-ε and k-ω models.

6.9. Solution > Solution Controls. Change Parameters as per below. (If you have problems with the
solution converging, you can decrease the Under –Relaxation Factors.)

Make sure you

scrolled down

26
6.10. Solution > Monitors > Residuals –Print, Plot > Edit…. Change the convergence limit to 1e-05
for all five equations.

6.11. Solution > Solution Initialization. Change the parameters as per below and click Initialize.

k-ε model k-ω model

27
6.12. Solution > Run Calculation. Change the Number of Iterations to 10000 and click Calculate.

Iteration history should look similar to the one below.

6.13. File > Save Project.

28
7. Results
7.1. Plotting and Saving Residual History

Solution > Monitors > Residuals –Print, Plot > Edit… > Plot > Ok.

File > Save Picture… > Save… Make sure the parameters are as per below and click Save… Name the file
“CFD Lab 2 Residual History” change the file directory to the CFD Lab 2 file you created on the H: drive and
click OK.

29
7.2. Plotting Pressure Coefficient Distribution with CFD and EFD Data

Results > Plots > XY Plot > Set Up… > Load File… Select “Pressure-coef-attack16.xy”. Change the
parameters as per below and click Plot. Save the picture the same way as you did for Residual History but in
this case, name it “CFD Lab 2 Pressure Coefficient Distribution”.

30
7.3. Plotting Contour of Pressure

Results > Graphics > Contours > Set Up… Change the parameters as per below and click Display. Save the
picture the same way as you did for Residual History but in this case, name it “CFD Lab2 Contour of
Pressure”.

31
7.4. Plotting Contour of Velocity Magnitude

Results > Graphics and Animations > Contours > Set Up… Change the parameters as per below and click
Display. Save the picture the same way as you did for Residual History but in this case, name it “CFD Lab 2
Contour of Velocity Magnitude”. Zoom in where you can see the airfoil clearly and the change in contour
levels around the airfoil.

32
7.5. Plotting Velocity Vectors at Trailing Edge

Results > Graphics and Animations > Vectors > Set Up… Change parameters as per below and click
Display. Save the picture the same way as you did for Residual History but in this case, name it “CFD Lab 2
Vectors of Velocity at Trailing Edge”. Zoom in on the trailing edge.

33
7.6. Plotting Streamlines Close to Airfoil Surface

Results > Graphics and Animations > Contours > Set Up… Change parameters as per below and click
Display. Save the picture the same way as you did for Residual History but in this case, name it “CFD Lab 2
Streamlines Close to Surface”. Adjust Min and Max values until streamlines cover the airfoil.

34
7.7. Printing Lift and Drag Coefficients

Results > Reports > Forces > Set Up… Change parameters as per below and click Print. This prints out the
drag force. If you change the X parameter to zero and the Y parameter to 1, this prints out the lift force. Save
the coefficients by clicking Write. This creates a text file of what was printed on the screen. Name the file
“Drag Coefficient” or “Lift Coefficient”.

35
8. Exercises

Parametric Studies of Turbulent Flow around an Airfoil

1. Turbulence model k-.

• Figures need to be saved: 1. Time history of residuals (residual vs. iteration number); 2.
Pressure coefficient distribution, 3. Contour of pressure, 4. Velocity vectors, and 5.
Streamlines
• Data need to be saved: lift and drag coefficients

2. Turbulence model k-ω

• Figures need to be saved: 1. Time history of residuals (residual vs. iteration number); 2.
Pressure coefficient distribution, 3. Contour of pressure 4. Velocity vectors, and 5. Streamlines
• Data need to be saved: lift and drag coefficients

3. Questions need to be answered:

Using the figures obtained in exercises 8.1 and 8.2 and those figures you created in this CFD Lab to
answer the following questions and present your answers in your presentation.

1. Which angle of attack produces higher lift/drag coefficients? Why?

2. What is the effect of angle of attack on lift and drag coefficients?
3. Describe the differences of streamline distributions near the trailing edge of airfoil surface for these
two different angles of attack. Do you observe separations for both? If so, does the separation occur
at the same location?
4. Do the two different turbulence models have the same convergence path? If not, which one requires
more iterations to converge.

36