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

Similar documents
Heat Transfer Enhancement of Solar Flat Plate Collector by Using V Corrugated Fins and Various Parameters

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

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

Solar Flat Plate Thermal Collector

UNIT FOUR SOLAR COLLECTORS

Simplified Collector Performance Model

International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN

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

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

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

Exergy Analysis of Solar Air Collector Having W Shaped Artificial Roughness

Numerical Study of the Three-dimensional Flow in a Flat Plate Solar Collector with Baffles

Chapter 5 MATHEMATICAL MODELING OF THE EVACATED SOLAR COLLECTOR. 5.1 Thermal Model of Solar Collector System

Vertical Mantle Heat Exchangers for Solar Water Heaters

OPTIMIZATION of the GEOMETRY & MATERIAL of SOLAR WATER HEATERS.

1D and 3D Simulation. C. Hochenauer

Numerical Analysis of Heat Transfer Performance of Flat Plate Solar Collectors

A CFD Analysis Of A Solar Air Heater Having Triangular Rib Roughness On The Absorber Plate

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

Optimization of Shell And -Tube Intercooler in Multistage Compressor System Using CFD Analysis

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

Performance Assessment of PV/T Air Collector by Using CFD

HEAT LOSS CHARACTERISTICS OF A ROOF INTEGRATED SOLAR MICRO-CONCENTRATING COLLECTOR

DETAILED MODELING OF SOLAR FLAT-PLATE COLLECTORS WITH DESIGN TOOL KOLEKTOR 2.2

International Journal of Advanced Engineering Technology E-ISSN

CFD ANALYSIS FOR INTEGRATED COLLECTOR STORAGE SOLAR WATER HEATERS (ICSSWH) BASED ON CUBOID TYPE GEOMETRIES

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

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER

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

Numerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led Application by using CFD

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

ANALYSIS OF FLAT PLATE PHOTOVOLTAIC-THERMAL (PVT) MODELS

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE

CFD ANALYSIS OF TRIANGULAR ABSORBER TUBE OF A SOLAR FLAT PLATE COLLECTOR

Project #1 Internal flow with thermal convection

STUDY OF FLAT PLATE COLLECTOR SOLAR WATER HEATER

PERFORMANCE ANALYSIS OF PARABOLIC TROUGH COLLECTOR TUBE WITH INTERNAL INTERMITTENT FINS

A Computational Fluid Dynamics Investigation of Solar Air Heater Duct Provided with Inclined Circular Ribs as Artificial Roughness

Heat Transfer Enhancement in Vertical Narrow Plates by Natural Convection

Comparative Analysis of Heat Transfer and Friction Characteristics in a Corrugated Tube

Comparative study of Different Geometry of Ribs for roughness on absorber plate of Solar Air Heater -A Review

Heat Transfer Enhancement of Solar Air Heater By Using Artificial Roughness double inclined ribs

Effect of Periodic Variation of Sol-air Temperature on the Performance of Integrated Solar Collector Storage System

Exam questions: HEAT. 2. [2003 OL][2004 OL][2005 OL][2006 OL][2007 OL][2008 OL][2009] Name two methods by which heat can be transferred.

The energy performance of an airflow window

CFD Analysis of High Temperature and High Velocity System: Plasma Torch

COMBINATION OF TAGUCHI METHOD AND ARTIFICIAL INTELLIGENCE TECHNIQUES FOR THE OPTIMAL DESIGN OF FLAT-PLATE COLLECTORS

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

C ONTENTS CHAPTER TWO HEAT CONDUCTION EQUATION 61 CHAPTER ONE BASICS OF HEAT TRANSFER 1 CHAPTER THREE STEADY HEAT CONDUCTION 127

Autumn 2005 THERMODYNAMICS. Time: 3 Hours

Simulation of Free Convection with Conjugate Heat Transfer

International Journal of Engineering Research and General Science Volume 3, Issue 6, November-December, 2015 ISSN

Numerical Analysis of Fe 3 O 4 Nanofluid Flow in a Double Pipe U-Bend Heat Exchanger

SIMULATION OF A DOUBLE CAVITY TYPE SOLAR TROUGH COLLECTOR

COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES OF DIFFERENT ASPECT RATIOS

Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel

Parametric Effect on Performance Enhancement of Offset Finned Absorber Solar Air Heater

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

STUDY OF A PASSIVE SOLAR WINTER HEATING SYSTEM BASED ON TROMBE WALL

The Electrodynamics of a Pair of PV Modules with Connected Building Resistance

Faculty of Technology and Science Department of Physics and Electrical Engineering. Helena Johansson. Nocturnal cooling

IJSRD - International Journal for Scientific Research & Development Vol. 4, Issue 05, 2016 ISSN (online):

EXPERIMENTAL AND NUMERICAL STUDIES OF A SPIRAL PLATE HEAT EXCHANGER

Computational Fluid Dynamics Based Analysis of Angled Rib Roughened Solar Air Heater Duct

CFD model to estimate the Effect of Tilt and Height on the Natural Air Flow inside a Solar Chimney

Flow and Temperature Analysis inside Flat Plate Air Heating Solar Collectors

COVENANT UNIVERSITY NIGERIA TUTORIAL KIT OMEGA SEMESTER PROGRAMME: MECHANICAL ENGINEERING

PHYSICAL MECHANISM OF CONVECTION

COMPUTATIONAL FLUID DYNAMICS ON DIFFERENT PASSAGES OVER A PLATE COIL EVAPORATOR FOR 40 LITER STORAGE TYPE WATER COOLER

HEAT TRANSFER THERMAL MANAGEMENT OF ELECTRONICS YOUNES SHABANY. C\ CRC Press W / Taylor Si Francis Group Boca Raton London New York

Performance Investigation of Cavity Absorber for Parabolic Dish Solar Concentrator

FLOW AND HEAT TRANSFER ANALYSIS OF VARIOUS RIBS FOR FORCED CONVECTION HEAT TRANSFER

AR/IA 241 LN 231 Lecture 4: Fundamental of Energy

Numerical Study of Convective Heat Transfer for Flat Unglazed Transpired Solar Collectors

Investigation of Heat Transfer Enhancement in Laminar Flow through Circular Tube Combined Wire Coil and Wavy Strip with Central Clearance

Chapter 1 INTRODUCTION AND BASIC CONCEPTS

Validation, Optimization and Simulation of Solar Thermoelectric Generator Model

Thermal Analysis of Shell and Tube Heat Ex-Changer Using C and Ansys

Applied CFD Project 1. Christopher Light MAE 598

Comparative Performance Analysis of Conventional Solar Air Heater with Longitudinal Rectangular and Triangular Fins Attached Absorber

Tick the box next to those resources for which the Sun is also the source of energy.

Simulation of a linear Fresnel solar collector concentrator

EFFECT OF STAGGERED RIB LENGTH ON PERFORMANCE OF SOLAR AIR HEATER USING V-RIB WITH SYMMETRICAL GAP AND STAGGERED RIB

Keywords: Spiral plate heat exchanger, Heat transfer, Nusselt number

Analysis of the Cooling Design in Electrical Transformer

Coolant. Circuits Chip

CFD study for cross flow heat exchanger with integral finned tube

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

Enhancement of Heat Transfer Effectiveness of Plate-pin fin heat sinks With Central hole and Staggered positioning of Pin fins

NUMERICAL HEAT TRANSFER ENHANCEMENT IN SQUARE DUCT WITH INTERNAL RIB

Heat Loss from Cavity Receiver for Solar Micro- Concentrating Collector

Simulation Application for Optimization of Solar Collector Array

HEAT TRANSFER AND FRICTION FACTOR CHARACTERISTICS OF COMPOUND TURBULATORS IN RECTANGULAR DUCTS OF SOLAR AIR HEATER

INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY

Chapter 3 NATURAL CONVECTION

FLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study

arxiv: v1 [physics.app-ph] 25 Mar 2018

Numerical simulation of heat transfer and fluid flow in a flat plate solar collector with TIM and ventilation channel

CFD AND CONJUGATE HEAT TRANSFER ANALYSIS OF HEAT SINKS WITH DIFFERENT FIN GEOMETRIES SUBJECTED TO FORCED CONVECTION USED IN ELECTRONICS COOLING

Numerical investigation to study effect of radiation on thermal performance of radiator for onan cooling configuration of transformer

Transcription:

HEAT TRANSFER ENHANCEMENT OF SOLAR FLAT PLATE COLLECTOR BY USING V CORRUGATED FINS AND VARIOUS PARAMETERS Manoj S. Chaudhari 1, Mahesh Jagadale 2 1,2 Department of Mechanical Engineering, SPP University, GSMCOE Balewadi, Pune (India) ABSTRACT Solar flat plate collector can be defined as typical metal box covered with a coated glass (known as glazing) on the top and dark colored absorber plate at the bottom. It has insulation on sides and bottom in order to minimize the heat loss. In this project work, initially an experimental set up is made for studying the flow through domestic solar flat plate collector and record experimental observations. After completion of experiment a 3D mathematical model with single and double glazed covers of Solar Flat plate collector (SFPC) is prepared and model is analyzed for CFD flow and heat transfer. Moreover, in this project use of vortex generators (V corrugation over absorber plate) and heat enhancement setups such as double glazing, optimum air gap, use of different gases(air, carbon dioxide and argon) are considered so as to increase the heat transfer rate between the interacting fluids. The basic geometrical dimensions for the 3D mathematical model is based on domestically available SFPC on which experiment was performed and CAD model is prepared using Ansys Design Modeler. The prepared model is simulated in ANSYS Fluent 14.5 considering radiation, natural convection and conduction all modes of heat transfer in boundary conditions. A parametric study is done so as to understand the effect of air gap thickness and use of different gases on the heat transfer rate. A comparative study is carried on between models having V corrugation and model having double glazing and another similar model but without fins and glazing, it was concluded that Heat transfer Enhancement occurs for models with double glazing and also with the use of v corrugation over absorber plates. Keywords: Glazing, Temperature, Fin, ANSYS Design Modular, ANSYS Fluent, Heat Transfer Enhancement I. INTRODUCTION With increasing living standards, growing world population and rapid developments in technology, raises the need for energy every day, however, meeting the demand is harder due to the limited energy reservoirs existed around the world. The environmental pollution resulting from the usage of conventional energy sources is another serious problem. Due to these difficulties related to the energy issue, the world is seeking for alternative energy sources If we look back to the 19th century, people used stoves to heat water by burning pieces of coal or wood. Duffie J. and Beckmann W[1] suggested a remedy to avoid these problems, a much safer and easier 355 P a g e

way to heat water was invented. This was achieved by placing a black painted metal tank full of water outside to absorb solar energy. Klein and Duffie [2] studied effects of thermal capacitance while modelling SFPC.They found that capacitance effect is significant due to variable whether conditions. They investigated three SFPC models first model of zero capacitance, the second model accounts for effects of capacitance assuming collector as unit and third model dividing the collector in nodes and isothermal segments. Wang and WU [3] proposed discrete numerical model for flow distribution and in solar collector. They considered no uniform flow, conduction and effects of buoyancy and found good agreement with experimental results. Oliva et al. [4] introduced another numerical method to determine solar collector thermal behavior. The distributed character model considers multidimensional and transient properties of heat transfer that characterize the solar collector, while of heat transfer flux by free convection at the air gap zone has been calculated using empirical expressions and the solar irradiance was integrated to be constant per hour. Anderson [5] scrutinized different coloured solar collector performances. Based on the product of transmittance absorptance of collectors of different colours, theoretical performances of the used collectors were determined using the Hottel-Whillier-Bliss 1-D steady-state model presented by Duffie and Beckmann. Gertzos and Caouris [6] enquired into a different heat transfer enhancement approach which used a pump for circulating the storage water in the indirect heating flat plate collector. This resulted in an increase in the outlet service water temperature. For flat-plate solar collector, Cadaflach [7] has presented a detailed numerical model. He noticed that the transfer of heat through the collector is essentially one dimensional. 2-D and 3-D effects almost always occur due to the influence due the edges and the non-uniform effects. A polymer solar collector developed by Martinopoulos et al. [8], where the solar energy is directly absorbed by the black-colored working fluid. The model was examined both experimentally and with computational fluid dynamics (CFD). Kim and Kumar [9] showed the investigations in order to enhance the Nusselt number and pressure drop by the use of rib baffles of different shapes, sizes, and orientations to produce artificial roughness in oreder to create more friction. Kiran kumar and Ranga Babu [10] investigated numerical and experimental analysis of SFPC in order to increase its thermal efficiency. The study made use of nano fluids like Al 2 O 3 and Cu 2 O with water and predicted heat loss coefficients. From the above literature survey, it has been concluded that so much investigation has been done analytically and experimentally for enhancing the thermal performance of solar flat collector and these results are interpreted and validated against the CFD simulation. However, many of these studies are considered constant heat flux condition and/or 1D/2D model. Hence, to understand the heat transfer phenomenon thoroughly, solar flat plate collector is modeled three dimensionally and actual solar radiation (which depends on geographical location, day and time) falling on the collector is taken into account. The aim of this thesis is to carry out parametric study on various parameters on a proposed design of a simple flat-plate glazed collectors suggesting possible design modifications which could prove healthy in future for increasing the performance and overall efficiency. II. PROBLEM STATEMENT A 3D numerical model for flat-plate liquid solar collectors has been designed in ANSYS 14.5 design module. This model is used for enhancing heat transfer the collector, for which the different heat transfer mechanisms involved are simultaneously taken into account and used to reduce the overall heat loss: solar radiation 356 P a g e

absorption, transmission and reflection natural convection in the air cavity heat conduction across the tubeabsorber welded junction mixed convection flow in the risers and heat losses by convection and radiation to the ambient. To ensure the reliability of the model, the heat transfer results inside the risers and in the air cavity were contrasted with experimental results. III. OBJECTIVE Objectives of this thesis are as follows- 1)Experimental set up and recording observations from domestically used Solar Flat plate collector. 2) A 3D geometric model of flat-plate solar liquid collectors has been developed in ANSYS Design Modular and has been solved in ANSYS 14.5 (Fluent). 3) To validate the simulated results with experimental values. 4) To observe the effects of variation of air gap thickness, Use of different gases between absorber plate and glass cover like air, CO2, argon gas, use of V corrugated absorber plate instead of flat plate and comparison of single and double glazing surface. IV. EXPERIMENTATION In the part of experimentation, various types of steps or parts are included as: A. Material and Specifications The specification is mentioned below: Table 1: Solar Flat Plate Collector Specification Solar flat plate collector Solar flat plate collector cover plate made of toughened glass. components Sheet for absorber made of copper. Absorber made of copper sheet and copper tube a)glass cover Cover plate: cover plate is toughened glass and thickness of 5.0 mm b)absorber plate Sheet for absorber is made from copper,thickness -1.2mm c) Air gap The air gap between the glazing and the absorber is 8 mm d)inner tank or cold tank Material- Stainless steel, 2.4m from terrace level e)gasket for flanges 3mm thick gasket of rubber is used for ceiling between flanges f)collector Support Collector can withstand wind velocity of 100 km/hr frame g)storage Tank or Hot Made of Stainless steel water tank h) Insulation Material- Rockwool I)Piping Material- Medium class (B class) GI as per IS 1239 357 P a g e

B. Experimental Setup: Fig : Set up of SFPC with six series pipe The above figure shows setup of SFPC with cold water tank and hot water tank and six series pipe. The flow in SFPC is through natural convection that is through siphon action water in pipe flows from cold water tank to hot water tank, thus no requirement of pump. In this project our main aim is to increase the heat transfer to pipes for increasing water temperature at outlets. Table 1: Construction details of Solar FPC with their material properties based on observation Eleme Mate Thermal De Specif Dimension (mm) nt rial Conduct ivity (W/m- K) nsit y (Kg /m 3 ) ic heat (J/Kgk) Pipe Copp er 387.6 897 8 381 I.D.=11.58 mm O.D.=12.7 mm & Length =800 mm Cover Glass 1.14 223 0 750 Length=805mm, width = 800 mm& thickness= 5mm Absorb er Copp er 387.6 897 8 381 Length=805 mm, Width =800 mm& Thickness =1.2 mm Insulati on Rock wool 0.04 48 840 Bottom side=50 mm Edge side= 25 mm 358 P a g e

Experimental Observations S.No Total flux on surface Mass flow rate Outlet Time Inlet temperature(ºc) I T(W/m 2 ) (Kg/sec) temperature(ºc) 1 10 AM 782.461 0.0063 27 30.2 2 11 AM 814.3508 0.0063 28.3 32.7 3 12 PM 826.4204 0.0063 30.1 35.2 4 13 PM 821.8579 0.0063 29.8 34.9 V. FINITE ELEMENT ANALYSIS A numerical method has an important contribution to solving complex computational mechanics problems, quickly and correctly. There are many methods generally used to solve a variety of problems. Each method has their own advantages and limitations. Finite element method is one which used differential equations to solve most of the problems related to automobile, civil, aerospace industries. The CFD simulations are performed by the CFD package in the ANSYS-14.5. To gain simulated results using CFD software, a simulation procedure needs to be followed. The procedure requires setting the volume and boundary conditions of the simulated module. A sketch of the geometry which was investigated is depicted in figure 4.1. It shows a standard flat-plate solar thermal collector where four different zones are identified: (1) the sheet-and-tube assembly, which consists of six 11.58 mm ID tubes, separated 115 mm and welded to a copper absorber plate along their upper generatrix. (2) A glass sheet which is transparent, which covers the absorber plate and reduces the convection and radiation losses to the ambient. (3) The air filled cavity between the glass and the absorber plate at atmospheric pressure. (4) An insulating layer of rockwool that covers the lateral and bottom sides of the assembly to minimize the conduction and convection losses to the surroundings. Fig.4 Geometry of solar flat plate collector model 359 P a g e

Fig.5 Meshing Fig.6 Temperature Contour at 10 AM Fig.7 Temperature Contour at 12 PM Fig.8 Streamline flkow in pipe The maximum velocity recorded at pressure outlet was found to be 1.87m/sec, velocity at outlet reduces from the inlet velocity of 2m/sec due to frictional losses and turbulence. The streamline can be seen in the above figure in order to know the profile of flow of water through pipes from velocity inlet to outlet VI. RESULT AND DISCUSSION In order to verify the accuracy of the calculation results, the simulation values are compared with the values obtained from Experimental observation.outlet fluid temperature is a measure of net heat gain from incident solar heat flux. By absorbing the incident solar heat flux on particular day at different timings, the outlet fluid temperature is varying. Thus below figure shows variation of outlet temperature of SFPC experimentally and numerically. Fig.9 Experimental versus Simulation Results 360 P a g e

Use of V corrugation over absorber plate In order to compare normal solar flat plate collector with V corrugated absorber plate collector a model of new Solar flat plate collector has been generated with v corrugations over absorber plate. The below figure shows the model of SFPC designed with V corrugation over absorber plate. Fig.10 Solar FPC with V corrugated absorber plate In order to avoid high number of nodes and meshed elements, for saving computation time insulation has not been considered. The insulation wall boundary condition has been provided below the pipes. Based on the results obtained of outlet temperature it was found that SFPC with V corrugation over absorber plate are more efficient then normal FPC. The comparison of outlet temperature is shown below. Fig.11 Normal SFPC versus SFPC with Fins Effect of Air Gap thickness When the size of the air gap was very small, the heat transfer from glass to absorber plate is dominated by conduction. As the gap between the glass and absorber plate is increased, natural convection takes place which leads to an increase in the convective heat loss. From the figure 5.8 the top loss coefficient is increased with the air gap up to some extent but after certain limit of gap net heat gain in terms of outlet fluid temperature will decrease because when thickness of gap will increase, shadow of side wall will be created on the absorber plate so net heat flux will be concentrated on certain area of absorber plate. 361 P a g e

Fig.12 : Effects of gap between absorber plate and cover Effect of different gases used between absorber plate and glass cover: Fig.13 Overall loss coefficients using different fluids Figure 13 shows the effect of different gases between the gap of absorber plate and glass cover. Here, air, carbon dioxide and inert gas like argon have been used. We can observe that in case of air overall heat loss coefficient is more for carbon dioxide and air and argon gas has the lowest heat losses. Argon (inert gas) has more absorption capacity related to temperature at different wavelength that reduces the heat loss. Use of Double Glazing As we can observe from figure below in case of double glazing i.e. two glass cover outlet fluid temperature is more than single glazing. In case of double glazing air gap is 20 mm provided between top glass cover and second glass cover. In this case more heat is trapped inside the glass cover that will reduce the heat losses hence outlet fluid gets more heated 362 P a g e

Fig.14 Comparison of Double and Single Glazing By using double glazing more radiation is trapped and re-radiation losses would be minimized, hence the overall loss coefficient is more in case of single glazing is more compared to the double glazing VII. CONCLUSION From the given study it is observed that The maximum difference between numerical and experimental results is less than 4%. Based on comparative study of normal SFPC and SFPC with V corrugated absorber plate, it was observed that SFPC with V corrugation has much more heat transfer due to increased turbulence(enhancing heat transfer through convection) and reflectance(enhancing heat transfer through radiation reflection and absorption). The size of the gap between the absorber plate and glass cover is investigated and it is found that there is certain limit for gap and it is good to have small air gap. When different gases like air, carbon dioxide and argon have been used. We observed that in case of air overall heat loss coefficient is more for carbon dioxide and while carbon dioxide has more heat loss coefficient then argon and argon gas has the lowest heat losses. Argon (inert gas) has more absorption capacity related to temperature at different wavelength that reduces the heat loss. REFERENCES 1. Duffie J. and Beckmann W., Solar engineering of thermal processes, 2nd edition Wiley Interscience, New York 1991. 2. Klein S., Duffie J., and Beckman W., Transient considerations of flat-plate solar collectors, Trans. of ASME: Journal Engineering for Power vol. 1974 96A, pp. 109-113. 3. Wang X. and Wu L., Analysis and performance of flat-plate solar collector arrays, solar energy, 45, 1974, 2, 71-78. 4. Oliva A., Costa M., and Perez Segarra C., Numerical simulation of solar collectors: the effect of non uniform and unsteady state of boundary conditions, Solar Energy, 1991, 47, 5, 359-373. 363 P a g e

5. Anderson T., Duke M., and Carson J., The effect of color on the thermal performance of building integrated solar collectors, Solar Energy Materials & solar cells, 2010, 350-354. 6. K.P. Gertzos, Y.G. Caouris, Experimental and computational study of the developed flow field in a flat plate integrated collector storage (ICS) solar device with recirculation, experimental Thermal and Fluid Science 31(8) 2007, 1133-1145. 7. Cadafalch J., A detailed numerical model for flat-plate solar thermal devices, Solar Energy Vol. 83, 2009, pp. 2157-2164. 8. Martinopoulos G., CFD modeling of a polymer solar collector, Renewable Energy Vol. 35, 2010, pp. 1499-1508. 9. Kumar, A. and Kim, M.H., Convective heat transfer enhancement in solar air channels. Applied Thermal Engineering, 89, pp.239-261.2015. 10. Rangababu, J.A., Kumar, K. and Rao, S., 2015. Numerical Analysis and Validation of Heat Transfer Mechanism of Flat Plate Collectors. Procedia Engineering, 127, pp.63-70. 364 P a g e