Resistance

Storyboard

The flow around a wing leads to the formation of whirlpools that, depending on the shape and angle of the wing with respect to the flow, can cause whirlpools in a section of these. If volume elements around the wing are considered and it is assumed that energy conservation can be madly assumed, different speeds will have to lead to different pressures (Bernoulli) on the surface.

The sum of all pressures on the surface in the direction of flight, both in front of the wing (backward force) and behind the wing (forward force) leads to a total force that we call resistance. To be able to propel a body (airplane / bird) it is necessary for the propulsion to overcome this resistance force.

>Model

ID:(464, 0)



Mechanisms

Iframe

>Top



Code
Concept
Application of coefficient of resistance
Coefficient of resistance of a sphere
Lift force in the flow
Total resistance force
Wing generating lift

Mechanisms

ID:(15180, 0)



Wing generating lift

Description

>Top


Upon observing the average flow around a wing, one can notice that the lines above the wing are longer than those below it. In simplified terms, it is argued that due to this longer path, the speed at the top ($v_t$) is expected to be greater than the speed at the bottom ($v_b$), although both are higher than the speed with respect to the medium ($v$).



If Bernoulli's law is applicable, the difference in velocities would result in a difference in pressures acting on the wing. In particular, if the speed at the top ($v_t$) is greater, its corresponding the pressure on top of wing ($p_t$) would be lower than with the speed at the bottom ($v_b$) and its corresponding the pressure on the bottom of the wing ($p_b$). This would imply the existence of a the lift force ($F_L$) due to the effect of this pressure difference.

However, as seen towards the end of the wing profile (right side), turbulence forms, limiting the applicability of Bernoulli's principle. Specifically, it should be considered that in a certain portion of the wing's perimeter, it may not be applicable, and there will be no contribution to lift.

ID:(11075, 0)



Lift force in the flow

Description

>Top


The object not only generates lift but also creates resistance to the airflow in its surroundings. While the direct application of Bernoulli's principle may not be feasible, we can still understand the kind of effect to expect, even if it needs to be modeled differently. In this context, at the wingtip, the velocity is zero, leading to maximum pressure. Similarly, at the trailing end of the object, the velocity is at its maximum, resulting in minimum pressure. This generates pressure opposing the object's forward motion, corresponding to resistance.



However, it's crucial to note that this argument is only partially correct. In particular, turbulence is generated at the wing's trailing edge, which complicates the high-speed argument and challenges the application of Bernoulli's principle.

Following the modeling approach used for lift, we can assume that the Kutta-Joukowski theorem, in which the lift force ($F_L$) with the density ($\rho$), the speed with respect to the medium ($v$), and the aerodynamic circulation ($\Gamma$), is:

$ \displaystyle\frac{ F_L }{ L } = - \rho v \Gamma$

We can assume it to be vectorial, meaning it has a vertical component explaining lift and a horizontal component modeling resistance.

ID:(11076, 0)



Coefficient of resistance of a sphere

Description

>Top


The coefficient of resistance $C_d$ is often a function of the Reynolds number $Re$. In the case of a sphere, the coefficient of resistance takes the form:

In the range of low Reynolds numbers, the coefficient of resistance is inversely proportional to the velocity $1/v$, which means that in this range, the drag force is proportional to the velocity (Stokes' law).

In the range of high Reynolds numbers, the coefficient of resistance becomes constant, resulting in a drag force that is proportional to the square of the velocity. However, it is important to note that there is a sudden drop in the coefficient, indicating a situation where the turbulent zone decreases.

ID:(1169, 0)



Application of coefficient of resistance

Description

>Top


If the goal is to reduce the resistance of a ball, it is necessary to increase the Reynolds number to take advantage of the decrease in the coefficient of resistance. This can be achieved by introducing grooves, such as those found on a golf ball. These grooves act as small dividers that keep the layer of air attached to the ball, allowing the flow to remain more linear for a longer period of time and reducing the area that generates resistance.

ID:(1170, 0)



Total resistance force

Description

>Top


When the wing is tilted, it not only generates lift but also produces a component of drag, as the direction of lift is orthogonal to the wing's surface. If we graph the lift force ($F_L$), the resistance force ($F_W$), and the angle of attack of a wing ($\alpha$), we obtain:

ID:(9578, 0)



Model

Top

>Top



Parameters

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$m$
m
Body mass
kg
$C_W$
C_W
Coefficient of resistance
-
$\rho$
rho
Density
kg/m^3
$g$
g
Gravitational Acceleration
m/s^2
$c$
c
Proportionality constant coefficient sustainability
1/rad
$S_w$
S_w
Surface that generates lift
m^2
$S_p$
S_p
Total object profile
m^2
$L$
L
Wing span
m

Variables

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
$\Gamma$
Gamma
Aerodynamic circulation
m^2/s
$\alpha$
alpha
Angle of attack of a wing
rad
$\alpha_s$
alpha_s
Angle required for lift
rad
$\rho$
rho
Density
kg/m^3
$F_L$
F_L
Lift force
N
$F_W$
F_W
Resistance force
N
$C_L$
C_L
Simple Model for Sustainability Coefficient
-
$v$
v
Speed with respect to the medium
m/s
$F_R$
F_R
Total resistance force
N

Calculations


First, select the equation: to , then, select the variable: to

Calculations

Symbol
Equation
Solved
Translated

Calculations

Symbol
Equation
Solved
Translated

Variable Given Calculate Target : Equation To be used




Equations

#
Equation

$ C_L = c \alpha $

C_L = c * alpha


$ \displaystyle\frac{ F_L }{ L } = - \rho v \Gamma$

F_ L / L = - rho * v * Gamma


$ F_L =\displaystyle\frac{1}{2} \rho S_w C_L v ^2$

F_L = rho * S_w * C_L * v ^2/2


$ F_R = F_W \cos \alpha + F_L \sin \alpha $

F_R = F_W *cos( alpha )+ F_L *sin( alpha )


$ F_R = \displaystyle\frac{1}{2} \rho S_p C_w v ^2 + \displaystyle\frac{2 m ^2 g ^2}{ c ^2 S_w \rho }\displaystyle\frac{1}{ v ^2}$

F_R = rho * S_p * C_w * v ^2/2+2* m ^2* g ^2/( c ^2* S_w * rho * v ^2)


$ F_W =\displaystyle\frac{1}{2} \rho S_p C_W v ^2$

F_W = rho * S_p * C_W * v ^2/2


$\cos\alpha\sim 1$

cos(alpha)=1


$\sin\alpha\sim\alpha$

sin(alpha)=alpha

ID:(15185, 0)



KuttaJoukowski theorem

Equation

>Top, >Model


From the works of Kutta [1] and Joukowski [2], a theorem was developed that shows the association between the aerodynamic circulation ($\Gamma$) and the lift force ($F_L$) through the wing span ($L$), the density ($\rho$), and the speed with respect to the medium ($v$) as follows:

$ \displaystyle\frac{ F_L }{ L } = - \rho v \Gamma$

$\Gamma$
Aerodynamic circulation
$m^2/s$
10203
$\rho$
Density
$kg/m^3$
10182
$F_L$
Lift force
$N$
6120
$v$
Speed with respect to the medium
$m/s$
6110
$L$
Wing span
$m$
6337

[1] "Über die Aufgabe der Flügeltheorie und ein neues Verfahren zur Herleitung derselben." (On the task of wing theory and a new method for its derivation.), Martin Wilhelm Kutta, Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse (1902)

[2] "Über die Erhaltung des Luftkreises um ein Profil." (On the conservation of the air circle around a profile.), Nikolai Zhukovsky, Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse (1904)

ID:(1166, 0)



Lift force

Equation

>Top, >Model


To generate higher pressure below than above the wing and generate lift, Bernoulli's principle is employed, correcting for the lack of energy density conservation using ($$). The pressure over the wing, the lift force ($F_L$), can be estimated using the density ($\rho$), the surface that generates lift ($S_w$), the coefficient of lift ($C_L$), and the speed with respect to the medium ($v$) through the following formula:

$ F_L =\displaystyle\frac{1}{2} \rho S_w C_L v ^2$

$\rho$
Density
$kg/m^3$
5342
$F_L$
Lift force
$N$
6120
$C_L$
Simple Model for Sustainability Coefficient
$-$
6164
$v$
Speed with respect to the medium
$m/s$
6110
$S_w$
Surface that generates lift
$m^2$
6117

The lift force ($F_L$), along with the wing span ($L$), the density ($\rho$), the wing top speed factor ($c_t$), the wing bottom speed factor ($c_b$), the upper wing length ($l_t$), the bottom wing length ($l_b$), and the speed with respect to the medium ($v$), is found in

$ F_L = \rho L ( c_b l_b - c_t l_t ) v ^2$



If we consider the surface that generates lift ($S_w$), given by the wing span ($L$), the upper wing length ($l_t$), and the bottom wing length ($l_b$),

$ S_w = \displaystyle\frac{1}{2} L ( l_t + l_b )$



and for the coefficient of lift ($C_L$), defined as

$ C_L = 4\displaystyle\frac{ c_t l_t - c_b l_b }{ l_t + l_b }$



we obtain

$ F_L =\displaystyle\frac{1}{2} \rho S_w C_L v ^2$

ID:(4417, 0)



Resistance force

Equation

>Top, >Model


The resistance force ($F_W$) kann mit the density ($\rho$), the coefficient of resistance ($C_W$), the total object profile ($S_p$) und the speed with respect to the medium ($v$) entsprechend berechnet werden folgende Formel:

$ F_W =\displaystyle\frac{1}{2} \rho S_p C_W v ^2$

$C_W$
Coefficient of resistance
$-$
6122
$\rho$
Density
$kg/m^3$
5342
$F_W$
Resistance force
$N$
6124
$v$
Speed with respect to the medium
$m/s$
6110
$S_p$
Total object profile
$m^2$
6123

Similarly to how the equation for the lift force ($F_L$) was derived using the density ($\rho$), the coefficient of lift ($C_L$), the surface that generates lift ($S_w$), and the speed with respect to the medium ($v$)

$ F_L =\displaystyle\frac{1}{2} \rho S_w C_L v ^2$



in this analogy, what corresponds to the surface that generates lift ($S_w$) will be equivalent to the total object profile ($S_p$) and the coefficient of lift ($C_L$) to the coefficient of resistance ($C_W$), thus the resistance force ($F_W$) is calculated:

$ F_W =\displaystyle\frac{1}{2} \rho S_p C_W v ^2$

The drag coefficient is measured and, in turbulent flows over aerodynamic bodies, values are generally found around 0.4.

ID:(4418, 0)



Total resistance force calculation

Equation

>Top, >Model


The total force of resistance is composed of the horizontal components of the wing profile's resistance force $F_W$ and the lift force $F_L$, which can be calculated from the angle of attack $\alpha$:

$ F_R = F_W \cos \alpha + F_L \sin \alpha $

$\alpha$
Angle required for lift
$rad$
6167
$F_L$
Lift force
$N$
6120
$F_w$
Resistance force
$N$
6124
$F_R$
Total resistance force
$N$
8480

The horizontal component of the lift force corresponds to the force $F_L$ multiplied by the sine of the angle of attack $\alpha$:

$F_L \sin\alpha $



And the horizontal component of the drag force corresponds to the force $F_W$ multiplied by the cosine of the angle of attack $\alpha$:

$F_W \cos\alpha $



Therefore, the total resistance force can be calculated as:

$ F_R = F_W \cos \alpha + F_L \sin \alpha $

ID:(9579, 0)



Lift constant

Equation

>Top, >Model


From measurements, it is concluded that the lift coefficient $C_L$ is proportional to the angle of attack $\alpha$:

$ C_L = c \alpha $

$\alpha_s$
Angle required for lift
$rad$
6167
$c$
Proportionality constant coefficient sustainability
$1/rad$
6165
$C_L$
Simple Model for Sustainability Coefficient
$-$
6164

After a certain angle, the curve decreases until it reaches zero. This is because beyond that critical angle, the vortices fully cover the upper surface of the wing, leading to a loss of lift. This phenomenon is known as \"stall\".

ID:(4441, 0)



Sinus for small angles

Equation

>Top, >Model


For small angles, the sine function can be approximated by a straight line that passes through the origin. If the angle of attack of a wing ($\alpha$) is expressed in radians, the slope of this line is equal to one, and we get:

$\sin\alpha\sim\alpha$

$\alpha$
Angle of attack of a wing
$rad$
6121

The sine can be calculated using a polynomial of the form:

$\sin \alpha = \alpha - \displaystyle\frac{1}{3!}\alpha^3 + \displaystyle\frac{1}{5!}\alpha^5 \ldots$



For small values of the angle of attack of a wing ($\alpha$) ($\alpha \ll 1$), the terms with higher powers are negligible, and we obtain:

$\sin\alpha\sim\alpha$

ID:(9580, 0)



Cosine for small angles

Equation

>Top, >Model


For small values of the angle of attack of a wing ($\alpha$), the cosine function can be approximated by an inverted parabola that passes through the origin. If the angle is expressed in radians and is approximately zero, we get:

$\cos\alpha\sim 1$

$\alpha$
Angle required for lift
$rad$
6167

The cosine can be calculated using a polynomial of the form:

$\cos \alpha = 1 - \displaystyle\frac{1}{2!}\alpha^2 + \displaystyle\frac{1}{4!}\alpha^4 \ldots$



For small values of the angle of attack of a wing ($\alpha$) ($\alpha \ll 1$), the terms with higher powers are negligible, and we obtain:

$\cos\alpha\sim 1$

ID:(14473, 0)



Total force of resistance

Equation

>Top, >Model


To calculate the total resistance force ($F_R$), we assume small angles and consider a situation where the angle is such that it maintains the body mass ($m$). Using this approximation and the variables the coefficient of lift ($C_L$), the coefficient of resistance ($C_W$), the surface that generates lift ($S_w$), the total object profile ($S_p$), the gravitational Acceleration ($g$), the proportionality constant coefficient sustainability ($c$), the density ($\rho$), and the speed with respect to the medium ($v$), we obtain the following expression:

$ F_R = \displaystyle\frac{1}{2} \rho S_p C_w v ^2 + \displaystyle\frac{2 m ^2 g ^2}{ c ^2 S_w \rho }\displaystyle\frac{1}{ v ^2}$

$m$
Body mass
$kg$
6150
$C_w$
Coefficient of resistance
$-$
6122
$\rho$
Density
$kg/m^3$
5342
$g$
Gravitational Acceleration
9.8
$m/s^2$
5310
$c$
Proportionality constant coefficient sustainability
$1/rad$
6165
$v$
Speed with respect to the medium
$m/s$
6110
$S_w$
Surface that generates lift
$m^2$
6117
$S_p$
Total object profile
$m^2$
6123
$F_R$
Total resistance force
$N$
8480

Using the relationships of the total resistance force ($F_R$) with the lift force ($F_L$), the resistance force ($F_W$), and the angle of attack of a wing ($\alpha$):

$ F_R = F_W \cos \alpha + F_L \sin \alpha $



we can calculate using the resistance force with the density ($\rho$), the coefficient of resistance ($C_W$), the total object profile ($S_p$), and the speed with respect to the medium ($v$):

$ F_W =\displaystyle\frac{1}{2} \rho S_p C_W v ^2$



and the lift force with the surface that generates lift ($S_w$) and the coefficient of lift ($C_L$):

$ F_L =\displaystyle\frac{1}{2} \rho S_w C_L v ^2$



using the relationship for the coefficient of lift ($C_L$) with the proportionality constant coefficient sustainability ($c$):

$ C_L = c \alpha $



using the relationship for the sine of the small angle of attack $\alpha$:

$\sin\alpha\sim\alpha$



and the cosine:

$\cos\alpha\sim 1$



with the condition to balance the weight of the bird or aircraft for the body mass ($m$) and the gravitational Acceleration ($g$):

$ \alpha =\displaystyle\frac{2 m g }{ c \rho S_w }\displaystyle\frac{1}{ v ^2}$



we obtain:

$ F_R = \displaystyle\frac{1}{2} \rho S_p C_w v ^2 + \displaystyle\frac{2 m ^2 g ^2}{ c ^2 S_w \rho }\displaystyle\frac{1}{ v ^2}$

ID:(4546, 0)