Gas Turbine,finned Heat Exchanger

  • June 2020
  • PDF

This document was uploaded by user and they confirmed that they have the permission to share it. If you are author or own the copyright of this book, please report to us by using this DMCA report form. Report DMCA


Overview

Download & View Gas Turbine,finned Heat Exchanger as PDF for free.

More details

  • Words: 25,445
  • Pages: 162
Republic of Iraq Ministry of Higher Education and Scientific Research Al- Mustansiriya University College of Engineering Mechanical Engineering Department

The Enhancement of Two- Shaft Gas Turbine Performance Using Improved Air Temperature

A THESIS SUBMMITTED TO THE COLLEGE OF ENGINEERING OF AL- MUSTANSIRIYA UNIVERSITY AS PARTIAL FULLFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER IN MECHANICAL ENGINEERING SCIENECE

BY

Ali Ahmed Abdulrasool

Supervised by Dr.Ali Hussain Tarrad

Asst.Prof Dr. Fouad Alwan Saleh

2009

‫ﻦ اﻟ َّﺮﺣِﻴ ِﻢ‬ ‫ِﺑﺴْ ِﻢ اﻟَّﻠ ِﻪ اﻟ َّﺮﺣْ َﻤ ِ‬

‫ن‬ ‫ب وَاﻟْﻤِﻴﺰَا َ‬ ‫ت وَأَﻧْﺰَﻟْﻨَﺎ َﻣ َﻌ ُﻬ ُﻢ اﻟْ ِﻜﺘَﺎ َ‬ ‫ﺳ َﻠﻨَﺎ ﺑِﺎﻟْ َﺒ ِّﻴﻨَﺎ ِ‬ ‫َﻟ َﻘﺪْ أَرْﺳَﻠْﻨَﺎ ُر ُ‬ ‫ﺤﺪِﻳ َﺪ ﻓِﻴ ِﻪ ﺑَﺄْسٌ ﺷَﺪِﻳﺪٌ َو َﻣﻨَﺎ ِﻓ ُﻊ‬ ‫ﻂ وَأَﻧْﺰَﻟْﻨَﺎ اﻟْ َ‬ ‫س ﺑِﺎﻟْ ِﻘﺴْ ِ‬ ‫ِﻟ َﻴﻘُﻮ َم اﻟ َﻨّﺎ ُ‬ ‫ن اﻟَّﻠ َﻪ َﻗ ِﻮيٌّ‬ ‫ﺐ ِإ َّ‬ ‫ﺳ َﻠ ُﻪ ﺑِﺎﻟْ َﻐﻴْ ِ‬ ‫ﺼ ُﺮ ُﻩ َو ُر ُ‬ ‫س َو ِﻟ َﻴﻌْ َﻠ َﻢ اﻟَّﻠ ُﻪ َﻣﻦْ َﻳﻨْ ُ‬ ‫ﻟِﻠ َﻨّﺎ ِ‬ ‫ﻋَﺰِﻳﺰٌ‬ ‫ﺻﺪق اﷲ اﻟﻌﻈﻴﻢ‬

‫ﺳﻮرة اﻟﺤﺪﻳﺪ اﻵﻳﺔ )‪(25‬‬

Supervisor Certificate We certify that the preparation of this thesis entitled

(THE

ENHANCEMENT OF TWO SHAFT GAS TURBINE PERFORMANCE USING IMPROVED AIR TEMPERATURE)

was made under our

supervision at the Mechanical Engineering Department , College of Engineering , Al-Mustansiriya University in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE IN MECHANICAL ENGINEERING. (SUPERVISORS) Signature:

Signature:

Name:

Name:

Date:

Date:

In view of the available recommendation I forward this thesis for debate by the examining committee Signature: Name: Date: (CHAIRMAN OF MECHANICAL ENGINEEING DEPARTMENT)

EXAMINATION COMMITTEE CERTIFICATE We certify that we have read the thesis entitled (THE ENHANCEMENT OF TWO SHAFT GAS TURBINE PERFORMANCE USING IMPROVED AIR TEMPERTURE) and as an examining committee, examined the student (Ali Ahmed Abdulrasool) in its context and that in our opinion it is adequate as a thesis for the degree of MASTER OF SCIENCE IN MECHANICAL ENGINEERING. Signature: Name: Date: (Chairman) Signature:

Signature:

Name:

Name:

Date:

Date:

(Member)

(Member)

Signature:

Signature:

Name:

Name:

Date:

Date:

(Supervisor)

(Supervisor)

Approved by the Dean of College of Engineering Signature: Name: Date:

DEDICATION TO MY DEAR FAMILY WITH LOVE AND RESPECT

ALI

Acknowledgement

All praise is due to Allah, the lord of the worlds

I would like to express my deep sense of appreciation and thank to my supervisors (Dr.Ali Hussain Tarrad and Dr. Fouad Alwan Saleh) for their continuous encouragement and support during the study. Finally, I would like to thank for all of people who helped me and introduce their opinion and companion to end this thesis including: Mr. Ahmed Muneer Dein, Mr. zyaad Talal ,Mr. Wa'al Najam and Mr. Nehaad Hashim.

Ali

ABSTRACT In the present research, a study of the performance improvement of a twoshaft gas turbine engine is conducted. This has been accomplished by utilizing sensible load cooling system. It is related to the effect of reducing compressor inlet temperature on gas turbine performance. An experimental study has been done on an existing a two-shaft gas turbine with a (5 kW) power, taking into consideration the effect of reducing compressor inlet temperature. An instrumented experimental rig was built for this object by adding an aircooled heat exchanger in series with water supplying system. The experimental results for gas turbine performance showed that the percent of design increases by (15%) for the power output increasing overall efficiency by (25%).Moreover, heat consumption has been reduced by (10%) when reducing compressor inlet temperature from (30°C to 15 °C).The percent of design is defined by the ratio between the parameter difference at both temperatures to the design point. A simplified new numerical model based on the step by step technique has been developed for the design and predicting the air cooled heat exchanger performance. The numerical model was designed in a new form so that variation of all design parameters can be calculated. The model has been checked and validated using the experimental laboratory data. The maximum discrepancy between the experimental data and model predicted values of the overall heat transfer coefficient and heat load were about (5%). This percentage value was obtained for the given range of the simulated conditions. Furthermore, a computational model program to investigate theoretically the effect of compressor inlet temperature on the performance of the gas turbine has been developed. The model utilizes a non-dimensional approach. The discrepancy percentage between the experimental and theoretical predictions was about (18 %).

TABLE OF CONTENTS Title

Pages

List of Figures

VIII

List of Tables

XII

Nomenclature

XIII

CHAPTER ONE:

Introduction

1

1.1 General

2

1.2 Two Shaft Engine

2

1.3 Influence of External Factors on the Gas Turbine Performance

3

1.4 Gas Turbine Inlet Air Cooling

5

1.4.1 Evaporative Cooling

5

1.4.2 Cooling with Absorption Chiller

5

1.5 Compact Heat Exchanger

6

1.6 Aims of the Present Work

7

CHAPTER TWO :

8

Literature Survey

2.1 General

9

2.2 The Previous Work

9

2.2.1 Modeling and Simulation of Air Cooled Heat Exchanger 2.2.2 Performance Improvement of the Gas Turbine Engine

9 11

2.3 The Present Work

16

CHAPTER THREE: Experimental Work

18

3.1 Test Rig for Heat Exchanger

19

3.1.1 Heat Exchanger

19

3.1.2 Water Supplying System

19

3.1.3 Measuring Instrumentation

21

3.1.4 The Electrical Board

21

3.1.5 Air Circulation System

22

3.2 Experimental Setup

26

3.3 Gas Turbine

26

3.3.1 Design Information

28

3.3.2 Performance design

28

3.4 Operating Principles

28

3.5 Main Component Parts of Gas Turbine

29

3.5.1 Gas Generator:

29

3.5.1.1 Centrifugal Compressor

29

3.5.1.2 Combustion Chamber

30

3.5.1.3 Gas Generator Turbine

31

3.5.2 Power Turbine

31

3.5.3 Fuel System

32

3.5.4 Oil Lubricating System

32

3.5.5 Starting System

33

3.6 Dimensionless and Parameter Groups

33

3.6.1 Corrected Compressor Data

34

3.6.2 Corrected Basic Data

35

3.6. 3 Corrected Derived Data

36

3.7 Experimental Work Procedure

40

3.8 Test Procedure

40

CHAPTER FOURE : Mathematical Model

43

4.1 General

44

4.2 Heat Exchanger Thermal Design

44

4.2.1 A Comprehensive Design Procedure

44

4.2.2 Numerical Modeling of Cross Flow Compact Heat Exchanger

46

4.2.2.1 Grid System

46

4.2.2.2 Physical Characteristics of Heat Exchanger

49

4.2.2.3 Mass Conservation

49

4.2.2.4 Log-Mean Temperature Difference

50

4.2.2.5 Heat Load

51

4.2.2.6 Overall Heat – Transfer Coefficient (Uo)

51

4.2.2.7 Forced Convection Heat Transfer coefficient Inside Tube

52

4.2.2.8 Forced Convection Heat Transfer coefficient for Air Side

52

4.2.2.9 Power of Fan

54

4-2-2-10 The Computer Program

55

4.3 Gas Turbine

57

4.4 Basic Gas Turbine Cycles

58

4.5 A non-Dimensional Analyses of Gas Turbine

58

Performance 4.5.1 Component Performance

59

4.5.2 Graphical Plot

61

4.6 Computer Calculations for Two-Shaft Gas Turbine

61

4.7 Results and Discussion of Theoretical Calculations

62

4.7.1 The effect of compressor inlet temperature on maximum to

62

minimum ratio (Ø) 4.7.2 The effect of compressor inlet temperature on the expansion

64

ratio CHAPTER FIVE : Results & Discussion

70

5.1 General

71

5.2 Computational Model Results for Heat Exchanger

71

5.2.1 Heat Load for heat exchanger

72

5.2.2 Heat Transfer Coefficient for air side

72

5.2.3 Overall Heat Transfer Coefficient

73

5.2.4 Air Temperature Distribution

74

5.2.5 The Effect of Core Aspect Ratio (H/L) and Size (LXDXH)

74

5.3 Experimental Results for Gas Turbine

75

5.3.1 Effect of Compressor Inlet Temperature on Power Output

76

5.3.2 Effect of Compressor Inlet Temperature on Fuel Mass Flow

76

Rate 5.3.3 Effect of Compressor Inlet Temperature on Specific Fuel

77

Consumption 5.3.4 Effect of Compressor Inlet Temperature on Heat

77

Consumption 5.3.5 Effect of Compressor Inlet Temperature on Heat Rate

78

5.3.6 Effect of Compressor Inlet Temperature on Overall

78

Efficiency 5.3.7 Effect of Compressor Inlet Temperature on Air Mass Flow

78

Rate 5.3.8 Effect of Compressor Inlet Temperature on Pressure Ratio

78

5.3.9 Effect of Compressor Inlet Temperature on Power Input to

79

Compressor 5.3.10 Effect of Compressor Inlet Temperature on Turbine Inlet

79

Temperature 5.4 Percent of Design

80

5-5 Comparison between the Experimental and Theoretical

81

Predictions of the gas turbine engine 5.6 Conclusion

81

CHAPTER SIX : Conclusions & Recommendations

114

6.1 Conclusions

115

- 6.2 Recommendations

116

References

117

Appendix (A)Experimental Work ,Data Tables

122

Appendix (B) Flow Charts and Computer Program

134

Appendix (C)Gas Turbine parameter Groups

146

List of Figures Title

Pages

Figure (1.1) Two Shaft Gas Turbine Diagram

3

Figure (1.2) Influences of External Factors on the Gas Turbine Performance

4

Figure (1.3) Evaporative Cooler

5

Figure (1.4 )Air "Chilling Cooling" System, Based on Absorption

6

Figure (2.1)Effect of Compressor Inlet Temperature on Gas Turbine Performance

13

Figure (2.2) Temperature-Power-Speed Interrelationships

15

Figure 2.3 Net Output Power Versus Inlet Temperature for Gas Cycle

16

Figure (3.1a) Schematic Diagram of the Built Rig (Heat Exchanger)

23

Figure (3.1b) Configuration of the Built Rig (Heat Exchanger)

24

Figure (3.1c) Configuration of the Test Rig (Inlet Cooling System) Preparing to Gas Turbine Figure (3.3.a) Top View of Heat Exchanger Geometry

24

Figure (3.3.b) Front View of Heat Exchanger Geometry

25

Figure (3.3) Schematic Diagram of the Test Rig (Overall)

27

Figure.(3.4) Schematic Arrangement for Gas Turbine (GT-85)

32

Figure (4.1) Methodology of Heat Exchanger Design

45

Figure (4.2) Step by Step method with two Directions

47

Figure (4.3a) Slice for inlet single Tube

48

Figure (4.3b) Exit of one row inlet to next row

48

Figure (4.3c) Nodal Points Distributions with two directions

49

Figure (4.4) The Mean Temperature Difference Along a Single Pass

50

Figure (4.5) Basic Gas Turbine Engine

59

Figure (4.6 )(T-s) Diagram for Irreversible Two-Shaft Circuit Simple Plant

56

Figure (4.7)Specific Heats and Their Ratios for ‘Real’ Gases-Air and Products of Combustion

62

Figure (4.8) Theoretical Overall Efficiency as a Function of (T3/T1 Ratio) with Isentropic Efficiency (ηt,ηc =0.9) ,T3=699 0C

65

Figure (4.9) Theoretical Overall Efficiency as a Function of (T3/T1 Ratio) with Isentropic Efficiency (ηt,ηc =0.9) ,T3=628.5 0C

65

25

Figure (4.10) Theoretical Overall Efficiency as a Function of (T3/T1 Ratio) with Isentropic Efficiency (ηt,ηc =0.9) ,T3=588 0C

66

Figure (4.11) Theoretical Overall Efficiency as a Function of (T3/T1 Ratio) with Isentropic Efficiency (ηt,ηc =0.9) ,T3=547.5 0C

66

Figure (4.12)The Effect of Turbine Inlet Temperature (High Pressure Turbine)on

67

the Theoretical Overall Efficiency ,Expansion Ratio=1.12 Figure(4.13) The Effect of Turbine Inlet Temperature (High Pressure Turbine)on

67

the Theoretical Overall Efficiency ,Expansion Ratio=1.08 Figure (4.14) Theoretical Overall Efficiency as a Function of Expansion Ratio with Isentropic Efficiency (ηt,ηc =0.9) andT1=150C

68

Figure (4.15) Theoretical Overall Efficiency as a Function of Expansion Ratio with Isentropic Efficiency (ηηt,ηc =0.9) andT1=300C

68

Figure (4.16) The Effect of Compressor Inlet Temperature (with Variable Expansion Ratio)on the Theoretical Overall Efficiency, T3= 669 0C

69

Figure(4.17)The Effect of Compressor Inlet Temperature (with Variable Expansion Ratio) on the Theoretical Overall Efficiency, T3= 628.5 0C

69

Figure (5.1) Comparison between the Experimental and Present Model for the Effect of Air Velocity on Heat Load at Water Flow Rate 2000 (L/h)

93

Figure (5. 2) Comparison between the Experimental and Present Model for the Effect of Air Velocity on Heat Transfer Coefficient.ha (w/m2.c) at Water Flow Rate 2000 (L/h) Figure (5. 3) Variation Heat Transfer Coefficient ha.(w/m2.c) a long Heat Exchanger Height at Water Flow Rate 2000 (L/h), Air Flow Rate 2000 cfm

94

Figure(5. 4) Variation Heat Transfer Coefficient ha.(w/m2.c) a long Heat Exchanger Height at Water Flow Rate 2000 (L/h),Air Flow Rate2000 cfm

96

Figure (5. 5) Variation Heat Transfer Coefficient. ha ( w/m2. c) a long Heat Exchanger Depth at Water Flow Rate 2000 (L/h), Water Entering Temp. 10 0C, Air Flow Rate 2000 cfm

97

Figure (5. 6) Comparison between the Experimental and Present Model for the Effect of Air Velocity on Overall Heat Transfer Coefficient (w/m2. c) at Water Flow Rate 2000 (L/h)

98

Figure(5. 7) Variation Overall Heat Transfer Coefficient (w/m2.c) a long Heat Exchanger Height with Water Flow Rate 2000 (L/h), Air Flow Rate 2000 cfm

99

95

Figure(5. 8) Variation Overall Heat Transfer Coefficient (w/m2.c) a long Heat Exchanger Height at Water Flow Rate 2000 (L/h), Air Flow Rate 500 cfm

100

Figure(5. 9) Variation Overall Heat Transfer Coefficient ( w/m2. c) a long Heat Exchanger Depth at Water Flow Rate 2000 (L/h), Water Entering Temp. 10 0C, Air Flow Rate 2000 cfm

101

Figure (5. 10) Comparison between the Experimental and Present Model for the Effect of Air Velocity on Air Exit Temperature (0C) At Water Flow Rate 2000 (L/h)

102

Figure (5. 11) Variation Exit Air Temperature (0C) a long Heat Exchanger Height with Water Flow Rate 2000 (L/h), Air Flow Rate 2000 cfm

103

Figure(5. 12)Variation Exit Air Temperature (0C) a long Heat Exchanger Height atWater Flow Rate 2000 (L/h), Air Flow Rate 500 cfm

104

Figure(5. 13) Variation Air Exit Temperature ( 0C) a long Heat Exchanger Depth at Water Flow Rate 2000 (L/h), Water Entering Temp. 10 0C,2000 cfm

105

Figure (5. 14) The Effect of Aspect Ratio (H/L) with Different Core Size (L × D × H) on the Pressure Drop in Air Side Water Flow Rate 2000 (L/h) ,Air Flow Rate 500 cfm

105

Figure (5.15) The Effect of Turbine Inlet Temperature on the Power Output with Variable Compressor Inlet Temperature at Rang Gas Generator Speed (50000-65000 RPM)

106

Figure (5.16) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed ) on the Fuel Mass Flow Rate

106

Figure (5.17) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed ) on the Specific Fuel Consumption

107

Figure (5.18) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the Heat Consumption

107

Figure (5.19) The Effect of Turbine Inlet Temperature (with Variable Compressor Inlet Temperature) on the Heat Rate, at Rang Gas Generator Speed (50000-65000 RPM)

108

Figure (5.20) The Effect of Turbine Inlet Temperature (with Variable Compressor Inlet Temperature) on the Overall Efficiency (%),at Rang Gas Generator Speed (50000-65000 RPM)

108

Figure (5.21) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the on Air Mass Flow Rate

109

Figure (5.22) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the on Air Flow Rate

109

Figure (5.23) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the on Compression Ratio

110

Figure (5.24) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the on Compressor Work

110

Figure (5.25) The Effect of Turbine Inlet Temperature (°C )on the HP.Turbine Work (kW) with variable Compressor Inlet Temperature (°C) At Rang Gas Generator Speed (50000-65000 RPM)

111

Figure (5.26) The Influence Compressor Inlet Temperature (°C )on Gas Turbine Performance at (Gas Generator Speed 45000 RPM)

112

Figure (5.27) The Influence Compressor Inlet Temperature (°C )on Gas Turbine Performance at(Gas Generator Speed 55000 RPM)

112

Figure (5.28) Comparison between the Experimental and Theoretical Predictions of the gas turbine engine (GT-85)

113

List of Tables Title

Pages

Table (2-1) GE-Design Parameters

12

Table (2-2) Turbine Inlet Cooling Options

14

Table (3.1) Heat Exchanger Geometry

20

Table (5.1) Air Temperature Distribution Along H.EX.Depth,Vw=2000

83

0

(l/h),Va=2000 (cfm) , Tw (in)=10 C Table (5.2) Air Temperature Distribution Along H.EX.Depth,Vw=2000 0

(l/h),Va=2000 (cfm), Tw (in)=50

84

C

Table (5.3) Air Temperature Distribution Along H.EX.Depth,Vw=2000

85

0

(l/h),Va=500 (cfm), Tw (in)=10 C Table (5.4) Air Temperature Distribution Along H.EX.Depth,Vw=2000 (l/h),Va=500 (cfm), Tw (in)=50

0

86

C

Table (5-5) Air Temperature Distribution a long First Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37

87

Table (5-6) Air Temperature Distribution a long Second Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37

88

Table (5-7) Air Temperature Distribution a long First Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37

89

Table (5-8) Air Temperature Distribution a long Second Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37

90

Table (5-9) Air Temperature Distribution a long First Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37

91

Table (5-10) Air Temperature Distribution a long Second Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37

91

Table (5-11) Air Temperature Distribution a long First Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37

91

Table (5-12) Air Temperature Distribution a long Second Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37

92

Nomenclature Symbol

Description

Units

A

Total Heat Transfer Surface Area



Af

Fin Area



Acw

Crosse Sectional Area of Water Side

m2

Aca

Crosse Sectional Area of Air Side

m2

Aexp

Exposed Area of the Bare Tube

Cpa

Specific Heat Capacity of Dry Air

J/kg..k

Cpw

Specific Heat Capacity of Water

J/kg..k

Cpg

Specific Heat Capacity of Gas

J/kg..k

dh

Hydraulic Diameter

mm

Dt

Tube Depth

mm

Df

Fin Depth

mm

D

Heat Exchanger Depth

cm

ha

Convection Heat Transfer Coefficient of Air Side

W/m².k

hw

Convection Heat Transfer Coefficient of Water Side

W/m².k

Ht

Tube Height

cm

H

Heat Exchanger Height

cm

H.V

Heating Value

Lf

Fin Length

mm

L

Heat Exchanger Length

cm

m& a

Air Mass Flow Rate

Kg/s

m& f

Fuel Mass Flow Rate

Kg/s

Nu

Nusselt Number =

Nt

Number of Tubes

-

Nr

Number of Rows

-

N

Number of Slices

-

P

Pressure

bar

Q&

Rate of Heat Transfer Loss or Gain

W

kJ/kg

hd k

-

ρud h μ

Re

Reynolds Number =

R

Gas Cosatant

rC

Compression Ratio

-

rT

Expansion Ratio

-

sg

Specific Gravity of Fuel

-

T

Temperature

°C

T

Gauge Measurement Temperature

°C

tf

Fin Thickness

mm

tt

Tube Thickness

mm

Uo

Overall Heat Transfer Ccoefficient



.

J/kg.K

W/m².˚K

V

Volumetric Flow

m3/s

v

Velocity

m/s

W

Power

Watt

XT

Transverse Space

mm

XL

Longitudinal Space

mm

xc

Isentropic Temperature Ratio for Compressor =

xt

Isentropic Temperature Ratio for Turbine =

.

T3 T4 s

T2 s T1

-

Greek Letters

ηo

= 1−

ηf

Fin Efficiency

ρ

Density (kg/m³)

μ

Viscosity (kg/m.s)

Ø

Maximum to Minimum Temperature Ratio through the gas turbine engine

α

t*

ηo

Overall Efficiency of Gas Turbine

ƒ

Fuel / Air Ratio

λt

Perimeter of Tube Side (m3)

λf

Perimeter of Fin (m3)

γ

Specific Heat Ratio

Ŧ

Torque (N.m)

Af A

(1 − η ) Overall Surface f

Efficiency

Ø η c* η

Subscripts a

Air

C

Corrected Value

c

Compressor

f

Fin

g

Gas

i

Input

o

Output

s

Isentropic Process

t

Tube

T

turbine

w

Water

1

Compressor Inlet

2 3

Compressor Discharge 1stTurbine Inlet

4

1st Turbine Outlet

5

2nd Turbine Outlet

Abbreviations NDCW Non-dimensional Compressor Work NDTW Non-dimensional Turbine Work c.c

Combustion Chamber

ISO

International Standard Organization

HVAC

Heating, Ventilating and Air Conditioning

CBT

Compressor Burner Turbine simple cycle

LMTD

Logarithmic Mean Temperature Difference

HP

High Pressure

LP

Low Pressure

CHAPTER ONE

Introduction

Chapter One…………………………………………………………………….Introduction

1.1 General Gas turbines are used in a wide range of services; they power aircraft of all types and drive mechanical equipment such as pumps, compressors, and generators in electric utilities. They also generate power for peak loads and base-load duties. Recently, the interest in gas turbines has grown significantly in combined-cycle plants. These plants use combinations of gas and steam turbines in various configurations of turbines, heat recovery steam generators, and regenerators. Gas turbines have many advantages over steam plants. These are as followed: 1. They are smaller in size, mass, and initial cost per unit output. 2. Their delivery time is relatively short and they can be installed quickly. 3. Their starting is quicker (as low as 10 s) Philip (2002) [1], often by remote control. 4. Their running are smooth and have a capacity factor (percent of time the unit is operating at full power) of 96 to 98 percent. 5. They can be used in a wide variety of liquid and gaseous fuels including gasified coal and synthetic fuels. 6. They can be subjected to fewer environmental restrictions other than prime movers. 1.2 Two Shaft Engine A two shaft gas turbine Figure (1.1) consists of an air compressor, a combustor, a gas generator turbine, and a power turbine. The air compressor generates air at a high pressure, which is fed into to the combustor where the fuel is burned. The combustion products and excess air leave the combustor at high pressure and high temperature. This gas is expanded in the gas generator turbine, which has the sole task of providing power to turn the air compressor. After leaving the gas generator turbine, the gas still has a high pressure and a high temperature. It is now further expanded in the power turbine. The power

Chapter One…………………………………………………………………….Introduction

turbine is connected to the driven equipment. It must be noted at this point that the power turbine (together with the driven equipment) can run at a speed which is independent of the speed of the gas generator portion of the gas turbine. The gas generator is controlled by the amount of fuel which was supplied to the combustor. Its two operating constraints: the firing temperature and the maximum gas generator speed. If the fuel flow is increased, both firing temperature and gas generator speed will also increased until one of the two operating limits is reached.

Exhaust Combustion Chamber

Power Turbine

Turbine

Compressor

Driven Equipment

Gas Generator

Inlet Air

Figure 1.1 Two Shaft Gas Turbine Engine

1.3 Influence of External Factors on Gas Turbine Performance A gas turbine uses ambient air; therefore, its performance is greatly affected by all factors that influence the flow rate of air delivered to the compressor, in terms of weight and its physical conditions. These factors are: 1. Ambient Temperature 2. Ambient Pressure 3. Relative humidity In this regard, the reference conditions for the three above-mentioned factors are (15 °C, 1013 mbar, and 60 %) respectively, ISO (1973)

[2]

.When the

compressor inlet temperature increases, the specific work needed to compress

Chapter One…………………………………………………………………….Introduction

air will also be increased. However, the weight of the air delivered will be decreased (because of a decrease in specific weight). Consequently, the turbine efficiency and useful work (and, therefore, power) diminish as well. If compressor inlet temperature decreases, the reverse process occurs. This temperature depends on the air aspirated by the compressor. The power and efficiency varies from turbine to turbine, according to cycle parameters, compression and expansion output and air delivery rate…etc. And as a result, the variation ratio of gas turbine performance parameters is taken proportional to design point (manufactured levels). Figure (1.2) shows an example of how power, heat consumption, heat rate and the delivery rate of exhaust gases depend on ambient temperature. Design point performance is a central to the engine concept design process. The engine configuration, cycle parameters, component performance levels and sizes are selected to meet the given specification.

Figure 1.2 Influences of External Factors on the Gas Turbine Performance, Frank (2002) [3]

Chapter One…………………………………………………………………….Introduction

1.4 Gas Turbine Inlet Air Cooling The turbine inlet air cooling methods can be divided into two categories: 1.4.1 Evaporative Cooling In this process water is distributed over pads of fibers through which the passing air should be humidified. Spray intercoolers or fogging systems were also used to cool the inlet air. When the power and efficiency can be increased

by decreasing compressor inlet temperature ,the latter can be

reduced artificially by using an evaporative cooler located upstream of the suction filter. Water, fractioned into drops or in the form of a liquid film, cools the air by evaporating in the cooler as it flows in contrary direction.

Figure 1.3 Evaporative Cooler

1.4.2 Cooling with Absorption Chiller The absorption chiller works on the principle of vapor absorption refrigeration cycle. The main advantage of this chiller lies in the fact that the inlet air can be cooled down to a specific temperature for a wide range of ambient air temperatures. Thus, the power output of a gas turbine remains more or less constant, independent of ambient air conditions.

Chapter One…………………………………………………………………….Introduction

The low grade exhaust energy can be used to drive the chiller. The chilled water produced by the absorption system, is passed through the inlet air cooler, which is an indirect type air to water heat exchanger.

Figure 1.4 Air "Chilling Cooling" System, Based on Absorption, GE [4]

1.5 Compact Heat Exchanger Compact heat exchangers have been widely used in various applications in thermal fluid systems including automotive thermal management systems. Radiators for engine cooling systems, evaporators and condensers for HVAC systems, oil coolers, and intercoolers are typical examples of the compact heat exchangers which can be found in ground vehicles. Among the different types of heat exchangers for engine cooling applications, cross flow compact heat exchangers with plain fins are of a special interest. This is because of their higher heat performance capability with the lower flow resistance.

Chapter One…………………………………………………………………….Introduction

1.6 Aims of the present Work The drop in overall performance of gas turbine engines is believed to be due to the increasing of compressor inlet temperature above ISO condition (15°C).The present work is an attempt to improve the gas turbine performance and bringing it near to ISO condition. This is accomplished by applying sensible cooling technique consisting of air cooled heat exchanger and water supplying system. In the present work a theoretical and experimental study has been developed for both, air cooled heat exchanger and gas turbine engine. A predictive numerical model for the air cooled heat exchanger has been developed which is based on step by step technique method so that a design tool for the heat exchanger can be developed. Such work requires developing a versatile experimental facility to examine the air cooling effect on gas turbine performance. This aim can not be achieved without some important requirements which can be summarized as follows: 1- Studying the performance of an existing gas turbine system available in the laboratory. This was achieved by controlling the inlet air temperature to the compressor. 2- An experimental rig was built up to conditioned air at different temperatures controlled by the water supply temperature and then studying this effect on the gas turbine performance. 3- Making theoretical assessments for the heat exchanger performance and the effect of compressor inlet temperature on the gas turbine performance. 4-Developing a computational program for simulate the heat exchanger design and estimating its performance with variable inlet conditions. This is done by feeding with the experimental data to validate this simulation. A computational program to predict the effect of air inlet temperature on the gas turbine performance was built for this object depending on the engine design information only.

CHAPTER TWO

Literature Survey

Chapter Two……………………………………………………………….Literature Survey

2.1 General The gas turbines are generally used for large scale power applications. The basic gas turbine cycle has low thermal efficiency. So it is important to look for improved gas turbine based cycles. The inlet air cooling helps to increase the performance of gas turbines. The demand of energy in the developing regions of the world, particularly in Asia, has witnessed pronounced increase in the recent past. According to a report of International Energy Outlook (2004)[5], the world net power consumption is expected to be doubles nearly over the next two decades. Much of the growth in new electricity demand is expected to come from countries of the developing world. Therefore, it is important to find improved technologies for power generations that have a high efficiency and specific power output, low emissions of pollutants, low investment, and low operating and maintenance cost for a sustainable use of available fuels. 2.2 Pervious Work The previous work can be classified into two categories: The first one is concerned with modeling and simulation of air cooled heat exchanger. The other related with the performance improvement of the gas turbine engine. The improvement of two-shaft gas turbine performance by making modifications or addition of some parts to the main components. This will lead to decrease in the compressor inlet temperature and enhance the power output and thermal efficiency of gas turbine. 2.2.1 Modeling and Simulation of Air Cooled Heat Exchanger Ganapathy (1979) [6] concluded that for air- cooled condensers, the ambient air is the most important variable in the design. Since ambient temperature in a location varies throughout the year. Using higher value, would result in over sizing the unit. Where as A lower value would give poor performance. Current practice is to use a design temperature that exceeds (2 to 5%) of the annual period.

Chapter Two……………………………………………………………….Literature Survey

Hedderich and Kellehere (1982)[7] developed a computer code for the analysis of air cooled heat exchangers and was coupled with a numerical optimization program to produce an automated air cooled heat exchanger design leading to optimization procedure. A general iteration free approximation method was used for the analysis which calculates the mean overall heat transfer coefficient and the overall pressure drop for many arrangements. The analysis takes into account the variation of heat transfer coefficients and the pressure drop with temperature and length of flow path .The capability is demonstrated by the design of an air cooled finned tube heat exchanger and is shown to be useful tool for the heat exchanger design. Zhang (1994) [8] proposed three dimension numerical model predicting the performance for large power plant condensers. He compared his predicted results with the experimental data. The prediction was achieved by solving the governing mass, momentum and air concentration using semi implicit consistent control-volume for simulation model with different conditions in work of condenser. Matthew and Joseph (2002) [9] developed a conceptual designs for wet and dry cooling systems applied to a new ,gas –fired, combined cycle 500-MW plant at four sites chosen to represent conditions in California. The requirements for cooling dry systems are four to six times those for wet systems. Dry systems, which are limited by the ambient dry bulb temperature, cannot be achieved as low a turbine back pressure as wet systems, which are limited by the ambient wet bulb temperature. Dohoy and Dennis (2006)[10] divided the heat exchanger core into small control volumes along the tube. Finite Difference Method (FDM) with staggered grid system was utilized in study. FDM can take into account the significant air temperature increase as well as the local variations of the properties and the heat transfer coefficient. The maximum difference between

Chapter Two……………………………………………………………….Literature Survey

the experimental data and calculated results was noticed to be 5% for the given range of the simulated conditions. Tarrad, et al. (2008)

[11]

investigated the performance prediction of the

cross flow air-cooled heat exchanger. They developed a new simplified correlation for the air side heat transfer coefficient which depends on the dimensional analysis with Buckingham-pi theorem. The discrepancy between the predicted and their own experimental values of the overall heat transfer coefficient and heat duty were within 2% and 4% respectively for both of the tested tube banks. 2.2.2 Performance Improvement of the Gas Turbine Engine There are two basic methods available for inlet air-cooling evaporative cooling and chilling cooling. The most widely accepted system is evaporative air-cooling. Evaporative coolers make use of the evaporation of water, and are the most cost-effective way to improve machine capacity during warm weather. Mostly percent of design concept is used to examine gas turbine performance. Of the two cooling methods of inlet air, namely, evaporative cooling and the absorption cooling, the absorption cooling technique demonstrated a higher gain in power output and efficiency than evaporative cooling for a simple cycle gas turbine. De Lucia, et al.(1995)

[12]

reported that evaporative inlet-cooling is

economical and simple, but only suitable for dry hot climates. they concluded that evaporative inlet cooling could enhance power out put by (2–4)% depending on the weather. Saleh (1996)[13] presented that water can be injected in a simple two shaft gas turbine (GT-85) to improve the performance. The studied cases were water injection prior to the combustion chamber and water injection in intake of the compressor. Maximum increasing in performance data was obtained when water injection prior to the combustion chamber where, the increase of

Chapter Two……………………………………………………………….Literature Survey

the power output up to (40.4٪) and thermal efficiency up to (34.8٪) .In addition, specific fuel consumption was reduced by (32.12٪ ). Ait-Ali (2001)

[14]

presented the concept of inlet air refrigeration to boost

the power output from the gas turbine. Chillers can increase the gas turbine power output by 15-20%.The absorption chiller works on the principle of vapor absorption refrigeration cycle. The main advantage of this chiller lies in the fact that the inlet air can be cooled down to a specific temperature for a wide range of ambient air temperatures. Therefore, the power output of a gas turbine remains more or less constant, independent of ambient air conditions. A typical absorption chiller with a capacity of 3000 refrigeration tons and a COP of 0.70. This absorption system uses the waste heat to produce required steam by the chiller. Nuovo Pignone (2002)[15] publishing curves of the compressor inlet temperature effect on gas turbine performance as shown in fig.(2.1).It is obvious that power output and heat rate are improved as compressor inlet temperature was decreased. Lowering the compressor inlet temperature can be accomplished by installing an evaporative cooler or inlet chiller in the inlet ducting downstream of the inlet filters. General electric estimates theses performance curves with respect to the design point and maximum speed as shown below in table (2-1). Table (2-1) GE-Design Parameters

Chapter Two……………………………………………………………….Literature Survey

Figure 2.1 Effect of Compressor Inlet Temperature on Gas Turbine Performance [15]

Donald and Icksoo (2003)[16] presented the various types of turbine inlet cooling applicable to small to mid-size turbines .These have been described along with their comparative benefits. The greatest benefit was shown to be obtainable from an exhaust heat-powered ammonia absorption cycle. An ammonia absorption cycle was especially designed for this application. A 300-refrigeration ton aqua ammonia refrigeration unit is required to cool the inlet of a (5 MW) gas turbine from (35°C to 5°C). This cooling will increase

Chapter Two……………………………………………………………….Literature Survey

the power output by 1 MW. The added power was at a marginal efficiency of 39%, compared to 29% for the base turbine power. The cooling option is listed in table (2-2). Table (2-2) Turbine Inlet Cooling Options 4.7 MWe Simple Brayton Cycle - 30% efficiency at ISO

Hameed (2004)[17] concluded that water injection in

the air intake is

strongly effecting the performance parameters of the two-shaft gas turbine cycle (GT-85) . the power output has been increased up to (23.15%) for simple cycle and the thermal efficiency is higher than that of normal cycle by (29%). Benjalool (2006)

[18]

concluded that, in September, the range of ambient

temperature in the Nafoora oil field varies typically between (29-36 °C). The temperature variation leads to change the maximum engine power output from (5.1 MW) to (4.85 MW). Tony Giampaolo (2006)

[19]

concluded that at a constant gas generator

speed, and ambient temperature decreases, turbine inlet temperature will be decreased slightly, and power output will be increased significantly, fig. (2.2).This increase in gas horsepower results from the increase in compressor pressure ratio and aerodynamic loading. Therefore, the control must protect

Chapter Two……………………………………………………………….Literature Survey

the gas turbine on cold days from overloading the compressor airfoils and over-pressurizing the compressor cases. Sensing ambient inlet temperature helps insure that engine internal pressures are not exceeded, and sensing turbine inlet temperature insures that the maximum allowable turbine temperatures are not exceeded. Sensing gas generator speed enables the control to accelerate through any critical speed points (gas turbines are typically flexible shaft machines and, therefore, have a low critical speed).

Figure 2.2 Temperature-Power-Speed Interrelationships

Kuamit (2006) [20] concluded that the effect of compressor inlet temperature has an important role on the power output as shown in fig.(2.3). It may be seen that the power out put is influenced by compressor inlet temperature due to the change of air density and compressor work .Since a lower compressor inlet temperature leads to a higher air density and a lower compressor work that in turn gives a higher gas turbine output.

Chapter Two……………………………………………………………….Literature Survey 160

Output Power, (M W )

150

140

130

120

110

100 260

270

280

290

300

310

320

330

340

Inlet Temp, (K)

Figure 2.3 Net Output Power Versus Inlet Temperature for Gas Cycle

2.3 Summary and Motivation From the review of previous works it is obvious that industrial gas turbines performance is different from turbine to turbine. It depends on the type and what full performance at ISO conditions .Also, what equipment that used to improve the performance. Clearly, the performance is not only affected by compressor inlet temperature, but also by other parameter such as relative humidity, inlet pressure, maximum temperature in the cycle, and speed of the shaft. However, operating at constant speed has constant volumetric flow rate. Since the specific volume of air is directly proportional to temperature, cooler air has a higher mass flow rate. It generates more power in the turbine. Cooler air also requires less energy to be compressed to the same pressure as warmer air. Thus, gas turbines generate higher power output when incoming air is cooler. A gas turbine inlet air cooling system is a good option for applications where electricity prices increase during the warm months. The power output increases by decreasing the compressor inlet temperature of the incoming air. The aim of this work is to study the approach used in the enhancement of gas turbine performance by making the compressor inlet temperature supplied to the gas turbine closed to ISO conditions. This was achieved by building a

Chapter Two……………………………………………………………….Literature Survey

cooling system for entering air to the compressor consisting of air cooled heat exchanger and water supply system. In this regard, a theoretical and experimental study for both air cooled heat exchanger and gas turbine engine will be conducted. A supporting computer program to simulate a new technique of heat exchanger design has been developed. In this model, heat exchanger was described in two directions, the height and depth to form horizontal slices which will described later. A computer program prepared for this purpose has the ability to analyze heat exchanger performance for any slice to be located in two heat exchanger dimension. The heat exchanger type used was a type of air cooled heat exchanger (finned-tube surfaces, flat tubes, continuous fins). It is part of cooling system accommodating different components like valves, pipe fittings, supply pump, and two reservoirs (hot reservoir and cold reservoir).

CHAPTER THREE

EXPERIMENTAL WORK

Chapter Three…………………………………………………………Experimental Work

3.1 Test Rig for Heat Exchanger The test rig layout was built as shown in Figure (3.1). It consists of inlet air cooling system having air cooled heat exchanger, tank, and water circulation system, pump and control panel. The water circulation system is a modification of an existing laboratory system, which had already built by Tarrad and Mohammed (2006) [21]. 3.1.1 Heat Exchanger The test section is made of a compact heat exchanger which is of 19 80Chevette radiator type using, Cross-flow exchanger with one flow mixed. Edges of flat vertical tubes, having dimensions of (55 cm) length,(3.5 cm) depth and (37 cm) height, heat exchanger geometry is illustrated in table (3.1). The Compact heat exchanger configuration is shown in Fig. (3.2). The water flows in the tubes in cross direction to the air flowing normal to the tubes. Thermometers are connected to the heat exchanger. The gauges are fixed in the specially prepared pockets mounted on the required locations. 3.1.2 Water Supplying System The cold water is supplied by a constant tank head of (200 liters) capacity. The water is pumped by a single stage centrifugal pump from the tank through the test section, and then it returns back to the tank. The water in the hot tank is heated by four electrical heaters, which have a total heating electrical power of (12 kW). The four heaters are fixed in the same level at 150 mm from the bottom of the hot tank separated at (90◦) apart. The hot tank is opened in cold tank by gating valve which be found in the pipe that linked the two tanks, as shown in Fig. (3.1). This arrangement of cold and hot tanks enable as to control manually the temperature in the tank .It is will help in providing a good control of the water temperature during tests. The mixing of water has been accomplished in the cold tank to obtain the required

Chapter Three…………………………………………………………Experimental Work

temperature. A piece of ice was immersed into cold tank so that the mixture becomes homogenous. Table (3.1) Heat Exchanger Geometry

Parameter Core

Fin

Tube

Dimension Length (L)

550 mm

Depth (D)

35 mm

Height (H)

370 mm

No.of Tubes

110

No.of tube/row

55

No.of fins/ tube

256

Normal Distance (XT)

9.92 mm

Longitudinal Distance (XL)

20.46 mm

Pitch (Pf)

1.46 mm

Length (Lf)

7.52 mm

Depth (Df)

15.88 mm

Thick

0.24 mm

(Tf)

Height (Ht)

2.4 mm

Depth

(Dt)

15.88 mm

Thick

(Tt)

0.28 mm

A special transparent glass tube level is fixed on the outer shell of the tanks in order to monitor the water level. Temperature gauge (thermometer) is connected to the shell of tank to monitor the water temperature in the tank during tests. Both hot and cold water flow is controlled by gate valves and the flow is measured by using a vertical variable area rotameters. Piping system was made of carbon steel metal, insulated by glass wool to minimize the heat loss.

Chapter Three…………………………………………………………Experimental Work

3.1.3 Measuring Instrumentation The parameters which are to be measured during the test as follow: (1) The inlet and outlet temperatures of the heat exchanger. (2) The air flow rate across the heat exchanger. (3) The water flow rate of the tube side. Thermometers: Thermocouples having a range of (0 ◦C – 120 ◦C), are used to measure the temperatures at inlet and outlet of heat exchanger. The accuracy rang is about (0.02). Pressure Gauges: Two pointer pressure gauges installed on both sides of heat exchanger to measure the pressure of the water and have a range of (0 – 2.5 bar), which are connected to the heat exchanger the accuracy rang about (0.03). Rotameters : Two vertical variable area rotameters are used to measure the flow rates of the water. A rotameter of (200 – 3000) L/h range is used the accuracy rang about (0.02). Pitot Tube: Pitot tube is located at end of air flow duct is used to measure air flow rates by calculating pressure difference for two points.

3.1.4 The Electrical Board The electrical board contains the main circuit breaker and other secondary switches, supplying power to the whole system components of the test rig. The board consists of an electrical contactor of (4×16 Amp), which is connected to a thermostat and four heaters. This contactor is controlled and receives the electrical signal from the thermostat switching on and off the electricity to the four heaters. There are four separate switches each of which controls one heater. In addition for safety, there is a switch controling the

Chapter Three…………………………………………………………Experimental Work

operation of the contactor (i.e., controls the operation of the four heaters and the thermometer), Tarrad and Mohmmed [21]. 3.1.5 Air Circulation System The air was supplied to the test heat exchanger through a fan. A forced draught arrangement was selected for the test object by variable fan speed. Three volumetric flow rate was prepared of capacity of (2000) cfm ,(1000 cfm) and (500 cfm) .The fan was close enough to the test section avoiding leakage of air to the surrounding. The air volumetric flow rate is measured by using (Pitot-Tube) .The pressure drop was calculated for two points through air flow duct.

Chapter Three…………………………………………………………Experimental Work

Heat Exchanger Duct

T

Fan

P

By Pass

Rotameter

T

Rotameter

P

Valve

Hot Tank

Cold Tank

Heaters

Ice

Valve

Pump

Valve

Figure (3.1a) Schematic Diagram of the Built Rig (Heat Exchanger)

Chapter Three…………………………………………………………Experimental Work

Fan

Cold Tank

Figure (3.1b) Configuration of the Built Rig (Heat Exchanger)

Pitot-Tube Rotameter

Heat Exchanger Hot Tank

Gas Turbine

Figure (3.1c) Configuration of the Test Rig (Inlet Cooling System) Preparing to Gas Turbine

Chapter Three…………………………………………………………Experimental Work

XL Dt XT Ht

Figuer (3.2.a) Top View of Heat Exchanger Geometry

Pf Lf H

D

L Figure (3.2.b) Front View of Heat Exchanger Geometry

Chapter Three…………………………………………………………Experimental Work

3.2 Experimental Setup A schematic diagram of the experimental setup is shown in Figure (3.3). The test loop was designed to allow easy control parameters such as upstream air flow rate to heat exchanger, and compressor inlet temperature to gas turbine. As mentioned in the previous section, the loop consists of the heat exchanger set prior to gas turbine apposite inlet box. To achieve this it is required

building

up

the

main

experimental

rig

with

measuring

instrumentations obtaining the suitable circumstances regarding air dry bulb temperatures, and flow rate. The air leaving the heat exchanger is fed directly to the compressor intake of the gas turbine engine.

3.3 Gas Turbine The chosen test engine for this present work is a two-shaft gas turbine type (GT-85) of 5 kW two shaft machine for industrial drive (Pumps, Compressors).The design strategy of a gas turbine unit is divided into two fields. These are; design information and performance information. Design information such as the compressor whether it is a centrifugal or axial flow, pressure ratio, the number of fuel nozzle in combustion chamber, a single shaft or two shaft. The performance information such as output power, heat rate, exhaust flow, exhaust temperatures. GT-85 performance parameter measurements have been obtained over variable gas generator speed. The speed variation is accomplished by manual fuel flow control valve. The measuring parameters obtained for a particular gas generator speed (N1) are, Vf,P2,T1,T2,T3,P3,T4,P4,T5,N2.

Pump

Valve

Rotameters

Cold Tank

Cold Tank

Valve

Hot Tank Hot Tank

Compressor

Duct

Heat Exchanger

Figure (3.3) Schematic Diagram of the Test Rig (Overall)

Ambient Air

Combustor

Chapter Three…………………………………………………………Experimental Work

Dynomometer

LP.Turbine

HP.Turbine

Chapter Three…………………………………………………………Experimental Work

3.3.1 Design Information Compressor: • Centrifugal flow • Pressure ratio , 2:1 Combustion: • Annular combustion, single chamber • Liquid fuel capability Turbine: • High Pressure turbine 1 stage , 1.69:1 expansion ratio • Low Pressure turbine 1 stage 1.18:1 expansion ratio Package: • Gas Generator, Power Turbine and auxiliary system mounted on a single base plate • Control system (instrumentation, sensors ,electronic panel, mechanical regulation) to monitor temperature and pressure 3.3. 2 Performance design • Power turbine output: 5 KW • Maximum cycle temperature: 750 (°C) • Maximum Gas generator speed : 90000 RPM • Maximum Power turbine speed: 35000 RPM The performance power, fuel consumption, temperatures, shaft speeds... etc. of a gas turbine engine is crucially dependent upon its inlet and exit conditions. The environmental envelop is impartment item for any gas turbine plant performance. 3.4 Operating Principles A gas turbine works in the following way: • It aspirates air from the surrounding environment • It compresses it to a higher pressure • It increases the energy level of compressed air by the addition of fuel gas

Chapter Three…………………………………………………………Experimental Work

which undergoes combustion in a combustion chamber • It directs high pressure and high temperature air to a turbine section, which converts thermal energy into mechanical energy allowing shaft to revolve. This serves on one hand, to supply useful energy to the driven machine, coupled to the machine by means of a coupling and, on the other hand, to supply energy necessary for air compression, which takes place in a compressor connected directly with the turbine section itself. • The remaining energy is supplying through power turbine, finally • Gas turbine sensors are designed to withstand the vibration and high temperatures found in these engines ,Watlow (2001) [22]. 3.5 Main Component Parts of Gas Turbine 3.5.1 Gas Generator: The gas generator consists of a centrifugal compressor; a combustion chamber and a radial flow turbine as shown in Fig. (3.3).Air enters the centrifugal compressor through an intake silencer and a bell mouth air flowmeter. The outlet of the compressor is directed towards a vertically oriented combustion chamber which in turn is connected to the compressor turbine. The compressor and turbine run at a maximum speed of 90000 RPM. Such a compressor speed can achieve a compression ratio of 2:1 for this particular compressor design. The speed of the compressor and turbine is measured by means of an opto-electronic technique. 3.5.1.1 Centrifugal Compressor : The centrifugal compressor draws in air at the center or the eye of the impeller and accelerates it around and outward. It consists of an impeller, a diffuser and compressor manifold. The diffuser is bolted to the manifold, and often the entire assembly is referred to as the diffuser. The impeller may be either single entry or dual entry. The principal differences between a single entry and the dual are the size of the impeller and the ducting arrangement. The single entry impeller permits ducting directly to the inducer vanes, as

Chapter Three…………………………………………………………Experimental Work

opposed to the more complicated ducting needed to reach the rear side of the dual-entry type .Although it is slightly more efficient in receiving air, the single-entry impellers must be of a greater diameter to provide sufficient air. The compressor draws in air at the hub of the impeller and accelerates it radially outward by centrifugal force through the impeller. It leaves the impeller at high speed and low pressure flowing through the diffuser. The diffuser converts the high speed, low–pressure air to low-speed, highpressure air. The compressor manifold diverts the low-speed, high-pressure air from the diffuser into the combustion chamber. In this design, the manifold has one outlet port for each combustion chamber. The outlet ports are bolted to an outlet elbow on the manifold. The outlet ports ensure that the same amount of air is delivered to each combustion chamber. The outlet elbows change the airflow from radial to axial flow. The diffusion process is completed after the turn .Each elbow contains from two to four turning vanes that perform the turning process and reduce air pressure losses by providing a smooth turning surface. 3.5.1.2 Combustion Chamber: In a gas turbine, the hydrocarbon fuel is burnt in the combustion chamber as a continuous process at essentially constant pressure. The combustion system of the (GT- 85) consists mainly of two components: the swirl unit incorporating a spray nozzle and the combustion chamber. The swirl unit produces a swirling air flow which mixes with the finely atomized fuel sprayed from the nozzle. The air/fuel mixture issuing from the swirl generator is fed into the combustion chamber as an expanding swirling flow. Schematic drawing Fig.(3.4), shows the flow issuing from swirl generator into the combustion chamber and the subsequent flow patterns. As it will be seen from the figure, the flow initially on entering the combustion chamber should form a toroidar vortex. This vortex being formed as a result of the expanding swirl flow and the air entry through the primary holes. Combustion takes place here

Chapter Three…………………………………………………………Experimental Work

under an approximately stoichiometric air/fuel ratio, i.e., A/F =15:1. The second row of holes, defined as the secondary zone admits more air to cool the combustion products and gives an overall A/F at this point of about 35:1 The dilution zone contains larger holes than any of the preceding sections and admits a larger quantity of air producing an exit gas flow with an overall a/f ratio of approximately 70:1 .For all sections in the combustion chamber, series of small holes is provided around the circumference of the combustion chamber. These holes provide a cooling flow for the walls and maintain the wall temperature within the safe operating temperature. 3.5.1.3 Gas Generator Turbine: A radial flow turbine consists essentially of a stationary casing containing a rotating impeller which is rotated as a result of the high velocity flow leaving the stationary nozzles. The function of the casing is to accelerate the flow smoothly producing uniform high velocity to the impeller tip. The turbine impeller is of radial vane type. Both the compressor and generator turbine run at the maximum speed of 90000 RPM. 3.5.2 Power Turbine: The hot gas issuing from the exhaust of the gas generator turbine is passed through a flexible circular duct to the inlet of the power turbine. The power turbine is larger than the gas generator turbine and is of inward radial flow design. The turbine operates at a maximum speed of 35000 RPM .The maximum power output is (5 kW) with an inlet temperature of 700C and a speed of 35000 RPM .The power turbine is directly coupled to the eddycurrent brake by means of a special high speed coupling .The speed of the power turbine is measured by electromagnetic sensor and a digital frequency meter.

Chapter Three…………………………………………………………Experimental Work

Fig.(3.4) Schematic Arrangement for Gas Turbine (GT-85)

3.5.3 Fuel System: The gas turbine utilizes kerosene as the operating fuel system. During the starting sequence, methylated alcohol (starting fuel) is ignited by the ignition system to provide a large flame in front of for the initiation of the kerosene combustion process. The fuel is pumped by an electrically driven gear pump at a pressure of 6 bars. The fuel flow is regulated by a manuall valve located on the control panel. Fuel flow rate is measured by a ‘Rotometer’. For safety operating process, a solenoid valve is incorporated in the fuel feed piping system to the combustion chamber. This will automatically stop the supply of fuel in the event of malfunction. 3.5.4 Oil Lubricating System: As mentioned before, the gas generator turbine and the power turbine run at very high speed. Therefore, a lubricating system, Fig.(3.4) is incorporated to ensure a safe running to turbines . Oil pressure of minimum 3 bars should

Chapter Three…………………………………………………………Experimental Work

be provided to lubricate the journal bearings by the gear type oil pump. Oil temperature and pressure are automatically monitored and the gas turbine will be shutdown in case of a non-safe running. 3.5.5 Starting System: In order to start the gas turbine, it is required to drive the compressor to a certain speed to achieve the required air flow rate to the combustion chamber, so that the starting fuel (methylated alcohol) can be ignited. This is done by blowing air through the compressor inlet duct using three electrical fans. The flow rate of air supplied by these fans is regulated by a gate valve. High energy spark ignition system is used to ignite the starting fuel. The sustained flame generated by starting fuel burning will raise the flow gases temperature to the required temperature to ignite the kerosene fuel. Consequently, self sustained gas generator running is obtained and the air blowing fans will be switched off Gilbert and Gordon LTD (1978) [23]. 3.6 Dimensionless and Parameter Groups The importance of dimensionless, referred and scaling parameter groups to all aspects of gas turbine performance cannot be over emphasized. Understanding and remembering the form of the parameter group relationships allows judgments concerning the performance effects of changing ambient conditions, scaling an engine, a change of working fluid, see Appendix C (Engine Parameter Groups) The parameter group for mass flow is then a function of , Philip (2004) [24]: 1. Ambient temperature 2. Ambient pressure 3. Engine rotational speed 4. Engine diameter (scale factor) 5. Gas constant of working fluid 7. Gamma for working fluid 8. Viscosity of working fluid

Chapter Three…………………………………………………………Experimental Work

in order to rationalize the performance data obtained from the gas turbine engine operating with a wide range of inlet conditions, it is necessary to reduce the data to known standard intake conditions. The correction formula used to achieve this rationalization can be derived by a non-dimensional analysis of the components in the gas turbine cycle. Non-dimensional analysis leads to various dimensionless parameters which are based on the dimension's mass (M) , length (L) , and time (T) . Based on these elements, one can obtain various independent parameters. These parameters will lead to form various dimensionless groups. By using non-dimensional groups as applied to the following basic equation for compressor non-dimensional analysis, Cohen,et al. (1996) [25]: f ( Nc, ma , P1 , P2 , R × T1, R × T2 ) =constant ……………….…….…..(3-1)

3.6.1 Corrected Compressor Data From the equations are illustrated in Appendix (C), and that referring to equation (3-1) for compressor non-dimensional analysis, and since :-

Ncc Ts

=

Nc T1

.

……………………………………………………………...…(3-2) .

M a Ta Pa

=

M ac Ts Ps

……………………………………………………..…(3-3)

T2 T2c = …………………………………………………………………(3-4) Ta Ts

where: Ts= 15 (°C) =288.16 (k) Ps= 1.0133 bar For other values of Ta different from the standard value a correction should be made as follows:

Chapter Three…………………………………………………………Experimental Work

N cc = N c

288.16 ………………………………………………….….….(3-5) Ta .

.

M ac =

1.0133 288.16

×

M a Ta

……………………………………………..…..(3-6)

Pa

⎛ 288.16 ⎞ ⎟⎟T2 ………………………………………………………….(3-7) T2c = ⎜⎜ ⎝ Ta ⎠

3.6.2 Corrected Basic Data Corrected basic data for the effects of variation in both, ambient pressure and temperature should be done. In the case of the actual test data, the pressure measurements taken are all static values, whilst the temperature measurements will be essentially total values. For the formula below, no differentiation between total and static is made because of the small differences will be small between the total and the static value of pressure or temperature. The following formulas are the corrected basic data, which are suggested by the manufacturer. −

1.

Ta = Ta ( o C ) + 273.16

2.

Pa =

3.

T1C = (T1 + 273.16)(

4.

P1 = Pa − 97.9 Χ10 −6 (ΔP)

5.

P2c = P2 + Pa

Barometer Re ading (mbar ) 1000 −

288.16 ) Ta





T2 c = (T2 + 273.16)(

7.

P3C = P3 + Pa

288.16 ) Ta





T3c = [T3 + 273.16](

9.

P4 c = P4 + Pa



bar K bar bar

6.

8.

K

K bar

288.16 ) Ta

K bar

Chapter Three…………………………………………………………Experimental Work −

288.16 ) Ta

K



288.16 ) Ta

K

10

T4 c = (T4 + 273.16)(

11.

T5c = (T5 + 273.16)(

12.

P5c = P4 c − 0.001333(ΔP4 5 )

bar

3.6.3 Corrected Derived Data The measured variables during the tests on the (GT-85) are presented in tables shown in Appendix (A).When the temperature and pressure at each point around the cycle, the following formulas may be used for the performance calculation in the experiments. The following represents the data reduction method applied in the present work. 1. Compression Ratio ( rc ) rc =

P2c …………………………………………………………………..(3-8) P1c

2. Compressor Isentropic Efficiency (η c ) ηc =

T1c ( γ −1) γ (rc − 1) × 100 0 0 ……………………………….………….(3-9) T2 c − T1c

3. Air Mass Flow Rate (kg/s) m& a

T1 P1

= 0.3005

m& a c = m& a(

ΔP P1

Ta ) 288.16

............................................................................(3-10)

where : ΔP in mm wg T1 in K

Chapter Three…………………………………………………………Experimental Work

4. Fuel Flow Rate (kg/s) m& fc = ρ f

v f × 10 −3

288.16 Ta

3600

…………………………………..……….…(3-11)

where:

ρ f = ρ w ∗ sg

The specific gravity of the fuel (sg) is 0.774 5. Air/Fuel Ratio (A/F) A/F = m& ac m& fc ……………………………………………………………(3-12) 6. Combustion Chamber Temperature Rise (K) −

ΔTcct

43740 − 10 T2 c = +2 ………………………………………..…….(3-13) 1.0078 A F + 6.6

Where:



T2C = T2 c − 273.16

7. Combustion Chamber Pressure Loss (٪) ΔPcc = (

P2 c − P3c ) × 100 ……………………………………………….…..(3-14) P2 c

8. Gas Generator Turbine Expansion Ratio rt1 =

P3c …………………………………………………………………(3-15) P4c

9. Theoretical Power Input to Compressor ( w& ) W& C = m& aCpa(T2 c − T1c ) ……………………………………………………(3-16)

Where: cpa = 1000 J kgK

Chapter Three…………………………………………………………Experimental Work

10. Theoretical Power Output Of Gas Generator Turbine ( w& ) W& t 1 = (m& ac + m& fc )Cp g (T3c − T4c ) ……………………………………………(3-17)

where: cp g = 1150 J kgK 11. Compressor Mechanical Efficiency (η mc ) η mc = (

2 w& c ) × 100 ………………………………………………………(3-18) w& c + wt′1

12. Gas Generator Turbine Mechanical Efficiency (η mt1 ) η mt1 = (

W& C + W& t1 ) × 100 …………………………………………….….……(3-19) 2W& t1

13. Compressor Overall Efficiency (η oc ) η oc = η c ×

η mc 100

……………………………………………………..………(3-20)

14. Power Turbine Expansion Ratio P4 c …………………………………………………………………(3-21) P5c

r12 =

15. Power Turbine Isentropic Efficiency (η12 )

η12 =

T4C − T5C × 100 ………………………………………………(3-22) 1 (γ −1) γ ] T4C [1 − ( ) r12

where: γ = 1.33 16. Power Turbine Power Output ( w& ) 2πN 2τ W& t 2c = 60

Ta ……………………………………………………(3-23) 288.16

Where: N 2 power turbine speed (RPM) and τ is the torque in (N.m)

Chapter Three…………………………………………………………Experimental Work

17. Power Turbine Theoretical Power Output ( w& ) W& t 2 = (m& ac + m& fc )Cp g (T4c − T5c ) ………………………………….………….(3-24)

18. Power Turbine Mechanical Efficiency (η mt 2 )

η mt 2 = (

W& t 2 c ) × 100 …………………………………………………………..(3-25) W& t 2

19. Power Turbine Overall Efficiency (η ot 2 )

η ot 2 = η t 2 ×

η mt 2 100

…………………………………………………………….(3-26)

20. Overall Thermal Efficiency (η th ) η th = (

W& t 2 c ) × 100 …………………………………………………..…(3-27) H .V × m& fc

Where: H.V= L.C.V × 4.1868 kJ kg for kerosene used in the tests (GT-85), and L.C.V= 10300 Kcal kg the lower calorific value of the fuel used, Gilbert and Gordon LTD [23]. 21. Specific Fuel Consumption ( kg Kw.s ) sfc = (

m& fc ) × 1000 ……………………………………………………..…..(3-28) W& t 2c

22. Heat Consumption (kW) Heat Consumption = fuel mass flow * heating value = mf (kg) × H.V (kW/kg)……………………………...(3-29) 23. Heat Rate Heat Rate = specific fuel consumption × heating value = sfc (kg/kW) × H.V (kW/kg)………………………………….(3-30)

Chapter Three…………………………………………………………Experimental Work

3.7 Experimental Work Procedure: Since the main object of present work is the enhancement of gas turbine performance by improved intake air temperature therefore, the test procedure has been separated to classify for heat exchanger once and gas turbine again another time. This was because of the following reasons : 1. Parameters effected on the design and performance calculations for heat exchanger are needed high air flow rate values, which are negatively effect on gas turbine performance. 2. The high inlet velocity into compressor may be increasing pressure drop for what of which negatively effect on gas turbine performance. 3. Difficulty to make sure of compressor inlet temperature at desired value for present work requirements. 3.8 Test Procedure: Two types of experiments were performed: 1. Measurements for the heat exchanger: On commencing the tests, all valves around the constructed rig are closed without water circulation through the heat exchanger. For the chilled water circulation tests, a piece of ice was added to the water tank continuously to keep a constant feeding temperature to the test heat exchanger section. The hot water tests were conducted by switching on the immersion heaters of the hot tank and controlling the temperature by setting the thermostat at the required water temperature. A check should be made to the air supplying fan prior to the experiments. After completing checking above steps, the test process begins by switching on the circuit breaker that supplies power to the system. The water pump will start and open the gate valve for the cold tank that controls the flow rate of water which circulates on the tube side of the heat exchanger. At particular time, when the temperature inside the cold tank is reaching (10°C), the operation conditions for heat exchanger were recorded when the air flow

Chapter Three…………………………………………………………Experimental Work

rate was fixed at (2000 cfm).Since the air flow passes through the duct that connected the heat exchanger to gas turbine, the air flow rate (cfm) can be measured by using (pitot-tube) fixed at compressor inlet box. After a few minutes, the gate valve of the hot tank is opened, so that the hot water is mixed with cold water to increase the temperature gradually until reaching particular values such as (20°C) , (30°C) , (40°C) and (50°C) .The process was repeated for other air flow rate such as (1000 cfm) and (500 cfm) for different circulated water flow rates.The following operating conditions were measured during the tests for each air flow rate: • The inlet and exit temperature of water side across the heat exchanger. • The circulated water flow rate. • The air temperature on both sides of the heat exchanger. The experimental data collected for the heat exchanger are listed in Appendix (A). 2. Measurements for the gas turbine: The Kerosene fuel is supplying in fuel tank about (20 Lit.), for lubricating shafts oil (SAE-10) is supplying in oil tank. The gas turbine engine is connected to a water supplying source to cooling purpose. Starting operation begins by firing of about (0.3) liter of (methylated alcohol) to ensuring the appropriate flame for kerosene combustion. The gas turbine was allowed to

operate until it

reaches the required

conditions at a particular speed; at this time the engine takes the appropriate air flow rate .The compressor inlet temperature is fixed at indicated value such as (15°C) ,(20°C), (25°C) and (30°C) by utilizing from: • Variation of ambient temperature along the day and season. • Variation of air exit temperature from heat exchanger. The time was too long to take the measurement variations of the ambient temperature (during the season) due to the fact that the temperature gradient across heat exchanger (Single-Pass Elliptical Tube, two Row Radiator) is

Chapter Three…………………………………………………………Experimental Work

small. It was decided to run the rig for (5 to 10) times at each specified ambient temperature with different times along the year. When the compressor inlet temperature is fixed at (15°C),the tests were conducted at different gas generator speed starting from (45000 RPM) . This speed was increased to (55000 RPM) by burning more fuel. The data were collected for each specified speed and compressor inlet temperature. The same procedure was repeated for other air intake temperature such as (20°C),( 25°C)and (30°C). At all tests, the ambient temperature was changing from (18°C) to (30°C).To show the relationship between the gas generator speed ,compressor inlet temperature and turbine inlet temperature, the turbine inlet temperature was fixed at a particular value such as (586 °C) and compressor inlet temperature was taken at different time such as (15 °C) and (22°C).The data collected during the tests are shown in Appendix (A).

CHAPTER FOURE

Chapter four…………………………………………………………………………Theory

4.1 General This chapter deals with the mathematical modeling of the core design of heat exchanger (Length, Depth, Height) corresponding to the intake air box of gas turbine and mathematical modeling of gas turbine. The step by step technique will be used to simulate core design of the air cooled heat exchanger. In this method, the heat exchanger is divided into two dimensions (Depth, Height) ,as it will be described later. On other hand, Modeling of the computational program to predicate the effect of compressor inlet temperature on gas turbine performance by using the non-dimension method was also considered. 4.2 Heat Exchanger Thermal Design: 4.2.1 A Comprehensive Design Procedure The methodology of arriving at an optimum heat exchanger design is a complex one. Not only because of the arithmetic involved, but it is more particularly because of the many qualitative judgments that must be introduced. The design procedure in a schematic presentation is shown in Fig (4.1) Kays and London (1984) [26]. The design theory procedure can be set-up on a computer program. The inputs to the design theory procedure include: 1. Surface Characteristics: flattened tubes, surfaces with flow normal to banks of smooth tubes, Finned-tube surfaces, normal distance (XT), longitudinal distance (XL) …etc. 2. Problem Specifications: The problem statement may specify a consideration of different exchangers. For instance, periodic-flow and directtransfer types. Like cross flow, inline tubes, both fluids unmixed flow …etc. 3. Physical Properties: Some options may be allowed in the physical properties the matrix material to be used in a periodic flow type exchanger.

Chapter four…………………………………………………………………………Theory

Surface Characteristics

Problem Specifications

Design Theory Procedure

Physical Properties

Optional Solutions

Evaluation Criteria

Evaluation Procedure

Optimum Solution

Figure 4.1 Methodology of Heat Exchanger Design

4. Optional Solutions: optional solutions may represent an estimate of what a competitor may offer, others may represent customer's suggestions. For example, what are exit parameters from heat exchanger, heat transfer coefficient, heat load …etc.

Chapter four…………………………………………………………………………Theory

5. Evaluation Procedure: The theoretical design must be furnished with evolutions criteria to obtain optimum solutions for design theory.

4.2.2 Numerical Modeling of Cross Flow Compact Heat Exchanger In order to develop a numerical model with the predictive capability for various design parameters of the heat exchanger, step by step method with two dimensions was employed in this study. Each tube row was divided into a number of horizontal slices occupying the total length of the heat exchanger and for air side the exit of one row considering inlet to next row as shown in Fig. (4.2).Step by step method enables us to take into account the significant air temperature increase as well as the local variations of the properties and the heat transfer coefficient. Forms of the mass and the energy conservation equations were derived for two dimensional grid systems. 4.2.2.1 Grid System In the present study the following assumptions were assumed : 1. Homogenous temperature distribution of air all over the frontal face area of the heat exchanger and hence for each slice. 2. Uniform mass flow distribution for both stream sides of heat exchanger. 3. The exit air condition for each row represents a mean value for all of the slices of the considered row. This will be the inlet condition for the next row. 4. The inlet air velocity for each row was assumed to be uniform represented by a mean representative value. 5. The water temperature variations between the rows were also assumed to be negligible. 6. The air velocity stream and maximum air velocity difference assumed to be negligible. In the present study the design requirements for the theoretical model are: 1. The velocity in tube side (vw) 2. The velocity of air side (va.)

Chapter four…………………………………………………………………………Theory

Figure 4.2 Step by Step method with two Directions

Chapter four…………………………………………………………………………Theory

3. The air temperature difference (∆Ta) per row 4. The air mass flow for each slice. Therefore only a single fin tube of one horizontal slice assembly was explained as shown in Fig. (4.3a) and the heat exchanger performance for any locate can be calculated the performance for each slice and calculating performance for each row Fig.(4.3b). The nodal points for the calculation of the variables of the air and the coolant were defined as illustrated in Fig. (4.3c). Nodal points for air temperature were assigned on the east and the west sides of the control volume and the nodal points for the coolant were assigned on the north and the south sides of the control volume. The nodes were evenly distributed in two directions Ni × Nj. The heat transfer coefficient and heat load were calculated at each control volume.

Twi

Water Flow Tao

Lf Slice 1

Dt Tar2

Tar1

Slice 2

Tai

Two Air Flow

Figure (4.3a) Slice for inlet single Tube

Figure (4.3b) Exit of one row inlet to next row

Chapter four…………………………………………………………………………Theory

2 ………………Nr

N-1.………………….……2

I=1

J=1

Figure (4.3c) Nodal Points Distributions with two directions 4.2.2.2 Physical Characteristics of Heat Exchanger From the selected velocity of the tube side of the water stream and the known tube cross section dimension, the total number of tubes of the heat exchanger can be calculated from: .

Vw (Nt) = [ N r × v w × Acw

] ………………………………………………….(4-1)

This value will be used for the estimation of the length and depth of the heat exchanger as: L = X T × ( N t + 1) ………………………………………………….(4-2) D = [ X L × ( Nr − 1)] + Dt ………………………………………………….(4-3)

4.2.2.3 Mass Conservation Mass conservation of the water flow through the tube is simply .

.

Σ m in −Σ m out = 0 ………………………………………………….(4-4) and the mass conservation equation for the water at each slice can be written as:

Chapter four…………………………………………………………………………Theory .

.

m w(i+1,j)= m w(i,j) ………………………..…..………………………….(4-5)

In case of the heat exchangers with plain continuous fins, the mass flow rate in the air flow direction can be calculated by the following equation. .

.

.

m a(j)= ρ a ( j ) × V a ( j ) ……………………………………………………...(4-6)

4.2.2.4 Log-Mean Temperature Difference To estimate the true mean temperature difference ( ΔTm ) between the two fluids Fig. (4.4) shows the possible flow direction of both streams in compact heat exchanger.The following relations may be used for the estimation of the logarithmic mean temperature difference according to counter flow directions, Smith(1997) [26]: LMTD =

(Ta (i, j ) − Tw (i, j + 1)) − (Ta (i + 1, j ) − Tw (i, j )) ………………………….(4-7) (Ta (i, j ) − Tw (i, j + 1)) ln (Ta (i + 1, j ) − Tw (i, j ))

Ta,(i,j) ∆Tm Ta,(i+1,j)

Tw,(i,j+1)

∆Tm Tw,(i,j) 0

Atotal area

Figure 4.4 The Mean Temperature Difference Along a Single Pass

Chapter four…………………………………………………………………………Theory

The actual temperature difference of a cross flow compact heat exchanger is obtained by applying a correction factor (F) to the (LMTD) value as, Hewitt (1998) [27]: ΔTm = F × LMTD

…………………………………………………….(4-8)

where: F=

( R 2 + 1) ln (1 − S ) (1 − RS ) ( R − 1) ln[

2 − S[ R + 1 − ( R 2 + 1) ]

……………………………………….(4-8a)

2 − S [ R + 1 + ( R 2 + 1) ]

R=

(Th (i, j ) − Th (i + 1, j )) ……………………………………….(4-8b) (Tc (i. j + 1) − Tc (i, j ))

S=

(Tc (i, j + 1) − Tc (i, j )) ……………………………………….(4-8c) (Th (i, j ) − Tc (i, j ))

4.2.2.5 Heat Load The heat load passes through a control volume on the water side is: .

Q (i,j)= m (i,j) × Cp (i,j) × ΔT (i,j) …………………………………..…….(4-9) It can be expressed with the overall heat transfer coefficient, U(i,j), as follows: Q(i,j)= Uo (i,j) × A × (i,j) × ΔTm (i,j) ……………………………...…….(4-10) The above equation is applying for each slices per one row.

4.2.2.6 Overall Heat – Transfer Coefficient (Uo): The local overall heat transfer coefficient for each slice based on the outside tube area can be written as follows , Holman(1989) [28]:

U o (i, j ) =

1

tt 1 1 1 [( )+( )+( )+( )] ha (i, j ) × η o K t ( Ain Aout ) hw (i, j ) × ( Ai n Ao ut ) h f × ( Ai n Ao ut )

....(4 − 11)

Chapter four…………………………………………………………………………Theory

4.2.2.7 Forced Convection Heat Transfer Coefficient Inside Tube: Numerous relations have been proposed for predicting fully developed turbulent flow in tubes. The popular Dittus-Boelter equation, Dittus-Boelter (1930) [30] is usually given in the form: h w = 0 . 023 Re 0.8 Pr n (

k d hw

………………………………………... (4-12)

)

and cross flow area for water side (rectangular tube with semi circular ends ):

[

]

Ac w = [( H t − Tt ) × ( Dt − H t )] + ( H t − Tt ) 2 × π ……………………………..(4-13)

also Re =

dh =

Pr =

ρud h w ……………………………………..…………………….(4-14) μ 4 Acw

λt

μCp k

……………………………………..………………….….(4-15) ……………………………………..…………………....….(4-16)

where the ranges of Re and Pr are: 6000 ≤ Re ≤ 10 7 0.5 ≤ Pr ≤ 120

where the coefficient (0.023) is recommended by McAdams (1954)

[31]

in

place of (0.0243) originally given by Dittus-Boelter. Also, n= 0.4 for heating n= 0.3 for cooling 4.2.2.8 Forced Convection Heat Transfer Coefficient for Air Side: The entrance region for the development of the longitudinal velocity profile and the temperature profile is about 10 times the hydraulic diameter. This criterion is particularly valid for calculating the time-averaged coefficient for fluids (air and water) Adrian [32] . There are several empirical relationships for heat transfer between the duct surface and the fully developed flow, the

Chapter four…………………………………………………………………………Theory

analytical form of these relationships is based on exploiting the analogy between momentum and heat transfer. The popular Dittus-Boelter equation is also used for fully developed turbulent flow (air side): h a = 0 . 023 Re 0 .8 Pr n (

k d ha

………………………………………….. (4-17)

)

For laminar flow, the Sieder and Tate (1930) [33] correlation can be used. ha = 1.86(Re Pr )

0 .3

0.3

k ⎛ d ha ⎞ ( ) ⎟ ⎜ D f ⎠ d ………………………………….……..….….(4-18) ⎝ ha

where cross flow area for air side: Pf

Ac a = (

2

)( X T − H t ) ………………………………………………..….(4-19)

The hydraulic diameter for the air side: dha =

4 Aca

……………………………..………………….….(4-20)

λf

λ f =2 Lf + pf

The number of slices (N) depends on the water temperature difference along the tube bare; by assuming the air mass flow rate across frontal area is divided equally on the slices (N), therefore the Reynolds number for air side is represented by: .

ma ( j ) ma(i, j ) = …………………………………………………………..(4-21) N .

.

m (i, j ) ……………………………………..………………..……. (4-22) Ga = a ( Aca ) Ra (i, j ) =

(d h a Ga ) ………………………………..………………….… (4-23) μ a (i, j )

By assuming that there is no temperature gradient between the tubes that share the fins then the fin can be considered to be insulated at the center. Therefore fin efficiency (ηf) can be calculated as that for the case with an adiabatic tip, Briggs and Young (1963) [34]:

Chapter four…………………………………………………………………………Theory

ηf =

[tanh(m × ( (m × (

Lf 2

Lf 2

))]

………………………………..………………… (4-24)

))

where:

m= [

(ha (i, j ) × ( Pf ) ( K t × Aca )

………………………………..………………..… (4-25)

The overall surface fin efficiency is: η o = 1 − [(

Af A

) × (1 − η f )] ………………………………..………………… (4-26)

where: Af H

Aexp H

= (4 × L f × D f )(

1 ) ………………………………..……………….. (4-27) Pf

= [((π × H t ) + (2 × ( Dt − H t )) − (2 × T f )](

1 ) ………..………………… (4-28) Pf

Af Aexp A =( )+( ) ………..…………………………………. ………….. (4-29) H H H

The log- mean temperature difference and overall heat transfer coefficient is calculated to obtain the height for one slice as follow:

H (i, j ) =

Q(i, j ) (U o (i, j ) × A × F × LMTD (i, j ))

…………………………………….(4-30)

4.2.2.9 Power of Fan: The pressure drop across heat exchanger can be calculated from the general pressure drop relationship, London(1983)[35] is most often written in terms of hydraulic diameter, 2 Df va ΔP = ( ρ )× f × 2 d ha

…………………………………………………..(4-31)

Chapter four…………………………………………………………………………Theory

To calculate ΔP, the friction factor ƒ must be known and it can be derived from the flow solution.The friction factors derived from the Colburn [35] flows described by eqs. (3.32) for fully developed laminar flow: ƒ=

64 …………………………………………………………………...(4-32) Ra

For isosceles triangular ducts, Bhatti and Shah(1987)[36] recommend, for fully-developed turbulent flow: ƒ=

0.078 ……………………………………………………………….(4-33) ( Ra) 0.25

For the power of the fan the equation is: .

Pfan =

V air Δp

η fan

……………………..………………………………….….(4-34)

And, .

V air =

m& air

ρ air

……………………..………………………………..…….(3-35)

Thus,

Pmotor , fan =

Pfan

η motor , fan

……………………..…………………………….….(3-36)

4.2.2.10 The Computer Program The computer program was built in this study to establish the thermal design of the compact heat exchanger incorporating the quick-basic computer language. The flow chart of the program (CPHE) is shown in appendix (B).The following procedure describing the calculation technique of the present model : 1. Input surface characteristics and inlet operation conditions for both fluids. 2. Calculating the cross sectional area of tube side. 3. Choosing the velocity in tubes and input water flow rate V&w 4. Calculating (No. of Tubes / Row) ,Eq.(4-1)

Chapter four…………………………………………………………………………Theory

5. Calculating the length and depth of heat exchanger Eq.(4-2) and Eq.(4-3). 6. Calculating the No. of slices from water temperature conditions. 7. Loop (j=1…….j=no. of row) 8. Assume the air exit temperature for first row. 9. Calculating the air mass flow rate for first row, Eq (4-6). 10. Assume the air velocity over tubes. 11. Calculating the height of the heat exchanger. 12. Loop (I=.1……i=no. of slices). 13. Choosing the water temperature difference for each slice. 14. Calculating the heat load for first slice and first row, Eq.(4-9). 15. Assume the air mass flow rate for each slice, Eq.(4-21). 16. Calculating the air exit temperature for first slice. 17. Correct the fluids properties. 18. Calculating the correct heat load for first slice. 19. Calculating the correct air exit temperature for first slice. 20. Calculating the cross sectional area, hydraulic diameter, Reynolds k No., Nusselt No, Eq.(4-13) to Eq.(4-20). 21. Calculating the overall heat transfer coefficient, Eq.(4-11). 22. Calculating the height of heat exchanger, Eq.(4-30). 23. Repeat the calculation with the iterated value of air exit temperature until the error percent calculated from: ξ% =

(Varible cal. − Varible assu.)

Variable cal.

is converged to a value within 1 × 10 −3 .hence,

the variable here is air exit temperature or air velocity over tubes. 24. Repeat process for all slices (Ni). 25. Calculating the mean air exit temperature and check it with assumed value. 26. Repeat the process for all rows (Nj). 27. Calculating heat exchanger performance (heat transfer coefficient for air side (ha), overall heat transfer coefficient (Uo), heat load)

Chapter four…………………………………………………………………………Theory

28. Calculating height of heat exchanger (H). The step by step method is formulated in this program to establish the following characteristics and operating conditions of the heat exchanger: 1. Air exit temperature 2. Air mass flow 3. Heat transfer coefficient for both sides 4. Overall heat transfer coefficient 5. Heat load 6. The height of heat exchanger (H) 7. The length of heat exchanger (L) 8. The depth of heat exchanger (D) 9. Number of fins per tube 10. Number of tubes per row 11.In addition, heat exchanger performance can be predicted for any required point in heat exchanger in two direction (Ni X Nj). 4.3 Gas Turbine .

A conventional power plant receiving fuel energy (Qadd), producing work ( w ) and rejecting heat to a sink at low temperature (Qrej) is shown in Fig.( 4.5) as a block diagram.

Figure. 4.5 Basic Gas Turbine Engine

Chapter four…………………………………………………………………………Theory

The objective is to achieve the least fuel input for a given work output as this will be economically beneficial in the operation of the power plant. Usually, a gas turbine plant operates on ‘open circuit’, with internal combustion Horlock(1987)

[37]

. Air and fuel pass into the compressor and

combustion chamber, respectively, and the combustion products leave the gas turbine after expansion through the turbine. 4.4 Basic Gas Turbine Cycles In power plant thermodynamics for high thermal efficiency led us to emphasis on raising the maximum temperature T3 and lowering the minimum Temperature T1. Thus, the efficiency will be increased with the ratio (T3/T1). In a gas turbine plant, this search for high maximum temperatures is limited by material considerations and cooling of the turbine is required. This is usually achieved in ‘open’ cooling systems, using some compressor air to cool the turbine blades and then mixing it with the mainstream flow. In practical open circuit gas turbine plants with combustion, real gas effects are present (in particular the changes in specific heats, and their ratio, with temperature), together with combustion and duct pressure losses.

4.5 A non- Dimensional Analyses of Gas Turbine Performance: In practical open circuit gas turbine plants with combustion, real gas effects are present (in particular the changes in specific heats, and their ratio, with temperature), together with combustion and duct pressure losses. Now some modifications of the air standard analyses and their graphical presentations for such open gas turbine plants must be developed, as an introduction to more complex computational approaches.

Chapter four…………………………………………………………………………Theory

The Hawthorne and Davis(1956)

[38]

analysis is first generalized for the

[CBT] open circuit plant, with fuel addition for combustion, (ƒ) per unit air flow, changing the working fluid from air in the compressor to gas products in the turbine, as indicated in Fig.(4.6).

T

3 Combustion Fuel 4 2 2s

4s

(1+ƒ) 5

Gases Products

5s

1 S Figure 4.6 (T-s) Diagram for Irreversible Two-Shaft Circuit Simple Plant

Real gas effects are present in this open gas turbine plant; specific heats and their ratio are functions of ƒ and T, and allowance is also made for pressure losses. 4.5.1 Component Performance Before moving on to the air standard analyses of irreversible gas turbine cycle, need to be define various criteria for the performance of some components. In addition, to the irreversibilities associated with these components, pressure losses (Δp) may occur in various parts of the plant (in the entry and exit ducting, the combustion chamber, and the heat exchanger).

Chapter four…………………………………………………………………………Theory

These are usually expressed in terms of non-dimensional pressure loss coefficients, Ѕ= Δp/ (p)in ……………………………………………………………(4-37) where: (p)in=is the pressure at entry to the duct. ΔP23=

P2 − P3 ………………………………………………………..(4-37a) P2

ΔP5a =

P5 − Pa Pa

(Pa is ambient pressure)………………………….…(4-37b)

As alternatives to the isentropic efficiencies for the turbomachinery components, ηT, ηC, which relate the overall enthalpy changes, small-stage or polytropic efficiencies (ηpT and ηpC) are often used. The pressure-temperature relationship along an expansion line Fig.(3.6) is then P/ Tz = constant ………………………………………………………(4.38) where: Zg=γg / [(γg-1) × ηpT ] …………………………………………………..…(4.39) and the entry and exit temperatures are related by T4 / T5s =( rT)1/Zg = xt …………………………………………………..(4.40) Along a compression line, P/ Tz = constant …………………………………………………….….(4.41) where: Za=[γa × ηpC]/ (γa-1) ……………………………………………………...(4.42) and exit and entry temperatures are related by T2s / T1 =( rC)1/Za = xc ………………………………………………….(4.43)

Chapter four…………………………………………………………………………Theory

4.5.2 Graphical Plot In graphical interpretation, using isentropic rather than polytropic efficiencies,Hawthorne and Davis [38] plotted the following non-dimensional quantities, all against the parameter x = r (γ-1)/γ as follows:(i) Non-dimensional compressor work, Wc= cpa × T1 × (xc-1)/ηc …………………………………………………(4.44) (ii) NDCW = Wc/cp(T3 – T1) = =(Xc-1)/(ηc × (Ø-1)) ………………….....(4.45) Non-dimensional turbine work, WT=(1+ƒ) × cpg × T3 × ηt(1-1/xt) …………………………………………(4.46) (iii) NDTW = WT/cp (T3 – T1) ==((1+f) × ηt × (1-1/Xt))/((1-1/Ø) × n) .…(4.47) Non-dimensional overall efficiency,

ηo=NDTW /(H. V × (1+ƒ)) …………………………………...…………(4.48) where: Ø=T3/T1 ratio ……………………..……………………………………(4.49) n= cpa/ cpg ……………………………………………………………(4.50) 4.6 Computer Calculations for Two-Shaft Gas Turbine The analytical approach outlined above for the two-shaft gas turbine plants is that used in modem computer codes. However, gas properties, taken from tables such as those of Keenan and Kaye(1945) [39] as shown in fig. (4.7).The flow chart of the program (CPTGTP) is illustrated in Appendix [B].

Chapter four…………………………………………………………………………Theory

1

Figure 4.7 Specific Heats and Their Ratios for ‘Real’ Gases-Air and Products of Combustion

The purpose of the computer calculations is to examine the possibility of increasing the overall performance of two-shaft gas turbine depending on inlet conditions only. 4.7 Results and Discussion of Theoretical Calculations The dimensionless parameter such as overall efficiency, expansion ratio, and maximum to minimum ratio (Ø) can be used for predicting the effect of compressor inlet temperature on gas turbine performance. The prediction is focused on the following concepts: 4.7.1 The effect of compressor inlet temperature on maximum to minimum ratio (Ø): One of two operating limits for gas turbine it is the turbine inlet temperature. It is not often when the engine runs at topping temperature will produce higher overall efficiency, but it depends on the minimum temperature. Therefore, wherever the minimum temperature is reduced, the engine will be at higher efficiency. This is due to the fact that the equilibrium condition between the compressor power requirement (which increases at high minimum temperatures) and the power produced by gas generator

Chapter four…………………………………………………………………………Theory

turbine (which is not directly influenced by the minimum temperature) will be satisfied at a lower value. Figures (4.8) to (4.11) show the effect of the Ø on the overall efficiency for inlet high pressure turbine temperature of (699°C),(628.5°C),(588°C) and (547.5°C) respectively at expansion ratios of 1.08 and 1.12.It is clear from these figures that for a given expansion ratio the overall efficiency shows a gradual increase with Ø .Further, the trend shows a raise in the overall efficiency as the inlet turbine temperature increases for a given expansion ratio. It is obvious that the increasing percentage varies with expansion ratio or turbine inlet temperature. When turbine inlet temperature was (699°C) the increasing percentage of overall efficiency was about (18%) at expansion ratio (1.08), while the increasing percentage of overall efficiency fell to (15%) at expansion ratio (1.12).Also, for the turbine inlet temperature (547.5°C) the increasing percentage of overall efficiency was about (20%) at expansion ratio (1.08),while the increasing percentage of overall efficiency was reduced to about (17%) at expansion ratio (1.12). However, there is relationship between (T3), (T1), and engine speed, at full load. Two shaft engines will run either at temperature topping or at speed topping. At speed topping, the engine will not reach its full firing temperature, while at temperature topping; the engine will not reach its maximum speed. The effect of the compressor inlet temperature on the overall efficiency is shown in figures (4.12) and (4.13) for expansion ratios of 1.12 and 1.08 respectively. When reducing the compressor inlet temperature from (30°C) to (15°C), the overall efficiency shows an increase by (17%) at expansion ratio (1.12).Also, the overall efficiency shows an increase by (19%) at expansion ratio (1.08).Noting that the higher the expansion ratio, the higher the overall efficiency is obtained as shown above.

Chapter four…………………………………………………………………………Theory

4.7.2 The effect of compressor inlet temperature on the expansion ratio: Changes in the minimum temperature have less impact the component efficiencies than the overall cycle output.Therefor, the impact of minimum temperature is usually less pronounced for the expansion ratio than for the overall efficiency. It is known from gas turbine cycle when expansion ratio increases, the overall efficiency will be increased as fig.(4.14),and fig.(4.15). But as both turbine inlet temperature (T3 ) and Ø (maximum to minimum ratio) increased it was found that the overall efficiency will be increased too. From the mentioned conditions, it can be concluded that any gas turbine cycle shows an improvement in thermal efficiency as long as increasing the turbine inlet temperature (T3) or lowering compressor inlet temperature. As, the expansion ratio increases, the work consumed by the compressor is reduced. Figure (4.16) and figure (4.17) showed the effect of reduced compressor inlet temperature on the overall efficiency. It is obvious that this effect can be clearly seen when the engine runs at high speed. So, when reducing compressor inlet temperature from (30°C) to (15°C), the increasing percentage in overall efficiency is fluctuated among (15 %), (17%) and (20%) at expansion ratio (1.12), (1.1) and (1.08), respectively.

Chapter four…………………………………………………………………………Theory

Theoretical Overall Efficiency (%)

8

7

6

5

4 Expansion Ratio (rT)rT=1. =1.12 Expansion Ratio 12 3 Expansion Ratio (rT)rT=1. =1.08 Expansion Ratio 08 2 3.08

3.1

3.12

3.14

3.16

3.18 T3/T1 (

Ø

3.2

3.22

3.24

3.26

3.28

)

Fig (4.8) Theoretical Overall Efficiency as a Function of (T3/T1 Ratio) with Isentropic Efficiency (ηt,ηc =0.9) ,T3=699 °C

7

Theoretical Overall Efficiency (%)

6.5 6 5.5 5 4.5 4 3.5 Expansion Ratio =1.12 Expansion Ratio(rT) rrT=1. T=1. 12 12

3

Expansion Ratio(rT) rrT=1. T=1. 08 08 Expansion Ratio =1.08

2.5 2 3.08

3.1

3.12

3.14

3.16

3.18 T3/T1 (

3.2

3.22

3.24

3.26

Ø)

Fig (4.9) Theoretical Overall Efficiency as a Function of (T3/T1 Ratio) with Isentropic Efficiency (ηt,ηc =0.9) ,T3=628.5 °C

3.28

Chapter four…………………………………………………………………………Theory

Theoretical Overall Efficiency (%)

6.5 6 5.5 5 4.5 4 3.5 3

Expansion Ratio r(rT) =1.12 T=1. 12

2.5

T=1. 08 Expansion Ratio r(rT) =1.08

2 3.08

3.1

3.12

3.14

3.16

3.18

T3/T1 (

3.2

3.22

3.24

3.26

3.28

Ø)

Fig (4.10) Theoretical Overall Efficiency as a Function of (T3/T1 Ratio) with Isentropic Efficiency (ηt,ηc =0.9) ,T3=588 °C 6.5

Theoretical Overall Efficiency (%)

6 5.5 5 4.5 4 3.5 3 Expansion Ratio (rT) rT=1.=1.12 12 2.5

Expansion Ratio (rT) =1.0 rT=1. 08 8

2 3.08

3.1

3.12

3.14

3.16

3.18 T3/T1 (

3.2

3.22

3.24

3.26

3.28

Ø)

Fig (4.11) Theoretical Overall Efficiency as a Function of (T3/T1 Ratio) with Isentropic Efficiency (ηt,ηc =0.9) ,T3=547.5 °C

Chapter four…………………………………………………………………………Theory

7.5

Theoretical Overall Efficiency (%)

7 6.5 6 5.5 5 4.5 T1= 15 C

4 3.5

T1= 30 C

3 500

520

540

560

580

600

620

640

660

680

Turbine Inlet Temperature (T3)

Fig (4.12) The Effect of Turbine Inlet Temperature (High Pressure Turbine)on the Theoretical Overall Efficiency ,Expansion Ratio=1.12

Theoretical Overall Efficiency (%)

5.5 5 4.5 4 3.5 3 T1= 15 C 2.5 T1=30 C 2 500

520

540

560

580

600

620

640

660

680

Turbine Inlet Temperature (T3)

Fig (4.13) The Effect of Turbine Inlet Temperature (High Pressure Turbine)on the Theoretical Overall Efficiency ,Expansion Ratio=1.08

Chapter four…………………………………………………………………………Theory

7.5

Theoretical Overall Efficiency (%)

7 6.5 6 5.5 5 4.5 T3=669 C ,

Ø=3.2

4 T3=628.5 C,Ø =3.1 3.5

T3=588 C, Ø =2.9

3 1.075

1.08

1.085

1.09

1.095

1.1

1.105

1.11

1.115

1.12

1.125

Expansion Ratio (rT)

Fig (4.14) Theoretical Overall Efficiency as a Function of Expansion Ratio with Isentropic Efficiency (ηt,ηc =0.9) and T1=15°C

Theoretical Overall Efficiency (%)

6

5.5

5

4.5

4 T3= 699 C, T3=628.5 C,

3.5

T3=588 C,

Ø =3.10 Ø=2.97 Ø=2.84

3 1.075

1.08

1.085

1.09

1.095

1.1

1.105

1.11

1.115

1.12

1.125

Expansion Ratio (rT)

Fig (4.15) Theoretical Overall Efficiency as a Function of Expansion Ratio with Isentropic Efficiency (ηηt,ηc =0.9) and T1=30°C

Chapter four…………………………………………………………………………Theory

Theoretical Overall Efficiency (%)

8

7

6

5

4 Expansion Ratio (rT) =1.12 Expansion Ratio rT=1. rT=1.1212

Expansion Ratio rT=1. 1 Expansion Ratio (rT) =1.1

3

Expansion Ratio rrT=1. Expansion Ratio (rT) = 1.08 T=1. 08 08 2 10

15

20

25

30

35

Compressor Inlet Temperature (C)

Fig (4.16) The Effect of Compressor Inlet Temperature (with Variable Expansion Ratio) on the Theoretical Overall Efficiency, T3= 669 °C 7

Theoretical Overall Efficiency (%)

6.5 6 5.5 5 4.5 4 Expansion Ratio (rT) rT= =1.12

3.5 3

Expansion Ratio (rT) r =1.1 T=

2.5

rT= =1.08 Expansion Ratio (rT)

2 10

15

20

25

30

35

Compressor Inlet Temperature (C)

Fig (4.17) The Effect of Compressor Inlet Temperature (with Variable Expansion Ratio) on the Theoretical Overall Efficiency, T3= 628.5 °C

CHAPTER FIVE

Results & Discussion

Chapter Five .......................................................................................Results &Discussion

5.1 General Since gas turbine is an air-breathing engine, its performance is changed by anything that affects the density and/or mass flow of the air intake to the compressor. The volumetric flow is constant with any shaft speed, it is possible to increase mass flow rate by increasing air density. Air-cooled heat exchanger has been sitting prior to intake of compressor .In order to find out the effect of compressor inlet temperature on the gas turbine performance. These important variables are: power output, fuel mass flow rate, heat consumption, heat rate, overall efficiency, air mass flow. The design and thermal performance of air-cooled heat exchanger are also studied. A new design technique is suggested which divided the heat exchanger in two directions to get more accurate thermal performance: heat load, heat transfer coefficient (air side), overall heat transfer coefficient, air mass flow, air exit temperature from heat exchanger. The present study concentrates on getting variable air exit temperature, air mass flow rate, size (aspect ratio), with known inlet operation conditions for both fluids only. 5.2 Computational Model Results for Heat Exchanger The computational model results of the developed program shall be discussed, including heat load, heat transfer coefficient for air side (ha), and overall heat transfer coefficient for each row of the heat exchanger. Air exit temperature and air mass flow rate discharged from heat exchanger are to be entered to gas turbine which must be controlled according to the requirements of present study. The computational model has been fed with the same operating conditions as those of the experimental test rig. For the object of validity of the theoretical prediction it was decided to use two set of experimental tests for the present heat exchanger in this work. These were conducted at entering water temperature of (10 °C) and (50 °C) for water flow rate of (2000 l/h) at two different air flow rates of (500 cfm) and (2000 cfm), as shown in tables (5-1) to (5-4).For the performance simulation, these tests

Chapter Five .......................................................................................Results &Discussion

were achieved at the same entering air temperature to the heat exchanger of (32 °C) .In these tables, three different heat exchanger core sizes were obtained for the same operating conditions on both sides, air and tube, sides of the test section. 5.2.1 Heat Load for heat exchanger Figures (5.1a) and (5.1b) show the variation of heat exchanger load with the air flow velocity through the tube bank for water entering temperature of (10°C) and (50°C) respectively. The air velocity was ranged between (1.2 m/s) and (4.6 m/s) corresponding to (500 cfm) and (2000 cfm) air flow rates respectively. It is obvious that the heat load experience an increase as the air flow rate increase .This is due to the improvement of the overall heat transfer coefficient of the heat exchanger by increasing the air side heat transfer coefficient. The heat load was increased by 3-4 times when the velocity was raised from (1.2 m/s) to (4.6 m/s) respectively, at entering water temperature of 10°C (Cooling Mode) and 50°C (Heating Mode). The heat exchanger load was ranged between (1 kW) and (4.5 kW) for the whole range of air flow rate. The step by step simulation model shows a good agreement with the experimental data as shown in figure (5.1) for the whole range of air velocity. The predicted heat exchanger performance (Heat Load, Overall Heat Transfer Coefficient ...etc) fell within 5 % for most of the simulation range. The trend of the lines of heat load is the same for upper point for both cases, but the magnitude of heat load is having different values at low air velocities. For example, in figure (5.1a) at air velocity (1.5 m/s), the heat load is (1 kW) at entering water temperature of (10 °C) and figure (5. 1b) for entering water (50 °C), the heat load is (1.5 kW).Which explains the heat load behavior is affected by the fluid properties at high temperature. 5.2.2 Heat Transfer Coefficient for air side Since the heat transfer performance of the heat exchanger for (gas-liquid) type is dependent on heat transfer coefficient of air side .Figures (5.2a), and

Chapter Five .......................................................................................Results &Discussion

(5.2b), show a comparison for the heat transfer coefficient for air side between the theoretical prediction and the experimental data. These figures represent the behavior of the heat transfer coefficient variation with water flow rates of (2000 l/hr). It is obvious from these data that, the heat transfer coefficient values decrease at low air velocity. And as a matter of fact, it has lower value when water entering at (10 °C).The prediction of the present model shows that it is possible to obtain the higher temperature difference for the air side at the same heat exchanger load when using a fat heat exchanger. This is due to the higher overall heat transfer produced at the same air flow rate as shown in tables (5-1) to (5-4). The predicted values from the computational model for the heat transfer coefficient for air side along its depth (D) and along heat exchanger height (H) are shown in figures(5-3) and (5-4) for air flow rates of (2000) and (500) cfm at water entering temperature (10 0C) and (50 °C) respectively. It is obvious that the heat transfer coefficient of the air side nearly stayed unchanged ,and it is essentially a constant value. However, it is more pronounced when the heat exchanger geometry was deeper in the depth direction as shown in figure (5.5).This is for a heat exchanger having an overall dimension of (18 × 11 × 18) cm3, table (5-1).It shows that the air heat transfer coefficient is in the range between (94) to (96) W/m2.°C. 5.2.3 Overall Heat Transfer Coefficient Figures (5.6a), (5.6b) show the variation of the overall heat transfer coefficient with air flow velocity through heat exchanger at water entering temperature of (10 °C) and (50 °C) respectively. The predicted values from computational model for the overall heat transfer coefficient along its depth (D) and along heat exchanger height (H) are shown in figures(5-7) and (5-8) for air flow rates of (2000) and (500) cfm at water entering temperature (10 °C) and (50 °C) respectively. It is obvious that the overall heat transfer coefficient of the air side did not change much and it is essentially a constant

Chapter Five .......................................................................................Results &Discussion

value. However, it is more pronounced when the heat exchanger geometry was deeper in the depth direction as shown in figure (5.9).This is for a heat exchanger having an overall dimension of (18 × 11 × 18) cm3, table (5-1).It shows that the air overall heat transfer coefficient is in the range between (84) to (85) W/m2.°C.The overall heat transfer coefficient of the heat exchanger approaches (47) and (45) W/m2.°C at the maximum air velocity falling to a minimum at (1.2 m/s) where Uo is (15 W/m2.°C) and (30 W/m2.°C) for cooling and heating modes respectively. The simulation prediction and experimental data showed a good agreement with maximum descripancy of (5%) for the whole range of test conditions. 5.2.4 Air Temperature Distribution Figure (5.10) shows the variation of the experimental and present model of the exit air temperature out of the heat exchanger and the air velocity. For both of the cooling and heating modes of the air passing through the heat exchanger, the exit air temperature shows a trend of decreasing as the air velocity increases. This is due to the fact that when the air velocity increases causes an increase in the air side heat transfer coefficient (ha ) and the (Uo )value which in turn produces higher heat exchanger load. However, the exit air temperature per any row and slice is one of the dominant parameters which were focused in the present work. Figures (5.11) and (5.12) show the air temperature distribution along its depth (D) and along heat exchanger height (H) for air flow rates of (2000 cfm) and (500 cfm) at water entering temperature (10 °C) and (50 °C). Figure (5.13) shows the temperature variation with row number with the leaving side for each row for a heat exchanger of 6 rows having (18 × 11 × 18) cm3. 5.2.5 The Effect of Core Aspect Ratio (H/L) and Size (L × D × H) The objective of this study is to develop a heat exchanger model which has the capability to predict the heat exchanger performance depending on the design parameters without relying on experimental data. The predictive

Chapter Five .......................................................................................Results &Discussion

capability of the model was demonstrated by studying two different cases. For the first case, the aspect ratio of the heat exchanger core was changed from (0.67 - 1) when air velocity (1.2 m/s). The water flow rate was kept constant at (2000 l/h). The aspect ratio is defined as the ratio of the height (H) to the length (L) of the core for the heat exchanger. Since the model is based on the down flow type with the water tanks on top and bottom of the core, larger aspect ratio means relatively longer tubes. For a given range of the aspect ratio, the results showed that a heat exchanger with smaller aspect ratio can perform better than that of larger aspect ratio case as shown in Figure (5.14). Pressure drop in air side (Δpa ) was increased by changing the aspect ratio from (0.67 - 1) with the different core size (L × D × H). This is because at higher aspect ratio, the frontal area decrease over air side experiencing more abstraction in flow. For the second case, pressure drop variation with the change of the core size was examined. The pressure drop was calculated for a given water flow rate (2000 L/h) and air velocity (1.2 m/s) . The effect of core downsizing on the pressure drop is presented in tables (5.1-5.4). As the core size is reduced the pressure drop rate is increased. This is another evidence of the importance of the predictive model which can properly reveal the effect of the core size (L × D × H) variation on the pressure drop in air side. 5.3 Experimental Results for Gas Turbine Changes in ambient temperature have an impact on full-load power and heat rate, and, also on part-load performance and optimum power turbine speed. Manufacturers typically provide performance maps that describe these relationships for International Organization for Standardization (ISO) conditions. The Design Point of the gas turbine engine is concerned with the following concepts: • Standard ambient conditions • Improved fuel type • Full gas turbine shaft speed

Chapter Five .......................................................................................Results &Discussion

• Full power output

• Minimum heat rate • Minimum heat consumption The percent of design considers the change ratio of performance information to gas turbine by changing location or operation conditions with respect to original design point, as shown in figure (1.1). Detailed discussion of the effect of compressor inlet temperature will be carried out in the following sections: 5.3.1 Effect of Compressor Inlet Temperature on Power Output The power output can be increased by increasing either air mass flow or fuel mass flow. However, increasing fuel mass flow is used according to the application but the power output can be increased by improving ambient conditions which attempts to make it near standard conditions. Power output versus turbine inlet temperature at various compressor inlet temperature is shown in fig.(5-15).The increasing of power output is due to increasing turbine inlet temperature and also because of reducing the compressor inlet temperature ,which increases the mass flow rate across the power turbine. An increase in power output about (22 %) was obtained due to reducing compressor inlet temperature from (30 - 15 °C) at gas generator speed of (65000 RPM). 5.3.2 Effect of Compressor Inlet Temperature on Fuel Mass Flow Rate Burning hydrocarbons (Kerosene-RT 10) with air leads to combustion gases that have practically the same gas constant as dry air. Thus the assumption that R = 287 J/kg.K is valid for both compressor and turbine. But, the isentropic exponents γC and γT are in reality not constant and they change significantly with temperature. Fuel mass flow rate versus compressor inlet temperature at various gas generator speed is shown in fig.(5-16).The increasing of fuel mass flow rate is due to increasing both ,gas generator speed and compressor inlet temperature. A decrease in fuel mass flow rate

Chapter Five .......................................................................................Results &Discussion

about (10 %) was obtained due to reducing compressor inlet temperature from (30 - 15 °C) at gas generator speed of (55000 RPM). 5.3.3 Effect of Compressor Inlet Temperature on Specific Fuel Consumption The fuel consumption is depended not only on the combustion chamber design, but also on required power output. The high pressure turbine work is accelerated by forced gas mass on it's blades. That means in decreasing compressor inlet temperature will increase gas generator speed without increasing in fuel mass flow. Also, specific fuel consumption is down to any power required. Specific fuel consumption versus compressor inlet temperature at various gas generator speed is shown in fig.(5-17).The increase in specific fuel consumption was due to the increase of both gas generator speed and compressor inlet temperature. A decrease in specific fuel consumption about (44 %) was obtained due to reducing compressor inlet temperature from (30 - 15 °C) at gas generator speed of (55000 RPM). 5.3.4 Effect of Compressor Inlet Temperature on Heat Consumption Work from a gas turbine can be defined as the product of mass flow, heat energy in the combusted gas (Cp), and temperature differential across the turbine. The mass flow in this equation is the sum of compressor airflow and fuel flow. The heat energy is a function of the elements in the fuel and the products of combustion. The power output from gas turbine can be increases by burning more fuel which means more heat consumption. Heat consumption versus compressor inlet temperature at various gas generator speed is shown in fig.(5-18).The increasing of heat consumption is due to increasing gas generator speed and compressor inlet temperature. A decrease in heat consumption of about (8 %) was obtained due to reducing compressor inlet temperature from (30 - 15 °C) at gas generator speed of (55000 RPM).

Chapter Five .......................................................................................Results &Discussion

5.3.5 Effect of Compressor Inlet Temperature on Heat Rate Heat Rate is the inverse of efficiency, which is an indication of the ratio between thermal energy, resulting from the combustion process, and mechanical energy, obtained on the power shaft. Heat rate versus turbine inlet temperature at various compressor inlet temperature is shown in fig.(519).The increasing of heat rate is due to increasing of both turbine inlet temperature and compressor inlet temperature. A decrease in heat rate of about (37 %) was obtained due to reducing compressor inlet temperature from (30 - 15 °C) at gas generator speed of (65000 RPM). 5.3.6 Effect of Compressor Inlet Temperature on Overall Efficiency Overall efficiency versus turbine inlet temperature at various compressor inlet temperature are shown in fig.(5-20).The increase in overall efficiency is due to the increase in turbine inlet temperature and also because of reduction of compressor inlet temperature. An increase in overall efficiency about (40 %) was obtained due to reducing compressor inlet temperature from (30 15 °C) at gas generator speed of (65000 RPM). 5.3.7 Effect of Compressor Inlet Temperature on Air Mass Flow Rate In general, for lowering compressor inlet temperature the air mass flow rate will be increased because of increasing air specific weight at low temperature. Air mass flow rate versus compressor inlet temperature at various gas generator speed are shown in fig.(5-21) and fig.(5-22).The increasing of air mass flow rate is due to the increase of gas generator speed and also because of reduction of compressor inlet temperature. An increase in air mass flow rate and air volumetric flow rate about (15 %) and (6 %) respectively were obtained due to reducing compressor inlet temperature from (30 - 15 °C) at gas generator speed of (55000 RPM). 5.3.8 Effect of Compressor Inlet Temperature on Pressure Ratio The pressure ratio of the compressor at constant speed becomes smaller with increasing compressor inlet temperature. There will be more work (or

Chapter Five .......................................................................................Results &Discussion

head) required to achieve a certain pressure rise. The increased work has to be provided by the gas generator turbine, and is thus a loss for the power turbine. Pressure ratio versus compressor inlet temperature at various gas generator speed are shown in fig.(5-23).The increasing of pressure ratio is due to the increase gas generator speed and also because of reduction of compressor inlet temperature. An increase in pressure ratio about (8 %) was obtained due to reducing compressor inlet temperature from (30 - 15 °C) at gas generator speed of (55000 RPM). 5.3.9 Effect of Compressor Inlet Temperature on Power Input to Compressor It is known the theoretical power input to compressor is equal to the theoretical power to high pressure Turbine. The compressor overall efficiency from isentropic efficiency and mechanical efficiency is increasing by decreasing compressor inlet temperature together with compression ratio is also increasing. Therefore, by reducing compressor inlet temperature the energy consumed to compressor decreasing and theoretical power work (same spool with high pressure gas turbine) increasing at the same gas generator speed as showed in Fig (5.24). On other hand the high air weight is needed little work to compress. Power input to compressor versus gas generator speed at various compressor inlet temperature is shown in fig.(5-24).The increasing of power input to compressor is due to the increase gas generator speed and also because of reduction of compressor inlet temperature. An increase in power input to compressor about (30 %) was obtained due to reducing compressor inlet temperature from (30 - 15 °C) for at gas generator speed (55000 RPM). 5.3.10 Effect of Compressor Inlet Temperature on Turbine Inlet Temperature The two parameters that play a main role in gas turbine design are pressure ratio and turbine inlet temperature for both gas generator and power generator

Chapter Five .......................................................................................Results &Discussion

engines.

When

decreasing

compressor

inlet

temperature,

both

the

compression ratio will be increased and air mass flow increasing too. On other hand turbine inlet temperature will be decreased. Therefore, the power output is increased by burning more fuel with increasing gas turbine speed as shown in the figure (5.25). Also an increase in high pressure turbine work of about (22 %) was obtained due to reducing compressor inlet temperature from (30 - 15 °C) for constant turbine inlet temperature at (650 °C) . 5.4 Percent of Design Design point performance is vital to the engine concept design process. The engine configuration, cycle parameters, component performance levels and sizes are selected to meet a given specification. This section describes this performance input, which cannot be divorced from component design. Design point performance must be defined before analysing of any of the other possible operating conditions. A number of key parameters that define overall engine performance are utilised to assess the suitability of a given engine design to the application. These engine performance parameters are described below to the present study which was taken from design performance calculations for gas turbine (GT-85) : • The full power output is 5 (kW)

• The maximum overall efficiency is 3.5 (%) • The minimum heat consumption The maximum heat consumption can be add for gas turbine (GT-85) as follow: heat consumption =

5 =142.8 (kW) 0.035

• The minimum heat rate The minimum heat rate can be calculated for (GT-85) as follow: heat rate=

1 =28.56 0.035

Chapter Five .......................................................................................Results &Discussion

Percent of Design =

Calculated Value at (T1 = 15 ) − Calculated Value at (T1 = 30 ) Design Value

× 100%

One main concept from present study is to attempt to work or operate near the above engine performance parameters of the design point as shown in figures (5.26) and (5.27). The change in ratio relation to original design (Percent of Design) is about about (15 %) increase in power output, (25 %) increase in overall efficiency and (10 %) reduction in heat consumption, these results were obtained with reducing compressor inlet temperature from (30 15 °C) at gas generator speed (55000 RPM). 5.5 Comparison between the Experimental and Theoretical Predictions of the gas turbine engine Figure (5-28) shows a comparison between experimental and theoretical predictions for overall efficiency of gas turbine engine (GT-85) with variable turbine inlet temperature. Both of the experimental data and the predicated values have the same trend of variation showering an increase in overall efficiency with raising turbine inlet temperature .The predicated values were higher than that the experimental data for the whole temperature range in the field between (550-675°C) about (18%) for the expansion ratio of (1.03). When the expansion ratio is raised then, the discrepancy percentage will be increased .However, the effect of irreversibility and operating with different conditions makes this discrepancy between the experimental and theoretical predictions for overall efficiency with the same turbine inlet temperature increasing. 5.6 Conclusion From literature survey and pervious discussion , it can be concluded that there is a relationship among compressor inlet temperature, turbine inlet temperature and gas generator speed. It is obvious that the effect of compressor inlet temperature on gas turbine performance is reduced at high

Chapter Five .......................................................................................Results &Discussion

gas generator speed. Therefore, it can be noticed that for a particular turbine inlet temperature such as (650 °C) fig. (5.15), the power output is increased of about (22%) due to reduction the compressor inlet temperature from (3015°C) at gas generator speed up to (65000 RPM). It can be also that an increase of about (44%) when reduction compressor inlet temperature from (30-15°C) at gas generator speed (55000 RPM).It is worth while mentioning that the industrial gas turbine is carried out by two types: power output with constant speed and power output with variable speed. The first type is popular using in the power plant (Power Generation, Co-Generation), the speed is important parameter which rely on it the current frequency, therefore the plant worked at full speed for all time. The second type is used to feed power for pumps and compressors (Oil, Gas). Therefore, the volumetric flow rate is incorporated with the speed of device which is operating at higher speeds.

47.28972

62.55761

62.80377

63.05259

63.30411

84.24339

84.56969

84.89945

85.56966

85.23277

85.91026

2.328

1.164

1.164

1.164

1.164

.776

.776

.776

.776

.776

.776

96.19064

95.78031

95.37463

94.97345

94.57677

21.49

23.59

25.69

27.69

29.79

19.39

21.49

23.59

25.69

27.69

29.79

23.59

10

10

10

10

10

10

10

10

12

12

12

12

12

12

12

12

12

12

12

12

.235

.237

.238

.239

.241

.242

.238

.239

.241

.242

.241

.242

6

5

4

3

2

1

4

3

2

1

2

1

Heat Exchanger (55 3.6 37) cm3 Heat Exchanger (28.5 7.7 27) cm3

Heat Exchanger (18 11 18) cm3

×

32

25.69

25.69

10

10

10

10

×

94.18447

69.22968

27.69

27.69

29.79

27.69

29.79

Core Size (L D H)

×

68.93851

29.79

32

29. 79

32

No. of ROW

×

68.65057

68.36579

50.81485

50.60406

(mbar)

∆pa

×

47.10219

Uo(w/m2. °C) ha w/m2. °C) Ta (in) °C Ta (out) °C Tw (in) °C Tw (out) °C

×

2.328

Q (KW)

Table (5-1) Air Temperature Distribution Along H.EX.Depth ,Vw=2000 (L/h),Va=2000 (cfm)

Chapter Five .......................................................................................Results &Discussion ×

×

45.31548

60.24699

60.69621

60.46152

60.46152

81.55405

81.84826

81.53701

81.82931

82.11871

81.81271

2.284

1.142

1.142

1.142

1.142

.761

.761

.761

.761

.761

.761

90.22915

90.59176

90.24413

89.89323

90.2618

40.49

38.39

36.29

34.19

32.09

42.59

40.49

38.39

36.29

34.19

32.09

38.39

50

50

50

50

50

50

50

50

48

48

48

48

48

48

48

48

48

48

48

48

.250

.248

.247

.246

.245

.243

.247

.246

.245

.243

.245

.243

6

5

4

3

2

1

4

3

2

1

2

1

Heat Exchanger (55 3.6 37) cm3 Heat Exchanger (28.5 7.7 27) cm3

Heat Exchanger (18 11 18) cm3

×

32

36.29

36.29

50

50

50

50

×

89.90861

65.76105

34.19

34.19

32.09

36.19

34.09

×

65.25092

32.09

32

34.09

32

Core Size (L D H)

×

65.51846

65.00219

48.59657

47.79097

No. of (mbar) ROW

∆pa

×

45.04536

Ta (out) °C Tw (in) °C Tw (out) °C

×

2.284

Q (KW) Uo (w/m2. °C) ha (w/m2. °C) Ta (in) °C

Table (5-2) Air Temperature Distribution Along H.EX.Depth ,Vw=2000 (L/h),Va=2000 (cfm)

Chapter Five .......................................................................................Results &Discussion ×

×

20.13704

18.50125

18.57706

18.65369

18.73119

25.65959

25.76418

25.86992

25.97685

26.08497

26.1943

.456

.2329

.2329

.2329

.2329

.155

.155

.155

.155

.155

.155

27.07296

26.95747

26.8433

26.7304

26.61876

21.57

23.67

25.78

27.89

29.99

19.46

21.57

23.67

25.78

27.89

29.99

23.67

10

10

10

10

10

10

10

10

10.4

10.4

10.4

10.4

10.4

10.4

10.4

10.4

10.4

10.4

10.4

10.4

6.15

6.18

6.22

6.25

6.29

6.32

6.22

6.25

6.29

6.32

6.29

6.32

6

5

4

3

2

1

4

3

2

1

2

1

Heat Exchanger (55 3.6 37) cm3 Heat Exchanger (28.5 7.7 27) cm3

Heat Exchanger (18 11 18) cm3

×

32

25.78

25.78

10

10

10

10

×

26.50834

19.21763

27.89

27.89

29.99

27.89

29.99

×

19.1368

29.99

32

29.99

32

×

19.05688

18.97783

20.68979

20.64679

Core Size (L D H) ×

20.09592

E-02

No. of (mbar) ROW

∆pa

×

.456

Q (KW) Uo (w/m2. °C) ha (w/m2. °C) Ta (in) °C Ta (out) °C Tw (in) °C Tw (out) °C

Table (5-3) Air Temperature Distribution Along H.X.Depth ,Vw=2000 (L/h),Va=500 (cfm)

Chapter Five .......................................................................................Results &Discussion

×

×

22.79861

25.26978

26.01909

25.5097

26.04461

35.53004

34.93548

35.29956

35.51754

35.57734

36.40247

.684

.342

.342

.342

.342

.228

.228

.228

.228

.228

.228

37.36954

36.50051

36.43759

36.2082

35.82522

78.37

69.88

54.46

43.02

34.12

88.12

78.37

69.88

54.46

43.02

34.12

69.88

50

50

50

50

50

50

50

50

49.5

49.5

49.5

49.5

49.5

49.5

49.5

49.5

49.5

49.5

49.5

49.5

7.58

7.25

6.98

6.75

6.57

6.41

6.95

6.78

6.57

6.41

6.57

6.41

6

5

4

3

2

1

4

3

2

1

2

1

Heat Exchanger (55 3.6 37) cm3 Heat Exchanger (28.5 7.7 27) cm3

Heat Exchanger (18 11 18) cm3

×

32

54.46

54.46

50

50

50

50

×

36.45073

26.53511

43.02

43.02

34.12

43.02

34.12

×

25.98009

34.12

32

34.12

32

×

26.0863

25.73128

23.4331

23.84225

Core Size (L D H) ×

23.18965

E-02

No. of (mbar) ROW

∆pa

×

.684

Q (KW) Uo (w/m2. °C) ha (w/m2. °C) Ta (in) °C Ta (out) °C Tw (in) °C Tw (out) 0C

Table (5-4) Air Temperature Distribution Along H.X.Depth ,Vw=2000 (L/h),Va=500 (cfm)

Chapter Five .......................................................................................Results &Discussion ×

×

Chapter Five .......................................................................................Results &Discussion

Table (5-5) Air Temperature Distribution a long First Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37 Tw. in = 10(°C) Tw.out =12(°C)

Ta .in =32(°C)

ha (w/m2.°C)

Q (watt)

Uo (w/m2.°C)

Ta out (°C)

A (m2)

50.60411

116.4624

47.09801

29.89904

.1187113

50.6041

116.4564

47.0985

29.89915

.1192767

50.6041

116.4504

47.09899

29.89925

.1198476

50.6041

116.4474

47.09924

29.89931

.1201254

50.6041

116.4474

47.09924

29.89931

.1201254

50.60409

116.4385

47.09997

29.89947

.1210062

50.60408

116.4326

47.10045

29.89958

.1215941

50.60408

116.4267

47.10094

29.89968

.1221879

50.60407

116.4208

47.10142

29.89979

.1227876

50.60406

116.4149

47.10191

29.89989

.1233933

50.60406

116.409

47.0239

29.9

.1240051

50.60405

116.4032

47.10287

29.9001

.1246231

50.60405

116.3974

47.10334

29.90021

.1252474

50.60405

116.3945

47.10359

29.90026

.1255509

50.60404

116.3858

47.10431

29.90042

.1265153

50.60403

116.3801

47.10478

29.90052

.1271591

50.60403

116.3743

47.10526

29.90063

.1278096

50.60403

116.3686

47.10573

29.90073

.1284669

50.60402

116.3629

47.1062

29.90083

.1291311

50.60402

116.3572

47.10667

29.90093

.1298023

Vw = 2000 l/h Va=2000cfm

Chapter Five .......................................................................................Results &Discussion

Table (5-6) Air Temperature Distribution a long Second Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37 Tw. in= 10(°C) Tw.out =12(°C)

Ta .in =32(°C)

ha (w/m2. °C)

Q (watt)

Uo (w/m2. °C)

Ta out (°C)

A (m2)

50.81489

116.4624

47.28552

27.799

. 129798

50.81488

116.4564

47.28601

27.7991

.131507

50.81488

116.4504

47.2865

27.79921

.132205

50.81488

116.4474

47.28675

27.79926

.1329107

50.81488

116.4474

47.28675

27.79926

.133253

50.81487

116.4385

47.28749

27.79943

.133253

50.81486

116.4326

47.28798

27.79953

.1343451

50.81486

116.4267

47.28846

27.79964

.1350742

50.81486

116.4208

47.28896

27.79975

.1358114

50.81485

116.4149

47.28944

27.79985

.1365567

50.81485

116.409

47.28992

27.79996

.1373105

50.81484

116.4032

47.29041

27.80006

.1380727

50.81483

116.3974

47.29089

27.80017

.1388436

50.81483

116.3945

47.29113

27.80022

.1396233

50.81482

116.3858

47.29185

27.80038

.140001

50.81482

116.3801

47.29234

27.80048

.1412096

50.81482

116.3743

47.29281

27.80058

.1420166

50.81481

116.3686

47.29329

27.80069

.142833

50.8148

116.3629

47.29376

27.80079

.1436589

50.8148

116.3572

47.29423

27.80089

.1444946

Vw = 2000 l/h Va=2000cfm

Chapter Five .......................................................................................Results &Discussion

Table (5-7) Air Temperature Distribution a long First Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37 Tw. in= 50(°C) Tw.out =48(°C) Ta .in =32(°C)

Vw = 2000 l/h Va=2000cfm

ha (w/m2. °C)

Uo (w/m2. °C)

Ta out (°C)

Q (watt)

A (m2)

48.69674

114.1267

45.37363

34.09845

.1494603

48.69672

114.1362

45.37796

34.09862

.1503539

48.69672

114.1457

45.38214

34.0988

.1512587

48.69672

114.1504

45.38418

34.09888

.1516965

48.69672

114.1504

45.38418

34.09888

.1516965

48.6967

114.1645

45.39007

34.09914

.1531029

48.69669

114.1739

45.39383

34.09932

.1540425

48.69669

114.1833

45.39747

34.09949

.1549942

48.69668

114.1926

45.40099

34.09966

.1559581

48.69667

114.2018

45.40438

34.09983

.1569344

48.69666

114.2111

45.0767

34.1

.1579234

48.69665

114.2203

45.41086

34.10017

.1589252

48.69664

114.2295

45.41394

34.10034

.1599401

48.69664

114.234

45.41545

4.10042

.1604302

48.69663

114.2477

45.41983

34.10067

.1620101

48.69662

114.2567

45.42263

34.10084

.1630657

48.69662

114.2657

45.42537

34.10101

.1641354

48.69661

114.2747

45.428

34.10117

.1652195

48.6966

114.2837

45.43056

34.10133

.1663181

48.69659

114.2926

45.43305

34.1015

.1674317

Chapter Five .......................................................................................Results &Discussion

Table (5-8) Air Temperature Distribution a long Second Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37 Tw. in = 50(°C) Tw.out =48(°C)

Ta .in =32(°C)

ha (w/m2. °C)

Q (watt)

Uo (w/m2. °C)

Ta out (°C)

A (m2)

48.49659

114.1226

45.19305

36.19837

.1674317

48.49654

114.1267

45.19515

36.19855

.1714795

48.49653

114.1362

45.19946

36.19872

.172653

48.49652

114.1457

45.20361

36.19881

.1738432

48.49652

114.1504

45.20564

36.19881

.1744162

48.49652

114.1504

45.20564

36.19907

.1744162

48.4965

14.1645

45.21149

36.19924

.1762752

48.4965

114.1739

45.21524

36.19941

.1775177

48.49649

114.1833

45.21885

36.19958

.1787783

48.49648

114.1926

45.22234

36.19976

.1800573

48.49648

114.2018

45.22573

36.19992

.1813552

48.49647

114.2111

45.229

36.2001

.1826723

48.49646

114.2203

45.23216

36.20026

.1840092

48.49645

114.2295

45.23523

36.20035

.1853661

48.49644

114.234

45.23672

36.2006

.1860176

48.49644

114.2477

45.24109

36.20076

.1881417

48.49643

114.2567

45.24387

36.20093

.1895614

48.49642

114.2657

45.24658

36.2011

.191003

48.49641

114.2747

45.2492

36.20126

.192467

48.49641

114.2837

45.25175

36.20142

.1939538

48.4964

114.2926

45.25422

36.20128

.195464

Vw = 2000 l/h Va=2000cfm

Chapter Five .......................................................................................Results &Discussion

Table (5-9) Air Temperature Distribution a long First Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37 Tw. in = 10(°C) Tw.out=10.4(°C) Ta .in =32(°C)

Vw =2000 l/h

Va=500cfm

ha (w/m2. °C)

Q (watt)

Uo (w/m2. °C)

Ta out (°C)

A (m2)

20.64677

116.4624

20.09579

29.89989

.2781911

20.64677

116.4564

20.09587

29.9

.2795178

20.64677

116.4504

20.09595

29.90011

.2808575

20.64676

116.4474

20.09598

29.90016

.2815096

Table (5-10) Air Temperature Distribution a long Second Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37 Tw. in = 10(°C) Tw.out=10.4(°C) Ta .in =32(°C)

Vw =2000 l/h

Va=500cfm

ha (w/m2. °C)

Q (watt)

Uo (w/m2. °C)

Ta out (°C)

A (m2)

20.73277

116.4624

20.17804

27.79996

.3081484

20.73277

116.4564

20.17811

27.80007

.3097861

20.73277

116.4504

20.17819

27.80018

.3114415

20.73276

116.4474

20.17823

27.80023

.3122447

Table (5-11) Air Temperature Distribution a long First Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37 Tw. in = 50(°C) Tw.out=49.5(°C) Ta .in =32(°C)

Vw =2000 l/h

Va=500cfm

ha (w/m2. °C)

Q (watt)

Uo (w/m2. °C)

Ta out (°C)

A (m2)

19.86854

114.1267

19.32309

34.09974

.3509701

19.86854

114.1362

19.32381

34.09991

.353089

19.86854

114.1457

19.32451

34.10009

.3552339

19.86854

114.1504

19.32485

34.10017

.356272

Chapter Five .......................................................................................Results &Discussion

Table (5-12) Air Temperature Distribution a long Second Row Height (No. of Slice in i-Direction), Core Size=55 × 3.6 × 37 Tw. in = 50(°C) Tw.out=49.5(°C) Ta .in =32(°C)

Vw =2000 l/h

Va=500cfm

ha (w/m2. °C)

Q (watt)

Uo (w/m2. °C)

Ta out (°C)

A (m2)

19.52759

114.1267

18.99886

43.02019

.3620837

19.52758

114.1362

18.99919

43.02028

.3668006

19.52758

114.1457

18.99919

43.02028

.3668006

19.52758

114.1504

19.00015

43.02053

.3632253

Chapter Five .......................................................................................Results &Discussion

6 Simulation Experiment

Heat Load (kW)

5

4

3

2

1

0 0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

Air Velocity (m/s)



5

Entering water at 10 °C Fig (5. 1a)

Simulation Experiment

Heat Load (kW)

4

3

2

1

0 0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

Air Velocity (m/s)



Entering water at 50 °C Fig (5. 1b) Fig (5.1) Comparison between the Experimental and Present Model for the Effect of Air Velocity on Heat Load at Water Flow Rate 2000 (L/h)

Chapter Five .......................................................................................Results &Discussion

Heat Transfer Coefficient .ha (W/m 2.C)

60

Simulation Experiment

50

40

30

20

10

0 0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

Air Velocity (m/s)



Heat Transfer Coefficient. ha (W/m 2.C)

60

Entering water at 10 °C Fig (5. 2a)

Simulation Experiment

50 40 30 20 10 0 0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

Air Velocity (m/s)



Entering water at 50 °C Fig (5. 2b) Fig (5. 2) Comparison between the Experimental and Present Model for the Effect of Air Velocity on Heat Transfer Coefficient ha, at Water Flow Rate 2000 (L/h)

Chapter Five .......................................................................................Results &Discussion

Heat Transfer Coefficient ha.(w/m2.c) 50.45

50.5

50.55

50.6

50.65

50.7

50.75

50.8

50.85

Slice Number (Bar Tube)

1

2 ROW 1 ROW 2 3

4



Entering water at 10 °C Fig (5. 3a)

Heat Transfer Coefficient ha.(w/m2.c) 48.35

48.4

48.45

48.5

48.55

48.6

48.65

48.7

48.75

Slice Number (Bar Tube)

1

2 ROW 1 ROW 2 3

4



Entering water at 50 °C Fig (5. 3b)

Fig (5. 3) Variation Heat Transfer Coefficient (ha) a long Heat Exchanger Height at Water Flow Rate 2000 (L/h), Air Flow Rate 2000 cfm

Chapter Five .......................................................................................Results &Discussion

Heat Transfer Coefficient ha.(w/m2.c) 20.6

20.62

20.64

20.66

20.68

20.7

20.72

20.74

Slice Number (Bar Tube)

1

2 ROW 1 ROW 2 3

4



Entering water at 10 °C Fig (5. 4a)

Heat Transfer Coefficient ha.(w/m2.c) 19.4 19.4 19.5 19.5 19.6 19.6 19.7 19.7 19.8 19.8 19.9 19.9

Slice Number

1

2

ROW 1 ROW 2

3

4



Entering water at 50 °C Fig (5. 4b)

Fig(5. 4) Variation Heat Transfer Coefficient( ha )a long Heat Exchanger Height at Water Flow Rate 2000 (L/h), Air Flow Rate 500 cfm

Chapter Five .......................................................................................Results &Discussion

Heat Transfer Coefficient (W/m2.c)

96.5 96 95.5 95 94.5 94 93.5 93 1

2

3

4

5

6

Row Number

Fig(5. 5) Variation Heat Transfer Coefficient (ha) a long Heat Exchanger Depth at Water Flow Rate 2000 (L/h), Water Entering Temp. 10 °C, Air Flow Rate 2000 cfm

Chapter Five .......................................................................................Results &Discussion

Overall Heat Transfer Coefficient (W/m 2.C)

50 Simulation 45

Experiment

40 35 30 25 20 15 10 5 0 0

1

2

3

5

4

Air Velocity (m/s)

Overall Heat Transfer Coefficient(W/m 2.C)



Entering water at 10 °C Fig (5. 6a)

50 Simulation 45

Experiment

40 35 30 25 20 15 10 5 0 0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

Air Velocity (m/s)



Entering water at 50 °C Fig (5. 6b) Fig (5. 6) Comparison between the Experimental and Present Model for the Effect of Air Velocity on Overall Heat Transfer Coefficient at Water Flow Rate 2000 (L/h)

Chapter Five .......................................................................................Results &Discussion

Overall Heat Transfer Coefficient (w/m2.c) 47

47.05

47.1

47.15

47.2

47.25

47.3

47.35

Slice Number (Bar Tube)

1

2

ROW 1 ROW 2

3

4



Entering water at 10 °C Fig (5. 7a)

Overall Heat Transfer Coefficient (w/m2.c) 45.1

45.15

45.2

45.25

45.3

45.35

45.4

Slice Number

1

2

ROW 1 ROW 2

3

4



Entering water at 50 °C Fig (5. 7b)

Fig(5. 7) Variation Overall Heat Transfer Coefficient a long Heat Exchanger Height at Water Flow Rate 2000 (L/h), Air Flow Rate 2000 cfm

Chapter Five .......................................................................................Results &Discussion

Overall Heat Transfer Coefficient (w/m2.c) 20.04

20.06 20. 08

20.1

20.12

20.14

20.16 20.18

20.2

Slice Number

1

2

ROW 1 ROW 2

3

4



Entering water at 10 °C Fig (5. 8a)

Overall Heat Transfer Coefficient (w/m2.c) 18.8

18.9

19

19.1

19.2

19.3

19.4

Slice Number

1

2 ROW 1 ROW 2 3

4



Entering water at 50 °C Fig (5. 8b)

Fig(5. 8) Variation Overall Heat Transfer Coefficient a long Heat Exchanger Height at Water Flow Rate 2000 (L/h), Air Flow Rate 500 cfm

Chapter Five .......................................................................................Results &Discussion

Overall Heat Transfer Coefficient(W/m2.c)

86.5 86 85.5 85 84.5 84 83.5 83 1

2

3

4

5

6

Row Number

Fig(5. 9) Variation Overall Heat Transfer Coefficient a long Heat Exchanger Depth at Water Flow Rate 2000 (L/h),Water Entering Temp. 10 °C, Air Flow Rate 2000 cfm

Chapter Five .......................................................................................Results &Discussion

27.95

Air Exit Temperature (C )

27.9 27.85 27.8 27.75 27.7 27.65

Simulation 27.6

Experiment 27.55 0

1

0.5

1.5

2

2.5

3

3.5

4

4.5

5

Air Velocity (m/s)



Entering water at 10 °C Fig (5. 10a)

44

Air Exit Temperature (C )

43 42 41 40 39 38 37 Simulation 36 Experiment 35 0

0.5

1

1.5

2

2.5

3

3.5

4

4.5

5

Air Velocity (m/s)



Entering water at 50 °C Fig (5. 10b) Fig (5. 10) Comparison between the Experimental and Present Model for the Effect of Air Velocity on Air Exit Temperature at Water Flow Rate 2000 (L/h)

Chapter Five .......................................................................................Results &Discussion

Exit Air Temperature (c) 26.5

27

27.5

28

28.5

29

29.5

30

30.5

Slice Number (Bar Tube)

1

2

ROW 1 ROW 2

3

4



Entering water at 10 °C Fig (5. 11a) Exit Air Temperature (c)

33

33.5

34

34.5

35

35.5

36

36.5

Slice Number (Bar Tube)

1

2 ROW 1 ROW 2 3

4

Entering water at 50 °C Fig (5. 11b) Fig (5. 11) Variation Exit Air Temperature (°C) a long Heat Exchanger Height at Water Flow Rate 2000 (L/h), Air Flow Rate 2000 cfm •

Chapter Five .......................................................................................Results &Discussion

Exit Air Temperature (c) 26.5

27

27.5

28

28 .5

29

29.5

30

30.5

2 ROW 1 ROW 2 3

4



Entering water at 10 °C Fig (5. 12a)

Exit Air Tem perature (c) 30

32.5

35

37.5

40

42.5

45

1 Slice Number (Bar Tube)

Slice Number (Bar Tube)

1

2

Row 1 Row 2

3

4



Entering water at 50 °C Fig (5. 12b)

Fig(5. 12) Variation Exit Air Temperature) a long Heat Exchanger Height at Water Flow Rate 2000 (L/h), Air Flow Rate 500 cfm

Chapter Five .......................................................................................Results &Discussion

35

25

20

15

10

5

0 1

2

3

4

5

6

Row No.ofNumber Row Fig(5. 13) Variation Air Exit Temperature a long Heat Exchanger Depth atWater Flow Rate 2000 (L/h), Water Entering Temp. 10 °C,2000 cfm

43 Water Inlet Temperature 10 C

40 Pressure Drop in Air Side (mbar)*10-2

Air Exit Temperature (c)

30

37 Water Inlet Temperature 50 C

34 31 28 25 22 19 16 13 10 0.6

0.7

0.8

0.9

1

Aspect Ratio (H/L)

Fig (5. 14) The Effect of Aspect Ratio (H/L) with Different Core Size (L × D × H) on the Pressure Drop in Air Side, Water Flow Rate 2000 (L/h) Air Flow Rate 500 cfm

1.1

Chapter Five .......................................................................................Results &Discussion

2.5

Power Output (kW)

2

1.5

1

0.5

T1=15 C T1=30 C

0 610

620

630

640

650

660

Turbine Inlet Temperature (C )

Fig (5.15) The Effect of Turbine Inlet Temperature on the Power Output With Variable Compressor Inlet Temperature, at Rang gas Generator Speed (50000-65000 RPM)

2.6

Fuel Mass Flow Rate (g/s)

2.5

2.4

2.3

2.2

45000 RPM

2.1

55000 RPM 2 10

15

20

25

30

35

Compressor Inlet Temperature (C )

Fig (5.16) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the Fuel Mass Flow Rate

Chapter Five .......................................................................................Results &Discussion

Specific Fuel Consumption (g/kW.s)

2 1.8 1.6 1.4 1.2 1 45000 RPM

0.8 55000 RPM

0.6 10

15

20

25

30

35

Compressor Inlet Temperature (C )

Fig (5.17) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed ) on the Specific Fuel Consumption

110

Heat Consumption (kW)

108 106 104 102 100 45000 RPM 98 55000 RPM 96 10

15

20

25

30

35

Compressor Inlet Temperature (C )

Fig (5.18) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the Heat Consumption

Chapter Five .......................................................................................Results &Discussion

100 90

Heat Rate (kj/kW.s)

80 70 60 50 40 30 20

T1=15 C

10 0 610

T1=30 C 620

630

640

650

660

Turbine Inlet Temperature (C)

Fig (5.19) The Effect of Turbine Inlet Temperature (with Variable Compressor Inlet Temperature) on the Heat Rate, at Rang Gas Generator (50000-65000 RPM)

2 1.8

Overall Efficiency (%)

1.6 1.4 1.2 1 0.8 0.6 0.4

Ti=15 C

0.2 0 610

T1=30 C 620

630

640

650

660

Turbine Inlet Temperature (C )

Fig (5.20) The Effect of Turbine Inlet Temperature (with Variable Compressor Inlet Temperature) on the Overall Efficiency (%),at Rang Gas Generator Speed (50000-65000 RPM)

Chapter Five .......................................................................................Results &Discussion

0.14 0.13 Air Mass Flow Rate (kg/s)

0.12 0.11 0.1 0.09 0.08 0.07 0.06

45000 RPM

0.05

55000 RPM

0.04 10

15

20

25

30

35

Compressor Inlet Temperature (C )

Fig (5.21) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the Air Mass Flow Rate

450

400

cfm

350

300

45000 RPM

250

55000 RPM

200 10

15

20

25

30

35

Compressor Inlet Temperature (C )

Fig (5.22) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the Air Flow Rate

Chapter Five .......................................................................................Results &Discussion

1.44

Compression Ratio

1.39 1.34 1.29 1.24 45000 RPM

1.19 55000 RPM

1.14 10

15

20

25

30

35

Compressor Inlet Temperature (C )

Fig (5.23) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the Compression Ratio

6

Compressor Work (kW)

5 4 3 2 45000 RPM

1

55000 RPM

0 10

15

20

25

30

35

Compressor Inlet Temperature (C )

Fig (5.24) The Effect of Compressor Inlet Temperature (with Variable Gas Generator Speed) on the Compressor Work

Chapter Five .......................................................................................Results &Discussion

9

HP.Turbine Work (kW)

8 7 6 5 4 3 2 T1=15 C 1 0 610

T1=30 C 620

630

640

650

660

Turbine Inlet Temperature (C )

Fig (5.25) The Effect of Turbine Inlet Temperature (°C) on the HP.Turbine Work with variable Compressor Inlet Temperature (°C) , at Rang Gas Generator Speed (50000-65000 RPM)

Chapter Five .......................................................................................Results &Discussion

80 70

Percent of Design

60 50 40 30 POWER OUTPUT 20 OVERALL EFFICIENCY 10

HEAT CONSUMPTION

0 10

15

20

25

30

35

Compressor Inlet Temperature (C)

Fig (5.26) The Influence Compressor Inlet Temperature (°C )on Gas Turbine Performance at (Gas Generator Speed 45000 RPM)

80 70

Percent of Design

60 50 40 30 POWER OUTPUT

20 OVERALL EFFICIENCY

10

HEAT CONSUMPTION

0 10

15

20

25

30

35

Compressor Inlet Temperature (C )

Fig (5.27) The Influence Compressor Inlet Temperature (°C )on Gas Turbine Performance at(Gas Generator Speed 55000 RPM)

Chapter Five .......................................................................................Results &Discussion

2.8 2.6

Overall Efficiency (%)

2.4 2.2 2 1.8

Theoretical,Expansion Ratio=1.03

1.6

Experimental,Expansion Ratio=1.03

1.4 1.2 500

525

550

575

600

625

650

675

700

HP.Turbine Inlet Tem perature (C)

Fig (5.28) Comparison between the Experimental and Theoretical Predictions of the gas turbine engine (GT-85)

CHAPTER SIX

Conclusions & Recommendations

Chapter Six ..................................................................... Conclusions &Recommendations

6.1 Conclusions: The main goal of this work was to examine the possibility of increasing the power output and overall efficiency of gas turbine engine. Also, it was aimed to reduce heat rate. This was accomplished by reducing compressor inlet temperature. The following major conclusions are drawn from this work: 1. A computational model for design of air-cooled heat exchanger has been developed. For validation of the model, the heat load performance of a typical air-cooled heat exchanger was simulated over wide ranges of the air velocity and water flow rate .These were compared with the experimental data provided by experimental work. The specifications of the test heat exchanger are given in [Table (3.1)] shows a good agreement with the predicted operating conditions of the present model. The maximum discrepancy between the experimental data and calculated results for overall heat transfer coefficient and heat load was about (5%) for the given range of the simulated conditions. The model evaluation is obtained by checking its validity against experimental results obtained in the present work. 2. Increasing the ambient temperature lowers the density of the compressor inlet air. Thus reducing the mass flow through the turbine, and therefore will be reduced the power output. When the volume flow remains approximately constant, the mass flow will increase with decreasing temperature and will be decreased with increasing temperature. Therefore, when the compressor inlet temperature reduced from (30 °C) to (15 °C), the percent of design increase up to (15%) of the power output ,and ,the overall efficiency increased up to (25%) ,while the heat consumption reduction was about (10%). 3. The prediction of the gas turbine performance by using non-dimensional analysis is a good method which has been used in modern computer codes. The results revealed that the discrepancy percentage between theoretical and experimental prediction for gas turbine performance was about (18%).

Chapter Six ..................................................................... Conclusions &Recommendations

4. There is a relationship among the compressor inlet temperature, gas generator speed and turbine inlet temperature. At the higher speeds the increase in percentage of power output is lower than that of lower speeds when reducing the compressor inlet temperature. 6.2 Recommendations: 1. At extreme high climate temperatures, the cooling of the ambient temperature technique is preferable and worth doing where the benefits are high. On the other hand, the benefits are low in regions with mild or lowround temperatures, and that what has been detected through the experimental data. 2. Using the computational model program to design and study the effect of varying different design parameters: such as (number of rows, fin spacing, air velocity, velocity in tube, number tubes. For air-cooled heat exchanger and entering these parameters as data to intake of gas turbine is good idea, therefore by using different of inlet cooling system like absorption refrigeration cycle, evaporative cooling especially in low relative humidity regions is benefits. 3. Modifying the computational model program to investigate the gas turbine performance and study the effect of compressor inlet temperature on the performance using non-dimensional analysis, this requires further future investigation by other researchers.

References

1-Philip, K., "Power Generation Hand Book", 1st edition, McGraw-Hill Professional, ISBN: 0071396047, pp.115-117, 2002. http:// www.digitalengineeringlibrary.com 2-ISO , Gas Turbine Acceptance Tests, ISO 2314, International Organisation for Standardisation, Geneva,(1973). 3-Frank ,J. B.," GE Gas Turbine Performance Characteristics", GE Power Systems Schenectady, NY, GER-3567H, pp.6-8, 2002. 4-General Electric Energy (Oil & Gas)," MS 5001 Gas Turbine" Technical Training. 5-International Energy Outlook, DOE/EIA-0484, 2004. 6-Ganapathy,V., "Process-Design Criteria of Air Cooled Heat Exchangers" ,Chemical Engineering,pp.418-425,McGraw-Hill Publication Book,Co.,Newyork,1979. 7-Hedderich, C.P., and Kellerher, M.D. "Design and Optimization of Air Cooled Heat Exchangers", Journal of Heat Transfer, vol. 104 pp683-890, Nov. (1982). 8-Zhang, C., "Numerical Modeling Using a Quasi-Tree Dimensional Procedure for Large Power Plant", Trans. of ASME Journal of Heat Transfer, v.116, p.841, 1994. 9-Matthew, S., Layton and Joseph O'Hagan, "Comparison of Alternate Cooling Technologies for California Power Plants", Electric power research institute (2002). 10-Dohoy, J. , Dennis, N.," Numerical Modeling of Cross Flow Compact Heat Exchanger with Louvered Fins using Thermal Resistance Concept" SAE.Paper , 2006-01-0726, The University of Michigan. 11-Tarrad, A. H., Khudor, D. S., and Wahed, M.A.," A Simplified Model for the Prediction of the Thermal Performance for Cross Flow Air Cooled Heat Exchangers with a New Air Side Thermal Correlation ", Journal of Engineering and Development, vol. 12,No.3, 2008.

12-De Lucia ,M., Lanfranchi, C., and Boggio, V., “Benefits of compressor inlet air cooling for gas-turbine cogeneration plants”, Proceedings of the International Gas Turbine and Aero-engine Congress and Exposition, Houston Texas, 5–8 June 1995. 13-Saleh, F. A., "The Influence of Water Injection on Two Shaft Gas Turbine Performance", M.Sc.Thesis, AL-Mustansiria University College of Engineering, 1996. 14-Ait-Ali,M. A., “Optimum power boosting of gas-turbine cycles with compressor inlet air refrigeration.”, Journal of Engineering for Gas Turbines and Power, Vol. 119, pp. 124–129, 2001. 15-GE Nuovo Pignone Internal DT-'N',"Gas Turbine Performance Curves",Hyundi Engineering and Construction CO, Commessa-Job 170552021,2002. 16- Donald C. Erickson and Icksoo Kyung, “Aqua Absorption Turbine Inlet Cooler”, ASME International Mechanical Engineering Congress and Exposition, Draft IMECE2003-42870,pp.113-115. 17- Hameed,N.,"The Influence of Water Injection on Two-Shaft Gas Turbine Performance with Regeneration", M.Sc.Thesis, AL-Mustansiria University College of Engineering, 2004. 18-Benjalool,A.,"Evaluation of Performance Deterioration on Gas Turbines due to Compressor Fouling",MSc Thesis,School of Engineering,Academic Year, pp.46-48,2006,. 19-Tony Giampaolo, " Gas Turbine Handbook: Principles and Practices" 3rd Edition, Published by the Fairmont Press, Inc.700 Indian Trail, 2006. 20-Kuamit, A.A.," Design of Combined Cycle Power Plant and Air Cooling System", Ph.D,Thesis, AL-Mustansiria University College of Engineering, 2006.

21-Tarrad, A. H., and Mohammed, A. G., " A Mathematical Model for Thermo-Hydraulic Design of Shell and Tube Heat Exchanger Using a Step by Step Technique" Engineering and Development Journal, Vol. 10, No.4, December 2006. 22-Watlow products and technical support delivered worldwide," Diesel and Gas Turbine Temperatures sensors" Watllow electric Manufacturing Company,2001. 23-Gilbert Gilkes and Gordon LTD, "Technical Handbook of Two-Shaft Gas Turbine " , Gilkes, Kendal, England, 1978. 24-Philip P., Walsh," Gas Turbine Performance", Blackwell Publishing, 2nd Edition, 143 pp, 2004. 25-H., Cohen, G., F., C., Rogers and H., I., H., Saravanamuttoo," Gas Turbine Theory", 4th ed, 1996. 26-Kays, W.M., and London, A.L., "Compact Heat Exchangers", 3rdedition, McGraw-Hill Book Company, 1984. 27-Smith, E.M.," Thermal Design of Heat Exchangers ,a Numerical Approach", John Wiley and Sons, New York, 1997. 28-G.,F.,Hewitt, " Heat Exchanger Design Hand book",Begell House, New York, 1998. 29- Holman, J.,P., "Heat Transfer" 7th edition McGraw-Hill, New Yourk,1989. 30- Dittus, F.W., and Boelter, L.M.K., Univ. Calif. (Berkeley) Pub. Eng., Vol.2, pp.443, 1930. 31-McAdams, W.H., "Heat Transmission", 3rd, Ed., McGraw-Hill, New York, 1954. 32- Adrian, B.," Forced Convection: Internal Flows", Mechanical Engineering and Materials Science, Chap. 5,420 pp. 33- Sieder, E.N., and Tate, C.E., "Heat Transfer and Pressure Drop of Liquids in Tubes", Ind. Eng. Chem., Vol.28, pp.1429, 1930.

34- Briggs D.E. and Young, E.H., "Convection Heat Transfer and Pressure Drop of Air Flowing Across Triangular Pitch Banks of Finned Tubes", Chem. Eng. Prog. Symp. Ser. Vol.59, No.41, pp1-10, (1963). 35- London,A. L.,"Compact Heat Exchanger-Design Methodology", ed.Kakac,S.,Shah,R.K.and Bergles,A.E.,Hemisphere,New York,(1983). 36-Bhatti,M.S. and Shah,R.K.,"Turbulant and Transition Flow Forced Convective Heat Transfer in Ducts",Single Phase Convective Heat Transfer,Wiley,1987. 37-Horlock,J.H.,"Co-generation:Combined Heat and Power", Pergamon Press, Oxford, See also 2nd edn, Krieger, Melbourne, FL, 1987. 38-Hawthorne,W.R., and Davis, G.de V."Calculating gas turbine performance", Engineering 181, 361 -367, 1956. 39- Keenan, J.H. and Kaye, J.," Gas Tables", Wiley, New York, 1945.

,

Appendix (A) Experimental Work Data Tables

The experimental Data for gas turbine: Table (A-1) Ambient Temperature: 19.5 (°C) Ambient Pressure: 1.0133 bar 15 (°C) Compressor Inlet Temperature:

Engine Parameters

Gas Generator Speed

Gas Generator Speed

(45000 RPM)

(55000 RPM)

T1 (°C)

15

15

Δp (mmwg)

32

43

Vf

5

6

P2 (bar)

0.23

0.32

T 2 (°C)

46

52

T3 (°C)

561

577

P3 (bar)

0.22

0.31

T4 (°C)

519

531

P4 (mmbar)

35

45

T5 (°C)

485

492

Δp 4/5 (mmHg)

25

30

Torque (N.m)

0.92

0.93

Power Turbine

14400

16300

(L/h)

Speed (RPM)

Table (A-2) Ambient Temperature: 25 (°C) Ambient Pressure: 1.0133 bar Compressor Inlet Temperature: 20 (°C)

Engine Parameters

Gas Generator Speed

Gas Generator Speed

(45000 RPM)

(55000 RPM)

T1 (°C)

20

20

Δp (mmwg)

30

40

Vf

4.9

5.7

P2 (bar)

0.25

0.33

T 2 (°C)

57

63

T3 (°C)

588

604

P3 (bar)

0.23

0.31

T4 (°C)

546

558

P4 (mmbar)

37

50

T5 (°C)

506

517

Δp 4/5 (mmHg)

25

35

Torque (N.m)

1.3

1.36

Power Turbine

17100

21900

(L/h)

Speed (RPM)

Table (A-3) Ambient Temperature: 29 (°C) Ambient Pressure: 1.0133 bar Compressor Inlet Temperature:25 (°C) Engine Parameters

Gas Generator Speed

Gas Turbine Speed

(45000 RPM)

(55000 RPM)

T1 (°C)

25

25

Δp (mmwg)

40

50

Vf

4.7

5

P2 (bar)

0.29

0.34

T 2 (°C)

56

63

T3 (°C)

618

635

P3 (bar)

0.25

0.31

T4 (°C)

551

565

P4 (mmbar)

40

52

T5 (°C)

508

523

Δp 4/5 (mmHg)

25

33

Torque (N.m)

0.66

0.85

Power Turbine

17900

22000

(L/h)

Speed (RPM)

Table (A-4) Ambient Temperature: 29(°C) Ambient Pressure: 1.0133 bar Compressor Inlet Tempertaur:30 (°C) Engine Parameters

Gas Turbine Speed

Gas Turbine Speed

(45000 RPM)

(55000 RPM)

T1 (°C)

30

30

Δp (mmwg)

39

59

Vf

5

6

P2 (bar)

0.31

0.45

T 2 (°C)

64

73

T3 (°C)

631

650

P3 (bar)

0.28

0.35

T4 (°C)

552

578

P4 (mmbar)

38

59

T5 (°C)

511

531

Δp 4/5 (mmHg)

42

57

Torque (N.m)

0.71

1.9

Power Turbine

11800

21600

(L/h)

Speed (RPM)

Table (A-5) Ambient Pressure:1.0133 bar Turbine Inlet Temperature: 586 (°C) Running with

Running with

Compressor Inlet

Compressor Inlet

Temperature 15 (°C)

Temperature 22 (°C)

Engine Parameters

T3 (°C)

586

586

Gas Generator ( RPM) 55000

50000

Δp (mmwg)

43

47

Vf

6.2

5.8

P2 (bar)

0.32

0.35

T1 (°C)

15

22

T2 (°C)

52

60

P3 (bar)

0.3

0.32

T4 (°C)

539

534

P4 (mbar)

50

46

T5 (°C)

499

496

Δp 4/5 (mmHg)

32

30

Torque (N.m)

1.09

0.74

Power Turbine(RPM)

18800

16500

(L/h)

Table (A-6) Ambient Pressure:1.0133 bar Turbine Inlet Temperature: 590 (°C) Running with

Running with

Compressor Inlet

Compressor Inlet

Temperature 18 (°C)

Temperature 22 (°C)

Engine Parameters

T3 (°C)

590

590

Gas Generator ( RPM) 50000

45000

Δp (mmwg)

39

33

Vf

5.5

5

P2 (bar)

0.3

0.25

T1 (°C)

18

22

T2 (°C)

52

53

P3 (bar)

0.27

0.21

T4 (°C)

544

549

P4 (mbar)

45

40

T5 (°C)

504

509

Δp 4/5 (mmHg)

30

27

Torque (N.m)

1.16

1.33

Power Turbine(RPM)

19800

17900

(L/h)

Table (A-7) Ambient Pressure:1.0133 bar Turbine Inlet Temperature: 620 (°C) Running with

Running with

Compressor Inlet

Compressor Inlet

Temperature 20 (°C)

Temperature 30 (°C)

Engine Parameters

T3 (°C)

620

620

Gas Generator ( RPM) 55000

50000

Δp (mmwg)

44

40

Vf

5.2

4.7

P2 (bar)

0.34

0.29

T1 (°C)

20

30

T2 (°C)

58

60

P3 (bar)

0.31

0.25

T4 (°C)

550

554

P4 (mbar)

50

40

T5 (°C)

500

510

Δp 4/5 (mmHg)

32

25

Torque (N.m)

0.69

0.66

Power Turbine(RPM)

21200

17900

(L/h)

Table (A-8) Ambient Pressure: 1.0133 bar Turbine Inlet Temperature: 625 (°C) Running with

Running with

Compressor Inlet

Compressor Inlet

Temperature 15 (°C)

Temperature 30(°C)

Engine Parameters

T3 (°C)

625

625

Gas Generator ( RPM) 65000

50000

Δp (mmwg)

53

38

Vf

5.6

5.2

P2 (bar)

0.41

0.28

T1 (°C)

15

30

T2 (°C)

61

64

P3 (bar)

0.34

0.26

T4 (°C)

550

549

P4 (mbar)

55

35

T5 (°C)

508

500

Δp 4/5 (mmHg)

36

22

Torque (N.m)

0.59

0.36

Power Turbine(RPM)

21500

13100

(L/h)

Table (A-9) Ambient Pressure: 1.0133 bar Turbine Inlet Temperature: 635 (°C) Running with

Running with

Compressor Inlet

Compressor Inlet

Temperature 25 (°C)

Temperature 29 (°C)

Engine Parameters T3 (°C)

635

635

Gas Generator ( RPM) 57200

50000

Δp (mmwg)

41

39

Vf

5.2

4.7

P2 (bar)

0.34

0.31

T1 (°C)

25

29

T2 (°C)

62

61

P3 (bar)

0.31

0.26

T4 (°C)

564

566

P4 (mbar)

52

45

T5 (°C)

522

523

Δp 4/5 (mmHg)

33

29

Torque (N.m)

0.85

0.79

Power Turbine(RPM)

22200

19500

(L/h)

The experimental Data for heat exchanger Table (A-10) Vw=2000 L/h , Va=2000 Cfm

Twi (°C)

Two (°C)

Tai (°C)

Tao (°C)

50

48

32

36.7

40

38

32

35

20

22

32

27.7

15

17

32

27.5

10

12

32

28

Table (A-11) Vw=2000 L/h , Va=1000 Cfm

Twi (°C)

Two (°C)

Tai (°C)

Tao (°C)

50

49

32

36.2

40

38.9

32

35.9

20

21

32

28.2

15

15.9

32

28

10

10.9

32

27

Table (A-12) Vw=2000 L/h , Va=500 Cfm

Twi (°C)

Two (°C)

Tai (°C)

Tao (°C)

50

49.5

32

38.4

40

39.4

32

36.5

20

21.4

32

28

15

15.4

32

27.9

10

10.4

32

26.9

Table (A-13) Vw=3000 L/h , Va=2000 Cfm

Twi (co)

Two (co)

Tai (co)

Tao (co)

50

48.5

32

36.3

40

39.7

32

33

20

21.5

32

27.9

15

16.5

32

27.7

10

11.5

32

27

Table (A-14) Vw=3000 L/h , Va=1000 Cfm

Twi (co)

Two (co)

Tai (co)

Tao (co)

50

49.4

32

36

40

39.3

32

35

20

20.8

32

27

15

15.9

32

26.5

10

11.7

32

26.8

Table (A-15) Vw=3000 L/h , Va=500 Cfm

Twi (co)

Two (co)

Tai (co)

Tao (co)

50

49.8

32

35.8

40

39.7

32

35.1

20

20.4

32

27

15

15.4

32

26.9

10

10.4

32

26.4

Appendix (B) Flow Charts and Computer Programs

Flow Chart for Theoretical heat exchanger Performance Prediction(CPHE)

Start

Input Surface Characteristics & Inlet Operation Conditions for Both Fluids

Choose the Velocity in Tubes Calculate the No. of Slices from Water Temperature conditions

Calculate (No. of Tubes / Row) ,EQ.(3-1) & Calculate Length of H.EX, EQ.(3-2) & Calculate Depth of H.EX, EQ.(3-3)

YES ∆Twater = 0.1, nw=0.3,na=0.4,F=1

IF (Twi >Two)

NO ∆Twater= 0.1, nw=0.4 ,na=0.3,F=-1

Calculate the air mass flow rate for first row, Eq (4-6) Calculate the Heat Load for First Slice & First Row,EQ.(3-9)

Assume the Air Exit Temperature for First Row

B

A

B

A Calculate the Air Mass Flow Rate for First Row

J=1, Calculate the No. of Slices from Water Temp. Conditions

I=.1, Divided the Air Mass Flow Rate per No. of Slices

Calculate the Fluids Properties from Inlet Conditions

Calculate the Heat Load for First Slice

Calculate the Air Exit Temperature for First Slice

Correct the Fluids Properties

Calculate the Correct Air Exit Temperature for First Slice

Assume the Air Velocity Over Tubes.

Calculate the Height of Heat Exchanger

D

C

E

C

D

E

Calculate the Cross Sectional Area Temperature, Hydraulic Diameter, Reynolds No., Nusslat No. .

Calculate the Overall Heat Transfer Coffi.,

Calculate the Height of Heat Exchanger Calculate Error from Calculating Value & Assuming Value.

Repeat Process for all Slices (Ni) NO IF Error <=0.001 YES Calculate the Mean Air Exit Temp. Calculate the Error from Calculating Value & Assuming Value

NO

IF Error <=0.001 YES Repeat the process for all Rows (Nj)

Calculate Heat Exchanger Performance

Calculate Heat Exchanger Physical Characteristics

End

Flow Chart for Theoretical Gas Turbine Performance Prediction(CPTGTP)

STRAT

INPUT DESIGN INFORMATION Tmax,T1,Rc,Rt,f,H.V,ηc,ηt

Ømax=Tmax/T1,M=0 For Ø=Ømax to 1 step -1 T3=Ø*T1 Za=γa/(γa-1), Xcm=(Rc)^(1/Za) For Xc=Xcm to 1.1 step -.05 NDCW=(Xc-1)/(ηc*(Ø-1))

Zg=γg/(γg-1), Xtm=(Rt)^(1/Zg),n=cpa/cpg

NDTW=((1+f)*ηt*(1-1/Xt))/((1-1/Ø)*n) ηo=NDTW/(H .v*(1+ƒ))

PRINT T3,Ø,Xc,Xt,ηo,NDCW NDTW

NEXT Xc ,

NEXT Xt , NEXT Ø

END

Appendix (C) Gas Turbine parameter Groups

(C-1)

(C-2)

Related Documents