 Open Access
 Total Downloads : 147
 Authors : Mudhafar A. Ali, Salih E. Najim, Qais A. Rishack
 Paper ID : IJERTV6IS080092
 Volume & Issue : Volume 06, Issue 08 (August 2017)
 Published (First Online): 11082017
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
CFD Analysis of a 200MW Tangential Gas Fired Furnace Boiler
Mudhafar A .ALi Salih E.Najim
Dep. Mechanical Engineering Dep. Mechanical Engineering University of Basra University of Basra
Basra ,Iraq Basra ,Iraq
Qais A. Rishack
Dep. Mechanical Engineering University of Basra
Basra ,Iraq
Abstract Three dimensional numerical cases to study the combustion occur in furnace working by gas fuel for a 200MW boiler capacity. The idea of this study has an important role in industry and various engineering applications such as steam power station to produce electricity. Finite volume method can be adopted to perform by the gaseous fuel in furnace with tangential fuel firing system efficiently . Hence, it is planning that this research can be achieve by using ANSYS Fluent ( CFD ) software tool . The non – premixed combustion model and probability density function (PDF) table are used to analyze the processes ,species transport model is applied for prediction of pure methane and air mixture combustion . Design Modeler ( DM ) of Fluent is used for drawing and designing all the parts such as walls, celling , bottom, outlet and wind boxes in four corners with three levels . Initially this work was solved to get acceptable simulation results of temperature and velocity distribution inside the furnace. Finally, the theoretical results of the simulation proved an agreement with the design results of the manufacturer, moreover this study aims to further understanding of combustion process hoping to enhance designing and manufacturing of boilers.
Keywords Furnace, Non premix combustion, Temperature, Turbulent flow, CFD, Fluent.

IINTRODUCTION
Recently boilers and furnaces are important subjects had been studied for many industries applications to provide steam or hot water. There are two types of boiler, first fire tube boiler by which the fire passes the through the tubes ,second water tube boiler by which the fire surrounding the tubes that contain the water. The size of water – tube furnace is very large depending on the designer view, demand of power ,fuel type and burner design . [1]. Herein, there is interest is to study the aerodynamic of the flow of fuel and air entering the furnace in tangential directions firing , the combustion processes and effects on the temperature and velocity distribution, by using the finite volume analysis and FLUENT programme. Some of variables such as the mean velocity of fuel and air are estimated from the mass flow rates and also the temperatures of the furnace, provided by the contracting company for the power plant. The use of CFD to visualize parametric effects in furnaces which potentially can be an accurate and cost effective tool. The furnace boiler configuration is called ''Box type'' as shown in Fig (1)that has
been specified for production steam with high pressure and it is combustion by natural gas fuel. [2].
Azazi [3] presented study for power plant furnace with tangential corner firing utilize (2 – D) aerodynamics and thermal aspects by using FORTRAN computer program. In his study deduced that inside the fire ball and in the furnace the pressure amount with the design value is agreement. That losses were 14.5% in the furnace to get1500oC inside temperature and the tangential velocity played a great role for keeping the stability of the fire ball.
Fande and Joshi[4] Studied the effects of operating conditions in a tangential fired furnace using commercial code FLUENT. 3D combustor model was used to determine the temperature, velocity and other thermal characteristics like O2 mass fraction, CO2 mass fraction for a typical 660 MW utility boiler firing medium volatile coal. Simulations based on real cases at different operating conditions. Results appeared relations among the temperature, O2 mass fraction and CO2 has been clearly demonstrated based on the calculated distributions. Hence, for the above, they concluded, by their results showed good agreement with measured and reported data.
Fig 1 Furnace boiler configuration
Munisamy et al. [5] presented a CFD investigation in comparing between the( 30Â° )offdesign burner angle and (0Â°)design condition of burner firing angle with for 700 MW Tangentially fired furnace boiler using coal fuel , then can by boiler operation be reduce rear pass temperature unbalance by
tilting burner angle 30Â° ,hence their conclusion illustrate, this effect back to reduce the turbulence fire ball in mixing inside the furnace to reduce the pass temperature shifted the fire ball position. However reducing the burner tilt angle, lead to decreasing the efficiency of the furnace. Their results for tilting angle 30Â° for burner, hence, velocity profile explain that fire ball is shifted below than the plane so as the condition design for fire follows not clear as that of 0Â° burner firing angle.
Azeez et al. [6] used computational fluid Dynamic simulations on combustion process in 3D tangential furnace of 120 MW. Gas and air fired in to furnace from four corners at three different levels The effects of various parameters and configurations on the thermal and hydraulic behavior inside the boiler were examined. Temperature profiles and distributions , velocity, and combustion gas flow patterns and prediction process are analyzed and they concluded that temperature profile along the furnace height exhibited a lower temperature at the bottom furnace , a drastic increase at the combustion region, and remained approximately constant at the upper furnace.
Himachandra and Prasad [7] investigated CFD analysis on 210 MW, type of boiler tangential fire pulverize tower utilized coal fuel for combustion. At different positions temperature are estimated. In their results the tangential fired system generates a swirling fire ball at the center of the combustion zone which is shown in Fig10. The swirling flow is stronger at the lower level approximately13 m from ground level as compared at 30m from the ground level .The maximum temperature of 1850 K is found to exist at the central part of the furnace.

CASE STUDY SPECIFICATION
The present case study, is parallelogram furnace boiler which is suitable for capacity around, 200MW. Tangential corner firing is adopted with following dimensions: X=10.81m,Y=22.5m, Z =9.4860 m, with volume capacity 2307.23 cubic meter, schematic furnace chamber shown in Fig.( 2).The type of gas burners are flat head nozzle as shown in figure( 3a )are mounted with entrances of primary and secondary air. In wind box arrangement and angle of nozzle position in horizontal plane for corners 1,3 are 45o and for corners 2,4 are 36o, 54o respectively as shown in fig. (3b).

BOUNDARY CONDITION
The boundary conditions can be define by properties of entering fluid (Gas) such as temperature, velocity, mass flow rate, density and viscosity and air properties and viscosity and air properties for simulation is listed in Table 1 ,which has been done under design condition for different load operation [2]. The burners in present study with = 0o.

Boundary Conditions at the wall
The fluid velocity in direction, in the immediate vicinity of furnace wall, must be equal to that of wall itself, thus Z equal zero at all points on the surface. The same is for axial velocity U and radial velocity V both of these equal zero at wall.

Boundary Conditions at the inlet
In most problems the conditions of the entering fluid is known the temperature, velocity , viscosity and density may be given or calculated, this just when flow coming out from nozzles. Velcity is given from the equation of mass flow rate for thermal power station.
= (1)
Where :
fluid; is the mass flow rate of fluid gas for gas fuel, oil for liquid air ;mass flow rate for air (either primary air m p.a , or secondary air m s.a)
:density of fluid(gas ,oil, air)in (kg/m3),A: area of nozzle(gas, oil ,air ) in (m2).
= (V gas , V oil , V p.a , V s.a) = V nozzle , V nozzle =
, n: No. of nozzles
The analysis of velocity components illustrated in Fig. 4 Where ; : Angle of nozzle position in horizontal plan ,
: Tilt angle in present case = 0o .
(b) wind box corners
Fig. 2 Furnace case study
(a) Flat head nozzle
Fig. 3 Gas burner type and wind box corners arrangement[2].
V. THEORTICAL ANALYSIS
In the analysis of flow and combustion inside furnace, the following equations must be considered[8],[9] :
Continuity equation
+= 0 (2)
Equation of Motion (Momentum)[9,10] In z – direction:
( ) –
(3)
Where : Source=
In r direction :
( ) –


COMBUSTION AND CHEMICAL REACTIONS
The combustible of the fuel elements and compounds in the fuel with all the oxygen as combination that requires high temperatures enough to ignite the constituents, mixing or turbulence to provide intimate oxygenfuel contact, and sufficient time to complete process [10]. The equation of stoichiometric equation of combustion for fossil fuels are mainly compounds of carbon and hydrogen (hydrocarbons – CmHn).The reaction of its oxidation can be written by the equation of stoichiometry. It is important for complete combustion one mole of fuel (CmHn) necessary exactly
(m + mole of oxygen.
Cm H n+(m+ )O2 m CO+ H2O (7)
This is illustrated in the following as balanced equation[11]: CH4+2O2 CO2+2H2O (8)
1 mole of CH4 +2 mole of 1 mole CO2 +2H2O
The total mass of oxygen required per kg of fuel as shown in
Where; source=:
(4)
(r )
table II, given by following equation:
Item
Unit
Gas firing
Load(MCR)%
–
100
70
60
40
30
CH4
Mass flow rate ()
Nm3/h
4870
Kg/s
0.90
0.63
0.54
0.36
0.27
Temperature (T)
k
300
–
–
–
–
Density(
Kg/m3
0.67
–
–
–
–
LHV
J/kg
50Ã—106
–
–
–
–
Primary Air
Mass flow rate( )
Kg/s
4.55
3.18
2.74
1.83
1.36
Temperature (T)
K
300
–
–
–
–
Density(
Kg/m3
1.18
–
–
–
–
Secondary Air
Mass flow rate
( )
Kg/s
1.52
1.06
0.90
0.61
0.35
Temperature (T)
K
300
–
–
–
–
Density(
Kg/m3
1.18
–
–
–
–
TABLE I. Boundary Conditions[2]
(4)
In direction :
( ) –
Where; source=
Equation of Energy
(5)
Where ; Source =
(6)
Fuel Content
%Content by weight of the fuel
Combustion Equation
Theoretical O2 Required per kg of Fuel
1
C
0.855
C+O2 CO2
C%
2
H
0.115
2H2 +O2 2H2O
H%
3
S
0.025
S+O2 SO2
S%
TABLE 2.Oxygen total mass required per kg of fuel[2]
X = U = sin cos
Y = V= sin
Z = W = cos cos
Fig. 4 Analysis of velocity components
mox=C%+H2%+S% (9)
Most gaseous involved in combustion calculations can be approximated as ideal gases, the reacting substance can be modeled by equation of state [10].Analysis process neglecting ash value because it is not entering combustion. High heat value11,144 kcal/Nm, Low heat value10,085kcal/Nm3.

MATHMATICAL MODEL
The mathematical mod el implemented here is totally based on the commercial CFD code, FLUENT , where the gas flow is described by the of mass, momentum and energy equations. Firstly the geometry has been made by design modular for case study and basically made fluid domain implemented into model then edit mash and name the fuel primary air , secondary air injectors and out let. The standard k turbulence model, nonpremix model as species model
,probability density function (PDF) and the P1 radiation models are used .

k turbulence model
The turbulence model [11] is used to simulate the turbulence flow and combustion process. In present study standard kepsilon model and radiation model p1 and non premixed species model have been chosen, the turbulence viscosity is expressed as:
+ + =
(10)

Energy dissipation ()
+ + =
(11)
where (G)is referred to the generation term and is given by [12]and[13].
G=2[( + + +( )]+
+( (12)
=Col (13)
Where :Co is a constant ,when assuming high Reynolds number, the value to be proportional to 3/2/l) ,so the previous equation becomes ,t is a hypothetical property of the flow and must be modeled, and it varies with position, so the pervious equation become[14].
= (14)
Where C : is a constant .
Some adjustable constants are included in equations also consist of. By numerous iterations of data fitting for a wide range of turbulent flows, the values of these constants have been arrived .The values of coefficients illustrated in table


NUMERICAL SOLUTION
The present study deal with combustion process ,there is heat transfer between the gas and furnace walls, while the gas have been experiencing few transitions during entry from nozzle flowing to furnace region .Numerical analysis is perform by using commercial software ANSYIS FLUENT 15.The set of governing equations represent in present study
are convection equations .Theoretical analysis for turbulent combustion flow is done to solve the case study. The continuity, momentum and energy equations are numerically analyzed by using finite volume method, also the SIMPLE algorithm which is an acrnym for SemiImplicit method for Pressure Linked Equations, uses to get relationship between velocity and pressure correction to impose mass conservation and obtain the pressure field . Combustion regime choose in this process is nonpremixed where gas and air entered the reaction zones separately [15].

Design Modular
ANSYS version 15 was used to for drawing all parts of domain for case study with high accuracy using data documents given by manufacture company in manual [2] such design represented in (Fig.5 a).

The Mesh
The discrete representation of the geometry that is involved in the problem,in computational solutions of partial differential equations can be defining as meshing , needed for the following:(1)Designates the cells or elements on which the flow is solved,(2)To discrete representation of the geometry of the problem, (3)Has celss grouped into boundary zones where b.c.'s are applied , (4) Significant impact on rate of convergence ,solution accuracy and CPU time required . Finally the mesh generation for solve present case study shown in (Fig5.b).
TABLE 3 Turbulent Model constants
1.44
1.92
0.09
1
1.3
0.7
0.4
0.9
(a) Design model
(b)Mesh model
Fig..5 Design and mesh model


RESULTS AND DISSCUSION
The computational model which has been applied to the furnace of a 200MW boiler fired with operational conditions is selected in correspondence with the design data [2].Ycutting planes shown in table 4 within the furnace are chosen to analysis the results of temperature distribution and velocity inside furnace.
Y=2.5m
Y=10.265m
Y=11.782m
Y=22.5m
Y=4.5m
Y=8.748m

Temperature Distribution
The contours of static temperature distribution along Y axis direction shown in figure 6,the temperature increase rapidly with elevation especially in wind box region up to exit. Six planes (present study) located at the vertical center plane of the furnace in direction of y axis were chosen for analysis, these horizontal planes are located diagonally and in contact with the fuel and air nozzles, as shown in figure 6, thus the gas temperature distribution can be clearly analyzed and studied at full load operation condition. The fuel and oxidizer are separated before entering the reaction zone where they can be mixed and burned. The combustion reactions in such cases are called nonpremixed flames, where the fuel and air that are supplied through the nozzles in the wind box region combust as they enter into the furnace. This can be shown by the temperature gradually increases as it flows toward the center of furnace. The lower temperature region is located at the area below the combustion region near the furnace bottom high dynamic range (HDR) level then increase gradually .In gas recirculation level as effect of hot gases which introduced by pass from line of chimney, then gas temperature increases rapidly with the increase of elevation specially in the wind box region, but when reach the top near the location of primary super heater platen tubes(SH level), the temperature remains high in center of plane with a little decreasing near furnace walls.
Table 4 Locations of Y planes
Yplanes
Elevation (m)
Temperature
(K)
Bottom HDR level
2.5
1990
GAS RECIRC.FAN Level
4.5
2120
Wind box(level1)
8.748
2139
Wind box(level2)
10.265
2139
Wind box(level3)
11.782
2135
Super heater(SH)level
22.5
2060
Fig .6 Contours of static temperature for different horizontal planes
Figure 7 shows the temperature variation along Yelevation indifferent levels across the furnace. The real picture contours for static temperature shown in Fig. 8.
Fig..7 Temperature distribution along the elevation of combustion zone
Fig..8 Static temperature distribution along Yplane furnace

Velocity Distribution
The general flow where natural gas (methane) and air undergo combustion process inside the boiler furnace, the mix of fuel and air have been intensive in the wind box region ,so intense magnitude velocity appeared at all the corners of the furnace where the fuel and air injected while this intensity reduced slightly in levels under and upper the wind box region. The real picture for Yplane in Fig.9, illustrated the contours velocity magnitude for the whole domain.
Fig .9 Contours of velocity magnitude along Yplane furnace
Fig..10 Trajectories of velocity magnitude vector for gaseous along Yplane
As well as the trajectories of velocity magnitude vector for gaseous along elevated of the furnace shown by Fig.10.
Horizontal planes for levels that chosen in table 4 used to describe intensity of velocity magnitude by figures (1116).
Fig.11Intensity of Magnitude velocity at Y=2.5
Fig.12Intensity of Magnitude velocity at Y=4.5
Fig.13 Intensity of Magnitude velocity at Y=8.748m
Fig.14 Intensity of Magnitude velocity at Y=10.265m
Fig.15 Intensity of Magnitude velocity at Y=11.782m
Fig.16 Intensity of Magnitude velocity at Y=22.5m
Path lines are the trajectories that individual fluid particles follow. These can be thought of as "recording" the path of a fluid element in the flow over a certain period. Figures.(17 22),how swirling started from bottom of furnace up to exit,
Fig.17 Path lines at Y=2.5m
Fig.18 Path lines at Y=4.5m
Fig.19 Path lines at Y=22.5m
The swirling in different magnitude and direction have been occurred, include randomness movement near the furnace bottom
Fig.20 Path lines at Y=8.748.5m
Fig.21 Path lines Y=10.265 m
Fig.22 Path lines Y=11.782 m
HDR level gradually become highly as elevation increase especially inn region of wind box where combustion occurred due to movement of gaseous from four corners with high velocity and the picture show the high swirling occurred in center of 3 level for wind box ,up to exit the magnitude decreases and swirling vanish especially in the SH level .


CONCLUTIONS
The conclusions drawn from the present analysis are summarized below:

The CFD simulations of furnace boiler using fluent are found to be in close agreement with manufacturer data.

The temperature degree increases when furnace elevation increase especially at wind box region .

High magnitude velocity affected by intense mixing of air and fuel in four corners of wind box region.

High swirling magnitude effected by high temperatures .
REFERENCES

F. William Payne and Richard E. Thompson "Efficient Boiler Operations" Sourcebook, Fourth Edition, 1996 .

The Mitsubishi Heavy Industries Ltd . Manual for the AL – Hartha Power Station "Boiler and Auxiliaries "Vol.1 Pg. 12

F.A. Azazi Aerodynamics and Thermal Prdiction of Al Hartha Power Generation Plant Furnace, Ph.D. Thesis, College of Engineering , Basrah University, Iraq 2001.

A.Sh. Fande and P.V.Joshi "CFD Investigation of Combustion in 660 Tangential Boiler "International Journal of Engineering Research& Technology (IJERT)ISSN:22780181,Vol.3 Issue 1 ,January, 2014.

K.M.Munisamy, M.Z. Yusoff, K.S Thanggaraju., H. Hassan, A. Ahmed "Burner titlting angle effect on velocity profile in 700 MW Utility Boiler"7thInternationalConference on Cooling & Heating Technologies (ICCHT ),2014.

W. Azeez. , H. Harsil , Z.Yusoff " CFD Investigation Of Effect Of Tilt Angle To Flow And Combustion In 120MW Gas Fired Full Scale Industrial Boiler " International Journal of Scientific & technology Research,Vol.4,issue12,December,2015.

J. J. Himachandra and L. L. Prasad, "CFD Analysis of a 210 M210210MWTangential Fired Boiler '' International Journal of Scientific Engineering and Technology Volume No.6 Issue No.1pp:33 44,2017.

R. B . Bird , W . E . Stewart and Lightfoot W . E., Transport Phenomena, Wiley International Edition, USA, 1980.

H. Schlichting, Boundary Layer Theory , 6th edition, Mc Graw
Hill , New York, USA, 1968. Singer, Joseph G. (Ed.)," Combustion / Fossil Power Systems". Windsor, Conn., Combustion Engineering Inc.,1981.
[11] V.S. Arpaci, Micro scales of Turbulent Combustion, Progress in Energy and Combustion Science Vol.21 , P.P. 153 – 171 ,Pergamon,Press,1995. B.E.Launder, and Spalding D.B, Lectures in Mathematical Models of Turbulence, Academic press , London and New York, 1972

I. Kubo and F. Gouldin, Numerical Calculations of Flow Fluid Mechanics of Combustion , Joint Fluid Engineering and CSME Conference , Montreal , Quebec , May , 1974.

D.B. Taulbee,"An improved algebraic Reynolds stress model and corresponding nonlinear Stress model",Phys,Fluid , 28,pp 25552561,1992

T.J.Poinsot and D,P, Veynante "Combustion" Encyclopedia
\of Computational Mechanics. Edited by Erwin Stein, RenÂ´e
de Borst and Thomas J.R. Hughes John Wiley & Sons, Ltd. 2004.