1SOLUTIONSHomework #3AAE 334Assigned: Tuesday, June 30, 2020Due: Tuesday, July 7, 2020 by midnightSubmit as a single PDF file in GradescopeProblem 1. [15 points]Assuming incompressible flow, consider an uncambered, rectangular wing with span ?? andchord ?? as shown in the figure below. Model this wing using a single horseshoe vortex, with itsbound vortex at the quarter-chord line of the wing (?? =
...[Show More]
1
SOLUTIONS
Homework #3
AAE 334
Assigned: Tuesday, June 30, 2020
Due: Tuesday, July 7, 2020 by midnight
Submit as a single PDF file in Gradescope
Problem 1. [15 points]
Assuming incompressible flow, consider an uncambered, rectangular wing with span ?? and
chord ?? as shown in the figure below. Model this wing using a single horseshoe vortex, with its
bound vortex at the quarter-chord line of the wing (?? = ??⁄4), as shown in the figure.
(a) Use the Biot-Savart law to calculate the velocity induced by the horseshoe vortex at the
control point, which is located at �3 4 ??, 0,0�. Remember that the bound vortex is finite and
the trailing vortices are semi-infinite.
SOLUTION
2
3
(b) Determine the strength of the vortex required to enforce the flow tangency condition (no
flow through the airfoil) at the control point for a given angle of attack ??. Assume that the
thickness of the wing and the angle of attack are small.
SOLUTION
4
(c) Find the lift on the wing for an angle of attack ??.
SOLUTION
(d) Find the slope of the lift curve ??????⁄???? for a wing with aspect ratio AR = 10. Compare the
computed lift curve slope to that given by the lifting line theory for an elliptically loaded
wing with the same aspect ratio.
SOLUTION
5
6
Problem 2. [10 points]
The Supermarine Spitfire shown in the figure had elliptical cord distribution and no twist, resulting
in elliptical wing loading. The airplane had a maximum velocity of ??∞ = ???????? = 362 mph at an
altitude of 18,500 ft. Its weight was 5820 lb, wing area was 242 ft2, and wing span was 36.1 ft.
Given: ??
∞ = ???????? = ?????? mph = 531 ft/s at an altitude H = 18,500 ft; W= 5820 lb, S = 242
ft2, and b = 36.1 ft.
Engine power = 1050 hp, and the propeller efficiency can be assumed to be 0.9. Since 1 hp
= 550 ft∙lb/s, the thrust power is P = 1050 x 550 x 0.9 ft∙lb/s = 519750 ft∙lb/s.
For Part (c), taper ratio =0.4. For Part (d), landing velocity ??∞ = ???? mph.
(a) Calculate the induced drag coefficient of the Spitfire at ???????? = 362 mph and 18,500 ft
altitude.
SOLUTION
For level steady flight at maximum altitude, T=W. Density at 18,500 ft is ?? = ??. ????????????????
slug/ft3. Using these values, we can calculate the lift coefficient and then the induced drag
coefficient:
7
(b) What fraction of the total drag is the induced drag at these conditions? (To calculate the
total drag, note that in steady, level flight the drag ?? equals the thrust ?? which, in turn, is
related to the power ??: ?? = ????∞. The Spitfire had a supercharged Merlin engine that
produced 1050 hp at 18,500 ft, and the propeller efficiency can be assumed to be 0.9. Since
1 hp = 550 ft∙lb/s, the thrust power is ?? = 1050 × 550 × 0.9 ft∙lb/s = 519750 ft∙lb/s.)
SOLUTION
Total drag:
Thus, the induced drag is only 4.5% of the total drag in these cruise conditions.
(c) If the elliptical wing of the Spitfire were replaced by a tapered wing with a taper ratio of
0.4, and everything else remained the same, what would the induced drag coefficient be?
Compare this with the actual Spitfire induced drag coefficient obtained in Part (a) and
comment on the effect of planform shape on the drag at high speed. For the tapered wing,
you can use the figure on slide 14 of Lecture 15 to estimate the induced drag factor ??.
SOLUTION
For the taper ratio ????⁄???? = ??. ??, we can estimate the induced drag factor ?? from the figure
on slide 74 of the Week 3 notes: ?? = ??. ??????. Therefore,
This differs by less than 2% from the induced drag coefficient for elliptical loading. Thus,
tapered wings can be used instead of difficult-to-make elliptical ones.
8
(d) Assume that the Spitfire’s landing velocity at sea level is 70 mph. Calculate the induced
drag coefficient in this case and compare the result with that in high-speed case of Part (a).
Comment on the relative importance of the induced drag coefficient at low speed versus
that at high speed.
SOLUTION
For landing, ??∞ = ???? mph = 102.67 ft/s, and the density at sea level is ?? = ??. ????????????
slug/ft3. We also still have L=W.
which is much higher than for the high-speed cruise. Then:
This induced drag coefficient is more than 2 orders of magnitude higher than that in cruise
flight. Thus, in contrast with cruise conditions, the induced drag is a major component of
the total drag during landing.
9
Problem 3. [45 points]
You will run XFLR5 to perform lifting line analysis on the wing from the general aviation aircraft
shown on the attached sheet. The instructions for performing a wing analysis are given in the
instructions file in the Homework 6 folder in Brightspace. The airplane has a wingspan of 36 ft. 1
in. The root chord length is ???? = 5 ft. 6 in. The wing has a constant chord length over the center
section, which is 16 ft. 9 in. long. Then the wing is tapered to the wing tip, where the chord length
is 4 ft. 0 in. We will assume the leading edge of the entire wing is straight (no sweep), except as
noted below. The airfoil section is a NACA 2412 for the entire wing. We will consider the design
cruise conditions as follows:
• Steady level flight at an altitude of 8,500 ft. (assume a standard atmosphere and use
Sutherland’s law to determine the viscosity of air at the relevant temperature)
• Air speed of 124 kts.
• Weight of 2550 lbs.
(a) Determine the wing area (in ft2 and m2), wing aspect ratio, freestream Mach number, and
wing lift coefficient needed to support the weight of the plane (using the wing planform
area as the reference area and ignoring the tail). Is the freestream Mach number below 0.3?
(It should be. Therefore, we will analyze the flow as incompressible, setting the Mach
number to 0.) Determine the chord length Reynolds number of the wing root airfoil section
and the chord length Reynolds number of the wing tip airfoil section.
SOLUTION
10
11
(b) We first need to run a batch airfoil analysis, as outlined in the instructions. (Note: use a
refined geometry with 201 panels.) Since we will want to run the wing analysis at other
flight conditions, we need a wider range of conditions than just the two Reynolds number
found in Part (a). Run a batch analysis for Reynolds numbers ranging from 1,000,000 to
7,500,000 with an increment of 500,000. Set the boundary layer trip locations on the upper
and lower surface to ??????⁄?? = 0.20. This will provide nearly fully turbulent boundary layers
at all conditions. (While this increases viscous drag, it helps avoid or reduce boundary layer
separation.) Run a sequence of angles of attack from ?? = -18° to 24° with an increment
of ∆?? = 0.1°. Save the airfoil polars as described in the instructions. Plot the lift curve,
drag polar, and lift to drag ratio from the batch analysis. How do the maximum lift
coefficient, minimum drag coefficient and maximum lift to drag ratio change with
Reynolds number?
SOLUTION
12
13
14
(c) Follow the instructions and create the wing geometry using the NACA 2412 airfoil. Print
a planform view of your wing and a perspective view.
SOLUTION
(d) Run a Type 2 wing analysis, entering the mass of the airplane and the other needed flight
parameters. (Use the atmospheric properties you determine from the Standard Atmosphere
Tables in Anderson and Sutherland’s law for viscosity (available on-line) rather than the
conditions available from XFLR5.) When you enter the mass, set the center of gravity to
?? = 18.25 in. downstream of the wing root leading edge. Use a range of angles of attack
from ?? = 0° to 12° with an increment of ∆?? = 0.1°. (There may be angles of attack for
which the airfoil data is not available and the solution will not converge, but we will have
the range of angles of attack we need.) Plot the wing lift curve, drag polar and lift to drag
ratio. Determine the maximum lift to drag ratio of the wing. How does this value compare
15
to that of the airfoil? If one accounts for the drag of the fuselage, horizontal and vertical
tails, and other appendages, such as antennae, the drag coefficient would more than double.
Assuming the drag of the entire airplane is twice that of the wing, determine the maximum
lift to drag ratio of the airplane? (This is a more realistic value for the lift to drag ratio of
an airplane.) At what angle of attack does the maximum lift to drag ratio of the wing occur?
What is the wing lift coefficient at this angle of attack? How does this lift coefficient
compare to the one computed in Part (a)? Use a finite difference approximation to compute
the lift curve slope near the angle of attack for maximum lift to drag ratio. (Do this by
computing ??????⁄???? ≈ ∆????⁄∆?? using values of ?? on either side of the desired angle of
attack.) How does the value of the lift curve slope compare to the theoretical lift curve
slope for an elliptically loaded wing with NACA 2412 airfoil sections and the aspect ratio
found in Part (a)? (The theoretical lift curve slope of an elliptically loaded wing is given by
where ??0 is the airfoil lift curve slope, which we will take to be ??0 = 6.588 per radian.)
SOLUTION
16
17
(e) On the right side of the toolbar (near the top of the window) is a dropdown list of the angles
of attack for which results are available. Select the angle of attack for the maximum lift to
drag ratio of the wing. Then click on <OpPoint> <Current OpPoint> <Properties>. A
window will pop up with tabulated data. The variables of interest to us are the following:
CD, the drag coefficient
VCD, the viscous drag coefficient
ICD, the induced drag coefficient
Cm, the moment coefficient about the moment reference center (the c.g.)
Tabulate the values of these quantities at this angle of attack and at ?? = 10°. Compute
the ratio of the induced drag coefficient to the total drag coefficient
for the two angles of attack. Compute the span efficiency factor ?? from ?????? = ????2⁄(????????).
Discuss how angle of attack affects the ratio of induced drag to total drag and the value of
the span efficiency factor. Compute and tabulate the total dimensional drag for these two
angles of attack. How do they compare? What are the corresponding flight speeds for these
two angles of attack? How do they compare? The Breguet range equation relates the range
of an aircraft to its velocity, lift to drag ratio, specific impulse of the propulsion system,
and the initial and final weights of the airplane (due to the consumption of fuel). (See
http://web.mit.edu/16.unified/www/FALL/thermodynamics/notes/node98.html). From
this equation velocity times lift to drag ratio is important. What is the ratio of the product
of velocity times lift to drag ratio for the two angles of attack (with the product for ?? =
18
10° divided by the product for the angle of attack corresponding to the maximum lift to
drag ratio)?
SOLUTION
19
(f) Select the OpPoint view from the toolbar. Choose <Graphs> <Graph 4> for the lift
distribution (“Local Lift”), ??′(??). Choose the angle of attack that corresponds to maximum
lift to drag ratio of the wing. Add the elliptic load distribution curve following the
instructions. Plot the comparison and discuss similarities and differences between the load
distribution for this wing and an elliptic load distribution. Look at the spanwise load
distribution for ?? = 10°. Is there much of a difference due to the change in angle of attack?
SOLUTION
(g) Assume the horizontal tail is made of thin symmetric airfoil sections. Then the center of
pressure is the same as the aerodynamic center, which is at the tail quarter chord location.
Assume the quarter chord locations of the tail airfoils are lined up at the same axial position,
although that is not exactly the case from the drawing. The tail aerodynamic center is
located 15 ft. 6 in. behind the wing leading edge. What dimensional value of lift is needed
on the tail to balance the moment about the c.g. at the angle of attack that corresponds to
the maximum lift to drag ratio? What value is needed at ?? = 10°.
SOLUTION
Problem 4. [15 points]
(Make sure your project from problem 1 is saved.)
Edit the wing and save it with a different name. To do this click on <Plane> <Current Plane>
<Duplicate...>. In this problem you will include a small amount of twist (with negative values
corresponding to washout of the wingtips). Entering a value of twist introduces a linear distribution
of wing twist from one end of a wing section to the next. You should keep the wing root at zero
twist angle. Set the wing tip to have -1.5° of twist. (Set the twist at the junctions of the middle
and outer wing sections so that the twist is proportional to the spanwise location, i.e., so that the
twist varies linearly and smoothly from the wing root to the wing tip.) Once you have created the
new wing, run the analysis. How does the maximum lift to drag ratio compare to the wing from
Problem 3? At what angle of attack does this occur? Does this make sense, given that the wing is
twisted? Plot the spanwise load distribution at this angle of attack, with the elliptic curve. Is the
load distribution closer or further away from the elliptic distribution? Compute the span efficiency
factor ?? from ?????? = ????2⁄(????????) at this angle of attack and compare it to the value found in
Problem 3 at the angle of maximum lift to drag ratio. What is the ratio of the induced drag to the
total drag at this angle of attack and how does that compare to what was found in Problem 3?
Because the induced drag coefficient varies with lift coefficient, it is better to compare the wing
performance at the same lift coefficient. Use linear interpolation of the data that is available for
the twisted wing to determine the angle of attack where the lift coefficient equals the value from
the table in Problem 3 for the non-twisted wing at its maximum lift to drag ratio. Compute the ratio
of the induced drag to the total drag using the interpolated data for the twisted wing at this angle
of attack where the ???? values match and describe how this value of the ratio compares to the value
found in Problem 3.
SOLUTION
Problem 5. [15 points]
Create a new wing without twist but with rearward wing sweep. Do this by defining the offset
values of the wing sections so that the leading edge is swept back by 30°. Since LLT (lifting-line
theory) does not account for wing sweep, choose to run the analysis with the vortex lattice method
(use VLM1). Do the same tasks and answer the same questions as in Problem 4, and discuss the
differences rearward wing sweep makes on the results.
SOLUTION
[Show Less]