Study of flow through combustion swirler with the effect of diffuser on the recirculation zone

Similar documents
The Effect of Radial Swirl Generator on Reducing Emissions from Bio-Fuel Burner System

EPA Method 1 Field Procedures Example:

ME 239: Rocket Propulsion. Real Nozzles. J. M. Meyers, PhD

Pumps, Turbines, and Pipe Networks, part 2. Ch 11 Young

[4163] T.E. (Mechanical) TURBO MACHINES (2008 Pattern) (Common to Mech. S/W) (Sem. - II)

The effect of hydrogen containing fuel blends upon flashback in swirl burners. Nick Syred Cardiff School of Engineering Wales, U.K.

The Study of Interference Effect for Cascaded Diffuser Augmented Wind Turbines

Numerical Analysis on the CO-NO Formation Production near Burner Throat in Swirling Flow Combustion System

Design and Simulation of Very Low Head Axial Hydraulic Turbine with Variation of Swirl Velocity Criterion

FLAME AERODYNAMICS COMBUSTION AND FUELS

Detached Eddy Simulation of High Turbulent Swirling Reacting Flow in a Premixed Model Burner

Performance of an Open Ducted Type Very Low Head Cross- Flow Turbine

Changes in combustion properties of Natural gas when mixed with Hydrogen

Numerical Investigation of the Combustion of Methane Air Mixture in Gas Turbine Can-Type Combustion Chamber

On the topic of gas turbine combustion

Simulation of Engine Flow with Swirl Using Re-Stress Turbulence Model in KIVA Code

Comparison Between PIV Measurements and CFD Simulation on a Model of GT Annular Burner

THE FORMATION OF CO-NO IN A CONFINED CAN COMBUSTOR DUE TO SWIRLING FLOW

Experimentelle Untersuchung des Einflusses eines axialen Freistrahls auf drallstabilisierte, mager vorgemischte Wasserstofflammen

INTERNATIONAL JOURNAL OF ADVANCED RESEARCH IN ENGINEERING AND TECHNOLOGY (IJARET)

Interference of Wind Turbines with Different Yaw Angles of the Upstream Wind Turbine

EFFECTS OF LoD AND PoD IN COMBINED-HOLE FILM COOLING

Design and Testing of a Rotating Detonation Engine for Open-Loop Gas Turbine Integration

FLOW VISUALIZATION IN A SWIRLING FLUIDIZED BED Sulastri Sabudin 1*, Shahzada Khan Nawab Khan 1 and Mohd Faizal Mohideen Batcha 1

EXPERIMENTAL AERODYNAMICS FLOW TROUGH AN ORIFICE LAB SHEET

WATER JET CONTROL TECHNIQUE FOR SWIRLING FLOWS IN FRANCIS TURBINES DIFFUSER

COMSOL Multiphysics Simulation of Flow in a Radial Flow Fixed Bed Reactor (RFBR)

Performance Study of Solar Air Heater with Baffled Duct

THE AMERICAN SOCIETY OF MECHANICAL ENGINEERS 345 E. 47th St., New York, N.Y

Equipment for Engineering Education

Theory Comparison between Propane and Methane Combustion inside the Furnace

CFD on Small Flow Injection of Advanced Accumulator in APWR

CYLINDRICAL VESSEL LIMIT LOAD ESTIMATION FOR OBLIQUE NOZZLE BY USING ANSYS AS ANALYSIS TOOL

PERFORMANCE ANALYSIS OF NATURAL DRAFT WET COOLING TOWER AT OPTIMIZED INJECTION HEIGHT

Effect of temperature on atomization in gas centered coaxial injection systems

PAPER-I (Conventional)

Heat transfer modelling of slot jet impinging on an inclined plate

ISO 3966 INTERNATIONAL STANDARD. Measurement of fluid flow in closed conduits Velocity area method using Pitot static tubes

Fabrication and Study of the Parameters Affecting the Efficiency of a Bladeless Turbine

Design, Modeling & Analysis of a Submersible Pump and to improve the Pump Efficiency

Tutorial letter 101/0/2016

CFD RE-DESIGN OF A GAS TURBINE CAN-TYPE COMBUSTION CHAMBER HYDROGEN FIRED

Three-Dimensional Numerical Simulation of a Model Wind Turbine

Design of Radial Turbo-Expanders for Small Organic Rankine Cycle System

Analysis of numerical models for the formation of NOX for combustion of gases in the strong swirl flow

Effect of Twisted Tape Inserts On Heat Transfer in A Concentric Tube Heat Exchanger

Effects of non axisymmetric endwall contouring and film cooling on the passage flowfield in a linear turbine cascade

Investigation of Rotating Instability in the Last Stage of Low Pressure Turbine during Low volume flow Operation

Optimization Of Pattern Factor Of The Annular Gas Turbine Combustor For Better Turbine Life

Performance Analysis for Natural Draught Cooling Tower & Chimney through Numerical Simulation

Analysis of Flue Gas Flow Behavior in Economiser Duct Using Cfd

Aerodynamic Design of 2.5 MW Horizontal Wind Turbine Blade in Combination with CFD Analysis

MATHEMATICAL MODELING OF DRYING KINETICS OF CORN IN ELECTRON FIRED FLUIDIZED BED DRYER

Simulation of the Flow in a Packed-Bed with and without a Static Mixer by Using CFD Technique

Characteristics of fans used in low-power boilers

CRHT VII. Design and CFD analysis of Pico- hydro Turgo turbine. Paper no. CRHT17-11

Solar Flat Plate Thermal Collector

ASME PTC Flow Measurement

CFD Analysis of Pelton Runner

The Design and Analysis of First Stage Gas Turbine Blade with a Modification on Cooling Passages Using ANSYS

The Effect of Swirl Strength on the Flow Characteristics and Separation Efficiency of Supersonic Gas Separator

Recent Studies in Turbine Blade Cooling

An Investigation of Oxide Layer Impact on Heat Transfer in a Fuel Element of the MARIA Reactor

Numerical Prediction of Turbulent Combustion Flows for 1700 C Class Gas Turbine Combustor

The influence of the lengths of turbine blades on the power produced by miniature wind turbines that operate in non-uniform flow fields

CFD Analysis of a Low Head Propeller Turbine with Comparison to Experimental Data By: Artem Ivashchenko, Mechanical Solutions, Inc.

PEMP M.S. Ramaiah School of Advanced Studies, Bengaluru 1

Analysis of capillary diffuser for micro jet engine

HYDRAULIC DEMONSTRATION CHANNEL

MECHANICS OF PASTE FLOW IN RADIAL SCREEN EXTRUDERS

Design of a Centrifugal Compressor with Low Specific Speed for Automotive Fuel Cell

Aerothermal Study of Intermediate Turbine Ducts

Latest Simulation Technologies for Improving Reliability of Electric Power Systems

The University of Liège Wind Tunnel Facility

Numerical Investigation of the Flow Dynamics of a Supersonic Fluid Ejector

CFD analysis of coolant flow in the nuclear reactor VVER440

Fluid Mechanics. The Energy Equation [7] Dr. Mohammad N. Almasri

Advanced Gas Turbine Concept, Design and Evaluation Methodology. Preliminary Design of Highly Loaded Low Pressure Gas Turbine of Aircraft Engine

Performance Test and Structural Analysis of Cross-Flow Turbine

FSI ANALYSIS OF TRAILING EDGE REGION COOLING IN HP STAGE TURBINE BLADE

Experimental Testing of the Transverse Horizontal Axis Water Turbine

Design and Analysis of Multilayer High Pressure Vessels and Piping Tarun Mandalapu 1 Ravi Krishnamoorthy. S 2

Mathematical Modelling of Wind Turbine in a Wind Energy Conversion System: Power Coefficient Analysis

Effect of Fuels on Gas Turbine Can-Type Combustor using CFD Code

DECELERATED SWIRLING FLOW CONTROL IN THE DISCHARGE CONE OF FRANCIS TURBINES

Porous Wall Fed Liquid Fuel Nonpremixed Swirl-Type Tubular Flames

Numerical investigation of a GTM-140 turbojet engine

Effect of air flowrate on particle velocity profile in a circulating fluidized bed

Solid Circulation and Gas Bypassing in an Internally Circulating Fluidized Bed with an Orifice-Type Draft Tube

System Identification and Performance Improvement to a Micro Gas Turbine Applying Biogas

flow work, p. 173 energy rate balance, p. 174 nozzle, p. 177 diffuser, p. 177 turbine, p. 180 compressor, p. 184 (4.4b) p. 166

CFD ANALYSIS OF THERMAL INFLUENCE ZONE OF EARTH AIR TUNNEL HEAT EXCHANGER UNDER TRANSIENT CONDITIONS

FE ANALYSIS OF RUNNER BLADE FOR WELLS TURBINE

Analysis of the corrosion test process for heat exchangers in corrosion test chambers by CFD simulation

The Paradigm Case Method of Flowmeter Selection

COMPUTATIONAL AND EXPERIMENTAL STUDY OF OPTIMAL DESIGN ON HYBRID VERTICAL AXIS WIND TURBINE

Flow simulation and efficiency hill chart prediction for a Propeller turbine

Advanced Gas Turbine Concept, Design and Evaluation Methodology. Preliminary Design of Highly Loaded Low Pressure Gas Turbine of Aircraft Engine

ORCHID Turbine. Fluid-dynamic design and characterization of a mini-orc turbine for laboratory experiments

Hydro Electric Power (Hydel Power)

Transcription:

International Journal of Engineering Research and Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 3, Issue 6 (September 2012), PP. 68-73 Study of flow through combustion swirler with the effect of diffuser on the recirculation zone Ashok K.Dhakiya 1, Bhavin K.Shah 2, Arvind.S.Mohite 3 1,2,3 Department of mechanical engineering, Faculty of technology & engineering, M.S. University, Vadodara. Abstract Swirling flows offer an interesting field of study for aerospace and mechanical engineers in general and for combustion engineers in particular since it involves complex interaction of recirculation and turbulent mixing which aid flame stabilization in combustion systems. Swirling flow generates rotating flows, turbulence and free jet wakes at the downstream of swirler in combustion chamber. So there is complex interaction between pressure gradients and fluid flow. Swirling jets are used as a means of controlling flames in combustion chambers. The main objective of this study is to present the details of the experimental procedure and measurements made in an axisymmetric swirler model with and without incorporating diffuser after swirler. This study mainly focuses on the effect diffuser on the recirculation zone in the expansion chamber. The flow through 8 blade 10 inlet angle and45 outlet angle vane swirler is studied experimentally in two cases. A hemi spherical shape is provided on the hub in order to overcome the inlet flow field distortion. Swirler is placed coaxially in the inlet pipe of 120 mm inner diameter. Diffuser is placed after swirler in the first case. In the second case diffuser is not used, inlet pipe is directly connected to the expansion chamber. Expansion chamber and tail pipe is attached in the sequence. The recirculation velocities at the center core and outer wall regions are well captured with the five-hole probe of cobra type. Five-hole probe measurements along the entire length of the model at different stations, allow the determination of axial mean velocity component and pressure field parameters, which provide comprehensive information to aid and understand such complex recirculation flows. Flow field characteristics downstream of the swirler at different axial locations are examined in detail. The recirculation length and width are found to be strongly dependent on the amount of turbulence created by the vane swirler. Keywords Swirler, Adverse pressure gradient, Vane angle, Recirculation region, Five-hole probe. I. INTRODUCTION Swirling flow is a main flow produced by air swirler. such flow is the combination of swirling and vortex breakdown. swirling flows are widely employed in combustion systems for flame stabilization purpose. swirl flows offer an interesting field of study for aerospace and mechanical engineering in general and for combustion engineers in particular since it involves complex interaction of recirculation and turbulent mixing which aid flame stabilization in combustion system. swirling flows, which are highly complex, have the characteristic of both rotating motion and free turbulence phenomenon encountered as in jets and wake flows. the presence of swirl results in setting up of radial and axial pressure gradients, which in turn influence the flow fields. in case of strong swirl the adverse axial pressure gradient is sufficiently large to result in reverse flow along the axis and generating an internal circulation zone. generating an internal circulation zone. The degree of swirl is usually characterized by means of a non dimensional number called the swirl number S, which is defined as the ratio of the axial flux of swirl momentum (GØ) divided by the axial flux of axial momentum (Gx) times the equivalent nozzle radius. S = GØ / (Gx R) However for vane swirlers, the tan (θ) alone may be taken as a measure of the degree of swirl where θ is the vane angle. Based on swirl number S Swirling flows can be classified into weak swirl, medium swirl and strong swirl. If swirl number is less than 0.3 it is usually classified as weak swirl and if it is between 0.3 and 0.6 it is called medium swirl and if the swirl number is greater than 0.6, it is called strong swirl. Swirling flows are highly three dimensional and it is quite complex to obtain enough details experimentally to fully understand it and to comprehend the mechanisms involved. The main objectives of the paper are: (i) to provide a systematic and experimental measurements of axial, radial and tangential velocity profiles downstream of the swirler in order to understand the flow field characteristics, (ii) to provide a base line experimental data for the validation of numerical results, and (iii) to understand the effect of turbulence and recirculation zone created by vane swirler. II. BASIC GEOMETRY AND ITS PARTS The swirler model with vane angle 45 is made in transparent, Perspex material as shown in Fig. 1. The geometry consists of an inlet pipe of length 350 mm, 120 mm outer diameter and hub diameter of 40 mm. The inlet pipe is provided with a 10mm hole to measure the mean velocity at the inlet with the help of a static pitot tube. The inlet pipe is followed by a conical diffuser with inlet dia of 120 mm and 250 mm outlet dia. Diffuser is followed by sudden expansion round chamber 68

of 250 250 mm cross-section and a length of 1100 mm. The holes of 10 mm dia. are drilled in diffuser and chamber to measure the velocity by five hole probe. A tail pipe of 120 mm diameter and length of 1300 mm is provided to prevent the atmospheric disturbances in the development of the flow. The swirler is placed in the inlet pipe with the vane tip made to coincide exactly with the outer plane of the inlet pipe as stated by Mathur and Macallum7. III. DESIGN DETAILS OF 45 ANNULAR VANE SWIRLER The design of vane swirler is based on Mathur and Macullum1. The design of 45 vane swirler is shown in Fig. 1. There are eight vanes 2 mm in thick. The vanes are symmetrical, and the trailing edges of the vanes do not lie in the plane of the hub exit. The angle subtended by a vane at the axis, when viewed in the axial direction is 75, giving an overlap of 30 between adjacent vanes. The length of the hub is 175 mm and a hemispherical bluff body is attached upstream of the hub in order to smoothly guide the fluid particle circumferentially, impinging on the hub. IV. EXPERIMENTAL SET UP WITH TEST FACILITY Experiments have been carried out in a low speed vane blower with cascade tunnel test facility. Blower is used to provide the flow of air through combustion swirler set up. Blower has a backward flow type of vanes. Cascade wind tunnel is of 75*75*60 in the cross section and placed at the exit of the blower with the help of convergent type of nozzle and small pipe. In wind tunnel the condition of the flow is of stagnation type. The three-dimensional swirler geometry is made from transparent Perspex material and it consists of an inlet pipe of length 350 mm, outer diameter 120 mm and hub diameter of 40 mm. The inlet pipe is attached to the wind tunnel as shown in Fig. 3. The vane swirler is placed in the inlet pipe with vane tip made to coincide with the exit plane of the inlet pipe. fig.1 swirler portion geometry fig.2 geometry of combustion swirler with diffuser The inlet pipe is followed by a sudden expansion chamber of 250 250 mm round cross section of length 1100 mm. Holes of 10 mm diameter are drilled downstream of the swirler at various stations (A to P) in diffuser and chamber. The round chamber is followed by a tail end pipe of length 1300 mm and diameter of 120 mm. As already mentioned, the tail end pipe is provided to prevent the back flow of fluid affecting the swirl flow development in the round chamber. Based on the axial velocity profile at the inlet of the test section the bulk mean velocity is found to be 20 m/s for free jet. Table 1 gives the non-dimensional position of the measuring locations with respect to the axial length of the swirler downstream region, considering the exit of the swirler as the reference location to measure the positive X values as shown in fig. 3. 69

V. PROBE CALIBRATION The present experiment is related with the axial, tangential and radial velocity of the flow. So here five hole probe as shown in fig.4 is used to find out these velocity components. For this, the calibration of five hole probe is done and then it is used for the experiment. The five hole probe was calibrated in a separated fabricated calibration section of 150cm x 220cm attached at the end of blower shown in fig.5. It was calibrated for zero setting and for dynamic pressure and pitch angle constant. Fig.3 Experimental set up with five hole probe measurement procedure fig.4 details of five hole probe The probe was calibrated against a standard Pitot tube for dynamic pressure and pitch angle. Pressure signals from the Pitot tube and five tubes of the probe were read on multi tube manometer. Tunnel flow velocity was varied in steps and for each velocity readings of true dynamic pressure from standard pitot tube and the pressure of each tube of the probe after null position of side tube were recorded. These were repeated for pitch angle which was varied from -20 o to -20 o in steps of 5 o. The total pressure is then read directly from the centre tube when the pressure of side tubes is equal. VI. CALIBRATION FOR DYNAMIC PRESSURE AND PITCH ANGLE The different values of true dynamic head and probe head were plotted. After plotting all the values, average tread line is derived which give the value of dynamic head constant K = 2.32.At different pitch angle the value of (P 3 - P 1 ) / (P 5 - P 1-4 ) plotted and get average value of pitch angle constant which comes C = 8.72. VII. Experiment performed in two different cases (1) Case-I With diffuser TWO CASES OF EXPERIMENT Fig.5 Swirler with diffuser 70

(2) Case-II Without diffuser VIII. Fig.6 Swirler without diffuser MEAN VELOCITY MEASUREMENT The instrument used for measuring velocities at various stations is the five-hole probe with a spherical diameter of 6.5 mm. The details of the probe are shown in Fig.4. The sensing head is a cobra shaped to allow probe shaft rotation without altering the probe tip rotation. The spherical head is turned in its guide about its axis until the pressures on the apertures 4 and 5 are equal and the torsional angle β can be read off the dial on the fixed marker. The readings in the pressure transmitter s p1, p2, p3 and p4 are noted down. The calibration curve is used to get the remaining properties such as the angle of flow and speed of flow Ū, with all the components of velocity (Eqs (1)-(4)) with the pressure values measured from the above. The traversing mechanism is employed for linear traversing. Measurements of axial, radial and tangential components of velocities are obtained for every 10 mm radial distance at various X/D axial stations downstream of the swirler. A total of 16 axial locations downstream of the swirler axially are considered for the measurement. ----------------(1) U = Ū cos β cos α ---------------(2) V = Ū sin α ----------------------(3) W = Ū cos α sin β ---------------(4) Where U, V and W are the axial, radial and tangential components of velocity U, k24 is the constant obtained from calibration curve, p1, p2, p3 and p4 are the static pressures at respective holes and r is the density of air. IX. RESULT AND DISCUSSIONS The axial, radial and tangential velocities are measured at various axial locations. At each axial location the fivehole pitot sphere is moved for measuring the various velocity components. The fig.7 &8 and fig.9 & 10 shows the axial velocity profiles from experimental values for case-i and case-ii respectively. At station A in fig.7, it is found that the negative velocity exists only at the centre point which is mainly due to the presence of hub and the axial velocity increases gradually as we move from the centre to 20 mm for eight vanes swirler and 15 mm for six vanes swirler radially and thereafter the velocity reduces and almost reaches a constant value near the wall. The increase in the axial velocity is due to impart of kinetic energy to the moving fluid by the vane swirler. For eight vanes swirler, reverse velocity occurs from station A to station M and maximum reverse velocity occurs is of 4.92 m/s at station J in the central portion. While for six vanes swirler, reverse velocity occurs from station A to K and maximum reverse velocity occurs is of 4.68 m/s at station F in the central portion. It can also be seen from the profiles that the values of axial velocities are positive from station N to P for eight vanes swirler and from station L to P for six vanes swirler. Recirculation Zone for 45 o Vane The radial velocity, which is very small and negligible as compared to the axial and tangential aids in flame stabilization by providing a hot flow of recirculated combustion products and a reduced velocity region. The flame speed and flow velocity is matched in the recirculation zone. A good mixing of the reactants with air occurs due to the high turbulence generated by high shear stresses. 71

Fig.7 Recirculation zone for case-i Fig.8 Recirculation zone for case-i by using MATLAB The boundary of the recirculation eddy is determined by the radial points at which the forward mass flow equals the reverse mass flow at the axial stations. A torroidal vortex may be formed within the recirculation zone, which is sustained by the momentum exchange into the main stream. The region which is enclosed where the axial velocity reaches zero magnitude in each axial plane is connected to obtain the recirculation region. This recirculation zone encloses a large torroidal vortex reverse flow situated in the centre of the jet. It is found to produce highly stable flames and also enables matching of zones of high turbulence intensity with those of high fuel concentration, resulting in higher combustion efficiency. From fig.7 it can be seen that the central recirculation occurs up to the station M and ends at station N in eight blade swirlerwith diffuser. while for case-ii the central recirculation occurs up to the station H. From the recirculation zone for case-i axial velocity, the recirculation expands rapidly and width of the liner of the recirculation becomes maximum nearer about 45 mm at station J in the radial direction from the central axis. The total length of recirculation zone for case-i is 4.33D mm, starting from station A and ends at station N. maximum width of recirculation is 0.75Dmm at station J, 320mm downstream side of the swirler Fig 8 & 10 shows graph plotted by using MATLAB for both the cases gives more clear idea about recirculation zones. As shown in fig 9 the total length of recirculation zone for case-ii is 2.75Dmm, starting from station A and ends at station H. maximum width of recirculation zone is 0.66Dmm at station D, 100mm downstream side of the swirler Fig.9 Recirculation zone for case-ii Fig.10 Recirculation zone for case-ii by using MATLAB 72

X. PRESSURE LOSS FACTOR The pressure loss factor is a measure of the flow resistance imparted to the air stream when it passes through the vane swirler. This pressure loss should be as small as possible for better fuel consumption in a combustion system. It is evident that the introduction of the vane swirlers necessarily introduces a certain amount of pressure loss in the system. The swirlers of varying vane angle are optimized when the pressure loss factor of each is known. The minimization of the total pressure drop is necessary in order to have proper recirculation zone for flame stabilization and sufficient turbulence and mixing. Total pressure loss factor (PLF) is a dimension less number and is defined as Pressure loss factor (PLF) = Where the inlet stagnation pressure (P O1 ) is measured in the inlet pipe upstream of the swirler and exit stagnation pressure (P O2 ) is obtained at the exit of the expansion chamber and leaving to the atmosphere, ρ is density of the fluid in kg/m 3 and U O is mean bulk inlet velocity to swirler (m/s). XI. CONCLUSIONS For the case of eight blades swirler with diffuser we get positive velocities from 450 mm to the downstream of the swirler while for without diffuser we get positive velocities from 260 mm to the downstream of the swirler. Due to the tailpipe of 120 mm diameter attached at the end of the expansion chamber less low pressure region occurs in the expansion chamber. Tail pipe will avoid the atmospheric disturbance and detain the atmospheric entering into the chamber. Recirculation zone of case-i is batter than case-ii. A 45 8 vane swirler with diffuser provides a reasonably good recirculation zone as compared to swirler without diffuser which may help flame stabilization more efficiently. The pressure loss factor is less for the 8 vane swirler. The tailpipe provided at the downstream of the swirler will drop the pressure loss by avoiding the atmospheric disturbance. The measurement of axial, radial and tangential velocity helps to determine the mean flow field characteristics downstream of the swirler, which can be used for reacting flow studies in burners and combustors. REFERENCES [1]. Agrawal, A.K., Kapat, J.S., & Yang, T. 1996. Flow Interactions in the Combustor Diffuser System of Industrial Gas Turbines. ASME Paper 96-GT-454R [2]. Thundil Karuppa Raj & V Ganesan, Experimental study of recirculating flows induced by vane swirler, Indian Journal of Engineering & Material Sciences, Vol. 16, February 2009, pp. 14-22. [3]. R. Thundil Karuppa Raj, V. Ganesan, Study on the effect of various parameters on flow development behind vane swirlers, International Journal of Thermal Sciences 47 (2008) 1204 1225 1205. [4]. E.Kilik, Better Swirl Generation by using Curved Vanes, California State University, Long Beach, California.. [5]. M.L. Mathur, N.R.L. Maccallum, Swirling air jets issuing from vane swirlers. Part 1: free jets, Journal of the Institute of Fuel 214 (1967) 214 225. [6]. Mathur, M.L., & Maccallum, N.R.L. 1967. Swirling Air Jets Issuing from Vane Swirlers. Part 2: Enclosed Jets. Journal of the Institute of Fuel, June, 238 245. [7]. H.S.Choksi & Ambalal Patil, Design aspect & performance analysis of Gas Turbine Combustion Chamber,January-2002,Faculty of Technology & Engineering, M.S.University, Vadodara. [8]. Combustion Aerodynamics by J.M.Beer and Chigier. [9]. Combustion fundamentals by Fischier. [10]. T.H.Chang, Experimental Study on turbulent swirling flows in a cylindrical annuli using the PIV techniques, Division of Mechanical and Automation Engineering, Kyungnam University, Gyeongnam 631-701, Korea International Journal of Automotive Technology, Vol. 5, No. 1, pp. 17-22 (2004). 73