Flow Over An Airfoil - Pre-Analysis & Start-Up
Flow Over An Airfoil - Pre-Analysis & Start-Up
Flow Over An Airfoil - Pre-Analysis & Start-Up
To begin, we need to tell ANSYS what kind of simulation we are doing. If you look to the left of the start up window, you will see the Toolbox Window. Take
a look through the different selections. We will be using FLUENT to complete the simulation. Load the Fluid Flow (FLUENT) box by dragging and dropping it
into the Project Schematic.
Once you have loaded FLUENT into the project schematic, you are ready to create the geometry for the simulation.
Go to Step 2: Geometry
Go to all FLUENT Learning Modules
MAE 4230
MAE 5230
MAE 5230
MAE 5230
(Advanced)
MAE 4272
Bifurcating Artery
MAE 6690
3D Bifurcating Artery
MAE 6690
MAE 7340
Older Tutorials
These tutorials are for older version of FLUENT and are no longer being updated. They use GAMBIT for meshing.
Supersonic Flow Over a Wedge
Geometry
For users of ANSYS 15.0, please check this link for procedures for turning on the Auto Constraint feature before creating sketches in
DesignModeler.
and select
to launch the Design Modeler. When prompted, select Meters as the unit of measurement.
Airfoil
First, we will create the geometry of the airfoil. In the menu bar, go to Concept > 3D Curve. In the Details View window, click Coordinates File and select
the ellipsis
click
to browse to a file. Browse to and select the geometry file you downloaded earlier. Once you have selected the desired geometry file,
to create the curve. Click
Next, we need to create a surface from the curve we just generated. Go to Concepts > Surfaces from Edges. Click anywhere on the curve you just
created, and select Edges > Apply in the Details View Window. Click
Click
to create a new
to generate the new coordinate system. In the Tree Outline Window, select the new coordinate system you created (defaulted to
For the end points of the arc, first select a point on the vertical axis above the origin (a C symbol will show), then select a point on the vertical axis
below the origin. You should end up with the following:
Now, we need to get rid of necessary lines created by the rectangle. Select Modify in the Sketching Toolboxes window, then select
. Click
the lines of the rectangle the are collinear with the positive and negative vertical axises. Now, select the Dimensions toolbox to dimension the C-Mesh
domain. Click
, followed by the arc to dimension the arc. Assign the arc a value of 12.5. Next, select
vertical axis and the vertical portion of the rectangle in the right half plane. Also assign the horizontal dimension a value of 12.5.
. Click the
Select the layer that corresponds to the airfoil and the airfoil will be highlighted.
Create Quadrants
In the final step of creating the geometry, we will break up the new surface into 4 quadrants; this will be useful for when we want to mesh the geometry.
To begin, select Plane 4 in the Tree Outline Window, and click
. Open the sketching menu, and select
. Draw a line on the vertical
axis that intersects the entire C mesh. Trim away the lines that are beyond the C-Mesh, and you should be left with this:
Next, go to Concepts > Lines from Sketchs. Select the line you just drew and click Base Objects > Apply , followed by
. Now that you
have created a vertical line, create a new sketch and repeat the process for a horizontal line that is collinear to horizontal axis and bisects the geometry.
Now, we need to project the lines we just created onto the surface. Go to Tools > Projection . Select Edges press Ctrl and select on the vertical line we
drew (you'll have to select both parts of it), then press Apply . Next, select Target and select the C-Mesh surface, then click Apply .
Once you click
, you'll notice that the geometry is now composed of two surfaces split by the line we selected. Repeat this process to
create 2 more projections: one projection the line left of the origin onto the left surface, and one projecting the right line on the right surface. When
you're finished, the geometry should be split into 4 parts.
The geometry is finished. Save the project and close the design modeler, as we are now we are ready to create the mesh for the simulation.
Go to Step 3: Mesh
Go to all FLUENT Learning Modules
Mesh
Mapped Face Meshing
First, we will apply a mapped face meshing control to the geometry. In the Outline window, click on Mesh to bring up the Meshing Toolbar. In the Meshing
Toolbar, select Mesh Control > Mapped Face Meshing. Making sure the face selection filter is selected
, select all four faces by holding down the
right mouse button and dragging the mouse of all of the quadrants of the geometry. When all of the faces are highlighted green, in the Details view Window
select Geometry > Apply. Next, select
Edge Sizing
Next, we will apply edge sizing controls to all of the edges of the mesh. To begin, go to Mesh Control > Sizing. Next, click the edge selection filter
Select the following 4 edges buy holding Ctrl and using the left mouse button:
Notice that the element sizes get smaller towards the airfoil. This will give us a better resolution around the airfoil where the flow gets more complicated.
Create a new edge sizing with the same parameters, but choose the 4 remaining straight edges (see figure below). The number of divisions will still be 50,
but now will be selecting a different biasing type by selecting the second Bias option: Bias > - - -----. Again, set the Bias Factor to 150
Edge Bias
It is important to make sure that the edge divisions to this point are biased towards the center of the mesh: otherwise you may run into some
problems later. If your mesh does not match the pictures in the tutorial, make sure to change the parameters of the mesh to make sure that
they do: this might mean choosing different edges for the different biasing types than those outlined in this tutorial.
Finally, create a third edge sizing, and select the rounded edges as the geometry. Again, select Type > Number of Divisions, and change Number of
Divisions to 100. Select Behavior > Hard. This time, we will not bias the edges.
Now, select Mesh > Generate to generate the mesh. It should look like this.
Named Selections
Now will assign names to some of the edges to make creating boundary conditions for the mesh easier. Let's recall the boundary conditions we planned in
the Pre-Analysis Step:
The edges highlighted blue are the inlet, the edges highlighted red are the outlet, and the airfoil is highlighted white in the middle. Now we are ready to
name the sections. In the Outline window, select geometry - this will make seeing the edges a little easier. Again make sure the edge selection tool
is selected. Now, select the two vertical edges on the far right side of the mesh. Right click, and select Create Named Selections. Name the edges outle
t. Next, select the edges that correspond to the inlet of the flow as defined by the picture above. Again, right click and select Create Named Selections
and this time name the selection inlet. Finally, select the two edges making up the airfoil, and name the selection airfoil.
Go to Step 4: Physics Setup
Go to all FLUENT Learning Modules
Physics Setup
Launch the Solver
In this step, we will open fluent and define the boundary conditions of the problem. If you haven't already,
close the meshing window to return to the Project Outline window. Now, click
. This
will load the mesh into FLUENT. Now, double click Setup. The Fluent Launcher Window should open.
Check the box marked Double Precision. To make the solver run a little quicker, under Processing
Options we will select Parallel and change the Number of Processes to 2. This will allow users with a
double core processor to utilize both.
Under Properties, ensure that density is set to Constant and enter 1 kg/m^3 as the density. Click
Change/Create to set the density.
Boundary Conditions
Inlet
Now that the fluid has been described, we are ready to set the boundary conditions of the simulation.
Bring up the boundary conditions menu by selecting Problem Setup > Boundary Conditions. In the
Boundary Conditions window, look under Zones. First, let's set the boundary conditions for the inlet.
Select Inlet to see the details of the boundary condition. The boundary condition type should have
defaulted to velocity-inlet: if it didn't, select it. Now, click Edit to bring up the Velocity-Inlet Window. We
need to specify the magnitude and direction of the velocity. Select Velocity Specification Method >
Components. Remember, we want the flow to enter the inlet at an angle of 6 degrees since the angle of
attack of the airfoil is 6 degrees; thus, the x velocity will be
, and the y velocity will be
.
Specify X-Velocity as 0.9945 m/s and Y-Velocity as 0.1045 m/s. When you have finished specifying
the velocity as entering the inlet at 6 degrees (the same thing as having an angle of attack of 6 degrees),
press OK
Outlet
In the Boundary Conditions window, look under Zones. Select Outlet to see the details of the boundary
condition. The boundary condition type should have defaulted to pressure-outlet: if it didn't, select it.
Click Edit, and ensure that the Gauge Pressure is defaulted to 0. If it is, you may close this window.
Airfoil
In the Boundary Conditions window, look under Zones and select airfoil. Select Type > Wall if it hasn't
been defaulted.
Reference Values
The final thing to do before we move on to solution is to acknowledge the reference values. Go to
Problem Setup > Reference Values. In the Reference Values Window, select Compute From > Inlet.
Check the reference values that appear to make sure they are as we have already set them.
Go to Step 5: Numerical Solution
Go to all FLUENT Learning Modules
Numerical Solution
Methods
First, go to Solution > Solution Methods. Everything in this section should have defaulted to what we want, but let's make sure that under Flow the
selection is Second Order Upwind. If this is the selection, we may move on.
Monitors
Now we are ready to begin solving the simulation. Before we hit solve though, we need to set up some parameters for how Fluent will solve the simulation.
Let's begin by going to Solution > Monitors. In the Monitors Window, look under Residuals, Statistic, and Force Monitors. Select Residuals - Print,Plot and
press Edit. In the Residual Monitors Window, we want to change all of the Absolute Criteria to 1e-6. This will give us some further trust in our solution.
Initial Guess
Now, we need to initialize the solution. Go to Solution > Solution Initialization. In the Solution Initialization Window, select Compute From > Inlet. Ensure
the values that appear are the same values we inputted in Step 5. If the are, initialize the solution by clicking Initialize.
Solve
Once the solution has been initialized, we are ready to solve the simulation. Go to Solution > Run Calculation. Change Number of Iterations to 3000, then
double click Calculate. Sit back and twiddle your thumbs until Fluent spits out a converged solution.
Go to Step 6: Numerical Results
Go to all FLUENT Learning Modules
Numerical Results
Velocity
First, we will look at the velocity vectors of the solution to see if they make intuitive sense. To plot the velocity vectors, go to Results > Graphics and
Animations. In the Graphics and Animations Window, select Vectors and click Set Up.... This will bring up the Vectors Menu.
Make sure the settings of the menu match the figure above: namely Vectors of > Velocity, Color by > Velocity, and set the second box as Velocity
Magnitude. To see the velocity vectors, press Display.
Pressure Contours
To view the pressure contours over the entire mesh, go to Results > Graphics and Animations again, and in the Graphics and Animations Window, select C
ontours.
Click Set Up... to bring up the Contours Menu. Check the box next to Filled. Under Contours Of, ensure that the two boxes that are selected are Pressur
e... and Static Pressure.
Once these parameters are set, press Display to see the pressure contours.
Streamlines
To view the streamlines, keep the Contours window open, and change the Contours Of box to Velocity, and the box below to Stream Function. Change L
evels to 100. Also, uncheck the box marked Auto Range, and set Min(kg/s) to 13.11, and Max(kg/s) to 14.16
Pressure Coefficient
Next, we will plot the pressure coefficient along the surface of the airfoil. Click on Results > Plots to open up the Plots Window. Under Plots, select XY
Plot, and click Set Up.... In the window that pops up, change the settings Y-Axis Function > Pressure, and change the second box to Pressure Coefficient.
Ensure X-Axis Function > Direction Vector. Under Surfaces, select airfoil. See the figure below for help.
When all the settings are correct, press Plot to plot the data to the command window. To save the data to a text file, check the box next to Write to File.
You'll notice that the Plot button has been replaced by a button marked Write..., click it. Change the file type to All Files and save the file name as Pres
sure_Coefficient.txt
We need to set the parameters so drag across the airfoil (keep in mind, which is at an angle) will be displayed. In the Force Reports window change the Di
rection Vector such that X > .9945 and Y > .1045. Click Print to print the drag coefficient to the command window. To print the lift coefficient, in the
Force Reports window change the Direction Vector such that X > -.1045 and Y > .9945. Again, press Print.
Go to Step 7: Verification & Validation
Go to all FLUENT Learning Modules
Exit out of the mesher. First, right click Setup and select Reset. Then click
in the project schematic. Open up the solver, and
solve the simulation using the same solver and boundary conditions (you'll have to input them again), but this time change the number of iterations to 50
00. Again, calculate the force coefficients and graph the pressure coefficient.
Validation
To validate our data, we will take a compare the data from actual experiment.
Unrefined Mesh
Refined Mesh
Experimental Data
Lift Coeffient
0.6315
0.6670
0.6630
Drag Coefficient
0.0122
0.0063
0.0090
Below is a graph displaying the comparing Coefficient of Pressure along the airfoil for the experimental data and the CFD simulation. The data is from
Gregory & O'Reilly, NASA R&M 3726, Jan 1970.
Exercises
Consider the low-speed airflow over the NACA 0012 airfoil at low angles of attack. The Reynolds number based on the chord is Rec = 2.88 10^6. This
flow can reasonably be modeled as incompressible and inviscid.
c. Coefficient of Pressure
Run a simulation for the NACA 0012 airfoil at angles of attack at 0 degrees and 10 degrees for two cases: a mesh with 15000 elements and a mesh with
40000 elements. Plot the pressure coefficient obtained from FLUENT on the same plot as data obtained from experiment The experimental data is from
Gregory & OReilly, NASA R&M 3726, Jan 1970 and is provided here Follow the aeronautical convention of flipping the vertical axis so that negative Cp
values are above and positive Cp values are below. This can be done in MATLAB using set(gca, YDir, reverse);
Conclusions
Comment on the comparison with experiment for the two angles of attack. Also,comment on the effect of mesh refinement. How does the pressure
distribution over the airfoil change on increasing the angle of attack?
Go to Comments
Go to all FLUENT Learning Modules