Microsoft Word - 0_Vol_16_No_2_2015_Special Issue on Energy_Compiled.docx IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 13 ANALYSIS AND MODELING OF SOLAR EVAPORATOR-COLLECTOR ZAKARIA MOHD. AMIN1, M.N.A. HAWLADER2 AND YE SHAOCHUN3 1School of Chemistry, Physics and Mechanical Engineering, Queensland University of Technology, Australia. 2 Department of Mechanical Engineering, Faculty of Engineering, International Islamic University Malaysia, Jalan Gombak, 53100, Kuala Lumpur, Malaysia. 3Department of Mechanical Engineering, National University of Singapore, 9 Engineering Drive 1, Singapore 117576. zakaria.amin@qut.edu.au, mehawlader@iium.edu.my, and ye_shaochun@nmc.a- star.edu.sg (Received: August 10, 2015; Accepted: October 30, 2015; Published on-line: November 30, 2015) ABSTRACT: Solar energy is considered a sustainable resource that poses little to no harmful effects on the environment. The performance of a solar system depends to a great extent on the collector used for the conversion of solar radiant energy to thermal energy. A solar evaporator-collector (SEC) is basically an unglazed flat plate collector where refrigerants, such as R134a is used as the working fluid. As the operating temperature of the SEC is very low, it utilizes both solar irradiation and ambient energy leading to a much higher efficiency than the conventional collectors. This capability of SECs to utilize ambient energy also enables the system to operate at night. This type of collector can be locally made and is relatively much cheaper than the conventional collector. At the National University of Singapore, the evaporator-collector was integrated to a heat pump and the performance was investigated for several thermal applications: (i) water heating, (ii) drying and (iii) desalination. A 2-dimensional transient mathematical model of this system was developed and validated b y experimental data. The present study provides a comprehensive study of performance. ABSTRAK: Tenaga suria merupakan sumber lestari yang tidak memudaratkan alam sekitar. Keupayaan sistem ini bergantung sebahagian besarnya kepada kebolehan pengumpul untuk menukar tenaga suria kepada tenaga haba. Pengumpul penyejat suria (SEC) pada asasnya merupakan plat datar tanpa glasair di mana bahan penyejuk seperti R134a digunakan sebagai bendalir kerja. Memandangkan suhu kendalian SEC adalah amat rendah, ia boleh menggunakan penyinaran suria dan tenaga ambien agar memberikan efisiensi yang lebih tinggi daripada pengumpul konvensional. Keupayaan SEC dalam mempergunakan tenaga ambien membolehkan sistem ini beroperasi pada waktu malam. Pengumpul merupakan keluaran tempatan dan ianya lebih murah berbanding pengumpul lazim. Di Universiti Nasional Singapura, pengumpul penyejatnya diintegrasikan pada pam haba dan prestasinya dikaji untuk beberapa aplikasi terma: (i) pemanasan air, (ii) pengeringan dan (iii) penyahgaraman. Model matematik transien dua dimensi bagi sistem ini dikembang dan divalidasi berdasarkan data eksperimental. Selidikan terkini memberikan kajian prestasi yang menyeluruh. KEYWORDS: heat pump; evaporator-collector. 1. INTRODUCTION In view of the growing global energy needs and concern for environmental degradation, the possibility of running thermal systems using energy from the sun is IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 14 receiving considerable attention in recent years. Solar energy is clean and is the most inexhaustible of all known energy sources. The low temperature thermal requirement of a heat pump makes it an excellent match for the use of solar energy. The combination of solar energy and heat pump system can bring about various thermal applications for domestic and industrial use, such as water heating, solar drying, space cooling, space heating, and refrigeration. Unlike thermosyphon solar water heaters, solar heat pump systems offer the opportunity to upgrade low-grade energy resources from the surroundings as well as solar energy and makes use of it for domestic and industrial applications [1]. The concept of direct expansion solar-assisted heat pump (DX-SAHP) was first proposed in an experimental study by Sporn and Ambrose [2]. They used double glazed collector which also act as evaporator for the heat pump with R-12 being the working fluid. They did not demonstrate the full potential of the concept due to a mismatched, oversized compressor. They also did experiments by removing the glazing and found that the removal did not affect the performance significantly. They found that removal of glazing or back insulation does not affect performance of two-phase flow solar collector significantly. The collector operating temperature in a solar-assisted heat pump system can be lower than the ambient temperature. In this case, an un-glazed solar evaporator collector is used. The Simple structure of the un-glazed solar collector makes it an economical type of solar collector. However, its performance is highly dependent upon the environment, because its surface is exposed directly to the ambient. To improve the performance of un- glazed solar collectors, a clear understanding of how the environment affects it is required. An uninsulated, bare two-phase collector is firstly used in Franklin et al.’s [3] heat pump system. Chaturvedi et al.[4, 5] found a variation of the evaporator temperature from 0°C to 10°C above the ambient temperature under favorable solar conditions. Many authors [6-8] reported that for ambient temperatures above 25°C, the evaporator could be operated at an elevated temperature. Hawlader et al.[1] performed analytical and experimental studies on a solar assisted heat pump using unglazed, flat plate solar collectors, known as an evaporator-collector, using a steady-state one-dimensional mathematical model. Two-phase flow in collectors was first considered by Sporn and Embrose [2]. They used a double glazed collector which also act as evaporator for the heat pump with R-12 being the working fluid. They did not demonstrate the full potential of the concept due to a mismatched, oversized compressor. They also did experiments by removing the glazing and found that the removal did not affect the performance significantly. They found that removal of glazing or back insulation does not significantly affect performance of the two- phase flow solar collector. Literature on two-phase collectors report that two models of calculating pressure drop are most widely used. These are known as the Martinelli Nelson's method for separated flows and Owen's homogeneous equilibrium model for misty or bubbly flow [9]. The homogeneous equilibrium model is the simplest method determining the pressure drop in two-phase flow and makes the basic assumption that the two-phases have the same velocity. By adopting an equilibrium, homogeneous, two-phase model as mentioned above, Chaturvedi et al. [10] carried out preliminary theoretical performance studies concerning a solar-assisted heat pump that uses a bare collector as the evaporator. The analysis was IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 15 subject to limitation of a constant temperature evaporator with no superheating or sub cooling. It is the first time wherethe single phase flat plate collector theory was extended to the two-phase regime. The first comprehensive work about two-phase solar collector was carried out by Soin et al. [6]. The collector efficiency was found to increase linearly with the liquid level. The researchers investigated the thermal performance of a thermosyphon collector containing boiling acetone and petroleum ether, and developed a modified form of the Hottel- Whillier-Bliss equation (HWB)[11] to account for the fraction of liquid level in the collector with glazed. The steady-state thermal efficiency ηB containing a modified heat- removal factor FR, in the HWB equation, for a two-phase solar collector is firstly defined by Al-Tamimi and Clark [7], accounting for the boiling and the sub-cooled portions of the collector. Researchers conducted a detailed study of boiling fluid solar collector in 1984 [8], in which an analytical model was developed to investigate of sub cooling the fluid entering the collector and the level of fluid in the collector on the collector efficiency. The term Z* was defined as reference to the fraction of the collector required to heat the fluid to its boiling temperature. Mainly based on the results of Al-Tamini and Clark's research, an ASHRAE 109-A test standard for two-phase solar collectors was established in 1984. This standard defined five sets of conditions for two-phase collector thermal efficiency [12]. Ahmed et al. [13] developed and defined a new generalized heat removal factor, Fs and a new overall thermal loss coefficient, UL for the two-phase flat-plate collector. Further modification on UL was also done in 1989 [14]. First steady-state system simulations made with TRNSYS was done by Freeman et al. [15]. The first thermodynamic model to analyze two-phase solar collector is developed by Chaturvedi et al. [4]. The equilibrium homogeneous theory is used to model the two-phase flow in solar collectors. O’Dell et al. developed a design method for heat pumps with refrigerant filled solar collectors [16]. The researchers obtained the heat gain at condenser and the COP as a function of evaporation temperature. Numerical calculations of the collector efficiencies for a double-glazed two-phase flat-plate collector employing R-11 were carried out by Kishore et al. [17]. They firstly express the collector efficiency as a function of the saturation temperature and liquid level by combining their experimental data. Ramos et al.[18] also carried out theoretical investigation on two-phase collectors assuming laminar homogeneous flow and experimentally confirmed by them. Mathur et al. [19] developed a method calculating boiling heat transfer coefficient in a two-phase thermosyphon loop. The first detailed model for boiling solar collector using TRNSYS was developed by Price et al. [20] All the modeling or analysis mentioned above assume a homogeneous two-phase flow. The first theoretical model concerning non-homogenous for two-phase flow thermosyphon in the collector and condenser was carried out by Yilmaz [21]. His results show that the homogenous model is not sufficient to describe the two-phase flow in the collector. Variation of the properties of the working fluid and water with temperature are taken into account. Considering the effect of long-wave radiation and wind speed, a modified form of the Hottel-Whillier-Bliss solar collector efficiency function for unglazed solar collector characterization is investigated by Morrison and Gilliaert [22]. Using the approach of element analyze of the steady one-dimensional two-phase energy conservation equations in finite, a mathematical model for natural circulation two- phase solar collector with cover was developed by Radhwan et al. [23]. Torres Reyes et al. conducted the first exergy analysis on a solar-assisted heat pump with two-phase collector- IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 16 evaporator [24]. They theoretically analyzed and experimentally confirmed that 80% irreversibilities in SAHP system occurs in collector-evaporator because a significant fraction of the total solar radiation absorbed is poorly employed. Similar investigation and results were obtained by Cervantes [25]. Most of the researches before 1999 were for steady-state in two-phase flow. Hussein firstly investigated theoretically and experimentally a two-phase closed thermosyphon flat plate solar collector under transient conditions [26]. He improved the modified the model in 2002 [27]. However no one has investigated the effect of the condensation phenomenon caused by high relative humidity (RH) on the surface of un-glazed solar collectors. The effect of the condensation phenomenon on the collector performance was investigated in this study. A transient two-dimensional mathematical model of the evaporator-collector has been developed to predict temperature distribution and useful energy gain. A series of experiments were performed under the meteorological conditions of Singapore to validate the model. 2. SYSTEM DESCRIPTION A solar assisted heat pump system for air-conditioning, water heating and drying is shown in Fig. 1. This system consists of five major components, namely the solar collector, evaporator, condenser (water cooled condenser and air cooled condenser), expansion valve, and compressor. Refrigerant R-134a was used as the working fluid due to environmental and thermodynamic considerations [4]. In the present study, two evaporators connected in parallel can increase the system coefficient of performance (COP) significantly. One of the evaporators performs as a solar collector, absorbing heat from solar radiation and ambient air; while another evaporator absorbs heat from a room for cooling purposes. The energy from these two heat sources, plus the electrical energy supplied to the compressor, is used for water heating and air heating for drying purposes. Fig. 1: Schematic overview of solar assisted heat pump air-conditioner, water heater and dryer. The superheated refrigerant vapor leaving the compressor enters the water-cooled condenser and performs water heating. The refrigerant, a mixture of vapor and liquid, leaves the water-cooled condenser and enters an air-cooled condenser to ensure complete condensation of the refrigerant vapor. The hot air obtained from the air condenser is used for cloth drying purposes. Once complete condensation of the refrigerant vapor has taken place, the refrigerant path splits into two: one goes to evaporator-collector to absorb solar energy and the other goes to air-conditioned space to absorb thermal energy from the room IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 17 to provide cooling. Superheated refrigerant vapor from the two evaporators mix together and enters into the compressor for next cycle. Two evaporator-collectors each with an area of 1.5 m2, made of 1mm copper plates were connected in series, as shown in Fig. 2. A 9.52 mm diameter copper tube was soldered at the back of each absorber plate in a serpentine configuration (Figure 2a). Adequate insulations were provided at the back of the collector but no glass cover was used on the top surface, i.e. the unglazed collector (Figure 2b). The collector surface was coated with a black paint with absorptivity of 90%. There is a bypass line from the exit of the first evaporator-collector to the exit of the second evaporator-collector, which remains closed or open depending on the solar irradiation. (a) (b) Fig. 2: Photograph of the un-glazed solar evaporator-collector. 3. EVAPORATOR-COLLECTOR MODEL A transient two-dimensional mathematical model of the evaporator-collector has been developed to predict temperature distribution and useful energy gain. Fig 3: Geometry and coordinate system of the unglazed solar evaporator collector. The function of an evaporator collector is to deliver the heat absorbed on the collector surface to the fluid in the tube. Because the copper plate is thin and its onductivity is high, the temperature variation in the z-direction can be considered as negligible. However, the heat transfer in the y-direction, as shown in Fig. 3, cannot be considered negligible because of the fact that the fluid temperature is not constant along the length of the tube. For this reason, the governing equation is derived to predict the performance of the evaporator collector by considering heat transfer in both x and y directions. IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 18 Fig. 4: Energy balance in a control volume in the interior area of collector. An energy balance on the control volume is shown in Fig. 4 gives: ( )x y u x x y y p T Q Q Q Q Q c x y t             (1) The term ( )p T c x y t       represents the change of internal energy in the elemental control volume. The terms , and u x x y yQ Q Q   are given by:  u uQ q x y   ; 2 2x x x T Q Q k x y x        ; 2 2y y y T Q Q k x y y        Substituting these three terms into equation (1) gives the partial differential equation: 2 2 2 2 pu cq T T T k x y k t            (2) Where ( )u L L aq S Q S U T T     This partial differential equation (2) requires one initial condition and four boundary conditions for its solution. The initial condition is 0, for all values of and aAt t T T x y  3.1 The boundary at   0 and 0 2 W D y x     Because the boundary at 0y  is exposed to the ambient, so  1 0, y L y aAt y q U T T    x S xQ x xQ  y yQ y yQ  LQ IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 19 Where  1 0 and 0 2 W Dy y x T T      Fig. 5: Energy balance on the control volume at y = 0. At 0y  , an energy balance on the control volume shown in Fig. 5 gives: 1 2 ( ) 2 y x y u x x y p y Ty Q Q Q Q Q c x t             (3) The terms 2 y y Q   are given by: 1 2 y y y T Q k x y        Substituting into equation (3):   2 1 1 1 12 2 2y y p yu L y a T T c Tq U T T k x k y y y k t                (4) Similarly, the other three partial differential equations for boundaries are derived. 3.2 The boundary at   and 0 2 W D y L x       2 2 2 2 22 2 2y y p yu L y a T T c Tq U T T k x k y y y k t                3.3 The boundary at x = 0 and 0 < y < L 2 1 1 1 2 2 pu x x xcq T T T k y x x k t             3.4 The boundary at ( ) and 0 2 W D x y L      22 22 1 1 x f pu x x fi i b T T cq T T k D x k t k D h D C                  x S xQ x xQ  2 y yQ 2 y y Q   LQ IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 20 Where    2 1 2 and 0 0 and 0 and 0 2 2 ; ;W D W Dy x xy L x x y L x y LT T T T T T            The four corner points are special conditions of the boundary conditions. The PDE for them are derived by the same energy balance method as for the boundary conditions. 3.5 The Overall Heat Transfer Coefficient UL As shown in Fig. 6 the plate on the tube has five parts of heat gain. Fig. 6: The cross-section of plate and the tube. The energy released to the fluid is " , and '' '' '' '' '' ''[ ]u b r w m cq q q q q q     (5) Where, " is decided by the flow condition inside the tube. '' ( ) 1 1 b f u fi i b T T q h D C          Where              L hxJ D k h fgDi i f fi 1Re0082.0 2 The effect of condensation is built into the overall heat transfer coefficient, which can be separated into three parts: wind heat transfer coefficient ℎ , the radiation heat transfer coefficient ℎ and the condensation heat transfer coefficient ℎ l w r condU h h h   (6) IIUM Engineering Journal, Vol. 16, No. 2, 2015 Special Issue on Energy Amin et al. 21 Where: 2.8 3w ah V  (7)     4 4 p sky r p p a T T h T T     (8) An empirical equation has been developed by Swinbank in [8], which relates, Tsky to ambient temperature, Ta, as shown below: 4 1.50.0552sky aT T (9) where, they are both in Kelvin. For the condensation heat transfer coefficient ℎ , if Tp is greater than Tdew, condensation will not occur, then hcond = 0 If Tp