Performance Assessment of an NH3/LiNO3 Bubble Plate Absorber Applying a Semi-Empirical Model and Artificial Neural Networks

In this study, ammonia vapor absorption with NH3/LiNO3 was assessed using correlations derived from a semi-empirical model, and artificial neural networks (ANNs). The absorption process was studied in an H-type corrugated plate absorber working in bubble mode under the conditions of an absorption chiller machine driven by low-temperature heat sources. The semi-empirical model is based on discretized heat and mass balances, and heat and mass transfer correlations, proposed and developed from experimental data. The ANN model consists of five trained artificial neurons, six inputs (inlet flows and temperatures, solution pressure, and concentration), and three outputs (absorption mass flux, and solution heat and mass transfer coefficients). The semi-empirical model allows estimation of temperatures and concentration along the absorber, in addition to overall heat and mass transfer. Furthermore, the ANN design estimates overall heat and mass transfer without the need for internal details of the absorption phenomenon and thermophysical properties. Results show that the semi-empirical model predicts the absorption mass flux and heat flow with maximum errors of 15.8% and 12.5%, respectively. Maximum errors of the ANN model are 10.8% and 11.3% for the mass flux and thermal load, respectively.


Introduction
Absorption cooling systems appear to be promising energy-saving technologies if they use residual/available heat or solar energy for their operation [1]. However, actual absorption equipment with traditional working mixtures has a number of limitations. To overcome the inconvenience and restrictions of traditional working mixtures, one option suggested in the open literature is the use of NH 3 /LiNO 3 [2,3]. The most attractive aspect of absorption systems using a NH 3 /LiNO 3 mixture is that they can be cooled by air because crystallisation does not occur at high condensation temperatures. In addition, systems with a NH 3 /LiNO 3 mixture do not require rectification of the refrigerant vapor emitted from the generator because the absorbent is a salt, and these systems can operate at lower activation temperatures compared to ammonia/water (NH 3 /H 2 O) systems. Nevertheless, former experimental investigations have shown that the major limiting aspect of the NH 3 /LiNO 3 mixture is its elevated viscosity, which restricts heat and mass transfer, primarily in the absorber, in comparison to the NH 3 /H 2 O mixture [4][5][6][7]. This limitation has motivated studies aimed at enhancing the absorption process with NH 3 /LiNO 3 by employing passive techniques [8][9][10]. A Based on the literature review, no validated semi-empirical heat and mass transfer models exist in the open literature to predict the NH 3 bubble absorption in an H-type corrugated PHE with NH 3 /LiNO 3 mixture. The present study was part of an R&D project investigating the characterization of the absorption process with NH 3 /LiNO 3 . The aim was to improve the information available about the performance of NH 3 /LiNO 3 bubble absorbers and contribute to their technological development. Therefore, this paper aims to provide a description and validation of dimensionless correlations and two different models developed to predict the performance of the NH 3 absorption process in a bubble plate absorber. The first model developed is based on discretized heat and mass balances in addition to empirical heat and mass transfer correlations fitted from experimental data, whereas the second model is based on trained artificial neural networks. Initially, each model is assessed and validated individually. Then, a performance comparison between both models is conducted. Finally, the advantages of each model are presented.

Methodology
This study involved the proposal of correlations generated from the experimental characterization of the absorption process of ammonia vapor, development of two model options to predict the NH 3 absorption process in a bubble plate absorber, validation of both models, and comparison of the results. The key details of the experimental test setup are presented in the next sub-section. The first model is a semi-empirical model detailed in Section 2.2, while the second model is an artificial neural networks-based model described in Section 2.3. Figure 1 shows the experimental setup used to assess and correlate the ammonia vapor absorption in a plate heat exchanger. In this facility, measurable variables, such as inlet flows, inlet temperatures, inlet solution concentration in ammonia, and inlet solution pressure in the test section, can be monitored and controlled. The facility consists of 3 circuits: the NH 3 /LiNO 3 mixture circuit, the heating circuit, and the heat dissipation/cooling circuit. In the NH 3 /LiNO 3 mixture circuit, the solution (weak in refrigerant) is pumped (B1) from the storage tank 1 (ST1) to the heat exchanger (SHE) where it reaches the required solution temperature at the test section inlet. After the absorption process takes place in the test section, the solution (strong in refrigerant) leaves the test section to flow to a vapor-liquid separator (VLS) and finally to the storage tank 2 (ST2). The heating circuit includes a 5 kW heater (R1), a water-water heat exchanger (HX2), a pump (B2), and a magnetic flow meter (F). The heat dissipation/cooling circuit is used to dissipate the excess of heat released during the absorption process in the test section. This circuit includes a 5 kW heater (R2), a magnetic flow meter (F), a pump (B3), and a heat exchanger (HX1). The test section is also connected to an ammonia bottle which supplies the ammonia vapor at a pressure slightly higher than that in the absorber inlet.

Experimental Setup
Various measuring devices were used for data monitoring and collection. Figure 1 shows the location of instrumentation such as sensors for temperature (T) and pressure (P), and magnetic (F) and Coriolis flowmeters (C). The data collection began after the ammonia vapor injection and once the steady-state was reached in the thermal conditions of the flows. Data was collected every 5 s for around 15 min.
The absorber used for this study was a plate heat exchanger (PHE) with four braze plates (three channels). The PHE is made of stainless steel AISI 316 and provided by Alfa Laval. It has Chevron-H type corrugations which correspond to an angle of 30 • from the horizontal axis.
The operation of the absorber is in counter-current flow. The NH 3 /LiNO 3 mixture and ammonia vapor flow upward in the inner channel while the cooling-water flows downward in the side channels. For more details of the PHE, please see Table 1 and Figure 2. More details about the experimental procedure and test facility are given by Amaris in [41].       The nominal operating conditions of the present application are shown in Table 2. These conditions of operation correspond to those of a simple-effect absorption chiller driven by low-temperature heat and a chilled-water temperature of 5 • C.

Semi-Empirical Model
A semi-empirical model was developed to predict the performance of the NH 3 /LiNO 3 bubble plate absorber. Modelling of the absorption process was carried out by the discretization of the PHE in small volumes of control. In each volume of control, energy and mass balances were solved at the same time as the simulation advanced progressively from the inlet to the outlet of the PHE (see Figure 3).
The NH 3 /LiNO 3 thermodynamic properties reported by Libotean et al. [42,43] and Cuenca et al. [44] were implemented in the model. The enthalpy values were determined according to the method proposed by Haltenberger [45] and implemented by McNeely [46].
The main considerations for the development of the model are presented below: The nominal operating conditions of the present application are shown in Table 2. These conditions of operation correspond to those of a simple-effect absorption chiller driven by low-temperature heat and a chilled-water temperature of 5 °C.

Semi-Empirical Model
A semi-empirical model was developed to predict the performance of the NH3/LiNO3 bubble plate absorber. Modelling of the absorption process was carried out by the discretization of the PHE in small volumes of control. In each volume of control, energy and mass balances were solved at the same time as the simulation advanced progressively from the inlet to the outlet of the PHE (see Figure 3).
The NH3/LiNO3 thermodynamic properties reported by Libotean et al. [42,43] and Cuenca et al. [44] were implemented in the model. The enthalpy values were determined according to the method proposed by Haltenberger [45] and implemented by McNeely [46].
The main considerations for the development of the model are presented below: • Model in one dimension and steady-state conditions, • Pressure drop is negligible, • Equal heat and mass transfer areas, • Refrigerant temperature equal to the inlet solution temperature, • Adiabatic exterior border.

Discretized Governing Equations and Correlations
The numerical model relies basically on heat and mass balances, heat and mass transfer correlations, and equations for the thermodynamic properties of the solution and water flows.

Discretized Governing Equations and Correlations
The numerical model relies basically on heat and mass balances, heat and mass transfer correlations, and equations for the thermodynamic properties of the solution and water flows. The heat and mass balances and correlations are described as follows: Solution-side global mass balance: Solution-side partial mass balance: where d • m Sol is the differential mass transfer in each control volume, F AB is the NH 3 absorption mass flux, and dTa is the differential transfer area. dTa is determined as dTa = 2·w·dl, where w is the width of the plate while dl is the control volume differential length. F AB is determined from Equation (5) where k m is the mass transfer coefficient and LMXD is the logarithmic mean concentration difference. k m is determined from Equation (6) where the Sh Sol is approached by Equation (16). In Equation (6), β Sol is the diffusivity of the NH 3 /LiNO 3 which is estimated as in [47].
Solution-side energy balance: The thermal load (d . Q AB ) is approached as follows: The temperature of the water flow at the outlet of each control volume is determined from: The overall heat transfer coefficient (U) is estimated considering the resistances for the heat flow in each volume of control: where α Sol represents the solution heat transfer coefficient and α Cw is the water-side heat transfer coefficient. The heat transfer coefficients are estimated as follows: α Cw = Nu Cw · λ Cw D h_Cw_Side (12) where λ represents the thermal conductivity of each fluid and D h refers to the hydraulic diameter in each side, estimated as 4 × Ap/Pw and simplified to 2 times the plate separation. In this case, the hydraulic diameter is the same for both sides. The solution-side Nusselt number (Nu Sol ) and water-side Nusselt number (Nu Cw ) were estimated using the correlations proposed in Section 3.1.

Model Resolution
The numerical model follows a sequential resolution of heat and mass balances in every control volume beginning from the first control volume to the last.
Because the PHE absorber operates in counter-current configuration, the outlet water-side temperature is initially assumed by the model, while the outlet bulk solution conditions and inlet water-side conditions in each control volume are estimated from the heat and mass balances. The simulation finishes when the estimated inlet water-side temperature reaches the actual inlet water-side temperature. The code of the numerical model was written in the Engineering Equation Solver (EES) software.
The model simulation procedure is as follows: • Input of the absorber geometrical characteristics.

•
Input of inlet operating conditions such as: The model guesses an initial outlet water-side temperature T Cw (n). • Evaluation of the first control volume: Estimation of thermodynamic properties. Estimation of Sh Sol , k m , F AB and mass transfer balances in the control volume. Estimation of α Sol , α Cw , T Sol (n + 1), T Cw (n + 1), and energy balances.
The model guesses T Sol (n + 1), carries out energy balances, and calculates T Cw (n + 1). If balances are close, the simulation runs on the next control volume, else, a new T Sol (n + 1) is assumed. T Sol (n + 1) is set by applying the secant method.

•
When the simulation reaches the last control volume, the calculated T Cw (n + 1) is compared to the actual T Cw . The simulation finishes if values are equal, else, a new T Cw (n) is guessed, and the simulation starts again from the first control volume. T Cw (n) is estimated by applying the secant method.

Artificial Neural Networks Model
The ANN requires input and output information so that training and therefore, the learning process, can take place. The ANN consists of 3 layers: an input layer, a hidden layer, and an output layer. The hidden layer is the heart of the ANN and consists of processing units called neurons. These neurons are grouped and linked to other layers by interconnections called weights. The input layer collects the input data while the output layer reproduces the required data from the network as an output. The required neurons in the hidden layer are defined by the prediction accuracy of the required parameters. Consequently, there are no explicit guidelines to define the design of the ANN structure. The network responses to the input data are the transfer functions at the outlet of each neuron. The hyperbolic tangent sigmoid function (tansig) is the most commonly used transfer function (Equation (13)) in the hidden layer while the linear transfer function (purelin) is used in the output layer (Equation (14)): Energies 2020, 13, 4313 8 of 20 Normalization of the input data was required to avoid any influence of the magnitude of input values in the defined ranges of the experimentation. Therefore, inputs were normalized within the range from 0.1 to 0.9 by Equation (15) [48]: The neurons in the input and output layers correspond to the input and output parameters. The back-propagation method employing the Levenberg-Marquardt optimization algorithm was used to compare the experimental data feeding the network and the target output. Obtaining an accurate and solid ANN depends on the input parameters, the information available, and the method employed for training [48]. An experimental set of data was previously collected from the absorption process of ammonia vapor in an H-type plate absorber with NH 3 /LiNO 3 [41].
The inlet parameters of the PHE, P Sol , T Cw , V Cw , m Sol , T Sol , and x Sol , were set as input data to the ANN. Previously, experimental data collection was conducted at a x Sol of 0.45 with a total of 30 experiments according to the procedure described in Section 2.1. From the available information, 17 experiments corresponded to a heat dissipation temperature of 35 • C, and 13 experiments to a heat dissipation temperature of 40 • C. In every set of experiments, some parameters were kept constant while others were varied. The following information was employed for the design of the ANN: Inlet values from 0.01389 to 0.00289 kg·s −1 forṁ Sol , from 0.450 to 0.1312 m 3 ·h −1 for V Cw , and from 40.8 to 34.9 • C for T Cw . The ANN design was developed using the neural network application of MATLAB software.

Results and Discussion
This section presents the correlations proposed, and results from the semi-empirical and ANN models developed to quantify the performance parameters of NH 3 bubble absorption with a NH 3 /LiNO 3 mixture. The validation and comparison of both models with experimental data are also reported.

Heat and Mass Transfer Correlations
This sub-section presents the heat and mass transfer correlations developed and proposed considering the data collected by Amaris [41] from an experimental study of the H-type PHE absorber. Correlations in Equations (16) and (17) correspond to the solution Nusselt (Nu Sol ) and Sherwood (Sh Sol ) numbers that characterize the absorption process. Validations of these correlations are presented in Figure 4. Moreover, correlation in Equation (18) refers to the water-side Nusselt number (Nu Cw ) also obtained from the performance of the H-type PHE using water in both the inner and side channels. This means that prior to the experimental investigation of the ammonia absorption process, a preliminary heat transfer experimental study was conducted using water in the inner and side channels. Results from this study are reported Equation (18) and Figure 5.
The solution Nusselt, Sherwood, and water-side Nusselt numbers were correlated by reducing the mean square error (MSE) when contrasting the correlated and experimentally estimated dimensionless parameters. Variables such as the inlet heat dissipation water temperature, inlet solution temperature, and solution equilibrium temperature were found to have a significant effect on the heat and mass transfer. Therefore, they were included in the correlation to account for these effects [49].
Nu Sol = 1.24 · Re 0.6 Sol · Pr 1/3 Sol · T Sol,Eq,in T Sol,in Energies 2020, 13, 4313 9 of 20 Figure 4b demonstrates that Equation (16) can adequately estimate 90% of the solution mass transfer coefficient points with relative differences below 15%. Moreover, Equation (17) can estimate 90% of the solution heat transfer coefficient data with relative differences below 15% (see Figure 4a). These correlations estimate the experimental data with average relative differences of 6.6% and 7.3%, for the heat and mass transfer coefficients, respectively. The fit parameter R 2 for each correlation is shown in Figure 4. The correlations proposed in Equations (16) and (17) were developed for Reynolds numbers ranging from 10 to 55. In this study, the Prandtl number ranged between 46 and 51, and the Schmidt number in Equation (16) ranged from 5194 to 5947. Figure 5a,b shows the overall heat transfer coefficient and water-side Nusselt number (Nu Cw ) obtained from the experimental study in the H-type PHE using water in both channels in counter-current configuration and correlated in Equation (15). The black line represents the trend of the results. In these experiments, temperatures were changed in a range from 30 to 60 • C while the mass flow was varied from 50 to 352 kg·h −1 . The correlation proposed in Equation (18) is recommended for a water turbulent flow with a Reynolds number ranging from 400 to 1400. The Prandtl number for the water flow ranged between 4.20 and 4.75. Figure 5a shows that Equation (18) can satisfactorily estimate 82.4% of the overall water heat transfer coefficient data with relative differences below 2% and an average relative difference of 1.17%.   Figure 5a,b show the overall heat transfer coefficient and water-side Nusselt number (NuCw) obtained from the experimental study in the H-type PHE using water in both channels in counter-current configuration and correlated in Equation (15). The black line represents the trend of the results. In these experiments, temperatures were changed in a range from 30 to 60 °C while the mass flow was varied from 50 to 352 kg.h −1 . The correlation proposed in Equation (18) is recommended for a water turbulent flow with a Reynolds number ranging from 400 to 1400. The Prandtl number for the water flow ranged between 4.20 and 4.75. Figure 5a shows that Equation (18) can satisfactorily estimate 82.4% of the overall water heat transfer coefficient data with relative differences below 2% and an average relative difference of 1.17%.   Figure 5a,b show the overall heat transfer coefficient and water-side Nusselt number (NuCw) obtained from the experimental study in the H-type PHE using water in both channels in counter-current configuration and correlated in Equation (15). The black line represents the trend of the results. In these experiments, temperatures were changed in a range from 30 to 60 °C while the mass flow was varied from 50 to 352 kg.h −1 . The correlation proposed in Equation (18) is recommended for a water turbulent flow with a Reynolds number ranging from 400 to 1400. The Prandtl number for the water flow ranged between 4.20 and 4.75. Figure 5a shows that Equation (18) can satisfactorily estimate 82.4% of the overall water heat transfer coefficient data with relative differences below 2% and an average relative difference of 1.17%.

Semi-Empirical Model Results
This subsection presents the estimation of the temperatures and concentration at the inlet and outlet of each partial control volume of the PHE using the semi-empirical model (Mod). Modelling of the plate absorber was performed at a mass flow of 40 kg.h −1 in the solution side. Inlet operating conditions were set as in Table 2. For comparison purposes, measured inlet and outlet conditions for temperature and solution concentration in the absorber are also depicted.
Simulation results were also contrasted to those from the theoretical model described by Cerezo et al. [14] (Mod [14]) at the same conditions. The model described in the literature considers the interface area between the NH3/H2O and the ammonia bubble vapour being absorbed. In addition, it

Semi-Empirical Model Results
This subsection presents the estimation of the temperatures and concentration at the inlet and outlet of each partial control volume of the PHE using the semi-empirical model (Mod). Modelling of the plate absorber was performed at a mass flow of 40 kg·h −1 in the solution side. Inlet operating Energies 2020, 13, 4313 10 of 20 conditions were set as in Table 2. For comparison purposes, measured inlet and outlet conditions for temperature and solution concentration in the absorber are also depicted.
Simulation results were also contrasted to those from the theoretical model described by Cerezo et al. [14] (Mod [14]) at the same conditions. The model described in the literature considers the interface area between the NH 3 /H 2 O and the ammonia bubble vapour being absorbed. In addition, it employs several correlations reported in the literature from configurations where separated heat and mass transfer processes were considered. Moreover, it also included the correlations reported by Herbine and Perez-Blanco [15] to estimate the heat transfer between the NH 3 /H 2 O mixture and the wall. Figure 6 illustrates the profiles of temperature and concentration along the PHE absorber. Results show that the measured and simulated outlet temperatures and concentrations in the PHE were reasonably close. For instance, Figure 6a shows that at a heat dissipation temperature of 40 • C at the PHE inlet, the temperature of the solution at the PHE outlet using the present model is 0.5 • C lower than that of the measured value, whereas the outlet temperature obtained with the model described by Cerezo et al. [14] is 2.2 • C lower. Moreover, at a water temperature equal to 35 • C, Figure 6d shows an immediate drop in the hot-side temperature when employing the semi-empirical model because of the greater potential for heat transfer due to the larger temperature difference. In this case, the solution-side temperature at the PHE outlet using the present model is 2.3 • C lower than that of the measured value, whereas the outlet temperature obtained with the model described in [14] is 3.2 • C below the measured value. In addition, Figure 6a,d shows that the water-side temperatures at the PHE absorber outlet follow a linear trend in both models, however, the present model approach was closer to the measured value.
As regards Figure 6b,e, the outlet concentration estimated with the model described in [14] is closer to the measured values in comparison to that of the semi-empirical model. It is then observed that the semi-empirical model tends to overestimate the solution concentration at the PHE absorber outlet, however, errors do not surpass 1.3%. Furthermore, Figure 6c,f presents the thermal load distribution along the control volumes of the PHE. Results show that both models predicted similar thermal load trends along the PHE, however, the thermal loads at the PHE outlet estimated with the semi-empirical model are closer to those of the measured values.
The total absorption mass flux obtained from the semi-empirical model and that calculated from the experiments present relative differences of 1.4% and 17.6%, at heat dissipation temperatures of 35 and 40 • C, respectively. For the thermal load, these relative differences are around 8.7% and 4.6%, at heat dissipation temperatures of 35 and 40 • C, respectively. Regarding the procedure described in [14], the relative differences of the total absorption mass flux with respect to the experimental values are around 11.2% and 16.8%, at heat dissipation temperatures of 35 and 40 • C, respectively. For the thermal load, these relative differences are around 16.7% and 24.9%, at heat dissipation temperatures of 35 and 40 • C, respectively.
In general, the estimations of the absorption mass flux and thermal load with the semi-empirical model were sufficiently close to the measured values. Although the liquid-vapor bubble interface area was not considered, approaches achieved were satisfactory. Figure 6d shows an immediate drop in the hot-side temperature when employing the semi-empirical model because of the greater potential for heat transfer due to the larger temperature difference. In this case, the solution-side temperature at the PHE outlet using the present model is 2.3 °C lower than that of the measured value, whereas the outlet temperature obtained with the model described in [14] is 3.2 °C below the measured value. In addition, Figure 6a,d show that the water-side temperatures at the PHE absorber outlet follow a linear trend in both models, however, the present model approach was closer to the measured value.

Artificial Neural Network Results
In this subsection, the neural network model configuration and results are presented. As mentioned in Section 2.3, six input parameters at the absorber inlet were selected: absorber pressure (P Sol ), solution temperature (T Sol ), solution mass flow rate (ṁ Sol ), solution concentration (x Sol ), water-side temperature (T Cw ), and water-side flow rate (V Cw ). As output parameters, three performance parameters were selected to assess the NH 3 absorption phenomenon: the solution-side heat transfer coefficient (h Sol ), the absorption mass flux (F AB ), and the solution-side mass transfer coefficient (k m ). The ANN model was designed in three steps: learning, validation, and evaluation of the model. The learning step included training and design of the neural network. Moreover, the back-propagation technique was used for the ANN training. The training process finished when the relative difference between the ANN output data and the measured values (target) corresponded to the tolerance set. Different network designs were tested to obtain the best correlation between the output data and the targets, and the minimum number of neurons. Therefore, each ANN design was assessed using the linear regression coefficient (R 2 ) and the root mean square error (rmse). Figure 7 shows the final structure design of the ANN used to assess the H-type plate bubble absorber performance. Training results showed that a neural network employing a hidden layer with five neurons is a suitable structure for a satisfactory prediction of the selected parameters. As a result, a 6-5-3 ANN structure was defined as shown in Figure 7. Therefore, it consists of six neurons in the input layer, five neurons in the hidden layer, and three neurons in the output layer. When the number of processing neurons was reduced to fewer than five, the linear regression coefficient dropped. In addition, when the number of neurons was increased, no improvement in the linear regression coefficient was obtained. For instance, when the number of processing neurons was set to four, the linear regression coefficient dropped to 0.9096, 0.9429, and 0.8561, for the absorption mass flux, solution heat transfer coefficient, and mass transfer coefficient, respectively. By comparison, when the number of processing neurons was set to six, the linear regression coefficients were 0.9735, 0.9401, and 0.8968, for the absorption mass flux, solution heat transfer coefficient, and mass transfer coefficient, respectively. Equation (19) was used to predict the selected parameters for the evaluation of the NH3 absorption process in the H-type PHE. In this equation, i refers to the inputs (i = 6), j refers to the neurons defined in the hidden layer (j = 5), Pi refers to the input "i" normalized, IW is the weight matrix in the hidden layer, LW is the weight matrix in the output layer, b1 and b2 refer to the bias in hidden and output layers, respectively; a refers to the output parameters and o refers to the selected output. For instance, o = 1 stands for the output a1 = FAB, o = 2 to km, o = 3 to hsol. Table 3 shows the results from the ANN training process to model the selected parameters which allow assessment of the plate bubble absorber. It contains the weights and biases of the ANN finally resulting from the simulation study. Equations (20) to (25) show a sample of the calculations conducted for the output parameter 1 (o = 1). Parameter 1 refers to the absorption mass flux (FAB). Therefore, this output depends on the inputs, the weights, and the biases as shown below: Equation (19) was used to predict the selected parameters for the evaluation of the NH 3 absorption process in the H-type PHE. In this equation, i refers to the inputs (i = 6), j refers to the neurons defined in the hidden layer (j = 5), P i refers to the input "i" normalized, IW is the weight matrix in the hidden layer, LW is the weight matrix in the output layer, b 1 and b 2 refer to the bias in hidden and output layers, respectively; a refers to the output parameters and o refers to the selected output. For instance, o = 1 stands for the output a 1 = F AB , o = 2 to k m , o = 3 to h sol . Table 3 shows the results from the ANN training process to model the selected parameters which allow assessment of the plate bubble absorber. It contains the weights and biases of the ANN finally resulting from the simulation study. Equations (20) to (25) show a sample of the calculations conducted for the output parameter 1 (o = 1). Parameter 1 refers to the absorption mass flux (F AB ). Therefore, this output depends on the inputs, the weights, and the biases as shown below: 1+exp(n 4 ) + LW 1,5 1+exp(n 5 ) −(LW 1,1 + LW 1,2 + LW 1,3 + LW 1,4 + LW 1,5 ) + b2 1 (20) where n 1 , n 2 , n 3 , n 4 , and n 5 are determined from Equations (21) to (25).
n 5 = −2[IW 5,1 P 1 + IW 5,2 P 2 + IW 5,3 P 3 + IW 5,4 P 4 + IW 3,5 P 5 + IW 5,6 P 6 + b1 5 ] The experimental and estimated values using the present model were contrasted to establish the model prediction accuracy. Figure 8 shows the estimated or correlated data vs experimental data for each output parameter (F AB , h sol , and k m ). The black lines represent the trends of the results. The results evidenced a satisfactory correlation between the experimental and simulated data. This highlights that the ANN developed is suitable for estimating the selected parameters to assess the PHE absorber.
The ANN model can represent the collected data with average deviations of 4.26%, 4.2%, and 1.74%, for the absorption mass flux, solution heat transfer coefficient, and mass transfer coefficient, respectively. Moreover, the ANN can satisfactorily predict 97% of absorption mass flux data at relative differences below 10%. Similarly, the ANN model can predict 90% and 97% of the solution heat transfer coefficient and mass transfer coefficient data, respectively, at relative differences below 10%. Values of the corresponding fit parameter, R 2 , of the ANN model are indicated in Figure 8.
The experimental and estimated values using the present model were contrasted to establish the model prediction accuracy. Figure 8 shows the estimated or correlated data vs experimental data for each output parameter (FAB, hsol, and km). The black lines represent the trends of the results. The results evidenced a satisfactory correlation between the experimental and simulated data. This highlights that the ANN developed is suitable for estimating the selected parameters to assess the PHE absorber. The ANN model can represent the collected data with average deviations of 4.26%, 4.2%, and 1.74%, for the absorption mass flux, solution heat transfer coefficient, and mass transfer coefficient, respectively. Moreover, the ANN can satisfactorily predict 97% of absorption mass flux data at relative differences below 10%. Similarly, the ANN model can predict 90% and 97% of the solution heat transfer coefficient and mass transfer coefficient data, respectively, at relative differences below 10%. Values of the corresponding fit parameter, R 2 , of the ANN model are indicated in Figure 8.

Heat and Mass Transfer Rates Comparison
In this sub-section, the results of the absorption phenomenon in the H-type PHE using the semi-empirical model (Mod) and the ANN model (ANN) are contrasted with the available experimental results (Exp). This analysis was carried out under the same conditions; see Table 2. These conditions of operation correspond to those of a single effect absorption cooling system with NH 3 /LiNO 3 driven by low-temperature heat sources at a chilled-water temperature of 5 • C. Figure 9 shows that the semi-empirical and ANN models follow similar trends to those of experimental results. In addition, uncertainty errors for the experimental results are shown with bars with a resulting maximum value of ± 8% following the procedure recommended by the National Institute of Standards and Technology [50]. Experimental results in Figure 9a,c correspond to water-side temperatures of 35 and 40 • C, and solution-side mass flows ranging from 10 to 50 kg·h −1 . In these figures it can be observed that when the solution-side flow increases, the absorption mass flux rises from 0.0029 to 0.0061 kg·m −2 ·s −1 and from 0.0021 to 0.0044 kg·m −2 ·s −1 at the given temperatures, respectively. Moreover, the thermal load varied from 4.84 to 13.08 kg·m −2 and from 3.38 to 7.56 kg·m −2 for the same temperatures, respectively.
Experimental results in Figure 9b,d correspond to water-side temperatures of 35 and 40 • C, while the water-side flow ranges between 130 and 450 l.h −1 . In these figures can be noted that when the water-side flow increases, the absorption mass flux rises from 0.0056 to 0.0064 kg·m −2 ·s −1 , and from 0.0033 to 0.0044 kg·m −2 ·s −1 at the specified temperatures, respectively. Thermal load varies between 10.31 and 13.2 kW·m −2 , and between 6.02 and 7.97 kW·m −2 at the specified temperatures, respectively. experimental results. In addition, uncertainty errors for the experimental results are shown with bars with a resulting maximum value of ± 8% following the procedure recommended by the National Institute of Standards and Technology [50]. Figure 9 also presents the influence of the water-side temperature, solution-side mass flow, and water-side volumetric flow rate on the absorption mass flux and thermal load in the PHE. Furthermore, results in Figure 9a,c were collected at a water-side volumetric flow rate of 268 l.h −1 , and results in Figure 9b,d were obtained at a solution-side mass flow of 40 kg.h −1 . Experimental results in Figure 9b,d correspond to water-side temperatures of 35 and 40 °C, while the water-side flow ranges between 130 and 450 l.h −1 . In these figures can be noted that when the water-side flow increases, the absorption mass flux rises from 0.0056 to 0.0064 kg.m −2 .s −1 , and from 0.0033 to 0.0044 kg.m −2 .s −1 at the specified temperatures, respectively. Thermal load varies Results from the present study show that when the solution-side flow rises at the given values, the absorption process and heat transfer are positively affected in the corrugated PHE. This effect is obtained by the improvement in the solution and ammonia vapor mixing when increasing the solution flow regimen. Because the absorption process is exothermic, the higher the process, the more heat is released to the water-side and thus, the thermal load increases. It can also be observed that the vapor absorption improves as the water-side temperature at the PHE inlet decreases. This effect is associated with the improvement in the absorption and heat transfer potential; that is, the absorption capacity is enhanced when excess heat is released from the absorption process to the cooling water, allowing the solution to remain subcooled. As a result, this solution sub-cooling extends the difference between the current state and equilibrium state, and therefore, the solution can continue absorbing ammonia vapor. Regarding the effect of the water-side flow, only a slight increase in the heat and mass transport is observed. Therefore, the overall heat transfer coefficient is controlled by the solution-side heat transfer coefficient.
In general, the maximum and average relative differences between the semi-empirical and experimental models were 15.8% and 6.0%, respectively, for the absorption mass flux. In addition, 83.3% of the data showed a relative difference below 10.0%. In the case of the thermal load, the maximum and average relative differences were 12.5% and 6.4%, respectively. Results also showed that 86.7% of the data presented a relative difference below 10.0%. Figure 9 also shows that the relative differences between results from the experimental study and the semi-empirical models are more pronounced at a water-side temperature of 35 • C.
Predictions of the artificial neural network showed that the maximum and average relative differences between the ANN model and experimental results were around 10.8% and 4.3%, respectively, for the absorption mass flux. In addition, 90% of the data evidenced a relative difference below 9.0%. In the case of the thermal load, the maximum and average relative differences were 11.3% and 1.3%, Energies 2020, 13, 4313 16 of 20 respectively. Results also showed that 97% of the thermal load data presented a relative difference below 9.0%.
Results confirmed that the semi-empirical model is suitable for estimating temperatures and concentration at the inlet and outlet of each partial control of volume in the PHE absorber, in addition to heat and mass transfer. As in the case of most heat and mass transfer models, this model requires the thermodynamics properties of the fluids and dimensions of the bubble absorber under study as inputs. Alternatively, the ANN model can estimate the absorption process phenomenon in the PHE absorber with lower errors compared to the semi-empirical model. Although internal details of the absorption process phenomenon are not needed, the ANN model is not suitable for the prediction of internal parameters, such as temperature and pressure profiles, or local heat and mass transfer details. It is finally highlighted that the ANN is a more effective method for prediction and control of the absorption process performance parameters in the studied configuration compared to the semi-empirical method. Furthermore, the semi-empirical method is preferable depending on the levels of internal/local heat and mass transfer details requirements. Finally, it is worth noting that, because the semi-empirical model uses dimensionless correlations, it could be adapted to other working fluids and slight variations in the absorber geometry, whereas the ANN model is specific to the configuration investigated in the present paper and within the ranges provided for the training process (see Tables 1 and 2) and thermophysical properties of NH 3 /LiNO 3 .

Conclusions
This paper presented the development and comparison of two models, a one-dimensional semi-empirical model and an ANN model, designed and developed to quantify the heat and mass transfer in an H-type PHE in bubble mode employing a NH 3 /LiNO 3 mixture.
Correlations to estimate the solution-side heat transfer coefficient, mass transfer coefficient, and water-side heat transfer coefficient are proposed. Results from the numerical models, the semi-empirical and ANN models, were contrasted with experimental results and indicated the proposed approach was suitable. Moreover, the semi-empirical model provided a reasonable approach to internal temperature and concentration details; for instance, temperature and concentration at the PHE outlet were sufficiently close to the results obtained by experiment. Moreover, the ANN model allowed better estimation of the overall heat and mass transfer in the studied bubble absorber compared to that of the semi-empirical model. In addition, one of the advantages of the ANN is that there is no need for internal details of the absorption process and thermophysical properties, however, its application is limited to the configuration studied.
Results showed that the semi-empirical model was able to predict the absorption mass flux and thermal load with maximum differences of 15.8% and 12.5%, respectively, relative to the experimental results. In the case of the ANN model, the maximum relative difference was 10.8% and 11.3% for the absorption mass flux and thermal load, respectively. Moreover, the semi-empirical model is preferable when more details of the internal mass absorption of the absorber are needed.
Finally, results showed that the correlations developed and proposed provided a good prediction of local heat and mass transfer along the absorber. Analysis of the simulation results also demonstrated that both models can be practically applied to the design of bubble plate absorbers with NH 3 /LiNO 3 mixtures.  Greek Letters λ thermal conductivity, kW·m −1 ·K −1 α heat transfer coefficient, kW·m −2 ·K −1