A flat plate has a chord c = 0.75 m, negligible thickness and a span wise length b = 1 m. It is immersed in an air flow with a free stream velocity U = 40 m/s. The x-axis is parallel to the air flow, while the y-axis is vertical to the stream velocity. The x = 0 coordinate coincides with leading edge of the flat plate and x = c coordinate coincides with trailing edge of the flat plate.

First part

In the first part of the exercise, the flat plate is displaced horizontally, thus α = 0 as seen at cover Figure.

Question 1 – Calculate the vorticity vector at the upper surface at point (c,0,0) understanding that the local friction coefficient is cƒ = 0.003

Re = U∙x/υ = (40∙0.75)/(15∙10-6) = 2000000 = 2∙106 → Turbulent

ω = – ∂u/∂y ; τw = μ∙(∂u/∂y) → (∂u/∂y) = τw

cƒ = τw/q → τw = cƒ∙q

ω = – cƒ∙q/μ = – (0.003∙960)/(18∙10-6) = 160000 Pa/(Pa∙s) = 1.6∙105 s-1

As can be seen, the air flow has enough speed to be considered turbulent. The vorticity is known as – ∂u/∂y. However, it is obtained through the boundary layer friction since cƒ is already given and there is no velocity profile equation. Hence, it is possible to conclude that in this exercise the vorticity is basically proportional to the square of the free stream velocity.

Question 2 – Calculate the mass flow rate of the boundary layer in m²/s

Calculate q’ respective to (x = c , z = 0) and for y > 0 knowing that δ = 0.016 m and the boundary layer thickness is δ* = 0.002. This question is basically formula application, thus:

q’ = U∙(δ – δ*) = 40∙(0.016 – 0.002) = 0.56 m²/s

Question 3 – Develop the analytical expression of the mass flow rate at (x = c , y > 0 , z = 0) for a given velocity profile

The velocity profile is given by:

U(y) = U(y/δ)1/7

δ* = ∫0y≥δ[1 – (U(y)/U)]dy = ∫0δ[1 – (y/δ)1/7]dy

δ* = y|0δ – (1/δ(x)1/7)∙(7/8)∙y8/7|0δ = δ(x) – (1/δ(x)1/7)∙(7/8)∙δ(x)8/7 = δ(x) – (7/8)δ(x) = δ(x)/8

δ* = δ(x)/8 → Analytical expression

This question asked about the analytical expression for the mass flow rate due to the displacement thickness. Since the velocity profile is given, it is just necessary to describe the boundary layer displacement thickness equation. At the end, it is possible to conclude that the flow rate is only proportional to the boundary layer thickness. The velocity profile and the displacement thickness only vary the degree of proportionality.

Second part

In the second part of the exercise, the flat plate exhibits an inclination α = – 0.1 rad and a lift behavior with L = 180 N without the viscous effects.

FIGURE 1

Question 4 – Estimate the induced drag considering the elliptical wing hypothesis

αi = CL/(π∙AR) → CL = αi∙π∙AR

cos(αi) = L/Fi ; sin(αi) = Di/Fi ; Fi = L/cos(αi) = Di/sin(αi)

Di = L∙(sin(αi)/cos(αi)) = L∙tg(αi)

L = CL∙q∙S → CL = L/(q∙S) = αi∙π∙AR → αi = L/(q∙S∙π∙AR)

q =½∙ρ∙U² = ½∙1.2∙40² = 960 Pa

S = 1∙0.75 = 0.75 m²

AR = b/c = 1/0.75

αi = L/(q∙S∙π∙AR) = 180/(960∙0.75∙(1/0.75)∙π) = 0.059683104 ≃ 0.059

Di = L∙tgαi = 180∙tg(0.059) = 10.756 N

It is important to remember that the lift L given is the lift component without the lift generated by the viscous force τ (Figure 1), this is Lƒ = τ∙sin(α). As can be seen, the lift L, which was given, does not account the viscous force. If the total lift was asked, may the induced component would be useful.

Third part

The third section of this exercise, the flat plate is perpendicular to the air flow, α = -π/2 rad. The drag coefficient CD = 1.3 and the pressure coefficient at the base is CP,b = -1.1.

Question 6 – Calculate the average pressure coefficient at the front CP,f

Db = -q∫CP,b∙(n∙i)dS = -q∙CP,b∙c ; n∙i = 1

Df = -q∫CP,f∙(n∙i)dS = q∙CP,f∙c ; n∙i = -1

D = q∙CD∙c

D = Df + Db = q∙CD∙c = q∙CP,f∙c + (-q∙CP,b∙c)

qCD∙c = qCP,f∙c-qCP,b∙c

CD = CP,f – CP,b → CP,f = CD + CP,b = 1.3 – 1.1 = 0.2

Question 7 – Calculate the energetic content of the wake generated by the flat plate per unit length

E0‘ = D = q∙CD∙S = 960∙1.3∙1∙0.75 = 936 Nm/m

Conclusion

The main point of this exercise is the attack angle variation. Different from the other exercises, it was assumed that the flat plate behaves as an elliptic wing. This makes easier to find the induced angle αi, because elliptic wing are the ones which the distribution of the circulation is elliptical and symmetrical along the wing span. Hence, it is quite straight forward to find αi since L was given. What makes a bit confusing is that difference between αi and α, which is illustrated at Figure. Understanding this, it is just application of geometrical and trigonometrical rules. Usually, these exercises assumes the simplification for very small angles.