Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (2024)

Airfoil theory is largely governed by potential flow theory. Specifically, a key component of airfoil aerodynamics theory is the combination of:

  • Free-stream velocity elementary flow
  • Line vortex model elementary flow
  • Flow tangency
  • Kutta condition
  • Kutta-Joukowski theorem

These concepts are depicted in the figure below. Let us review all of these concepts to understand the core theory of airfoil aerodynamics.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (1)

Airfoil

Airfoils are the cross-sections of a wing or lifting surface (i.e., propellers and fins). These shapes drive the underlying performance of a lifting surface. As indicated in the Figure below, the shape is described as the cross-section of the lifting surface. These shapes can vary from a high-drag cylinder, used for support, to low-drag shapes (conventionally thought of as an airfoil) associated with wings and high-speed flight.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (2)

Airfoils can be characterized as indicated in the Figure. A few key terms are the angle of attack, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (3), which is defined as the incidence angle of the chord line with respect to the free-stream velocity. The chordline is the line connecting the leading and trailing edges. The chord length (c) is the length of the chordline and the common lengthscale for airfoils. Then the camberline is defined as the line through the mean thickness of the airfoil. Through these characteristics, airfoils can be generalized.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (4)

Aerodynamic loads generally include a lift, drag, and pitching moment. These forces, respectively relate to supporting the aircraft, power, and trim (or straight and level flight). The convention to describe these loads are provided in the Figure below.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (5)

Baseline Aerodynamic Scaling, Parameters, and Concepts

Aerodynamic scaling and nondimensional numbers are critical to aerodynamics. They enable us to scale loads from one condition to another. In general, aerodynamic forces relate to the integration of pressure. Because of this, it is common to define the dynamic pressure, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (6), as

(1) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (7)

This dynamic pressure relates to the overall aerodynamic energy in the flow. Two other critical parameters relate local flow quantities to aerodynamic forces. Recall the momentum equation, given as

(2) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (8)

For aerodynamic forces, we want to focus on the loads on the aerodynamic surfaces themselves. In Fig. 1, we can see the forces on the airfoil are reaction forces to moving the fluid. These are the pressure and viscous forces on the control volume defined on the airfoil surface.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (9)

Aerodynamic forces are reaction forces to the terms on the right-hand side of the equation, integrated over the aerodynamic surface. Note that the pressure/shear forces of the fluid on the airfoil (i.e., the aerodynamic forces) are equal and opposite to the pressure and shear forces the airfoil imposes on the fluid. Hence, the aerodynamic forces (i.e., – of fluid convention) can be described as

(3) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (10)

The force coefficients scale as follows

(4) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (11)

Here, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (12) is a reference area that is defined specifically to applications or conventions. For example, a wing normally uses planform area while a fuselage may use a projected frontal area. The reference A is really not important, what is important is consistent usage of Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (13). We can then compute a force coefficient as follows:

(5) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (14)

Which can be reduced as follows:

(6) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (15)

Note that the surface integral of constant pressure, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (16), is zero (recall from hydrostatics). Hence, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (17). This can be subtracted from Eq. 6, to yield

(7) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (18)

These terms yield local forces in terms of a local pressure coefficient,

(8) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (19)

which yields the surface-normal force. We can also obtain the local skin friction coefficient,

(9) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (20)

which is a normalized, surface-parallel force from viscous forces.

The underlying flow field drives the pressure and viscous forcing. In the context of aerodynamics, there are two key scaling parameters associated with the flow dynamics.The first is the Reynolds number,

(10) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (21)

Conceptually, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (22) defines how important viscous forces are with respect to the inertia of the flow. Here, density, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (23), free-stream velocity, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (24), a length scale (typically chord), Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (25), and dynamic viscosity, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (26), all relate to this nondimensional number.

The second key driving fluid scale is the Mach number,

(11) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (27)

Conceptually, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (28) relates to the speed of a vehicle in comparison to how fast pressure waves move. These two parameters govern kinematic similarity in fluid dynamics, which drive pressure and viscous forces, which ultimately relate to the aerodynamic force response.

Kutta condition

Simple in concept and critical (despite ad-hoc) component to aerodynamics. At its core, the Kutta condition implies either a smooth flow at the trailing edge of an airfoil or a stagnation point at the trailing edge. Conventional wisdom has in the past implied that the Kutta condition is a “viscous patch” that ties potential flow to a physical condition.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (29)

Alternatively, a recent theory derived based on variational theory was developed arguing that rather than a viscous patch, Kutta condition is essentially a result of the incorporation of the momentum equation (recall that the momentum equation is only used for pressure). The authors describe this physical as discussed by Gonzalez and Taha (2022):

“In words, lift is due to (i) momentum conservation, (ii) continuity (the fluid has to fill the given space), and (iii) the body’s hard constraint (i.e. presence of the body inside the fluid), which forces the flow to move around the body without creating void to maintain (ii).”

Regardless, Kutta condition is a physical observation that is used to fix a potential flow solution based on the conservation of mass to a physical condition relevant to airfoils. It is typically done by enforcing one of two conditions at the airfoil trailing edge: (1) a stagnation point, or (2) no lift production (or circulation).

Kutta-Joukowski Theorem:

The concept of the Kutta-Joukowski Theorom drives a correlation between the lift per unit span to be directly proportional to the circulation around the body. Specifically, it states that

(12) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (30)

where Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (31) is the lift per unit span, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (32) is the fluid velocity, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (33) is the circulation, and Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (34) is the free stream velocity. Recall that circulation is defined as

(13) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (35)

Where (+) is in the clockwise direction. What does this mean? First, let’s recall that a flow with Circulation is not the same as a flow with rotation. Rotational flows are described by viscous forces and imply the fluid elements rotate, see for example this link. In the later example in the link are two flows that both have circulation (as the flows rotate like a tornado, hurricane, or vortex), but only one is rotational (where the fluid elements themselves rotate). Circulation is a feature relevant to vortices, and via Kutta-Jukowski, is why line-vortices are critical for lift generaton in potential flow.

Examples: Kutta Joukowski Tho*rom Derviation

Now let’s explore thin airfoil models. The concept is depicted below. Note the distinct difference between our panel code and the model depicted in the first Figure. That is, the airfoil is treated to have no thickness and approximately look like a camberline. At the camberline, flow tangency is enforced which implies that the line vortex strength is established through satisfying flow tangency on the camberline. This model drives the basics of airfoil theory and will be explored in the context of (1) thin-airfoil theory, (2) numerical thin-airfoil theory, and (3) Wiessinger’s approximation. Normally theory is introduced in this order, but let’s flip the ordering to gradually bring in more complex concepts.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (36)

Wiessinger’s approximation

Wiessinger’s approximation is a simplified form of TAT. The overall model is described in the Figure below. Note that compared to the previous model, the line vortex is positioned at the quarter chord and on the chordline. The line vortex induces a velocity as follows
begin{equation}
v_theta=frac{Gamma}{2pi r}.
end{equation}
As indicated in the flow tangency constraint, this induced velocity is evaluated only at the 3/4 chord location and on the chordline (with the slope of the camberline). Note that positioning on the chordline (rather than the camberline) is consistent with a small angle assumption (i.e.,Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (37) and Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (38)). This certainly causes some confusion, but simplifies the math and is not essentially to understanding the model. The model ultimately solves for the vortex strength, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (39), that yeilds lift through Kutta-Joukowski Th. (Eq. (ref{eq:KuttaJouk})) using this value of Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (40).

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (41)

Weissinger’s approximation is much more physically intuitive than the formal thin airfoil theory, so please focus to understand and master this first. The result of this model is a very flexible and is rather simple model of airfoil lift aerodynamics. The model is also extendable to the aerodynamics of a cambered airfoil, flapped airfoils, and airfoils operating together. A sample of its application is provided in the video below.

Examples: Development of Wiessinger’s approximation

Numerical Thin-Airfoil Theory

Now consider expanding Wiessinger’s approximation into many elements that represent the airfoil camberline in detail. This approach is referred to as numerical thin-airfoil theory (and is the discretized representation of formal thin-airfoil theory). The overall model is described in the Figure below. There the extension of Weissinger’s approximation should be seen, that is, we basically break the camber line into individual segments (Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (42)) indicated by the purple sections. In this case, we discretize the camberline into Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (43)-chordal sections. Each of these sections has a length Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (44) and is composed of a line vortex with strength Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (45) at its 1/4 length, while flow tangency is enforced at its 3/4 length. And, similar to Wiessinger’s approximation, the camberline slope is projecting to the chordline. We can then see the final Numerical thin-airfoil theory (TAT) model composed of N-line vortices with N-flow tangency conditions. Similar to the panel code, this develops an NxN algebraic and linear system that is solveable.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (46)

Let’s dig a little deeper into how we would solve this. As indicated in the blowup in the figure above, each panel has a line-vortex and corresponding strength. Each of these vortices induces a flow on each of the panels.The general equation for an induced velocity induced by a line vortex is given as

(14) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (47)

To generalize into the model above, we can describe an induced velocity (Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (48)) which is normal to the chordline by the vortex on panel Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (49) onto panel Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (50). The formal extension is given as:

(15) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (51)

Here, the induced velocity from panel Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (52) onto control point Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (53) is given as

(16) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (54)

Recall, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (55) is positive for a clockwise direction, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (56) is the distance from the vortex position on section Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (57) to the flow tangency point on section Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (58). Hence, the sign flips depending on if Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (59) is upstream or downstream of Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (60). The diagram below indicates the model formulation from the induced flow from the 3rd panel (Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (61)). Note that the induced flow is downwards on the downstream panels, and upwards on the upstream panels.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (62)

In this entire process, remember, we are utilizing an underlying assumption that flow tangency is enforced on the chordline. Without this assumption, there is a demand to expand the analysis to directions that are not chord normal which is much more complex. Ultimately we build a potential flow solution for velocity such that it is given by the free stream + all the contributions from each of the Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (63) sections, which is given as

(17) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (64)

Here, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (65). We then satisfy flow tangency at each control point via

(18) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (66)

which is a constraint that states that the slope of the camberline is equal to the slope of the velocity.

Perhaps the best path to understanding the theory involves working through a stepwise procedure as given below:
Step 1: Define the BCs

  • Free stream flight condition

    (19) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (67)

  • Airfoil discretization which demands (for all Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (68)) dividing the camberline into segments. Segment Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (69) is formed between two points on the camberline: (1) the starting point Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (70), and (2) the endpoint Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (71). Note that this needs to be done for all segments.
  • Airfoil shape which demands (for all Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (72)) by defining the slope for segment Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (73) as

(20) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (74)

Step 2: Compute the control points

  • Line vortex points are (1/4 length of each panel). Use the edge points (Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (75)) of each section:

    (21) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (76)

  • Flow tangency points are (3/4 length of each panel)

    (22) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (77)

    Why are Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (78) and Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (79) not collocated? What happens to the induced velocity at Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (80)? Why can we not handle this numerically?

    Step 3: Compute the panel length and slope or normal vector

  • Provides a scaling of the overall strength. Here, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (81) is per-unit length, hence it is important to think of these as “finite” source sheet segments that are assessed at the quarter chord. The source sheet size can be computed as:

(23) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (82)

Step 4: Now compute the influence coefficients
This is where the discussion of Eqn (15) is involved. Here, these induced velocities are solved to (as a group of all segments) satisfy flow tangency which forms a linear algebraic system of equations. These equations are therefore used to form coefficients for an Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (83) matrix. Recall from Eqn (15) that Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (84) is positive for a clockwise direction hence, then Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (85) (the vector from the control point at Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (86) to the vortex point Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (87)) is positive, the induced flow is downwards Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (88). When the sign flips, ie, the control point is upstream of the vortex, the induced velocity is upwards (Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (89)).Flow tangency is obtained by the requirement that:

(24) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (90)

We can rearrange to obtain:

(25) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (91)

Our A matrix elements come from the right-hand side of Eq. (ref{eq:linSys}). Specifically, each element can be described as

(26) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (92)

Our x vector elements represent the strengh of the vortex sheet segment, which

(27) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (93)

Finally, our Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (94) matrix elements come from enforcing flow tangency, from Eq. (25) this is

(28) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (95)

Recall the approximations and be sure to keep these in mind when utilizing the method/model.

Step 5: Solve for x (or the Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (96)s)

This step involves utilizing an algorithm to solve the linear system, recall Guass-Siedel, iterative methods, etc. Here we just assume we have software that can provide this for us as follows

(29) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (97)

Step 6:Compute the Load
The last step involves extracting the aerodynamic lift. Rather than force integration, here we utilize the Kutta Joukowski Theorom as follows

(30) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (98)

Examples: Numerical thin airfoil theory

Thin-airfoil theory (TAT)

Overview of method:

  • Collapse airfoil to its camberline projected to the chordline
  • Create a vortex sheet on the chordline
  • Find strengths that satisfy flow tangency on the camberline
  • Fully analytic

Below are two lectures on TAT.

Example: Thin Airfoil Tho*ry (TAT) Derivation

Examples: Thin Airfoil Tho*ry (Varient 2)

Potential Flow Represented in Complex Theory

(In progress, draft do not read yet)

Recall, the velocities in terms of Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (99) and Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (100) and their relationship to the Cauchy-Riemann conditions:

(31) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (101)

and

(32) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (102)

We seek to utilize the complex potential funciton

(33) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (103)

where z is the plane of interest. Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (104) hasderivatives given as

(34) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (105)

We can take conventional potential flow theory to formulate the flow over a rotating cylinder in cross flow that generates lift. Here, we sum the Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (106)‘s of the free-stream, cylinder, and vortex as follows (and described in detail here if you do not recall from undergrad):

(35) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (107)

(36) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (108)

Recall the characters from analytic lifting cylinder solutions in the figure below.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (109)

Joukowski transformation is a conformal mapping technique to derive an analytic solution of the flow around a Joukowski airfoil. First, consider the concept as indicated in the Figure below. The transformation maps the flow around a cylinder in cross flow (in the Z plane) onto the A plane where the circle is mapped to a Joukowski airfoil. This new curve preserves the conformal properties of the original circle, which means that it preserves angles (relevant to stream function and velocity potential) locally. As indicated by the lines, the following is mapped:

  • The cylinder surface Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (116) the Joukowski airfoil surface
  • The foremost stagnation point on the cylinder Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (117) the airfoil stagnation point near the leading edge.
  • The aft stagnation point on the cylinder Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (118) the airfoil trailing edge (Kutta Condition).
  • The flow is preserved in this mapping

Hence, at its basis, a Joukowsi airfoil is an airfoil created from “transforming” a cylinder (and its flow) to an airfoil using a special transformation.

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (119)

Let’s now move to the details of the method. First, you must consider the Z-plane, which is a complex plane where the y-axis is imaginary and the x-axis is real. In the context of the complex potential function in the Z-plane, we can define Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (120) as

(37) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (121)

Here, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (122) and Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (123) need to be defined for the flow over a rotating cylinder which is the composition of a uniform flow + doublet + line vortex and is described here. Let’s walk through this in steps.

Starting with a free stream velocity we have a Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (124) and Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (125). These are given as

(38) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (126)

Similarly, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (127) is given as

(39) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (128)

Constructing into the complex potential is as follows:

(40) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (129)

Recall that Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (130), hence,

(41) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (131)

This is a very simple flow and can be expanded to a flow with an angle of attack (Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (132)) with a free-stream velocity Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (133) as

(42) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (134)

Let’s now add the doublet and the line vortex (both at the origin, no angle of attack).

(43) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (135)

Note that the 2nd and 3rd terms on the right-hand side are the doublet and line vortex, respectively. The differential of this gives us velocity, that is, Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (136), which yields what looks like a normal differential as

(44) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (137)

yielding the velocity field.

Examples: Complex potential flow over a rotating cylinder

Try this example matlab code. But sure to use the auto text cleanup.

clear

Uinf=1; %Free stream velocity

R=1; %Cyl radius

Gamma=10;%Ciruclation

%create mesh

x = meshgrid(-5:.1:5);

y = x’;

%create z on imaginary

z = x + i*y;

%clean up inside of cyl (basically let’s focus on only the flow)

for ii = 1:length(x)

for jj = 1:length(y)

if abs(z(ii,jj)) <= R

z(ii,jj) = NaN;

end

end

end

%create our complex potential function F

%%%%%Note, don’t forget to use term-by-term operations on z (use .* and ./)

F = Uinf*(z) + Uinf*R^2./z – i*Gamma/(2*3.14)*log(z);

%create our W

%%%%%Note, don’t forget to use term-by-term operations on z (use .* and ./ and .^2)

W = Uinf – Uinf*R^2./z.^2 + i*Gamma/(2*3.14)./z;

%plot

figure(1)

subplot(2,2,1)

contourf(x,y,real(F))

%xlabel(‘x (m)’)

ylabel(‘y (m)’)

title(‘Velocity Potential (phi)’)

set(gca,’Fontsize’,18)

subplot(2,2,2)

contourf(x,y,imag(F))

%xlabel(‘x (m)’)

ylabel(‘y (m)’)

title(‘Stream Function (psi)’)

set(gca,’Fontsize’,18)

subplot(2,2,3)

contourf(x,y,real(W))

xlabel(‘x (m)’)

ylabel(‘y (m)’)

title(‘u (m/s)’)

set(gca,’Fontsize’,18)

subplot(2,2,4)

contourf(x,y,-imag(W))

xlabel(‘x (m)’)

ylabel(‘y (m)’)

title(‘v (m/s)’)

set(gca,’Fontsize’,18)

figure(2)

contourf(x,y,sqrt(real(W).^2+imag(W).^2).^.5)

xlabel(‘x (m)’)

ylabel(‘y (m)’)

title(‘vmag (m/s)’)

set(gca,’Fontsize’,18)

Next, consider shifting the doublet and the line vortex from the origin, by some distance Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (138). Note that Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (139). We will also recover an angle of attack. The complex potential becomes

(45) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (140)

Note you can simply subtract Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (141) from z. Similar to before, we obtain this differential

(46) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (142)

yielding the velocity field. See the example Matlab code below.

Examples: Lifting cylinder in flow with incidence and offset

clear
Uinf=1; %Free stream velocity
R=1; %Cyl radius
Gamma= 0.25*4*3.14*R*Uinf;%Ciruclation
aoa=25*pi/180; %aoa in rad
mu=0.5 +0.5*i;
%create mesh
x = meshgrid(-5:.1:5);
y = transpose(x);
%create z on imaginary
z = x + i*y;
%clean up inside of cyl (basically lets focus on only the flow)
for ii = 1:length(x)
for jj = 1:length(y)
if abs(z(ii,jj)-mu) <= R
z(ii,jj) = NaN; %comment this and see the high speed if this is not done.
end
end
end
%create our complex potential function F
%%%%%Note, don’t forget to use term-by-term operations on z (use .* and ./)
F = Uinf*((z-mu))*exp(-i*aoa) + Uinf*R^2./(z-mu) – i*Gamma/(2*3.14)*log((z-mu))*exp(i*aoa);
%create our W
%%%%%Note, don’t forget to use term-by-term operations on z (use .* and ./ and .^2)
W = Uinf*exp(-i*aoa) – Uinf*R^2./(z-mu).^2 + i*Gamma/(2*3.14)./(z-mu)*exp(-i*aoa);
%plot
figure(1)
subplot(2,2,1)
contourf(x,y,real(F))
%xlabel('x (m)’)
ylabel('y (m)')
title('Velocity Potential (\phi)')
set(gca,'Fontsize',18)
subplot(2,2,2)
contourf(x,y,imag(F))
%xlabel('x (m)')
ylabel('y (m)')
title('Stream Function (\psi)')
set(gca,'Fontsize',18)
subplot(2,2,3)
contourf(x,y,real(W))
xlabel('x (m)')
ylabel('y (m)')
title('u (m/s)')
set(gca,'Fontsize',18)
subplot(2,2,4)
contourf(x,y,-imag(W))
xlabel('x (m)')
ylabel('y (m)')
title('v (m/s)')
set(gca,'Fontsize',18)
figure(2)
contourf(x,y,sqrt(real(W).^2+imag(W).^2).^.5)
xlabel('x (m)')
ylabel('y (m)')
title('vmag (m/s)')
set(gca,'Fontsize',18)

Now let’s consider the Joukowski transform. The idea here is to transform from the z-plane (where z=x+iy), to the Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (143)-plane (where Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (144)). The mapping is as follows:

(47) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (145)

Examples: Joukowski transform

Try changing the inpu variables (R and mu) and see how things change.

clear

R=1;

mu=-0.2+i*0.1;

%loop over various angles

angle = 0:.01:2*pi;

%define the circle in the z plane

z_circle = R*(cos(angle)+i*sin(angle)) + mu;

%define the Jooukowski airfoil in the w plane

eta_airfoil = z_circle+(R-mu)^2./z_circle;

%plot comparisons

plot(z_circle)

hold on

plot(eta_airfoil)

hold off

Now move to the actual method. The steps as as follows:

  • Identify the stagnation points in the flow over a cylinder.
  • Change Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (146) such that the aft stagnation point falls on the real axis (i.e., no imaginary component)

(48) Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (147)

Theory of Airfoil Lift Aerodynamics – Intermediate Aerodynamic Theory and Analysis (2024)
Top Articles
Basketball Bros On Crazy Games
Closest Chase Bank Branch Near Me
Ucsf Ilab
Express Pay Cspire
Elektrisch koffiezetapparaat Philips CSA240/61 1450 W Zwart 1450 W | bol
Coverwood Terriers For Sale
Paulding County Bus Stop Locator
Solo Player Level 2K23
North Station To Lowell Schedule
I Feel Pretty (2018) | Rotten Tomatoes
Dr. med. Dupont, Allgemeinmediziner in Aachen
Myvetstoreonline.pharmacy
Zenuwbeknelling in de voorvoet (Mortons neuroom)
Adt First Responder Discount
Spectrum Store Downey Photos
Real Estate Transfers Erie Pa
Post-Tribune Obits
Childersburg Busted Mugshots
Maryse Mizanin Nip Slip
Food Stamp System Down
Methodwow
The Quiet Girl Showtimes Near Amc Shirlington 7
How a 1928 Pact Actually Tried to Outlaw War
Cocaine Bear Showtimes Near Harkins Cerritos
Craigslist Gaming Chair
Wok Uberinternal
Toernooien, drives en clubcompetities
Closest Dollar Tree Store To My Location
Ullu Web Series 123
Usc Human Biology
A 100% Honest Review of M. Gemi Shoes — The Laurie Loo
Nsa Panama City Mwr
South Bend Tribune Online
Ssbbw Coomer
Louisiana Physical Therapy Jurisprudence Exam Answers
Chicken Coop Brookhaven Ms
Aerospace Engineering | Graduate Degrees and Requirements
10000 Blaulicht-Meldungen aus Baden-Württemberg | Presseportal
Balmorhea Fishing Resort & Rv Spaces
Basis Independent Brooklyn
Terraria Cement Mixer
Mvsu Canvas
Glyph Of The Trusted Steed
About My Father Showtimes Near Marcus Saukville Cinema
Alj Disposition Data
The Realreal Temporary Closure
Destep: 10 tips voor de scherpste destep-analyse ooit!
Roblox Mod Menu Platinmods
Kamzz Llc
2006 Ford E350 Startrans RV Conversion for sale by owner - Medford, OR - craigslist
XY6020L 6-70V CNC einstellbares stabilisiertes Spannungsnetzteil Konstantspannung Konstantstrom 20A/1200W Buck-Modul Bewertungen
Neuer Extraction-Shooter auf Steam will Escape from Tarkov Konkurrenz machen, wird von echten Militär-Veteranen entwickelt
Latest Posts
Article information

Author: Errol Quitzon

Last Updated:

Views: 6741

Rating: 4.9 / 5 (59 voted)

Reviews: 82% of readers found this page helpful

Author information

Name: Errol Quitzon

Birthday: 1993-04-02

Address: 70604 Haley Lane, Port Weldonside, TN 99233-0942

Phone: +9665282866296

Job: Product Retail Agent

Hobby: Computer programming, Horseback riding, Hooping, Dance, Ice skating, Backpacking, Rafting

Introduction: My name is Errol Quitzon, I am a fair, cute, fancy, clean, attractive, sparkling, kind person who loves writing and wants to share my knowledge and understanding with you.