Main Content

Calculating Best Glide Quantities

This example shows how to perform glide calculations for a Cessna 172 following Example 9.1 in Reference 1 using the Aerospace Toolbox software.

Best glide calculations provide values (velocity and glide angle) that minimize drag and maximize lift-drag ratio (also called the glide ratio).

Aircraft Specifications

The aircraft parameters are declared as follows.

W = 2400;% weight, lbfS = 174;% wing reference area, ft^2;A = 7.38;% wing aspect ratioC_D0 = 0.037;% flaps up parasite drag coefficiente = 0.72;% airplane efficiency factor

Conditions

Set the current aircraft conditions. The bank angle (phi) is zero for this case.

h = 4000;% altitude, ftphi = 0;% bank angle, deg

Convert altitude to meters usingconvlength. The atmospheric calculations in the next step require values in metric units.

h_m = convlength(h,'ft','m');

Calculate atmospheric parameters based on altitude usingatmoscoesa:

[T, a, P, rho] = atmoscoesa(h_m,'Warning');

Convert density from metric to English units usingconvdensity:

rho = convdensity(rho,'kg/m^3','slug/ft^3');

Best Glide Data

Best glide velocity is calculated using the following equation. TAS (true airspeed in feet per second) is the velocity of the aircraft relative to the surrounding air mass.

T A S b g = 2 W ρ S × [ 1 4 C D 0 2 + C D 0 π e A cos 2 ϕ ] 1 4

TAS_bg = sqrt( (2*W)/(rho*S) )...*( 1./( (4*C_D0.^2) + (C_D0.*pi*e*A*(cos(phi)^2)) )).^(1/4);% TAS, fps

Convert velocity from fps to kts usingconvvel. KTAS is true airspeed in knots.

KTAS_bg = convvel(TAS_bg,'ft/s','kts')';

Convert KTAS to KCAS usingcorrectairspeed. KCAS (calibrated airspeed in knots) is the velocity corrected for instrument error and position error. This position error comes from inaccuracies in static pressure measurements at different points in the flight envelope.

KCAS_bg = correctairspeed(KTAS_bg,a,P,'TAS','CAS')';

Best glide angle is calculated using:

sin γ b g = - 4 C D 0 π e A cos 2 ϕ + 4 C D 0

This is the angle between the flight path and the ground that provides the highest L/D ratio.

gamma_bg_rad = asin(-√(4。* C_D0”)。/ (pi*e*A*cos(phi)^2 + 4.*C_D0')) );

Convert glide angle from radians to degrees usingconvang:

gamma_bg = convang(gamma_bg_rad,'rad','deg');

Best glide drag is calculated using:

D m i n = D b g = 1 2 ρ ( T A S b g 2 ) S ( 2 C D 0 ) = - W sin γ b g

D_bg = -W*sin(gamma_bg_rad);

Best glide lift is calculated using:

L b g = L m a x = W cos γ b g = W 2 - D b g 2

L_bg = W*cos(gamma_bg_rad);

Calculate dynamic pressure usingdpressure:

qbar = dpressure([TAS_bg' zeros(size(TAS_bg,2),2)], rho);

Calculate drag and lift coefficients using:

C D b g = D b g q S

C_D_bg = D_bg./(qbar*S);

C L b g = L b g q S

C_L_bg = L_bg./(qbar*S);

Summary of Best Glide Values

Here are the best glide values:

K C A S b g = 7 1 . 9 K C A S

γ b g = - 5 . 3 8 deg

C D b g = 0 . 0 7 4

C L b g = 0 . 7 8 5 9

D b g = 2 2 4 . 9 l b f

L b g = 2 3 8 9 . 4 l b f

Verification

These plots show drag and lift-drag ratio plots for the aircraft as a function of KCAS. The plots are used to verify the best glide calculations.

Set range of airspeeds and convert to KCAS usingconvvelandcorrectairspeed:

TAS = (70:200)';% true airspeed, fpsKTAS = convvel(TAS,'ft/s','kts')';% true airspeed, ktsKCAS = correctairspeed(KTAS,a,P,'TAS','CAS')';% corrected airspeed, kts

Calculate dynamic pressure for new airspeeds usingdpressure:

qbar = dpressure([TAS zeros(size(TAS,1),2)], rho);

Calculate parasite drag using:

D p = 1 2 ρ S C D 0 ( T A S 2 )

Dp = qbar*S.*C_D0;

Calculate induced drag using:

D i = 2 W 2 ρ S π e A 1 ( T A S 2 )

Di = (2*W^2)/(rho*S*pi*e*A).*(TAS.^-2);

Calculate total drag using:

D = D p + D i

D = Dp + Di;

Approximate lift as weight (assuming small glide angle and small angle of attack). At this speed, assuming

C L = 2 π α

and using

C L b g

from above, the angle of attack is about 7 degrees. Adding the flight path angle (i.e. best glide angle) from above shows the fuselage pitch (attitude angle theta) to be about 2 degrees.

L = W;

Plot L/D versus KCAS

As expected, the maximum L/D occurs at approximately the best glide velocity calculated above.

h1 = figure; plot(KCAS,L./D); title('L/D vs. KCAS'); xlabel('KCAS'); ylabel('L/D'); holdonplot(KCAS_bg,L_bg/D_bg,'Marker','o','MarkerFaceColor','black',...'MarkerEdgeColor','black','Color','white'); holdofflegend('L/D','L_{bg}/D_{bg}','Location','Best'); annotation('textarrow',[0.49 0.49],[0.23 0.12],'String','KCAS_{bg}');

Figure contains an axes object. The axes object with title L/D vs. KCAS contains 2 objects of type line. These objects represent L/D, L_{bg}/D_{bg}.

Plot parasite, induced, and total drag curves

Notice the minimum total drag (i.e. D_bg) occurs at approximately the same best glide velocity calculated above.

h2 = figure; plot(KCAS,Dp,KCAS,Di,KCAS,D); title('Parasite, induced, and total drag curves'); xlabel('KCAS'); ylabel('Drag, lbf'); holdonplot(KCAS_bg,D_bg,'Marker','o','MarkerFaceColor','black',...'MarkerEdgeColor','black','Color','white'); holdofflegend('Parasite, D_p','Induced, D_i','Total, D','D_{bg}','Location','Best'); annotation('textarrow',[0.49 0.49],[0.23 0.12],'String','KCAS_{bg}');

Figure contains an axes object. The axes object with title Parasite, induced, and total drag curves contains 4 objects of type line. These objects represent Parasite, D_p, Induced, D_i, Total, D, D_{bg}.

close(h1,h2);

Reference

[1] Lowry, J. T., "Performance of Light Aircraft", AIAA(R) Education Series, Washington, DC, 1999.