**World Journal of Engineering and Technology**

Vol.07 No.01(2019), Article ID:90501,17 pages

10.4236/wjet.2019.71011

Exergy Analysis of a Double-Effect Solar Absorption Refrigeration System in Ngaoundere

Maurice Tenkeng^{1}^{,2}, Paiguy Armand Ngouateu Wouagfack^{3*}, Réné Tchinda^{2,3 }

^{1}Department of Physics, University of Dschang, Dschang, Cameroon

^{2}LISIE, University Institute of Technology Fotso Victor, University of Dschang, Dschang, Cameroon

^{3}Department of Renewable Energy, Higher Technical Teachers’ Training College, University of Buea, Kumba, Cameroon

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: September 2, 2018; Accepted: February 11, 2019; Published: February 14, 2019

ABSTRACT

Solar energy is replacing more and more traditional sources of energy because of the fact that it’s also fighting about global warming. This study is based on exergy analysis of a double-effect series flow absorption refrigeration system powered by solar energy in Ngaoundere. The simulation is done on the basis of a half hourly analysis for the first time, from 6.30 AM to 6.30 PM, using water-lithium bromide as working pair. The main parameters for the performance of an absorption cycle, which are the COP and the ECOP, have been analyzed and the results show that this two parameters increase while increasing the temperature of the main generator. The exergy loss of each component of the system and the total exergy loss of the system have been analyzed and their effectiveness calculated, using the first and second law of thermodynamics. The highest exergy loss occurs in the main generator GI and in the absorber, making these components more important in an absorption cycle. This analysis is based on a mathematical model using FORTRAN language. The results obtained may be useful for the optimization of solar absorption refrigeration systems.

**Keywords:**

Double Effect, Refrigeration, Half Hourly, Exergy Loss

1. Introduction

Refrigeration is one of the active cooling methods studied here to enable numerous villages from remote areas and why not towns, to benefit from welfare offered by nature; areas where electricity is in short supply are also concerned. Refresh air or refrigeration is becoming more and more important in the daily activities nowadays. The refrigerators are useful for making ice, for storing vaccines and food in areas where electricity is unavailable or high in cost that is the case in Africa in general and Cameroon in particular. Passive cooling systems, that is: good insulation, double glazing or use of thermal mass and ventilation are no longer sufficient. In inter-tropical regions, with the sun available all the yearlong, solar energy has a considerable interest because of the promising advantages. The First solar refrigeration system was introduced by Ferdinand Carré in 1895 [1] . Double-effect absorption refrigeration cycle was introduced for the first time in 1956 for developing the system performance within a heat source at high temperature [2] .

In order to improve the performance of a solar absorption refrigeration system, number optimization studies based on the energy analysis have been done. Energy analysis takes into account the first law of thermodynamics which deals with the conversion of energy; this cannot show where the irreversibility occurs in a system that has many sources of energy [3] .

The exergy analysis based on the second law of thermodynamics, is the only way to detect irreversibility in different components of the system that is the only rational basis for evaluating the system performance [4] .

Many researchers have carry out studies on exergy analysis of double-effect refrigeration systems. K. Talukdar et al. [5] performed an exergy analysis to show the performance variation of both the topping PC and the bottoming ARS with changing HPG temperature from 120˚C to 150˚C. S.A.M Said et al. [6] brought two modifications to the traditional solar absorption refrigeration system by recovering waste heat from dephlegmator and utilization of a refrigerant storage unit. The COP increases 18% because of these modifications. The exergy variation of a solar assisted double-effect absorption refrigeration system, was studied in detail and the parameters of the cycle analyzed by Ravikumar T.S et al. [7] on the basis of first and second law effectiveness and the results indicated that a better quality of the evaporator, has more effect on the system performance than the better quality of other components. The study of Nahla Bouazis et al. [8] reveals that, the performance of a novel configuration is better than that of the two stage conventional configuration and it allows a lower temperature, about 60˚C - 120˚C instead of 100˚C - 160˚C for conventional cycle. T. Avanessian et al. [9] analysed and compared different water-cooled LiBr-H_{2}O absorption systems under different operating and climatic conditions, and the effect of considering the chemical exergy of the LiBr-H_{2}O solution on the second law analysis of such a system is probed. Milad Ashouri et al. [10] examined the performance of a double pressure organic Rankine cycle couple with a solar collector, via a thermal tank during a day. Results showed that the system is capable of generating stable power during the day with a solar fraction of 100%. Kai Zhang et al. [11] discussed in detail the exergy loss for each component, in order to identify the inefficient components and then to indicate the way to improve the performance of a double-effect evaporating temperature chiller. Rabah Gomri [12] investigated the potential of application of a single effect and multiple effect absorption cooling systems for the production of chilled water. The simulation results were used to study the influence of the varios operating parameters on the performance of single, double and triple effect absorption cooling cycles. Omer Kaynakli et al. [13] investigated the equation temperature of the overall system and its components by conducting a parametric study on heat capacity and exergy destruction of the HPG, coefficient of performance (COP) of the system and mass flow rate of heat sources. R. Maryami et al. [14] done a comparative study between the five classes of LiBr/water absorption refrigeration systems to investigate the influence of various operating parameters on the coefficient of performance. Zeyu Li et al. [15] give the yearly performance analysis of a solar air cooled double effect LiBr/H_{2}O absorption cooling system in subtropical city, including total efficiency, cooling capacity per area of collector and solar fraction. Ahmed Hamed et al. [16] optimized the coefficient of performance of a double-effect absorption system and predicted the coefficient of performance and exergetic efficiency of the system with two equations developed. Rabah Gomri et al. [17] simulated a double-effect vapor absorption cooler system to calculate exergy losses, the coefficient of performance and the number of exergy of each component. A. Allouhi et al. [1] presented a review of the available Technologies to provide cooling from solar energy for both thermal and photovoltaic ways. A market study is also carried out to evidence the potentiel of these solar cooling technologies with overview of the various solar cooling installations in Europe, Egypt and China was also presented. S.W Sharshir et al. [18] highlighted design methods to allow the researchers to optimize the solar still further development and proposed some recommendations.

In this study, in order to analyse the system with a great accuracy, we decided to look more closely by the simulation of a half hour analysis, from 6:30 AM to 6:30 PM of a double effect series flow absorption refrigeration system, using water-lithium bromide solution as working fluid with an Evacuated tube solar collector type. The data, which are solar insolation and environmental temperature are taken from M. Tenkeng et al. [19] and were collected on the 15th of January 2014, in Ngaoundere, Cameroon. The main focus of this study is concentrated on the exergy analysis of each component of the system, precisely the exergy loss [20] . The coefficient of performance (COP) and the exergetic coefficient of performance (ECOP) of the system are also investigated.

2. Description of the Model and Assumptions

2.1. Description of a Series Flow Solar Absorption Cycle

Double-effect consists of components similar as in single-effect absorption system with an additional generator, see Figure 1. The weak solution at state 1 is pumped from the absorber to the first-effect solar assisted generator (GI) through two heat exchangers (HEX I and HEX II). The LiBr-H_{2}O solution in the first-effect generator is heated at high temperature. This high temperature heat comes from a solar collector to supplies to the first-effect generator to boil out the water vapour from the solution. The vapour refrigerant generated is condensed at high pressure in the second-effect generator. Thus, the primary vapour, from the first-effect generator goes to the secondary condenser incorporated in the second-effect generator (GII). The heat of condensation from this condenser is utilized by the second-effect generator. The strong solution, so called because it’s poor in water, which is released by the first-effect generator at state 5 flows to the second-effect generator, through HEX II where the solution is cooled to some extent by exchanging heat with the weak solution. The heat rejected is used to produce addition refrigerant vapour from the solution coming from the first-effect generator. At the outlet of the second-effect generator vapour is produced, which together with the condensed water vapour from the secondary condenser, flows into the main condenser through a throttle valve where the heat of condensation is rejected to the atmosphere. Thus, the overall amount of liquid refrigerant leaving the main condenser is the sum of refrigerant coming from the first and second-effect generators. The refrigerant in a liquid phase from this condenser, on expansion, flows to the evaporator through a throttle valve (REV) where it is evaporated at low pressure. After extracting heat from the medium to be cooled, the refrigerant evaporates and then passes to the absorber where it gets absorbed by the strong solution coming from the second-effect generator, through HEX I and a throttle valve (SEV). The resulting weak solution in the absorber is then pumped to the first-effect generator and the cycle completes.

Figure 2 shows the P-T-X diagram of our chiller. Idealy, this system has three pressure levels and five operating temperature levels which are: T_{e}; T_{a}; T_{c}; T_{gI} and T_{gII}. Low pressure prevails in the evaporator and absorber (
${P}_{1}={P}_{a}={P}_{e}$ LPL), high pressure in the first-effect generator (GI) (
${P}_{3}={P}_{GI}$ HPL) and intermediate pressure in the condenser and the second-effect generator (GII) (
${P}_{2}={P}_{c}={P}_{GII}$ IPL).

2.2. Assumptions

For the simulation of the cycle, following assumptions were made [21] [22] [23] :

1) The ambient temperature is assumed to be 298.15 K;

2) The refrigerant completely evaporates in the evaporator;

3) The thermal and pressure losses can be ignored;

4) The solution is saturated after is passes through the absorber and the generator;

5) The refrigerant vapor is completely condensed in the condenser;

6) For the temperature range considered, the absorbent is non-volatile;

7) The refrigerant leaves the condenser and evaporator at saturated state;

8) Changes occur only at pump and pressure-reducing valves for the pressure;

9) Because thermodynamic work required change in specific volume liquid, the solution pump work is neglected;

10) In the valves, the expansion process is isenthalpic.

3. Thermodynamic Analysis

3.1. Exergy Analysis of the System

The exergy of fluid stream can be defined as [24] :

$\phi =\left(h-{h}_{0}\right)-{T}_{0}\left(s-{s}_{0}\right)$ (1)

where $\phi $ is the exergy of the fluid at temperature T, h and s are respectively enthalpy and entropy of the fluid, ${h}_{0}$ and ${s}_{0}$ are respectively the enthalpy and entropy of the fluid at environmental temperature ${T}_{0}$ (298.15 K).

The specific exergy content of a mixture with m components is defined by A. Sözen et al. 2007 [25] as:

$\phi ={\displaystyle \underset{n=1}{\overset{m}{\sum}}{x}_{n}{h}_{n}-{T}_{0}{\displaystyle \underset{n=1}{\overset{m}{\sum}}{x}_{n}{s}_{n}-{\displaystyle \underset{n=1}{\overset{m}{\sum}}{x}_{n}{\mu}_{n0}}}}$ (2)

where
${T}_{0}$ is the reference temperature (25˚C),
${x}_{n}$ is mass fraction and
${\mu}_{n0}$ the chemical potential of the n^{th} component of the mixture at
${T}_{0}$ and
${P}_{0}$ .

The exergy loss or the availability loss in each component is given by [21] :

$\Delta \phi ={\displaystyle \sum {m}_{i}{\phi}_{i}}-{\displaystyle \sum {m}_{0}{\phi}_{0}}-\left[{\displaystyle \sum Q{\left(1-\frac{{T}_{0}}{T}\right)}_{i}}-{{\displaystyle \sum Q\left(1-\frac{{T}_{0}}{T}\right)}}_{ref}\right]+{\displaystyle \sum W}$ (3)

where $\Delta \phi $ is the lost exergy or irreversibility that occurs in the process. The first two terms of the right hand side are the exergy of the inlet and outlet streams of the control volume. The third and fourth terms are the exergy associated with the heat transferred from the source maintained at a temperature T. The last term is the exergy of mechanical work added to the control volume. The term is negligible for absorption systems as the solution pump has very low power requirements.

The total exergy loss of absorption system is the sum of exergy loss in each component and is written as:

$\Delta \phi ={\displaystyle \underset{n=1}{\overset{m}{\sum}}\Delta {\phi}_{n}}$ (4)

De Vos (1992) [26] established equations for calculating the exergy loss to the collector per unit area as function of the exergy emitted from the sun minus the albedo of the earth and the radiation emitted from the solar collector. The exergy analysis of the collector can be obtained with Pridasa was et al. 2004 [27] :

Exergy (radiation) input:

$\phi =f\sigma {T}_{sun}^{4}+\left(1-f\right)\sigma {T}_{p}^{4}-\sigma {T}_{sc}^{4}$ (5)

where sun is for sun, p for planet and sc for solar collector, f is the sunlight dilution factor equal to $2.16\times {10}^{-5}$ on earth. The solar radiation that reaches the solar collector is transformed into heat. This heat is partly absorbed by thermal fluid and the surrounding equipment and partly lost to the environment. The available solar radiation is transformed into available heat for the process, however, the second law of thermodynamics hinders the transformation of all heat into exergy [26] [27] . The exergy of the solar heat input to the solar collector is given by the relationship:

Exergy (heat) input:

${\phi}_{s,h}={Q}_{ava}\left(1-\frac{{T}_{ref}}{{T}_{sc}}\right)$ (6)

Exergy loss during the transformation process for each component is as follows:

${I}_{sc,r}={\phi}_{s}-{\phi}_{s,h}$ (7)

The useful exergy gained by the solar collector is:

${\phi}_{su}={Q}_{u}\left(1-\frac{{T}_{ref}}{{T}_{sc}}\right)$ (8)

where Q_{u} is the useful steady state energy gain to the solar collector. The exergy loss from the input of solar collector to the working fluid can be calculated using the equation:

$\Delta {\phi}_{sc}={\phi}_{s}-{\phi}_{su}-\Delta {\phi}_{en}-\Delta {\phi}_{sc,r}$ (9)

where Δψ_{sc}_{,r} is the exergy lost during the transformation from solar energy radiation to heat on the solar collector and written as [26] :

$\Delta {\phi}_{sc,r}={\phi}_{s}-{\phi}_{s,h}$ (10)

And Δψ_{en} is the average loss caused by the wind effect on the solar collector. It’s impossible to calculate a priori, since the speed, the direction and the sense of the wind are unknown [27] . Moreover, as the plate temperature is a function of Δψ_{en}, Δψ_{en} cannot be determined either. For these reasons, it will be assume for simplicity Δψ_{en} = 0 and therefore:

$\Delta {\phi}_{sc}={\phi}_{s}-{\phi}_{su}-\Delta {\phi}_{sc,r}$ (11)

The exergy efficiency for cooling is the ratio of the chilled water exergy at the evaporator to the exergy of the heat source at the generator:

${\phi}_{cooling}=\frac{{m}_{22}\left({\phi}_{22}-{\phi}_{23}\right)}{{m}_{18}\left({\phi}_{18}-{\phi}_{19}\right)}$ (12)

For the thermodynamic analysis, we have equations adapted from O. Kaynakli et al. (2007) [28] :

The mass balance:

${\stackrel{\dot{}}{m}}_{w}={\stackrel{\dot{}}{m}}_{s}+{\stackrel{\dot{}}{m}}_{{H}_{2}O}$ (13)

${\stackrel{\dot{}}{m}}_{w}{X}_{w}={\stackrel{\dot{}}{m}}_{s}{X}_{s}$ (14)

where X is the LiBr concentration, w (weak solution) and s (strong solution).

The flow rate of the strong and weak solutions can be determined from equation below:

${\stackrel{\dot{}}{m}}_{s}=\frac{{X}_{w}}{{X}_{s}-{X}_{w}}{\stackrel{\dot{}}{m}}_{{H}_{2}O}$ (15)

${\stackrel{\dot{}}{m}}_{w}=\frac{{X}_{s}}{{X}_{s}-{X}_{w}}{\stackrel{\dot{}}{m}}_{{H}_{2}O}$ (16)

The double circulation ratios have been defined in the cycle as:

${f}_{1}=\frac{Xw}{{X}_{s}{}_{1}-{X}_{w}}=\frac{{X}_{1}}{{X}_{5}-{X}_{1}}$ and ${f}_{2}=\frac{{X}_{s1}}{{X}_{s2}-{X}_{s1}}=\frac{{X}_{5}}{{X}_{14}-{X}_{5}}$ (17)

Using circulation ratios, the heat capacity of each component is:

For the generators

${Q}_{GI}={h}_{10}+{f}_{1}{h}_{5}-\left({f}_{1}+1\right){h}_{17}$ (18)

${Q}_{GII}=\left({f}_{2}\left({f}_{1}+1\right)\right){h}_{13}-{f}_{2}{h}_{9}+{f}_{1}\left({f}_{2}/\left({f}_{1}+1\right)\right){h}_{14}$ (19)

For the condenser

${Q}_{C}=\left({f}_{2}/\left({f}_{1}+1\right)\right){h}_{13}-\left({f}_{2}/\left({f}_{1}+1\right)+1\right){h}_{15}$ (20)

For the evaporator

${Q}_{E}=\left(\left({f}_{2}/\left({f}_{1}+1\right)\right)+1\right)\left({h}_{17}-{h}_{16}\right)$ (21)

For the absorber

${Q}_{A}=\left({f}_{2}/\left({f}_{1}+1\right)\right){h}_{17}+\left(\left({f}_{21}{f}_{1}\right)/\left({f}_{1}+1\right)\right){h}_{7}-\left({f}_{2}+1\right){h}_{1}$ (22)

For the heat exchangers

${Q}_{HEXI}=\left({f}_{2}+1\right)\left({h}_{3}-{h}_{2}\right)$ (23)

${Q}_{HEX\text{}II}=\left({f}_{2}+1\right)\left({h}_{4}-{h}_{3}\right)$ (24)

The pump work can be expressed as:

${W}_{p}={\nu}_{ws}\left({P}_{GII}-{P}_{E}\right)\left({f}_{2}+1\right)/{\eta}_{p}$ (25)

The coefficient of performance for the system is given by:

For the exergy destruction, the following equations are determined as:

For the generators

${I}_{GI}=\stackrel{;}{{m}_{10}}\left({\phi}_{18}\left({f}_{2}+1\right)-{\phi}_{5}{f}_{2}-{\phi}_{10}\right)+\stackrel{.}{{m}_{18}}\left({\phi}_{18}-{\phi}_{19}\right)$ (26)

${I}_{GII}=\stackrel{;}{{m}_{10}}\left({\phi}_{9}{f}_{2}-{\varphi}_{13}\left({f}_{2}/\left({f}_{1}+1\right)\right)-{\phi}_{14}\left({f}_{1}\left({f}_{2}/\left({f}_{1}+1\right)\right)\right)\right)+\stackrel{.}{{m}_{10}}\left({\phi}_{10}-{\phi}_{11}\right)$ (27)

For the condenser

${I}_{C}=\stackrel{;}{{m}_{10}}\left({\phi}_{13}\left({f}_{2}/\left({f}_{1}+1\right)\right)+{\phi}_{12}-{\phi}_{15}\left({f}_{2}/\left({f}_{1}+1\right)\right)\right)+\stackrel{.}{{m}_{20}}\left({\phi}_{20}-{\phi}_{21}\right)$ (28)

For the evaporator

${I}_{E}=\stackrel{;}{{m}_{10}}\left(\left({f}_{2}/\left({f}_{1}+1\right)\right)+1\right)\left({\phi}_{16}-{\phi}_{17}\right)+\stackrel{.}{{m}_{22}}\left({\phi}_{22}-{\phi}_{23}\right)$ 29)

For the absorber

${I}_{A}=\stackrel{;}{{m}_{10}}\left(\left({\phi}_{17}\left({f}_{2}/\left({f}_{1}+1\right)\right)+1\right)+{\phi}_{7}\left(\left({f}_{2}{f}_{1}\right)/\left({f}_{1}+1\right)\right)-{\phi}_{1}\left({f}_{2}+1\right)\right)+\stackrel{.}{{m}_{24}}\left({\phi}_{24}-{\phi}_{25}\right)$ (30)

For the throttle valves

${I}_{REV}=\stackrel{;}{{m}_{10}}\left(\left({f}_{2}/\left({f}_{1}+1\right)\right)+1\right)\left({\phi}_{15}-{\phi}_{16}\right)$ (31)

${I}_{SEV}=\stackrel{;}{{m}_{10}}{f}_{2}\left({\phi}_{6}-{\phi}_{7}\right)$ (32)

For the heat exchangers

${I}_{HEXI}=\stackrel{;}{{m}_{10}}\left(\left({f}_{1}+1\right)\right)\left({\phi}_{2}-{\phi}_{3}\right)+\left({f}_{1}{f}_{2}\right)/\left({f}_{2}+1\right)\left({\phi}_{6}-{\phi}_{14}\right)$ (33)

${I}_{HEXII}=\stackrel{;}{{m}_{10}}\left({f}_{2}\left({\phi}_{5}-{\phi}_{8}\right)+\left({f}_{2}+1\right)\left({\phi}_{3}-{\phi}_{4}\right)\right)$ (34)

For the pump

${I}_{PUMP}=\stackrel{;}{{m}_{10}}\left({f}_{1}+1\right)\left({\phi}_{1}-{\phi}_{2}\right)$ (35)

The total exergy destruction for all the system is:

${I}_{S}={I}_{GI}+{I}_{GII}+{I}_{C}+{I}_{A}+{I}_{E}+{I}_{REV}+{I}_{SEV}+{I}_{HEXI}+{I}_{HEXII}+{I}_{PUMP}$ (36)

A. Sözen et al. 2007 [27] defined the exergetic coefficient of performance on which depends the effectiveness of an ARS:

$ECOP=\frac{{\stackrel{\dot{}}{q}}_{e}\left(1-\frac{{T}_{0}}{{T}_{e}}\right)}{{\stackrel{\dot{}}{q}}_{g}\left(1-\frac{{T}_{0}}{{T}_{g}}\right)+{\stackrel{\dot{}}{W}}_{pe}}=\frac{{m}_{17}\left({\phi}_{17}-{\phi}_{16}\right)}{{m}_{4}{\phi}_{4}+{m}_{10}{\phi}_{10}-{m}_{5}{\phi}_{5}}$ (37)

And the coefficient of performance of the system for cooling purpose is:

$COP=\frac{{\stackrel{\dot{}}{q}}_{e}}{{\stackrel{\dot{}}{q}}_{g}}=\frac{{m}_{17}\left({h}_{17}-{h}_{16}\right)}{{m}_{4}{h}_{4}+{m}_{10}{h}_{10}-{m}_{5}{h}_{5}}$ (38)

3.2. Input Parameters for the Simulation

To simulate the above equations, a simulation program written in FORTRAN was developed. Properties of the LiBr-H_{2}O solution and typical overall heat transfer coefficient were adopted from Gur Mittelman 2007 [29] . Also, the solar insolation and the ambient temperature, the solution mass flow rate and the solution concentration are given in Table 1.

The cooling load of the system is Q_{E} = 10 kw and the evaporator temperature is T_{E} = 6˚C, the reference temperature is T_{0} = 25˚C. The temperature of the cooling water in the absorber and the condenser is 3˚C. The solar collector is a flate plate collector with 60˚ angle of incidence and 100 m^{2} of area. Properties of the LiBr-H_{2}O solution and typical overall heat transfer coefficient were adopted from Gur Mittelman [29] . For the heat transfer area or total conductance, we adopted the values from M. Tenkeng et al. [19] : (UA)a = 1.800 kW/K, (UA)e = 2.250 kW/K, (UA)g = 1.000 kW/K, (UA)c = 1.200 kW/K.

Figure 1. Doubleeffectserie flow solarabsorptionrefrigerartionchiller.

Figure 2. P-T-X diagram of double effect series flow absorptionregfrigerationcycle.

4. Results and Discussion

To have a better understanding of the results and their interpretations, thermodynamics properties of each point of the entired system are given in Tables 1-4.

Table 1. Ambient temperature and Solar insolation in Ngaoundere on the 15th of January 2014 [19] .

Table 2. The four level temperature at 12.30 PM.

Table 1 gives the solar insolation and the ambient temperature collected in the Meteological Station of Ngaoundere on the 15th of January 2014 by M. Tenkeng

Table 3. Thermodynamic properties of the cycle at 12.30 PM.

Table 4. The effectiveness for each component at 12.30 PM.

et al. [19] . Table 2 presents the temperature of the components of the system at 12.30 PM. Table 3 shows the thermodynamic properties of each of the points within the cycle, enthalpy, entropy and exergy at 12:30 PM can be slightly compared to that obtained by Omer Kaynakli et al. [13] for the different points of the system. The values of the exergy for the two studies are similar. Table 4 shows the effectiveness of the two generators, the absorber, the evaporator and the condenser. It is noted that the first effect generator GI is a component with a low effectiveness. This indicates that the generator needs a particular attention and must be optimized. Figure 1 shows the variation of insolation during the day chosen for our study. This curve has the same behavior with that of the ambient temperature shown in Figure 3. The variation of the exergy loss of the whole system with the temperature of the first effect generator GI is presented in Figure 4. While the temperature of GI increases, the exergy loss decreases. The lowest exergy loss is obtained for a generator temperature of 150˚C, the maximum temperature for the second effect generator powered by solar energy [13] . Figure 5 shows the variation of the ambient temperature with the maximum value at 12.30 PM. The temperature of the solar collector has the same behavior. Figure 6 and Figure 7 show the variation of the COP and ECOP respectively with the temperature of the first effect generator GI. The COP and ECOP depend on the generator load Qge and the evaporator load Qev. The COP decreases almost nonlinearly while the temperature of GI increases. This is a general

Figure 3. Variation of solar insolation during the day.

Figure 4. Variation of the exergy loss of GI with the temperature of GI.

Figure 5. Variation of the ambient temperature during the day.

Figure 6. Variation of the COP with the temperature of GI.

Figure 7. Variation of ECOP with the GI temperature.

behavior of Absorption Refrigeration Systems between the COP and the temperature of the generator receiving the main energy (GI). The ECOP is a very

Figure 8. Variation of the COP with the temperature of GII.

important tool for the evaluation of the performance of a system while applying the second law of thermodynamics. The ECOP decrease when the temperature of the first effect generator increase with a minimum value of ECOP which is 0.8, result in accordance with the litterature. Figure 8 shows the variation of the COP with the temperature of GII, the COP increases when the temperature increases. This result is in adequation with the need of the second generator in the double effect series flow absorption refrigeration system. So, the main impact of the second generator is to increase the value of the COP of the system.

5. Conclusions

By applying the first and the second laws of thermodynamics, the exergy analysis of a double-effect series flow solar absorption refrigeration system, with the most commun working pair, that is LiBr-H_{2}O, has been performed under weather conditions of Ngaoundere, in Cameroon, of the 15th January 2014. The study was made throughout the middle of the dry season. The exergy loss of the absorption cooling device and half hourly exergy destruction values of components were determined. The COP and the ECOP of the system were also investigated. The main conclusions obtained from the present study are as follows:

1) The maximum exergy destruction occuring in the solar energy assistance system is in the generator and the evaporator.

2) The exergy loss in the generator indicated that the generator is the component that needs to optimise.

3) The COP and ECOP variation with the temperature of the generators brings us to the conclusion that the second generator has a great contribution to the double effect refrigeration system. The system can therefore be implemented in Cameroon and in countries with high solar availability.

Our study can serve as preliminary investigation to solar water with cooled double effect absorption refrigeration system in tropical cities, and it is necessary to provide a foundation for further researches.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Tenkeng, M., Wouagfack, P.A.N. and Tchinda, R. (2019) Exergy Analysis of a Double-Effect Solar Absorption Refrigeration System in Ngaoundere. World Journal of Engineering and Technology, 7, 158-174. https://doi.org/10.4236/wjet.2019.71011

References

- 1. Allouhi, A., Kouskso, T., Jamil, A., Bruel, P., Mourad, Y. and Zeraouli, Y. (2015) Solar Driven Cooling Systems: An Updated Review. Renewable and Sustainable Energy Reviews, 44, 159-181. https://doi.org/10.1016/j.rser.2014.12.014
- 2. Ullah, K.R., Saidur, R., Ping, H.W., Akikur, R.K. and Shuvo, N.H. (2013) A Review of Solar Thermal Refrigeration and Cooling Methods. Renewable and Sustainable Energy Reviews, 24, 499-513. https://doi.org/10.1016/j.rser.2013.03.024
- 3. Sozen, A. and Yücesu, H.S. (2007) Performance Improvement of Absorption Refrigeration System Using Triple-Pressure-Level. Renewable Energy, 32, 267-284.
- 4. Lior, N. and Zhang, N. (2007) Energy, Exergy, and Second Law Performance Criteria. Energy, 32, 281-296. https://doi.org/10.1016/j.energy.2006.01.019
- 5. Talukdar, K. and Gogoi, T.K. (2016) Exergy Analysis of a Combined Vapor Power Cycle and Boiler Flue Gas Driven Double Effect Water-LiBr Absorption Refrigeration System. Energy Conversion and Management, 108, 468-477. https://doi.org/10.1016/j.enconman.2015.11.020
- 6. Said, S.A.M., El-Shaarawi, M.A.I. and Siddiqui, M.U. (2015) Analysis of a Solar Powered Absorption System. Energy Conversion and Management, 97, 243-252. https://doi.org/10.1016/j.enconman.2015.03.046
- 7. Ravikumar, T.S., Suganthi, L. and Samuel, A.A. (1998) Exergy Analysis of Solar Assisted Double Effect Absorption Refrigeration System. Renewable Energy, 14, 55-59. https://doi.org/10.1016/S0960-1481(98)00047-0
- 8. Bouaziz, N. and Lounissi, D. (2015) Energy and Exergy Investigation of a Novel Double Effect Hybrid Absorption Refrigeration System for Solar Cooling. International Journal of Hydrogen Energy, 40, 13849-13856. https://doi.org/10.1016/j.ijhydene.2015.05.066
- 9. Avanessian, T. and Ameri, M. (2014) Energy, Exergy, and Economic Analysis of Single and Double Effect LiBr-H2O Absorption Chillers. Energy and Buildings, 73, 26-36. https://doi.org/10.1016/j.enbuild.2014.01.013
- 10. Ashouri, M., Ahmadi, M.H., Pourkiaei, S.M., Astaraei, F.R., Ghasempour, R., Ming, T. and Hemati, J.H. (2017) Exergy and Exergo-Economic Analysis and Optimization of a Solar Double Pressure Organic Rankine Cycle. Thermal Science and Engineering Progress, 6, 72-86.
- 11. Zhang, K., Zhu, Y.T., Liu, J.X., Niu, X.F. and Yuan, X.L. (2017) Exergy and Energy Analysis of a Double Evaporating Temperature Chiller. Energy & Buildings, 165, 464-471.
- 12. Gomri, R. (2010) Investigation of the Potential of Application of Single Effect and Multiple Effect Absorption Cooling Systems. Energy Conversion and Management, 51, 1629-1636. https://doi.org/10.1016/j.enconman.2009.12.039
- 13. Kaynakli, O., Saka, K. and Kaynakli, F. (2015) Energy and Exergy Analysis of a Double Effect Absorption Refrigeration System Based on Different Heat Sources. Energy Conversion and Management, 106, 21-30. https://doi.org/10.1016/j.enconman.2015.09.010
- 14. Maryami, R. and Dehghan, A.A. (2017) An Exergy Based Comparative Study between LiBr/Water Absorption Refrigeration Systems from Half Effect to Triple Effect. Applied Thermal Engineering, 124, 103-123. https://doi.org/10.1016/j.applthermaleng.2017.05.174
- 15. Li, Z., Ye, X. and Liu, J. (2014) Performance Analysis of Solar Air Cooled Double Effect LiBr/H2O Absorption Cooling System in Subtropical City. Energy Conversion and Management. Energy Conversion and Management, 85, 302-312. https://doi.org/10.1016/j.enconman.2014.05.095
- 16. Hamed, A., Kaseb, S.A. and Hanaf, A.S. (2015) Prediction of Energetic and Exergetic Performance of Double-Effect Absorption System. International Journal of Hydrogen Energy, 40, 15320-15327. https://doi.org/10.1016/j.ijhydene.2015.05.150
- 17. Gomri, R. and Hakimi, R. (2008) Second Law Analysis of Double Effect Vapour Absorption Cooler System. Energy Conversion and Management, 49, 3343-3348. https://doi.org/10.1016/j.enconman.2007.09.033
- 18. Sharshir, S.W., Elsheikh, A.H., Peng, G., Yang, N., El-Samadony, M.O.A. and Kabeel, A.E. (2017) Thermal Performance and Exergy Analysis of Solar Stills—A Review. Renewable and Sustainable Energy Reviews, 73, 521-544. https://doi.org/10.1016/j.rser.2017.01.156
- 19. Tenkeng, M., Wouagfack, P.A.N., Lissouck, D. and Tchinda, R. (2017) Exergy Analysis of a Solar Absorption Refrigeration System in Ngaoundere. Journal of Power and Energy Engineering, 5, 1-18. https://doi.org/10.4236/jpee.2017.510001
- 20. Atmaca, I. and Yigit, A. (2003) Simulation of Solar-Powered Absorption Cooling System. Renewable Energy, 28, 1277-1293. https://doi.org/10.1016/S0960-1481(02)00252-5
- 21. Xu, Z.Y., Wang, R.Z. and Wang, H.B. (2015) Experimental Evaluation of a Variable Effect LiBr-Water Absorption Chiller Designed for High-Efficient Solar Cooling System. International Journal of Refrigeration, 59, 135-143.
- 22. Lizarte, R. and Marcos, J.D. (2016) COP Optimisation of a Triple-Effect H2O/LiBr Absorption Cycle under Off-Design Conditions. Applied Thermal Engineering, 99, 195-205. https://doi.org/10.1016/j.applthermaleng.2015.12.121
- 23. Anand, S., Gupta, A. and Tyagi, S.K. (2013) Simulation Studies of Refrigeration Cycles: A Review. Renewable and Sustainable Energy Reviews, 17, 260-277. https://doi.org/10.1016/j.rser.2012.09.021
- 24. Sencan, A., Yakut, K.A. and Kalogirou, S.A. (2005) Exergy Analysis of Lithium Bromide/Water Absorption Systems. Renewable Energy, 30, 645-657. https://doi.org/10.1016/j.renene.2004.07.006
- 25. De Vos, A. (1993) The Endoreversible Theory of Solar Energy Conversion: A Tutorial. Solar Energy Materials and Solar Cells, 31, 75-93. https://doi.org/10.1016/0927-0248(93)90008-Q
- 26. Pridasawas, W. and Lundqvist, P. (2004) An Exergy Analysis of a Solar-Driven Ejector Refrigeration System. Solar Energy, 76, 369-379. https://doi.org/10.1016/j.solener.2003.11.004
- 27. Sozen, A. (2001) Effect of Heat Exchangers on Performance of Absorption Refrigeration Systems. Energy Conversion and Management, 42, 1699-1716. https://doi.org/10.1016/S0196-8904(00)00151-5
- 28. Kaynakli, O. and Kilic, M. (2007) Theoretical Study on the Effect of Operating Conditions on Performance of Absorption Refrigeration System. Energy Conversion and Management, 48, 599-607. https://doi.org/10.1016/j.enconman.2006.06.005
- 29. Mittelman, G., Kribus, A. and Dayan, A. (2007) Solar Cooling with Concentrating Photovoltaic/Thermal (CPVT) Systems. Energy Conversion and Management, 48, 2481-2490. https://doi.org/10.1016/j.enconman.2007.04.004

Nomenclature

COP: coefficient of performance

ECOP: exergy coefficient of performance

f: circulation ratio

h: enthalpy (kJ/kg)

s: entropy (kJ/kgK)

m: mass flow rate (kg/s)

Q: heat flow rate (kW)

T: temperature (˚C or K)

P: pressure (kPa)

ψ: exergy (kJ/kg)

x: mass fraction of lithium bromide (%)

Subscripts

i: component i

ref: reference

g: ge: generator

e: ev: evaporator

a: ab: absorber

c: co: condenser

hex: heat exchanger

ev: expansion valve

pump: p: pump

p: planet

ws: w: weak solution

sc: solar collector

s: solution or system

sc, r: solar collector radiation

h: heat

p: planet

HEX I: heat exchanger

HEX II: heat exchanger

SEV: solution expansion valve

S1: weak solution at the inlet of GI

S2: weak solution at the inlet of GII

W: work

REV: refrigerant expansion valve

u: useful

ava: available.

sun: sun