Open access peer-reviewed chapter

Steady State Thermal Analysis of a Tri-Wing Solar Chimney

By R. S. Bello, C. N. Ezebuilo, T. A. Adegbulugbe and K. A. Eke

Submitted: October 22nd 2014Reviewed: October 24th 2014Published: May 6th 2015

DOI: 10.5772/59721

Downloaded: 1065

1. Introduction

The solar chimney is one of the technologies which work on the principle of buoyancy, in which air is heated through the greenhouse effect generated by solar radiation (heat energy) at low cost. The solar chimney is a passive solar ventilation system (non-mechanical) that can be installed on roofs or in walls. The heat is transferred based on the convective cooling principle which works on the fact that hot air rises upward; these chimneys reduce unwanted heat during the day by displacing interior (warm) air with exterior (cool) air. Solar chimneys are mainly made of a black, hollow thermal mass with an opening at the top as an exit for the hot air. The air in the room exits from the top of the chimney (Figure 1). The process can also be reversed for room heating.

Figure 1.

Solar chimney [1]

Solar chimney performance analysis for natural circulation dryers reported that simple air heater increases ventilation up to some extent but not sufficiently [1]. Also, Duffie, and Beckman, (2003) experimentally analyzed an inclined window-sized roof solar chimney and found its summer performance to increase. Duffie, and Beckman, (2003) equally studied the effects of various performance parameters like chimney width and height and solar radiation on flat-plate collectors. The thermal analysis on tri-wing collectors both at steady and transient states is an entirely new area of collector configuration research with limited research outcomes.

1.1. Project objective

For effective design of tri-wing solar chimney, solar parametric equations were utilized to model drying conditions at chimney inlet and outlet to maximize the differentials between system temperatures and air densities and their effects on drying applications. The air density depends on the temperature; hence, it also implies that the maximum differentials between the chimney air temperature and the ambient temp should give the best chimney performance.

The first mathematical model for the solar chimney design (Trombe wall) was given by [3, 16], and they also reported the concept of increasing airflow by increasing solar irradiation. This theoretical study also reported an air change per hour with change in the coefficient of fluid (air) discharge. Alter, (2011) reported the mathematical model of a conventional vertical chimney which operates under natural convention conditions where the temperature of the air inside the chimney is warmer than outside. Shiv et al., (2013) reported that solar chimney as a solar air heater may be represented by position (as vertical or horizontal chimney), or as a part of a wall (in the form of Trombe wall) or as a roof solar collector [6]. The roof solar chimney is the most convenient and mature technology used for buoyancy-driven natural ventilation systems [7, 8, 9].

A complete analysis of the tri-wing collector with a mathematical model is cumbersome because of its distinct features compared to an ordinary flat-plate model comparison of its performance effectiveness with experimental design data using high-precision apparatus and equipment offers a realistic solution. The objective of this study is to use analytical method to derive expressions for modelling drying effects of buoyant airflow created by solar heating of a tri-wing collector.

Advertisement

2. Experimental setup

The solar chimney used in modelling drying conditions is shown in Figure 2. The experimental solar chimney is a hollow cylindrical channel of glass glazing. The walls of the chimney were made as smooth as possible to reduce pressure losses due to wall friction. The absorbing surface (collector) is a tri-wing multi-flapped selective absorber plate draped inside the glazed glass. This chimney is mounted above the room space (drying chamber) through which dry air passes.

Figure 2.

Schematic diagram of solar chimney

3. Solar incident radiation on each wing

The incident angle of radiation reaching any collector surface is expressed by the general expression

Cos θ=sin δ cosφ cosγ+cos δ sin φ cosγcosω+cos δ sinγ  sinωE1

This expression was used to model the following equations for each of the wings in terms of ωas

Cos θ=(0.1761+0.053 cosω+0.81 sinω) for wing 1 E2
Cos θ=(0.353+0.1106 cosω) for wing 2 E3
Cos θ=(0.17640.053 cosω.17640.0ω) for wing 3E4

Previously the expression for all the wings has been derived as

Cosθz=0.9288Cosθ.9288E5

4. System heat transfer mechanism in the model

A complex differential analysis is required to fully analyze the multi-flapped collector due to shadows and the transient heat flow of flaps at different potentials (Figure 3). For instance, when the sun is between the azimuth of east and south, wings 1 and 2 receive full-area radiation at different incident angles, while wing 3 receives partial-area direct radiation due to the shadow of wing 2 cast on it. The case is reverse when the sun is between the azimuths of south and west, the critical time of changeover. When the sun is over the azimuth of 0°, wings 1 and 3 receive full-area direct radiation at different incident angles, and wing 2 receives partial direct radiation. This only happens intermittently (Figure 3).

Figure 3.

Azimuth positions of the solar insolation with respect to the three absorber surfaces

The analysis of a tri-axially configured surface is cumbersome because of its distinct features compared to an ordinary flat-plate model; to overcome this, an assumption is made that the collector absorber material is a highly conductive metal such that heat due to absorbed radiation is evenly distributed in all the flaps. This implies that solar radiation incident on the collector is such that there is no temperature gradient anywhere in the collector (i.e. the absorber is isothermal at any point in time).

4.1. Assumptions on analytical methods

In the course of analysis of this solar collector, a number of overriding assumptions were made which made the correlation of analyzed values comparable to experimental data values. Experimental data for this type of solar collector at the study location (Nsukka, Enugu state, Nigeria) was not available; however, such correlations provide a complete analytical method. Such assumptions considered a vertical tri-wing absorber collector heated by insolation to a temperature Tp with free convection boundary layer effects and also considered a zero air stream velocity which could be increased to some maximum value and then decreased to zero again due to free stream conditions, provided the gap between the plate and the cover is much greater than the boundary layer thickness.

To effectively model flow within the solar chimney, the analytical method employed required the following steps:

  1. First, determination of available solar radiation data in the study location.

  2. Second, evaluation of optical properties of glazing materials (assumed to be same as that of flat glazed collectors).

  3. Third, consideration of the electrical analogy of the system (from which the collector efficiency factor F’expression and the overall collector heat loss ULexpression were derived).

  4. Last, the heat transfer analysis and collector performance.

4.2. Solar insolation data at experimental setup location

The experimental setup was located at Nsukka, Nigeria, longitude 7° 23’ 45” E and latitude 6° 51’ 24” N [10], and the data acquisition method was adopted from the works of [11]. Solar data of the year and declination of that day were taken. This data was used as a case study (data presented in another chapter of this book) with the assumption that insolation recorded for a particular month is approximately the same for the same month for every other year. Also obtained from the literature were the measurement of instantaneous values of total insolation IT against diffuser Id and beam Ib insolations. From the instantaneous values, an expression relating to the total and beam insolations was deduced by the method of least squares for a third-order polynomial that fits the distribution of data. The polynomial is deduced as

Ib=0.01524IT30.0152IT2+0.988IT0.9882(MJm2)for IT>0.22(MJm2)E6

And with linear regression of natural logarithm of data distribution,

Ib=e(1.761lnIT1.761(KJm2)for IT>220(KJm2) E7
Id=ITIb E8

where IT=total insolation, Ib= direct beam insolation and Id=diffuse insolation.

These parameters were used to compute instantaneous values of hourly insolation for the day of the year when data were taken.

Solar radiation on a tilted collector is made up of the beam radiation component, the diffuse radiation component and the ground-reflected diffuse radiation components. Hence, the total hourly solar radiation on a tilted surface is the summation of the three components:

Ic=Ibc+Idc+IdgcE9

where subscript c denotes solar radiation on the collector. Furthermore,

Ic=IbcRb+Id[(1+cosβ)/2]+(Id+Ib)ρ[1/2(1cosβ)]E10

For the case study, considering the collector being a vertical wall and no vegetation covering the ground, β=90° and  ρ=0.2. Substituting these conditions into the equation of the total hourly radiation on the collector gives

Ic=IbcRb+0.5Id+0.1(Id+Ib)(Rb+0.1)Ib+0.6Id E11

Both sides of each wing are exposed to diffuse radiation, while only one side of the wings is exposed to direct beam radiation depending on the area factor receiving this direct radiation. Hence, the total radiation aimed at each wing of the collector is

Ici=IbRbi+2[0.5Id+0.1(Id+Ib)]=(Rb+0.2)Ib+1.2IdE12

Thus, the total radiation aimed at the whole collector in the chimney is

Ic=Ic1+Ic2+Ic3E13

Numbered subscripts 1, 2 and 3 denote wing positions.

5. Optical properties of glazing and absorption of solar radiation

The determination of collector performance requires the evaluation of the amount of solar radiation actually absorbed by the collector through glazing. Optical properties (transmittance, reflectance and absorptance) of the glaze cover, affects the absorption of solar radiation. These optical properties depend on the thickness of glazing and the refractive index and extinction coefficient of the glaze material. At a particular zenith angle of the sun, there are different incident angles of radiation on the circumference of the cylindrical glazing of the collector, unlike the flat-plate collector which has one uniform incident angle of radiation throughout the surface of glaze. Due to this feature, the evaluation of optical properties of a cylindrical glaze is complex. Though the glaze is cylindrical, the absorber plate is flat; hence, it was assumed that optical properties of a flat glass glazing are approximate to that of a cylindrical glass glazing. According to [12], expressions for reflection of unpolarized radiation passing from a medium with refractive index n1 to another with refractive index n2 are

rr=sin2(θ2θ1)sin2(θ+θ1)and r1=tan2(θ2θ1)tan2(θ2+θ1)E14

where  θ1and θ2are the angles of incidence and refraction. The average of reflectance of unpolarized radiation [12] is then

r=0.5(r+r)E15

where r. and r represents perpendicular and parallel reflectance components of the unpolarized radiation. The angle of refraction, θ2, is related to the refractive indices and the incident angle by Snell’s law expressed as

θ2= sini[n1n2(sinθ1)]E16

For air, n1 is approximately equal to 1. Consequently, the expressions for transmission of unpolarized radiation passing through a slab of glass from one interface to another interface and neglecting absorption decomposed to perpendicular and parallel components are as follows:

For a perpendicular component,

τ=[1τ1+τ]E17

For a parallel component,

τ=[1τ1+τ] E18

where τ and τ represent the perpendicular component and the parallel component of transmittance of unpolarized radiation for single cover glaze.

The average transmittance of these two components is

τr=(τ+τ)2 E19

In order to account for absorption, a loss factor was introduced by Bourguer’s law:

τ2= e[kL/cosθ2]E20

where k is the extinction coefficient (which varies from 4m-1 for white glass to 32m-1 for green glass). L is the thickness of the glass. In this study, the thickness of cover glass was taken to be 2.5mm, the extinction coefficient 5m-1 and the refractive index 1.526. The product k.L = 0.0025m x 5m-1 = 0.0125.

Some of the radiation passing through the glaze and striking the absorber plate is reflected back to the cover system. However, not all of this radiation is lost, since some of it is reflected back to the plate. The multiple reflections and absorptions between the plate and the cover is the greenhouse effect. To account for this greenhouse phenomenon, the actual fraction of incident radiation absorbed by the plate is called the transmittance-absorptance product, (τα). tThis is reasonably approximated as

(τα)=A(τ)αp E21

The constant A ranges from 1.01 to 1.02, but for conservativeness, 1.01 is preferred:

(τα)=1.01(τ)αp  E22

For transmittance of diffuse radiation, the effective incident angle θ1for vertical collectors is 59.5° for both ground and sky diffuse radiation. Hence, the effective refractive angle θ2of diffuse radiation from Snell’s law is

θ2= sini[n1n2(sinθ1)] E23

Likewise, the reflection components are

r= sin2(θ2 θ1)sin2(θ2+ θ1)E24
r= tan2(θ2 θ1)tan2(θ2+ θ1) E25

The transmittance components are

τ=(1τ)(1+τ) E26
τ=(1τ)(1+τ) E27

Cover absorptance τawas evaluated with the relation

τa= ewas evθ2E28

Hence, the resultant transmittance for diffuse radiation is

τd=τr.τa E29

Since all necessary optical properties of the glass glazing have been evaluated, the amount of solar radiation actually absorbed by the collector can be easily deduced. The total incident radiation aimed at each wing of the collector is

Ici=IbRbi+2[0.5Id+0.1(Id+Ib)]E30

the total absorbed solar radiation of each wing is

Si=(τα)bIbRbi+2(τα)b[0.5Id+0.1(Id+Ib)] E31
=(τα)bIbRbi+(τα)b[Id+0.2(Id+Ib)]E32

The mean absorbed solar radiation, S, of the whole collector plate is evaluated by the relation

S=1/3[S1+S2+S3]E33
Advertisement

6. The collector efficiency factor and the collector loss coefficient

Considering the unusual configuration of the tri-wing absorber plate, a section of it was reduced to a vertical flat-plate collector over both sides of the plate without back insulation. This implies that a wing of the absorber can be treated as a flat-plate air heater with flow over both sides of the plate. The determination of the collector efficiency factor F′ and the collector loss coefficient UL provided the basis for the analytical prediction of collector performance when used in the Hottel-Whillier-Bliss equation [12, 13]:

qu=F'[SUl(TrTa)] E34

The collector efficiency factor F'is defined as the ratio of the actual useful energy gain to the useful gain that would result if the collector absorbing surface had been at the local fluid temperature, that is

F'=qu[SUl(TrTa)] E35

where qu is the actual useful heat collector rate, S is absorbed radiation and Tris average fluid temperature.

The collector loss coefficient Ulis defined as the lumped overall heat loss value of the whole collector to the operating temperature difference of the collector, that is

Ul=Eloss(TpTa)E36

where Eloss= heat loss, Tp=absorber plate temperature and Ta= ambient temperature of surrounding air.

To conform to the performance equation of a flat-plate collector, the configuration of the chimney collector is transformed to suit the equation. Criteria for the transformation are as follows:

  1. The total area of the tri-wing absorber plate is equal to the total area of the flat absorber plate.

  2. The height of the chimney collector is the same as the height of the resulting flat absorber plate.

  3. The area of the circular inlet column of the chimney collector and the rectangular area inlet column of the resulting flat-plate collector are equal throughout the height of the collector.

  4. The resulting flat-plate absorber is positioned within the rectangular column such that the flow is halved.

  5. Size difference of the glass glazing of the two cases is neglected. The resulting flat-plate collector has an absorber plate 5.3m in height and 2.4m in breadth enclosed by a rectangular channel of glass glazing measuring 2.46m by 0.8587m in length and breadth. Thus, the flow width on either side of the plate is about 43cm.

6.1. Electrical analogy

According to [12], the usual procedure of deriving  F'and Ulis by converting a schematic diagram of the collector system to a thermal network of electrical analogy and analyzing the resulting circuit. A schematic diagram with equivalent thermal network is shown in Figure 4.

Figure 4.

Schematic diagram with equivalent thermal network

At the leading edge of the flat absorber plate, airflow is halved to flow on both sides of the plate. Absorbed solar energy heats up the plates to a temperature Tp, and energy is transferred to fluid at temperature Trthrough the convection heat transfer coefficient h2and to the cover glass through the linearized radiation heat transfer coefficient  h1. Energy is also transferred to the cover glass from the fluid though the convective heat transfer coefficient ht, and finally energy is lost to the ambient air through the combined convection and radiation coefficient Ut from the cover glass.

7. Energy balance equations

Energy balance equations on the cover, the plates and the fluid are derived as follows:

2Ut(TaTc)+2h1(TpTc)+2h2(TrTc)=0 E37
S+2h1(TcTp)+2h2(TrTp)=0 E38
2h1(TcTr)+2h2(TpTr)=qu E39

Simplifying and expanding the expressions as a function of heat transfer coefficients and temperatures yields

UtTaUtTc+hrTphrTc+h1Trh1Tc=0E40
0.5S+h1Tch1Tp+h2Trh2Tp=0 E41

Equations (40) and (41) can be expressed as

A=BTc=TpE42
D=ETp=Tc E43

where A=UtTa+h1Tr, B=Ut+hr+h1,C=hr, D= 0.5S+h2Tr and E=hr+h2.

Solving equations (42) and (43) simultaneously and re-substituting values into equations (40) and (41),

Tp=AC+BD[BEC2]=(UtTa+h1Tr)hr+(Ut+hr+h1)0.5S+h2Tr(Ut+hr+h1)(hr+h2)hr2 E44
Tp=AE+CD[BEC2]=(UtTa+h1Tr)(hr+h2)+(0.5S+h2Tr)hr(Ut+hr+h1)(hr+h2)hr2 E45

Subtracting Trfrom both sides of equations (44) and (45),

TpTr=Uthr(TaTr)+0.5S(Ut+hr+h1)(Ut+hr+h1)(hr+h2)hr2 E46
TcTr=(Uthr+Uth2) (TaTr)+0.5Shr(Ut+hr+h1)(hr+h2)hr2 E47

Substituting equations (46) and (47) into equation (39),

0.5qu=0.5S(h1hr+Uth2+hrh2+h1h2)+Ut(h1hr+hrh2+h1h2)+ (TaTr)(Ut+hr+h1)(hr+h2)hr2 E48

Rearranging equation (46),

qu=S(h1hr+Uth2+hrh2+h1h2)(Ut+hr+h1)(hr+h2)hr22Ut(h1hr+hrh2+h1h2)+ (TaTr)(Ut+hr+h1)(hr+h2)hr2 E49

Comparing equation (49) with the Hottel-Whillier-Bliss equation,  qu=F'[SUL(TrTa)], it can be deduced that

F'= h1hr+Uth2+hrh2+h1h2(Ut+hr+h1)(hr+h2)hr2 E50
UL= 2Ut(h1hr+hrh2+h1h2)+ (TaTr)(h1hr+Uthr+hrh2)+h1h2 E51

It has been demonstrated experimentally that h1=h2=haccording to Duffie and Beckman; therefore, equations (47) and (48) reduce to

F'= h(Ut+2hr+h)(Ut+hr+h1)(hr+h2)hr2 and UL= 2Ut(2hr+h)+ (TaTr)(Ut+2hr)+h E52

It is necessary to outline expressions of heat coefficients as functions of F'and  UL. These heat coefficients are expressed as

 Ut=(hwind+hrea) E53

where hreais the radiative loss coefficient between the cover and ambient air expressed as

hrea=εcα(Tc2+Ts2)(TcTs)E54

where glass cover emittance εc=0.88, Stefan-Boltzmann constant α=5.6 x 101w/m2k4and sky temperature Ts=0.0552Ta1.5.

The convective atmospheric heat loss coefficient due to wind, hwind, is dependent on the prevailing wind velocity given by [14] as

hwind=5.8+3.8E55

where (ms1)is the wind velocity (2  ms1) for the study location); thus,

hwind, location=5.8+3.82=13.4 w/m2°E56

The radiative heat loss coefficient between the plate and the cover after linearization is given as

hr=σ(εch+εphh)(Tc2+Tp2)(TcTp) E57

where absorber plate emittance εp=1αp=0.1.

8. Heat transfer and collector performance analysis

A review of the principles and theories governing natural free convection heat transfer [4] is done for the determination of other collector performance parameters. Thus, considering a vertical absorber plate of one wing of the collector heated by insolation to a temperature Tp, a free convection boundary layer is formed. The boundary layer is such that at the wall of the plate, the velocity of air stream is zero. This increases to some maximum value and then decreases to zero due to free stream conditions, provided the gap between the plate and the cover is much greater than the boundary layer thickness. The equation of motion for such a system is derived from the Navier-Stokes momentum equation of fluid flow given by [15]:

ρ(ux+vuy)=ρg+μ2uy2  E58

where ρis the density of air; xis the displacement in the vertical direction, the direction of buoyant flow; yis the displacement perpendicular to the direction of buoyant flow; Uis the velocity component in the x-direction; vis the velocity component in the y-direction; pis air pressure; μis the dynamic viscosity of air and g is acceleration due to gravity of value 9.18 (ms2).

Equation (58) can be expressed as

(Uux+vuy)=pxρg+μ2uy2 E59

where px=g. ρand ρis the free stream air density.

The density difference can be expressed in terms of the volume coefficient of expansion, β, defined by

=1V(VT)p=1V[VVTT]=[ρρTT] E60

Substituting equation (59) into (58) brings the momentum equation to

ρ(Uux+Vuy)=ρgβ(TT)+μ2uy2 E61

integrating equation (61) with respect to yfrom the boundary layer thickness δto the plate surface; since the velocity in the flow direction uis much greater than the velocity component perpendicular to the flow direction v, hence vtends to zero compared to  u. The integral momentum equation becomes

ddx[0δρu2dy]=0δρβg(TT)dyμuy=0 E62

The functional relation between the temperature distribution and the boundary layer thickness is needed to solve equation (61). From the prevailing conditions of the system, temperature T=Tpat  y=0, T=Tat y=δand dtdy=0 at y=δ,a parabolic temperature function can be assumed to represent the temperature profile of the system with respect to  y.

Therefore, let the function be of the form

T=Ay2+ By+E63

Substituting the conditions of flow into equation (61) gives

TTp=Ay2Anδy and TpT=2Aδ22δ2=Aδ2E64

Also the expression for the velocity profile is necessary to solve the integral momentum equation which can be deduced from the following flow conditions: u=0at y=0, u=0at y=δand 2uy2=βg(TT)/vat y=0.These four flow conditions can be fitted into a cubic equation of the form

uux=Ay3+By2+Cy+E65

where uxis a fictitious velocity which is a function of displacement in the  x-direction. Applying the conditions of flow into the cubic equation gives

D=0,C=[βg(TpT)δ4v],B=[βg(TpT)2v],A=[ββg(TpT)4vδ]}E66

Substituting these expressions into equation (65) and evaluating yields

uux=[βgδ2(TpT)4uxv]y(1yδ)2 E67

The terms involving  (TpT), δ2and uxmay be incorporated into the function uxso that the velocity profile can be reduced to

uux=y/δ(1yδ)2E68

Now the momentum equation of equation (61) when solved has its terms reduced to

0δ:u2dy=0δux2δ2[(y24y3δ+6y4δ24y5δ3+y6δ4]dy=ux2δ2(y33y4δ+65y6δ3+y7δ4)yoδ =ux2δ(13 =6523+17)E69
0δ:u2dy=ux2δ1.05 E70

Solving equation (69) in terms of temperature variables,

0δ:(TT)dy=0δ:(TpT)(12yδ+y3δ2)0δ=(TpT)δ(11+13)E71
0δ:(TT)dy=13(TpT)δ E72
δuδy|y=0=Uxδ(y2y2δ+y3δ2)|y=0=Uxδ(14yδ+3y2δ2)|y=0δuδy|y=0=Uxδ E73

Substituting equations (67), (68) and (70) into equation (61),

11.05(ddxux2)=13gβ(TpT)δvxδ E74

where v=μρis known as kinematic viscosity.

The energy equation for free convection of the collector system is expressed as

ρCp(UTx+VUTy)=k2Ty2, or (UTx+VUTy)=α2Ty2 E75

where α=kρCpis called the thermal diffusivity. The integral form of equation (72) with velocity tending to zero is

ddx0δ:U(TT)dy=αdTdy|y=0E76

From equation (67), it can be deduced that the relation of the order of magnitude between uxand δis  uxδ2. Inserting this relation into equation (68) and solving yields  δx1/2. Hence, it can be assumed that

ux=Ax1/2 and δ=Bx1/4E77

Introducing these expressions into equations (67) and (74) yields

5240A2Bx1/4=gβ(TT)(Bx1/43)(AB)vx1/4 E78

Solving for A and B in the above equation,

A=5.17v(0.952+Pr).95[gβ(TT) x3v2]x95[Pr]r95 E79
B=3.93Pr14[gβ(TT) x3v2]x3.[Pr]r3. E80

Substituting the expression for A into equation (77) gives the fictitious velocity uxat any point xas

ux=5.17v(0.952+Pr)12(Grx)r E81

where Grx is the Grashof number at point x, this is given by

Grx=gβ(TT) x3v2 E82

The resultant expression for boundary layer thickness from equation (75) is

δx=BxBe =3.93Pr12(0.952+Pr)14(Grx)14 E83

The heat transfer coefficient may be evaluated using equation (62)

qp=kAdTdy=kA(2δ)(TpT)=hA(TpT)E84

where it was derived that  h=2kδ. Hence, the dimensionless equation for the heat transfer coefficient becomes the expression for Nusselt number  Nux:

 Nux=0.508Pr12(0.952+Pr)14(Grx)14 E85

In the case study, it was assumed that solar heating of the collector is one of the constant heat flux conditions. In such a case, a modified Grashof number is introduced:

Grx*=Grx.Nux=gβ x4qpkv2E86

where  qpis the plate heat flux. The local heat transfer coefficient is correlated as

Nux=hxk=0.6(GrxPr)0.2,  105<Grx*< 1011,qp=constant E87

For the turbulent range, the average heat transfer coefficient for the constant heat flux cause is obtained from

h=1L01hxdx=1.25hxE88

It has been established that during distribution the local heat transfer coefficient is essentially with xwhen turbulent free convection is encountered. In such a case, the average heat transfer coefficient is h=hx. The mean velocity of air stream due to natural convection was obtained from evaluating the mean value of velocity distribution of equation (64) as

u=ux01[yδ(1yδ)]d(yδ)=0.0833UxE89

The volumetric flow rate is

Vo=CdAoUE90

The coefficient of discharge Cdis taken as 0.6 as adapted from [1]. Then the mass flow rate is obtained from the expression

m˙=ρVo E91

In the previous section, expressions for the overall collector heat loss coefficient UL and the collector efficiency factor F', which are very important performance parameters, were derived. These parameters are applied in the evaluation of useful energy rate extracted from the collector. To incorporate the flow rate and express the system energy equation using the collector inlet temperature, another performance parameter is introduced as the collector heat removal factor. The collector heat removal factor FRis defined as the ratio of actual useful heat collector rate to useful heat collector rate attainable with the entire collector surface at the inlet fluid temperature. This is stated mathematically as

FR=mCp(ToT1)Ac[SUL(T1Ta)]=[mCpAcU][1eAcUFmCp] E92

Thus, the energy equation of the system becomes

Qu=2AcFR[SUL(T1Ta)] E93

The efficiency of the whole collector system ηdefined as the ratio of the useful heat extracted from the collector Qu to the total incident solar radiation on the collector is mathematically expressed as

η=QuAcI=2FR[(τα)cULI(T1Ta)] E94

where (τα)is the effective transmittance-absorptance product of the collector glazing,

(τα)c=STI E95

where temperatures T1and  Taare the collector’s inlet and outlet temperatures, respectively.

9. Results and discussions

The modelled parameters from the analysis of the momentum equation provide a platform to analyze the thermal behaviour on each wing based on the use of a selective surface. This was used to graphically illustrate the heat flow pattern in the chimney and the radiative effects on the emitting and receiving surfaces. At any position of the sun, only two wings of the collector receive full-area direct radiation. At a smaller zenith angle, a larger portion of the farthest wing is lightened, while a smaller portion is lightened at a higher zenith angle. Figure 5 gives the direct beam area ratio ∀ for the three wings of the collector with respect to midpoint solar time.

Figure 5.

Direct beam area ratio ∀ for the three wings

The shading of the wings was taken into account to evaluate the actual direct beam radiation absorbed by each wing so that the average radiation absorbed by the whole collector can be deduced. To investigate the shape of the lightened area of a wing at any time interval, the zenith angle and the hour angle are required. The shape of the lightened region of the shadowed wing is always triangular depending on the zenith angle of the sun. Five variants of the shadowed wing model are shown in Figure 6.

Figure 6.

Variants of the shadowed wing model

Figure 7 shows the graph of temperature elevation and mass flow rate and insolation which points out that increase in temperature elevation is not always consequent on increase in the mass flow rate as illustrated by the negative gradient of mass flow rate between 11.30 am and 1.30 pm compared to the positive gradient of temperature elevation within the same insolation range.

Figure 7.

Graph of temperature elevation and mass flow rate and insolation

Advertisement

10. Conclusion

The result shows that modelling buoyant airflow within the chimney and the variation of air temperature elevation with insolation is minimal while the chimney has better efficiency at lower values of solar radiation. Thus, the method of analysis of thermal performance for this type of solar chimney has been accomplished. Nevertheless, a method of complete analytical evaluations that will give satisfactory results can be achieved by obtaining a table of correlation factor from experimental data values to analytical values of this kind.

© 2015 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

R. S. Bello, C. N. Ezebuilo, T. A. Adegbulugbe and K. A. Eke (May 6th 2015). Steady State Thermal Analysis of a Tri-Wing Solar Chimney, Solar Radiation Applications, Segun Raphael Bello, IntechOpen, DOI: 10.5772/59721. Available from:

chapter statistics

1065total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Performance Characteristics of Modelled Tri-Wing Solar Chimney and Adaptation to Wood Drying

By R. S. Bello, C. N. Ezebuilo, K. A. Eke and T. A. Adegbulugbe

Related Book

First chapter

Potential of the Solar Energy on Mars

By Dragos Ronald Rugescu and Radu Dan Rugescu

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us