The Comparison between the Effects of Using Two Plane Mirrors Concentrator and that without Mirror on the Flat- Plate Collector

Similar documents
Simplified Collector Performance Model

Scholars Research Library

Chapter 1 Solar Radiation

Construction and performance analysis of a three dimensional compound parabolic concentrator for a spherical absorber

Research Article Study on Effect of Number of Transparent Covers and Refractive Index on Performance of Solar Water Heater

Honeycomb Encapsulated Atmospheric Solar Collector Along with Single Basin Solar Still in Highly Energy Absorbing Weather Condition

ISESCO Science and Technology Vision

OPTIMIZATION of the GEOMETRY & MATERIAL of SOLAR WATER HEATERS.

Mohammed Awwad Al-Dabbas. Mutah University, Assistance Prof, Mechanical Engineering department, Mutah University, Karak,Jordan

UNIT FOUR SOLAR COLLECTORS

Thermal Analysis of a Flat-Plate Solar Collectors in Parallel and Series Connections Huseyin Gunerhan

Sunlight and its Properties II. EE 446/646 Y. Baghzouz

Page 1. Name:

Effect of Solar Angles on Incident Energy of the Flat Collectors

Experimental study on heat losses from external type receiver of a solar parabolic dish collector

Calculating equation coefficients

The optimal tilt angle of a solar collector

INTERNATIONAL JOURNAL OF MECHANICAL ENGINEERING AND TECHNOLOGY (IJMET)

Estimation of Hourly Solar Radiation on Horizontal and Inclined Surfaces in Western Himalayas

CFD ANALYSIS OF HEAT TRANSFER AND FLUID FLOW IN FLAT PLATE NATURAL CONVECTION SOLAR AIR HEATER

THERMAL PERFORMANCE OPTIMIZATION OF A FLAT PLATE SOLAR WATER HEATER COLLECTOR USING MATLAB

Thermal Analysis of Solar Collectors

Feasibility study of one axis three positions tracking solar PV with low concentration ratio reflector

ME 476 Solar Energy UNIT THREE SOLAR RADIATION

An Experimental Study of Flat Radial Fresnel Len Solar Collector for High Temperature Thermal System

CHAPTER 3 PROBLEM DEFINITION AND OBJECTIVE

4. Solar radiation on tilted surfaces

FLATE Hillsborough Community College - Brandon (813)

Chapter 2 Available Solar Radiation

Principles of Energy Conversion Part 11A. Solar Energy Insolation

EFFECT OF SOME PARAMETERS ON LINEAR FRESNEL SOLAR CONCENTRATING COLLECTORS

Experimental Study of the Performance of Compound Parabolic Concentrating Solar Collector

Basic information about developed Calculator for solar / wind hybrid power supply

ME 430 Fundamentals of Solar Energy Conversion for heating and Cooling Applications

Simulation Application for Optimization of Solar Collector Array

Measuring the Solar Constant

EXPERIMENTAL INVESTIGATION OF DIFFERENT TRACKING MODES OF THE PARABOLIC TROUGH COLLECTOR

Evaluation the impact tilt angle on the sun collectors

Thermal conversion of solar radiation. c =

Solar Radiation 1.1 THE SUN

Motion of the Sun. View Comments

Theoretical Design and Analysis of Gravity Assisted Heat Pipes

Theoretical Analysis of Overall Heat Loss Coefficient in a Flat Plate Solar Collector with an In-Built Energy Storage Using a Phase Change Material

Double-Skin Facade in Low-Latitude: Study on the Absorptance, Reflectance, and Transmittance of Direct Solar Radiation

EE Properties of Sunlight. Y. Baghzouz Professor of Electrical Engineering

AVAILABLE SOLAR RADIATION THEORETICAL BACKGROUND

Development and Validation of Flat-Plate Collector Testing Procedures

Determination of Optimum Fixed and Adjustable Tilt Angles for Solar Collectors by Using Typical Meteorological Year data for Turkey

LECTURE 3 - SOLAR ENERGY AND SOLAR RADIATION -

Solar Energy Fundamentals COURSE CONTENT

Chapter Seven. Solar Energy

Heat Transfer F11-ENG Lab #5 Photovoltaic Solar Cell School of Engineering, UC Merced.

COMPUTER PROGRAM FOR THE ANGLES DESCRIBING THE SUN S APPARENT MOVEMENT IN THE SKY

Design of Concentrating Parabolic Reflector System to Generate Energy

MEASUREMENTS OF TIME-SPACE DISTRIBUTION OF CONVECTIVE HEAT TRANSFER TO AIR USING A THIN CONDUCTIVE-FILM

NUMERICAL SIMULATION OF THE AIR FLOW AROUND THE ARRAYS OF SOLAR COLLECTORS

Optimum Collector Tilt Angles For Low Latitudes

Estimation of insolation of flat plate solar collectors in selected cities in Nigeria

Observer-Sun Angles. ), Solar altitude angle (α s. ) and solar azimuth angle (γ s )). θ z. = 90 o α s

1. SOLAR GEOMETRY, EXTRATERRESTRIAL IRRADIANCE & INCIDENCE ANGLES

Department of Energy Science & Engineering, IIT Bombay, Mumbai, India. *Corresponding author: Tel: ,

Monthly performance of passive and active solar stills for different Indian climatic conditions

Experimental investigation of the performance of a Parabolic Trough Collector (PTC) installed in Cyprus

Selection and evaluation of different tracking modes performance for parabolic trough solar collector

Solar Energy Cooking with the Sun

Flow and Temperature Analysis inside Flat Plate Air Heating Solar Collectors

Optical analysis of low concentration evacuated tube solar collector

Performance Assessment of PV/T Air Collector by Using CFD

Performance Investigation of Cavity Absorber for Parabolic Dish Solar Concentrator

Modeling of the Optimum Tilt of a Solar Collector to Receive Maximum Radiation

Simulation of a linear Fresnel solar collector concentrator

EFFECT OF NON-UNIFORM TEMPERATURE DISTRIBUTION ON SURFACE ABSORPTION RECEIVER IN PARABOLIC DISH SOLAR CONCENTRATOR

Effect of Selective Coatings on Solar Absorber for Parabolic Dish Collector

Theoretical Analysis of Solar Thermal Collector with a Flat Plate Bottom Booster Reflector

StudyOn The Thermal Performance Of Flat Plate Solar Air Collectors With Dust Deposition On The Transparent Covers

Numerical study of optical performance of a parabolic-trough concentrating solar power system

Solar Flat Plate Thermal Collector

ISO/TC 180/SC. Date: ISO/xxxx yyyy. Secretariat: SA. Solar energy Collector fields Check of performance. Énergie solaire xxx yyy

SURFACE ORIENTATIONS AND ENERGY POLICY FOR SOLAR MODULE APPLICATIONS IN DHAKA, BANGLADESH

Questions you should be able to answer after reading the material

Development and Validation of Flat-Plate Collector Testing Procedures

Optimizing the Photovoltaic Solar Energy Capture on Sunny and Cloudy Days Using a Solar Tracking System

Performance Evaluation of Low Inertia Multi- Stage Solar Still

ENHANCEMENT OF THE HEAT TRANSFER RATE IN FREE CONVECTION SOLAR AIR HEATER USING PIN SHAPED ARTIFICIAL ROUGHNESS ON ABSORBER PLATE

EXPERIMENTAL STUDY ON A CASCADED PCM STORAGE RECEIVER FOR PARABOLIC DISH COLLECTOR

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 6, Issue 11, May 2017

Chapter 2. Heating Earth's Surface & Atmosphere

The Sun. Fabio Peron Università IUAV - Venezia. Earth-Sun relationships. The Sun. Photosphere (Emits much of the solar radiant power)

16. More About Polarization

Energy Efficiency, Acoustics & Daylighting in building Prof. B. Bhattacharjee Department of Civil Engineering Indian Institute of Technology, Delhi

XI. DIFFUSE GLOBAL CORRELATIONS: SEASONAL VARIATIONS

Thermal Analysis of a Compound Parabolic Collector

EGEE 437: HWK #2. Brownson. Yaqdaan Alkayyoomi, Mackenzie Ailes, Nicholas Minutillo, Sheel Vora. Group 17. Due on Thursday, Feb.

Optimization of the Air Gap Spacing In a Solar Water Heater with Double Glass Cover

EAS 535 Laboratory Exercise Solar Radiation

CONSTRUCTION AND CALIBRATION OF A LOCAL PYRANOMETER AND ITS USE IN THE MEASUREMENT OF INTENSITY OF SOLAR RADIATION

BMM4753 RENEWABLE ENERGY RESOURCES

Lecture Notes Prepared by Mike Foster Spring 2007

HEAT TRANSFER ENHANCEMENT OF SOLAR FLAT PLATE COLLECTOR BY USING V CORRUGATED FINS AND VARIOUS PARAMETERS

Basic Solar Geometry. Contents

Transcription:

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 The Comparison beteen the ffects of Using To Plane Mirrors Concentrator and that ithout Mirror on the Flat- Plate Collector S. Himran, B.Sudia,.H. Piarah Universitas Hasanuddin, Indonesia Kampus Tamalanrea, Makassar 90245, Indonesia Phone: +062 0411 588400, Fax: +062 0411 588400,.mail: syukri_h@yahoo.com ABSTRACT: The aim of the study is to compare beteen the effects of using to planes mirror concentrator and that ithout mirror on the flat plate collector. The ork is performed on solar ater heating system ith flat plate collector installed in it. Measuring of the data are taken beteen 9.00 to 16.00 h over a day. Based on the global and diffuse radiations ith the average of 594 and 119 /m 2, the average efficiencies of the system over a day are 41.90 % and 22.80 % for collector ith and ithout mirror respectively Keyords: plate collector, mirrors concentrator NOMNCLATUR A p - absorbing area [m 2 ] I b - hourly beam radiation [/m 2 ] I d - hourly diffuse radiation [/m 2 ] I g - hourly global radiation [/m 2 ] I bm - reflected beam radiation [/m 2 ] I T - flux falling on the surface of collector [/m 2 ] r - tilt factor [-] S - incident solar flux absorbed in the absorber plat[/m 2 ] T - temperature [ o C] U - loss coefficient [/m 2 -K] Greek letters - solar altitude angle [ o ], absorptivity [-] - slope angle of collector [ o ] - declination angle [ o ] - latitude angle [ o ] - incidence angle of reflected radiation on collector surface [ o ] - incidence angle on collector surface [ o ] () - transmissivity-absorptivity product [-] - hour angle [ o ] ξ - azimuth angle [ o ] - mirror inclination angle [ o ] 080-1

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 Subscripts a - absorption b - beam d - diffuse g - global l - overall coefficient p - absorber plate r - reflection and refraction u - useful - all, - ast M - mirror - est 1. INTRODUCTION From many studies, it has been proved that a flat-plate collector for heating liquids and gasses gives only to a maximum temperature around 100 o C. hen higher temperatures are required, it becomes necessary to concentrate the incoming solar radiation. Concentrating collectors are of various types and can be classified in many ays i.e.: reflecting type utilizing mirrors and of the reflecting type utilizing Fresnel lenses. The reflecting surfaces used may be flat, parabolic or spherical. In this paper e analyzed the performance a flat-plate collector ith adjustable mirrors at the edges of collector to reflect radiation to the absorber plate. An advantage associated ith this type of concentrating collector is that the diffuse component of the incoming solar radiation is not entirely asted. A flat-plate collector ith plane reflectors is a simple non-imaging concentrating collector and represents an effective means of getting slightly higher temperatures than are obtainable ith a flat-plate collector alone. 2. XPRIMNTAL APPARATUS AND PROCDURS: The major components that make the flat-plate collector are: (1) absorber plate made of Alsheet metal painted black, 0.4 mm thickness, size 1.60.9 m, (2) tubes fixed to absorber plate ith nine tubes ½ inch diameter ith spacing 9.9 cm, made of brass. (3) to transparent cover glasses 4 mm thick ith spacing 4 cm. The spacing beteen absorber plate and first cover glass is 4 cm, (4) collector box made of plyood. The bottom and sides of the box are insulated by foam cork 4 cm thick. A schematic diagram of the test set-up is shon in Fig.1. It is a closed lop consisting of the solar ater heater is to be tested. ater can be circulated by loering the outlet pipe from collector 10 cm belo the ater surface in the tank. The experiments ere performed for collector ith mirror and that ithout mirror, during 8 hours over a day, from 9.00 AM to 16.00 PM. The solar radiations in to days give nearly the same value. The collector is equipped by to mirror-reflectors (mirrors), one of hich facing est and the other ast. 080-2

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 2 4 3 N 1 S 1. ater flat-plate collector 2. ater tank 3. circulating pipes 4. concentrating mirrors Figure 1 Schematic diagram of a closed loop set-up for testing flat-plate collector ith concentrating mirror 3. ANALYSIS: From Figure 2 and 3, the angle of incidence of direct solar radiation on the collector can be related by a general equation by 3 position angles of the sun in the sky i.e.: declination angle, latitude angle, hour angle ω and 2 orientation angles of the plane of the collector: slope of the collector, and surface azimuth angle ξ. The governing equation of incidence angle on the surface is given by: cos sin sin cos cos cos cos sin cos cos cos cos sin cos sin (1) cos sin sin sin The flat-plate collector is facing due south and is inclined at an angle = 0 o. Since Makassar located at latitude 5.1 o S and longitude 119.42 o. The angle of incidence on the collector, based on equation (1) becomes: cos sin sin cos cos cos (2) The inclination angle of the reflector is adjusted manually every hour. In order to get the optimum solar radiation reflected in to collector, the reflector mirrors are adjusted in such a ay that the sun s rays striking the top edge of the mirrors are reflected to the top edge of the collector. This needs the reflecting mirrors have the same length L as the collector in order to boost the output, Fig. 4. The inclination angle can be calculated from equation (1) by changing to, the incident angle on top edge of the reflector by changing to r, see Fig. 5. The incident angle on top edge of the reflector is given by: 080-3

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 cos r sin sin cos cos cos cos sin cos cos cos cos sin cos cos sin sin sin sin (3) Figure 2 The celestial sphere = solar altitude; = declination angle; = latitude angle; = hour angle; = solar azimuth angle Latitude for Makassar: = 5.1 o Declination: o 360 23.45 sin 365 284 N Day of the year on 24 October 2009, N = 298 days, gives declination angle = 12.786 o Surface azimuth angle: = 0 o for facing south; = 180 o for facing north Hour angle: o = (12 t) 15 o, t = local apparent time 080-4

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 Normal to tilted plane z N Normal to horizontal plane (zenith) Tilted plane Horizontal plane S Figure 3 Diagram illustrating the angle of incidence z =zenith angle, =solar altitude angle, = slope of the plane, = solar azimuth angle, = all azimut angle L Mirror reflector Mirror reflector Flat-plate collector L 1 2 Figure 4 Flate-plate collector ith reflectors Based on equation (3) and Figure 4, the incident angles of beam radiation on reflector θ r are given by : 080-5

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 for reflector facing est, surface azimuth angle: = 90 o cos r sin sin cos cos cos cos cos for reflector facing ast, surface azimuth angle: = 90 o cos r sin sin cos cos cos cos cos From Fig: 5, e found the relations for,, r are: cos sin sin..(4.a) cos sin sin (4.b) = 90 o r = 2 Then: = 2 r (5) From equations (4.a), (4.b) and (5), e found by trial-and-error, the inclination angle of the reflector and for facing est and ast respectively. r r Reflector= mirror r Flat-plate collector Flat-plate collector r Reflector =mirror (a) N (b) S Figure 5 Flat-plate collector and reflectors facing due est and ast, = 0 o The flux incident on the top cover of collector I T and the flux absorber in the absorber plate S are: Solar radiation flux incident on the top cover of flat-plate collector I T and flux absorbed in the absorber plate S are sum of the beam and diffuse radiation falling directly on the surface and the radiation reflected onto the surface from the surrounding ith the factors r b, r d, and r r. Fluxes I T and S are given by: I T = I b r b + I d r d + (I b + I d ) r r (6) S = I b r b () b + {I d r d + (I b + I d ) r r }() b (7) Tilt factor r b, diffuse radiation factor r d and reflected radiation factor r r are given as: r cos, r 1 cos 2, and 1 cos 2 respectively. b z d r r here, z,, are incidence angle on the surface, zenith angle, surface slope angle, reflectivity respectively. Since = 0 o, then r b and r d = 1, and r r = 0. 080-6

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 Based on equations (6) and (7), solar radiation incidents on the existing flat-plat collector become: (i). Incident flux on the top of the collector: ithout mirror: I I r I r (6.a) T b b d d ith mirror: I I I r I r (6.b) Tm b bm b d d I bm = beam radiation reflected by the to-mirrors = I b + I b I b I b b sin, I sin b cos I cos (ii). Flux S absorbed on the absorber plate: ithout mirror: S I b r b b I d r d d (7.a) ith mirror: Since the reflected radiations are not influenced by r b and r d then, S I I r I r (7.b) M b bm b b d d d The factors () b and () d are transmisivity-absorptivity product for beam and diffuse radiations. The procedures to calculate () can be found in Ref. [2]. The useful heat gain = the rate of heat transfer to the orking fluid: q u =A p S q l (8) The instantaneous collection efficiency is given by: u (9) Radiation Useful incident heat gain on the collector 4. RSULTS AND DISCUSSIONS A q p I T The global (I g ) and diffuse (I d ) solar radiations ere measured by using digital pyranometer over a day from 9.00 to 16.00 h and the variations are as shon on Fig.6. For performing the calculations, data are taken for the temperatures at: plate-absorber surface, to glass cover surfaces, and ater at the inlet and outlet ports of the collector. The ambient temperature and the hourly ind speed over the surface of the collector are also measured. Based on the equations (1) to (5), e can find the variation of hour angle ω, incident angle on the collector surface θ, incident angle on mirror facing est θ r and facing east θ r, inclinations of the mirrors facing est and facing east. The results obtained are given in Table 1. By considering the calculated inclination angles and for every hour during the day, the slopes of the to mirrors are adjusted. All the data calculated are available in the files. Figure 6 shos the variations of global I g and diffuse I d radiations during a day. It is seen that the values increase a little bit sharp from 0900 to 1100 h, touch a peak around noon and then drop sharply to 1600. It is also interest to study the performance of collector h for collector ith and ithout mirrors over a day. The variation of the important items are tabulated and presented on graphs are shon in Table 2 and Figures 7, 8, and 9. The variation obtained indicates the strong independence of these items on the radiation incident on the collector. From Figure 7, it seen 080-7

Ig, Id, (/m 2 ) ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 that the flux incident I T for collector increase from 0900 to 1100 h and has a maximum value at noon then decrease to 1600 h. It is also shon that collector using mirrors give higher value compared to that ithout mirrors. The flux incident I T using mirror facing est has higher values at the 0900 h and the decrease to 1600 h and inversely for collector ith mirror facing ast gives loer values at 0900 h and higher at afternoon 1600 h. The overall heat loss coefficients do not vary much during the day but have a higher value for collector ith mirrors compared to that ithout mirrors. The average efficiencies over a hole period, during hich considering the useful energy q u, flux absorbed S, heat loss from collector q l, flux incident on the top cover I T, could also be calculated. Making the approximation, then the efficiencies average over 8 hours from 0900 to 1600 h for collector ith and ithout mirrors ork out to be 41.9 % and 22.8 % respectively. 900 800 700 600 500 400 300 200 100 0 Ig Id 9 10 11 12 13 14 15 16 Time (h) Figure 6 Variations of global I g and diffuse I d radiations over a day Table 1 The variation of the inclination angles over a day on 24 October. N = 297, δ = 12.79 o, ξ. = 90 o, ξ. = +90 o Time 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 ω( o ) 45 30 15 0-15 -30-45 -60 ( o ) 45.74 31.27 17.57 9.25 17.57 31.27 45.74 60.43 r ( o ) 74.829 69.904 64.978 60.054 55.133 50.220 45.320 40.447 r ( o ) 45.475 50.328 55.208 60.10 65.007 69.918 74.835 79.758 ( o ) 30.342 40.193 50.044 59.893 69.734 79.561 89.360 99.106 ( o ) 89.050 79.345 69.585 59.801 49.987 40.164 30.330 20.485 080-8

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 Table 2 Performance of a flat-plate collector ith and ithout mirrors over a day, T a 32 o C LST(h) 09.00 10.00 11.00 12.00 13.00 14.00 15.00 16.00 I b (/m 2 ) 382.2 403.4 509.8 613 544.2 496.4 427.4 388.6 I d (/m 2 ) 90.8 108.6 127.2 159.0 137.8 123.6 108.6 97.4 I bm (/m 2 ) 543.58 461.26 523.86 608.86 599.21 567.60 595.50 757.07 I T (/m 2 ) 473.00 512.00 637.00 772 692 620.00 536.00 486.00 I TM (/m 2 ) 891.74 930.21 1010.10 1130.11 1100.21 1066.46 1060.12 1054.60 S (/m 2 ) 300.0 345.0 380.0 430.0 404.54 376.0 345.0 294.0 S M (/m 2 ) 563.05 644.44 699.58 790.0 751.03 720.0 685.0 615.0 T p ( o C) 65.2 70.6 76.1 80 77 74 71.5 65 T pm ( o C) 68.0 73.5 83.0 95.0 87.0 83.0 74.5 66.0 U l (/m 2 -K) 4.96 5.55 5.62 5.67 5.45 5.34 5.33 5.04 U l. M (/m 2 -K) 5.60 5.77 5.85 5.91 5.87 5.76 5.73 5.38 q u () 170.5 188.0 235 288.86 255 226 194 172 q u.m () 552.00 570.00 633.00 715.58 685.65 655.00 620.12 571.23 (%) 22.50 22.92 23.03 23.36 23.00 22.75 22.59 22.10 M (%) 38.25 38.64 39.12 39.53 38.56 38.34 36.52 33.81 1200 IT_C, IT_, IT_, IT_tot (/m 2 ) 1000 800 600 400 200 IT_ Collector IT Mirror-est IT_Mirror ast IT_Total 0 9 10 11 12 13 14 15 16 Time (h) Figure 7 Variations of flux incident on the top cover of the collector over a day 080-9

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 900 S_C, S_, S_, S_tot (/m 2 ) 800 700 600 500 400 300 200 100 S_total S_C S_ S_ 0 9 10 11 12 13 14 15 16 Time (h) Figure 8 Variations of Insident solar flux in the absorber plate over a day ff M, ff (%) 50 45 40 35 30 25 20 15 10 5 0 ff. ith mirrors ff. ithout mirrors 9 10 11 12 13 14 15 16 Time (h) Figure 9 Variation of efficiencies ith mirror and ithout mirror over a day 5.CONCLUSIONS From the calculations of useful energy q u, flux absorbed S, heat loss from collector q l, flux incident on the top cover I T, and making approximation, the efficiencies average over 8 hours from 0900 to 1600 h for collector ith and ithout mirrors ork out to be 41.9 % and 22.8 % respectively. 080-10

ICCHT2010 5 th International Conference on Cooling and Heating Technologies, Bandung, Indonesia 911 December 2010 RFRNCS [1]. Sudia B., Himran S., Piarah. H., 2009, Flat-Plate Solar Collector Performance ith to-mirrors Concentrators, S2 Thesis, Mech. ngineering Dept., Unhas, Makassar [2]. Foniman H, et.al, Himran S., Jalaluddin, 2007,2008, The Comparison among the Performance of Flat Plate Collector ith One, To, and Three Covers, S1 Minithesis, Mech. ngineering Dept., Unhas, Makassar [3]. Sukhatme, S. P., Nayak J. K., 2008, Solar nergy, Principles of Thermal Collection and Storage, 3 rd d, Tata McGra Hill Private Limited, Ne Delhi [4]. Duffie J. A., Beckman. A., 1991, Solar ngineering of Thermal Process, J.. & Sons, Inc, Ne York 080-11