Experimental Evaluation and Theoretical Optimization of an Indirect Solar Dryer with Forced Ventilation under Tropical Climate by an Inverse Artiﬁcial Neural Network

: In this theoretical–experimental study is presented a hybridization strategy based on the application of an inverse artiﬁcial neural network model (ANNi) coupled with metaheuristic optimization algorithms to optimize the drying velocity ( v d ) of an active indirect solar dryer for plantain and taro ( Colocasia antiquorum ). In the experimental stage, both fruits were evaluated in periods from 9:00 a.m. to 5:00 p.m. under a humid tropical climate region, varying the voltage of the air extractor fan (at 6 V, 9 V, and 12 V) to control the fan velocity. The experimental results showed that the maximum drying velocities were reached at 9 V, achieving a drying velocity of 1.5, 0.9, and 0.55 g/min, with a total drying time of 465 min for the taro, and 1.46, 1.46, and 0.36 g/min, with a total drying time of 495 min, for the plantain. As a result of the drying curves, it was observed that the drying velocity is higher in taro than in plantain. Subsequently, an artiﬁcial neural network (ANN) architecture was trained using the database generated from the solar dryer’s experimental stage. Six environmental variables and one operational variable were considered as the model’s inputs, feeding the ANN to estimate the drying velocity ( v d ), obtaining a linear regression coefﬁcient R = 0.9822 between the experimental and simulated data. A sensitivity analysis was performed to determine the impact of all the input variables. A hybrid strategy based on ANNi was developed and evaluated with three metaheuristic optimization algorithms; the best result was obtained by genetic algorithms (ANNi-GA) with an error percentage of 0.83% and an average computational time of 4.3 s. The scope of this optimization strategy was to obtain a theoretical result that allows predicting the behavior of the dryer and improving its performance for its practical application in future work through the implementation in development boards. Lastly, the optimization strategy presented is not limited to indirect solar dryers and opens a research window for evaluating other solar drying technologies.


Introduction
In Latin America, the Mexican agricultural industry is one of the most significant, establishing itself as the largest producer of vegetables in the region and occupying the second place in fruit exports worldwide. Agricultural activities play an essential role in the Mexican economy, with a contribution to gross domestic product (GDP) of more than 2% and an economic spillover from foreign direct investment (FDI), which in 2019 exceeded 50 million USD [1]. Among the country's growing regions, the south-southeast zone (made up of the States of Tabasco, Chiapas, Campeche, and southern Veracruz) stands out due to its favorable climatic and hydrographic conditions, where 96% of the agricultural surface has rainwater irrigation. In this region, 50% of the entire country's water runs off, and the prevailing tropical rainy climate allows for harvests practically all year round [2]. From 2011 to 2017, the south-southeast region sustained the highest agricultural growth, with an annual rise in the area, in real terms, of 35.7%, making this region very important for the national agricultural use [3]. However, the fact of having harvests almost all year round brings about an economic problem related to the placement and sale of the products. This process is subjected to the postharvest lifetimes of the crop since it can generate an excess of production, saturating the market, or generate economic losses for farmers, due to wastage. Therefore, alternatives are required to preserve seasonal food, helping the development of the rural economy.
One of the most attractive options for preserving fruits, vegetables, and herbs is drying and dehydration, which also provides other benefits such as lowering transport costs, increasing lifetime, and reducing crop waste [4]. Specifically, in regions with a tropical climate, solar dryers represent a feasible option for the development of such activity from an ecological approach, presenting lower capital and operating costs than conventional dryers [5]. In solar dryers, thermal energy from near-infrared wavelengths is harnessed to supply heat for the drying process and evaporate the required moisture content from the product. The effectiveness in the use of this technology depends on the radiation levels in the region. In the case of Mexico, the annual average daily irradiation on a horizontal surface is just over 4.4 kWh/m 2 ·day. Specifically, in the south-southeast region, there is average daily irradiation of 5.5 kWh/m 2 ·day due to its geographical location between the tropics of Cancer and Capricorn, the area where the most significant amount of solar radiation is captured in the world [6]. Therefore, the use of solar thermal technology in Mexico contributes significantly to the energy requirements for postharvest processing, with the adequate technical feasibility to be used as a technological improvement in agricultural areas of tropical climate regions (like the south-southeast region of Mexico).
Due to the agricultural biodiversity of this tropical region, in recent years, the practical potential of using solar thermal technology drying has been analyzed for various products with good commercial value for the south-southeast of Mexico [7]. M. Castillo Téllez et al. [8] developed an indirect tunnel type solar dryer for application in red chili. They evaluated the drying kinetics obtaining drying times of 2.75, 3.0, and 6.25 h, for temperatures of 55 • C, 65 • C, and 45 • C, respectively. They also evaluated the velocity ranges, between 0.7 and 2.6 m/s, reaching drying times and final moisture contents of 16 h, with 0.057 to 0.90 kg water /kg dry mass, and 21 h, with 0.0611 to 0.109 kg water /kg dry mass, respectively. López-Vidaña et al. [9] developed a passive solar dryer to operate as a mixed or indirect dryer for tomatoes in a tropical region; the drying efficiency was in a range between 5.47% and 4.48%, respectively. For their part, M. Castillo Téllez et al. [10] presented an experimental study for a Stevia leaf solar dryer. A comparative study was carried out between direct and indirect drying by natural and forced convection. The heat supplied by a solar water heating system was used through a water/air heat exchanger. Natural convection produced a drying time of 360 min, while forced convection produced a drying time between 550 min and 600 min. Without a shade mesh, these values were 250 min and 300 min, respectively. Likewise, it was shown that an indirect solar dryer can be controlled better than a direct one.
In the southeastern region of Mexico, there are two crops that have excellent economic viability due to their characteristics. One is the plantain, a product highly commercialized nationally and internationally, which can expand its economic potential if its postharvest processing is diversified and its production is maximized. Another crop is taro, a product with good nutritional properties, which has been underutilized; hence, its processing and use make it a market niche in the agro-industrial and food sector. To dry these products and to preserve their organoleptic properties, it is necessary to achieve drying temperatures that range between 60 • C and 75 • C. This makes it necessary to use low-temperature solar thermal technology capable of establishing limits to the UV radiation that may damage the samples. Thus, the best option for this purpose is an indirect solar dryer.
However, an indirect solar drying system comprises processes where different interconnected physical phenomena take place, such as energy interactions, mass exchanges, fluid movements, and heat flows, to mention a few. Therefore, it is a complex system whose analytical solution cannot be found, which implies finding numerical solutions [11]. On the other hand, during the drying process, it is crucial to monitor the different variables such as environmental conditions, temperatures inside the solar dryer, relative humidity, air velocity, incident solar radiation, the type of product to be dried, the initial moisture contents, and total mass of the product, since they are variables that affect the drying of the product [12]. Being able to determine the best operating conditions for the drying system implies cost reduction, higher thermal performance, reduced drying time, better use of solar energy, and better dehydration. Therefore, it is necessary to develop and implement tools that help solve the optimization problem [13].
Therefore, it is important to establish an optimization criterion that, as a function of the design variables, as well as the modeling, simulation, and validation stages, allows us to minimize an objective function to satisfy the system's operating conditions and the physical restrictions of the system [14]. Furthermore, many control and optimization problems can be solved by simple regulators; however, some control loops require advanced optimization techniques, which allow us to achieve the desired optimal operating conditions and, at the same time, reduce the complexity of their solution [15,16]. A feasible solution is the use of artificial intelligence modeling techniques, combined with reverse-engineering strategies and heuristic optimization algorithms, such as artificial neural network (ANN), genetic algorithm (GA), particle swarm optimization (PSO), ant colony optimization (ACO), and deep learning [17,18]. Table 1 shows a qualitative analysis of the most recent studies for indirect solar dryers. As can be seen, these studies are categorized into different aspects, such as type of research, the modeling technique applied, the purpose of the research, variables of interest, type of ventilation, and type of analysis. Concerning the type of research, it was detected that few works simultaneously addressed the modeling and experimental approaches of indirect solar dryers [19,20]. In addition, it was identified that there is only one work that used artificial intelligence tools (diffuse modeling, Takagi-Sugeno type) for modeling the thermal phenomena that occur during the drying of the product [21]. With regard to the variables of interest, it was identified that the drying time was the most analyzed, followed by the drying efficiency, the pick-up efficiency, and the exergetic efficiency [22]; however, there are no reports of the drying velocity, despite it being an important parameter for the development of control strategies focused on improving production. Lastly, there are no works that analyzed the indirect solar dryers from the perspective of sensitivity analysis.
According to the literature review, there is no report of a study focused on optimizing and controlling a solar drying device; thus, a hybridization strategy with artificial intelligence and metaheuristic optimization techniques has not been used either. In the same sense, none of the reported works considered the drying velocity as a variable of interest to control. Furthermore, despite the knowledge of the variables that influence the operation of solar dryers, there are no studies that quantitatively describe the level of impact that environmental and operational variables have on the final product.
For this reason, in this work, an experimental study of an indirect solar dryer for tropical fruits of high economic potential was developed, using samples of plantain and taro, and a hybrid strategy was subsequently implemented based on an artificial intelligence model and a metaheuristic optimization algorithm. The optimization process considered the voltage of the air extractor inside the dryer as an operating parameter in order to improve the drying velocity. The objective was to theoretically establish the optimal control values for the drying operation from an inverse neural network solved with three metaheuristic optimization algorithms, in addition to determining the theoretical optimal drying speed for the process. An experimental study of a solar dryer for plantain and taro is presented to establish the optimal operating conditions during drying in this work An artificial neural network inverse was used, and three metaheuristic optimization strategies were evaluated: genetic algorithm, particle swarm optimization, and ant colony optimization, to maximize the drying velocity of the process

Description of the Experimental System
The experimental system used in the present study corresponds to a prototype of a solar dryer that operates by forced ventilation, designed mainly for drying plantain and taro. The prototype was evaluated in the Municipality of Cunduacán Tabasco, México, at the facilities of the Chontalpa Campus, Academic Division of Engineering and Architecture of the Juárez Autonomous University of Tabasco (18 • 4 24.27 N latitude and 93 • 9 52.72 W longitude). Figure 1 shows the schematic diagram of the dryer used. It consists of a suspended absorber plate of stainless steel painted with a selective coating of nanostructured black nickel (α = 0.91; = 0.19), with a thickness of 4 mm and an area of 0.7 m 2 . The dimensions of the solar collector are 1.0 m × 0.7 m, with a fixed inclination angle of 18 • . In order to minimize the heat losses to the environment, the bottom and sides of the collector were built with a central piece of thermal insulation using polyurethane panels coated with galvanized steel, with a thickness of 3.81 cm and a coefficient of thermal conductivity (k) of 0.018 W·m −1 ·K −1 . The transparent cover is 3 mm thick clear glass (t c = 0.84, c = 0.9 and α c = 0.06), to reduce heat loss at the entrance of solar radiation. The dimensions of the drying chamber are 0.7 × 0.5 × 0.8 m 3 , whereas its walls are also made with polyurethane panel covered with galvanized steel. The inside of the chamber has a metal guide to place the drying tray with dimensions of 0.6 × 0.45 m 2 , constructed with galvanized steel mesh. Above the drying chamber, there is a 1 m high chimney, with an opening in the outlet to allow the circulation of humid air outside the dryer. Lastly, to generate forced ventilation, a direct current (DC) fan with nominal values of 60 W, 12 V, 5.0 A, and 1800 rpm was installed on top of the chimney.
The solar dryer was instrumented to measure the environmental parameters that affect its performance and the operational parameters of interest. For the measurement of solar irradiance, an Amprobe Solar 100 pyranometer mounted on the solar dryer at an inclination of 18 • was used. The temperature and relative humidity of the environment were sensed using an Extech RH101 Hygro-Thermometer hygrometer, while an Extech 451,126 fin anemometer was used for wind velocity. Concerning the operational parameters, the temperature inside the dryer was measured using type K thermocouples, distributed across the solar collector, the drying chamber, and the chimney. On the other hand, the air velocity inside the collector was measured using an Extech 451,126 anemometer, while the mass of the product to be dried was measured using a Torrey L-PCR 20 balance. Table 2 presents the technical information with the precision and operation ranges of these sensors. The sensors were connected to an ADAMS-4016 8-channel data acquisition card for the automated recording of information storage. The data were sent by network cable to a data logger, interconnected to a PC. The sampling period was 1 min with averages of 15 min. The selection of the sampling time at 15 min was carried out according to the design criteria for solar trackers, considering that the solar angle moves 1 • every 15 min [12]. The schematic diagram of the experimental system is shown in Figure 2.   To measure the mass of evaporated water, a kilogram of fruit was sliced, placed in trays, and subsequently introduced to the solar dryer. The samples were left in the dryer for a period of 8 h. During this time, the set of all samples was weighed every 15 min and returned to the tray, identifying the evaporated water from the samples through weight loss. For practical terms, in each measurement, the mass of all the slices was quantified each measurement interval. In this period, the mass of water evaporated from the product was recorded to calculate the drying rate. This procedure was performed separately for both plantain and taro.
The experimental drying velocity was estimated as the ratio between the mass of evaporated water of the product and the drying time interval of the same. These values were obtained by measuring the loss mass of the fruit and the time elapsed between each sampling. This relationship is expressed by the following equation [34]: where m w is the mass of evaporated water of the product (g), ∆t is the drying time interval (min), and v d the drying velocity (g/min). Figure 3 presents a diagram of the computational methodology used in the present work to specify the optimization process of the drying velocity for the proposed solar dryer. As can be seen, this is made up of three stages. The first stage consisted of creating an experimental database from the drying of plantain and taro by the photothermal device. This database was used in the second stage for the training and creation of a multivariate regression model based on ANN that correlated the environmental and operational measurements with the drying velocity of the samples. Once the optimal ANN model was obtained, it was used as an objective function to evaluate various metaheuristic optimization algorithms and determine which provided the best numerical results.

Experimental Database Construction
After implementing the experimental system, the measured data were collected and arranged in a spreadsheet according to the variables of interest. In total, a database made up of 3168 data was formed from eight experimental variables. For the development of the correlational model, the following variables were considered: solar radiation, ambient temperature, relative humidity, wind velocity, input temperature dryer, output temperature dryer, and fan voltage as the independent variables, and the velocity drying as the output variable. Table 3 summarizes the ranges of the variables considered in the present study. Before training the ANN, the input variables were normalized because of the active functions' nonlinearity and the operating range fluctuating between 0 and 1. These are special sigmoidal functions for activating the hidden layer. This normalization allows the surface of the error function to be reduced, making training more efficient, considering the initialization of the variables to zero. Each input variable was normalized using the following equation [35]: where u i,N is the normalized variable, u i is the real value of the variable before normalization, and u max is the maximum value of the variable before normalization.

Artificial Neural Network Modeling
ANN was used as a computational tool [36] for the development of a mathematical model that correlates the environmental and operational parameters of the solar dryer to the drying velocity. The main reason for its choice is that it is one of the multivariate regression techniques with the greatest adaptability in complex thermal processes such as those exposed in this work [37]. In short, an ANN is a set of mathematical subfunctions called neurons, which are distributed in layers and interconnected with each other through constant values (weights). The ANN architecture used in this work consisted of an input layer, a single hidden layer, and an output layer. The numbers of neurons in the input and output layers were defined by the experimental variables in Table 3, while the hidden layer was made up of sigmoid neurons whose number was difficult to specify and required an iterative process for identification. In this work, two ANN architectures were evaluated with two different activation functions in the hidden layer, a hyperbolic tangent function and a logistic sigmoid, whose mathematical expressions with respect to their network architectures are presented in Equations (3) and (4).
where W i and W o correspond to the connection weights of the neurons between the inputhidden layer and the hidden-output layer, respectively. The elements b 1 and b 2 represent factors of adjustments of the neural network associated with the neurons of the hidden layer and output layer, respectively. S and K represent the total number of neurons in the hidden layer and the input layer, respectively. Lastly, In represents the k-th input variable of the mathematical model. For the development of the ANN model, the normalized input parameters contained in the database were used ( Table 3). The data were randomly divided into subsets of training (80%), validation (10%), and testing (10%). By using a backpropagation algorithm to train the ANN through an iterative trial and error process, different numbers of neurons in the hidden layer were evaluated to find the optimal values of W i , W o , b 1 and b 2 for Equations (3) and (4) that minimized the error between the experimental data (x out,exp ) and those estimated by the ANN (x out,sim ). In this work, the backpropagation algorithm applied was Levenberg-Marquardt. This is a modification of the gradient descendent algorithm that performs numerical minimization processing to solve nonlinear least squares problems from an estimated initial parameter vector. The main advantage of using this algorithm is its fast convergence time, which is why it is widely used in ANN training [38].
To quantify the performance of the neuronal model, the following statistical criteria were used: root-mean-square error (RMSE), linear regression coefficient (R), coefficient of determination (R 2 ), and mean absolute percentage error (MAPE), synthesized in Table 4. Determination coefficient

Optimization Approach by Inverse ANN
In order to identify the fan voltages that guarantee optimal forced convection conditions to improve the operation of the solar system and achieve the maximum drying velocity, a computational strategy known as ANNi [39,40] was implemented, which is based on the algebraic reduction of Equations (3) and (4). This approach uses the coefficients of the weights and bias of the ANN obtained after their training (W o , W i , b 1 , and b 2 ) and inverts the ANN to form the objective function defined by Equation (5).
where y is the function of the experimental variable to be optimized (y = f (In(v d )).
Equations (6) and (7) present the mathematical expressions for the case of ANN architectures with hyperbolic tangent and logistic sigmoid activation functions in the hidden layer neurons, respectively.
where X(s) is determined using Equation (8), from which K acquires all In input values except for the one corresponding to v f (In v f = In (7)).
For the solution of both objective functions, the evaluation of three metaheuristic optimization algorithms is proposed to determine the most efficient alternative dryer operation with the least computation time. The metaheuristic algorithms evaluated were genetic algorithm (GA), particle swarm optimization (PSO), and ant colony optimization (ACO). A detailed description of their operating principle is provided in Appendix A.

Experimental Evaluation of the Solar Dryer
The forced convection solar dryer prototype was experimentally evaluated during March, April, and May 2019, considering a drying schedule from 9:00 a.m. to 5:00 p.m. The experiments were carried out by applying different operating voltages to the fan (6, 9, and 12 V). Figures 4 and 5 show the curves of moisture content on a dry basis and the velocity of drying referred to solar radiation for a test day, for taro and plantain, respectively.

Taro (Colocasia esculenta)
Due to cloudiness, solar radiation had significant variations, which caused fluctuations in the air temperature at the entrance to the drying chamber. Three important aspects are combined in the phenomenology corresponding to the drying process: forced ventilation, the impossibility of maintaining constant the temperature of the air surrounding the samples, and the porosity of the taro. Due to its porosity, the critical moisture content in taro is high, which implies that the samples have a greater amount of water in the last drying stage (decreasing drying stage) and higher drying velocity than the plantain. The drying curves reflect the problem of migrating from the free moisture content zone and reaching the equilibrium moisture content, causing disturbances in the drying velocity curves. This is because the heat transfer rate is different from the mass transfer rate, which does not allow constant drying rates, and sudden changes occur since the concentration gradient decreases, sharply decelerating the mass transport and delaying the filling of water in the interstitial spaces of the taro by capillarity. However, the temperature gradient between the inlet of the solar collector and the outlet of the drying chamber could be maintained, and, despite the disturbances during the test, the evaporative drying conditions and behaviors were established. In Figure 4, the drying curves for taro are observed, which were sectioned according to the deviations that occur in the drying velocity. The choice of the time ranges indicated in the subfigures was given by observing changes in the drying velocity curves, which showed different stages within the dryer. On the other hand, the derivative of the moisture content on a dry basis was considered because the drying speed was evaluated through the variations of the mass per unit of time, for each sampling period. In the case of the fan voltage at 6 V Figure 4(a1,b1), the drying time ranges were 0 to 200 min, 200 to 300 min, 300 to 450 min, and 450 to 525 min, with average drying velocities for each region of 1.8, 1.5, 0.9, and 0.3 g/min, respectively. These correspond to the changes in the dry basis moisture curve in the same ranges. According to the figures, despite the environmental conditions, the drying process stages were present, and a more stable drying velocity in the decreasing drying region was achieved, with a decrease of 33% compared to the previous region. In this case, the solar radiation and the ambient temperature conserved almost the same morphology during the entire test time, except for the period between 3:00 and 5:00 p.m. The temperature rose from 42 • C to 44 • C. This temperature rise helped sustain the drying process at the end of the curve. In the case of the 9 V fan voltage (Figure 4(a2,b2)), the drying time ranges were 0 to 300 min, 300 to 400 min, and 405 to 465 min, with average drying velocities for each region of 1.5, 0.9, and 0.55 g/min respectively, corresponding to the changes in the dry basis moisture. The drying velocity curve in the first range presented a constant and prolonged rate of variation, while the drying velocity in the last time range decreased by 61% and stably maintained its decreasing velocity. According to the variations in solar radiation and ambient temperature, their behavior was at times similar in morphological terms, even more so from 3:30 to 5:00 p.m.; solar radiation decreased gradually, while the temperature remained oscillating between 40.7 • C and 41.5 • C, helping the drying process in its last stage.
In the case of the 12 V fan voltage (Figure 4(a3,b3)), the observed drying time ranges were from 0 to 270 min, from 270 to 405 min, and from 405 to 465 min, with average drying velocities for each region of 1.6, 1.17, and 0.53 g/min, respectively. These results correspond to the changes in the rate of variation of the dry basis moisture. Except for a disturbance at the beginning of drying, solar radiation remained with a favorable rate of change, maintaining its morphology concerning ambient temperature, including during the last hour and a half of the test where the decrease in solar radiation and ambient temperature agreed with the drastic decrease in drying speed. The average drying velocity decreased by 45.3% in the last time range, corresponding to the stability zone of the drying velocity.
From the conditions generated by the drying chamber and forced ventilation, and from the behavior of the drying curves, it was identified that the fastest drying velocities occurred when operating the fan at 9 V. It was observed that the average drying rates for the equilibrium moisture content drying region for the three voltage conditions were similar, with a value of 1.5 g/min, as well as for the first drying region, decreasing by 0.9 g/min. In contrast, in the second part of the decrease in drying, there was disparity, since the decrease percentages varied for each voltage at 33% for 6 V, 61% for 9 V, and 45.3% for 12 V, revealing the best drying conditions the 9 V experiment. From the above, we could perform a search to find the optimal voltage for drying, and this also allowed us to reduce the optimal response surface of the search.

Plantain
Similar to the case of taro, solar radiation had variations, albeit less significant in this case due to the presence of low cloudiness in the experiment period, improving the air temperature conditions at the entrance to the drying chamber. However, it was still not possible to keep a constant temperature of the air surrounding the fruit.
Contrary to taro, plantain porosity is low, also reducing the critical moisture content in the plantain. Due to this, plantain samples have less water in the last drying stage, and their drying velocity is lower than that of taro. Similar to the case of taro, there were disturbances, but they were decreased due to better conditions of solar radiation and a decrease in the porosity of the samples. The phenomenology was similar to that which occurs in taro. The temperature gradient for hot air circulation was still achieved, establishing the conditions and behaviors of drying by evaporation and improving the drying curves, as shown in Figure 5.
In the case of 6 V fan voltage, (Figure 5(a1,b1)), three drying time ranges were identified from 0 to 225 min, 225 to 350 min, and 350 to 465 min, with average drying velocities for each region of 1.27, 1.25, and 0.255 g/min respectively. This corresponded to the rate of variation of the dry basis moisture, which, at the beginning of drying, was observed to decrease, but not remain constant. The same behavior occurred for solar radiation, which only suffered from some disturbances at the beginning of the drying time. In Figure 5(a1), after 150 min, a behavior derived from a combination of effects was observed. On the one hand, there was a measurement error, which could be subject a change in the relative humidity of the environment. On the other hand, in the solar radiation curve, a sudden decrease of 0.735 kW/m 2 occurred between 2:00 and 3:00 p.m., due to the presence of cloudiness. The temperature developed in a range between 32.3 • C and 44 • C, with a morphology similar to solar radiation, except for the sudden decrease in radiation, where the slow and delayed nature of the heating process was revealed, as there was no change in temperature until after the sampling period. In the last drying stage, a decrease of 20.4% in the drying velocity was achieved, 12.6 percentage points below the taro for the same supply voltage.
For the fan operating at 9 V ( Figure 5(a2,b2)), two ranges of drying time were observed from 0 to 270 min and from 270 to 495 min, with average drying velocities of 1.46 and 0.36 g/min, respectively. The rate of variation of the dry base humidity ratio during the start of drying decreased and remained almost constant, thus covering a wide range of time (4.5 h). This was due to a reduction in disturbances of solar radiation, as well as the ambient temperature, reaching values between 38 • C and 28 • C, which also allowed a stabler and broader decreasing drying region of almost 3.75 h, with an average drying velocity of 0.36 g/min and a 25% reduction in drying velocity, 36 percentage points below the taro for the same supply voltage.
In the case of the fan at 12 V, (Figure 5(a3,b3)), three ranges of drying time were presented from 0 to 200 min, 200 to 345 min, and 345 to 435 min, with average drying velocities for each region of 1.95, 0.79, and 0.12 g/min, respectively. It was observed that disturbances only occurred in the solar radiation curve in the second region of the drying time, with a negative peak of 0.232 kW/m 2 between 11:00 and 11:30 a.m.; with these exceptions, solar radiation and ambient temperature maintained their relative morphology, with a range between 27.9 • C and 36 • C, helping the drying process. On the other hand, in the last stage of a decrease in drying, the average drying velocity was 0.12 g/min, with a reduction in drying velocity of 15.4%, 29.9 percentage points below the taro for the same feed voltage.
The fan voltage with the highest detected drying velocities was 9 V, reaching maximum average drying velocities of 1.46 and 0.36 g/min for each drying region, with a total drying time of 495 min. It was observed that the average drying velocities for the initial drying region and the three voltage conditions were consistent, in a range from 1.27 to 1.952 g/min. In contrast, there were significant deviations for the final drying region (0.255, 0.36, and 0.12 g/min for 6 V, 9 V, and 12 V, respectively) with percentage decreases for each voltage of 20.4%, 25%, and 15.4%, respectively. Therefore, it was possible to perform a faster and narrower search for the optimal voltage for plantain drying.

Artificial Neural Network Model
A mathematical model based on ANN was developed to obtain a numerical expression that relates the environmental conditions of the experiment and the operating variables of the solar dryer. All the numerical calculations were executed through the MATLAB software and its ANN toolbox complement, due to the ease and flexibility of the platform for the application of artificial intelligence codes [38]. This model considers G T , T a , RH, w v , T i , T o , and V f as the input variables and the drying velocity (v d ) as the output variable.
The ANN model was trained through a supervised learning process in which 80% of the data were used to train the model. This process was aimed at minimizing the error between the measured dry velocity and that generated by the model. In parallel, 10% of the data assigned to validation were used to avoid overtraining the model and guarantee its generalizability. Lastly, the remaining 10% of data were used to verify performance using unknown data. During this supervised learning process, the maximum number of iterations and the validation check (deviation index between the training and validation results) were considered as stopping criteria. Furthermore, in order to enable the general predictive capacity of the ANN, the k-validation approach suggested by Paneiro and Rafael [41] was conducted. Table 5 presents the evaluation of the different network architectures (using the hyperbolic tangent and logistic sigmoid activation functions) to find the optimal ANN model. According to the analysis, for both cases, although the network architecture with 10 hidden neurons presented the best performance, the differences obtained with respect to the architecture with eight hidden neurons were minimal. Thus, to obtain a shorter computing time and better use regarding computational resources, the latter was chosen; the architecture of the artificial neural network can be observed in Figure 6.   Figure 7 shows the regression curves between the simulated (v d,sim ) and experimental (v d,exp ) data of v d , for the models with the hyperbolic tangent ( Figure 7a) and logistic sigmoid (Figure 7b) transfer functions. As can be seen, the model based on hyperbolic tangent (v d,sim = 0.9865 v d,exp + 0.1000) was the one with the least dispersion of data and was more in line with the values measured experimentally. The above is supported by the statistical information in Table 5, where, for the model with eight neurons in the hidden layer and with the hyperbolic tangent transfer function, the statistical parameters were as follows: RMSE = 0.1912, R 2 = 0.9647, R = 0.9822, and MAPE = 31.9%. For the same neural model, but with the logistic sigmoid transfer function in the hidden layer, we had the following parameters: RMSE = 0.2139, R 2 = 0.9521, R = 0.9757, and MAPE = 25.4%. In relation to results with this same approach reported in the literature, Cervantes-Bobadilla, et al. [42] reported acceptable RMSE values of 0.4052 and 2.9178, while, in this work, the corresponding values were 0.1992 and 0.2139. Furthermore, Ajbar, et al. [40] performed modeling with a similar ANNi approach, obtaining R values of 0.9753 and 0.9678, which were considered acceptable to proceed with the optimization phase. Compared to this, this study presented better correlation values (R of 0.9822 and 0.9757). In the case of MAPE, high values were presented, although the model was adequately adjusted, as proven above and as shown by the regression curves (See Figure 7). This is due to the fact that the drying curves presented regions close to zero, due to the phenomenology of the drying process itself, where the drying rates were close to zero toward the end of each experiment. This fact is important because the MAPE divides the absolute error from the real data; hence, real values that are close to 0 can significantly increase the MAPE [43]. Thus, the mathematical function that describes the behavior of the solar drying process under study was given by Equation (3). Table 6 summarizes the weights and biases of the neuronal model for the hyperbolic tangent transfer function. The neural model used these parameters to reproduce the values of v d during the optimization process. Table 6. Weighting and bias parameters obtained for the best ANN model.  a s is the number of neurons in the hidden layer, k is the number of neurons in the input layer, and l is the number of neurons in the output layer (l = 1).

Garson Sensitivity Analysis
Once we obtained the artificial neural network model synthesizing the complex thermal processes inside the solar dryer in a single mathematical expression, a sensitivity analysis was performed using the Garson equation [44]. This determined the relative importance of the input variables on the drying velocity of the agricultural products evaluated. For this purpose, Garson's equation takes into consideration the values of the weights (presented in Table 6), estimating the relative importance of each variable from the following mathematical expression: where I j is the relative importance of the input variable j over the output variable, N i and N h are the numbers of input and hidden neurons, respectively, and W is the weight of the connection. The superscripts "i", "h", and "o" refer to the input, hidden, and output layers, respectively, and the subscripts "k", "m", and "n" refer to the input, hidden, and output neurons, respectively. Figure 8 presents the results of the sensitivity analysis. It was determined that, except for the wind velocity, all the variables significantly impacted the drying velocity. In quantitative terms, the percentage of relevance was as follows: solar radiation (G T ) (17.6%), ambient temperature (T a ) (16.2%), relative humidity (RH) (15.1%), wind velocity (w v ) (0.3%), drying chamber inlet temperature (T i ) (17.4%), drying chamber outlet temperature (T o ) (17.7%), and fan voltage (V f ) (15.7%). From these results, it can be seen that, due to it being a forced convection system, the effect of natural wind circulation was diminished. On the other hand, the wind velocity had a very low influence because there was not much variation in the database during the experiment (See Table 3), and the drying system was closed. The major influence regarding the wind velocity was represented by losses via convection on the surface, thus presenting a lower influence in the drying process. Likewise, the other parameters presented a similar significance level, differing by less than 2%. The fan voltage was shown to have an important impact, validating the application of control strategies as it is the only manipulable element.
Derived from the decoupling of wind velocity by forced ventilation, the other input variables were uniformly distributed, with an average contribution value of 16.61%. This means that they were strongly linked, due to the energetic interactions of the heat transport phenomena that occur in the drying process, such as the supply of heat by solar radiation to the solar collector. This heat is further injected into the air stream, raising its temperature to generate a gradient and provide the hot air flow, needed to remove any excess of water.

Optimal Operation Conditions
In this section, the inverse artificial neural network (ANNi) approach was implemented, considering as an objective function the ANN model developed (architecture 7-8-1; hyperbolic tangent activation function in the hidden layer). Based on the sensitivity analysis results and the discussion in Section 3.3, the ventilator's operating voltage (V f ) was considered as the control variable. To establish a strategy in the theoretical operation that allows maximizing the drying velocity of both fruits (plantain and taro), the optimization of the objective function was expressed as follows in order to evaluate various operating scenarios considering the range used in the experimental process: maximize : min(y) = −v d,sim V f ,subject to : V f j,min < V f j < V f j,max where j = 6to12. (12) To execute the optimization by ANNi, a Matlab script, aided by the Global Optimization Toolbox [45], was created in which the feasibility of the solution was evaluated using three different metaheuristic algorithms: genetic algorithm (GA), particle swarm optimization (PSO), and ant colony optimization (ACO) (described in detail in Appendix A). These algorithms were selected because they have shown successful results in the optimization of complex problems in diverse fields of research [46][47][48][49].
Using Equation (13), the error percentage (ε) was calculated to identify the ANNmetaheuristic optimization tuple that best approximates its calculations. In the same way, the computation time was recorded due to its importance during the possible execution of the optimization strategy in the operation of the fan voltage in real time. Table 7 shows the results of implementing ANNi, considering a set of random samples of both plantain and taro that were not part of the training or validation model stages. The execution of the algorithms was carried out on a personal computer with an Intel i7 processor, eight cores, 16 threads, and 32 GB of RAM. Because the algorithms work based on stochastic methods, a total of 20 runs were performed to identify the average computation time to reach the solution. It is observed that, for this function, the smallest ε for each optimization strategy was as follows: 0.83% for ANNi-GA, with a computation time of 4.3 s; 1.11% for ANNi-PSO, with a computation time of 8.7 s; 3.00% for ANNi-ACO, with a computation time of 9.2 s. As can be seen, ANNi-PSO and ANNi-ACO presented the slowest calculation times of the three techniques evaluated; however, the percentage error for the optimal fan voltage for both algorithms was the highest. In comparison with ANNi-PSO and ANNi-ACO, ANNi-GA generated lower error rates, lower by 0.28 and 2.17 percentage points, respectively. Moreover, it presented better calculation times, with differences of 4.4 s and 4.9 s, respectively. The results of the table indicate that ANNi is an appropriate tool for application in indirect solar dryers, obtaining, for taro, error rates in the fan voltage lower than 3.33%, 5.56%, and 6.67% for ANNi-GA, ANNi-PSO, and ANNi-ACO, respectively, according to the experimental data. In the case of plantain, the errors were 5.33%, 4.58%, and 4.89% for ANNi-GA, ANNi-PSO, and ANNi-ACO respectively. By applying three different optimization algorithms to the ANNi, it was possible to obtain the input values associated with the solar dryer. As for the computational time, ANNi-GA applied to the plantain (5.3 s) required less time than for taro (5.5 s). When comparing the computation time with the different algorithms for the case of taro, a maximum difference of 4.4 s was obtained between ANNi-GA and ANNi-PSO, while 2.9 s was measured between ANNi-GA and ANNi-ACO, and 1.5 s was obtained between ANNi-PSO and ANNi-ACO. For plantain, the maximum computation time difference was between ANNi-GA and ANNi-ACO with 14.6 s, whereas the difference was 10.8 s between ANNi-PSO and ANNi-ACO and 3.8 s between ANNi-GA and ANNi-PSO. From the above, it was concluded that ANNi-GA is more likely to generate reduced computation time values. Table 7. Comparison between experimental data and those obtained from the ANNi model applying the hyperbolic tangent transfer function for taro and plantain.

Experimental Data
ANNi-GA ANNi-PSO ANNi-ACO From the optimization results, the small error margins imply that the fan voltage values can be reproduced to allow higher drying velocities and reduce drying times. In addition, the reduced computation time will make it possible to take full advantage of the sampling times of the drying process, as well as define better strategies for the use of control algorithms that allow fine adjustment, not only to reduce drying times. In terms of drying velocity, the results suggest establishing drying ramps to achieve drying with specific characteristics for other agricultural products, as shown by the development for taro and plantain.
In conclusion, the hybrid strategy of ANNi applied to metaheuristic optimization obtained good precision with an acceptable error range. The computation times are reasonable for a drying system since the control and optimization strategies can be implemented online for the applied sampling time, establishing an automatic working regime for lowtemperature drying systems.

Conclusions
In this work, the design, instrumentation, experimental evaluation, and modeling of an indirect forced convection solar dryer with a chimney-type drying chamber were presented for its application in taro and plantain crops. A hybridization between artificial intelligence (artificial neural networks) and metaheuristic optimization techniques was established to obtain the operating conditions of the solar thermal system that maximize the drying velocity of the fruits. The system was evaluated under tropical climate conditions of the Mexican southeast by varying the fan voltage between 6 V and 12 V. The experimental results showed that, by applying a voltage of 9 V to the fan, the highest drying velocities were achieved for both fruits, which helped narrow down the search surface for the optimal answer.
From the point of view of both regional fruits, it was observed that, as the plantain has less porosity than the taro, it presented the best drying curves, taking less time to reach the critical moisture content. At the same time, the higher critical moisture content in taro allowed it to have a greater amount of water in its last drying stage, thus leading to higher drying velocities than in plantain.
An artificial neural network model was developed using the artificial intelligence application approach to estimate the drying velocity of both fruits, considering as independent variables the fan voltage, solar irradiance, ambient temperature, wind velocity, relative humidity, and the inlet and outlet temperatures of the solar dryer. The best network architecture consisted of eight hidden neurons, obtaining statistical values of RMSE of 0.1912, R 2 of 0.9647, R of 0.9822, and MAPE of 31.9%, for a hyperbolic tangent transfer function in the hidden layer. Although the aim was not to replace existing theoretical models, the results indicate that the ANN is an efficient alternative to synthesize the thermal phenomena inside the solar dryer and estimate the drying velocity with reduced error margins. The first contribution of the ANN model was to perform a sensitivity analysis to determine the degree of influence of the input variables on the drying velocity of the evaluated fruits and its relationship with the phenomenology of the drying process. Except for the wind velocity, all the variables had a significant impact on the drying velocity, and the fan voltage could be determined as the control variable allowing us to achieve the optimal theoretical response.
The second theoretical contribution of the ANN model was its use as an objective function to optimize the drying velocity. The application of the ANNi method allowed finding the optimal voltage for a given drying velocity. The ANNi was evaluated with three different types of metaheuristic optimization algorithms, obtaining the following hybridization results: genetic algorithm (ANNi-GA) with an error percentage of 0.83% and a computation time of 4.3 s, followed by particle swarm particle optimization (ANNi-PSO) with an error of 1.11% and a computation time of 8.7 s, and ant colony optimization (ANNi-ACO) with an error of 3.00% and 9.2 s of computation time.

Data Availability Statement:
The data used to support the findings of this study are available from the corresponding author upon request.

Appendix A
Algorithm for genetic algorithm (GA) optimization. move the best fitness to population g + 1 7 for i = 2 to Pobn do 8 Selection of two individuals 9 Apply crossover to create two new individuals 10 Apply mutation to the remaining population 11 end for 12 move new individual, mutate and crossed to population (g + 1) 13 end for Algorithm for particle swarm optimization (PSO). Table A2. Pseudocode and parameters used to develop particle swarm optimization [51].

Line Type Algorithm Simulation Parameters
Particle Swarm Parameter Value 1 Initialization for particle 1 to N Initialize the position x i (0) ∀ i ∈ 1: N Initialize particle best position to initial position p i (0) = x i (0) Calculate the fitness of each particle and, if ƒ(X j (0)) ≥ ƒ(X i (0)) ∀ i = j, initialize the global best as g = x j (0) Steps (n s ) 100 2 while stopping criteria is not reached Particles (N) 250 3 Update velocity for particle using the equation: v i (t + 1) = v i (t) + c i (p i − x i (t)) R 1 + c 2 (g − x(t)) R 2 Cognitive parameter (c 1 ) 1 4 Update the particle position using the equation: x i (t + 1) = x i (t) + v i (t + 1) Social parameter (c 2 ) 2 5 Evaluate fitness of the particle ƒ (X i (t + 1)) Minima Inertia weight (w 1 ) 0. 9 6 if ƒ (X i (t + 1)) ≥ ƒ(pi), update particular best: pi = xi (t + 1) Maxima Inertia weight (w 2 ) 0. 2 7 if ƒ (X i (t + 1)) ≥ ƒ(g i ), update global best: g i = x i (t + 1) 8 end. The best solution is presented at the end of the iterative process by g Algorithm for ant colony optimization (ACO).