Modeling and Optimization of a Latent Heat Thermal Storage with Solar Absorption Chiller in Iran
Mohammad Mehdi Keshtkar*, Amin Sheibani
Department of Mechanical Engineering, Kerman Branch, Islamic Azad University, Kerman, Iran
To cite this article:
Mohammad Mehdi Keshtkar, Amin Sheibani. Modeling and Optimization of a Latent Heat Thermal Storage with Solar Absorption Chiller in Iran. Science Journal of Energy Engineering. Vol. 4, No. 5, 2016, pp. 44-49. doi: 10.11648/j.sjee.20160405.12
Received: October 20, 2016; Accepted: November 3, 2016; Published: November 25, 2016
Abstract: The research in this paper focuses on the theoretical modeling and optimization of a solar absorption chiller system with a latent heat thermal storage using phase change material (PCM). At first for modeling the absorption chiller cycle, we use the thermodynamic principles. After that with considering the technical limitations we change this model to a standard mathematical planning model by using exergy analyses. At the second step, the mathematical model will be optimized by a genetic algorithm which is inspired by nature. The optimized parameters of cycle will be calculated accurately by doing this step. These parameters consists of the temperature of essential parts i.e. generator, evaporator, absorber, and condenser. Designing a storage system is necessary for resolving the important problem of sun systems which is not absorbing the light of sun at night. So that in the final step, a cylindrical system for storage of energy sun by using storage technique of thermal energy and PCM is designed. This system is analyzed by the enthalpy and Finite Volume Method which their parameters are, dimensions and type of the PCM.
Keywords: Exergy, Genetic Algorithm, Optimizing, Storage System, Finite Volume Method
Nowadays, development and increase explosion of new and renewable energy systems is one of main needs of every society. Two reasons can be broached for this increasing of energies: 1) decreasing the amount and level of fossil fear’s production. 2) Today’s being discussed fossil fuels have more environmental results and new energies are excluded of this rule. One of these kinds of energy is solar energy of photovoltaics that currently have been used in small scales in different country. Main reason of run-used this energy on large industrial scales is high cost of structure against low efficiency. Researchers obliged to studying wide investigate of using this energy on large industrial scales in addition to small scales, because the freedom of input energy to solar systems.
Misra et al. (2002, 2003) have used first concept to optimize a LiBr-water absorption refrigeration system. They determined average cost per unit exergy of all internal flows. Kizilkan et al. (2007) used second concept to optimize a LiBr-water absorption refrigeration system thermoeconomically. Kim et al. (2011) used a micro genetic algorithm for optimizing the long-term performance of an active-indirect solar hot water plant in a flight kitchen at Changi international airport services. Kulkarni et al. (2007) proposed the "design space methodology" to design and optimization of an active-direct solar water heating system with one storage tank. A dynamic simulation tool specially devoted to the typical climates of South Italy for both winter heating and summer cooling load were presented by Calise et al. (2011) based on a rigorous mathematical algorithm. He also prepared a complete dynamic model of a solar-assisted heating and cooling system and used to develop a case study, including a detailed energy and economic analysis and optimization for different types of Italian school buildings. Prakash et al. (1985) applied a rectangular storage unit in water heater that contained a layer of PCM placed at the bottom. Sharma et al. (2003) developed latent heat storage unit for providing hot water in the evening and morning. Their used PCM was paraffin wax (m. p. 54°C) in a box type solar collector. Ghoneim and Klein (1989) investigated theoretically phase change heat storage for hot air and hot water. Sodium sulfate decahydrate and paraffin were used as PCMs. Buddhi and Sahoo, 1997 used Sytric Acid as a PCM in a box type solar cooker. Tan et al. (2009) investigated heat transfer during constrained melting of PCM inside a spherical capsule. A smooth circular tube with effective specific heat technique was developed for evaluating the heat transfer performance of a PCM by Kandasamy et al. (2008). Also, a report was presented about the variation of thermal resistance in PCM and showed the effect of natural convection will not be as significant (Tay et al., 2012a). Moreover, they modeled shell and tube heat exchanger with PCM using CFD method that was developed by ANSYS code (Tay et al., 2012b). Solar energy is usable of two aspects: 1) photo voltaic systems produce electric energy. 2) Thermal system that operated in heating and cooling fields of building. By focus to this subject that there isn’t any research in the field of cooling by absorption chiller along phase change material and optimization of these cycles, so this aims follow in the present work. Absorption system, operate heat of energies such as solar energy, geothermal energy and loss energy of industrial instruments for heating or cooling special environment like interior space of building. This system has high potential to save energy consumption and low-pollutant heat of environment.
By attention to the mention notes, aim of performing research and process of solving problem is as follow. In the first step, thermodynamic modeling and exeagy analyze done in single effect absorption chiller cycle with lithium-bromide and then principal parameters of cycle optimized by genetic algorithm. In the next step, type of collector, diameter and mass of phase change material are assigned. Finally, by focus on the time of using absorption chiller, temperature variation of PCM and HTF are being discussed on term of time.
2. Thermodynamic Modeling and Exergy Analyze
2.1. Energy Analysis
Schematic process of single effect absorption chiller cycle is shown in Fig. 1. Simple absorption system is composed from four main components. This system is composed from one evaporator, and absorber tower, on the low-pressure side of system and one generator and condenser on the high-pressure side of system. Working fluids are refrigerant (water) and absorbent (LiBr). The path circulation of refrigerant is from condenser to evaporator, tower, generator and the condenser, respectively.
Some assumptions are considered from the point of view of system modeling:
1) System operates in the permanent conditions.
2) Output solution of generator and absorber are considered at saturation temperature.
3) Output fluid of condenser is considered at saturation temperature.
4) Output steam of evaporator is considered at saturation temperature.
5) Valves operate at constant enthalpy process.
The coefficient of performance of cooling system (COPcool) can be defined for system as the ratio of heat extracted from evaporator to the total work of pump and input heat to generator and the coefficient of performance of heating system (COPheat) can be defined for system as the ratio of total extracted heat of absorber and condenser to total work of pump and input heat to generator that are defined as follow:
Equilibrium equation of energy for total cycle of absorption as follow:
2.2. Exergy Analysis
According to exergy analysis, aim is obtaining destruction exergy for all components of system (Ek) that k index is referred to each of system component separately: evaporator, condenser, absorber tower, heat exchanger of solution and pump. Exergy equation balance of system can be defined as follow:
By attention to above analysis, the effect of charging evaporator temperature, condenser, absorber and generator and also amount efficiency of heat exchange on COP and exergy efficiency of system are investigated (Fig. 2). After investigating the effect of each parameters, highest COP and exergy efficiency obtained by genetic algorithm. Then, operating point of system obtained from the point of view of COP and exergy efficiency.
2.3. Determination of Type of PCM, and Monthly Solar Radiation for Cooling System Along by Thermal Storage System
For solar absorption chiller effective and economic daily thermal energy storage is unavoidable, because solar energy is unavailable during certain periods of the day and during the night. Among the various methods of energy storage, the latent heat thermal energy system using PCM is quite appealing, because of its high energy storage density and its ability to provide heat at a constant temperature.
Main system (chiller) is assumed to be used in hours from 9 AM to 9 PM. This system operates from 9 AM to 5 PM by solar energy that is provided by the collector, but Support Energy System (SES) supports it during period 5 PM to 9 PM that is during the time that collector does not receive sufficient solar radiation. The PCM to be used in the design of thermal storage systems should possess desirable thermophysical, kinetic and chemical properties and to have a melting temperature near to working temperature of main system. Working temperature of main system in this study to be provided by collector and storage system is 88°C. So, the PCM used for this investigation is considered "Erythrol". The problem of phase change falls into the category of moving boundary problems. Solving this problem is done by enthalpy method. In this method the governing equation is divided into three states (liquid, mushy and solid) and each of state has specified governing equation. Up to now, after energy and exergy analyzing, here, by using thermal storage system with phase change material (PCM), need heat of generator provided by solar collector.
Location of Iran on north of 25-45 degrees latitude, is caused average solar energy about 19.23 MJ/m2 (5.3 kWh/m2. day). In this country amount of sunny hours is more than 2800 hours in year and according to government’s policy to develop produce of energy by new technique indicated the high potential of Iran to operating solar energy in different purpose. In present study, design of solar system is considered for 700m2 building area and positive space of ventilation is 600m2 in Tehran. Operation time of system is 16 hours, 10 hours relate to day and 6 hour relief to night. The period of operation time of system is to consider for 5 month of year between 15 Ordibehsht (May) to 15 Mehr (October). So, we refer to NASA data that is one of most valuable information resource to calculate solar radiation intensity and absorption.
According to this resource in July (Tir-Mordad) at Tehran, the average of solar radiation intensity is equal to, 7.5 kWh/m2 day, on horizontal surface. By assumption of about 9 hour radiation for this month of summer, so the amount of solar radiation monthly on horizontal surface is equal to G=834 w/m2. This power of solar-radiation of monthly is used to calculate of input thermal to collector.
By attention to gathering information of articles, we install collector on south side with angle of 36 degree which is about equal to Tehran’s latitude.
2.4. Selecting Collector Which Is Connected to Storage System
Main section and the heart of solar system is solar collector of system which transfers thermal energy to the working fluid inside the collector by absorbing solar radiation energy and exchange of it to heat.
It is noticeable that, the most suitable collector for solar absorption chiller (single-effect) is vacuum pipe collector type.
Need heat of generator is obtained by specified parameters in previous section by attention to synchronic optimization of COP and exergy efficiency.
By focus to below equation one can be calculate transferred heat by collector:
For vacuum pipe collector value of FRUL, is about 0.7 and is about 0.689.
Mass flow rate of nano fluid inside the collector is calculated by using optimum amount of Qgen and so we can assign need surface of solar collector.
The collector efficiency, is given by below equation:
Finally, by attention to the collector area and absorbent area of each vacuum pipe which is equal to 0.0804 m2, so, the number of vacuum pipe is calculated.
By focus to this note, that each 22 pipe compose one collector, so the number of need collector is calculated.
In this section, system calculations have been done by attention to previous condition along thermal storage system in two state of day and night.
Through the day, 40 percent of need thermal of generator and thermal energy to melting PCM is provided by solar energy and 60% remain of this thermal energy is provided by auxiliary heating system (AXHS) (Fig. 3a). Absorbed solar energy of collector (qcollector) is used to supply need thermal energy of generator and need thermal energy of melting PCM at during of the day. By assumption of 50 percent of mass flow rate of nano fluid is for the supply of 40 percent of need thermal energy of generator, so, amount of mass flow rate will be equal to 2.7 kg/s. 50 percent of remain mass flow rate of HTF will be arrived into the thermal storage system. HTF in thermal storage section provides need heat for melting of solid PCM. Solid PCM is liquefied and then prepared to giveback heat and supply need thermal of generator during the night.
During the night, elective valves SV2, SV1 are closed (Fig. 3b). At this case, input temperature of HTF is equal to 30°C, that this temperature increases by PCM until 40°C and remain thermal energy is provided by auxiliary heating system. Now, by attention to this subject that the mass flow rate of HTF provides need thermal energy for melting of PCM, so, accumulated thermal energy inside the PCM released and the temperature of HTF increases from 30°C to 40°C during the night. We can assign time of discharging of PCM by considering of charging time and mass flow rate of HTF during the night. Accumulated energy inside the PCM in the period of charging time is equal to removal heat of PCM in the period of discharging time and absorbing heat by Nano fluid. According to Fig. 4 time of attaining HTF to 40°C is equal to 100 minutes after operating thermal discharge. It should be mentioned that in the time of discharging, accumulated thermal energy inside the PCM transmits to nano fluid material. According to this subject, need thermal of HTF in this system will be 40°C until about 5 hours. Need thermal energy to bring HTF temperature from 40°C to 88°C which is the temperature of operating single-effect absorption cooling system is provided by auxiliary heating system.
By attention to optimization of cooling cycle in previous section, amount need thermal of generator in optimum state calculated. According to the optimization results by genetic algorithm, maximum of COP will be equal to 0.778, temperature of evaporator, condenser, generator and absorber in COP mode is 10, 30, 80 and 30°C, respectively. Maximum exergy efficiency is 0.423 and temperature of evaporator, condenser, generator and absorber in exergy efficiency mode is 10, 30, 80, 30°C, respectively.
By synchronic optimization of exergy efficiency and COP, these parameters will be obtained 0.417 and 0.778, respectively. Temperature of evaporator, condenser, generator and absorber in multi object optimization are 8.4, 30.2, 80.2 and 30.33°C, respectively. Table 1 shows range of COP difference according to various temperature range for components at TE=10°C
|Parameter||Temp Range||COP Diff|
|Absorber (°C)||30< T2<40||0.032|
|Generator (°C)||80< T5<95||0.015|
|Condenser (°C)||30< T9<40||0.014|
Table 2 shows range of exergy efficiency difference according to different temperatures for components at TE=10°C
|Parameter||Temp Range||exergy efficiency|
|Absorber (°C)||30< T2<40||0.016|
|Generator (°C)||80< T5<95||0.008|
|Condenser (°C)||30< T9<40||0.01|
Finally, by attention to modeling and optimization of absorption cooling system of liBr-water, concluded that:
1) By increasing evaporator temperature, COP and exergy efficiency will be increase.
2) By increasing condenser temperature COP and exergy efficiency will be decrease.
3) By increasing generator temperature, COP and exergy efficiency will be increase and then decrease.
4) By increasing absorber temperature, COP and exergy efficiency will be decrease.
According to positive area of 600m2, cooling capacity is 84.48 kW and generator power is 118.45kW. So 47.38 kW of this thermal is provided by solar collector and 71.07kW of it, is provided by AXHS. By synchronic optimization of exergy efficiency and COP, these parameters will be obtained 0.417 and 0.778, respectively. Temperature of evaporator, condenser, generator and absorber in multi object optimization are 8.4, 30.2, 80.2 and 30.33°C, respectively. Material of HTF provides need thermal energy to melting 62 kg of PCM (Erythrol) which is inside 50 storage pipes and melting temperature of it is equal to 135°C. Total transfer thermal by HTF is equal to 132.51 kW. By attention to this subject that output power to collector area is equal to 541.35 W/m2, so collector need area will be equal to 242.76m2 and then we will need 3014 discharge pipes. Each 22 pipe constitute one separate collector and it is clearly that we will need 138 collectors. Also, amount efficiency of it, is equal to 64.9%. By focus to above chart, it is clearly that. Time of charging is assumed 10 hours and time of attaining HTF to 40°C is equal to 100 minutes. According to this subject, need thermal of HTF in this system will be 40°C until about 5 hours.
P: Solution Pump
This work supported by Research Program supported by the Department of Mechanical Engineering Azad university of kerman, Iran.