draft

Climbing and Gliding Caveat

Since this course comprises both aircraft performance and flight mechanics, the nomenclature can be a little confusing when trying to keep things consistent across the disciplines.

Aircraft performance uses simplified models of aircraft mechanics to determine basic relationships for flight. Using these to predict climbing and gliding flight is possible, and will be completed in the following, but some explanation of the nomenclature will be described first.

Flight Mechanics Angles

To describe the aircraft attitude a set of angles are used called Euler angles - these will be fully described in Module 3. The aircraft pitch is defined as the angle between the horizontal and the earth axis, and is given the symbol θ\theta.

To describe the aircraft flight path a set of angles describing the orientation of the aircraft velocity vector with respect to the earth are used. These are the flight path angle γ\gamma, and the track angle, τ\tau.

To describe the aerodynamic incidence, the aerodynamic angles α\alpha and β\beta are used. α\alpha is angle of attack, and β\beta is sideslip.

The angles above will be fully explored and utilised in the later modules, but an appreciation of the longitudinal angles is required here:

Fig. 1: Relationship between longitudinal angles

Aircraft Performance

In aircraft performance, the climb or glide angle is denoted by θ\theta in numerous texts - for example [Anderson:1999AP] - and this angle represents the angle between the aircraft flight direction and the horizontal plane.

Looking at Fig. 1, you will see that this should really be denoted by γ\gamma. In the following, the glide or climb angle will be denoted by γ\gamma, but you should be aware that in most texts you will see it denoted by θ\theta.

In reality, γ=θ\gamma=\theta if and only if (iff) α=0\alpha=0.

The distinction is made because this course later requires us to use θ\theta, and I don’t want you getting confused about what the angles mean.

Since the flight path angle, γ\gamma, is defined as positive in a glide and negative in a climb (positive γ\gamma gives positive WeW_e which is velocity towards the earth), it leads to requiring a negative angle continually in the climb expressions. [Yechout:2014vg] uses a good approach and introduces

γˉ=γ\bar{\gamma}=\left|\gamma\right|

to avoid the negative expression. For the following, γˉ\bar{\gamma} will be used.

Gliding Flight (Unpowered Descent)

The parameters determine above, the lift and hence speeds for minimum thrust and minimum power are important for determination of an aircraft’s glide characteristics.

With no thrust, the aircraft has T=0T=0 but D>0D>0, so there is no longer force equilibrium. The aircraft will descend at the glide angle γˉ\bar{\gamma}. In a steady glide, the horizontal component of weight is equal to the drag.

Fig. 2: Forces on aircraft in a glide

Glide Angle

In Fig. 2, the forces are parallel and perpendicular to the flight path as before, but you should note that this is no longer parallel with the ground. In the horizontal direction, equilibrium gives Wsinγˉ=DW\cdot\sin\bar{\gamma}=D and in the vertical direction, equilibrium gives L=WcosγˉL=W\cdot\cos\bar{\gamma}. Hence the glide angle can be determined from simple trigonometry:

tanγˉ=DL=CDCL\begin{aligned}\tan\bar{\gamma}&=\frac{D}{L}\\ &=\frac{C_D}{C_L}\end{aligned}

So it can be seen that the best/shallowest glide (the smallest γˉ\bar{\gamma}) occurs at the minimum drag to lift ratio - i.e., at CL,mdC_{L,md}.

The glide angle dictates the glide distance, which is very important to be aware of in the cockpit in a glider as this dictates the furthest distance that can be covered along the ground.

The aircraft speed equation can be used with L=WcosγˉL=W\cdot\cos\bar{\gamma} in place of L=WL=W - which shows that the speed that the shallowest glide occurs is not V=Vmd=[BA]14V=V_{md}=\left[\frac{B}{A}\right]^\frac{1}{4} as VmdV_{md} as calculated previously is only valid for cruise.

V=L12ρSCL,md=Wcosγˉ12ρSCL,md=W12ρSCL,mdcosγˉ=Vmdcosγˉ\begin{aligned} V&=\sqrt{\frac{L}{\tfrac{1}{2}\rho\cdot S\cdot C_{L,md}}}\\ &=\sqrt{\frac{W\cdot\cos\bar{\gamma}}{\tfrac{1}{2}\rho\cdot S\cdot C_{L,md}}}\\ &=\sqrt{\frac{W}{\tfrac{1}{2}\rho\cdot S\cdot C_{L,md}}}\cdot\sqrt{\cos\bar{\gamma}}\\ &=V_{md}\cdot\sqrt{\cos\bar{\gamma}}\end{aligned}

Small angle assumption?

The speed for minimum slide is therefore slower than the minimum drag speed in cruise. Compare the ratio of the two speeds for a range of glide angles:

Interactive demo

It can be seen that for glide angles up to around 11 degrees, the best glide speed is 99% of VmdV_{md}. So for most glide performance work, it is assumed that the shallowest glide occurs at VmdV_{md}.

Since most performance gliders have great glide ratios (CLCD\frac{C_L}{C_D}), the above small angle approximation is suitable for gliders and aircraft with reasonable glide ratios. It may not hold for some other aircraft, however, in which case the shallowest glide speed of VmdcosγˉV_{md}\sqrt{\cos\bar{\gamma}} must be used.

Sink rate

Whilst the glide angle dictates the distance along the ground (hence giving the range), the best endurance is given by the slowest overall descent. Thus the sink rate needs to be minimised.

The sink rate, VsinkV_{sink}, is the vertical component of flight speed - Vsink=VsinγˉV_{sink}=V\cdot\sin\bar{\gamma}. From the reasoning used previous regarding the exchange of potential energy for kinetic energy, you should be able to intuit that the minimum sink rate occurs at the minimum power speed.

The sink rate is the vertical component of flightspeed, Vsink=VsinγˉV_{sink}=V\cdot\sin\bar{\gamma} - see Fig. 3

Fig. 3: Glide Velocities

so

Vsink=Vsinγˉ=VDW\begin{aligned}V_{sink} &= V\sin\bar{\gamma}\\ &=V\frac{D}{W}\end{aligned}

from equilibrium and making the small angle assumption for γˉ\bar{\gamma}

VDL=VCDCL\simeq V\frac{D}{L}=V\frac{C_D}{C_L}

using the aircraft speed equation and using the small angle assumption to substitute LWL\simeq W

Vsink=2LρSCLDL2WρSCDCL3/2V_{sink}=\sqrt{\frac{2\,L}{\rho\,S\,C_L}}\frac{D}{L}\simeq\sqrt{\frac{2\,W}{\rho\,S}}\frac{C_D}{C_L^{3/2}}

since the other terms are just constants, or outside of our control, then the minimum sink rate is found at the minimum of CDCL3/2\frac{C_D}{C_L^{3/2}} which can be found from some elementary calculus by substituting CD=CD0+KCL2C_D=C_{D0}+K\cdot C_L^2 and differentiating by CLC_L:

CDCL3/2=CD0+KCL2CL3/2=CD0CL3/2+KCL1/2\frac{C_D}{C_L^{3/2}} = \frac{C_{D0} + K\,C_L^2}{C_L^{3/2}} = \frac{C_{D0}}{C_L^{3/2}} + K\,C_L^{1/2} CDCL3/2CL=32CD0CL5/2+12KCL1/2=0\frac{\partial\frac{C_D}{C_L^{3/2}}}{\partial C_L} = -\frac{3}{2}\frac{C_{D0}}{C_L^{5/2}} + \frac{1}{2}K\,C_L^{-1/2}=0     CL,minsink=3CD0K\implies C_{L,minsink} = \sqrt{\frac{3\,C_{D0}}{K}}

Substitution of this into the aircraft speed equation, Eq. (1), yields the speed for minimum sink

Vminsink=2WρSCL,minsink=2WρS3CD0K=[K3CD0]1/4[2WρS]1/2V_{minsink}=\sqrt{\frac{2\,W}{\rho\,S\,C_{L,minsink}}}= \sqrt{\frac{2\,W}{\rho\,S\,\sqrt{\frac{3\,C_{D0}}{K}}}}=\left[\frac{K}{3\,C_{D0}}\right]^{1/4}\left[\frac{2\,W}{\rho\,S}\right]^{1/2}

which is equal to Vminsink=[B3A]14V_{minsink}=\left[\frac{B}{3\,A}\right]^\frac{1}{4}, or the minimum power condition, which is what was expected from intuition.

Glider Polar

For the pilot, as has been mentioned previously, the actual lift coefficient or angle of attack for a given condition is not usually of relevance. A glider pilot will be keenly aware of their best sink rates, and will desire to know the airspeed(s) at which these occur.

Glider manufacturers tend to provide a chart of sink rate vs. airspeed, which is proportional to the DL\frac{D}{L} ratio against the inverse of the lift coefficient (with some additional exponents thrown in for good measure).

You can see that an increase in the aircraft weight moves the curve along so its tangent remains the same.

Note: The weights for the plot are likely unfeasible for this aircraft but they make for a nice plot

Powered Descent

The procedure for powered descent and determination of descent slope and sink rate is the same but with D=T+WsinγD=T + W\sin\gamma