Journal of Engineering Science and Technology

Vol. XX, No. Y (Year) PPP - QQQ

© School of Engineering, Taylor’s University

THIS IS THE TEMPLATE YOU USE TO

FORMAT AND PREPARE YOUR MANUSCRIPT

ANNA A. LEE*, FRANK B. KONG

School of Engineering, Taylor’s University, Taylor's Lakeside Campus,

No. 1 Jalan Taylor's, 47500, Subang Jaya, Selangor DE, Malaysia
*Corresponding Author: [email protected]

Abstract
An investigation has been made to predict the effects of forebody and afterbody
shapes on the aerodynamic characteristics of several projectile bodies at
supersonic speeds using analytical methods combined with semi-empirical
design curves. The considered projectile bodies had a length-to-diameter ratio
of 6.67 and included three variations of forebody shape and three variations of
afterbody shape. The results, which are verified by comparison with available
experimental data, indicated that the lowest drag was achieved with a cone-
cylinder at the considered Mach number range. It is also shown that the drag
can be reduced by boattailing the afterbody. The centre-of-pressure assumed a
slightly rearward location for the ogive-cylinder configuration when compared
to the configuration with boattailed afterbody where it was the most forward.
With the exception of the boattailed afterbody, all the bodies indicated inherent
static stability above Mach number 2 for a centre-of-gravity location at about
40% from the body nose.
Keywords: Projectile, Forebody and afterbody, Supersonic speed, Aerodynamics,
Next keyword.

1. Introduction
The shape of a projectile is generally selected on the basis of combined
aerodynamic, guidance, and structural considerations. The choice of seeker,
warhead, launcher, and propulsion system has a large impact on aerodynamic
design [1]. Consequently, various configurations have evolved, each resulting
from a series of design compromises. During supersonic flight, the drag
component that results from the change of the cross section of the projectile
is referred to as wave drag and it is attributed to the shock waves formed. This

1
2 A. A. Lee and F. B. Kong

Nomenclatures

CBT Boattail factor

Do Ratio of diameter of nose blunting to cylinder diameter (2ro/d)
d Body diameter, m
Kb Correction factor for base drag
LN Original nose length of pointed conical nose (Fig. 1), m
M Freestream Mach number
pb Base pressure coefficient for cylindrical boattail
Re Reynolds number
ro Radius of nose blunting (Fig. 1), m
Sref Reference area (d2/4), m2
xcp Centre-of-pressure location measured from the nose apex, m

Greek Symbols
 Angle of attack, deg.
 Mach number parameter, M 2  1
BT Ratio of base diameter to cylinder diameter (dBT/d)

 Cyl Ratio of cylindrical part length to nose length of pointed cone shape
N Fineness ratio of original pointed cone ( LN / d )
 Semi-vertex angle of the conical nose (Fig. 1), deg.

Abbreviations
ISA International Standard Atmosphere
NACA National Advisory Committee for Aeronautics
WHO World Health Organization
normally happens at the forebody (nose) and afterbody (tail). Since the wave may
be the prevailing drag form at supersonic speeds, careful selection of the nose and
tail shapes is mandatory to ensure performance and operation of the over-all system.
Shahbhang and Rao [2] conducted an experimental investigation to determine
aerodynamic characteristics of cone-cylinder and ogive-cylinder bodies of
different fineness ratios at Mach number of 1.8. Their results indicated that the
normal force for ogive-cylinder body is slightly higher than that for cone-cylinder
body of nose fineness ratio 3 and lower than that for cone-cylinder body for nose
fineness ratio 7 and there is crossing of normal force curves for nose fineness
ratio equal to 5. This interesting phenomenon requires further investigation.
The objectives of this paper is to predict the aerodynamic characteristics of
projectiles using analytical and semi-empirical methods and study the effect of
body shape; forebody and afterbody, on the aerodynamic characteristics of
projectiles at supersonic speeds. For this purpose five widely used projectile shapes
are investigated. The geometry and full dimensions of these projectile shapes are
shown in Fig. 1. The models are: (a) cone-cylinder, (b) ogive-cylinder, (c) blunted
cone-cylinder, (d) cone-cylinder boattail (4o), and (e) cone-cylinder boattail (8o). All
the models have a fineness ratio of 6.67 and a centre-of-gravity location at about the
40% body station. The supersonic Mach number range considered is from 1.6 to 5
for zero-angle of attack.

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

 This is the Template You Use to Format and Prepare Your Manuscript 3

(a) Pointed-cone cylinder. (d) Cone cylinder boat-tail (4o).

(b) Pointed-ogive cylinder. (e) Cone cylinder boat-tail (8o).

(c) Blunted-cone cylinder. (f) Blunted-cone geometry.

Fig. 1. Investigated shapes of projectiles (Geometry and dimensions).

2. Prediction of Aerodynamic Coefficients

Analytical methods and design charts used for the prediction of zero-lift drag
coefficient CD0, normal-force-curve slope CN, and centre-of-pressure location xcp
of body of revolution at supersonic speeds are presented in this section. The
analytical methods are based on supersonic linearised theory and thus they are
limited to slender bodies and low angles of attack, i.e., in the linear range of the
normal-force-curve slope while the design charts are produced from semi-
empirical results. The design charts are coupled with the analytical methods to
improve the accuracy of the results. These design charts used for the prediction of
aerodynamics characteristics are adapted from Refs. [3,4] and converted to
numerical data, as outlined in Appendix A.

2.1. Zero-lift drag coefficient CD0

The total zero-lift drag coefficient of the body is usually considered to be of three
components; friction drag, wave drag, and base drag as shown in Eq. (1). These
different components are further discussed in the following sub-sections.
C D0  C D fr  C Dw  C Db (1)

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

4 A. A. Lee and F. B. Kong

2.1.1. Friction drag coefficient

For fully-turbulent and compressible flow, the friction coefficient is given by Eq. (2)
[5, 6]

0.455 (log 10 Re) 2.58 S wet

C D fr  (2)
(1  0.21M 2 ) 0.467 S ref

2.1.2. Wave drag coefficient

The main contribution to the wave drag arises from nose and afterbody. The
magnitude of the wave drag depends primarily on the Mach number, the shape
and dimensions of the nose or afterbody. Therefore, the total wave drag of the
body is simply the summation of the nose and afterbody wave drags
C Dw  C Dw  N  C Dw  BT (3)

The wave drag of pointed cone-cylinder (CDw)cone and pointed ogive-cylinder

(CDw)ogive can be obtained from Figs. A-1 and A-2 (Appendix A) as a function of
nose fineness ratio N, and Mach number. For blunted cone-cylinder the wave
drag can be determined as a function of N, Mach number, and diameter of nose
bluntness D0 using Eq. (4) [4]

(C Dw ) blunted cone  C Dw cone (1  D02 cos 2  )  C Dw sphere D02 (4)

where (CDw)cone is the wave drag of the original pointed cone with N  0.5 / tan  ,
and (CDw)sphere is the wave drag of hemispherical nose, which can be determined
from Fig. A-3 as function of N and Mach number.
The wave drag of conical boattail (CDW)BT can be evaluated from Fig. A-4 as a
function of BT, BT and Mach number.

2.1.3. Base drag coefficient

At supersonic speeds the base drag of the body, caused by a large negative
pressure, results in a substantial increase in the body drag. The base drag
coefficient of the body is related to the base pressure coefficient as follows [4]

C Db   p b K b BT
2
(5)

where pb is the base pressure coefficient for cylindrical base (determined from Fig.
A-5 as a function of Mach number), and Kb is a correction factor, which depends on
the geometry of boattail Kb = f (CBT, M), obtained from Fig. A-6 [4], where
1   BT
C BT  (6)
2 BT  BT
2

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

 This is the Template You Use to Format and Prepare Your Manuscript 5

2.2. Normal-force-curve slope CN

The total normal-force-curve slope of nose-cylinder-boattail body is determined
by the summation of the normal-force-curve slopes of the nose (with the effect of
cylindrical part) and afterbody.
C N a  C N a N  C N a BT (7)

At supersonic speeds design charts are presented for estimating the normal-
force-curve slope of bodies of revolution composed of conical or ogival noses and
cylindrical afterbodies.
The normal-force-curve slope of blunted cone-cylinder can be evaluated from
[3] as a function of N, Mach number, and D0
C N N  C N cone (1  D02 )  C N sphere D02
  
(8)

where (CN)cone is the normal-force-curve slope of pointed cone with  / N and
.curve slope of hemispherical nose-force-the normalis sphere)NC, 
and (  Cyl

2.3. Location of the centre-of-pressure xcp

The identification of the location of the centre of pressure of a projectile body is
motivated by the need for calculating aerodynamic moments, stability and
structural analyses. The centre-of-pressure location of bodies composed of conical
noses and cylindrical afterbodies is determined as follows [5]
( xcp ) N  (0.667 CCyl ) LN (9)

3. Computer Programme: Validation and Verification

To ensure the validity and accuracy of the calculations, the results are compared
to available experimental wind tunnel data. Normal force coefficient and base
drag coefficient are compared as a function of Mach number and angle of attack.
Two typical projectile configurations (as shown in Figs. 2 and 3) are selected for
this purpose. The specifications of the models and test conditions are shown in
Table 1.
Table 1. Test model Specifications and test conditions.
Model No. 1 [2] Model No. 2 [7]
Configuration Type Cone-cylinder Ogive-cylinder
Body fineness ratio 13 12.2
Nose fineness ratio 3 7.5
Body Diameter, d (inches) 1 6
Reference length, Lref d d
Reference area, Sref d2/4 d2/4
Testing Mach number 1.77 1.5, 1.6, 1.79, 1.99
Testing angle of attack (deg.) 0–6 0

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

6 A. A. Lee and F. B. Kong

Figures 2 and 3 show the comparison between the current results and the wind
tunnel experimental data. Figure 2 shows that at low angles of attack the normal
force coefficients are in excellent agreement with the experimental data. This is to
be expected due to the assumption of small angle of attack. The figure also shows
that the current results are closer (average percentage error less than 0.5%) to the
experimental data than those predicted analytically (average percentage error
about 6%) by Shahbahang and Rao [2]. This is expected as the analytical methods
contained a number of simplifying assumptions that limit their accuracy and range
of use.

A comparison of base drag coefficient as a function of Mach number is shown

in Fig. 3. The comparison shows that at low supersonic Mach numbers the
average percentage error is 12%, while at higher Mach number the accuracy is
excellent (error less than 2%). At low supersonic Mach numbers the base drag
contribution is greater than the wave drag of nose, but with increasing Mach
number the wave drag contribution is the largest [8]. However the obtained
accuracy is still within the range of 10% error, which is considered sufficient to
be used in the preliminary design of projectiles and missiles [9].

Pointed-Cone Cylinder [2]. Pointed-Ogive Cylinder [7].

Theory (present work) 0.2

0.4
Theory [2]
Experiment
0.3
CDb

0.15
CN

0.2

Theory (present work)

0.1
Experiment [4]
0
0.1 Linear (Theory (present
0 1 2 3 4 5 6
1.4 1.6work)) 1.8 2 2.2
 (deg.) M

Fig. 2. Variation of normal force vs. Fig. 3. Variation of base drag

angle of attack at M = 1.77. coefficient vs. Mach number.

4. Results and Discussion

The prediction of the aerodynamic coefficients of the investigated projectiles
shown in Fig. 1 was carried using the methods and the computer programme
described above. The effects of forebody and afterbody shapes on the
aerodynamics at supersonic speeds are analysed in this paper.

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

 This is the Template You Use to Format and Prepare Your Manuscript 7

4.1. Effect of forebody

4.1.1. Zero-lift drag CD0
Figure 4(a) shows the effect of nose shape on CD0 with cylindrical afterbody as a
function of Mach number. The drag of cone-cylinder combination was the lowest
at the considered Mach numbers. It is clear that the bluntness of nose causes the
drag to increase.
4.1.2. Normal-force-curve-slope CNand location of centre-of-pressure xcp
The effect of forebody on the normal force curve slope and centre-of-pressure
location for the cylindrical afterbody is shown in Figs. 4(b) and (c) as a function
of Mach number. For all three shapes the variations of xcp are reasonably similar
and in general indicate the most rearward location with the ogival nose. The
centre-of-pressure locations are apparently a result of the normal force
distribution over the bodies with the blunted cone producing the least lift forward,
thus resulting in a more rearward xcp.

0.4 PCC
0.4
CCB4
0.3 0.3
CCB8
CD0
CD0

0.2 0.2

0.1 0.1

0 0
1 2 3 4 5 6 1 2 3 4 5 6
M M

(a) Zero-lift drag coefficient CD0. (a) Zero-lift drag coefficient CD0.
3.5 3.5

3
CN (1/rad)

3
CN (1/rad)

2.5 2.5

2 2
1 2 3 4 5 6 1 2 3 4 5 6
M M

(b) Normal-force-curve slope CN. (b) Normal-force-curve slope CN.

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

8 A. A. Lee and F. B. Kong

70
70

60
60

xcp /L (%)
xcp /L (%)

50 50

x cg x cg
40 40

30 30

20 20
1 2 3 4 5 6 1 2 3 4 5 6
M M

(c) Centre-of-pressure location xcp/L. (c) Centre-of-pressure location xcp/L.

Fig. 4. Effect of forebody shape Fig. 5. Effect of afterbody shape

on the aerodynamic characteristics. on the aerodynamic characteristics.
For conical, ogival, and blunted cone forebody shapes, an inherent static
stability occurs for a centre-of-gravity location of about 40% body length at Mach
number above around 1.6, 1.8 and 2 respectively. Such a centre-of-gravity
location may not be difficult to achieve with a projectile [10].

4.2. Effect of afterbody

4.2.1. Zero-lift drag CD0
For the projectile configuration comprising conical forebody and boattail, the effect
of boattail shape on the drag is shown in Fig. 5(a) as a function of Mach number.
For all the investigated configurations, the zero-lift drag coefficient, decreases as
Mach number increases. This is a typical behaviour for this curve for all bodies
flying at supersonic speeds. The high drag for the cone-cylinder combination was
primarily a result of higher base drag than the bodies with boattail. It is also seen
from this figure that the higher the angle of boattail the lower is the drag.

4.2.2. Normal-force-curve-slope CNand location of centre-of-pressure xcp

The effect of afterbody shape on the normal-force-curve slope and centre-of-
pressure location for conical nose are shown in Figs. 5(b) and (c) as a function of
Mach number. The variations with Mach number are reasonably similar with the
most forward location of xcp occurring with the boattail. The increasing of the
angle of boattail results in decreasing of the projectiles static stability.
Accordingly the projectile with cone-cylinder is statically stable above the
entire Mach number range, while the projectile with boattail (4o) is stable above
Mach number 1.8. The cone-cylinder boattail (8o) projectile is shown to be stable
only within the Mach number range 2 to 4.

5. Conclusions
An investigation has been made of the effects of forebody and afterbody shapes of a
series of projectiles on the aerodynamic characteristics at Mach numbers from 1.6 to
5. This is done using analytical methods combined with semi-empirical design
curves. Some concluding observations from the investigation are given below.

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

 This is the Template You Use to Format and Prepare Your Manuscript 9

 A pointed cone-cylinder produced the lowest drag at the considered Mach

number range, and the highest drag was produced by the blunted cone-cylinder.
 The shape of forebody slightly affects the normal force and centre-of-
pressure location. The farthest aft centre-of-pressure locations were
obtained with the ogive-cylinder and the most forward locations with a
boattailed afterbody.
 With the exception of the boattail afterbody, all the considered projectile
shapes indicated inherent static stability above a Mach number of about 2
with the centre-of-gravity location of about 40% body length.
 Configurations with boattail have higher wave drag but appreciably lower
base drag with a resultant decrease of total drag. The decrease of the boattail
angle increases the base drag but reduced the projectile wave drag with a
resultant decrease of the total drag.
References
1. Vukelich, S.R.; and Jenkins, J.E. (1982). Evaluation of component buildup
methods for missile aerodynamic prediction. Journal of Spacecraft and
Rocket, 19(6), 481-488.
2. Shahbhang, V.V.; and Rao, R.U. (1970). Normal force characteristics of cone-
cylinder & ogive-cylinder bodies at Mach number 1.8. Technical memorandum
No. TM-PR.235/69-70, National aeronautical laboratory, Bangalore.
3. Lebedev, A.A.; and Chernobrovkin, L.S. (1973). Dinamika poleta
bezpilotnykh letatelykh apparatov. "Flight dynamics of unmanned vehicles".
Mashinostroenie, Moscow.
4. Jankovic, S. (1979). Aerodynamika projektila. "Aerodynamics of projectiles"
Masinski Fakultet, Univerziteta u Beogradu, Belgrade.
5. Krasnov, N.F. (1985). Aerodynamics. Part 2 (Methods of Aerodynamic
Calculations). Translated from the Russian language by G. Leib, Mir
Publisher, Moscow.
6. Hoerner, S.F. (1965). Fluid-dynamic drag. Horner Fluid Dynamics, Brick
Town N.J.
7. Cohen, R.J. (1951). Aerodynamic characteristics of four bodies of
revolution showing some effects of afterbody shape and fineness ratio at
free-stream Mach numbers from 1.5 to 1.99. NACA research
memorandum, RM E51C06.
8. Al-Obaidi, A.S.M. (2005). The influence of Configuration Design on the
Aerodynamics and Stability of SSM at Supersonic Speeds. Ph.D. Thesis. Al-
Rasheed College of Engineering and Science, University of Technology,
Baghdad, Iraq.
9. Brebner, G.G. (1971). A brief review of air flight weapons. AGARD Lectures
No. 98 'Missile Aerodynamics', NATO. Aerodynamics Department, Royal
Aircraft Establishment. Farnborough, Hants GU14 6TD, UK.
10. Author, A.B. (2000). This is how to cite an internet reference. Retrieved
October 5, 2000, from http://www.author.com.
11. Roberson, J.A.; and Crowe, C.T. (1997). Engineering fluid mechanics (6th ed.).
New York: John Wiley and Sons Inc.

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

10 A. A. Lee and F. B. Kong

Appendix A
Representation and Figures of Design Charts
In the present work a number of empirical and semi-empirical design charts are
used for the prediction of the aerodynamic characteristics (Figs. A-1 and A-2).
These figures are adapted from the design charts given by Lebedev et al. [3] and
Jankovic [4]. The curves of those charts are read and converted to numerical data
and then stored in a separated subroutine in a computer programme described by
Al-Obaidi [8]. A simple linear interpolation is used to find the value of the
parameters used in the calculations.

0.42 0.24

0.35 0.18  N = 2.0

 N = 1.5
0.28 = 2.5
= 2.0
0.12

(CDW )N
= 3.0
(CDW )N

0.21 = 2.5
= 4.0
= 3.0 0.06
0.14 = 5.0
= 4.0
= 5.0 0
0.07

0 -0.06
0 1 2 3 4 5 6 0 1 2 3 4 5 6
M M

Fig. A-1. Wave drag coefficient of conical Fig. A-2. Wave drag coefficient of
noses with different fineness ratios [3]. ogival noses with different fineness
Ratios [3].

Appendix B
Computer Programme
B. 1. Introduction
A computer code, for the prediction of projectile aerodynamic characteristics as a
function of projectile geometry, Mach number and altitude of flight, is developed
in the present work. This programme is based on the analytical and semi-
empirical methods presented in section 2.
The computer programme can serve two main purposes: firstly, in the design
stage, a rapid parametric study of configuration can be performed to allow the
optimum configuration compatible with the requirements to be found and
secondly, by calculating the forces acting on a projectile at a range of speeds, the
programme is used in conjunction with both trajectory and stability calculations to
provide a complete picture of the projectile over its whole flight.

B. 3. Programme Structure and Description of Subroutines

Fortran-77 language is used in programming the prediction methods. Each
estimation method is programmed in a separate subroutine for case of
modification or addition to the programme. Many of design charts are used in this
programme. These charts are converted to numerical data (Appendix A) and stored
in a separate subroutine for convenience. The main flow chart of the programme
is shown in Fig. B-1.

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)

 This is the Template You Use to Format and Prepare Your Manuscript 11

START

Read Input Data

Geometry, Dimensions
and Flight Conditions

Calculate Projectile
Geometrical Parameters
Independent of Mach Number

For a certain Altitude Calculate

Air Properties (density,
acoustic speed, kinematic
viscosity)

Begin Mach Number Loop

I=1

Compute Zero-Lift Drag Coefficient

A of a Projectile
Sec. 2.1, Eqs. (1)-(6), Figs. A-1 to A-6

B Compute Normal-Force-Curve Slope of a

Projectile, Sec. 2.2, Eqs. (7)-(9), Figs. A-7 to A-10

Compute Centre-of-Pressure Location

of a Projectile
Sec. 2.3, Eqs. (10)-(14), Figs. A-11 to A-13

I=I+1
End of
No
Mach Number
Loop?
Yes
Print out Results

END

Fig. B-1. Main Flow Chart of the

Computer Programme Used in this Study.

Journal of Engineering Science and Technology Month Year, Vol. XX(Y)