Engineering, Technology & Applied Science Research Vol. 7, No. 4, 2017, 1741-1745 1741 www.etasr.com Fakhim: An Investigation of the Effect of Different Nanofluids in a Solar Collector An Investigation of the Effect of Different Nanofluids in a Solar Collector Hamed Agabalaie Fakhim Department of Mechanical Engineering, Tabriz Branch, Islamic Azad University, Tabriz, Iran hamedfakhim@iaut.ac.ir Abstract—In this article, we examine the use of different nanofluids in a solar collector in a parabolic form. Temperature, thermal efficiency and outlet average temperature for a conventional parabolic collector and a collector with nanofluid are compared. The effect of various parameters (concentration ratio, volume fraction of nanoparticles, absorption and fluid velocity) are studied. Results are discussed and it is shown that nano fluid increases the efficiency of the collector. Keywords- nanofluid; water; solar collector; thermal efficiency; Numerical Analysis I. INTRODUCTION The importance of solar energy is widely recognized and is of increased significance for countries where favorable weather conditions permits increased exploitation through PVs and solar collectors [1-2]. However, the potential is often larger than the amounts of energy actually used. In Iran, the potential is calculated to be about 4,000 times its energy consumption [3]. Researchers have made different experiments and simulations to evaluate the performance of collectors. An experimental investigation can be in conjunction with the plate coverage [4]. The two-dimensional finite element method has been employed in order to investigate the collector’s performance [5]. Another numerical simulation of the flat plate solar collector in which water has streamed showed good agreement with experimental results [6]. An effort to optimize a parabolic solar collector was conducted in [7] that included a comprehensive mathematical model of thermal and optical performance conditions of the collector. A considerable improvement in reducing radiation convection losses was reported [7]. The natural flow inside a waveform solar collector with the use of two nanofluids, using the finite element method, was investigated in [8]. It was concluded that the silver-water nanofluid had a better effect on increasing the heat transfer rate [8]. In this paper, the effect of nanofluids (water-silver and silver-copper) on the performance of a parabolic concentrator collector is examined. To improve the efficiency of the system, the nanoparticles are added to the fluid in the absorbing tube and as a result direct heat transfer is possible. This study numerically models the flow of nanofluid in the solar collector and describes its difference with a common collector. II. DESIGN OF SOLAR SYSTEMS Solar collectors work out as a heat exchangers in which a form of energy, solar energy, converts to another form of energy [9]. In general, it can be said that solar collectors absorb the solar radiation energy and convert it to heat. An essential part of the set up is anti-corrosion and the use of water and/or water with some additives such as antifreeze [10]. There are various different schemes that can be roughly categorized as active and passive, direct and indirect, full and back gravity systems. Passive systems include the use of very few moving and no electronic equipment. In these systems, storage and heat circulation are supplied with Thermosiphon or cold input water pressure. In active systems, electricity is used to supply the pump and the controllers. These systems are also known as mandatory rotation and have better management potential and performance on absorbed solar heat. Thermal, economic, environmental efficiency and lifetime of solar water heaters of Thermosiphon type have been examined in various works (e.g. [11]). An additional issue is the choice of solar collectors. Flat plate solar collectors are widely used due to their lower price, however vacuum tube collectors gain increased attention due to their improved characteristics. III. MODEL PARAMETERS Figure 1 shows the basic scheme of a solar collector, a water storage tank, heat exchangers and back tubes. Collector includes the absorber plate, glass coverage and side plates. Cube is with 1.2 length and 120 mm lateral distance. Height of absorbent to glass is 240 mm and insulator thickness under the absorbent is 30 millimeters. Pipe diameter inside the collector is 1.2 mm which passes from the adsorbent. Tubes on the parts outside the collector and the water tank are all assumed to be adiabatic. In the exact modeling, production of a precise computational grid is required for the geometric range. The control volume forms should be in a way to reduce computational error. The accuracy of the computational results in fluid dynamics depends heavily on the grids. Further, regular Engineering, Technology & Applied Science Research Vol. 7, No. 4, 2017, 1741-1745 1742 www.etasr.com Fakhim: An Investigation of the Effect of Different Nanofluids in a Solar Collector computing grids that include cubes controls volume increase the accuracy of the results and reduce the number of grids. The level of radiation in experiments has been obtained by averaging the amount of radiant energy on the top and bottom of the absorber tube. To evaluate the thermal performance of vacuum solar water heaters, a closed system instead of an open system is used. The thermal performance of a vacuum tube solar water heater is obtained via (1) GA TCm c p    (1) In (1), m represents the rate of mass flow of water heater inside the tank (via incoming and outgoing tank), Cp represents the thermal capacity of the working fluid, ΔΤ represents the average increase in the temperature of the working fluid between the tank inlet and the outlet, Ac represents the tube absorber level and G represents the average solar radiation. For closed systems, the mass of the fluid inside the tank is used instead of m. 3D components related to the workspace are considered to only appear in the time of activity performance in simulation. Temporary activities were corresponding to the construction resources in three categories: materials, machinery or equipment and manpower. For each of the resources, an attributed space was designed that covered the main resource in the form of transparent and geometric volumes. The workspace of the construction resource was defined to produce the building components of one day and it is appeared daily and continuously besides building components and is hidden in the end of that day. Actually, after the construction of one day's workspace were hidden and then the corresponding volumes with 3D building components were appeared and remained to the end. In (2), T represents the increase in the average temperature of the fluid inside tank, G represents the average solar radiation per tube and time represents the time of the experiment. Gtime TCm pk , tan   (2) Fig. 1. Display of solar water heating system IV. MATHEMATICAL MODELING To evaluate and analyze the theory of parabolic collector, the radiation source, nanofluids and concentrator have been numerically modeled. The spectral distribution of radiation by a black object is a function of wavelength and temperature which has been defined by the Planck distribution, in (3) [12].              1exp 2 ),( 0 2 0 sumB sumb Tk hc hc TI     (3) In this equation, Tsum represents the solar surface temperature, h is Planck's constant, kB represents Boltzmann's constant, c0 represents the speed of light in a vacuum and λ the wavelength. Spectral radiation intensity is calculated via (3). Energy emission wavelength is considered between 0.2-6 μm. In this research, the adsorbent is modeled as geometric and in form of a circular channel which has placed along linear parabolic concentrator. Water tank has been designed 10 liters regarding single-tube collector. To simplify the model, it is assumed that the absorber has a transmittance coefficient equaled to 1. Nanofluid is modeled as a half-transparent environment, moving with constant speed. In this modeling, the radiation landing and then the energy transmission are considered through the absorption and Rayleigh scattering methods [13]. Received radiation in the water is due to absorption but received radiation caused by the interaction of nanoparticles is due to both scattering and absorption. The concentrator in Figure 1 converges the radiations caused by solar energy on adsorbent and a parameter called concentration ratio is defined to determine the convergence size [14]. For ease of severity, received radiation caused by solar energy is assumed in the radial direction. To receive and calculate radiation at various radiuses, the corrected form of Bayes equation in form of (4) is used. rIKrIKK dr rId esa ...).( ).(     (4) In (4), Kaλ represents the absorption coefficient of the spectrum, Ksλ represents the spectral scattering coefficient and Keλ represents the coefficient of received scattering. The tube inside collector has been made of copper with a thickness of one millimeter and a diameter of 1.27 cm. It has been assumed that the inner surface of the collector is lined with black tubes to increase the absorption of radiation; tube as opaque surface has the absorption coefficient of 0.94. The intensity of the radiation received in the nanofluid is caused by the effect of base fluid and nano-particles. When radiation passes through the water, received radiation takes place through absorption; received radiation takes place through scattering and absorbing incoming radiation in nanofluids. In the proposed model, the average diameter of the nanoparticles is about 8 nm thus the overall received coefficient is expressed in (5) via Riley equation. D Q K meesnanopartice 2 3 ),( )(      (5) Engineering, Technology & Applied Science Research Vol. 7, No. 4, 2017, 1741-1745 1743 www.etasr.com Fakhim: An Investigation of the Effect of Different Nanofluids in a Solar Collector In (5), ϕ represents the volume fraction of nanoparticles and Qeλ represents the received efficiency and D represents particle diameter. α represents particle size and m represents the relative refractive index of particles obtained in (6) [7-15].    D   (6) fluid particles n m m   (7) Received efficiency in computing is obtained via (8). 2 2 2 4 2 22 2 22 2 2 1 1 3 8 } 32 3827 1 1 15 1 1 1 {4                              m m m mm m m m m IQ me     (8) The ratio of the total received radiation for the nanofluid is obtained by the superposition of received base fluid and Nanofluid coefficient. Thus, in order to obtain the relationship between the radiation flows in the flow system, the placement of this amount in the Bayes equation is used. In the software, time-dependent equations in terms of both time and space are discrete. In (9), the amount of energy production (Egn) in the overall energy balance equation is expressed. .... stgenoutin EEEE   (9) In (9), Ėin, Ėout, Ėgen, Ėst represent the energy transmission to system, power transmission from system, power generation and energy storage respectively. To find the temperature field for the two-dimensional steady-state from the overall equation of equilibrium temperature, (9) is used. After obtaining the overall nanofluid temperature, average outlet temperatures, optical efficiency, thermal efficiency, etc. are calculated. A. Equations governing the flow To study and analyze the behavior of flows, the conservation of mass equations are used. For compressible flows or flows of heat transfer, energy conservation equations are solved. Under turbulent flow, turbulence modeling equations should be used. The purpose of calculating flows, especially turbulent flows, is to determine the terms such as Reynolds stress, turbulent mass flux or turbulent heat flux by linking the mentioned amounts to the average flow quantities and especially existing gradients at average flow [9]. Mass conservation equation is written in (10). 0)()( 11      uu xt i    (10) by considering (10), the equation of momentum is written in (11). VPguu xDt VD j i 2)(       (11) Energy equation is written in (12).                                  11 1 1 1 1 1 2 x u x u x u x u TuC x T k xDt TD C jj j p i p    (12) V. NUMERICAL SOLUTION To solve this set of equations, finite difference method is used, so that the total volume is decomposed to the initial controls volume with grids at their centre [10]. Absorbent with radius and length of 40 mm and 2.5 mm decomposes along radial direction. The size of the nodes along radial direction of axis r equals 0.5 mm and 1.2 mm. The size of the grids is selected so that solutions are stable and convergent in iterations from simulation. Fluid temperature incoming the environment and the collector temperature are considered 25°C. Further, convective heat transfer coefficient is equal to 6.75W/m2K. Incoming solar radiation flux to system is also considered be 1000W/m2K [11]. In this calculation, the diameter of the nanoparticles and the volume fraction (ϕ) have been considered 8 nm and 0.03 wt%, respectively. To obtain effective thermo- physical properties of nanofluid, the rule of mixture is used [12]. (1 )nanofluid particle basedfluidP P P     (13) VI. STUDY ON RESULTS The results of two-dimensional field for collectors with Nanofluid and conventional collectors are displayed in Figure 2 and the comparison of the two collectors is shown in Table I. The results of the simulation in Table I show that the nanofluid collector outperforms the conventional collector under the same working conditions. TABLE I. COMPARISON OF TWO COLLECTORS IN THE PROPOSED MODEL Optical efficiency (%) Thermal efficiency (%) Average outlet temperature (K) Working fluid 322.5 27.8 33.2 pure water 331.72 95.2 97.8 Nanofluid To optimize the proposed model, working parameters such as concentration ratio, volume fraction, the adsorbent and the flow rate were varied and their impact on the performance index is evaluated. The impact of nanofluid volume fraction (Ag/Water, Al/Water) on the temperature of the fluid in the tube is shown in Figures 3 and 4. As shown, the temperature Engineering, Technology & Applied Science Research Vol. 7, No. 4, 2017, 1741-1745 1744 www.etasr.com Fakhim: An Investigation of the Effect of Different Nanofluids in a Solar Collector difference will be about 2 degrees. If different percentage of nanoparticles is added to the base fluid, inlet and outlet temperature difference with respect to the used volume percent varies. Fig. 2. Results from the two-dimensional field for collectors with nanofluid and conventional collectors Fig. 3. The impact of the volume fraction of Nanofluids (Ag/Water) on the temperature of the fluid in the tube Fig. 4. The impact of the volume fraction of Nanofluids (Al/Water) on the temperature of the fluid in the tube Considering Table I, it is specified that independent optical and thermal efficiency are the ratio of concentration, but the outlet temperature changes linearly with concentration ratio as shown in Figure 5. A change of the volume fraction has a significant impact on all three performance indicators. As shown in Figure 6, the average outlet temperature increases rapidly with increasing volume fraction to reach a certain speed and then it is almost constant. Indeed, there is a certain amount of volume fraction, after which certain progress does not occur in the average temperature outlet. Figure 7 shows the comparison of the temperature of the absorber plate and the working fluid along the tube using nanofluid with the same volume percentage. In case of using water as the fluid, the difference in the inlet and outlet temperature of water has the least temperature, and also the temperature difference between the beginning and end of the absorber plate has the highest amount, so that we will witness less efficiency using water as fluid. In case of nanofluid, there will be the highest difference on inlet and outlet temperature and lowest difference on the absorber plate. The Black lines display fluid temperature with nanoparticles and red lines display the temperature of the flowing water. Fig. 5. Average outlet temperature of nanofluid based on collector concentration ratio Fig. 6. Average outlet temperature in terms of volume fraction of nanoparticles Fig. 7. Comparison of the temperature of the absorber plate and the working fluid along the tube As shown in Figure 8, the average outlet temperature increases by increasing the absorber. This behavior is because of the fact that the duration at which a specific amount of working fluid passes absorber increases by constant flow speed and increased absorber, the fluid has more time to warm up along tube. On the contrary to Figure 8, thermal efficiency reduces by increasing the absorption. This increases due to increased absorber and thereby increased the time spent in the working fluid and heat loss increases through convection. On the other hand, by increasing the flow rate, the elapsed time in the adsorbent decreases, or in other words, less time is spent to increase temperature through radiation. Thus, with increasing speed, the average temperature of outlet reduces. By increasing the flow rate of heat loss, convection reduces and thermal efficiency increases (Figure 9). Engineering, Technology & Applied Science Research Vol. 7, No. 4, 2017, 1741-1745 1745 www.etasr.com Fakhim: An Investigation of the Effect of Different Nanofluids in a Solar Collector Fig. 8. Average outlet temperature in terms of absorber Fig. 9. Heat efficiency in terms of flow speed Fig. 10. Average outlet temperature in terms of collector concentration ratio. VII. CONCLUSION In this research, the effect of different Nanofluids (Ag/Water, Al/Water) on a parabolic solar collector performance is discussed. The results of the simulations demonstrate that adding a small amount of nanoparticles, not only improves its absorption characteristics, but also improves the efficiency of the collector. An increase in the volume fraction will increase collector efficiency and fluid speed to a certain range. Further, by increasing the absorption, average outlet temperature increases but thermal efficiency decreases. On the contrary, this governs the speed increase so that the average temperature was reduced by increasing flow rate, but the heat efficiency increases. Overall, using nanofluids as the base fluid is shown to increase collector efficiency and outlet temperature. VIII. REFERENCES [1] A. W. Badar, R. Buchholz, F. Ziegler, “Experimental and theoretical evaluation of the overall heat loss coefficient of vacuum tubes of a solar collector”, Solar Energy, Vol. 85, pp.1447–1456, 2011 [2] I. Budihardjo, G. L. Morrison, M. Behnia, “Natural circulation flow through water in glass evacuated tube solar collectors”, Solar Energy, Vol. 81, pp. 1460–1472, 2007 [3] F. Jorge, O. Armando, “Numerical simulation of a trapezoidal cavity receiver for a linear Fresnel solar collector concentrator”, Journal of Renewable Energy, Vol. 36, pp. 90–96, 2011 [4] K. Sopian, M. Syahri, S. Abdullah, M. Y. Othman, B. Yatim, “Performance of a non-metallic unglazed solar water heater with integrated storage system”, Renewable Energy, Vol. 29, pp.1421–1430, 2004 [5] R. S. R. Gorla, “Finite element analysis of a flat plate solar collector”, Finite Elements in Analysis and Design, Vol. 24, pp. 283-290, 2010 [6] M. Selmi, M. J. Al-Khawaja, A. Marafia, “Validation of CFD simulation for flat plate solar energy collector”, Renewable Energy, Vol. 33, pp. 383-387, 2008 [7] A. Kahrobaian, H. Malekmohammadi “Exergy Optimization Applied to Linear Parabolic Solar Collectors”, Journal of Faculty of Engineering, Vol. 42, No. 1, pp. 131-144, 2008 [8] T. P. Otanicar, P. E. Phelan, R. S. Prasher, G. Rosengarten, R. A. Taylor, “Nanofluid-based direct absorption solar collector”, Journal of Renewable and Sustainable Energy, Vol. 2, No. 3, p. 033102, 2010 [9] T. Yousefi, F. Veysi, E. Shojaeizadeh, S. Zinadini, “An experimental investigation on the effect of pH variation of MWCNT–H2O nanofluid on the efficiency of a flat-plate solar collector”, Solar Energy, Vol. 86, No. 2, pp. 771-779, 2012 [10] M. F. Modest, Radiative Heat Transfer, Academic Press, 2003 [11] A. Kasaeian, A. Mohmmadkarim, A. Kaabinejadian, “Performance investigation of solar evacuated tube collector using TRNSYS in Tehran”, International Journal of Renewable Energy Research, Vol. 4, No. 2, pp. 497-503, 2014 [12] F. P. Incropera, D. P. De Witt, Fundamentals of Heat and Mass Transfer, Wiley and Sons, New York, 2007 [13] H. Tyagi, P. Phelan, R. Prasher, “Predicted efficiency of a Low- temperature Nanofluid-based direct absorption solar collector”, ASME Journal of Solar Energy Engineering, Vol. 131, No. 4, pp. 0410041- 0410047, 2009 [14] S. K. Das, S. U. Choi, W. Yu, T. Pradeep, Nanofluids: Science and Technology, Wiley and Sons, New York, 2008 [15] Z. Li, C. Chen, H. Luo, Y. Zhang, Y. Xue, “All-glass vacuum tube collector heat transfer model used in forced-circulation solar water heating system”, Solar Energy, Vol. 84, pp. 1413-1421, 2010