Report Final PDF

School of Innovation, Design and Engineering
BACHELOR THESIS IN
AERONAUTICAL ENGINEERING
15 CREDITS, BASIC LEVEL 300
Conceptual Design of a
Supersonic Jet Engine
by
Joakim Kareliusson
Melker Nordqvist
Vsters, September 2014

Advisor at MDH: Konstantinos G. Kyprianidis
Examiner: Mirko Senkovski
ii
Abstract
This thesis is a response to the request for proposal issued by a joint collaboration between the AIAA
Foundation and ASME/IGTI as a student competition to design a new turbofan engine intended for a
conceptual supersonic business jet expected to enter service in 2025. Due to the increasing
competition in the aircraft industry and the more stringent environmental legislations the new
engine is expected to provide a lower fuel burn than the current engine intended for the aircraft to
increase the chances of commercial success.
The thesis has one main and one secondary objective, the main objective is to perform a preliminary
design of a jet engine, complying with the specifications stated in the request for proposal. The
secondary objective is to evaluate whether the knowledge gathered from the BSc program in
aeronautical engineering at Mlardalen University is sufficient to perform this kind of study and to
recommend courses for future students who wish to perform a similar project. Due to time
restrictions a full conceptual design has not been performed and the work has mainly been focused
on the thermodynamic and aerodynamic design phases.
The thermodynamic analysis and optimization has been carried out using the Numerical Propulsion
System Simulation (NPSS) code, where the cycle parameters such as fan pressure ratio (FPR), overall
pressure ratio (OPR), turbine inlet temperature (TIT) and bypass ratio (BPR) have been optimized for
cycle overall efficiency. With the cycle selected, and the fluid properties at the different flow stations
known, the component aerodynamic design, sizing and efficiency calculations were performed using
MATLAB, with equations and formulas found in the open literature together with literature
developed at Chalmers University of Technology in Gothenburg, Sweden. Furthermore a number of
research papers have been used for various parts of the thesis. Several aspects of the
turbomachinery components have been evaluated to assure satisfactory performance. The result is a
two spool low bypass axial flow engine of similar dimensions as the current engine but with
increased efficiencies. A weighted fuel flow comparison of the two engines at the key operating
conditions shows a fuel burn improvement of 11,8% for the new engine.
The conclusions drawn by the authors on the secondary objective is that even though the knowledge
might be slightly insufficient to undertake this type of task, with proper guidance and determination
its certainly not impossible.
iii
Sammanfattning
Detta examensarbete r ett bidrag till en studenttvling skapad genom ett samarbete mellan the
AIAA Foundation och ASME/IGTI dr en ny motor avsedd fr en business jet med
verljudsegenskaper som frvntas tas i bruk r 2025 efterfrgas. P grund av den kande
konkurrensen inom flygindustrin och de strngare miljkraven frvntas den nya motorn frbrnna
mindre brnsle n den nuvarande motorn avsedd fr flygplanet och drmed ka chanserna p
marknaden.
Arbetet har tv syften, det frsta r att utfra en preliminr design av en jetmotor som uppfyller de
krav som finns i specifikationen. Det andra r att utvrdera huruvida frkunskaperna frn
flygingenjrsprogrammet p Mlardalens hgskola r tillrckliga fr att utfra den hr typen av
arbete och att rekommendera kurser fr framtida studenter som vill utfra ett liknande projekt. P
grund av tidsbegrnsningar har en full preliminr design inte utfrts, utan arbetet har huvudsakligen
fokuserats p den termodynamiska och aerodynamiska designen.
Den termodynamiska analysen och optimeringen har utfrts med hjlp av programmet the
Numerical Propulsion System Simulation (NPSS) dr vergripande cykelparametrar s som fan
pressure ratio (FPR), overall pressure ratio (OPR), turbine inlet temperature (TIT) och bypass ratio
(BPR) har optimerats fr hga verkningsgrader i motorn. D den termodynamiska analysen och
optimeringen frdigstllts och egenskaperna fr gasen vid de olika stationerna i motorn var knda,
berknades komponenternas aerodynamiska design, storlekar och effektiviteter i MATLAB med hjlp
av formler och ekvationer frn ppen litteratur tillsammans med material som tagits fram p
Chalmers tekniska hgskola i Gteborg. I vrigt har ett antal forskningsuppsatser anvnts till olika
delar av arbetet. Flera aspekter har utvrderats fr kompressorer och turbiner fr att skerstlla
fullgod prestanda. Resultatet r en tvspolig lg-bypass axialmotor med liknande dimensioner som
den nuvarande motorn men med kade effektiviteter. En jmfrelse av brnslefrbrukning ver ett
typiskt uppdrag visar en brnslebesparing p 11,8% fr den nya motorn.
P arbetets andra syfte r det frfattarnas sikt att ven om frkunskaperna r aningen otillrckliga
fr den hr typen av projekt, s r det med god handledning och en gedigen arbetsinsats definitivt
ingen omjlig uppgift.
iv
Acknowledgements
We would like to thank Konstantinos Kyprianidis for proposing this thesis, guiding us through the
process and always responding quickly to questions.
We also want to thank Mirko Senkovski for considering us when the thesis proposal came, as well as
the other teachers involved with the aeronautical program.
v
Contents
Abstract ...................................................................................................................................... iii
Sammanfattning...........................................................................................................................iv
Acknowledgements ...................................................................................................................... v
List of Figures ............................................................................................................................. viii
List of Tables ................................................................................................................................ix
Nomenclature .............................................................................................................................. x
1 Introduction ........................................................................................................................ 12
1.1 Problem area ............................................................................................................................. 12
1.2 Aim and objectives .................................................................................................................... 12
1.3 Work scope and limitations ....................................................................................................... 13
1.4 Thesis outline ............................................................................................................................. 13
2 Literature review ................................................................................................................. 14
2.1 Reducing fuel consumption ....................................................................................................... 14
2.2 Computer aided design.............................................................................................................. 15
2.3 Cost of development ................................................................................................................. 15
3 Methodology....................................................................................................................... 16
3.1 The preliminary engine design process ..................................................................................... 16
3.2 Thermodynamic analysis ........................................................................................................... 17
3.2.1 Design point and off design performance calculations ........................................................ 17
3.2.2 Software description ............................................................................................................ 17
3.2.3 Cycle modeling and optimization ......................................................................................... 17
3.3 Aerodynamic design .................................................................................................................. 21
3.3.2 Turbomachinery ................................................................................................................... 21
3.3.3 Combustion chamber ........................................................................................................... 24
3.3.4 Internal ducts........................................................................................................................ 25
3.3.5 Inlet....................................................................................................................................... 25
3.3.6 Nozzle ................................................................................................................................... 26
3.3.7 Mixer..................................................................................................................................... 26
3.3.8 Iterations .............................................................................................................................. 27
3.4 Mechanical design ..................................................................................................................... 28
3.4.2 Discs ...................................................................................................................................... 28
vi
3.4.3 Shafts .................................................................................................................................... 28
4 Results and discussion ......................................................................................................... 29
4.1 Baseline engine .......................................................................................................................... 29
4.2 Cycle optimization ..................................................................................................................... 30
4.2.1 TIT and OPR optimization ..................................................................................................... 30
4.3 Component design..................................................................................................................... 34
4.3.1 Inlet....................................................................................................................................... 34
4.3.2 Internal ducts........................................................................................................................ 34
4.3.3 Compressors ......................................................................................................................... 35
4.3.4 Combustion chamber ........................................................................................................... 38
4.3.5 Turbines ................................................................................................................................ 39
4.3.6 Mixer..................................................................................................................................... 42
4.3.7 Nozzle ................................................................................................................................... 43
4.3.8 Discs ...................................................................................................................................... 44
4.4 Off-design .................................................................................................................................. 45
4.5 Sensitivity analysis ..................................................................................................................... 46
4.6 Comparison to baseline engine ................................................................................................. 47
5 Conclusions ......................................................................................................................... 48
5.1 Future work ............................................................................................................................... 48
6 References .......................................................................................................................... 49
7 Appendices ......................................................................................................................... 52
A Formulas ....................................................................................................................................... 52
B NPSS outputs ................................................................................................................................ 55
C MATLAB program output files ...................................................................................................... 58
D T-AXI DISK output files ................................................................................................................. 69
vii
List of Figures
Figure 1: Preliminary engine design process [10] ................................................................................. 16
Figure 2: Different cycle efficiencies over velocity ratio in the mixer ................................................... 19
Figure 3: over OPR for different types of cycles [11] .................................................................... 20
Figure 4: Velocity triangles for a typical compressor stage [19] ........................................................... 22
Figure 5: Different types of flow paths.................................................................................................. 23
Figure 6: Typical burner [22] ................................................................................................................. 25
Figure 7: Theoretical thrust gain due to mixing [11] ............................................................................. 26
Figure 8: Different types of disc designs [28] ........................................................................................ 28
Figure 9: Section view baseline engine ................................................................................................. 29
Figure 10: Engine layout with flow stations .......................................................................................... 30
Figure 11: Optimum OPR at different TIT: s. ......................................................................................... 31
Figure 12: Efficiencies over OPR ............................................................................................................ 32
Figure 13: The MJ Haran S14 section view ......................................................................................... 33
Figure 14: The MJ Haran S14 section view with supersonic inlet ...................................................... 33
Figure 15: Supersonic inlet with flow stations ...................................................................................... 34
Figure 16: Ducts ..................................................................................................................................... 35
Figure 17: Stage load over flow coefficient for the FAN and HPC ......................................................... 35
Figure 18: Fan component .................................................................................................................... 35
Figure 19: Velocity triangles for the fan at mean blade ........................................................................ 36
Figure 20: HPC component.................................................................................................................... 37
Figure 21: Velocity triangles for the HPC at mean blade ...................................................................... 38
Figure 22: Combustor component ........................................................................................................ 38
Figure 23: Stage load over flow coefficient for the HPT and LPT .......................................................... 39
Figure 24: HPT component .................................................................................................................... 39
Figure 25: Velocity triangles for the HPT at mean blade....................................................................... 40
Figure 26: LPT component..................................................................................................................... 41
Figure 27: Velocity triangles for the LPT at mean blade ....................................................................... 42
Figure 28: Mixer component ................................................................................................................. 42
Figure 29: Nozzle component................................................................................................................ 43
Figure 30: Rate of change in area over distance and area over distance for the nozzle ...................... 43
Figure 31: First stage HPC disc output file from T-AXI DISC .................................................................. 44
Figure 32: Typical mission for the engine.............................................................................................. 45
viii
List of Tables
Table 1: Key data for the baseline engine ............................................................................................. 29
Table 2: Required HPT cooling flow and LPT inlet temperature ........................................................... 31
Table 3: Key parameters for design point cycle .................................................................................... 32
Table 4: Inlet data.................................................................................................................................. 34
Table 5: Duct data ................................................................................................................................. 34
Table 6: Key fan data ............................................................................................................................. 36
Table 7: Fan velocity triangles data at mean blade ............................................................................... 36
Table 8: Key HPC data............................................................................................................................ 37
Table 9: HPC velocity triangles data at mean blade .............................................................................. 38
Table 10: Key combustor data ............................................................................................................... 39
Table 11: Key HPT data .......................................................................................................................... 40
Table 12: HPT velocity triangles data at mean blade ............................................................................ 40
Table 13: Key LPT data........................................................................................................................... 41
Table 14: LPT velocity triangles data at mean blade ............................................................................. 42
Table 15: Key mixer data ....................................................................................................................... 42
Table 16: Key Nozzle data ..................................................................................................................... 43
Table 17: Key data for HPC and HPT discs ............................................................................................. 44
Table 18: M 1.5 alt 51000ft ................................................................................................................... 45
Table 19: M 0 alt 0ft .............................................................................................................................. 45
Table 20: M 0,98 alt 38500ft ................................................................................................................. 45
Table 21: M 1,15 alt 42000ft ................................................................................................................. 45
Table 22: Design point exchange rates for efficiencies and pressure losses ........................................ 46
Table 23: Design point exchange rates for power off take, HPT metal temp and customer bleed ...... 46
Table 24: Comparison between baseline engine and MJ - Haran S14 .................................................. 47
Table 25: Primary contributors of increased efficiencies ...................................................................... 47
ix
Nomenclature
Parameter Symbol Units
Absolute Velocity /
Angularity Coefficient
Aspect Ratio
Axial Mach number
Axial Velocity /
Blade Stress Level 2 2 2
Blade Velocity /
Bypass Ratio
Chord
Convergent Divergent
Corrected mass flow /
Core Efficiency
Entry into Service
Fan Pressure Ratio
Flow Coefficient
Fuel to Air Ratio
Gross Thrust
High Pressure Compressor
High Pressure
High Pressure Turbine
Inlet Guide Vane
Intercompressor Duct
Interturbine Duct
Isentropic Efficiency
Loss Coefficient for Stator Blades
Lower Heating Value /
Low Pressure
x
Low Pressure Turbine
Mach number
Mass flow /
Momentum Drag
Net Thrust
Nozzle Guide Vane
Overall Efficiency
Overall Pressure Ratio
Pitch s
Polytropic Efficiency
Pressure Loss Factor
Pressure Recovery
Propulsive Efficiency
Relative Blade Mach number
Relative Velocity /
Revolutions per Minutes
Revolutions per Seconds
Specific Fuel Consumption
Specific Thrust /
Stage Load
Static Pressure
Static Temperature
Thermal Efficiency
Thrust Coefficient
Total Pressure 0
Total Temperature 0
Transfer Efficiency
Turbine Inlet Temperature
Water to Air Ratio
xi
1 Introduction
1.1 Problem area
The request for proposal is a joint AIAA Foundation and ASME/IGTI student design competition
where a new turbofan engine for a conceptual supersonic business jet aircraft expected to enter
service in 2025 is requested. It will be able to travel from North America to Europe and back again
within one business day. It can cruise at Mach 1.15 over land without creating a sonic boom on the
ground. Over water, it can cruise at Mach 1.5. Transonic cruise at Mach 0.98 should offer similar cost
per distance to subsonic private jets. [1]
Todays aircraft industry is very competitive. The commercial success of an aircraft highly depends on
its ability to provide light weight, low noise, competitive fuel consumption, and in the end lower
operating costs and passenger fares. A contributing factor to achieve this is the engine design.
Furthermore, due to the increase of the environmental concerns over the impact of the growing civil
aviation, the emission legislations set by the International Civil Aviation Organization (ICAO) are
becoming more stringent. The airlines need to continuously reduce their operating costs, which
combined with stringent emission legislations introduces new challenges for the aero-engine
industry: new engines need to be developed for low operating costs as well as reduced
environmental impact. A reduction of fuel consumption, which is directly proportional to CO2
emissions, can be achieved primarily by increasing component efficiencies, increasing OPR and TIT
reducing SFN and engine size and weight. A number of novel engine concepts addressing these issues
have been proposed over the years, some of which will be presented briefly in chapter two of this
thesis. Implementations of such concepts however, involve a great financial risk for the
manufacturers [2]. In this thesis, a more conventional cycle is considered.
Engine development is a multidisciplinary process covering thermodynamics, aerodynamics and solid

mechanics, where the requirements of one discipline may very well contradict the requirements of
another, and trade-offs between performance, size, weight, and costs continuously need to be made.
Engine specifications are provided in the RFP, including thrust requirements, maximum dimensions
etc. A generic model of a baseline engine is also provided. The new engine is expected to have
improved fuel consumption and overall performance, thus increasing the chances of commercial
success for the aircraft. The new engine has been given the designation MJ Haran S14.
1.2 Aim and objectives

There is one primary objective and one secondary objective for this thesis. The primary is to carry out
a preliminary design for an engine complying with the requirements specified in the RFP. This can be
broken down in the following steps:
1. Model the baseline engine and perform a thermodynamic study of the cycle performance.
2. Optimize the engine cycle for low fuel consumption, and high cycle overall efficiency.
3. Carry out preliminary sizing, efficiency calculations and aerodynamic design for the main
engine components.
4. Create a schematic design of the engine including dimensions.
12
The secondary objective is to evaluate if the knowledge acquired from the courses at the
aeronautical program at MDH is sufficient to carry out a preliminary engine design process, and to
recommend courses for future students who wish to perform a similar task.
1.3 Work scope and limitations

Due to time restrictions a full preliminary design has not been carried out in more detail. The
mechanical design is limited to the disc design of the HPC and HPT. No other material choices have
been made and hence no overall weight estimations.
1.4 Thesis outline

Chapter 2 contains a literature review where published literature regarding the problem areas of the
thesis is discussed. It provides an overview of relevant past and current research on the topic.
Chapter 3 provides the methodology used to achieve the results of the thesis, it motivates why and
how certain parameters are chosen. It also contains a description of the software used.
In chapter 4 the results of the simulations and choice of parameters is presented and discussed.
Some comparisons to the baseline engine are also provided.
Chapter 5 concludes the thesis, looking back and evaluating the objectives. Recommendations for
future work are also presented.
13
2 Literature review
A large number of research papers have been published regarding the aero engine industry over the
years. In this chapter, a few of these regarding important aspects of the engine design process will be
discussed.
2.1 Reducing fuel consumption

Today one major concern in the industry is producing engines with better fuel economy, partially to
reduce the operating costs for the airlines leading to lower fares for paying passengers, but also to
address todays environmental concerns. A few approaches to achieve this has been proposed.
A summary of some of the current research being done in the search for more fuel efficient engines
is provided in [2]. The author starts by discussing the limitations of increasing OPR and TIT in order to
increase thermal efficiency and reduce SFC, which has been the trend for conventional cores over the
past decades. Increasing OPR further than current designs is limited by HPC delivery temperature at
take-off. Increasing TIT is limited by the maximum HPT rotor temperatures at take-off and top of
climb. The increase of cooling flows for this reason is also a limited strategy since this represent
losses in the cycle and will eventually lead to a loss in thermal efficiency.
In [3] the authors investigate the benefits of introducing Intercooled and intercooled recuperated
cores. An intercooled core can be designed for a significantly higher overall pressure ratio, with
reduced cooling requirements, providing a higher thermal efficiency than could be practically
achieved with a conventional core. In an intercooled recuperated core a high thermal efficiency can
be achieved at a low overall pressure ratio. The use of a variable geometry auxiliary nozzle for the
intercooled engine, and a variable geometry LPT for the intercooled recuperated engine is also
evaluated, showing further benefits. The use of recuperation is also investigated in [4], where a
recuperated engine for a UAV is compared to a conventional design. The novel engine is expected to
maintain the same weight as the conventional engine and extend the mission endurance by 15%.
In [5], a geared open rotor configuration engine is compared to an ultra-high bypass geared turbofan
engine, with focus on minimizing fuel consumption and hence engine emissions. The open rotor
concept is not a new idea, it was evaluated in the 80s in several flight tests, but most of them never
got past that phase and lately the concept has resurfaced. In an open rotor configuration the
propulsive efficiency may be increased without increasing nacelle drag, and losses due to transfer
efficiency are comparable to those of a turbofan. The results show that the open rotor engine is
heavier but the reduced SFC and nacelle drag makes up for this and a mission fuel burn improvement
of approximately 15% compared to the turbofan is presented.
In [2] the author concludes by highlighting the fact that even though research of more efficient
engine concepts exists, the aero engine design industry is primarily driven by economic
considerations. Introducing these novel concepts to the market involves an economical and
technological risk, and whether the potential reduction in fuel consumption outweighs these risks
remains to be seen.
14
2.2 Computer aided design
Over the years, computers have become an invaluable tool for engineering purposes. In [6] the
authors discuss the impact computers have had on the industry during the 1980s. Since then the use
of computers have increased exponentially and software is constantly being developed and updated.
The authors of [7] discuss the development of a multidisciplinary optimization code at the time being
developed at NASA Lewis research center for the design of structural components, subsonic and
supersonic aircraft configuration design and airbreathing propulsion engine design. The design tool
can optimize a system which can be defined in terms of fifty optimization subproblems. The system
first formulates the design as a nonlinear mathematical problem, and solves the resulting problem
based on data from specified input files.
In [8] another piece of software is described, meant to be used in the preliminary design phase of gas
turbines. The software was created to allow for a quick preliminary design process allowing the
engine suppliers the ability to evaluate numerous concepts to fulfill the market requirements. The
system is built upon three major features:
i. All major engine components and their interrelations are assessed.
ii. Several relevant disciplines are considered.
iii. Designing is done on several operating points and the off design characteristics are
evaluated.
2.3 Cost of development

In [9] the relationship between development and cost is discussed. As mentioned earlier the industry
is driven by economic consideration and the cost of achieving a performance target must be taken
into consideration. The paper evaluates the tools used at Rolls Royce for the preliminary design
phase of military engines and emphasizes a paradox in the preliminary design stage; there is little
product knowledge but there is a powerful impact on the final design!
15
3 Methodology
3.1 The preliminary engine design process
The engine design process usually starts with an engine specification. The specification can either be
a result of the requirements presented by a customer, or by the manufacturer itself to meet a market
need.
The first step in the design process is the thermodynamic cycle analysis. The configuration of the
engine, the cycle parameters and the performance of the components are selected and evaluated to
meet the given specification.
When the major parameters are set by the thermodynamic analysis, the aerodynamic design of the
turbomachinery and other components can begin. The number of stages, overall sizes, rotational
speeds, efficiencies and other aerodynamic parameters can be determined.
Once the aerodynamic design is complete, the mechanical design of the components can begin, in
this stage the mechanical properties such as stress and vibrations are evaluated.
Figure 1: Preliminary engine design process [10]
16
This process needs constant feedback between the different disciplines, and changes in one may lead
to changes in another. Figure 1 shows a schematic diagram representing a general design procedure.
This thesis will cover the thermodynamic, aerodynamic and a small part of the mechanical design.
The full specification of the engine requirements can be found in the RFP [1].
3.2 Thermodynamic analysis

3.2.1 Design point and off design performance calculations
At the design point, the thermodynamic performance of the engine is evaluated at a fixed operating
condition. In the case of the MJ Haran S14 the operating condition chosen is supersonic cruise at M
1.5 at 51000ft, where the engine is expected to spend most of its time. Changing the operating
condition in the design point stage would result in a different engine with a different geometry. This
conditions used as a reference point.
The next step is to run off design performance calculations. The engine geometry is now fixed and
the performance of the engine at different operating conditions is evaluated. Many off design cases
may be run to ensure the engine performance over the entire operating range. Both design point and
off design performance calculations are highly iterative processes [11].
3.2.2 Software description

For the thermodynamic analysis of the engine, the Numerical Propulsion System Simulation (NPSS)
code has been used. NPSS was developed by NASA in cooperation with U.S aerospace industries with
the aim to develop a state of the art simulation tool for complex systems with an open architecture
for user flexibility. The code is written in C++ and allows users to add their own unique objects and
calculations. The program has a built in solver with the ability to make one parameter dependent on
another, and hence varying the user input of these parameters to satisfy a certain condition specified
by the user [12].
3.2.3 Cycle modeling and optimization

A replica of the baseline engine was modelled as a starting point for the cycle optimization. The cycle
was then modified to achieve the required thrust levels at design point. The key cycle performance
data can be seen in chapter 4. All comparisons herein are to the modified baseline engine. For the
initial cycle optimization, turbomachinery efficiencies and pressure losses were assumed to be the
same as for the baseline engine.
3.2.3.1 Cycle efficiency

A number of different types of efficiencies can be used to characterize the performance of an aero-
engine:
Propulsive efficiency the ratio of useful propulsive energy. This is the product of thrust and flight
velocity divided by the wasted kinetic energy of the jet.
Thermal efficiency the increase of the kinetic energy of the gas stream divided by the energy of the
fuel, which is the product of the fuel mass flow and the fuel LHV. For a turbofan engine, because of
17
the two gas streams the thermal efficiency is divided into two terms, the core efficiency and the
transmission efficiency.
Core efficiency the energy available after all the power requirements of the core stream are
satisfied, which is the energy available at the core exit, divided by the energy of the fuel. This is
evaluated assuming an isentropic expansion from the state at the core exit to ambient pressure.
Transfer efficiency the quality of the energy transfer from the core stream to the bypass stream.
This is the energy available at the nozzle divided by the energy at the core exit.
Overall efficiency the ratio of useful work done in overcoming the drag of the airplane, and the
energy of the fuel. For a turbofan engine, this is the product of the core efficiency, the transfer
efficiency and the propulsive efficiency [13].
SFC The mass flow of burnt fuel per unit of time per unit of output thrust [11]
Detailed formulas are provided in appendix A.
3.2.3.2 Selection of cycle parameters

The primary parameters for optimizing the cycle of a turbofan engine are the following: FPR, BPR,
OPR and TIT. The definitions can be found in appendix A. These are limited by a number of
mechanical design constraints such as creep, oxidation, casing rupture and vibrations. Herein the
following secondary design parameters have been evaluated: Cooling air, pressure losses,
turbomachinery efficiencies, power off-take, mixer efficiency and exhaust efficiency. It should be
noted that the engine mass flow was kept constant in this study in order to maintain a constant
engine diameter.
Fan pressure ratio

The fan pressure ratio was selected by letting NPSS solver vary the FPR to always achieve the
required specific thrust. The required FPR to satisfy this condition is highly dependent on the
remaining parameters in the engine, and will vary as other parameters vary. For the MJ Haran S14,
the core stream FPR and the by-pass stream FPR is assumed to be the same, which in reality may be
hard to achieve due to the different blade speeds at the fan tip and the fan hub.
Bypass ratio
It can be shown that for every FPR there is an optimum BPR and vice versa. At this condition the
overall energy conversion is maximized giving minimum SFC, maximum specific thrust and maximum
engine overall efficiency. In [14] it is shown that for a mixed flow turbofan this occurs when the

velocity ratio between the two streams ~0,8 . The BPR has been selected using this

relationship, letting the NPSS solver vary the BPR keeping the velocity ratio fixed at any FPR and
overall condition.
, , and SFC over the velocity ratio are illustrated in figure 2. SFN, TIT and OPR are kept
constant. It can be seen that overall efficiency is at its maximum and SFC at its minimum when the
velocity ratio is approximately 0,8.
18
Figure 2: Different cycle efficiencies over velocity ratio in the mixer
Turbine inlet temperature

In principle, increasing the TIT and OPR will increase thermal efficiency and hence SFC. However, due
to material temperature limitations, a higher TIT will need higher cooling flows for the turbines.
Cooling flows essentially represent losses in the cycle and at a certain point these losses will
overcome the gain in efficiency [2]. The MJ Haran S14 uses an uncooled LPT to keep cooling flow
losses to a minimum and production costs down, thus limiting the TIT.
Overall pressure ratio

For a turbofan engine, there is an optimal OPR for a given TIT where the core and thermal efficiencies
are at their maximum and hence the SFC at its minimum [15]. This was found by creating diagrams of
core efficiency vs OPR, using different temperatures and constant SFN, letting the NPSS solver
optimize the remaining parameters. A generic diagram showing this relationship for different types
of cycles at a specific TIT can be seen in figure 3.
19
Figure 3: over OPR for different types of cycles [11]
Cooling flows
The necessary cooling flows have been calculated using the method described in [16]. The required
cooling flow will depend on the cooling flow efficiency, which varies depending on the technology
being used. The calculations roughly correspond to a film/cross flow impingement technology. Two
cooling flows have been considered for the HPT: one for the NGV and one for the rotor. An upper
limit of 18% of the HPC inlet flow was decided to avoid excessive losses. The cooling flows were
added as a variable to the NPSS solver to keep constant blade metal temperatures based on the
blade temperatures for modern designs under all overall conditions. The flow is taken from the last
stage of the HPC to assure sufficient pressure.
Power off-take
The power off-take from the HPT has been set to 100hp in accordance to the RFP [1].
20
3.3 Aerodynamic design
The aerodynamic design can commence when the flow properties are defined from the
thermodynamic analysis. In this process, the dimensions, pressure losses, efficiencies, and other
aerodynamic parameters for the different components can be determined and evaluated.

MATLAB is a widely used numerical computing environment and programming language suitable for
a range of applications. Its ability to read data from external sources, built in functions, advanced
plotting functions and extensive documentation made it an appropriate choice for this project [17].
3.3.2 Turbomachinery
The initial turbomachinery sizing requires a few assumptions from the designer in order to proceed
with further analysis of the component. Based on these assumptions together with the fluid
properties from the thermodynamic analysis the areas, rotational speeds, blade speeds, number of
stages, stage loading etc. can be determined. The initial sizing procedure can be found in [18], where
guidelines for initial assumptions of the necessary parameters are also given. With the initial sizing
procedure done, the parameters described in the following sections can be evaluated, and the
turbomachinery properties can be modified to assure adequate performance.
3.3.2.1 Velocity triangles

In order to design the blades, the power input per stage needs to be related to velocity triangles.
These give an indication of how the blades change from the hub to the tip by showing the variations
in velocity.
The velocity triangles for a typical stage at the mean radius can be seen in figure 4. The fluid
approaches the rotor with an absolute velocity of 1 at an angle 1 relative to the axial direction. By
combining 1 vectorially with the blade speed U, the velocity relative to the blade 1 at an angle 1
relative to the axial direction can be found. As the flow passes through the rotor, the absolute
velocity of the fluid increases leaving the rotor with a velocity 2 at an angle 2 . With the axial
velocity 1 constant over the stage, 2 can be found and the outlet velocity triangle constructed by
combining 2 and U vectorially to give 2 at the angle 2 . The fluid then passes to the stator where
its diffused to a velocity 3 at an angle 3 . Generally 3 approximately equals 1 and 3
approximately equals 1 so that the fluid is prepared for entry into another similar stage [19].
21
Figure 4: Velocity triangles for a typical compressor stage [19]
In this work the axial velocity has been determined assuming values of , and , based on
guidelines given in [18] with a linear interpolation within the component. The velocity triangles have
been calculated at the hub, mean, and tip of each blade of the turbomachinery components. With
regard to radial equilibrium, the exponential blading method was used to determine the angles and
the velocities at the hub and the tip. The necessary equations and formulas can be found in [19].
A number of aerodynamic properties have been evaluated to achieve an adequate turbomachinery

performance including: degree of reaction, diffusion factor, De Haller number, stage loading, flow
coefficient and deflection. Detailed formulas are provided in appendix A.
Degree of reaction
The degree of reaction is defined as the ratio between the static enthalpy change in the rotor and
that of the whole stage. It provides a measure to the extent of which the rotor contributes to the
overall pressure difference in the stage [19]. For a compressor, this should typically be in the range of
0,5 0,8 [20], while for a turbine 0,3 0,5 [11].
Diffusion factor
The air passing over an airfoil will accelerate on the convex side, leading to a drop in static pressure.
On the concave side, the air will decelerate. The losses in a blade row arise primarily from growth of
boundary layers on the two blade sides. At the blade trailing edge the boundary layers come
together forming a wake, causing a local drop in total pressure. Thick boundary layers, causing high
losses have been found in regions where rapid velocity changes occur. The diffusion factor concept
was developed by NACA based on cascade testing. For the rotor hub region and stators the losses are
unaffected up to a diffusion factor of 0,6. In the rotor tip region, the losses grow rapidly at a diffusion
factor above 0,4 [19].
22
De Haller number
The de Haller number is defined as the velocity ratio of the trailing edge and the leading edge of a
blade or a vane. This is a simple measurement of diffusion used in preliminary design work. The De
Haller number should exceed 0.72. For final calculations, the diffusion factor is preferred [19].
Stage loading
Stage loading is a measure of how much work is demanded by a turbomachinery stage and therefore
the required turning of the flow. It is defined as the enthalpy change per unit mass flow of air
divided by the blade speed squared. A lower stage loading leads to a higher efficiency, but in turn,
more stages are required [11].
Flow coefficient
The flow coefficient relates primarily to component size, and if it is non-optimal to the efficiency. It is
defined as the axial velocity divided by the blade speed. A typical value is 0,4 0,8 [20].
Deflection
The deflection is the difference between the angle of the relative velocity of the leading edge and the
trailing edge. It is a direct measure of flow turning. A high deflection would imply a high rate of
diffusion in a compressor or acceleration in a turbine.
3.3.2.2 Flow path

The flow paths for the turbomachinery components have been evaluated using three different design
approaches: constant mean diameter, constant outer diameter and constant inner diameter. The
flow path selection will influence some of the above mentioned parameters, and they all have some
advantages and disadvantages [19]. In modern designs an intermediate choice of constant inner and
constant outer is often made. Figure 5 shows the three designs of a compressor.
Figure 5: Different types of flow paths
23
3.3.2.3 Turbomachinery efficiencies
In an ideal compression or expansion process the assumption is made that there is no change in
entropy. In a real process however, several losses take place as the fluid flows through the
component resulting in an entropy generation. For a compressor this results in a higher temperature
and enthalpy, and hence work, at a given pressure ratio compared to what would be expected from
an ideal process. For a turbine, the outlet pressure will always be lower at a given power output
requirement [15]. The two most common ways to account for losses in the turbomachinery
components is the isentropic and polytropic efficiency.
Isentropic efficiency
The isentropic efficiency of a compressor is defined as the ratio between the change in enthalpy for
an ideal process, and the change in enthalpy for the actual process. For a turbine the isentropic
efficiency is defined as ratio between the change in enthalpy for the actual process, and the change
in enthalpy for the ideal process [15].
Polytropic efficiency
The isentropic efficiency can be misleading when comparing compressors and turbines with different
pressure ratios or at different inlet conditions. If the necessary calculations are performed for each
stage of a multistage component, with the assumption that each stage operates at the same
isentropic efficiency, the result would be that the efficiency of the whole component is lower than
the stage efficiency. Consequently isentropic efficiency is not a suitable parameter when performing
optimization studies. For such studies, the polytropic efficiency can be used. The polytropic efficiency
can be defined as the isentropic efficiency at an infinitely small stage in the compression or
expansion process, such that it can be assumed constant throughout the entire process [15].
Therefore the turbomachinery polytropic efficiencies have been used in this project.
The efficiency calculations are based on outputs from both the thermodynamic analysis and the
sizing procedure, and are based on the following parameters: Entry into service correction, Reynolds
number index effect, size variation and normalized efficiency. These are all based on empirical data.
The method used can be found in [21].
3.3.3 Combustion chamber

In order to commence the combustion chamber sizing, a few assumptions need to be made by the
designer. These assumptions are based on the guidelines in [22] where the full sizing procedure can
also be found. A typical burner can be seen in figure 6.
Loading
Combustor volume must be derived considering the combustor loading at a number of flight
conditions. At static sea level maximum power the combustor loading should be less than

10 1.8 3 . A critical point for the burner is the ability to relight during windmilling following an
engine flameout. In order to achieve this the designer must make sure that the combustor loading at
the highest required altitude and lowest flight Mach number for relight does not exceed

300 1.8 3 [23]. To minimize the volume of the burner, the loading should be close to 300 at
windmilling conditions. A smaller burner however, leads to a higher pressure loss. The definition can
be found in appendix A.
Residence time
The residence time depends on the length of the liner and the Mach number after the diffuser, and
should typically exceed 3ms. Long residence times however, can result in increased NOx.
24
Figure 6: Typical burner [22]
Pressure loss
A pressure loss occurs in the combustion chamber due to skin friction and the rise in temperature.
The skin friction loss is called the cold loss and accounts for approximately 90% of the total pressure
loss. The loss due to rise in temperature is called the fundamental loss and accounts for the
remaining 10% [23]. The pressure loss has been estimated using the method described in reference
[22] and is dependent on the maximum area of the burner, the PLF and the inlet 0 and 0 . The PLF is
an assumed value and is typically in the range of 20 to 25 [23].
Efficiency
The efficiency of the burner can be found from a chemical analysis of the combustion products.
Knowing the fuel to air ratio and the proportion of incompletely burnt constituents, it is possible to
calculate the ratio of the actual energy released to the theoretical quantity available. Due to a
number of reasons this is a very complicated approach. In practice, the combustion is complete at
high-power to mid-power conditions [19] and the burner efficiency has therefore been assumed to
be 0,9999.
3.3.4 Internal ducts

The ducts have been designed using cubic splines to allow for smooth surfaces between the
components trying to minimize the pressure loss. However, no pressure loss calculations have been
performed, the losses are assumptions based on the guidelines given in reference [11].
3.3.5 Inlet
Due to time limitations, only a basic 2D/axisymmetric inlet using oblique shocks to diffuse the
incoming air has been designed. Tabulated values of normal shock properties found in reference [24]
have been used to calculate the total pressure loss and new Mach numbers after the shocks. The
inlet has been designed to create an oblique shock at the inlet start, reducing the Mach number to a
value close to one, after which a normal shock at the venturi reduces the value to below one. After
this, the velocity is diffused in the divergent part using the continuity equation to an acceptable value
before entering the fan. These calculations assume an isentropic flow. To compensate for this the
total pressure loss has been multiplied by a factor of 0,997 for the subsonic region to account for
pressure losses due to skin friction. The necessary formulas and equations can be found in [24]. With
25
the angles and areas known, the geometry of the inlet has been calculated using standard
trigonometric relations.
3.3.6 Nozzle
A convergent divergent nozzle lets the flow expand to ambient pressure, which for an isentropic
process produces maximum thrust [11]. However, due to skin friction losses and the extra weight and
length, a con - di nozzle is only a viable option if the nozzle pressure ratio is greater than about three
[15] which tends to be the case for supersonic engines due to the high ram pressure ratio. The thrust
loss associated with flow angularity increases sharply at great exhaust angles, therefore the angle of
the divergent section should be less than 30 degrees [11].
To account for losses in the nozzle, different coefficients are used, the definitions of these may vary
between different manufacturers. In this thesis the thrust coefficient, CV and the angularity
coefficient are used. CV is defined accordingly to [15]. This coefficient gives the actual thrust from
the thrust that could be achieved by an ideal con di nozzle by accounting for friction and flow non-
uniformity. The value of CV has been assumed based on guidelines given in [11]. is defined
as 0, 5 (1 + cos()), where is the half angle of the divergent section [25]. The area of the nozzle
inlet, throat and exit is calculated in the thermodynamic analysis.
3.3.7 Mixer
The use of a mixer to combine the hot and cold stream prior to the exhaust has a few advantages
compared to using two separate exhausts: a small gain in net thrust and SFC can be achieved, the
optimum fan pressure ratio at a constant specific thrust is lower leading to a lower weight and cost
for both the fan and the LPT, the jet noise is lower due to the velocity of the mixed stream being
significantly lower than the core stream of a separate jets engine. In the decision whether to adopt a
mixer these advantages need to be balanced against the disadvantages of the extra weight and cost
[11]. For the thrust gain to be maximized the mixing needs to be complete (i.e. a mixer efficiency of
100%), this would require a large and heavy mixer, which in turn would require more thrust [1].
Therefore a compromise between mixer efficiency and thrust gain is necessary. Figure 7 shows the
theoretical gross thrust gain by introducing a mixer assuming that: the FPR is at its optimum, the
mixing efficiency is at 100% and no pressure loss is accounted for.
Figure 7: Theoretical thrust gain due to mixing [11]
26
An approximate length and mixer efficiency have been estimated based on guidelines in [11]. To
perform mixer calculations the designer must set the Mach number of the primary incoming stream,
after which the Mach number of the second stream is varied until the static pressures of the two
streams are equal in the mixer. The Mach number has been set based on guidelines given in [19]. The
mixer area for the two streams is calculated in the thermodynamic analysis.
3.3.8 Iterations
The efficiencies of the components highly depend on the choices made in the aerodynamic design,
once these are set, the new values need to be inputs to the thermodynamic analysis, which in turn
will generate new optimum values for the key parameters, and the aerodynamic design will need to
be evaluated again. As mentioned earlier constant feedback is required between the different
disciplines. Iterations were performed until the residual of the turbomachinery efficiencies were less
than 0,0005.
27
3.4 Mechanical design
The discs have been designed using T AXI DISK V2.5, a software developed primarily for education
purposes. The system can be used to design multistage compressors and turbines from a small
number of physical design parameters [26]
3.4.2 Discs
To carry out the disc analysis, the designer needs the aerodynamic properties of the component,
such as stage geometry, RPM and temperatures. Due to time limitation, the disc analysis has only
been carried out on the HPC and the HPT where the disc stresses are the highest due to high disc
speeds. The RPM of the HP shaft is often limited by disc stresses in the HPT rather than the first stage
, for the HPC [27]. The disc designs: hyperbolic, web, continuous slope and ring can be seen in
figure 8.
Figure 8: Different types of disc designs [28]
Target safety line

The target safety line is a safe margin until the material starts to yield, this has been kept at the
default 10%.
Optimization
The optimization was carried out with the T AXI DISK built in function, parametric disc, which
optimizes the disc in order to get a design that meets the target safety line and has the lowest
possible weight. Parameters such as rim width, bore width, bore radius etc. need to be set by the
designer before the optimization can start.
3.4.3 Shafts
The mechanical efficiencies of the shafts are input parameters to the thermodynamic cycle analysis,
in this work they are assumed to be 99,5% for both shafts [11]. No further shaft assumptions or
calculations have been made.
28
4 Results and discussion
4.1 Baseline engine

The modified baseline engine performance data, which served as a starting point for the cycle
optimization can be seen in table 1. A schematic 2D view of the engine can be seen in figure 9.
15545 5%
1,5 2,5%
82270 0,5%
42500 , 0,8914
90/ , 0,9180
472/ , 0,9136
1,068/( ) , 0,9151
16 0,579
1472 0,625
0,3 0,943
3,29 0,546
6% 0,341
Table 1: Key data for the baseline engine
Figure 9: Section view baseline engine
29
4.2 Cycle optimization
A two spool axial flow low by-pass turbofan configuration was selected for the MJ Haran S14. The
engine layout with flow stations can be seen in figure 10.
Figure 10: Engine layout with flow stations
4.2.1 TIT and OPR optimization

The optimum OPR for different TIT:s can be seen in figure 11. SFN and polytropic efficiencies have
been kept constant. Other parameters have been varied to keep the optimum interrelationships
between the components. Figure 11 shows that increasing TIT and OPR increases the core efficiency.
The design is however limited by the blade metal temperature in the turbines. At cruise condition
this temperature was set to 1150K which is representative of a single crystal metal with thermal
barrier coating in the HPT and Inconel 718 in the LPT. Table 2 shows the required cooling flow to
keep this temperature in the HPT, and the inlet temperature of the uncooled LPT for each of the TIT:s
at its optimum OPR.
30
Figure 11: Optimum OPR at different TIT: s.
TIT (K) Opt. OPR Required HPT LPT inlet

cooling flow (%) temperature (K)
35 9.8 1052
37 12 1082
38 16.2 1150
40 20.2 1218
41 24 1278
Table 2: Required HPT cooling flow and LPT inlet temperature
Table 2 shows that the maximum TIT in order to keep the desired LPT inlet temperature is 1644K, at
this temperature the cooling flow for the HPT is also within acceptable limits. Based on this, the TIT
has been set to 1644K and the OPR to 38.
Figure 12 shows the relationship between core efficiency, thermal efficiency, overall efficiency and
SFC at the chosen TIT. It can be seen that the maximum values occur at approximately the same OPR,
which also corresponds to the lowest value of SFC. The propulsive and transfer efficiencies are
relatively constant over this range, and are excluded in the diagram.
31
Figure 12: Efficiencies over OPR
With the OPR and TIT set, the remaining parameters were found using the interrelationships
described in chapter 3. The key parameters for the engine design point cycle can be seen in table 3. A
full output of the engine cycle performance data can be found in appendix B.
15545 9,7%
1,5 6,5%
82270 , 0,901
42500 , 0,941
90/ , 0,91
472/ , 0,934
0,98/( ) 0,632
38 0,624
1644 0,941
0,4 0,595
3,93 0,371
Table 3: Key parameters for design point cycle
32
A schematic 2D view of the complete engine can be seen in figure 13. The engine with inlet can be
seen in figure 14. Note that the discs are preliminary, the final HPC and HPT discs can be found in
appendix D.
Figure 13: The MJ Haran S14 section view
Figure 14: The MJ Haran S14 section view with supersonic inlet
33
4.3 Component design
4.3.1 Inlet
A 2D view of the inlet, with the shocks represented as dotted lines can be seen in figure 15. The inlet
has a total pressure recovery of 0,983 and a total length 2,088 m. The deflection angle was set to
9,26 to achieve maximum pressure recovery in combination with the following normal shock at
design point. The corresponding angle is 55. The remaining inlet data can be seen in table 4.
Figure4.4 Supersonic inlet.
Figure 15: Supersonic inlet with flow stations
Station Mach number PR Area( )

1,500
1,152 0,9897
0,874 0,9966 0,992
0,650 0,9970 1,211
Table 4: Inlet data
4.3.2 Internal ducts

The assumed pressure ratios and the lengths of the internal ducts can be seen in table 5 and a
schematic view can be seen in figure 16.
Duct Bypass Intercompressor Interturbine LPT Jet pipe

duct duct exhaust
PR 0,97 0,98 0,98 0,985 0,99
Length 1,163 0,149 0,107 0,150 0,5
Table 5: Duct data
34
Figure 16: Ducts
4.3.3 Compressors
In order to use an uncooled LPT, a two shaft no-booster design was chosen to maximize the work
done by the HPT, thus lowering the inlet temperature to the LPT. The tip clearance has been set to
1% of the blade height. The blade spacing has been set to 20% of the upstream chord [11]. The two
compressors can be seen on a Smith chart [29], based on the average stage load and flow coefficient
in figure 17.
Figure 17: Stage load over flow coefficient for the FAN and HPC
4.3.3.1 Fan
With a FPR of 3,93, a four stage fan was deemed necessary to achieve satisfactory performance. The
fan flow path is a result of a linear interpolation starting from the hade angle, gradually decreasing to
Figure 18: Fan component
35
zero at the last stage stator outlet. The hade angle was selected to 21 in order to achieve a decent
blade velocity at the first stage rotor hub exit. , has been set high to keep the area to a
minimum.The low pressure shaft RPM was set to reach a , close to 1,45. The fan can be seen
in figure 18.
Table 6 shows some key fan data. Figure 19 shows the fan velocity triangles at the mean blade and
table 7 shows the velocity triangles data at the mean blade. A full output of stage by stage
performance can be found in appendix C. It can be seen in the appendix that the first stage hub
velocity triangles are unrealistic. This is the result of a high stage loading at the hub due to a big
enthalpy change and a low blade speed which indicates that the fan pressure ratio has to vary from a
lower value at the hub region to a higher value at the tip.
, 0,65 4
, 0,38 . 0,82
2,40 . 0,63
2,30 . 42,3
/ 0,27 0,901
/ 0,69 ,, 1,43
2025 / 7000/455,15
/ 1,20 0,9876
21 1,2418
Table 6: Key fan data

53 52 56 58
34 27 34 39
0 12 13 22
34 45 47 50
222 212 197 177
222 212 197 177
290 321 340 351
303 331 346 353
365 348 356 331
269 239 238 226
222 217 202 190
269 300 286 274
0,77 0,66 0,58 0,50
0,94 0,89 0,79 0,65
18 25 22 19
0,74 0,60 0,63 0,60
0,34 0,46 0,47 0,47
0,74 0,69 0,67 0,68
. 39 45 45 39
Table 7: Fan velocity triangles data at mean blade Figure 19: Velocity triangles for the fan at mean blade
36
4.3.3.2 HPC
A five stage constant outer diameter design was chosen for the MJ Haran S14. This configuration
allows for high hub and mean blade speeds through all stages, which reduces the stage loading. The
first stage hub to tip ratio was chosen for an optimal radius of the compressor. Placing it lower would
lead to higher pressure losses in the duct following the splitter due to the steep angle, placing it
higher would lead to a lower RPM in order to keep , at reasonable values, which in turn could
lead to more stages. The final radius is a compromise between these conflicting requirements. To
keep the diffusion factor at the tip below 0,4, the reaction has been adjusted to reduce from 0,75 at
the first stage towards 0,5 at the last stage. With a 50% reaction the blades are symmetrical 1 =
2 = 3 , 1 = 2 giving the most efficient diffusion over the stage [19]. However, a 50% reaction
was not achievable with five stages. An IGV was implemented to increase the absolute velocity of the
first stage, allowing a higher rotational speed while maintaining the same , . The HPC
component is illustrated in figure 20. Note that the discs are preliminary, the final disc design can be
seen in appendix D.
Figure 20: HPC component
Table 8 shows some key compressor data. Figure 21 shows the velocity triangles at the mean blade
and table 9 shows the velocity triangles data at mean blade. A full output of stage by stage
performance for the HPC can be found in appendix C.
5 , 0,52
. 0,58 , 0,33
. 0,38 1,63
. 89,8 0,93
0,941 / 0,70
,, 1,34 2025
/ 13500/613,8 / 1,2
18 0,37
Table 8: Key HPC data
37

67 67 65 67 68
52 57 54 57 58
18 15 26 20 24
48 46 52 51 52
221 225 223 216 205
221 225 223 216 205
523 582 592 596 597
561 587 594 597 598
568 568 534 560 546
357 417 380 394 389
233 233 247 230 225
328 322 361 343 337
0,42 0,39 0,38 0,36 0,34
0,66 0,57 0,56 0,57 0,54
15 9 11 11 10
0,70 0,75 0,67 0,71 0,70
0,47 0,36 0,42 0,42 0,41
0,63 0,73 0,71 0,70 0,71
Figure 21: Velocity triangles for the HPC at mean blade
. 87 92 92 92 87
Table 9: HPC velocity triangles data at mean blade
4.3.4 Combustion chamber

An annular combustor was chosen for the MJ Haran S14 due to its compactness and clean
aerodynamic design. The combustor has been sized to maintain a reasonable pressure loss, while
allowing for relight at windmilling conditions. was set to allow for a residence time in excess of
3ms. A slope angle of 8 was chosen in order to increase the HPT radius, improving its stage loading.
The combustor is shown in figure 22 and data for the combustor can be seen in table 10.
Figure 22: Combustor component
38
/ 4,11 0,04
4 23
1,15 / 2
231 / 1,75
0,35 / 0,5
0,0593 0,999
0,0293 . 2300
8 0,195
Table 10: Key combustor data
4.3.5 Turbines
A constant hub design was chosen for the turbines. This configuration allows for a higher blade speed
due to the increasing blade tip radius, which reduces the stage loading. The loss coefficient was
set to 0,05, at this value the stator blades are convergent, which allows for good performance over a
variety of power conditions [19]. To avoid choking the turbines, the angle 3, has been adjusted

to satisfy the condition ( 01 ) < . The highest relative Mach number occurs at the inlet hub.
2
This value needs to be kept below 0,75 [19] to avoid high pressure losses due to shock waves. The tip
clearance has been set to 1% of the blade length. The blade spacing was set to 20% of the upstream
chord [11]. The two turbines can be seen on a smith chart [30], based on the average stage load and
flow coefficient for the LPT, in figure 23.
Figure 23: Stage load over flow coefficient for the HPT and LPT
4.3.5.1 HPT
A highly loaded one stage HPT was chosen to minimize the required cooling flow. The HPT inlet area
is derived from 2 to maintain stress levels based on material technology for 2025. , was then
derived from the area. To avoid high pressure losses in the upstream ducting and to assure that the
Figure 24: HPT component
39
gas accelerates at all points along the NGV, it has been kept below 0,2 [11]. The HPT is shown in
figure 24. Note that the disc is preliminary, the final disc design can be seen in appendix D.
Some key turbine data can be found in table 11. The velocity triangles for the mean blade can be
seen in figure 25 and the corresponding data can be seen in table 12. A full output of stage by stage
performance can be seen in appendix C.
7247,5 1
, 0,16 . 494
, 0,4 3,10
1,16 0,41
1,16 0.910
0,05 , 0,35
2025 0.097
Table 11: Key HPT data

52
68
75
4
259
630
419
699
994 Figure 25: Velocity triangles for the HPT at mean blade
260
0,41
3,10
16
0,25
, 0,35
. 494
Table 12: HPT velocity triangles data at mean blade
40
4.3.5.2 LPT
The first stage LPT mean radius was set to correspond with the last stage HPT mean radius in order to
achieve a decent stage loading and to approximately align with the fan to reduce flow
turning. , was set high to reduce the rotor outlet area, thus reducing 2 . Since 2 will be
lower on the first stage due to the smaller area, calculations are only performed on the second stage.
The first stage stator was replaced with an aerodynamically loaded strut to accomplish the same
function as a stator, while giving structural support to the engine. With two stages there is no exit
swirl entering the LPT exhaust at design point. Some swirl will however be present at off-design
conditions. The LPT is shown in figure 26.
Figure 26: LPT component
Some key turbine data can be found in table 13. The velocity triangles at the mean blade can be seen
in figure 27 and the corresponding data can be seen in table 14. A full output of stage by stage
performance can be found in appendix C.
, 0,35 2
, 0,5 . 2,06
1,84 . 0,81
6,72 0,934
0,05 ,,, 0,40
2025 ,, 0,44
. 99,7 0,089
5115
Table 13: Key LPT data
41

4 0
53 49
54 49
0 0
254 295
335 344
254 295
420 453
435 453
254 295
2,11 2,00
0,76 0,86
49 49
0,47 0,50
Figure 27: Velocity triangles for the LPT at mean blade
, 0,40 0,44
. 100 100
Table 14: LPT velocity triangles data at mean blade
4.3.6 Mixer
The engine uses a lobed annular mixer to improve thrust and SFC and to reduce noise. The mixer
length was chosen to attain / ~1 to achieve a good percentage of the theoretical thrust gain [11].
The mixer can be seen in figure 28 and the mixer data can be seen in table 15.
0, 7
0, 552
0, 5
0,124 2
0,461 2
0,2
/ 1.0767
Table 15: Key mixer data
Figure 28: Mixer component
42
4.3.7 Nozzle
Due to the high nozzle pressure ratio at supersonic speed, a convergent divergent nozzle was
chosen, which can be seen in figure 29. The areas are set to allow for the fluid to fully expand at
design point. The angles are a compromise between length and thrust, a lower angle will lead to
higher thrust, at the expense of a longer component. For very low angles, the losses due to skin
friction will take over. The nozzle data can be seen in table 16. The rate of change in area as a
function of distance together with the area as a function of distance in the nozzle can be seen in
figure 30.
0, 97
0,988
1,0992
0, 4532
1,042
1,18m
12,5
35,2
13,09
/ 0,44
/ 2,43
Table 16: Key Nozzle data
Figure 29: Nozzle component
Figure 30: Rate of change in area over distance and area over distance for the nozzle
43
4.3.8 Discs
The first stage HPC disc with some stress data can be seen in figure 31. The optimum design was a
hyperbolic disc with the material Inconel 718. However, the disc had to become very heavy and wide
to meet the target safety line due to high rotational speeds, which in turn caused a geometrical
problem with the second stage disc. This was solved by making the gap between the center axis of
the engine and the inner rim of the disc, also called bore radius, different for the two stages. The first
stage has a bore radius of 4cm while the second stage has a bore radius of 9cm.
All the HPC discs, as well as the HPT disc, were designed with a hyperbolic disc and the material
Inconel 718, with a bore radius of 4cm, except for the second stage HPC. Table 17 shows some key
data for all the discs. The full disc analysis can be found in appendix D.
Figure 31: First stage HPC disc output file from T-AXI DISC
() () ()

160,8 4 23,4
38,3 9 7
19,8 4 3
16,9 4 2,5
16,9 4 2,4
104,3 4 16,5
Table 17: Key data for HPC and HPT discs
44
4.4 Off-design
Figure 32 shows the key operating conditions of a typical mission for the aircraft with the thrust
requirements of one engine. Off-design simulations were carried out at these conditions. Some key
turbomachinery off-design data can be seen in table 18 21. The full off-design performance data
can be seen in appendix B.
18000 M = 1,5
16000 Fn = 43kN
M = 1,15
14000 M = 0,98
Fn = 44kN
Fn = 25kN
12000
Altitude (m)
10000
8000
6000 Climb
Fn = 96,5kN
4000
2000
0
0 1 2 3 4 5 6
Time (hr)
Figure 32: Typical mission for the engine
238 7000 0,901 184,6 5542 0,885

54,9 13500 0,941 52,2 12204 0,951
14 13500 0,91 13,7 12204 0,901
59,9 7000 0,934 60,9 5542 0,923
Table 18: M 1.5 alt 51000ft Table 19: M 0 alt 0ft
190,6 5348 0,891 241,2 6714 0,883

52,5 11622 0,951 55,2 12744 0,939
14 11622 0,909 14 12744 0,91
60,9 5348 0,923 59,8 6714 0,935
Table 20: M 0,98 alt 38500ft Table 21: M 1,15 alt 42000ft
45
4.5 Sensitivity analysis
The purpose of a sensitivity analysis is to investigate the effect that a change in efficiency or pressure
loss of a component will have on the efficiency of the engine cycle. This is useful for example when
determining whether the cost of developing more efficient components, will pay back in a more
efficient engine. Table 22 shows the change in efficiency and pressure loss of certain components
and the corresponding change is SFC at design point. Table 23 shows a change in power off-take, HPT
blade metal temperature, customer bleed and the corresponding change in SFC at design point.

, + % 0,47% + % +0,33%
, + % 0,59% + % +0.27%
, + % 0,41% + % +0.28%
, + % 0,24% + % +0,06%
, + , % 0.49% + % +0.28%
+ % 2,52% + % +0.34%
+ % 2,47% + % +0.28%
, + , % 0.24%
Table 22: Design point exchange rates for efficiencies and pressure losses

0,24%
. + 0,19%
+2%
( ) + %
+2,43%
( ) + %
Table 23: Design point exchange rates for power off take, HPT metal temp and customer bleed
46
4.6 Comparison to baseline engine
Table 24 shows a comparison of some key performance parameters for the MJ-Haran S14 and the
baseline engine. Table 25 shows the primary contributors of the increased efficiencies. The slight
decrease of is due to the higher BPR of the MJ-Haran S14, and the slight decrease in is
due to the higher jet velocity, both of these are however very small and are therefore excluded in
table 25.
MJ Haran S14 Baseline Difference (%)

0,632 0,579 + 5,3
0,595 0,546 + 4,9
0,941 0,943 0,2
0,624 0,625 0,1
0,371 0,341 +3
0,981 1,068 8,9
( )

() 4,23 4,71 11,3
() 1,24 1,26 1,6
Table 24: Comparison between baseline engine and MJ - Haran S14
Turbomachinery Increase of OPR Increase of TIT

increase (%) (%) (%)
+0,8 +2,6 +1,9
+0,8 +2,8 +1,3
+0,5 +1,8 +0,7
1,5 4,9 2,2
Table 25: Primary contributors of increased efficiencies
A weighted fuel flow comparison over the key operating points of a typical mission shows a fuel burn
improvement of 11,8% for the MJ-Haran S14 compared to the baseline engine.
47
5 Conclusions
In this thesis a preliminary design of a supersonic jet engine for a conceptual business jet expected to
enter service in 2025 has been carried out. The engine was expected to provide better performance
than that of the current engine chosen to propel the aircraft. Based on the provided engine
requirements, the thermodynamic cycle has been optimized and basic sizing and aerodynamic design
of the main components has been performed. The result is a shorter engine of similar diameter with
improved efficiency. A weighted fuel flow comparison of the two engines at the key operating points
shows a fuel burn improvement of 11,8% for the new engine, improving the economic aspects of
operating the aircraft and lowering the emissions for a reduced environmental impact.
The design process of a jet engine is a complex process which covers many different disciplines and
there is quite often no obvious solution as an improvement of one parameter often comes at the
expense of another.
A secondary objective of the thesis was to evaluate whether the knowledge obtained from the
courses of the aeronautical program at Mlardalen University are sufficient to perform this type of
study, and to recommend courses for students who in the future wish to perform a similar work. The
authors opinion here is that the knowledge obtained is in the lower part of the range of necessary
knowledge, but with proper guidance and determination, its certainly not an impossible task. To
cover all parts of the design process within the limits of the program may not be possible. Courses
covering: thermodynamics, mechanics, aerodynamics, mathematical programming, computer aided
design, aero-engine technology and general mathematics are recommended. Most of these are
available in the current curriculum of the program. A course going deeper in the theory of gas
turbine technology, with a greater focus on the turbomachinery components, could however be
developed to provide a better starting point for a project like this. To follow the method provided in
this thesis, good MATLAB skills are necessary.
5.1 Future work

Only the first stages of the conceptual design process have been carried out and some future work is
necessary to complete the engine.
Further off-design analysis need to be performed, including a turbomachinery analysis at the
off design points.
Further mechanical design including material selections, weight estimations and stress
calculations need to be carried out.
A variable geometry inlet and exhaust design may be necessary to assure good performance
over the entire flight envelope.
48
6 References
[1] Dr. I. Halliwell, Candidate Engines for a Supersonic Business Jet, Joint AIAA
Foundation/ASME.IGTI Student Design Competition 2013/14, Reston, VA, U.S, 2013.
[2] K. G. Kyprianidis, Future Aero Engine Designs: An Evolving Vision, Advances in Gas
Turbine Technology, InTech, Dr. Ernesto Benini (Ed.), ISBN: 978-953-307-611-9 , doi:
10.5772/19689, 2011
[3] K.G. Kyprianidis, T. Grnstedt, S.O.T. Ogaji, P. Pilidis, R. Singh, Assessment of Future
Aero - engine Designs With Intercooled and Intercooled Recuperated Cores , ASME
Journal of Engineering for Gas Turbines and Power, Vol. 133, pp. 011701-1 011701-
10, doi: 10.1115/1.4001982, January, 2011.
[4] T. Lindquist, R. Stieger, M. Hillel,Recuperated Small Gas Turbines Recuperator

Installation, SEAS DTC Technical Conference, Edinburgh, UK, 2008.
[5] L. Larsson, T. Grnstedt, K.G. Kyprianidis, Conceptual Design and Mission Analysis for
a Geared Turbofan and an Open Rotor Configuration, Proceedings of ASME Turbo
Expo, GT2011-46451, Vancouver, 2011.
[6] D.P Raymer, Use of Computers in the Conceptual Design Process, Aircraft Design
Short Course, Dayton, Ohio, 1984.
[7] S.N. Patnaik, D.A Hopkins, General-purpose optimization method for multidisciplinary
design applications, Advances in Engineering Software 31, pp. 57-63, June, 1999.
[8] P. Jeschke, J. Kurzke, R. Schaber, C. Riegler,Preliminary Gas Turbine Design Using the
Multidisciplinary, Design System MOPEDS, Journal of Engineering for Gas Turbines
and Power, Vol. 126, pp. 258 264, April, 2004.
[9] M.J. Jones, S.J. Bradbrook, K. Nurney, A preliminary Engine Design Process for an
Affordable Capability, RTO-MP-089, Paris, France, 2002.
[10] F. Ersavas, Multidisciplinary Conceptual Design of a Transonic High Pressure

Compressor, MSc thesis, Chalmers University of Technology, Gothenburg, Sweden,
2011.
[11] P.P. Walsh, P. Fletcher, Gas Turbine Performance, Second Edition, Blackwell Science,
2004. ISBN: 0-632-06434-X.
[12] http://www.wolverine-
ventures.com/index.php?option=com_content&view=article&id=2&Itemid=2,
[20140804].
49
[13] J. Kurzke, GasTurb12 Design and Off-Design Performance of Gas Turbines,
Technical Manual, Germany, 2012.
[14] A. Guha, Optimum Fan Pressure Ratio for Bypass Engines with Separate or Mixed
Exhaust Streams, Journal of Propulsion and Power, Vol. 17, No. 5, pp. 1117 - 1122,
September - October, 2001.
[15] K.G Kyprianidis, Introduction to Gas Turbine Performance, Chalmers University of

Technology, Gothenburg, Sweden, 2012.
[16] K.G Kyprianidis, MSc course on Gas Turbine Technology, Lecture 3, Chalmers
University of Technology, Gothenburg, Sweden, 2012.
[17] http://www.mathworks.se/products/matlab, [20140804].
[18] T. Grnstedt, Conceptual Aero Engine Design modeling Engine Sizing, Chalmers
[19] H.I.H. Saravanamuttoo, C.F.G. Rogers, H. Cohen, P.V. Straznicky,Gas Turbine Theory,
Sixth Edition, Pearson Prentice Hall, 2009. ISBN: 978-0-13-222437-6.
[20] S.L. Dixon, C.A. Hall, Fluid Mechanics and Thermodynamics of Turbomachinery
Sixth Edition, Elsevier, 2010. ISBN-13: 978-1856177931.
[21] T. Grnstedt, Conceptual Aero Engine Design modeling Efficiency modeling,

Chalmers University of Technology, Gothenburg, Sweden, 2011.
[22] L. Elbrandt, Conceptual design of an annular combustor for the Whittle engine data,
Chalmers University of Technology, Gothenburg, Sweden, 2011.
[24] J.D Anderson, Jr.,Fundamentals of Aerodynamics, Fifth Edition, McGrawer-Hill

Higher Education, 2007. ISBN: 9780071289085.
[25] L.E Stitt,Exhaust Nozzles for Propulsion Systems with Emphasis on Supersonic Cruise
Aircraft, NASA reference publication 1235, Cleveland, Ohio, 1990.
[26] A. Merchant, D. Bruna, M.G. Turner, A Turbomachinery Design Tool for Teaching
Design Concepts for Axial Flow Fans, Compressors and Turbines, Proceedings of
ASME Turbo Expo, GT2006-90105, Barcelona, 2006.
[27] K.G Kyprianidis, personal communication, 2014
[28] D. Gutzwiller, T-Axi Disk V2.2 User's Guide and Tutorial, Cincinnati, Ohio, 2009
50
51
7 Appendices
A Formulas
Propulsive efficiency
0
= 92 02
9 2
Thermal efficiency
1 1
9 92 0 02
2 2
=

Core efficiency
1 2 1
47 47, 02
2 2
=

Transfer efficiency
Overall efficiency
Specific fuel consumption
OPR
03
=
01
52
FPR
021
=
02
TIT
= 04
BPR
13
=
21
Thrust coefficient

=
,
Angularity coefficient
1 + ()
=
2
Gross thrust
= [( 9 92 ) + (9 0 ) 9 ]
Stage loading

= 2
Corrected mass flow
0
=

0

Flow coefficient

=

53
Degree of reaction
2 1
=
02 01
De Haller number
2
=
1
Diffusion factor
2
= 1 +
1 21
Deflection
= 1 2
Combustor loading

= 1.8
, 100.00145(0,400)
Residence time
54
B NPSS outputs
Design point
(All values are in imperial units)
55
Off design
*******************************************************************************
56
57
C MATLAB program output files
-----------------------------FAN-----------------------------
Entry Into Service: 2025
Stages: 4
Av. Stage loading 2dH/Umid2: 1.0223
1st stage Flow coefficient Cax/Umid: 0.76598
1st stage Mtip,rel: 1.4837
1st stage Utip,in: 455.1531
1st stage Utip,in,corr: 435.7497
1st stage Umid,in: 289.5909
1st stage Cax,in: 221.8209
1st stage Hade angle: 21
(rhub/rtip)in: 0.2725
(rhub/rtip)out: 0.69474
AR,in 2.4
AR,out 2.3
M,ax,in 0.65
M,ax,in 0.38
Total length [m]: 0.9876
RNI: 0.35315
Mcorr: 238.1475
dETAEIS: 0.016354
dETARe: -0.012639
dETAM: 0
ETA*: 0.89693
ETApol: 0.90065
-----------------------------HPC-----------------------------
Stages: 5
1:st stg Flow function HPC: 0.42313
Avg. Stage load HPC: 0.7074
1:st stg(hub/tip)-ratio: 0.70373
last stg(hub/tip)-ratio: 0.94776
1:st stg Mtip,rel no IGV: 1.5247
AR,in 1.6325
AR,out 0.9325
M,ax,in 0.517
M,ax,in 0.329
RNI: 0.80258
Mcorr: 54.9092
dETAEIS: 0.018391
dETARe: -0.0021124
dETAM: -0.0014644
ETA*: 0.9258
ETApol: 0.94062
-----------------------------Combustor-----------------------
At Cruise conditions
Pattern Factor: 0.35259
dp: 4.1152
burning time: 0.004021
Vcc [m^3]: 0.058889
Vcc [m^3]: 0.029445
At windmilling conditions
Loading: 230.9597
58
-----------------------------HPT-----------------------------
Stages: 1
1:st stg Flow function HPT: 0.41448
Avg. Stage load HPT: 3.1288
AN2: 7247.5
AR,in 1.16
AR,out 1.16
M,ax,in 0.15852
M,ax,in 0.4
RNI: 1.822
Mcorr: 9.0066
dETAEIS: 0.0084819
dETARe: 0.0056299
dETAM: -0.034234
ETA*: 0.95693
ETApol: 0.9099
-----------------------------LPT-----------------------------
Stages: 2
1:st stg Flow function LPT: 0.75935
Avg. Stage load LPT: 2.1106
Last stg AN2: 5114.8347
AR,in 1.835
AR,out 6.7225
M,ax,in 0.35
M,ax,in 0.5
RNI: 0.64113
Mcorr: 38.6614
dETAEIS: 0.016457
dETARe: -0.0045817
dETAM: -0.0082715
ETA*: 0.93012
ETApol: 0.93372
-----------------------------Engine--------------------------
OPR: 38.0
BPR: 0.39702
FPR: 3.9283
HPC P-ratio: 10.0381
Inlet Diameter [m]: 1.2418
Inlet Area [m]: 1.2112
Engine Length [m]: 4.2305
Core length: FANstart-LPTend [m]: 2.0971
LowSpeedShaft [rps]: 116.6667
HighSpeedShaft [rps]: 225
-----------------------------Ducts---------------------------
ICD length [m]: 0.14886
ITD length [m]: 0.1067
LPT exhaust length [m]: 0.15
BP duct length [m]: 1.1627
Jet pipe length [m]: 0.5
Mixer length [m] 0.2
-----------------------------Nozzle--------------------------
Nozzle length [m] 1.1834
Petal angle [deg]: 35.2059
Divergent angle [deg]: 12.5
Throat Diameter [m]: 0.75978
Exit Diameter [m]: 1.1515
59
Ain/Athroat: 2.4249
Athroat/Aexit: 0.43538
Athroat: 0.45338
Aexit: 1.0414
Ain: 1.0994
-----------------------------Inlet---------------------------
Total Pressure recovery Inlet: 0.98342
Total length Inlet: 2.0884
Total Engine length with Inlet: 6.3189
Venturi Area [m^2]: 0.99211
Oblique shock angle "beta" [deg]: 55
Deflection angle "theta" [deg]: 9.2563
Mach number after oblique shock: 1.1517
Mach number after normal shock: 0.87387
Mach number inlet FAN: 0.65
----------------------------Velocity Diagram FAN------------
Mtip,rel: 1.4297
Average Stage Load: 0.81845
Average Flow coefficient: 0.62704
Average Temp rise: 42.2583
Stage Hubradii Midradii Tipradii

_____ ________ ________ ________
beta1 1 51.128 52.549 63.069

beta2 1 -10.74 34.427 55.205
alpha1 1 -32.218 0 5.216
alpha2 1 41.505 34.311 34.031
Deflection 1 61.868 18.122 7.864
C-axial in 1 203.19 221.82 220.97
C-axial out 1 254.59 221.82 203.29
Uin m/s 1 124.03 289.59 455.15
Uout m/s 1 176.99 303.42 429.84
V1 m/s 1 323.77 364.78 487.89
V2 m/s 1 259.13 268.92 356.26
C1 m/s 1 240.17 221.82 221.89
C2 m/s 1 339.95 268.55 245.31
Reaction 1 0.60803 0.73863 0.82703
Diffusion 1 0.38456 0.33965 0.2986
de Haller 1 0.63877 0.73721 0.7937
Flow coefficient 1 1.6383 0.76598 0.48549
Stage load 1 5.1384 0.94255 0.38156
Stage Temp-rise K 1 39.258 39.258 39.258

_____ ________ ________ ________
beta1 2 44.776 52.446 59.916

beta2 2 2.009 27.241 50
alpha1 2 2.742 11.984 17.446
alpha2 2 43.977 44.945 46.873
Deflection 2 42.767 25.204 9.9159
60
C-axial in 2 218.11 212.15 203.43
C-axial out 2 248.13 212.15 182.76
Uin m/s 2 226.85 320.97 415.1
Uout m/s 2 248.13 330.52 412.92
V1 m/s 2 307.25 348.07 405.83
V2 m/s 2 248.28 238.62 284.32
C1 m/s 2 218.36 216.88 213.24
C2 m/s 2 344.81 299.74 267.34
Reaction 2 0.44926 0.6 0.68797
Diffusion 2 0.43885 0.45729 0.37939
de Haller 2 0.71031 0.68555 0.77984
Stage load 2 1.7776 0.88795 0.53092
Stage Temp-rise K 2 45.258 45.258 45.258

_____ ________ ________ ________
beta1 3 50.622 56.374 61.686

beta2 3 12.93 34.122 51.766
alpha1 3 6.5513 12.817 17.186
alpha2 3 45.27 46.494 48.377
Deflection 3 37.692 22.252 9.9208
C-axial in 3 202.08 196.95 190.49
C-axial out 3 224.99 196.95 172.61
Uin m/s 3 269.41 340.95 412.49
Uout m/s 3 278.78 346.06 413.33
V1 m/s 3 318.51 355.65 401.62
V2 m/s 3 230.85 237.91 278.91
C1 m/s 3 203.4 201.98 199.39
C2 m/s 3 319.7 286.08 259.86
Reaction 3 0.53541 0.63 0.69312
Diffusion 3 0.49431 0.47296 0.39589
de Haller 3 0.65094 0.66892 0.76638
Stage load 3 1.268 0.79169 0.54089
Stage Temp-rise K 3 45.258 45.258 45.258

_____ ________ ________ ________
beta1 4 51.145 57.728 63.452

beta2 4 17.476 38.639 55.698
alpha1 4 16.638 21.913 26.171
alpha2 4 47.237 49.873 53.432
Deflection 4 33.669 19.088 7.754
C-axial in 4 186.71 176.87 166.42
C-axial out 4 207.95 176.87 147.88
Uin m/s 4 287.56 351.22 414.88
Uout m/s 4 290.32 353.22 416.12
V1 m/s 4 297.62 331.25 372.35
V2 m/s 4 218.01 226.43 262.4
C1 m/s 4 194.87 190.64 185.43
C2 m/s 4 306.27 274.43 248.21
Reaction 4 0.51202 0.6 0.66118
Diffusion 4 0.49207 0.46819 0.39644
61
de Haller 4 0.65771 0.68358 0.79309
Stage load 4 0.97205 0.65163 0.46701
Stage Temp-rise K 4 39.258 39.258 39.258
62
----------------------------Velocity Diagram Compressor-----
IGV angle: 18
First stage Mtip,rel: 1.3421
Average Stage Load: 0.58116
Average Flow coefficient: 0.37835
Average Temp rise: 89.8067

_____ ________ ________ ________
beta1 1 59.191 67.066 68.265

beta2 1 46.582 51.759 61.289
alpha1 1 11.848 18 20.69
alpha2 1 47.102 47.582 49.396
Deflection 1 12.609 15.306 6.9764
C-axial in 1 228.95 221.24 212.68
C-axial out 1 238.25 221.24 205.13
Uin m/s 1 431.95 522.88 613.8
Uout m/s 1 508.18 560.99 613.8
V1 m/s 1 447 567.76 574.31
V2 m/s 1 346.63 357.44 427.01
C1 m/s 1 233.93 232.63 227.34
C2 m/s 1 350 327.99 315.19
Reaction 1 0.64761 0.69971 0.73964
Diffusion 1 0.38034 0.47122 0.34827
de Haller 1 0.74519 0.62956 0.77085
Stage load 1 0.96328 0.65739 0.47705
Stage Temp-rise K 1 86.807 86.807 86.807

_____ ________ ________ ________
beta1 2 65.382 66.703 67.942

beta2 2 54.619 57.398 60.589
alpha1 2 13.717 15.017 16.21
alpha2 2 45.398 45.763 46.162
Deflection 2 10.762 9.3045 7.3534
C-axial in 2 226.76 224.67 222.5
C-axial out 2 231.52 224.67 218.02
Uin m/s 2 550.22 582.01 613.8
Uout m/s 2 560.78 587.29 613.8
V1 m/s 2 544.35 568.06 592.48
V2 m/s 2 399.87 416.98 443.97
C1 m/s 2 233.42 232.61 231.72
C2 m/s 2 329.72 322.04 314.78
Reaction 2 0.73637 0.75 0.76235
Diffusion 2 0.37158 0.36042 0.33478
de Haller 2 0.71946 0.73405 0.76474
Stage load 2 0.63967 0.5717 0.51401
Stage Temp-rise K 2 91.807 91.807 91.807
63
_____ ________ ________ ________
beta1 3 64.154 65.349 66.49

beta2 3 51.159 54.142 57.217
alpha1 3 24.576 25.624 26.638
alpha2 3 51.163 51.895 52.671
Deflection 3 12.995 11.207 9.2734
C-axial in 3 226.09 222.66 219.18
C-axial out 3 231.05 222.66 214.32
Uin m/s 3 570.14 591.97 613.8
Uout m/s 3 573.93 593.87 613.8
V1 m/s 3 518.61 533.83 549.47
V2 m/s 3 368.4 380.1 395.81
C1 m/s 3 248.62 246.94 245.21
C2 m/s 3 368.43 360.8 353.43
Reaction 3 0.65764 0.67 0.68152
Diffusion 3 0.42826 0.41595 0.39781
de Haller 3 0.69513 0.71203 0.73672
Stage load 3 0.60805 0.56402 0.52461
Stage Temp-rise K 3 91.807 91.807 91.807

_____ ________ ________ ________
beta1 4 66.495 67.293 68.063

beta2 4 54.411 56.763 59.071
alpha1 4 19.353 20.196 21.006
alpha2 4 50.423 50.924 51.45
Deflection 4 12.084 10.53 8.9918
C-axial in 4 217.95 216.04 214.1
C-axial out 4 222.22 216.04 209.93
Uin m/s 4 577.69 595.75 613.8
Uout m/s 4 579.35 596.58 613.8
V1 m/s 4 546.48 559.66 573.09
V2 m/s 4 381.84 394.16 408.44
C1 m/s 4 231 230.19 229.34
C2 m/s 4 348.79 342.73 336.86
Reaction 4 0.70106 0.71 0.71843
Diffusion 4 0.4306 0.41636 0.39987
de Haller 4 0.68531 0.70429 0.72686
Stage load 4 0.60446 0.56837 0.53542
Stage Temp-rise K 4 91.807 91.807 91.807

_____ ________ ________ ________
beta1 5 67.134 67.912 68.663

beta2 5 55.841 58.116 60.285
alpha1 5 23.151 23.952 24.731
64
alpha2 5 51.863 52.449 53.064
Deflection 5 11.293 9.7957 8.378
C-axial in 5 207.58 205.4 203.2
C-axial out 5 211.7 205.4 199.14
Uin m/s 5 580.98 597.39 613.8
Uout m/s 5 581.61 597.71 613.8
V1 m/s 5 534.2 546.24 558.49
V2 m/s 5 377.03 388.87 401.75
C1 m/s 5 225.76 224.76 223.72
C2 m/s 5 342.8 337.02 331.39
Reaction 5 0.69157 0.7 0.70798
Diffusion 5 0.42568 0.4114 0.39632
de Haller 5 0.69205 0.71191 0.73404
Stage load 5 0.57551 0.54432 0.51561
Stage Temp-rise K 5 86.807 86.807 86.807
65
-----------------------------Velocity Diagram HPT------------
Avg. Stage load 3.101

Avg. Flow coefficient 0.41135

_____ ________ ________ ________
beta2 1 58.989 51.748 42.042

beta3 1 66.991 68.206 70.037
alpha2 1 75.772 74.874 73.983
alpha3 1 4.2634 4 3.6162
Deflection 1 8.0023 16.458 27.995
C-axial in 1 259.33 259.33 259.33
U m/s 1 591.35 630.43 669.52
V2 m/s 1 503.36 418.86 349.2
V3 m/s 1 663.47 698.49 759.58
C2 m/s 1 1055.1 993.8 939.85
C3 m/s 1 260.05 259.96 259.85
Relative Mach # 1 0.42712 0.35277 0.2923
Reaction 1 0.19342 0.25351 0.3572
Stage load 1 3.101 3.101 3.101
Stage Temp-drop K 1 494.16 494.16 494.16
P2/Pchoke < 1 1 0.85975 0.85975 0.85975
66
-----------------------------Velocity Diagram LPT------------
Avg. Stage load 2.056

Avg. Flow coefficient 0.80783

_____ ________ ________ ________
beta2 1 21.753 4.1707 -12.257

beta3 1 49.486 52.837 56.013
alpha2 1 57.493 54.309 51.36
alpha3 1 0 0 0
Deflection 1 27.732 48.666 68.27
C-axial in 1 253.6 253.6 253.6
U m/s 1 296.77 334.54 372.31
V2 m/s 1 273.04 254.27 259.51
V3 m/s 1 390.37 419.8 453.66
C2 m/s 1 471.9 434.68 406.13
C3 m/s 1 253.6 253.6 253.6
Relative Mach # 1 0.43013 0.3978 0.40403
Reaction 1 0.35946 0.47236 0.59675
Stage load 1 2.1106 2.1106 2.1106
Stage Temp-drop K 1 99.661 99.661 99.661
P2/Pchoke < 1 1 0.73794 0.73794 0.73794

_____ ________ ________ ________
beta2 2 18.955 0.046184 -16.582

beta3 2 45.207 49.383 54.209
alpha2 2 53.485 49.402 45.763
alpha3 2 0 0 0
Deflection 2 26.252 49.337 70.79
C-axial in 2 294.64 294.64 294.64
U m/s 2 296.77 343.55 390.32
V2 m/s 2 311.53 294.64 307.42
V3 m/s 2 418.19 452.59 503.79
C2 m/s 2 495.16 452.77 422.34
C3 m/s 2 294.64 294.64 294.64
Relative Mach # 2 0.47108 0.44216 0.45906
Reaction 2 0.36351 0.49965 0.68916
Stage load 2 2.0014 2.0014 2.0014
Stage Temp-drop K 2 99.661 99.661 99.661
P2/Pchoke < 1 2 0.73757 0.73757 0.73757
67
68
D T-AXI DISK output files
First HPC disc
69
&INPUT_OUTPUT
C1="*******T-AXI_DISK_INPUT_FILE:_VERSION_2.4*******"
C2="*******PARAMATERIZATION_TYPE********************"
TYPE= 2
C3="*****************PARAMETERS*********************"
DRQR= 1.00010000169277191E-002
WWEB= 0.44678959250450134
DRQB= 0.57624977827072144
WBOR= 0.23997600376605988
RBOR= 4.00040000677108765E-002
DSF= 0.16814558207988739
PR1= 0.50000000000000000
PR2= 0.50000000000000000
PT1= 1.00000000000000000
PT2= 0.25000000000000000
PT3= 0.25000000000000000
S12= 0.50000000000000000
C4="*****************DISK_MATERIAL******************"
MATNAME="INCONEL718.dmat"
C5="**********DEAD_WEIGHT_SPECIFICATION*************"
DEAD_WEIGHT_FLAG= 1
M_B= 1.1856651650167098
R_CG_B= 0.35306778472642797
C6="**************ANALYSIS_CONDITIONS***************"
RPMO= 13500.000000000000
TADDER= 0.0000000000000000
RIMT= 150.00000000000000
TAXIS= 145.00000000000000
SFACT= 1.1000000000000001
C7="*****************OTHER_VALUES*******************"
BNUM= 65.000000000000000
BIRR= 0.30900000000000000
BERR= 0.35899999999999999
WRIM= 7.79999999999999999E-002
BSPN= 0.12500000000000000
BSHR= 2
BRHO= 8000.0000000000000
BTHP= 7.0000000000000000
RTQB= 0.25000000000000000
WBOR_MAX= 0.23997599999999999
RBOR_MIN= 4.00039999999999979E-002
BFILENAME="none"
BSCALE= 1.00000000000000000 /
70
Second HPC disc
71
&INPUT_OUTPUT
TYPE= 2
C3="*****************PARAMETERS*********************"
DRQR= 1.00010000169277191E-002
WWEB= 0.10898753255605698
DRQB= 1.1600420475006104
WBOR= 5.99999986588954926E-002
RBOR= 9.00089964270591736E-002
DSF= 0.10001000016927719
PR1= 0.50000000000000000
PR2= 0.50000000000000000
PT1= 1.00000000000000000
PT2= 0.25000000000000000
PT3= 0.25000000000000000
S12= 0.50000000000000000
C4="*****************DISK_MATERIAL******************"
DEAD_WEIGHT_FLAG= 1
M_B= 0.11507518039625086
R_CG_B= 0.39494879717528664
RPMO= 13500.000000000000
TADDER= 0.0000000000000000
RIMT= 150.00000000000000
TAXIS= 145.00000000000000
SFACT= 1.1000000000000001
C7="*****************OTHER_VALUES*******************"
BNUM= 60.000000000000000
BIRR= 0.38900000000000001
BERR= 0.39700000000000002
WRIM= 2.99999999999999989E-002
BSPN= 4.49999999999999983E-002
BSHR= 2
BRHO= 8000.0000000000000
BTHP= 7.0000000000000000
RTQB= 0.20000000000000001
WBOR_MAX= 6.99930000000000135E-002
RBOR_MIN= 9.00089999999999918E-002
BFILENAME="none"
BSCALE= 1.00000000000000000 /
72
Third HPC disc
73
&INPUT_OUTPUT
TYPE= 2
C3="*****************PARAMETERS*********************"
DRQR= 9.99999974737875164E-005
WWEB= 8.00120010972023010E-002
DRQB= 1.2000399827957153
WBOR= 2.49975007027387619E-002
RBOR= 4.00040000677108765E-002
DSF= 0.18009099364280701
PR1= 0.50000000000000000
PR2= 0.50000000000000000
PT1= 1.00000000000000000
PT2= 0.25000000000000000
PT3= 0.25000000000000000
S12= 0.50000000000000000
C4="*****************DISK_MATERIAL******************"
DEAD_WEIGHT_FLAG= 1
M_B= 5.81338048435073584E-002
R_CG_B= 0.40494046444841159
RPMO= 13500.000000000000
TADDER= 0.0000000000000000
RIMT= 150.00000000000000
TAXIS= 145.00000000000000
SFACT= 1.1000000000000001
C7="*****************OTHER_VALUES*******************"
BNUM= 60.000000000000000
BIRR= 0.40300000000000002
BERR= 0.40600000000000003
WRIM= 2.29999999999999996E-002
BSPN= 3.09999999999999998E-002
BSHR= 2
BRHO= 8000.0000000000000
BTHP= 7.0000000000000000
RTQB= 0.20000000000000001
WBOR_MAX= 2.99969999999999994E-002
RBOR_MIN= 4.00039999999999979E-002
BFILENAME="none"
BSCALE= 1.00000000000000000 /
74
Fourth HPC disc
75
&INPUT_OUTPUT
TYPE= 2
C3="*****************PARAMETERS*********************"
DRQR= 1.00010000169277191E-002
WWEB= 6.81580975651741028E-002
DRQB= 0.86708587408065796
WBOR= 2.46679987758398056E-002
RBOR= 4.00040000677108765E-002
DSF= 0.15354484319686890
PR1= 0.50000000000000000
PR2= 0.50000000000000000
PT1= 1.00000000000000000
PT2= 0.25000000000000000
PT3= 0.25000000000000000
S12= 0.50000000000000000
C4="*****************DISK_MATERIAL******************"
DEAD_WEIGHT_FLAG= 1
M_B= 4.14624665236959580E-002
R_CG_B= 0.40971401261695956
RPMO= 13500.000000000000
TADDER= 0.0000000000000000
RIMT= 150.00000000000000
TAXIS= 145.00000000000000
SFACT= 1.1000000000000001
C7="*****************OTHER_VALUES*******************"
BNUM= 60.000000000000000
BIRR= 0.40899999999999997
BERR= 0.40999999999999998
WRIM= 2.19999999999999987E-002
BSPN= 2.29999999999999996E-002
BSHR= 2
BRHO= 8000.0000000000000
BTHP= 7.0000000000000000
RTQB= 0.20000000000000001
WBOR_MAX= 2.49975000000000024E-002
RBOR_MIN= 4.00039999999999979E-002
BFILENAME="none"
BSCALE= 1.00000000000000000 /
76
Fifth HPC disc
77
&INPUT_OUTPUT
TYPE= 2
C3="*****************PARAMETERS*********************"
DRQR= 1.00010000169277191E-002
WWEB= 6.48523047566413879E-002
DRQB= 0.90963256359100342
WBOR= 2.29748580604791641E-002
RBOR= 4.00040000677108765E-002
DSF= 0.16116099059581757
PR1= 0.50000000000000000
PR2= 0.50000000000000000
PT1= 1.00000000000000000
PT2= 0.25000000000000000
PT3= 0.25000000000000000
S12= 0.50000000000000000
C4="*****************DISK_MATERIAL******************"
DEAD_WEIGHT_FLAG= 1
M_B= 4.39481341670272668E-002
R_CG_B= 0.41129828749125069
RPMO= 13500.000000000000
TADDER= 0.0000000000000000
RIMT= 150.00000000000000
TAXIS= 145.00000000000000
SFACT= 1.1000000000000001
C7="*****************OTHER_VALUES*******************"
BNUM= 60.000000000000000
BIRR= 0.41099999999999998
BERR= 0.41099999999999998
WRIM= 2.29999999999999996E-002
BSPN= 2.29999999999999996E-002
BSHR= 2
BRHO= 8000.0000000000000
BTHP= 7.0000000000000000
RTQB= 0.20000000000000001
WBOR_MAX= 2.39976000000000009E-002
RBOR_MIN= 4.00039999999999979E-002
BFILENAME="none"
BSCALE= 1.00000000000000000 /
78
HPT disc
79
&INPUT_OUTPUT
TYPE= 2
C3="*****************PARAMETERS*********************"
DRQR= 1.00010000169277191E-002
WWEB= 0.15955232083797455
DRQB= 0.90005499124526978
WBOR= 0.16498349606990814
RBOR= 3.00030000507831573E-002
DSF= 0.14009299874305725
PR1= 0.50000000000000000
PR2= 0.50000000000000000
PT1= 1.00000000000000000
PT2= 0.25000000000000000
PT3= 0.25000000000000000
S12= 0.50000000000000000
C4="*****************DISK_MATERIAL******************"
DEAD_WEIGHT_FLAG= 1
M_B= 0.18910831139661863
R_CG_B= 0.42570518476269659
RPMO= 13500.000000000000
TADDER= 0.0000000000000000
RIMT= 600.00000000000000
TAXIS= 145.00000000000000
SFACT= 1.1000000000000001
C7="*****************OTHER_VALUES*******************"
BNUM= 68.000000000000000
BIRR= 0.41799999999999998
BERR= 0.41799999999999998
WRIM= 4.39999999999999974E-002
BSPN= 5.39999999999999994E-002
BSHR= 2
BRHO= 8000.0000000000000
BTHP= 7.0000000000000000
RTQB= 0.14999999999999999
WBOR_MAX= 0.16498350000000001
RBOR_MIN= 3.00029999999999984E-002
BFILENAME="none"
BSCALE= 1.00000000000000000 /
80

Report Final PDF

Uploaded by

Copyright:

Available Formats

Report Final PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Report Final PDF

Uploaded by

Copyright:

Available Formats

School of Innovation, Design and Engineering

Vsters, September 2014

Axial Mach number

Blade Stress Level 2 2 2

Corrected mass flow /

Entry into Service

Fan Pressure Ratio

Fuel to Air Ratio

High Pressure Compressor

High Pressure Turbine

Inlet Guide Vane

Loss Coefficient for Stator Blades

Lower Heating Value /

Nozzle Guide Vane

Overall Pressure Ratio

Pressure Loss Factor

Relative Blade Mach number

Revolutions per Minutes

Revolutions per Seconds

Specific Fuel Consumption

Turbine Inlet Temperature

Water to Air Ratio

Engine development is a multidisciplinary process covering thermodynamics, aerodynamics and solid

1.2 Aim and objectives

1.3 Work scope and limitations

1.4 Thesis outline

2.1 Reducing fuel consumption

2.3 Cost of development

Figure 1: Preliminary engine design process [10]

3.2 Thermodynamic analysis

3.2.2 Software description

3.2.3 Cycle modeling and optimization

3.2.3.1 Cycle efficiency

Detailed formulas are provided in appendix A.

3.2.3.2 Selection of cycle parameters

Fan pressure ratio

Turbine inlet temperature

Overall pressure ratio

3.3.1 Software description

3.3.2.1 Velocity triangles

A number of aerodynamic properties have been evaluated to achieve an adequate turbomachinery

3.3.2.2 Flow path

Figure 5: Different types of flow paths

3.3.3 Combustion chamber

3.3.4 Internal ducts

Figure 7: Theoretical thrust gain due to mixing [11]

Figure 8: Different types of disc designs [28]

Target safety line

4.1 Baseline engine

Figure 9: Section view baseline engine

Figure 10: Engine layout with flow stations

4.2.1 TIT and OPR optimization

TIT (K) Opt. OPR Required HPT LPT inlet

Figure 13: The MJ Haran S14 section view

Figure4.4 Supersonic inlet.

Figure 15: Supersonic inlet with flow stations

Station Mach number PR Area( )

4.3.2 Internal ducts

Duct Bypass Intercompressor Interturbine LPT Jet pipe