<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJEE</journal-id><journal-title-group><journal-title>Open Journal of Energy Efficiency</journal-title></journal-title-group><issn pub-type="epub">2169-2637</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojee.2022.113007</article-id><article-id pub-id-type="publisher-id">OJEE-119860</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Off-Design Modelling of a Turbo Jet Engine with Operative Afterburner
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>O.</surname><given-names>O. Ighodaro</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>E.</surname><given-names>I. Aburime</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>A.</surname><given-names>A. Erameh</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Chemical Engineering, University of Benin, Benin City, Nigeria</addr-line></aff><aff id="aff1"><addr-line>Department of Mechanical Engineering, State University of Benin, Benin City, Nigeria</addr-line></aff><aff id="aff3"><addr-line>Department of Mechanical Engineering, Igbinedion University, Okada, Nigeria</addr-line></aff><pub-date pub-type="epub"><day>06</day><month>09</month><year>2022</year></pub-date><volume>11</volume><issue>03</issue><fpage>88</fpage><lpage>107</lpage><history><date date-type="received"><day>11,</day>	<month>April</month>	<year>2022</year></date><date date-type="rev-recd"><day>16,</day>	<month>September</month>	<year>2022</year>	</date><date date-type="accepted"><day>19,</day>	<month>September</month>	<year>2022</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Gas turbines are considered as one of the leading internal combustion engines in modern air transportation due to its favourable power to weight ratio and its continuous combustion process. Recent research focus h
  as
   been concerned with performance improvements aimed at reduced fuel consumption and hence reduced impact on the environment. This study is aim
  ed
   at using theoretical and computational methods to model the operation and performance a turbojet gas turbine engine. The commercial software GasTurb13 was used for the theoretical simulation while Microsoft Excel was used for the analytical study. GasTurb13 solved the model using pseudo-perfect gas models i.e. component maps since the specific gas ratio could not be inputted into the solver. The effect of changes in the Mach number and altitude on the engine performance was studied. Also the effect of changes in the compressor pressure ratio, the turbine inlet temperature and the afterburner exit temperature were also studied. Results obtained showed the optimum pressure ratio at maximum thrust constraint to be 16.78 for the turbojet engine operating at Mach number (Ma)
   
  = 0.8 and altitude = 10
  ,
  000
   
  m, Turbine inlet temperature (TIT) = 1200
   
  K and Afterburner exit temperature = 1800
   
  K
  .
 
</p></abstract><kwd-group><kwd>Turbo Jet Engine</kwd><kwd> Modelling</kwd><kwd> Simulation</kwd><kwd> Off-Design</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Gas turbines are among the leading internal combustion engines in the transportation industry. They are used mainly for aviation and stationary electricity generation purposes. They play an important role in the aviation industry due to their high power to weight ratio, are lighter and smaller than comparable engines and are cost-efficient in operation [<xref ref-type="bibr" rid="scirp.119860-ref1">1</xref>]. More and more attention is been paid on optimizing fuel consumption of these engines while reducing their noise level.</p><p>Gobran [<xref ref-type="bibr" rid="scirp.119860-ref2">2</xref>] studied the off-design performance of solar centaur-40 gas turbine engine. A Simulink model was used in simulating the off-design running for the single shaft Centaur 40 gas turbine engine used for power generation. The off-design modelling was done in two phases, initially, the engine was operated at speeds ranging between 65% to 100% at no load conditions, subsequently, the engine was loaded at a constant speed of 100%. The results obtained from this test were compared with the actual operation results of the turbine to check the validity of the Simulink model. He then went to investigate the effect of ambient temperature on engine performance at engine design conditions.</p><p>Lazzarettoand Toffolo [<xref ref-type="bibr" rid="scirp.119860-ref3">3</xref>] used neural network models for their studies on gas turbine. Both design and off-design simulations were carried out. The results obtained from their analytical studies were used in training the neural network model, they further studied the capability of using the developed neural network model to predict the performance of the turbine.</p><p>The predicting of off design performance of gas turbine engines using graphical analysis was carried in studies conducted by Shapiro and Caddy [<xref ref-type="bibr" rid="scirp.119860-ref4">4</xref>] and Wittenberg [<xref ref-type="bibr" rid="scirp.119860-ref5">5</xref>]. In the absence of manufactures component maps, gas dynamic relations were used in the analysis. The results obtained were far from accuracy due to the many assumptions used in the analysis.</p><p>Sellers and Daniels [<xref ref-type="bibr" rid="scirp.119860-ref6">6</xref>] used the software DYNGEN for the dynamic simulations of both jet and turbofan engine dynamic, they obtained poor results due to the software’s numerical stability issues and its poor user interface.</p><p>The so called “hot end method” used in predicting the steady state performance of LM-600 gas turbine engine has been used by a number of authors [<xref ref-type="bibr" rid="scirp.119860-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.119860-ref8">8</xref>]. The method is based on using simplified matching equations on the condition that the low pressure turbine is chocked (fixed high pressure turbine operating point). The solution method starts from hot end down to the compressor entry point.</p><p>Najjar and Balawneh [<xref ref-type="bibr" rid="scirp.119860-ref9">9</xref>] carried out optimization studies on turbojet engines. Initial results showed that the specific thrust from the engine depends heavily on the turbine inlet temperature and consequently on the specific fuel consumption (SFC) of the engine. They work rigorously in achieving minimal specific fuel consumption for a maximum thrust and after optimizing the turbojet cycle, they found the optimal pressure ratio for minimal specific fuel consumption.</p><p>Bakalis &amp; Stamatis [<xref ref-type="bibr" rid="scirp.119860-ref10">10</xref>] developed a small model of turbojet engine which was well calibrated and tried to calculate for the efficiency which proved to be inadequate for producing reasonable values. They suggested the possible reasons for the inadequacy were due to ignoring the heat transfer occurring in the turbine casing and errors when taking measurements. The analytical calculations carried out for the turbine inlet temperature (TIT) and the turbine exit total pressures deviated significantly from the predictions. In trying to resolve the differences, the model was recalibrated to use static pressure instead of total pressure measurements in the turbine exit, the recalibration produced much more reasonable results both for predicting measured parameters and for the turbine efficiency, their research concludes that for small engines, static pressure measurements at hot regions produces more reliable results than total pressure measurements. Wood &amp; Pilidis [<xref ref-type="bibr" rid="scirp.119860-ref11">11</xref>] also investigated a jet engine based on a variable cycle for a supersonic short take-off vertical landing (ASTOVL) aircraft. They analysed the on-design and off-design parameters for the engine</p><p>Najjar &amp; Al-Sharif [<xref ref-type="bibr" rid="scirp.119860-ref12">12</xref>] carried out a study to develop and find the effect of using and combining four cycle design variables to reduce the specific fuel consumption (SFC) of a turbofan engine. In the study, the bypass ratio (B), the fan pressure ratio, the overall pressure ratio and turbine inlet temperature (TIT) were selected as the variables. The SFC was reduced without any constraint in the minimum thrust. Afterwards, a minimum specific thrust constraint was used. The study discovered that the first condition of no constraint in thrust resulted in a two-dimensional optimisation problem whereas when the thrust constraint was introduced, it became a three-dimensional problem. Sensitivity analysis was also carried out on the optimised cycle, results showed that there was little or no deviations in SFC from design values and that the SFC were more sensitivity to TIT and FRP. The study revealed that the By-pass ratio no longer the limiting factor, as it is now a three-dimensional (B, FPR, TIT) problem, however the overall pressure ratio remained the limiting factor in all the cases studied.</p><p>In order to study and overcome the shortcomings of energy loss caused by compression heating in compressed air energy storage technology, Xin He et al. [<xref ref-type="bibr" rid="scirp.119860-ref13">13</xref>] proposed a novel constant pressure pumped hydro system. Their study adopted an off-design model so as to make the study close to reality. The model was used to calculate and analyse the effect of key parameters on system thermodynamics performance</p><p>In the analysis carried out by Hiromasa Suzuki et al. [<xref ref-type="bibr" rid="scirp.119860-ref14">14</xref>] they opted for experimental analysis to study the effect of an under expanded jet radially discharging from a circular slit nozzle which consists of two cylinders, the aim was to study the frequency of the noise emitted from the jet when different nozzle pressure ratios and different cylinder diameters were used.</p><p>At higher altitudes, a wrong configuration of the operating parameters of a jet engine could result to a reduced optimal propulsive efficiency. This could cause the engine to consume more fuel than normal and produce a much lesser thrust than capable. Najjar and Al-Sharif [<xref ref-type="bibr" rid="scirp.119860-ref12">12</xref>] conducted a research study on the reduction of the specific fuel consumption on a turbofan engine while Najjar and Balawneh [<xref ref-type="bibr" rid="scirp.119860-ref9">9</xref>] focused on finding the optimal pressure ratio for various design variables. However, not much work has been done on determining these factors when using an afterburner. Thus, the aim of this study is to evaluate the performance of a turbojet engine with and without afterburner for different combinations of design choices and operating conditions and to determine an optimal configuration of the cycle for optimal operation.</p></sec><sec id="s2"><title>2. Methodology</title><p>The major objectives of this study are: 1) to model a turbojet engine with an afterburner and validate the model using realistic data 2) determine the influence of ambient temperature on performance of design point 3) optimize the cycle for minimum specific thrust with and without minimum thrust constraint, these are achieved by modelling the process using commercial softwares and analysing the results.</p><sec id="s2_1"><title>2.1. Modelling</title><p>Modelling of the engine was conducted in two parts: Analytical modelling using Microsoft Excel and Computer modelling using GasTurb 13.</p><sec id="s2_1_1"><title>2.1.1. Developing the Model in Microsoft Excel</title><p>Governing Equations</p><p>The different processes and Equations (1) and (7) in the engine are outlined below:</p><p>1) Intake/Diffuser</p><p>The total temperature at the diffuser outlet depends on the diffuser efficiency, η d .The outlet pressure and temperature will be given by</p><p>P 02 = P a ( 1 + η d γ c − 1 2 M a 2 ) γ c / γ c − 1 (1)</p><p>T 02 = T a ( 1 + γ c − 1 2 M a 2 ) (2)</p><p>2) Compressor</p><p>From the diagram in <xref ref-type="fig" rid="fig1">Figure 1</xref>(a), the state 2-3 represents the adiabatic compression process, the compressor isentropic efficiency, η c is a determinant factor. The outlet conditions are given by:</p><p>P 03 = P 02 π c (3)</p><p>T 03 = T 02 [ 1 + π c γ c − 1 / γ c − 1 η c ] (4)</p><p>3) Combustion Chamber</p><p>A reduction in stagnation pressure across the combustion chamber is expected due to the effect of fluid friction. The outlet pressure from the combustion chamber is defined as</p><p>P 04 = P 03 ( 1 − Δ P c c % ) (5)</p><p>The combustor outlet temperature is limited by the turbine metallurgical considerations. The fuel-to air ratio is defined as</p><p>f = C p g T 04 − C p a T 03 η b Q R − C p g T 04 (6)</p><p>4) Turbine</p><p>From the diagram in <xref ref-type="fig" rid="fig1">Figure 1</xref>(b), the state (4)-(5) defines the turbine expansion which accounts for the power required by the compressor and mechanical losses, thus the outlet temperature from the turbine is defined</p><p>W t = W c η m (7)</p><p>C p g ( T 04 − T 05 ) = C p a ( T 03 − T 02 ) / η m (8)</p><p>T 05 = T 04 − C p a ( T 03 − T 02 ) / η m C p g (9)</p><p>While its outlet pressure is expressed as</p><p>P 05 P 04 = [ 1 − 1 η t ( 1 − T 05 T 04 ) ] γ g − 1 / γ g (10)</p><p>5) Afterburner</p><p>a) Inoperative afterburner</p><p>A pressure drop within the afterburner is defined, the drop is due to presence of drag and skin friction in the flame holders. The expression is defined as</p><p>P 06 = P 05 − Δ P a b (11)</p><p>Since no afterburning is carried out, thus the temperature is constant in the afterburner duct.</p><p>T 06 = T 05 (12)</p><p>b) Operative afterburner</p><p>Additional fuel is burnt in the afterburner, this leads to a rise in the output temperature. The maximum temperature in the cycle is now given as</p><p>T 06 A = T max (13)</p><p>The afterburner fuel-to-air ratio is defined as</p><p>f a b = ( 1 + f ) ( C p g T 06 A − C p g T 05 ) η a b Q R − C p g T 06 A (14)</p><p>6) Nozzle</p><p>The critical pressure is used in determining if the nozzle is choked.</p><p>Inoperative afterburner</p><p>The critical pressure is defined as</p><p>P 06 P c = 1 [ 1 − 1 η n ( γ g − 1 γ g + 1 ) ] γ g − 1 / γ g (15)</p><p>where η n is the efficiency of the nozzle. For unchoked nozzles, the outlet pressure from the jet is equal to the ambient pressure, therefore the jet speed is defined as:</p><p>V 7 = 2 C p g ( T 06 − T 7 ) (16)</p><p>However, if the nozzle is choked, the outlet temperature is defined as</p><p>T 06 T 7 = γ h + 1 2 (17)</p><p>The jet speed is then defined as</p><p>V 7 = R γ g ( T 7 ) (18)</p><p>a) Operative afterburner</p><p>In the presence of an afterburner, the critical pressure is defined as</p><p>P 06 A P c = 1 [ 1 − 1 η n ( γ g − 1 γ g + 1 ) ] γ g − 1 / γ g (19)</p><p>For unchoked nozzles, the jet speed is defined as</p><p>V 7 a b = 2 C p g ( T 06 A − T 7 ) (20)</p><p>When it is choked, it is expressed as</p><p>V 7 a b = R γ g ( T 7 ) (21)</p><p>The specific thrust of the engine is therefore defined as</p><p>T m ˙ a = [ ( 1 + f + f a b ) V 7 − V ] + A 7 m ˙ a ( P 7 − P a ) (22)</p><p>While the thrust-specific fuel consumption (TSFC) is expressed as</p><p>TSFC = m ˙ f + m ˙ f a b T (23)</p><p>For an inoperative afterburner, the same equation was used but m ˙ f a b and f a b is equal to zero.</p></sec><sec id="s2_1_2"><title>2.1.2. Developing the Model in GasTurb13</title><p>Gasurb13 is incorporated with pre-designed jet engine models. The choice model for the research work to be done was selected and loaded. There was no need to input Equations because it has its own solver. GasTurb does not require user-defined specific heat ratios, which are necessary for the calculations rather it uses component maps for its design point calculation. GasTurb uses a pseudo-perfect gas model in which the specific heat ratios are a function of temperature and gas composition [<xref ref-type="bibr" rid="scirp.119860-ref15">15</xref>].</p></sec></sec><sec id="s2_2"><title>2.2. Validation of Model</title><p>The model developed in Excel and GasTurb13 was validated by live data of functioning engines gotten from Nate Meier’s jet engine (Meier, 2005b, 2005a) as shown in <xref ref-type="table" rid="table1">Table 1</xref>. These values were evaluated at sea-level static conditions.</p></sec><sec id="s2_3"><title>2.3. Assumptions and Design Parameters</title><p>The design parameters and assumed valves used in developing the model are presented in <xref ref-type="table" rid="table2">Table 2</xref>.</p></sec><sec id="s2_4"><title>2.4. On-Design Engine Performance</title><p>On-design performance analysis is defined at the flight conditions for which the engine is originally sized. The data used for the calculations are obtained at static sea level which means that the Mach number and altitude are zero. The performance analysis of the engine at the on-design point will be referred to as the reference point performance as it will constantly by referred to when carrying out more calculations at different atmospheric and flight conditions. These analyses will be carried out for both operative and inoperative afterburner. Sensitivity analysis is carried out on the on-design parameters to determine the effect of small changes in compressor pressure ratio and Turbine inlet temperature on thrust and specific fuel consumption. Further calculations was carried out using other component parameters with varying levels of technology between level 1 and level 3 as discussed in the book titled Elements of Gas propulsion textbook [<xref ref-type="bibr" rid="scirp.119860-ref14">14</xref>]. These values are presented below in <xref ref-type="table" rid="table3">Table 3</xref> and <xref ref-type="table" rid="table4">Table 4</xref>.</p></sec><sec id="s2_5"><title>2.5. Off-Design Performance Analysis</title><p>A turbojet engine is said to be operating at off-design conditions when it operates</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Validation models gotten from Nate Meier’s engine</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Model name</th><th align="center" valign="middle"  colspan="2"  >Thrust (KN)</th><th align="center" valign="middle"  colspan="2"  >TSFC</th><th align="center" valign="middle"  rowspan="2"  >Airflow</th><th align="center" valign="middle"  rowspan="2"  >Pressure Ratio</th></tr></thead><tr><td align="center" valign="middle" >Dry</td><td align="center" valign="middle" >Wet</td><td align="center" valign="middle" >Dry</td><td align="center" valign="middle" >Wet</td></tr><tr><td align="center" valign="middle" >SNECMA, Atar09C, Mirage 111E</td><td align="center" valign="middle" >41.95</td><td align="center" valign="middle" >60.802</td><td align="center" valign="middle" >0.1029</td><td align="center" valign="middle" >0.2068</td><td align="center" valign="middle" >68.03</td><td align="center" valign="middle" >5.7</td></tr><tr><td align="center" valign="middle" >Rolls-Royce, Avon RA.7R MK 114</td><td align="center" valign="middle" >33.406</td><td align="center" valign="middle" >42.258</td><td align="center" valign="middle" >0.0879</td><td align="center" valign="middle" >0.1938</td><td align="center" valign="middle" >58.51</td><td align="center" valign="middle" >7.0</td></tr><tr><td align="center" valign="middle" >Rolls-Royce, Avon RB.146 MK 301</td><td align="center" valign="middle" >53.378</td><td align="center" valign="middle" >69.903</td><td align="center" valign="middle" >0.0949</td><td align="center" valign="middle" >0.1888</td><td align="center" valign="middle" >77.11</td><td align="center" valign="middle" >8.4</td></tr><tr><td align="center" valign="middle" >GE, J47-GE-1, F-86A</td><td align="center" valign="middle" >21.573</td><td align="center" valign="middle" >25.888</td><td align="center" valign="middle" >0.112068</td><td align="center" valign="middle" >0.11716</td><td align="center" valign="middle" >40.82</td><td align="center" valign="middle" >4.3</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Custom engine parameter</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Parameter Type</th><th align="center" valign="middle"  colspan="2"  >Turbojet</th><th align="center" valign="middle" >References</th></tr></thead><tr><td align="center" valign="middle" >Parameter</td><td align="center" valign="middle" >Values</td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle"  rowspan="3"  >Ambient Parameters</td><td align="center" valign="middle" >Altitude</td><td align="center" valign="middle" >10,000 m</td><td align="center" valign="middle"  rowspan="3"  >Klein, [<xref ref-type="bibr" rid="scirp.119860-ref16">16</xref>]</td></tr><tr><td align="center" valign="middle" >Ambient Temperature</td><td align="center" valign="middle" >223.3 K</td></tr><tr><td align="center" valign="middle" >Ambient Pressure</td><td align="center" valign="middle" >0.2650 bar</td></tr><tr><td align="center" valign="middle"  rowspan="10"  >Parameters of the Gas Turbine</td><td align="center" valign="middle" >Mach number</td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle"  rowspan="10"  >Assumed for Design Point calculation</td></tr><tr><td align="center" valign="middle" >Compressor Pressure Ratio</td><td align="center" valign="middle" >9</td></tr><tr><td align="center" valign="middle" >Turbine Inlet Temperature</td><td align="center" valign="middle" >1200</td></tr><tr><td align="center" valign="middle" >Intake Efficiency</td><td align="center" valign="middle" >0.9</td></tr><tr><td align="center" valign="middle" >Compressor Efficiency</td><td align="center" valign="middle" >0.9</td></tr><tr><td align="center" valign="middle" >Combustion Efficiency</td><td align="center" valign="middle" >0.9</td></tr><tr><td align="center" valign="middle" >Propelling Nozzle Efficiency</td><td align="center" valign="middle" >0.9</td></tr><tr><td align="center" valign="middle" >Mechanical Efficiency</td><td align="center" valign="middle" >0.9</td></tr><tr><td align="center" valign="middle" >Combustion Pressure loss</td><td align="center" valign="middle" >6%</td></tr><tr><td align="center" valign="middle" >Nozzle inlet Temperature (Afterburner exit)</td><td align="center" valign="middle" >2000 K</td></tr><tr><td align="center" valign="middle"  rowspan="5"  >Other parameters</td><td align="center" valign="middle" >Specific Heat at constant Pressure (Air)</td><td align="center" valign="middle" >1.4</td><td align="center" valign="middle"  rowspan="4"  >Mattingly [<xref ref-type="bibr" rid="scirp.119860-ref17">17</xref>]</td></tr><tr><td align="center" valign="middle" >Specific Heat at constant Pressure (Combustion Gas)</td><td align="center" valign="middle" >1.157</td></tr><tr><td align="center" valign="middle" >Specific Heat Ratio of Combustion Gases</td><td align="center" valign="middle" >1.333</td></tr><tr><td align="center" valign="middle" >Gas constant</td><td align="center" valign="middle" >287 J/(KgK)</td></tr><tr><td align="center" valign="middle" >Fuel Calorific value</td><td align="center" valign="middle" >42,100 KJ/Kg</td><td align="center" valign="middle" >Klein [<xref ref-type="bibr" rid="scirp.119860-ref16">16</xref>]</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> On-Design parameters</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Mass flow (Kg/s)</th><th align="center" valign="middle" >30</th></tr></thead><tr><td align="center" valign="middle" >Pressure ratio</td><td align="center" valign="middle" >9</td></tr><tr><td align="center" valign="middle" >TIT (K)</td><td align="center" valign="middle" >1200</td></tr><tr><td align="center" valign="middle" >Afterburner Temperature (K)</td><td align="center" valign="middle" >1800</td></tr><tr><td align="center" valign="middle" >Altitude (m)</td><td align="center" valign="middle" >0</td></tr><tr><td align="center" valign="middle" >Mach Number</td><td align="center" valign="middle" >0</td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> On-design parameters based on level of technology</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >LEVEL 1</th><th align="center" valign="middle" >LEVEL 2</th><th align="center" valign="middle" >LEVEL 3</th></tr></thead><tr><td align="center" valign="middle" >π<sub>c</sub></td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >0.84</td><td align="center" valign="middle" >0.88</td></tr><tr><td align="center" valign="middle" >TIT</td><td align="center" valign="middle" >1110</td><td align="center" valign="middle" >1390</td><td align="center" valign="middle" >1780</td></tr><tr><td align="center" valign="middle" >π<sub>t</sub></td><td align="center" valign="middle" >0.9</td><td align="center" valign="middle" >0.92</td><td align="center" valign="middle" >0.94</td></tr><tr><td align="center" valign="middle" >η<sub>b</sub></td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >0.85</td><td align="center" valign="middle" >0.89</td></tr></tbody></table></table-wrap><p>at atmospheric and flight conditions different from the prescribed on-design conditions. The performance projections at off design is expected to different from those at on-design conditions since the turbojet engine’s geometry is designed for on-design performance. In carrying out off-design analysis, the following assumptions are made:</p><p>1) The gas is ideal and perfect;</p><p>2) The Flow at the turbine inlet and nozzle throat is choked;</p><p>3) The pressure ratio is constant;</p><p>4) Efficiencies of the various components are constant;</p><p>5) All the Power developed in the turbine is used in powering the compressor;</p><p>6) The Fuel-air ratio is constant.</p><p>The formulations for the Off-design analysis are similar to that for the on-design calculations, except for the compressor ratio, temperature values and mass flow rates. Off-design calculations were carried out for both varying Mach numbers and altitude to determine the effect of ambient temperature, pressure and flight speed on engine performance.</p></sec><sec id="s2_6"><title>2.6. Optimization of Cycle</title><p>In the design of a turbojet engine, the thermodynamic variables of Turbine Inlet Temperature (TIT) and Pressure ratio,rp, are used in determining the optimum performance of the engine. Performance optimization is intended to determine under what conditions the maximum Thrust (F) and minimum SFC of the engine can be obtained. This is carried out by determining the optimum pressure ratio (rp) and the corresponding turbine inlet temperature (TIT) for which the turbojet engine performance gives maximum thrust and minimum SFC.</p><p>The optimisation studies are carried out in two modules: 1) The SFC is to be minimized without any specific thrust constraints. The study is carried out for different operating conditions and maximum TIT as indicated in Table. 2) The SFC is minimized with specific thrust constraints, in which for each TIT there is pressure ratio that gives maximum thrust. Afterwards, the analytical model developed with Microsoft Excel was used in determining the optimum running line (ORL) over a range of operating conditions.</p></sec></sec><sec id="s3"><title>3. Results and Discussion</title><p>The values were used in the analytical calculations of the turbojet engine parameters.</p><sec id="s3_1"><title>3.1. Results Validation</title><p>The validation simulations are obtained at sea level static conditions that are zero altitude and Mach number. At this level, the design point thrust is maximum. The difference in values of analytical calculations and the GasTurb simulation as shown in <xref ref-type="table" rid="table5">Table 5</xref> is due to the configuration of the software as described in Section 2.1.2. “Dry” result represents the performance of the engine without</p><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Validation results</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  ></th><th align="center" valign="middle"  colspan="3"  >Dry</th><th align="center" valign="middle"  colspan="2"  >Wet</th></tr></thead><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >Analytical</td><td align="center" valign="middle" >GasTurb13</td><td align="center" valign="middle" >Analytical</td><td align="center" valign="middle" >GasTurb13</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Model 1</td><td align="center" valign="middle" >Thrust</td><td align="center" valign="middle" >41.94073</td><td align="center" valign="middle" >41.880</td><td align="center" valign="middle" >60.802</td><td align="center" valign="middle" >60.849</td></tr><tr><td align="center" valign="middle" >SFC</td><td align="center" valign="middle" >0.1029</td><td align="center" valign="middle" >0.0899</td><td align="center" valign="middle" >0.2068</td><td align="center" valign="middle" >0.18532</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Model 2</td><td align="center" valign="middle" >Thrust</td><td align="center" valign="middle" >33.406</td><td align="center" valign="middle" >33.427</td><td align="center" valign="middle" >42.258</td><td align="center" valign="middle" >42.373</td></tr><tr><td align="center" valign="middle" >SFC</td><td align="center" valign="middle" >0.0879</td><td align="center" valign="middle" >0.07963</td><td align="center" valign="middle" >0.1938</td><td align="center" valign="middle" >0.1350</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Model 3</td><td align="center" valign="middle" >Thrust</td><td align="center" valign="middle" >53.378</td><td align="center" valign="middle" >53.803</td><td align="center" valign="middle" >69.904</td><td align="center" valign="middle" >70.08</td></tr><tr><td align="center" valign="middle" >SFC</td><td align="center" valign="middle" >0.0949</td><td align="center" valign="middle" >0.08629</td><td align="center" valign="middle" >0.1888</td><td align="center" valign="middle" >0.1888</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Model 4</td><td align="center" valign="middle" >Thrust</td><td align="center" valign="middle" >21.576</td><td align="center" valign="middle" >21.72</td><td align="center" valign="middle" >25.888</td><td align="center" valign="middle" >26.02</td></tr><tr><td align="center" valign="middle" >SFC</td><td align="center" valign="middle" >0.112068</td><td align="center" valign="middle" >0.0924</td><td align="center" valign="middle" >0.11716</td><td align="center" valign="middle" >0.14317</td></tr></tbody></table></table-wrap><p>afterburner while the “wet” result represents the performance of the engine with afterburner effects.</p></sec><sec id="s3_2"><title>3.2. On-Design Simulation Comparison</title><p>Components maps are not used for the simulations, the difference in values between the analytical solution and that from the model as shown in <xref ref-type="table" rid="table6">Table 6</xref> are caused by the specific heat variations.</p><p>The T-S diagrams for the design point analysis without afterburning and with after burning are shown in <xref ref-type="fig" rid="fig2">Figure 2</xref> and <xref ref-type="fig" rid="fig3">Figure 3</xref> respectively.</p>Sensitivity Analysis<p>The effect of changes in burner exit temperature, reheat exit temperature and pressure ratio on the net thrust and specific fuel consumption are presented in <xref ref-type="table" rid="table7">Table 7</xref>.</p></sec><sec id="s3_3"><title>3.3. Off-Design Simulation/Parametric Study</title><p>Off design parameters are calculated in GasTurb13 using component maps. The custom engine parameter was used as a base of reference for the off-design simulation.</p><sec id="s3_3_1"><title>3.3.1. Effect of Varying Altitude at a Mach Number of 0.8</title><p>As seen in <xref ref-type="fig" rid="fig4">Figure 4</xref>, increasing altitude results in reduction in the specific fuel consumption, the same trend occurs in the fuel flow and the thrust at a constant Mach number. It is noticed that the engine’s maximum attainable thrust occurs at sea level. The results are similar with and without an afterburner.</p><p>The thrust also decreased as seen in <xref ref-type="fig" rid="fig4">Figure 4</xref> as the altitude increases from 0 to 10,000 m. This is due to the increased drag on the jet as it moves through the air. Definitely, there is an increase in fuel consumption as seen in <xref ref-type="fig" rid="fig5">Figure 5</xref>. The relationship between the altitude and fuel flow is presented in <xref ref-type="fig" rid="fig6">Figure 6</xref>, the higher the altitude, the less fuel is consumed.</p><table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> Results from on-design analysis</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="3"  >Dry</th><th align="center" valign="middle"  colspan="2"  >Wet</th></tr></thead><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >Analytical</td><td align="center" valign="middle" >GasTurb13</td><td align="center" valign="middle" >Analytical</td><td align="center" valign="middle" >GasTurb13</td></tr><tr><td align="center" valign="middle" >Thrust</td><td align="center" valign="middle" >22.287</td><td align="center" valign="middle" >19.78</td><td align="center" valign="middle" >33.057</td><td align="center" valign="middle" >28.40</td></tr><tr><td align="center" valign="middle" >SFC</td><td align="center" valign="middle" >0.1177</td><td align="center" valign="middle" >0.09016</td><td align="center" valign="middle" >0.1847</td><td align="center" valign="middle" >0.18123</td></tr></tbody></table></table-wrap><table-wrap id="table7" ><label><xref ref-type="table" rid="table7">Table 7</xref></label><caption><title> Sensitivity analysis table</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  ></th><th align="center" valign="middle" ></th><th align="center" valign="middle" ></th><th align="center" valign="middle" ></th><th align="center" valign="middle" >FN</th><th align="center" valign="middle" >SFC</th><th align="center" valign="middle" >WF Total</th></tr></thead><tr><td align="center" valign="middle" >Unit</td><td align="center" valign="middle" >Basis</td><td align="center" valign="middle" >Delta</td><td align="center" valign="middle" >%</td><td align="center" valign="middle" >%</td><td align="center" valign="middle" >%</td></tr><tr><td align="center" valign="middle" >Burner Exit Temperature</td><td align="center" valign="middle" >K</td><td align="center" valign="middle" >1200</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >+0.43</td><td align="center" valign="middle" >−0.51</td><td align="center" valign="middle" >−0.08</td></tr><tr><td align="center" valign="middle" >Reheat Exit Temperature</td><td align="center" valign="middle" >K</td><td align="center" valign="middle" >1800</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >+0.30</td><td align="center" valign="middle" >+0.56</td><td align="center" valign="middle" >+0.87</td></tr><tr><td align="center" valign="middle" >Pressure Ratio</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >9</td><td align="center" valign="middle" >0.2</td><td align="center" valign="middle" >+0.26</td><td align="center" valign="middle" >−0.23</td><td align="center" valign="middle" >+0.02</td></tr></tbody></table></table-wrap></sec><sec id="s3_3_2"><title>3.3.2. Effect of Varying Mach Number at an Altitude of 10,000 m</title><p>As seen in <xref ref-type="fig" rid="fig7">Figure 7</xref>, increasing the Mach number leads to an increase in total fuel flow in the turbojet engine this implies increased specific fuel consumption at altitude of 10,000 m as the Mach number increases. This is also the trend in the net thrust as it increases with an increase in Mach number.</p><p>The thrust also increases as seen in <xref ref-type="fig" rid="fig8">Figure 8</xref> as the Mach number increases from 0 to 2. This is due to the increased drag on the jet as it moves through the air. Definitely, there is an increase in fuel consumption as seen in <xref ref-type="fig" rid="fig9">Figure 9</xref>.</p></sec><sec id="s3_3_3"><title>3.3.3. Performance Carpet at Variable Pressure Ratio, Turbine Inlet Temperature and Afterburner Exit Temperature</title><p>This was calculated at design point and at custom engine parameters. The values obtained show the effect of pressure ratio and thrust on specific fuel consumption and Net Thrust. The carpet performance of the engine at static sea levels for the burner exit temperature conditions and reheat exit temperature conditions are presented in <xref ref-type="fig" rid="fig1">Figure 1</xref>0 and <xref ref-type="fig" rid="fig1">Figure 1</xref>1 respectively, while the carpet performance when at altitude 10,000 m and Mach number 0.8 for both burner exit and reheat exit conditions are presented in <xref ref-type="fig" rid="fig1">Figure 1</xref>2 and <xref ref-type="fig" rid="fig1">Figure 1</xref>3 respectively.</p></sec></sec><sec id="s3_4"><title>3.4. Optimum Running Line</title><p>An optimum point line was superimposed on the compressor map as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>4. The point at which the specific fuel consumption is lowest pressure ratio is minimum is called the cruise mode. Since there is no compressor surge, there was no need to introduce a control bleed to counter the surge. The relationship between the fuel consumption and net thrust is presented in <xref ref-type="fig" rid="fig1">Figure 1</xref>5.</p></sec><sec id="s3_5"><title>3.5. Optimization Results</title><sec id="s3_5_1"><title>3.5.1. Minimizing Specific Fuel Consumption without Net Thrust Constraint</title><p>The results of the optimization simulations when there was no net thrust constraint are summarized in <xref ref-type="table" rid="table8">Table 8</xref>. The values were evaluated between pressure</p><table-wrap id="table8" ><label><xref ref-type="table" rid="table8">Table 8</xref></label><caption><title> Optimisation results at different operating conditions</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="4"  >Operating Conditions</th><th align="center" valign="middle"  colspan="3"  >Optimum conditions</th></tr></thead><tr><td align="center" valign="middle" >Altitude (m)</td><td align="center" valign="middle" >Mach Number</td><td align="center" valign="middle" >TIT (K)</td><td align="center" valign="middle" >Reheat Temperature (K)</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >Thrust (KN)</td><td align="center" valign="middle" >SFC Kg/N hr</td></tr><tr><td align="center" valign="middle" >0 m</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >1450</td><td align="center" valign="middle" >1800</td><td align="center" valign="middle" >16.8504</td><td align="center" valign="middle" >7.10</td><td align="center" valign="middle" >0.1751</td></tr><tr><td align="center" valign="middle" >10,000 m</td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >1450</td><td align="center" valign="middle" >1800</td><td align="center" valign="middle" >19.6336</td><td align="center" valign="middle" >11.58</td><td align="center" valign="middle" >0.0991</td></tr></tbody></table></table-wrap><table-wrap id="table9" ><label><xref ref-type="table" rid="table9">Table 9</xref></label><caption><title> Optimisation results at different operating conditions</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="4"  >Operating Conditions</th><th align="center" valign="middle"  colspan="3"  >Optimum conditions</th></tr></thead><tr><td align="center" valign="middle" >Altitude (m)</td><td align="center" valign="middle" >Mach Number</td><td align="center" valign="middle" >TIT (K)</td><td align="center" valign="middle" >Reheat Temperature (K)</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >Pressure Ratio</td><td align="center" valign="middle" >SFC Kg/N hr</td></tr><tr><td align="center" valign="middle" >0 m</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >1450</td><td align="center" valign="middle" >1800</td><td align="center" valign="middle" >16.8504</td><td align="center" valign="middle" >7.8</td><td align="center" valign="middle" >0.1538</td></tr><tr><td align="center" valign="middle" >10,000 m</td><td align="center" valign="middle" >0.8</td><td align="center" valign="middle" >1450</td><td align="center" valign="middle" >1800</td><td align="center" valign="middle" >19.6336</td><td align="center" valign="middle" >16.78</td><td align="center" valign="middle" >0.0887</td></tr></tbody></table></table-wrap><p>ratio of 5 and 20 and with a maximum turbine inlet temperature constraint of 1450 K.</p></sec><sec id="s3_5_2"><title>3.5.2. Minimizing Specific Fuel Consumption with a Constraint on Net Thrust</title><p>The effect of constraining the thrust is shown in <xref ref-type="table" rid="table9">Table 9</xref>. A constraint of thrust at 15 KN was introduced and the optimal pressure ratio at which there is minimal fuel consumption was calculated. The pressure ratio at a maximum Thrust of 15 KN was determined to be 16.78.</p></sec></sec></sec><sec id="s4"><title>4. Conclusion</title><p>This study has shown that the altitude and Mach numbers are limiting factors when the specific fuel consumption is minimised without constraining the minimum thrust thereby resulting in a two-dimensional optimisation problem. Results from this work show that TSFC reduces as the Altitude and Mach number increases at optimised turbine inlet temperature and compressor pressure ratio. Practical considerations, however, limit the potential improvements. The significant trends observed in this study are: 1) At increased Altitude and Mach number the TIT increases; 2) Pressure ratio decreases at increasing Altitude; 3) Significant reduction in Thrust at increasing; 4) Presence of an afterburner improved the thrust but also increases the Specific fuel consumption. The sensitivity study of the on-design cycle showed that there is little or no effect on the TSFC due to small deviations from optimum design values. It further showed that the TSFC of the optimised cycle is quite sensitive to variations in TIT and pressure ratio. However, the overall pressure ratio remains the limiting factor in all the cases studied</p></sec><sec id="s5"><title>Declaration of Competing Interest</title><p>The authors declared that there is no conflict of interest.</p></sec><sec id="s6"><title>Cite this paper</title><p>Ighodaro, O.O., Aburime, E.I. and Erameh, A.A. (2022) Off-Design Modelling of a Turbo Jet Engine with Operative Afterburner. Open Journal of Energy Efficiency, 11, 88-107. https://doi.org/10.4236/ojee.2022.113007</p></sec><sec id="s7"><title>Nomenclature</title><p>P: Pressure</p><p>T: Temperature</p><p>Ma: Mach Number</p><p>η: Isentropic Efficiency</p><p>F: Fuel to air ratio</p><p>Cp: Specific heat capacity</p><p>W: Power</p><p>TIT: Turbine inlet temperature</p></sec></body><back><ref-list><title>References</title><ref id="scirp.119860-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Farokhi, S. (2012) Aircraft Propulsion. 2nd Edition, John Wiley and Sons, Kansas.</mixed-citation></ref><ref id="scirp.119860-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Gobran, M.H. (2013) Off-Design Performance of Solar Centaur-40 Gasturbine Engine Using Simulink. Ain Shams Engineering Journal, 4, 285-298.  
https://doi.org/10.1016/j.asej.2012.08.007</mixed-citation></ref><ref id="scirp.119860-ref3"><label>3</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Lazzarettoand</surname><given-names> T. </given-names></name>,<etal>et al</etal>. (<year>2001</year>)<article-title>Analytical and Neural Network Models for Gas Turbine Design and Off-Design Simulation</article-title><source> International Journal of Applied Thermodynamics</source><volume> 4</volume>,<fpage> 173</fpage>-<lpage>182</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.119860-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Shpiro, S.R. and Caddy, M.J. (1974) NECOMP: The Navy Engine Performance Program. Paper No. 74-GT-83, American Society of Mechanical Engineers, New York.  
https://doi.org/10.1115/74-GT-83</mixed-citation></ref><ref id="scirp.119860-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Wittenberg, H. (1978) Prediction of Off-Design Performance of Turbojet and Turbofan Engines. AGARD Conference Proceedings, 242, 4.1-4.31.</mixed-citation></ref><ref id="scirp.119860-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Sellers, J.J. and Daniels, C.J. (1975) DYNGEN—A Program for Calculating Steady State and Transient Performance of Turbojet and Turbofan Engines. NASA TN D-7901, National Aeronautics and Space Administration (NASA), Washington DC.</mixed-citation></ref><ref id="scirp.119860-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Yang, Y.H. (1991) A New Method of Predicting the Performance of Gas Turbine Engines. Journal of Engineering for Gas Turbines and Power, 113, 106-111.  
https://doi.org/10.1115/1.2906516</mixed-citation></ref><ref id="scirp.119860-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Zhu, P. and Saravanamutto, H.I.H. (1992) Simulation of an Advanced Twin Spool Industrial Gas Turbine. Journal of Engineering for Gas Turbines and Power, 114, 180-186. https://doi.org/10.1115/1.2906568</mixed-citation></ref><ref id="scirp.119860-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Najjar, Y. and Balawneh, I. (2015) Optimization of Gas Turbines for Sustainable Turbojet Propulsion. Propulsion and Power Research, 4, 114-121.  
https://doi.org/10.1016/j.jppr.2015.05.004</mixed-citation></ref><ref id="scirp.119860-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Bakalis, D.P. and Stamatis, A.G. (2012) Data Analysis and Performance Model Calibration of a Small Turbojet Engine. Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, 226, 1523-1533.  
https://doi.org/10.1177/0954410011425126</mixed-citation></ref><ref id="scirp.119860-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Wood, A. and Pilidis, P. (1997) A Variable Cycle Jet Engine for ASTOVL Aircraft. Aircraft Engineering and Aerospace Technology, 69, 534-539.  
https://doi.org/10.1108/00022669710186002</mixed-citation></ref><ref id="scirp.119860-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Najjar, Y.S.H. and Al-Sharif, S.F. (2006) Thermodynamic Optimization of the Turbofan Cycle. Aircraft Engineering and Aerospace Technology, 78, 467-480.  
https://doi.org/10.1108/00022660610707139</mixed-citation></ref><ref id="scirp.119860-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">He, X., Wang, H., Tao, F. and Ge, G. (2021) Performance Analysis of Constant-Pressure Pumped Hydro Combined with Compressed Air Energy Storage System Considering Off-Design Model of Compressor. Energy and Power Engineering, 13, 11-18. https://doi.org/10.4236/epe.2021.134B002</mixed-citation></ref><ref id="scirp.119860-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Suzuki, H., Kawasaki, K., Endo, M. and Sakakibara, Y. (2020) Structure and Screech Tone of Radially Underexpanded Jet Emitted by Facing Cylinders. Open Journal of Fluid Dynamics, 10, 52-62. https://doi.org/10.4236/ojfd.2020.101004</mixed-citation></ref><ref id="scirp.119860-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Sung, R. (2013) A Comparative Study of the Gas Turbine Simulation Program (GSP) 11 and GasTurb 11 on Their Respective Simulations for a Single-Spool Turbojet. Master’s Thesis, University of Tennessee, Knoxville.</mixed-citation></ref><ref id="scirp.119860-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Klein, D. and Abeykoon, C. (2015) Modelling of a Turbojet Gas Turbine Engine. 2015 Internet Technologies and Applications, Wrexham, 8-11 September 2015, 200-206.  
https://doi.org/10.1109/ITechA.2015.7317395</mixed-citation></ref><ref id="scirp.119860-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Mattingly, J.D. (2010) Elements of Gas Turbine Propulsion. 1st Edition, TATA McGraw-Hill, New Dehli.</mixed-citation></ref></ref-list></back></article>