A New Computational Tool for the Development of Advanced Exergy Analysis and LCA on Single Effect LiBr–H2O Solar Absorption Refrigeration System

A single effect LiBr–H2O absorption refrigeration system coupled with a solar collector and a storage tank was studied to develop an assessment tool using the built-in App Designer in MATLAB®. The model is developed using balances of mass, energy, and species conservation in the components of the absorption cooling system, taking into account the effect of external streams through temperature and pressure drop. The whole system, coupled with the solar energy harvesting arrangement, is modeled for 24 h of operation with changes on an hourly basis based on ambient temperature, cooling system load demand, and hourly solar irradiation, which is measured and recorded by national weather institutes sources. Test through simulations and validation procedures are carried out with acknowledged scientific articles. These show 2.65% of maximum relative error on the energy analysis with respect to cited authors. The environmental conditions used in the study were evaluated in Barranquilla, Colombia, with datasets of the Institute of Hydrology, Meteorology and Environmental Studies (IDEAM), considering multiannual average hourly basis solar irradiation. This allowed the authors to obtain the behavior of the surface temperature of the water in the tank, COP, and exergy efficiency of the system. The simulations also stated the generator as the biggest source of irreversibility with around 45.53% of total exergy destruction in the inner cycle without considering the solar array, in which case the solar array would present the most exergy destruction.


Introduction
Most of the air conditioning systems used globally are based on vapor compression refrigeration systems (VCRs), a thermal cycle that is also widely used in the refrigeration field. Nevertheless, despite its global dominance and extension, new concerns regarding the rising cost of primary energy, the increase in demand for space conditioning due to the increase in temperature from global warming, the limit on the installed capacity of the electricity networks, and a progressive awareness towards the decrease of environmental impact have brought new alternatives into consideration [1,2].
There are thermally driven air conditioning systems that can be powered by heat from industrial exhaust gases, energy from waste heat from processes, among others, which replace the compressor and electrical energy that drives it from the power grid. These systems are known as vapor absorption refrigeration systems (VARs), being the most 2 of 27 common temperature operated system, and typically use water-ammonium or lithium bromide-water as the refrigerant pairs. An emerging technology for VARs is solar cooling, which promotes an economically attractive energy source with zero or very low emissions and can help mitigate peak energy loads associated with refrigeration and space conditioning. Solar cooling can include the direct use of solar energy through the photovoltaic effect, as well as solar thermal collectors that use a thermally driven cooling device, such as an absorption chiller, or even a combination of both systems [3]. Both types can be adapted for space conditioning and can also be combined with waste heat, geothermal, cogeneration plants, among others [4,5].
Absorption solar thermal cooling technologies use the heat obtained from the sun (using different types of collectors) in a thermally driven cooling process. Solar thermal cooling systems using absorption cooling are so far relatively small technologies, but the market is growing as costs come down. By the end of 2015, an estimated 1350 solar thermal cooling systems had been installed worldwide, of which about 80% were in Europe, mainly in Spain, Germany, and Italy [1,[6][7][8]. Costs for air conditioning by absorption with solar collectors have been reduced by more than half since 2007, largely due to the standardization of equipment, although they have not decreased as rapidly as solar photovoltaics [9,10].
In addition to the above, the authors of this paper recognize the hard times all of us are living due to the COVID-19 pandemic [11]. In educational settings, face-to-face classes have been avoided in order to break the chain of contagion. Many researchers from around the globe have focused their efforts on developing content, creating strategies, and finding resources to continue the learning process and guarantee the minimum impact on the quality of education [12,13]. This strategy should be accompanied by students taking responsibility for their own learning. For this reason, the product of this article could be used as a learning environment for the assessment of first and second laws on a LiBr-H 2 O cooling absorption system.
Research and development in the field of renewable energy, together with energy efficiency, makes it possible to mitigate the damage generated by the emission of greenhouse gases. The main difference between the two branches is that energy efficiency devotes a large part of its research to the improvement of approaches that generate long-term energy savings with little or no investment. Conversely, renewable energy requires investment, which is usually not economical. Therefore, thermodynamic simulation for such systems is widely used to determine whether or not physical implementation is feasible [14,15]. Such research can be aimed at full use of renewable energy (which is somewhat similar to the aim of the solar VAR), partial use between renewable and non-renewable energy [16][17][18], or purely non-renewable use where a reduction in use is pursued [19,20].
The aim of this work was to present an efficient and easy to use computer tool that allows the simulation of a basic and advanced cooling system, by means of energy, exergy, and disaggregated energy parameters. The advanced model has a variation where it is coupled with a solar collection system, which allows evaluation of the performance of the absorption cooling system under climatic conditions reported by meteorological databases. The software proposed in this article was performed in the MATLAB App Designer application as a standalone application and was simulated with the input conditions established in other articles on the topic in order to measure the performance of the software by obtaining correspondence with those articles.

Description of the System
A vapor absorption refrigeration (VAR) system is a thermodynamic cycle similar to the popular vapor compression refrigeration (VCR) system that consists of four main components: evaporator, condenser, valve, and compressor. The main difference between both systems is that VAR replaces the compressor by a pseudo inner cycle where the Lubricants 2021, 9,76 3 of 27 absorption process takes place and a similar effect of the compressor is achieved, taking the refrigerant pressure from a low value from the evaporator to a high value for the condenser.
The (VAR) cycle, shown in Figure 1 coupled with the solar assembly, uses water and aqueous lithium bromide as refrigerants, which is also an important difference from the VCR, which uses hydrochlorofluorocarbons (HCFCs) as refrigerants [21]; water in the VAR cycle performs the entire system cycle, while LiBr only performs the inner cycle. The inner cycle includes five components, denoted as absorber (A1), where the water and LiBr are combined as an aqueous LiBr solution; a solution pump (P1); a solution heat exchanger (ITC1); a generator (G1), where most of the water is separate from LiBr; and a solution valve (V2), where a strong solution of aqueous LiBr is taken back to the absorber. both systems is that VAR replaces the compressor by a pseudo inner cycle where the ab-sorption process takes place and a similar effect of the compressor is achieved, taking the refrigerant pressure from a low value from the evaporator to a high value for the condenser.
The (VAR) cycle, shown in Figure 1 coupled with the solar assembly, uses water and aqueous lithium bromide as refrigerants, which is also an important difference from the VCR, which uses hydrochlorofluorocarbons (HCFCs) as refrigerants [21]; water in the VAR cycle performs the entire system cycle, while LiBr only performs the inner cycle. The inner cycle includes five components, denoted as absorber (A1), where the water and LiBr are combined as an aqueous LiBr solution; a solution pump (P1); a solution heat exchanger (ITC1); a generator (G1), where most of the water is separate from LiBr; and a solution valve (V2), where a strong solution of aqueous LiBr is taken back to the absorber.
The solar collector, coupled with a thermal storage tank, is used as an energy source for the VAR cycle, as shown in Figure 1. The fluid considered for the solar assembly is water; properties such as the collector area and type of collector can be modified within the simulation, as well as the thermal tank losses and volume. The desired hour of operation of the VAR, hourly irradiation at the solar collector, and hourly variation of the room temperature can be specified as an array or vector. The main objective of the simulation of the refrigeration cycle coupled to the solar system is to evaluate how the system will perform under user-specified conditions, so a MATLAB© App Designer application was developed in order to provide an intuitive but completed software to do so. However, the functionality of the software can be informative in educational courses such as thermodynamics or even in courses about how to develop software for engineering. A longer-term objective with software is to provide awareness. According to the International Agency of Energy (IAE), traditional refrigeration with VCR is a blind-spot in global energy politics due to the fact that as the economy grows, more VCR will be needed for better human comfort, but traditional refrigeration with VCR is slowly taking a more and more energy from the grid, causing a need for more energy, which leads to the construction of additional power plants and ultimately more The solar collector, coupled with a thermal storage tank, is used as an energy source for the VAR cycle, as shown in Figure 1. The fluid considered for the solar assembly is water; properties such as the collector area and type of collector can be modified within the simulation, as well as the thermal tank losses and volume. The desired hour of operation of the VAR, hourly irradiation at the solar collector, and hourly variation of the room temperature can be specified as an array or vector.
The main objective of the simulation of the refrigeration cycle coupled to the solar system is to evaluate how the system will perform under user-specified conditions, so a MATLAB© App Designer application was developed in order to provide an intuitive but completed software to do so. However, the functionality of the software can be informative in educational courses such as thermodynamics or even in courses about how to develop software for engineering. A longer-term objective with software is to provide awareness. According to the International Agency of Energy (IAE), traditional refrigeration with VCR is a blind-spot in global energy politics due to the fact that as the economy grows, more VCR will be needed for better human comfort, but traditional refrigeration with VCR is slowly taking a more and more energy from the grid, causing a need for more energy, which leads to the construction of additional power plants and ultimately more greenhouse effect [22,23], which in turns increases the need of refrigeration in warm regions. This effect can be confronted by renewable technologies [24,25] and energy-saving techniques [26,27].
The software is subdivided into three phases, the first one provides a simple cycle of VAR, where no external flow is considered. This model allows a quick calculation, whose output is the power per component, the properties per state, and a TS diagram. This model is called model 1 and its validation is presented in the results section. The second model (model 2) allows a more detailed study, where there are more inputs, such as pressure drops, temperature differences, flow and temperature of external streams, cooling loads, Lubricants 2021, 9, 76 4 of 27 etc., with outputs such as energy, exergy, and exergy disaggregation. Its validation is also presented in the results section. In the last stage, the coupling of model 2 with the solar array is carried out. In this stage, a case study of the irradiation and ambient temperature of the city of Barranquilla is analyzed.

Thermodynamics Analysis: Vapor Absorption Refrigeration System
Considering the mass conservation law and the steady-state consideration assumed to all components of the VAR, the mass (Equation (1)), species (Equation (2)), and energy (Equation (3)) conservation laws were applied.
Enthalpy and entropy in the states of pure water in the cycle are functions of two independent intensive properties, usually temperature and pressure, while for the states where LiBr and water are involved, the properties also depend on the water concentration in LiBr:

Pure Water
LiBr-Water Pure Water LiBr-Water s = s(T, P) s = s(T, P, X) It is widely accepted in the analysis of this cycle that, at the output of the generator to the condenser, the refrigerant does not contain lithium bromide [6], meaning that the refrigerant is found as a stream of pure water, so: which means that for a concentration X = 1, there is pure water, while a concentration X below 1 means that there is an amount 1 − X of LiBr in the state. The generator, absorber, condenser, and evaporator stream outputs are considered to be in saturated condition: m 15 X 9 = X 11 = X 13 = X 15 = 1 Cooling Load . Q eva (7) Due to the extensive design possibilities, there are many considerations regarding the steam state in which the refrigerant exits the generator. For this model, we adopt the consideration presented by Wonchala et al. [28], where the generator refrigerant outlet temperature is considered to be equal to the equilibrium temperature of the weak solution: The minimum temperature difference in the absorber, generator, condenser, and evaporator and the efficiency of the solution heat exchanger are input values: ∆T gen = T 13 − T 2 (10) Pressure drops by length and those due to the individual components of the absorption cooling system, with the exception of the pump and expansion valves, which are usually negligible. For the present modeling, following the recommendations presented by Morosuk and Tsatsaronis 2008 [29], the pressure drops in the generator and the absorber will be taken into account; this supports the determination of ideal, real, and inevitable operations that are useful for the disaggregation of exergy destruction: Solution pump efficiency is defined as: where h 1s is the outlet enthalpy of the weak solution if the pump will operate isentropically. From the properties calculated at each point in the absorption cooling system, the system performance (COP) is defined as: The maximum possible efficiency for an absorption cooling system (reversible COP) operating between the temperature of the heat source, the cooling fluid, and the cooled environment is defined as: (20) To determine the thermodynamic properties of water, it was implemented in the library of properties of CoolProp. Coolprop is a software that has an extensive database of correlations and equations needed to determine the properties of various pure and pseudo pure substances, and it is completely free-to-use [30]. CoolProp contained the properties of the LiBr-H 2 O solution. However, these properties did not correspond satisfactorily with the validation of different works carried out during the development of the simulation software, so it was decided to use the properties of the LiBr-H 2 O solution programmed manually according to the tables and correlations presented in the scientific literature.
The enthalpy h and entropy s of the LiBr-H 2 O solution are described by the parameters and constants resulting from the implementation of the Gibbs free energy equation on electrolyte substances by Kim and Infante for a concentration range between 0 and 70 wt% of lithium bromide in water and a temperature range between 0 and 210 • C [31].
Viscosity, concentration, density, specific heat, and conductivity were taken from the technical report of the thermophysical properties for the LiBr-H 2 O solution of the American Society of Heating, Refrigerating, and Air-Conditioning Engineers, ASHRAE [32]. Table 1 presents the properties of the LiBr-H 2 O solution, commented on before, which is attached to the admissible range of the independent input variables and literary sources where they are used: The total flow exergy ψ of a state, considering that potential and kinetic components are negligible, is calculated as the sum of the physical and chemical exergy: The physical exergy of a flow ψ ph is the maximum work available when the system is taken from an initial state with conditions (T, p) to a reference state (T 0 , p 0 ) by a reversible process, exchanging energy in the form of heat and working only with the reference environment; physical exergy can be calculated using the following equation: where h. is the enthalpy, s. is the entropy, and the sub-index zero indicates the property is in a dead state at reference temperature (T 0 ). For water, the above equation depends on the temperature and pressure conditions, while for the LiBr-H 2 O solution, it also depends on the specific concentration.
The chemical exergy of a flow ψ ch is the maximum work that can be achieved when the substance is brought from its equilibrium state at ambient pressure and temperature (restricted dead state) to the equilibrium state of equal chemical potential (unrestricted dead state) by a process involving heat, work, and mass transfer with the environment.
For those states where only water is present, the chemical exergy of flow is taken with a fixed value of 49.96 kJ/kg, according to [34,35]. Since the LiBr-H 2 O solution is not ideal, the exergy of the chemical flow for the solution is a function of the standard activities and exergies of the pure species, as represented below: As mentioned above, Equation (23) consists of two parts, the standard chemical exergy of pure species (Equation (24)) and chemical exergy destroyed due to the dissolution process (Equation (25)): where ε 0 H 2 O and ε 0 LiBr represent the standard chemical exergies of water and lithium bromide, whose values are 0.9 kJ/mol and 101.6 kJ/mol, respectively.
where a H 2 O and a LiBr represent water and lithium bromide activity, respectively. A flowchart of the process to calculate the total exergy for LiBr-Water is presented in Figure 2; for a detailed explanation of each term, check [35]. After the determination and addition of both components of exergy, the total exergy expressed in terms of rate is obtained by multiplying the total flow exergy with the flow rate of that state: m ·ψ (26) Lubricants 2021, 9, x FOR PEER REVIEW 7 of 27 where ̃ and represent the standard chemical exergies of water and lithium bromide, whose values are 0.9 kJ/mol and 101.6 kJ/mol, respectively.
where and represent water and lithium bromide activity, respectively. A flowchart of the process to calculate the total exergy for LiBr-Water is presented in Figure  2; for a detailed explanation of each term, check [35]. After the determination and addition of both components of exergy, the total exergy expressed in terms of rate is obtained by multiplying the total flow exergy with the flow rate of that state: The definition of fuel-product-destroyed is applied to traditional exergy analysis, where the fuel is the amount of exergy that enters a component to produce an amount of product. Similarly, the product is defined as the amount of exergy left by a component converted by the fuel that previously entered the same component. For the specific case of the component under study k, the exergy of fuel, product, and destruction is given by Equation (27). Table 2 shows in detail the structure of fuel-product-destroyed exergy balance in each component of the absorption refrigeration system:  The definition of fuel-product-destroyed is applied to traditional exergy analysis, where the fuel is the amount of exergy that enters a component to produce an amount of product. Similarly, the product is defined as the amount of exergy left by a component converted by the fuel that previously entered the same component. For the specific case of the component under study k, the exergy of fuel, product, and destruction is given by Equation (27). Table 2 shows in detail the structure of fuel-product-destroyed exergy balance in each component of the absorption refrigeration system: Table 2. Fuel-product-destroyed structure for each component in the refrigeration cycle.
The ratio between the definitions of COPs set out by Equations (19) and (20) allows the total exergetic efficiency of the system to be defined, which represents a measure of the thermodynamic irreversibilities associated with the absorption refrigeration cycle: Exergetic efficiency η I I for the k component, in particular, is defined as the relationship between the produced exergy . E P,k over the required exergy . E F,k : The ratio of irreversibility (IR) of each component of the absorption refrigeration cycle is defined as the ratio of exergy destroyed in the particular component and the total exergy destroyed in the system: The exergy depletion ratio (FDR) is defined as the ratio of exergy destroyed in the particular component and the total exergy supplied to the system, in this case the exergy supplied to the generator:

Exergy Disaggregation
An exergy disaggregation allows a more detailed investigation of the cause of the exergy destruction, with the aim of identifying opportunities for improvement in each component of the system [36,37]. In order to disaggregate the irreversibilities in the absorption cooling system, the theory of advanced exergy analysis proposed by [38] is implemented, which is based on two postulates. The first proposes that the exergy destroyed can be divided into two categories according to the relevance of the component in the exergy destruction in the said component. In the first category, the exergy destruction is due to the operation of the component under consideration, which is called endogenous exergy destruction, while in the second category, there is a destruction of exergy of the component under consideration but due to the operation of the other components, which is known as exogenous exergy destruction: .
Equation (32) determines whether efforts to make system improvements should focus on the component with the greatest exergy destruction, or whether the irreversibilities associated with the component with the greatest exergy destruction can be reduced more efficiently by improving operation in the other components of the cycle. The second postulate disaggregates or categorizes the irreversibilities of the components of the cycle as avoidable and unavoidable: .
Equation (33) indicates that efforts to make improvements to the system should be focused on components with appreciable avoidable irreversibilities and should not be focused on those components or processes where inevitable irreversibilities predominate.
The combined implementation of the two postulates of Equations (32) and (33) allows the destruction of exergy of the k-component to be described as: .
Equation (34)  respectively stand for the avoidable and unavoidable destroyed exergy of the analyzed component-not due to itself but because of the interactions of other components of the cycle over this one.
The thermodynamics analysis development in the previous section features temperature differences in the condenser, evaporator, absorber, and generator as well as pressure drops in the absorber and generator (Equations (9)- (12) and Equations (14)- (17)). These configurations allow the definition of real, unavoidable, and ideal cycle conditions that allow executing the cycle with theoretical values to measure the real performance against a theoretical one. Thus, by defining real temperature differences and pressure drops, a reference model is obtained. In contrast, it is possible to define a theoretical cycle with the same characteristics as the real cycle but with no temperature difference and no pressure drop in the components, in addition to operating with isentropic turbines instead of the valves, such as the Carnot refrigeration cycle. These configurations allow the definition of ideal cycle conditions. Finally, an intermediate configuration is defined with very small temperature differences and pressure drops in the components. This configuration allows the definition of a cycle condition with unavoidable irreversibilities or the unavoidable cycle. . With these data, it is possible to make further improvement by focusing on the component of the system that has the most potentially avoidable exergy destroyed and to check whether the improvement could be achieved by focusing only on that component or on the whole cycle. This is an advantage over conventional exergy analysis, which can create motivation to improve a component that has the total exergy destroyed; however, maybe all that exergy is of an inevitable nature and, together with concepts such as thermoeconomic analysis, where each exergy destroyed is computed against a calculated cost, it can become a powerful tool for determining energy and economic improvements [39,40].

Solar Collector and Storage Tank
The model of the flat-plate solar collector is based on the equations presented by [41]. The rate of useful energy transported by the collector is defined as: where . Q u is the useful energy collected by the solar system, A c is the area of the collector, I t is the total radiation incident on the collector, and η c is the instantaneous efficiency of the collector, which is defined as the ratio of useful heat over a specified period of time to the energy of the incident solar irradiation over the same period of time. Instantaneous collector efficiency is defined as: where F R U L is the collector heat removal factor, which is calculated as the ratio of the actual rate of heat transfer to the working fluid to the actual rate of heat transfer at the minimum temperature difference between the collector absorber element and the environment; while τα is the product of transmittance−absorbance, or optical efficiency, where the transmittance τ is the fraction of solar incidence radiation transmitted through the collector cover, and the absorbance α is the fraction of solar incidence radiation absorbed by the collector material. In addition, in Equation (39), U L is the overall transfer coefficient, T f ,in is the temperature of the fluid at the collector inlet, and T a is the ambient temperature. Usually, the products F R (τα) and F R U L are characterized by a particular solar collector and dictate the performance of the collector. Their values are tabulated and obtained from experimental information in graphic form of the collector's performance η c against (T f ,in −T a ) I T , where F R (τα) is the intercept, and F R U L is the slope [41].
Data for total hourly solar irradiation I T is taken from metrological sources, such as the Institute of Hydrology, Meteorology and Environmental Studies (IDEAM), which have the hourly irradiation per month for different locations in Colombia for different years [42]; but, as can be seen later in the software section, the irradiation and the room temperature can be manually modified to assess the solar collector gain and therefore to assess the solar absorption system performance anywhere.
A thermal storage tank is integrated with the collection system in order to collect solar energy and provide a continuous supply of heat when solar radiation is no longer sufficient. To simplify the model, it is assumed that the water in the insulated tank is completely mixed with the water returning to the tank from the collector and generator. By performing an energy balance in the tank without stratification (i.e., completely mixed and without temperature variations lengthwise), the following equation is obtained: where m is the mass of water in the storage tank, C p is the calorific capacity of water, and T a is the ambient temperature around the tank; whereas, . Q u and . Q load represent the useful thermal energy obtained from the solar collector and the thermal energy required by the absorption cooling system (thermal load required by the generator), respectively.
Additionally, (U A) s represents the product coefficient that quantifies the losses per storage tank area. Its value is taken as 11.1 W/ • C [43]. However, it is not fixed and can be editable in the software. If . Q u , . Q load and the tank losses over a period of time ∆t are assumed to be constant, Equation (37) can be rewritten for each interval as: where T s,new is the temperature of the storage tank at the end of that time interval ∆t, and T s is the temperature of the storage tank during the ∆t. At the end of each time interval ∆t, the value of T s,new obtained will be assigned as the T s for the next interval, and so on for the time extension of interest. The rate of exergy absorbed from solar radiation is expressed, according to [44], as: .
where η o represents the optical efficiency, which is equal to the transmittance-absorbance product τα, T b is the apparent temperature of the sun, which is equal to 75% of the black body temperature of the sun [45]. The definition of FDR from Equation (31) is modified to take into account the total exergy entering the system from the solar source, in the form:

LCA Solar Absorption Refrigeration System
The life cycle analysis developed is a tool used for the systematic assessment of the likely environmental impacts of products, services, and technologies. This is done by examining inputs and desired emissions, wastes, and/or outputs throughout the entire product life cycle (Figures 3 and 4). Through this analysis, environmental burdens are identified, and their propagation from one stage of the life cycle to another or from one environmental threat to another is prevented. The implemented LCA is a standardized method that provides an analytical basis for the evaluation and comparison of products, processes, materials, and service systems. where , is the temperature of the storage tank at the end of that time interval ∆ , and is the temperature of the storage tank during the ∆ . At the end of each time interval ∆ , the value of , obtained will be assigned as the for the next interval, and so on for the time extension of interest.
The rate of exergy absorbed from solar radiation is expressed, according to [44], as: where represents the optical efficiency, which is equal to the transmittance-absorbance product , is the apparent temperature of the sun, which is equal to 75% of the black body temperature of the sun [45].
The definition of FDR from Equation (31) is modified to take into account the total exergy entering the system from the solar source, in the form:

LCA Solar Absorption Refrigeration System
The life cycle analysis developed is a tool used for the systematic assessment of the likely environmental impacts of products, services, and technologies. This is done by examining inputs and desired emissions, wastes, and/or outputs throughout the entire product life cycle (Figures 3 and 4). Through this analysis, environmental burdens are identified, and their propagation from one stage of the life cycle to another or from one environmental threat to another is prevented. The implemented LCA is a standardized method that provides an analytical basis for the evaluation and comparison of products, processes, materials, and service systems. The LCA involved the specifications of the absorption refrigeration system and compiled all relevant energy inputs, materials, and environmental releases in a quantitative inventory. The inputs and outputs of the system were evaluated taking into account the environmental impact they may have on human health, ecosystems, and resource depletion. The interpretation of these results provides further information for decision making regarding the industrial application of this equipment.
ISO 14040: 2006 [46] is used to elaborate the LCA of the solar absorption refrigeration system in four phases: definition of objective and scope, where the system boundaries and functional unit are established; inventory analysis, in which the input and output data pertaining to the product life cycle are quantified; impact assessment, which associates the inventory data with environmental impact categories; and interpretation, where the

MATLAB App Designer Application
A graphical user interface was developed in MATLAB using the App Designer application for the energy, exergy, and disaggregation stages of the VAR cycle of LiBr-H2O, which consists of three tools that perform the analyses presented in the previous section. A flow chart schematic of the software is presented in Figure 5. Figure 6 shows a preview of some of the graphical user interfaces of the software packages.
This software was registered with the number 13-79-348 in the DNDA. In order to execute the program, it is necessary to have already installed MATLAB Runtime and Coolprop in the default folder and Windows operating system with a least 100 Mb of free space on the local disk.
Model 1 (Figure 6a) allows validation of the LiBr-H2O single-effect absorption refrigeration system for basic conditions, i.e., when only information about the temperature of the four internal components, mass flows, and heat exchanger efficiency is presented. It is important to note that the number of models and considerations presented in the literature for this system is vast, so the results may differ from one particular study to another. Since the model does not include external streams, the energy and exergy capacity of the calculation is limited to providing the values of the thermal power used by the components, the properties by states of the system, and a temperature-entropy diagram of the system. The development of this model is based on the study presented by [47]. This model has a Specific Analysis mode, where detailed results are presented, and a Parametric Analysis mode, where it is possible to vary the temperature of a particular component and see its impact on the performance variable for this particular model, COP.

Definition of goals and objectives
Inventory analysis The LCA involved the specifications of the absorption refrigeration system and compiled all relevant energy inputs, materials, and environmental releases in a quantitative inventory. The inputs and outputs of the system were evaluated taking into account the environmental impact they may have on human health, ecosystems, and resource depletion. The interpretation of these results provides further information for decision making regarding the industrial application of this equipment.
ISO 14040: 2006 [46] is used to elaborate the LCA of the solar absorption refrigeration system in four phases: definition of objective and scope, where the system boundaries and functional unit are established; inventory analysis, in which the input and output data pertaining to the product life cycle are quantified; impact assessment, which associates the inventory data with environmental impact categories; and interpretation, where the results are summarized and conclusions are drawn.

MATLAB App Designer Application
A graphical user interface was developed in MATLAB using the App Designer application for the energy, exergy, and disaggregation stages of the VAR cycle of LiBr-H 2 O, which consists of three tools that perform the analyses presented in the previous section. A flow chart schematic of the software is presented in Figure 5. Figure 6 shows a preview of some of the graphical user interfaces of the software packages.
This software was registered with the number 13-79-348 in the DNDA. In order to execute the program, it is necessary to have already installed MATLAB Runtime and Coolprop in the default folder and Windows operating system with a least 100 Mb of free space on the local disk.
Model 1 (Figure 6a) allows validation of the LiBr-H2O single-effect absorption refrigeration system for basic conditions, i.e., when only information about the temperature of the four internal components, mass flows, and heat exchanger efficiency is presented. It is important to note that the number of models and considerations presented in the literature for this system is vast, so the results may differ from one particular study to another. Since the model does not include external streams, the energy and exergy capacity of the calculation is limited to providing the values of the thermal power used by the components, the properties by states of the system, and a temperature-entropy diagram of the system. The development of this model is based on the study presented by [47]. This model has a Specific Analysis mode, where detailed results are presented, and a Parametric Analysis mode, where it is possible to vary the temperature of a particular component and see its impact on the performance variable for this particular model, COP.   The second tool or model 2 ( Figure 6b) allows a complete energy and exergy analysis of the absorption cooling system to be performed, to the point of being able to break down the destroyed exergy of the components into four types (endogenous, exogenous, avoidable, and unavoidable) by implementing the real, ideal, and unavoidable conditions due to pressure drops and temperature differences in the same components, as well as taking into account the effect of external streams over the characteristic components of the LiBr-H 2 O single-effect absorption cooling system. This tool also features specific analysis for detailed performance information for fixed conditions and parametric analysis for conditions where external streams temperatures vary. The modeling of this tool is based on the model presented [33,34].
The third tool or solar analysis (Figure 6c) allows a dimensioning of a solar collection system composed of a thermal storage tank and a solar collector coupled with a LiBr-H 2 O absorption cooling system similar to the one presented in model 2. The tool simulates operation over a 24 h period, where the storage tank will supply temperature to the cooling system as if it were a thermal battery. The tool allows indication of whether the cooling system operates continuously or is turned on for a particular period of time, simulating a more appropriate and conscious use of the system. It is important to note that, due to the nature of change, the utility of the tool can be interpreted as a model in a transient state. However, the model actually applies a stationary operation at each hourly interval, so the properties during one hour of execution will not change until the next hour begins.  The second tool or model 2 ( Figure 6b) allows a complete energy and exergy analysis of the absorption cooling system to be performed, to the point of being able to break down the destroyed exergy of the components into four types (endogenous, exogenous, avoidable, and unavoidable) by implementing the real, ideal, and unavoidable conditions due to pressure drops and temperature differences in the same components, as well as taking into account the effect of external streams over the characteristic components of the LiBr-

Properties Validation
For the development of the software, correlations from different sources were used, as shown in Table 1, referring to the LiBr-H 2 O. It is of interest to validate the accuracy of these properties since the thermal properties are a function of the type of model used by each author and since different thermodynamic modeling approaches have been found for absorption cooling systems [6,48].
Due the clear explanation of thermodynamical modelling of Gong and Boulama [49] and Wonchola [28], their model was adopted; however, as mentioned, the thermal properties are taken from different sources, so the aim of this section is to validate such properties in the particular case presented in [49,50]. Table 3 presents a comparison regarding the flow-rate, concentration of LiBr, and temperature, while Table 4 presents pressure, enthalpy, and entropy; the maximum relative error found was of 2.65% for two states referring to enthalpy, followed by 1.61% for two states of entropy and 1.45% for two states of flow-rate.

Software Validation
The tab of tool number 1 was validated with the article by [47]; in Figure 7, it can be seen how the performance of a system with two temperature conditions for the absorber can be replicated while varying the temperature of the generator.
The second tool of the program was validated with [49]. Figure 8 shows how the same input conditions, reported in Gong's system, are established in the software. Figure 9 shows the properties' outputs per state and the energy and exergy per component of the tool against the results reported by Gong. Figure 10 corresponds to the output of the parametric analysis option, where one can see the effect of the temperature variation of one component on another or on the global system, taking as a reference the performance, exergy efficiency, total exergy destroyed, or its disaggregation by component.

Software Validation
The tab of tool number 1 was validated with the article by [47]; in Figure 7, it can be seen how the performance of a system with two temperature conditions for the absorber can be replicated while varying the temperature of the generator.
The second tool of the program was validated with [49]. Figure 8 shows how the same input conditions, reported in Gong's system, are established in the software. Figure 9 shows the properties' outputs per state and the energy and exergy per component of the tool against the results reported by Gong. Figure 10 corresponds to the output of the parametric analysis option, where one can see the effect of the temperature variation of one component on another or on the global system, taking as a reference the performance, exergy efficiency, total exergy destroyed, or its disaggregation by component.

Solar Analysis
The last software tool allows making a solar analysis of the absorption cooling system for a period of 24 h of operation by entering information similar to model 2. For this purpose, the hourly array configuration is presented in Figure 11, the solar collector and storage tank configuration are presented in Figure 12, and the VAR configuration in Figure 13.

Solar Analysis
The last software tool allows making a solar analysis of the absorption cooling system for a period of 24 h of operation by entering information similar to model 2. For this purpose, the hourly array configuration is presented in Figure 11, the solar collector and storage tank configuration are presented in Figure 12, and the VAR configuration in Figure 13.   As can be seen in Figure 14, a large dataset is generated for each hour of operation for the coupled system, including the collector generation, each component's power, exergy destroyed, irreversibility ratio (IR), and exergy depletion ratio (FDR), along with a plot solar irradiance per hour against COP and exergy efficiency, and finally, parametrical plots of the exergy disaggregation in the selected component for each hour of operation previously defined in the logical state of operation.  As can be seen in Figure 14, a large dataset is generated for each hour of operation for the coupled system, including the collector generation, each component's power, exergy destroyed, irreversibility ratio (IR), and exergy depletion ratio (FDR), along with a plot solar irradiance per hour against COP and exergy efficiency, and finally, parametrical plots of the exergy disaggregation in the selected component for each hour of operation previously defined in the logical state of operation. As can be seen in Figure 14, a large dataset is generated for each hour of operation for the coupled system, including the collector generation, each component's power, exergy destroyed, irreversibility ratio (IR), and exergy depletion ratio (FDR), along with a plot solar irradiance per hour against COP and exergy efficiency, and finally, parametrical plots of the exergy disaggregation in the selected component for each hour of operation previously defined in the logical state of operation. This particular case indicates that when the VAR system is turned up for high-temperature hours, the solar irradiance will not be enough to keep a stable temperature in the This particular case indicates that when the VAR system is turned up for hightemperature hours, the solar irradiance will not be enough to keep a stable temperature in the storage tank water, and therefore the thermal quality of water exchanging heat with the generator will drop slowly, causing the need of more thermal power for the fixed cooling load required of 10 kW. This case can be easily resolved by implementing an electrical resistance as a way to supply the temperature difference to guarantee a stable cooling supply without variation of thermal power in the generator. However, this would make outdated the energy and exergy output of the model since it would always be the same and could just be modeled by model 1 or model 2, while the only new unknown to solve would be the power needed to supply in each case due to the electrical resistance. Executing a different scenario, shown in Figure 15, the same model but with a 3.33 kW of required cooling load, gives a more stable performance in the VAR. storage tank water, and therefore the thermal quality of water exchanging heat with the generator will drop slowly, causing the need of more thermal power for the fixed cooling load required of 10 kW. This case can be easily resolved by implementing an electrical resistance as a way to supply the temperature difference to guarantee a stable cooling supply without variation of thermal power in the generator. However, this would make outdated the energy and exergy output of the model since it would always be the same and could just be modeled by model 1 or model 2, while the only new unknown to solve would be the power needed to supply in each case due to the electrical resistance. Executing a different scenario, shown in Figure 15, the same model but with a 3.33 kW of required cooling load, gives a more stable performance in the VAR.

LCA Solar Absorption Refrigeration System Results
The contribution to the climate change impact category is mainly due to the CO2 emissions generated in the construction process of the refrigeration system equipment, more specifically in the steel used. The results of the life cycle analysis in the evaluated case of the solar absorption refrigeration system show that the largest generator of potential emissions that impact climate change is the energy used in the equipment where unit operations occur, such as heat exchangers, contributing 48.1% of total emissions by the energy line used. Figure 16 shows the results of the evaluated inventory, thus allowing us to observe the potential kg of CO2 equivalents that are sent to the environment throughout the entire life cycle. These results are established by processes ( Figure 16a) and substances ( Figure  16b) that are involved in this category. Figure 16a shows the kg of CO2 equivalent potential emitted by the processes, showing that the base case of the plant evaluated in this study emits 1 kg of CO2 more per kg processed than the reference plant for this study due to the processes that are cogenerated with the production of thermal energy through solar energy, which shows the characteristics of the resources used in this type of refrigeration equipment. In the analysis of the impacts generated by the substances as shown in Figure 16b, it can be observed that the carbon dioxide emitted into the atmosphere by the use of fossil fuels emits approximately 4.4 kg of CO2 equivalent more than the reference plant used by the SimaPro ® database, showing the inefficiencies of the processes that are generated by the equipment used in Figure 15. VAR cycle coupled with the solar assembly outputs.

LCA Solar Absorption Refrigeration System Results
The contribution to the climate change impact category is mainly due to the CO 2 emissions generated in the construction process of the refrigeration system equipment, more specifically in the steel used. The results of the life cycle analysis in the evaluated case of the solar absorption refrigeration system show that the largest generator of potential emissions that impact climate change is the energy used in the equipment where unit operations occur, such as heat exchangers, contributing 48.1% of total emissions by the energy line used. Figure 16 shows the results of the evaluated inventory, thus allowing us to observe the potential kg of CO 2 equivalents that are sent to the environment throughout the entire life cycle. These results are established by processes ( Figure 16a) and substances (Figure 16b) that are involved in this category. Figure 16a shows the kg of CO 2 equivalent potential emitted by the processes, showing that the base case of the plant evaluated in this study emits 1 kg of CO 2 more per kg processed than the reference plant for this study due to the processes that are cogenerated with the production of thermal energy through solar energy, which shows the characteristics of the resources used in this type of refrigeration equipment. In the analysis of the impacts generated by the substances as shown in Figure 16b, it can be observed that the carbon dioxide emitted into the atmosphere by the use of fossil fuels emits approximately 4.4 kg of CO 2 equivalent more than the reference plant used by the SimaPro ® database, showing the inefficiencies of the processes that are generated by the equipment used in the process. In addition, we can observe that there are fields such as P11 and S5 where negative values in emissions are generated, which is due to the fact that in the process of converting solar energy to thermal energy through these systems, these processes and substances have no impact on the category of climate change because no emissions are generated.
Lubricants 2021, 9, x FOR PEER REVIEW 23 of 27 the process. In addition, we can observe that there are fields such as P11 and S5 where negative values in emissions are generated, which is due to the fact that in the process of converting solar energy to thermal energy through these systems, these processes and substances have no impact on the category of climate change because no emissions are generated.

Conclusions
A computational tool that allows simulating a vapor absorption refrigeration system was created through the use of a new MATLAB GUI application, App Designer. A reliable accuracy was obtained in the context of previously published research, although using different properties for the lithium bromide-water solution. This software can be executed as a standalone application. Its developed aims are to help the implementation of a new renewable system in the Caribbean Colombian zone, especially a renewable system that is fueled by solar energy because the Caribbean region benefits from high solar irradiance.
The software shows a good correspondence when calculating the energy, exergy, and advanced exergy components for the particular set of inputs proposed by Gong and Boulama [49] and Wonchala [28]. A solar harvesting system model was coupled to the vapor absorption refrigeration system, and under the user-defined ambient temperature, irradiance, and hours of working vectors, it is possible to visualize how such a system would work under 24 h of operation.
Due to the lack technical and theoretical knowledge, vapor absorption refrigeration (VARs) systems coupled with solar assembly have not been widely taken into account as a possible air conditioning solution in those regions where there is a good solar irradiance. In addition, due to the same reasons, there are not many manufactures and local industries; however, it seems that this situation will only last for a few more years because countries from Europe and Asia have been promoting such technologies, achieving a good impact in national industry and among domestic customers. It is the right time for this technology to be adopted in the Americas and specifically in the Caribbean area, where there is high solar irradiation that can be usefully employed in this system and where clean energy production can ultimately benefit human comfort.

Conclusions
A computational tool that allows simulating a vapor absorption refrigeration system was created through the use of a new MATLAB GUI application, App Designer. A reliable accuracy was obtained in the context of previously published research, although using different properties for the lithium bromide-water solution. This software can be executed as a standalone application. Its developed aims are to help the implementation of a new renewable system in the Caribbean Colombian zone, especially a renewable system that is fueled by solar energy because the Caribbean region benefits from high solar irradiance.
The software shows a good correspondence when calculating the energy, exergy, and advanced exergy components for the particular set of inputs proposed by Gong and Boulama [49] and Wonchala [28]. A solar harvesting system model was coupled to the vapor absorption refrigeration system, and under the user-defined ambient temperature, irradiance, and hours of working vectors, it is possible to visualize how such a system would work under 24 h of operation.
Due to the lack technical and theoretical knowledge, vapor absorption refrigeration (VARs) systems coupled with solar assembly have not been widely taken into account as a possible air conditioning solution in those regions where there is a good solar irradiance. In addition, due to the same reasons, there are not many manufactures and local industries; however, it seems that this situation will only last for a few more years because countries from Europe and Asia have been promoting such technologies, achieving a good impact in national industry and among domestic customers. It is the right time for this technology to be adopted in the Americas and specifically in the Caribbean area, where there is high solar irradiation that can be usefully employed in this system and where clean energy production can ultimately benefit human comfort.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: Indicates the property is in reference conditions. k K-th component F, K Fuel of K-th component P, K Product of K-th component D, K Destruction of K-th component ch 0 Standard chemical property of pure species ch; dis Chemical property destroyed due to the dissolution process