Design of Optimized 3D Tip-To-Tail Scramjet Engines
Design of Optimized 3D Tip-To-Tail Scramjet Engines
Design of Optimized 3D Tip-To-Tail Scramjet Engines
Abstract
Worldwide, significant research efforts are in progress towards the development of scramjet
engines. These efforts are primarily focused on the integration of the scramjet to either ‘a blended body’
or ‘a waverider derived’ hypersonic vehicle configurations; one of which may become the candidate for
the first access to space air vehicle. The design of the scramjet is a complicated process due to many
technical constraints, such as, aerothermodynamics, materials and mechanical. Examples of aerodynamic
constraints include inlet-unstarts, mass capture, isolator characteristics, boundary layer separation, and
constraints on combustor entrance flow profiles and operating temperatures. Examples of mechanical
constraints include variable geometry flexibility and cooling system limits. All in all, the design of the
scramjet is influenced greatly by its flight constraints. If the flight conditions were fixed, then the scramjet
design problem became deterministic. As such, with a few assumptions, the engine configuration that will
produce minimum drag and maximum thrust may be determined directly. However, the design problem
under consideration for many practical applications requires extensive variability in altitude and Mach
number, and as such, its solution can only be obtained through an optimization process. Additionally, this
analysis assumed that the scramjet powered vehicle or missile is launched with an external propulsion
source to some ‘pre-determined’ Mach number and altitude. At that point, the scramjet engine is started
and propels the system until freestream flight conditions once again preclude further operation of the
engine. The research efforts described in this paper focused on the design of an optimized scramjet
configurations that are decouple for the hypersonic vehicle. Efforts are focused on the derivation of
efficient compression forebodies, isolators, combustion chambers, and nozzle after-bodies and their
integration into optimized scramjet configurations.
Nomenclature
= angle of attack
= shock wave angle
Cf = skin friction coefficient
Cp = pressure coefficient
D = Drag, force component parallel to the freestream velocity
= specific heats ratio
M = Mach number
= wedge angle
P = pressure
u = velocity component parallel to the freestream velocity
S = Surface area
T = temperature
OOP = Object Oriented Programming
1
Professor & Director, Center for Aerospace Research, and AIAA Senior Member.
2
Graduate Student, Center for Aerospace Research, and AIAA Member.
3
Graduate Student, Center for Aerospace Research, and AIAA Member.
Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.
I. Introduction
A. Background
Aircraft designers of the 21st century are focusing on a revolutionary engine technology that is
capable of not only propelling vehicles to hypersonic speeds, but also one that can facilitate integrated air-
to-space operations. The SCRAMJET, abbreviated from the words; Supersonic Combustion Ramjet, is
the latest evolution of the jet engine family. The scramjet, like its predecessor, the ramjet, is a natural
extension of the jet engine concept. However, unlike the ramjet, the scramjet uses no rotating parts.
Scramjets will enable three categories of hypersonic craft; namely, weapons, such as cruise missiles;
aircraft, such as those designed for global strike and reconnaissance missions; and space-access vehicles
that will take off and land like conventional airliners1.
While the physical concepts behind the scramjet are very simple, the practical ramifications of
constructing such an engine are quite formidable. A few of the challenges are supersonic fuel-air mixing,
the heat dissipation both from the air friction and the internal combustion, and the engine operating
temperatures. Consequently, the flow path of the incoming air needs to be extremely precise to minimize
hot spots. However, by far, the biggest challenges arise from the intense operational temperatures. Since
the air entering the engine is already heated by friction with the engine walls, combustion chamber
temperatures could exceed 5000 degrees Fahrenheit, if left unchecked1-2. At these temperatures most
metals melt, and air and fuel become ionized, and the physics of their behavior becomes unpredictable.
The objective of this paper is to identify the preliminary design variables of an idealized ram-scramjet
engine under optimum conditions.
A B C D E F G
Using the notations provided in Figure 3, the aerodynamic parameters, such as, density, pressure,
temperature and velocity can be computed as prescribed in Ref . 9. For instance, with a known cruise
Mach number, M and a prescribed shock angle, 1 , the initial shape of the forebody lower surface
angle, 1 , and the Mach number behind the shock, M1 , are determined by using the oblique shock
relations given in equations (1) and (2);
M 2 sin 2 1 1
tan = 2 cot (1)
M 2( + cos 2 1 ) + 2
1 1
and
M n,1
M1 = (2)
sin ( 1 )
Where the symbol, , is a constant equals to 1.4, and the symbols, M n,1 and M n, , are defined as
follows,
1 + [( 1) 2] M n2,
M n2,1 = (3)
M n2, ( 1) 2
and
M n, = M sin( 1) (4)
It is of interest to note that the value of 2 is dictated by the cowl location and a fixed point on the
forebody. When the results of 1 and M1 are coupled to an appropriate value for 2 , the flow properties
at the cowl inlet can be calculated using equations (1) and (2) in which the values of M 2 and 1 are
replaced with those of M1 and 2 .
[ ]
1
M 22 1 + (( 1) 2 )M 22
2
2
1
M3 = (5)
(1 M 22 P3 P2 )
2 2
Similarly, with the exit Mach number known, the length of the isolator can be evaluated based on the
following experimental relationship developed in Ref 11:
L
=
H {50((P
3 P2 ) + 1) + 170((P3 P2 ) 1)2 } (6)
(Re )
1
H 4 M 22 1
where Re is the inlet Reynolds number based on the momentum thickness, , and H is the height of the
isolator duct.
Now, consider a typical transition, combustion and nozzle combination of scramjet elements as
illustrated in Figure 9. A close look at Figure 9, indicates that the transition and combustor elements are
depicted with their external surfaces removed. This is done to highlight the geometry of the inner tube.
The aero thermodynamic evaluation of this part of the scramjet involves the quasi 1D calculations of the
flow field variables at each station in the length-wise direction of the duct illustrated in Figure 9. This task
is accomplished through the application of the appropriate conservation laws at the inlet and exit of each
element making up the scramjet, as described below:
( U A)inlet + m& f = ( U A)exit (7)
(u 2
A + pdA )
Inlet (
+ m& f + Fvis )Ductt = (u 2
A + pdA )
Exit (8)
$ V ! Vf $ V !
"( VA) h + + m& f h f + + Q& = "( VA) h + (9)
# 2 Inlet 2 # 2 Exit
The symbols, m& f hf and Q& , in equations (7) – (9) represent the rate of fuel added to the mass capture as it
travels through the combustor, hf, represents the enthalpy of the fuel and Q& represents the heat generated
in the combustor.
It is of interest to note that by choosing the design points, A – G, in the form of design parameters and
independent of each other, the resulting integrated configuration can have its E-point beyond its F-point
location with respect to the x-axis. This approach provides a great deal of flexibility in arriving at
scramjet candidates with desirable aerodynamic and thrust characteristics. In addition, when generating
the 3D configurations, the choice of the forebody is of great importance, since it provides the geometric
information needed for the third dimension. In this research, star-shaped, forebodies are generated with
the number of blades ranging from three to eight, with either circular or rectangular combustor c-sections,
as illustrated in Figures 1 and 10.
Figure 10: Illustration of ‘Four-Point Star-shaped Scramjet’ with Square Combustor C-Sections
III. Scramjet Configuration Design
A. Ramjet-Scramjet Analysis
The analysis of the scramjet system consists of the evaluation of the following six major
components; namely, the primary shock zone, the inlet, the isolator, the diffuser and mixing region
(transition element), the combustor and the nozzle. The evaluation of the primary shock zone, the inlet
and the nozzle is straight forward and is conducted in accordance with Ref. 3-4, 10. However, the
evaluation of the flowfield parameters in the transition, combustion and nozzle elements are somewhat
complicated and are conducted in accordance with the quasi 1D influence coefficients method described
in Ref. 14-16.
dp
pitotal
+1 = pitotal 1.0 + (13)
p i
Where the coefficients, C1 through C13, are evaluated at each spatial step in accordance Ref. 14-16.
(
where, the impulse function, I, is defined as, I = pA 1 + M 2 . , and the points, C and G, are described in )
Figure 1. In a similar manner the drag force, Fx, generated by the combined forebody, inlet and isolator
elements are computed. The performance parameters of this interest to this study are the net thrust, Tnet =
T – Fx, and the Thrust to Drag ratio, TD, where TD is defined as T/D.
where x = x1, x2, …, x10 is the vector of scramjet design variables, F(x) is the scalar cost function,
that represents either the thrust or the thrust to drag ratio, and G(x) is a vector-valued set of nonlinear
constraints, that represents the aerodynamic and geometric constraints, such as, the nozzle to inlet area
ratio, the isolator back pressure ratio and the combustion chamber maximum temperature. In general, the
function, F(x), is the fitness function characterizing how much the shape suits the optimization problem.
G(x) represents all the constraints of the problem. The design variables xi represent the geometric
parameters that control the general locations of points; A – G, that defines the 2D cross sectional
geometry illustrated in Figures 2 and 11. The goal of this optimization problem is to identify the ideal set
of geometric parameters that define the scramjet geometry, subsystems shape and aerodynamic
constraints for which the maximum thrust characteristic is observed.
Fig. 11: Scramjet 2D Cross Section with Optimization Design Parameters Highlighted.
In general, there are two families of optimization algorithms; namely a first family of ‘local
minimum solution’ based algorithms such as gradient methods, and a second family of ‘global minimum
solution’ based algorithms such as the genetic algorithm (GA), and simulated annealing (SA). In this
research effort, the optimization of the 3D scramjet configuration is conducted through the use of a
genetic algorithm. Starting from an initial scramjet configuration, the genetic algorithm is able to quickly
find an optimum. The method used is based on an OOP Fortran code.
Z ZoneD
ZoneC
ZoneB
ZoneA
Fig. 11: Isolator Shear stress Distribution Fig. 12: Isolator C-Section Pressure Distribution
Fig. 13: Isolator Internal Mach Number Fig. 14: Isolator Internal Pressure Distribution
Distribution
A.1: Preliminary Forebody-Isolator configuration Results
Figure 11 illustrates the shear stress distribution on the surface of the forebody-isolator element. In
Figure 11 the outer surface to the isolator is removed to reveal the shear stress distribution on the internal
surfaces of the isolator. Figure 12, on the other hand, shows the pressure distribution at cross-sections
located at stations, A, B, C and D, of the forebody-isolator element. These results indicated that the
pressures are constant at each location, however, the magnitude of the pressures increases as expected,
and are predicted as expected from the results illustrated in Ref. 3-4, 7-10. Attempts were made to
independently validate the forebody-isolator performance using CFD Euler and Navier –Stokes Solvers.
Preliminary results, indicating the Mach numbers and pressure behavior within the forebody-isolator
system from the Euler solutions for the 3D Mach 5 configurations, are illustrated in Figures 13 and 14.
In an effort to identify the design parameters necessary for optimization, the overall performance of
the forebody-isolator configuration was studied as a function of its geometric and freestream input
parameters. The captured mass flow rates, the Mach numbers at the exit of the isolator, the forebody-
isolator drag and the non-dimensional drag ratios were observed as the freestream mach number and the
wedge angle were allowed to vary. The Mach number was allowed to vary in the range of 2 to 12 while
the wedge was kept at 12 degrees. On the other hand, the wedge angle varied from 8 to 16 degrees, while
the freestream Mach number was fixed at 5. The results from these studies are illustrated in Figures 15
through 20.
Fig. 15: Isolator Mach Exit vs. Wedge Angle Fig. 16: Mach Exit vs. Flight Mach Number
Fig. 17: Isolator Mass Capture vs. Wedge Angle Fig. 18: Mass Capture vs. Flight Mach Number
Fig. 19: Isolator Drag Ratio vs. Wedge Angle Fig. 20: Drag Ratio vs. Flight Mach Number
B: Conclusion
This paper gives a preliminary report on efforts to design a reconfigurable propulsion system; morphing
from ramjet to scramjet mode of operation in flight. Performance requirements will almost certainly
dictate the use of a translating center body and associated outer ‘clamshell’. The design of the forebody-
inlet-isolator was accomplish through the use of 2D planar flowfields. In addition, realistic ramjet-to-
scramjet propulsion systems were derived and analyzed from Quasi 1D flowfields. Detailed analysis of
the ‘ramjet-to-scramjet propulsion systems’ performance at various Mach numbers is of interest to this
study and will form a significant element of this research project as it develops.
C. Conclusion
In general, the code developed as part of this research effort was used to conduct the following
studies:
1. Generate propulsion systems configurations from prescribed 2-D shock waves;
2. Evaluate the resulting engine geometric characteristics;
3. Evaluate the thrust performance of the engine, and;
4. Identify the design parameters that affect the engine’s overall performance and shape; and
develop the engine ‘optimization design space’.
At this stage of its development, the optimization process outline in this paper is yet to be tested and
debugged.
The preliminary outcomes of this research can be classified in the following two categories. First,
the propulsion system design and assembly process led to the discovery of engineering parameters that
directly influence the aerodynamic performance of the resulting configuration. These parameters were
manipulated to generate configurations with superior thrust and Isp characteristics. Second, design
information needed to develop routines for the follow-on morphing and optimization of the ‘ramjet-to-
scramjet configuration’ were generated and analyzed.
1500 1500
1.8 1.8
1.6 1.6
1400 1400
Local C-Sectional Area
Local Velocity
Local Pressure
1.2 1300 1.2 1300
1 1
0.6 0.6
1100 1100
0.4 0.4
3 4 5 6 3 4 5 6
x Location in Scramjet x - Location in Scramjet
Fig. 21: Scramjet Area & Velocity Distribution Fig. 22: Scramjet Area & Pressure Distribution
6000
Isp
4400
1.8
4200
5500
1.6
4000
Local Total Temperature
Local C-Sectional Area
1.4
5000
3800
1.2
3600
4500
1
3400
0.8
4000
3200
0.6
3000
3500
0.4
2800
3 4 5 6 4 6 8 10 12
x - Location in Scramjet Freestream Mach Number
Fig. 23: Scramjet Area & Total Temperature
Fig. 24: Scramjet Isp vs Mach Number
Distribution
VI: Acknowledgments
This work has been partially sponsored by the following agencies; WPAB, NAVAIR, NASA Glenn
and Langley Research Centers. In addition, special appreciation is extended to Dr. Datta Gaitonde and
Mr. Donnie Saunders of the Air Vehicles Directorate at Wright Patterson Air Force Base for their
encouragement and support of this project. Appreciation is expressed to Dr. Isaiah Blankson of NASA
Glenn Research Center and Dr. Reginal Williams of the Naval Airforce Base at Patuxent River in
Maryland.
VI: References
1
I. M. Blankson et al., “GLENN 20/20”. Results of the GRC Visioning Team. A presentation by the
Visioning Team to the DLT, OH, Nov. 2004.
2
I. M. . Blankson et al. “Anatomy of a Hypersonic Air-breather: A Turbine Based Combined - Cycle
Hypersonic Cruise Application-Vehicle.” International Symposium on Air Breathing Engines, ISABE
Paper 1999, Florence, Italy. Sept 1999.
3
Hydar Apdin and Frederick Ferguson, ‘A Design Concept for the Construction of Completed
Hypersonic Vehicles’, 13th AIAA/CIRA International Space Planes and Hypersonic Systems and
Technologies Conference, CIRCA, Italy, May 2005
4
J. Zhang and F. Ferguson, ‘Construction and Analysis of Hypersonic Vehicle Configurations’, 13th
AIAA/CIRA International Space Planes and Hypersonic Systems and Technologies Conference, CIRCA,
Italy, May 2005.
5
Bowcutt, Kevin G., “Optimization of hypersonic Waveriders Derived from Cone Flows Including
Viscous Effects”, Ph.D. Dissertation, Department of Aerospace Engineering, University of Maryland,
College Park, Maryland, 1986
6
Bowcutt, Kevin G., Anderson, John D. Jr., and Capriotti, D.P., “Viscous Optimized Hypersonic
Waveriders,” AIAA paper 87-0272, 1987.
7
Ferguson, F., “Expanding the Waverider Design Space Using Arbitrary Generating Flowfields”,
Ph.D. Dissertation, Department of Aerospace Engineering, University of Maryland, College Park,
Maryland, 1993.
8
Frederick Ferguson, Terry L. Corbett, Jr., Stephen Akwaboa, and Haile Lindsay, “The Development
of Waveriders From an Axisymmetric Flowfield”, AIAA 2007-847, 45th American Institute of
Aeronautics and Astronautics Aerospace Sciences Meeting and Exhibit, Reno, NV 2007.
9
John D. Anderson “Fundamentals of Aerodynamics” Third Edition, McGraw-Hill, New York, 2001.
10
William H. Heiser, David T. Pratt “Hypersonic Airbreathing Propulsion.” AIAA Education Series.
ISBN: 1-56347-035-7, Published by American Institute of Aeronautics and Astronautics, 1994.
11
E. T. Curran and S.N.B. Murthy, “Scramjet Propulsion”, Vol. 189, AIAA, 2000.
12
Waltrup, P. J., and Billig, F. S., ‘Prediction of Recompression Wall Pressure Distributions in
Scramjet Engines, “Journal of Spacecraft and Rockets, Vol., 10, No. 9, 1973, pp. 620-622.
13
Billig F. S.; “Research on Supersonic Combustion” Journal of Propulsion and Power, Vol. 9, No 4,
July-Aug 1993
14
Shapiro, H. Asher, ‘The Dynamics and Thermodynamics of Compressible Flow, Volume I, The
Ronald Press Company, 1953.
15
Markell, Kyle C., ‘Exergy Methods for the Generic Analysis and Optimization of Hypersonic vehicle
Concepts”, MS Thesis, Mechanical Engineering Department, Virginia Polytechnic Institute and State
University, 2005.
16
Jing, C., Shuo, T., and Bingnan, K., “Airframe/Scramjet Integrated Design of Hypersonic Cruise
Vehicle”, 45th AIAA Aerospace Sciences Meeting and Exhibit, AIAA Paper 2007-642, Reno, NV, 2007.