Mathematical Problems in Engineering
Volume 2020 (2020), Article ID 9098709, 11 pages
https://doi.org/10.1155/2020/9098709

Research Article

Modeling a Thermochemical Reactor of a Solar Refrigerator by BaCl₂-NH₃ Sorption Using Artificial Neural Networks and Mathematical Symmetry Groups

Onesimo Meza-Cruz,1,2 Isaac Pilatowsky,² Agustín Pérez-Ramírez,¹ Carlos Rivera-Blanco,² Youness El Hamzaoui,¹ Miguel Perez-Ramirez,³ and Mauricio A. Sanchez⁴

¹Universidad Autónoma Del Carmen, 24039 Cd Del Carmen Campeche, Mexico
²Instituto de Energías Renovables, Universidad Nacional Autónoma de México, 62580 Temixco, Morelos, Mexico
³Instituto Nacional de Electricidad y Energías Limpias, Reforma 113, Col. Palmira, C. P. 62490 Cuernavaca, Morelos, Mexico
⁴Facultad de Ciencias Químicas e Ingeniería, Universidad Autónoma de Baja California, Tijuana, Mexico

Correspondence should be addressed to Agustín Pérez-Ramírez; nitsug@hotmail.com

Received 26 April 2020; Revised 26 July 2020; Accepted 10 August 2020; Published 30 September 2020

Guest Editor: Qian Zhang

Copyright © 2020 Onesimo Meza-Cruz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The aim of this work is to present a model for heat transfer, desorbed refrigerant, and pressure of an intermittent solar cooling system’s thermochemical reactor based on backpropagation neural networks and mathematical symmetry groups. In order to achieve this, a reactor was designed and built based on the reaction of BaCl₂-NH₃. Experimental data from this reactor were collected, where barium chloride was used as a solid absorbent and ammonia as a refrigerant. The neural network was trained using the Levenberg–Marquardt algorithm. The correlation coefficient between experimental data and data simulated by the neural network was r = 0.9957. In the neural network’s sensitivity analysis, it was found that the inputs, reactor’s heating temperature and sorption time, influence neural network’s learning by 35% and 20%, respectively. It was also found that, by applying permutations to experimental data and using multibase mathematical symmetry groups, the neural network training algorithm converges faster.

1. Introduction

It is estimated that 15% of the electrical energy in the world is used for cooling [1]. One way to decrease this percentage is to use energy from alternative sources such as geothermal or solar energy, among others [2]. In [3], sorption cooling is proposed (adsorption) which uses low-intensity sources such as solar energy. Sorption cooling systems use thermal energy and apply thermodynamic equilibrium: liquid-vapor and solid-vapor. The best known liquid-vapor works are mentioned in [4, 5]. The solid-vapor system is divided into physical sorption and chemisorption. The main physical sorption works are mentioned in [6, 7].

The solids frequently used are metal salts, in particular, alkaline metals and alkaline earth metal halides, which, in anhydrous state, can absorb large amounts of ammonia or other refrigerants. The salts that are most used as sorbents are calcium chloride (CaCl₂) [8–10], strontium chloride (SrCl₂) [11], and barium chloride (BaCl₂) [12, 13]. In [14], 36 solid-gas reactions are collected, mainly chlorides reacting with ammonia. On the contrary, composite materials are a mixture between the aforementioned salts and inert solids such as expanded graphite [15], vermiculite [16, 17], carbon Sibunit, and Al₂O₃ [17]. Composite materials have also been used as solid sorbents as they increase thermal conductivity and prevent agglomeration and swelling. The disadvantage of these composite materials is the low solid-gas ratio and the reactor’s increased volume due to the addition of new materials, compared to the use of pure salts. On the contrary, the predominant refrigerant in solid-gas systems is ammonia due to its vaporization properties, high availability, and minimal environmental impact.

A thermodynamic solid-gas sorption cycle is composed primarily of four elements: an evaporator, an absorber/generator, a condenser, and an expansion device. This cycle works intermittently; when the absorber/generator generates, it cannot absorb and vice versa. The reactor is the device where the phenomena of heat transfer occur by mass and kinetic reactions. In [18], the phenomena of heat and mass transfer are modeled, as well as the kinetic reaction by means of differential equation systems, in order to have a better knowledge in solid-gas reactors.

Mathematical models for studying solid-gas reactors are classified as local, global, and analytical. Local models deal with the phenomena of heat and mass transfer and are numerically solved. Global models deal with thermal conductivity, specific heat, and permeability averages. Analytical models consider reaction times and geometrical configurations [19]. There are also approaches where phenomena within the reactor have been analyzed at different scales [20] and different materials are used as sorbents [21]. It is also reported that the kinetic reaction in thermochemical reactors is complicated and has not been precisely modeled yet [6].

As it has been mentioned before, very complete studies on the phenomena of heat and mass transfer in solid-gas systems in numerical and analytical modalities with different salts have been carried out. There are also reports on intermittent cooling systems that use artificial neural networks (ANNs) to determine the performance in thermochemical reactors during the solid-gas sorption process [22, 23] and show that they are a promising technique in this area [24]. In this work, artificial neural networks and mathematical symmetry groups were used to model a thermochemical reactor of a solid-gas cooling system, where the solid is barium chloride (BaCl₂) and the refrigerant is ammonia (NH₃).

A thermochemical reactor working with the BaCl₂-NH₃ pair was designed and built. Experimental data were collected from this reactor in order to obtain a mathematical prediction based on neural networks that model such a reactor. This provides an artificial intelligence tool that contributes to the design of reactors used in solid-gas sorption cooling.

The rest of this work is organized as follows: Section 2 describes the experimental setup and the influence of the heating fluid on the thermochemical reactor. Section 3 shows the modeling of the thermochemical reactor using neural networks and mathematical symmetry groups. Section 4 shows the obtained results. Section 5 shows the most relevant variables in ANN learning. Then, Section 6 shows the results of the uncertainty analysis followed by the main conclusions and a list of references.

2. System’s Description and Experiments

Figure 1 shows the intermittent solid-gas sorption refrigeration cycle used in this work; this is composed as follows: the thermochemical reactor that acts as a generator/absorber, the condenser, the refrigerant storage tank, an expansion valve, and the evaporator. The system was built in stainless steel [12].

Figure 1: Thermochemical refrigerator system: 1, generator/absorber, 2, condenser, 3, ammonia condensate tank, 4, expansion valve, 5, evaporator, 6, auxiliary systems for heat and cold water supply, 7, data acquisition system, and 8, ammonia supply tank [12].

In addition, the following subsystems were built: (1) water supply to heat/cool the thermochemical reactor and to remove heat from the condenser. Both temperatures are controlled by this subsystem consisting of a thermal bath (heat), a vapor compression cooler (cool), and a water pumping system (remove heat); (2) ammonia supply to the thermodynamic cycle; and (3) a data acquisition system.

Figure 2 depicts an intermittent cooling system. The system works in two stages: (a) generation-condensation and (b) evaporation-sorption. In the former, the VA and VC ball valves are closed, and the VB ball valve is open. Subsequently, the generator/absorber containing the BaCl₂-NH₃ pair is heated (QG) to a temperature where the refrigerant (NH₃) is desorbed. The refrigerant in the vapor state passes to the condenser where heat (QC) is extracted, becomes a liquid, and is stored in the condensate tank.

Figure 2: Solid-gas sorption refrigeration cycle.

In the latter, VB closes, and VA and VC open. Subsequently, the expansion valve opens and regulates the refrigerant to enter the evaporator where it absorbs heat from the medium to be cooled (QE). The refrigerant then returns to the generator/absorber where heat (QA) is removed during the absorption process.

Three J-type thermocouples with different radii were installed in the generator/absorber, taking the reactor centre as reference: r1 = 0 cm, r2 = 2 cm, and r3 = 4 cm. This was done in order to obtain the temperatures Tr1, Tr2, and Tr3 at the respective radii, when heat (QG) is added to the reactor by heating water at various temperatures (Thf). Thermocouples were calibrated with high-precision Ametek brand electronic equipment built by Jofra Instruments model D55SE and had a temperature range of 0 to 150°C. A Cole-Parmer pressure sensor with an operating range of 0 to 20.68 bar (g) and 4 to 20 mA output was also installed to measure the generated pressure (Pg). To measure the volume of condensed or desorbed ammonia (Vd), an external glass meter of 38.1 mm (1½ ”) diameter was installed, coupled to their respective valves as can be seen in Figure 1 (3). To measure the amount of desorbed ammonia (accumulated volume), readings were taken every 5 minutes from the start of the condensation until the desorption process was completed.

2.1. Influence of the Heating Fluid on the Reactor

The desorption process was carried out in order to determine the behavior of the temperature distribution within the reactor and quantify the volume of the desorbed refrigerant. To perform the desorption, the reactor was heated with water at different Thf temperatures (70°C, 80°C, and 95° C) and at a condensation temperature of 23°C. The variables recorded were temperatures Tr1, Tr2, and Tr3, generated pressure Pg, and the volume of desorbed ammonia in its liquid state (Vd).

Temperature distributions within the reactor (Tr1, Tr2, and Tr3) were obtained; however, only the result of Tr2 is shown (Figure 3(a)) as an illustrative case of the influence of Thf on the reactor’s radial temperatures. A linear behavior can be observed in Figure 3, that is, the higher the temperature of Thf, the higher the temperature of Tr2, and the same happens with the variable Vg (Figure 3(b)). The pressure Pg (Figure 3(c)) does not show an increasing behavior; this is because the pressure during the desorption process (phase change) must behave as an isobaric process.

(a)

(b)

(c)

Figure 3: Influence of the heating fluid on the reactor parameters. (a) Influence of Thf (70°C, 80°C, and 95°C) on Tr2. (b) Influence of Thf (70°C, 80°C, and 95°C) on Vg. (c) Influence of Thf (70°C, 80°C, and 95°C) on Pg.

3. Neural Network Proposal for Modeling the Thermochemical Reactor

An artificial neural network (ANN) is a set of interconnected nodes, trying to mimic the functioning of the human brain. An ANN is a function . Each of the nodes has an input function and an output function. Once the ANN is designed and trained, then given an input dataset, the ANN can predict the values of the output variables. ANNs are used especially to model the behavior of nonlinear systems [24]. To train the ANN, the Levenberg–Marquardt algorithm [25, 26] was used because of its convergence speed and precision.

Figures 4 and 5 show the network structure and its input and output variables, which was implemented in MATLAB. Table 1 shows the definition of the neural network’s input and output variables. The ANN’s precision was determined by using the root-mean-square error (RMSE).

Figure 4: Numerical procedure used in the ANN learning process and the iterative architecture used by the ANN model to estimate temperatures, pressures, and ammonia production, respectively.

Figure 5: Optimized architecture of the ANN model.

Table 1: Working range of the thermochemical reactor by barium chloride-ammonia and variables used for the ANN training.

Parameter
Nomenclature
Minimum
Maximum
Unit

Input

Radius 1 inside the reactor
r1
0
0
cm

Radius 2 inside the reactor
r2
2
2
cm

Radius 3 inside the reactor
r3
4
4
cm

Heating fluid’s temperature
Thf
22.512
95.065
°C

Volume of the solid (BaCl₂)
Vs
0.0002
0.0002
m³

Refrigerant volume (NH₃)
Vr
0.0003
0.0003
m³

Desorption time
t
0
12.4375
hr

Reactor height
H
0.1702
0.1702
m

Reactor nominal diameter
Dn
0.0732
0.073 2
m

Output

Temperature at radius 1
Tr1
23.162
75.309
°C

Temperature at radius 2
Tr2
22.866
73.992
°C

Temperature at radius 3
Tr3
22.772
75.437
°C

Volume of the refrigerant
Vg
0
110
ml

Pressure into the reactor
Pg
2.572
11.834
bar

For the ANN training, 2 databases were used. The first database has 7637 experimental data and was generated with the Thf fluid at 70°C. The second database has 8956 experimental data and was generated with the Thf fluid at 95°C. The neural network was trained with these two merged databases. For a better network training, the databases were merged obtaining a single database of 16,593 records, and an alternating group was applied to this database [27, 28]. Then, tuples or permutations of the base 2 multibase mathematical symmetry group (obtained from the multibase symmetry group) were applied to the resulting database. The multibase mathematical symmetry groups are obtained as follows [29].

Let G be a finite group of order N. Rename the elements of G by the set . Let be a number that will indicate a numerical base (for example, base 10 and base 2). Then, can be represented as , where is a cyclic group of order and .

On the contrary, each number in the set can be represented as

where y.

Each element in is given by (vector form of the number), where .

Let be the operation of group and the operation of the cyclic group . If , then

where the cyclical group operation is represented as

The inverse of the element is defined as

With the previous group operations, the matrix is formed, where is the number of columns or records in the database used to train the ANN, and each row of this matrix represents a permutation of this database. The matrix is defined as follows:

On the contrary, to determine the optimal ANN architecture, a set of various network configurations was tested, each one with different numbers of hidden layers and different numbers of neurons in each layer. It was found that the simplest and most accurate network model was 9-9-5 (nine neurons in the hidden layer and 5 neurons in the output layer) as shown in Figure 6. The transfer functions used in the network training were tansig function (hidden layer) and pureline function (output layer).

Figure 6: The effect of the neuron number of the hidden layer on the correlation coefficient “R” and the mean square error “MSE” for the ANN model.

According to Figure 6, it was concluded that the best ANN model was achieved with 09 neurons on the hidden layer. However, the performance in the training set is good with more than 09 neurons, but the performance of the test set is significantly worse, which could indicate an overfit. Therefore, in order to improve the results, the number of neurons in the hidden layer has been reduced to 09 neurons [30]. Consequently, the optimal network configuration was approximately 09-09-05. This model presents smaller values of RMSE (0.4%) and higher R² (99.97%).

Once the ANN architecture was defined, each output of the neural network is calculated by the following equation:

[/lcmath]

Explicitly, each output of the neural network is calculated by the following equations (7)–(11):

where N = total neurons of the hidden layer, M = total input, LWij = weights of the output layer from neuron j of the hidden layer to the output, IWjk = weights of the input layer from input k to neuron j of the hidden layer, Ink = k-th input, b1j = bias of the hidden layer of the j-th neuron, b2i = bias of the output layer at the i-th output, and I = i-th output.

4. Results and Discussion

After training the neural network, the results were obtained as follows (see Table 2). The ANN was trained with 75% of data, and the other 25% of data were used for validation and testing (see Section 3).

Table 2: Weights and bias of the neural network.

Input layer

r1
r2
r3
Thf
Vs
Vg
t
H
Dn

Neuron 1
3.26E + 00
1.63E + 00
0.00E + 00
4.81E − 02
1.40E − 04
2.26E − 04
−9.14E − 01
1.39E − 01
5.96E − 02

Neuron 2
−1.71E + 00
−8.56E − 01
0.00E + 00
1.48E − 01
−7.36E − 05
−1.19E − 04
7.18E − 01
−7.29E − 02
−3.13E − 02

Neuron 3
4.29E − 01
2.15E − 01
0.00E + 00
5.96E − 02
1.85E − 05
2.97E − 05
5.72E − 01
1.83E − 02
7.85E − 03

Neuron 4
3.90E + 00
1.95E + 00
0.00E + 00
−1.75E − 01
1.68E − 04
2.70E − 04
−1.39E − 01
1.66E − 01
7.13E − 02

Neuron 5
−4.39E − 01
−2.19E−01
0.00E + 00
1.65E − 02
−1.89E − 05
−3.04E − 05
−5.77E − 03
−1.87E − 02
−8.03E − 03

Neuron 6
3.60E + 00
1.80E + 00
0.00E + 00
−3.76E − 01
1.55E − 04
2.49E − 04
1.59E + 00
1.53E − 01
6.58E − 02

Neuron 7
3.03E − 01
1.52E − 01
0.00E + 00
−5.32E − 03
1.30E − 05
2.10E − 05
−1.18E − 01
1.29E − 02
5.55E − 03

Neuron 8
−6.36E − 02
−3.18E − 02
0.00E + 00
−8.44E − 02
−2.74E − 06
−4.41E − 06
1.37E − 01
−2.71E − 03
−1.16E − 03

Neuron 9
2.76E + 00
1.38E + 00
0.00E + 00
−2.16E − 01
1.19E − 04
1.91E − 04
−2.28E + 00
1.17E − 01
5.05E − 02

Hidden layer

Bias

1
2
3
4
5
6
7
8
9

b1
−1.13E + 01
−5.76E + 00
−9.33E + 00
−4.22E + 00
9.58E−01
−5.05E + 00
−1.03E + 00
4.56E + 00
4.26E + 00

Output layer

Inputs

1
2
3
4
5
6
7
8
9

Neuron 1
−5.86E − 02
4.22E + 00
−3.83E + 00
9.03E − 01
5.15E + 01
6.45E − 01
−4.80E + 00
−2.78E + 00
−1.64E + 00

Neuron 2
−4.09E − 01
3.27E + 00
−3.77E + 00
8.74E − 01
5.41E + 01
8.69E − 01
−1.03E + 01
−7.84E − 01
2.38E − 01

Neuron 3
−1.10E + 01
1.05E + 00
−9.41E − 03
8.09E − 01
−3.78E + 00
−3.65E + 00
−4.87E + 00
−1.28E + 01
−1.08E − 01

Neuron 4
3.88E + 00
−2.67E + 01
3.95E + 01
−3.04E + 01
−5.10E + 01
9.09E−02
−6.28E + 01
−4.04E + 00
−7.94E + 00

Neuron 5
8.20E − 02
4.34E − 01
−2.80E − 01
1.70E + 00
5.02E + 00
3.54E − 01
−2.67E + 00
−2.41E + 00
2.14E − 01

Output layer

Bias

b2
6.19E + 01
6.08E + 01
4.83E + 01
3.76E + 01
6.93E + 00

In order to verify learning of the network, a database with 7,420 experimental data was used. These data were the heating time (t) from 0 to 10.31 h, the heating temperature (Thf) from 24.7°C to an average maximum temperature of 81.2°C, and the reactor geometry values which were constant, as well as the amounts of NH₃ and BaCl₂ (see Table 1).

By training the network with the appended database and verifying its learning, a correlation coefficient of r = 0.9671 was obtained (see Section 3). Applying an alternating group to the ANN training database and checking its learning with this new database, a correlation coefficient of r = 0.9920 was obtained. Finally, applying an alternating group and then applying elements from the multibase mathematical symmetry group to the ANN training database and verifying their learning with the new database, a correlation coefficient of r = 0.9957 was obtained. For a fair comparison, in all ANN trainings with different permutations (alternating group and multibase mathematical symmetry group) from the training database, the network weights were initialized with a random seed of 1230.

Figure 7(a) shows that the behavior of simulated temperature Tr3 during the addition of heat to BaCl₂-NH₃ (sensible heat), against the experimental data, coincides with a correlation coefficient of r = 0.9996. During the phase change or desorption, the difference of the simulated data against the experimental data becomes minimal.

(a)

(b)

(c)

(d)

(e)

Figure 7: Testing the neural network (correlation coefficient r = 0.9957): (a) temperature near the reactor wall (Tr3); correlation coefficient r = 0.9996, (b) temperature at the midpoint between the centre and the reactor wall (Tr2); correlation coefficient r = 0.9924, (c) temperature at the reactor’s centre (Tr1); correlation coefficient r = 0.9953, (d) ammonia produced; desorption volume (Vg); correlation coefficient r = 0.99, and (e) pressure into the reactor (Pg); correlation coefficient r = 0.9503.

In the Tr2 case (Figure 7(b)), it is observed that the simulated and experimental Tr2 temperatures are consistent and have a correlation coefficient of r = 0.9924. It should be noted that, during desorption, there are some minimal differences between simulated and experimental values; however, from the experimental point of view, this is acceptable. The evolution of the Tr1 simulated temperature (Figure 7(c)) coincides with a correlation coefficient of r = 0.9953 which is acceptable from a practical point of view.

The simulated values of the refrigerant volume versus the experimental ones coincide with a correlation coefficient of r = 0.99 (Figure 7(d)). It should be noted that the experimental data were visually recorded every five minutes. However, it has an approximate value to the recorded final volume at the experimental desorption.

The simulated pressure inside the reactor (Figure 7(e)) is consistent with the sensible heating process, as it occurred with temperatures at different radii, and it has a correlation coefficient of r = 0.9503. During desorption, the pressure must have a constant behavior; however, in the experimental stage, there is a gradual increase. In the case of simulated pressure, the same increase occurs, with a maximum difference of 1 bar; it must be taken into account that the pressure values in Figure 3(c) do not show an ascending or descending pattern when increasing Thf.

5. Sensitivity Analysis

In order to determine the impact of each input variable on the outputs estimated by the neural network, a sensitivity analysis was developed using a genetic programming algorithm executed in Eureqa [31]. The sensitivity analysis showed that the most used characteristic is the heating fluid’s temperature (35%) followed by the desorption time (20%), that is, these variables are the most important in the ANN’s learning process.

The sensitivity analysis for each of the input variables was also performed by generalizing the equations of Dimopoulos et al. [32] for outputs, which is given as follows:

where is the sensitivity of the q-th output of the neural network, for , N is the total input data to the neural network, and is the total number of outputs of the network.

where is the total input parameters to the neural network.

where ni is the number of neurons of the hidden layer, , is the q-th activation function of the output layer, is the -th activation function of the hidden layer, are the weights of the output layer, are the weights of the hidden layer, is the weight from the entrance to the hidden layer, from the i-th network entrance to the e-th neuron of the hidden layer, is the weight that goes from the hidden layer to the output layer, from the i-th neuron of the hidden layer to the q-th neuron of the output layer, and and are those derived from the activation functions of the neurons of the output layer and the hidden layer, respectively.

In equation (14), the domains of the activation functions are defined as follows:

where is the p-th input to train or simulate the network.

where is the domain of the activation function of the q-th neuron of the output layer.

Equations (12)–(16) were implemented in MATLAB yielding Table 3, as in Eureqa; it is concluded that the most important variables for the learning process of the neural network are Thf and t.

Table 3: Result of sensitivity analysis of the input parameters of the neural network (all data are multiplied by 1e + 57).

Output
Input

r1
r2
r3
Thf
Vs
Vg
t
H
Dn

Tr3
0
0
0
0.01442
0
0
0.000000000000006
0
0

Tr2
0
0
0
1.971712
0
0
0.000000000000185
0
0

Tr1
0
0
0
0.163830
0
0
0.000000000000321
0
0

Vd
0
0
0
0.576825
0
0
0.000000000001319
0
0

Pg
0
0
0
0.202679
0
0
0.000000000000003
0
0

6. Uncertainty Analysis

In order to verify the accuracy improvement of the ANN, the uncertainties of the temperature sensors, pressure sensor, and desorption were analyzed. Some methodologies and standards for obtaining uncertainties are mentioned in [33–35]. All the uncertainties calculated for the sensors used in this research were based on evaluation of measurements described in [35], and a probability and a coverage factor k = 3 were proposed. The expanded uncertainty of the J-type class 2 temperature sensor of Ametek brand [36] is U = 2.9°C. The expanded uncertainty of the pressure sensor of Cole-Parmer brand [37] is U = 0.18 bar. On the contrary, the desorption of the refrigerant was recorded manually, so the fuzzy set type 2 technique [34] and the central limit theorem were used to obtain the expanded uncertainty, obtaining U = 4.47 ml. Using these uncertainties, the correlation coefficient of r = 0.9957 was obtained between the experimental data (those never used for the network training) and the data simulated by the neural network.

7. Conclusions

The proposed artificial neural network allowed modeling the behavior of a BaCl₂-NH₃ reactor globally with a precision of 0.9957 without using complex systems of differential equations. However, it is still necessary to improve the learning precision and the architecture of this neural network to model the behavior of these reactors. For example, it is necessary to improve the output of the neural network that simulates pressure, which has a correlation coefficient of r = 0.9503.

Applying an alternating group of mathematical symmetry to the input data of the neural network allowed obtaining better precision in the learning process than only appending the 70°C and 95°C databases. However, when applying the permutations of the multibase mathematical symmetry groups to the input data of the neural network, a better precision (r = 0.9957) was obtained in the mentioned learning process than applying the alternating symmetry group .

By applying the permutations of the mathematical symmetry group to the neural network input database, the training algorithm (Levenberg– Marquardt) converges faster and with greater precision compared to the appended database. However, given the high computational complexity of the permutations generated for the mathematical symmetry group, not all permutations were applied to the input data of the neural network.

As shown by the sensitivity analysis, the most important variables for the neural network training are heating fluid temperature (Thf) and desorption time (t). The analysis of the uncertainties of experimental measurements allowed obtaining a higher learning precision of the neural network (r = 0.9957).

Data Availability

The experimental data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

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

Acknowledgments

Special thanks are due to the Institute of Renewable Energies (IER) for their valuable help in the development of this work and to the Secretary of Public Education (SEP) for granting the National Conventional Scholarship for doctoral studies PROMEP/103.5/09/4363.

References

L. Wang, X. Bu, and W. Ma, “Experimental study of an adsorption refrigeration test unit,” Energy Procedia, vol. 152, pp. 895–903, 2018.

H. Z. Hassan, “Energy analysis and performance evaluation of the adsorption refrigeration system,” ISRN Mechanical Engineering, vol. 2013, Article ID 704340, 14 pages, 2013.

F. Martínez-Tejeda, I. Pilatowsky, R. Best et al., “Experimental barium chloride-ammonia cooling cycle study at low generation temperatures,” Applied Thermal Engineering, vol. 141, pp. 751–761, 2018.

K. R. Ullah, R. Saidur, H. W. Ping, R. K. Akikur, and N. H. Shuvo, “A review of solar thermal refrigeration and cooling methods,” Renewable and Sustainable Energy Reviews, vol. 24, pp. 499–513, 2013.

R. Novella, V. Dolz, J. Martín, and L. Royo-Pascual, “Thermodynamic analysis of an absorption refrigeration system used to cool down the intake air in an Internal Combustion Engine,” Applied Thermal Engineering, vol. 111, pp. 257–270, 2017.

G. L. An, L. W. Wang, J. Gao, and R. Z. Wang, “A review on the solid sorption mechanism and kinetic models of metal halide-ammonia working pairs,” Renewable and Sustainable Energy Reviews, vol. 91, pp. 783–792, 2018.

F. Meunier, “Adsorption heat powered heat pumps,” Applied Thermal Engineering, vol. 61, no. 2, pp. 830–836, 2013.

S. O. Enibe and O. C. Iloejei, “Transient analysis and performance prediction of a solid absorption solar refrigerator,” Journal of Solar Energy Engineering, vol. 61, no. 1, pp. 43–59, 1997.

L. L. Vasiliev, D. A. Mishkinis, A. A. Antukh, and L. L. Vasiliev, “Solar-gas sorption refrigerator,” Adsorption, vol. 7, no. 2, pp. 149–161, 2001.

R. G. Oliveira and R. Z. Wang, “A consolidated calcium chloride-expanded graphite compound for use in sorption refrigeration systems,” Carbon, vol. 45, no. 2, pp. 390–396, 2007.

K. NagaMalleswara Rao, M. Ram Gopal, and S. Bhattacharyya, “Analysis of a SrCl2-NH3solid sorption refrigeration system,” International Journal of Low-Carbon Technologies, vol. 10, no. 4, pp. 365–373, 2013.

C. Rivera, I. Pilatowsky, E. Méndez, and W. Rivera, “Experimental study of a thermo-chemical refrigerator using the barium chloride-ammonia reaction,” International Journal of Hydrogen Energy, vol. 32, no. 15, pp. 3154–3158, 2007.

C. Dueñas, I. Pilatowsky, R. J. Romero, A. Oskam, and P. A. Finck, “Dynamic study of the thermal behaviour of solar thermochemical refrigerator: barium chloride-ammonia for ice production,” Solar Energy Materials and Solar Cells, vol. 70, no. 3, pp. 401–413, 2001.

P. Neveu and J. Castaing, “Solid-gas chemical heat pumps: field of application and performance of the internal heat of reaction recovery process,” Heat Recovery Systems and CHP, vol. 13, no. 3, pp. 233–251, 1993.

T. X. Li, R. Z. Wang, L. W. Wang, and J. K. Kiplagat, ““Performance study of a consolidated manganese chloride–expanded graphite compound for sorption deep-freezing processes”,” Journal Applied Energy, vol. 86, no. 7-8, pp. 1201–1209, 2009.

Y. Zhong, R. E. Critoph, R. N. Thorpe et al., ““Isothermal sorption characteristics of the BaCl2–NH3 pair in a vermiculite host matrix”,” Journal of Thermal Science and Engineering, vol. 27, no. 14-15, pp. 2455–2462, 2007.

J. V. Veselovskaya, M. M. Tokarev, and Y. I. Aristov, “Novel ammonia sorbents porous matrix modified by active salt for adsorptive heat transformation 1. Barium chloride in various matrices,” Journal of Thermal Science and Engineering, vol. 30, no. 6-7, pp. 584–589, 2010.

H.-B. Lu, N. Mazet, and B. Spinner, “Modelling of gas-solid reaction-coupling of heat and mass transfer with chemical reaction,” Chemical Engineering Science, vol. 51, no. 15, pp. 3829–3845, 1996.

D. Stitou and G. Crozat, “Dimensioning nomograms for the design of fixed-bed solid-gas thermochemical reactors with various geometrical configurations,” Chemical Engineering and Processing: Process Intensification, vol. 36, no. 1, pp. 45–58, 1997.

P. Neveu and J. Castaing-Lasvignottes, “Development of a numerical sizing tool applied to a solid-gas thermochemical transformer-II. Influence of external couplings on the dynamic behaviour of a solid-gas thermochemical transformer,” Journal of Thermal Science and Engineering, vol. 17, no. 6, pp. 519–536, 1997.

H.-B. Lu and N. Mazet, “Mass-transfer parameters in gas-solid reactive media to identify permeability of IMPEX,” AIChE Journal, vol. 45, no. 11, pp. 2444–2453, 1999.

B. A. Escobedo-Trujillo, D. Colorado, W. Rivera, and F. A. Alaffita-Hernández, “Neural network and polynomial model to improve the coefficient of performance prediction for solar intermittent refrigeration system,” Solar Energy, vol. 129, pp. 28–37, 2016.

V. Baiju and C. Muraleedharan, “Exergy assessment of single stage solar adsorption refrigeration system using ANN,” ISRN Mechanical Engineering, vol. 2012, Article ID 915154, 2012.

Y. Liu, N. Dinh, Y. Sato, and B. Niceno, “Data-driven modeling for boiling heat transfer: using deep neural networks and high-fidelity simulation results,” Applied Thermal Engineering, vol. 144, pp. 305–320, 2018.

R. Roohia, M. Jafarib, E. Jahantab et al., “Application of artificial neural network model for the identification the effect of municipal waste compost and biochar on phytoremediation of contaminated soils,” Journal of Geochemical Exploration, vol. 208, Article ID 106399, 2019.

M. T. Hagan and M. B. Menhaj, “Training feedforward networks with the Marquardt algorithm,” IEEE Transactions on Neural Networks, vol. 5, no. 6, pp. 989–993, 1994.

Oriol Vinyals, S. Bengio, and M. Kudlur, “Order matters: sequence to sequence for sets,” 2016, Published as a conference paper at ICLR.

A. G. Kurosh, The Theory of Group, Chelsea Publishing Company, New York, NY, USA, 2nd edition, 1960.

E. Shelomov, J. Seguel, A. Ortiz, J. Ruz, and G. Shelomova, “Fourier analysis in commutative finite groups with matlab,” in Proceedings of the Memories of the VII Conference of Electrical Engineering CINVESTAV, pp. 90–116, Mexico city, Mexico, November 2002.

K. Gnana Sheela and S. N. Deepa, “Review on methods to fix number of hidden neurons in neural networks,” J. Math. Probl. Eng, vol. 2013, Article ID 425740, 11 pages, 2013.

Eureqa 1.24.0, Nutonian Inc., Cambridge, MA, USA, 2019, http://www.nutonian.com.

Y. Dimopoulos, P. Bourret, and S. Lek, “Use of some sensitivity criteria for choosing networks with good generalization ability,” Neural Processing Letters, vol. 2, no. 6, pp. 1–4, 1995.

Y. Liu, N. T. Dinh, R. C. Smith, and X. Sun, “Uncertainty Quantification of Two-Phase flow and boiling heat transfer simulations through a data-driven modular Bayesian approach,” International Journal of Heat and Mass Transfer, vol. 138, pp. 1096–1116, 2019.

J. E. Moreno, M. A. Sanchez, O. Mendoza et al., “Design of an interval Type-2 fuzzy model with justifiable uncertainty,” Information Sciences, vol. 513, pp. 206–221, 2020.

Guide to the Expression of Uncertainty in Measurement, ISO, Geneva, Switzerland , https://www.bipm.org/en/publications/guides/gum.html.

M. T. Thermocouple Assembly, Data Sheet 1500. Class 2 Type-J, Ametek, Berwyn, PA, USA, 2020, https://www.ametekcalibration.cn/-/media/ametekcalibration/download_links/temperature%20sensors/fp%20sensors/1500%20and%201600%20series/temperature-sensors-1500-series-datasheet-us.pdf.

Cole-Parmer, in High-Accuracy Pressure Transducers , Cole-Parmer, Vernon Hills, IL, USA, 2020, https://www.coleparmer.com/p/cole-parmer-high-accuracy-pressure-transducers/15143.visited.