Drying of Carrot Strips in Indirect Solar Dehydrator with Photovoltaic Cell and Thermal Energy Storage

: Minimizing agricultural losses, accompanied by employing technologies capable of taking advantage of solar energy, are the current challenges of the Antofagasta region of Chile due to having an average solar irradiance of 7.2 kWh/m 2 per day. With this objective, using an indirect solar dryer with storage of thermal energy in the form of sensible heat, the effect of drying on the quality of carrot strips was studied using chromatic coordinates CIEL*a*b*, the color difference ( ∆ E), the relationship between redness/yellowness (R = a*/b*), browning index (BI), whiteness index (WI), chroma (C), hue angle (h ◦ ) and drying kinetics. The experimental drying data were ﬁtted to 15 typically employed nonlinear regression models. The ∆ E = 14.11 ± 0.14 between the carrots in the dry and fresh conditions represented a detectable color change, the R ratio increased from 0.75 to 0.89, the BI index decreased from 209.82 ± 0.62 to 148.38 ± 0.26 and the WI index increased from 24.5 ± 0.11 to 31.8 ± 0.17, indicating color affectations due to the process. The coefﬁcients of determination, (R 2 ) close to 1 and the values closest to 0 of χ 2 , SSE and RMSE, led to the best ﬁt corresponding to the Weibull distribution model. In addition, it was found that the operation of the drying installation in the hours without incident solar radiation maintained the drying temperature values throughout the process.


Introduction
Among the 39 main fruits and vegetables, carrots occupy the tenth place in the ranking of nutritional value, being a good source of dietary fiber and trace element molybdenum, which is rarely found in vegetables [1]. Carrot is one of the important root crops cultivated worldwide, with annual production in 2017 of 42.8 million tons, reported as carrot + turnip. China has been the world's leading carrot-producing country for many years, with 45% of the global production. Other key carrot-producing countries are the United States and Ukraine. However, Chile reported a production of 168,484 tons in 2017 [2]. In many cases, carrot must be dried before its use. Many drying techniques have indeed been employed to dry carrots to maintain their appearance and their nutrients, especially β-carotene [3].
Vegetable cultivation in Chile stands out in second place; after corn, the carrot represents 2.8% of the hectares sown at the national level [4], being not only an attractive raw material for its designation of origin but also by its contribution in carotenoids and antioxidant capacity [5]. In addition, it is considered one of the most commonly used vegetables for human nutrition due to its high content of vitamins and fiber [6].
The agriculture of the Antofagasta Region in Chile, in addition to its potential contribution to regional food autonomy, constitutes the basis of the existence and reproduction of Atacama culture (Lickan Antay) since it is a form of agriculture but in the desert heights, associated with water scarcity and abundance of the sun [7]. The solar resource of the Safe and accurate estimates of the effective moisture diffusivity (D eff ) for food drying systems in solar drying conditions and the development of a simple methodology with general applicability for reliable simulation of the drying process are some of the major challenges faced by the researchers at present. Furthermore, the modeling and simulation of the mass transport phenomenon during solar drying of carrot slices has not been studied yet [27], as well as in the forms of cutting, cubes and strips.
With this objective, using an indirect solar dryer, the effect of drying on the color degradation of the carrot (Daucus carota L.) cut into longitudinal strips was studied, and the thermal behavior of the installation was also determined, which uses the recirculation of air and thermal storage in water based on sensible heat as characteristics that allow the use of energy for drying day and night.

Solar Drier Equipment and Complementary Equipment
The experimental equipment mainly consists of a pilot solar dryer of indirect operation using thermal energy storage in water and equipped with solar collectors (U-pipe) for converting the solar irradiance in heat energy and photovoltaic cells (monocrystalline) for converting the solar irradiance in electricity ( Figure 1). The description of the operation, the mode in which the air heating is performed and the complementary equipment that make up the solar drying installation are exhaustively described by Mezquita et al. [14]. (1) Three photovoltaic cells, (2) voltage controller, (3) electrical switches board, (4) battery, (5) heat exchanger air-water, (6) digital balance # 2, (7) drying chamber # 2, (8) window for entry and exit of air dryer, (9) drying chamber # 1; (10) digital balance # 1, (11), (12) solar collectors for water (U-pipe), (13) hydro-pneumatic tank water with inner membrane for thermal energy storage (with external insulation), (14) water pump, (15) variable speed drive for airflow, (16) water pump, (17) variable speed fan.

Experimental Procedure
Carrots used in the research were purchased from a central market of vegetables in the city of Antofagasta through a supplier from the community of Chiu-Chiu, a locality in El Loa Province, recognized as one of the main carrot producing areas in the region of Antofagasta, Chile.
Carrots were washed with water manually, rubbed with brushes, cut in strips (without removing the skin) 0.2 cm thick by 2 cm long and placed between 327 to 342 g of weight per tray, distributed evenly on three trays before starting the drying process; the thickness of the thin layer was 4 mm. During the drying period, the weight of the trays was determined with the strips of carrots every 20 min in each of the drying chambers using the coupled balances. The drying process was continued until constant weight. In all experiments, a constant airflow velocity of 1 m/s was set for the entire drying process. Experimental processes were performed in the spring (southern hemisphere) between September and October of 2016.
The initial moisture determination of the carrot sample was performed by drying on infrared balance, model ML-50 (A&D, Moisture Analyzer, Tokyo, Japan). The initial amount used in each analysis was 5 g, and it was not concluded until a constant weight was reached. The analyzes were carried out in triplicate at 100 • C.

Drying Curves and Mathematical Modeling Calculations for Making Drying Curves and Determining the Drying Ratios
For the preparation of the drying curves, a record of the data obtained experimentally was taken at constant time intervals (minutes), with which the moisture (g of water/g dry solid) of the product was calculated at each time, being able to obtain the graphs that relate the moisture loss concerning time. The moisture content (g water/g dry solid) was determined using the following equation [28]: where M: moisture content at any time (g water/g dry solid); W o : initial weight of the sample (g); W: the amount of evaporated water (g); W 1 : dry matter content of the sample (g).
The moisture ratio (MR) was simplified to M M 0 , (where M o is the content of the initial moisture at time = 0) instead because, in solar drying, the sample can be dried until very close to the equilibrium condition and reach approximately the content of the dry matter [28]. In such a way, the measurement of the equilibrium moisture content is not necessary (M e ).
Different mathematical models have been proposed to describe the characteristics of agricultural products (fruits and vegetables) during the conventional drying process and in solar dryers [14,29]. All the equations applied in this study to model the drying kinetics of carrot strips are shown in Table 1. Three of these models are: "models derived from Newton's law of cooling"; nine are: "models derived from Fick's second law of diffusion"; three are: "empirical models" [29].

Color Analysis
Color determination of the fresh and dried samples was measured by a Hunter laboratory tristimulus colorimeter (Color Flex, Hunter Lab Inc., Reston, VA, USA). The instrument, equipped with D 65 illuminant and 10 • observer optical position, was calibrated using black and white reference tiles provided. All measurements were made using the CIEL*a*b* system. The parameters L*, a* and b* indicate whiteness/brightness, redness/greenness and yellowness/blueness, respectively. For each sample, at least ten measurements were made at different positions of the sample, and each of the measurements were made in triplicate. The color change was determined by color difference (Equation (2)), the chroma (Equation (3)) and Hue angle (Equation (4)) [37].
where ∆L * = L * 2 − L * 1 : Luminosity difference; ∆a * = a * 2 − a * 1 : Redness difference; The subscripts 1 and 2 were used to designate the conditions (initial and final) in which the colors were determined, establishing the color differences between both points. The relationship of redness on yellowness was calculated as follows (Equation (5)) [38]: where: a* and b* are chromatic coordinates. MR: moisture ratio; k: kinetical parameter (min −1 ); n and a: empirical parameters (dimensionless); t: drying time (min); b, c, g, h, k o and k 1 : empirical constants in the drying models [27].
A color parameter that was also calculated is that which represents the purity of brown color (brown) and is reported as an important parameter in the processes in which enzymatic and non-enzymatic browning take place [39]. This calculation is performed mainly in the drying processes to indicate degradative aspects of the color, such as the browning Sustainability 2022, 14, 2147 6 of 18 index (BI) (dimensionless) and whiteness index (WI) (dimensionless) [40], represented by the Equations (6)- (8): where x is:

Calculation of Effective Diffusivity Coefficient
Fick's second law was used to determine the effective diffusivity coefficient, which is adequate for diffusive phenomena during the mass transfer in the food during drying until reaching equilibrium. Fick's second law, symbolized as a mass-diffusion equation for drying of agricultural products in a falling rate period, is shown in the following Equation (9): The analytical solutions of Fick's second law (Equation (9)) for infinite slab can be given as Equation (10) with the assumptions of neglecting shrinkage, constant temperature, diffusion coefficients and uniform initial moisture distribution [41].
where D eff is the effective moisture diffusivity in m 2 /s, t is the time (h), n is a positive integer and L is the half-thickness of samples (m). For long drying periods, Equation (10) can be further simplified to only the first term of the series. Equation (11) is written in a logarithmic form as follows: The effective moisture diffusivity is calculated from the slope of a straight line by plotting data in terms of ln (MR) versus drying time, for which slope is K (Equation (12)):

Statistical Analysis
The fit of the best mathematical model that describes the thin layer drying curve of carrot strips samples was done using the 15 models presented in Table 1, employing the nonlinear regression analysis through the Statgraphics-Centurion XVI.1 software (StatPoint Technologies Inc., Warrenton, VA, USA).
The determination coefficient (R 2 ) (Equation (13)) [17] was one of the main criteria to evaluate the fit quality of proposed mathematical models, complemented by the calculation of the statistical parameters: root means square error (RMSE) (Equation (14)) [17], sum square errors (SSE) (Equation (15)) [17] and the reduced chi-square (χ 2 ) (Equation (16)) [17]. These parameters were calculated with the following equations: Sustainability 2022, 14, 2147 where MR exp,i is the i th experimentally observed moisture ratio found in any measurement and MR pre,i is the i th predicted moisture ratio for the measurement under consideration. N and n are the number of observations and the number of constants, respectively. The statistical analysis of experimental data was determined by applying an analysis of variance (ANOVA) to estimate any statistically significant difference at a confidence level of 95% (p < 0.05). Differences between mean values were established using Duncan's multiple range tests. The highest values of R 2 (≈1.0), along with the closeness of RMSE, SSE and χ 2 values to zero, were selected as optimal criteria to evaluate the fit quality of the proposed drying models. All data were presented as mean values with standard deviations. All solar drying experiments were replicated thrice under similar conditions.

Drying Curve
The initial moisture content of the carrot samples was 9.07 g of water/g of dry solid, corresponding to 9.93 ± 0.08% of the total solids, and, in that condition, they were subjected to the solar drying process until the weight losses were made practically constant. This moisture content in the fresh carrot, an ecotype of the community of Chiu-Chiu, is within the ranges of percentages reported in different investigations, which are between 87.0 and 93.0% (bh) [6,21,42].
In Figure 2, the drying curve, moisture ratio M M 0 vs. drying process time (min), is presented, in which a total average time of the drying processes of 1600 min (26.7 h) was recorded. In addition, three stages stand out, the first consisting of the induction period, called "stage A", of 80 min (1.3 h), the second consisting of the constant drying rate period, called "stage B", of 380 min (6.3 h), which was defined through the best fit of the straight line shown in Figure 2 with R 2 = 0.9941, between the end of the induction period and the beginning of the third stage, which is made up of the falling drying rate period, called "stage C", of 1140 min (19 h).
Regarding the hours of incidence of solar irradiance, up to 591 min (≈10 h) profited on the first day of drying, and 200 min (3.3 h) on the second day of drying, equivalent to 47.4% of the total process time, where the maximum average temperatures of 56.1 • C (air in chamber 1) were reached ( Figure S1 and Table S1; Supplementary Materials), with an average relative humidity of 43.4 ± 1.5%.
Due to the thermal storage of heat carrier liquid (distilled water) in the drying installation, it was possible to continue with the drying process in the night period without the influence of solar radiation, which represented 46.1% of the total time of the drying process, 770 min (12.83 h), achieving the carrot strips moisture (bh) of 20% at end of the night. However, the final moisture content reached by the carrot strips was 0.131 ± 0.001 g of water/g of dry solid, which corresponded to 11.56% moisture (bh), a value lower than 14%, which is recommended for a dry agricultural product that remains without deterioration [43] and is within the range of 5.7 to 22.5%, which are the optimal values to ensure microbiological stability. Urrea et al. [44] reported that the loss of carotenoids depends on the drying temperature to which the carrot samples are exposed. In the carotenoids, where β-carotene, the main component of carrots, is found, it is usually in a concentration of 45 to 80% [45], and they are relatively stable when drying takes place at low temperatures (40 to 60 °C) but are very sensitive to higher temperatures (70 to 90 °C) [46,47]. Moreover, Sun et al. [48] reported, using a heat pump assisted convective dryer, and Zielinska & Markowski [49], employing a fluidized bed dryer, concluded that the highest retention of β-carotene in carrot drying was observed at temperatures of 40 °C, 50 °C and 60 °C, reaching 81.8%, 77.9% and 73.0%, respectively. This indicated that the retention of β-carotene was largely dependent on the drying temperature. Due to this approach, it is to be expected that, in the present investigation, the sensitive components have not been significantly altered during the indirect solar drying process by maintaining an average temperature of 37.2 °C with a maximum of 56.1 °C during the process, which can be observed in Figure S1 and Table S1 (Supplementary Materials) in chambers 1 and 2.

Drying Rate Curve
The drying rate curve observed in Figure 3 shows a typical curve of a convective drying process with its induction period (A), constant drying rate period (B) and falling drying rate period (C), similar to those obtained by Doymaz [27] for convective hot-air drying carrot slices. The maximum temperature of the drying air was obtained in the period of constant drying speed, being 56.1 and 40.3 °C in drying chambers 1 and 2, respectively. In this period, the highest average value of the drying speed was also reached, being 0.023 ± 0.001 g water/g dry solid * minute. This period can be seen in Figure S1 (Supplementary Materials) between 11:50 and 18:10 h (380 min), this last hour being where the critical moisture (Mc) of the process was established, 2.01 g of water/g of dry solid (Figure 3), which is defined by the end of the constant drying rate period and the start of the falling drying rate period [50], showing the fitted straight line MR = −0.0019 Dt + 1.1071 (R 2 = 0.9941) of the constant rate period (Figures 2 and 3).
Critical moisture (Mc) is not a property of the food but depends on the ease of movement of moisture through the solid, the porosity of the food structure [51] and the drying Urrea et al. [44] reported that the loss of carotenoids depends on the drying temperature to which the carrot samples are exposed. In the carotenoids, where β-carotene, the main component of carrots, is found, it is usually in a concentration of 45 to 80% [45], and they are relatively stable when drying takes place at low temperatures (40 to 60 • C) but are very sensitive to higher temperatures (70 to 90 • C) [46,47]. Moreover, Sun et al. [48] reported, using a heat pump assisted convective dryer, and Zielinska & Markowski [49], employing a fluidized bed dryer, concluded that the highest retention of β-carotene in carrot drying was observed at temperatures of 40 • C, 50 • C and 60 • C, reaching 81.8%, 77.9% and 73.0%, respectively. This indicated that the retention of β-carotene was largely dependent on the drying temperature. Due to this approach, it is to be expected that, in the present investigation, the sensitive components have not been significantly altered during the indirect solar drying process by maintaining an average temperature of 37.2 • C with a maximum of 56.1 • C during the process, which can be observed in Figure S1 and Table S1 (Supplementary Materials) in chambers 1 and 2.

Drying Rate Curve
The drying rate curve observed in Figure 3 shows a typical curve of a convective drying process with its induction period (A), constant drying rate period (B) and falling drying rate period (C), similar to those obtained by Doymaz [27] for convective hot-air drying carrot slices. The maximum temperature of the drying air was obtained in the period of constant drying speed, being 56.1 and 40.3 • C in drying chambers 1 and 2, respectively. In this period, the highest average value of the drying speed was also reached, being 0.023 ± 0.001 g water/g dry solid * minute. This period can be seen in Figure S1 (Supplementary Materials) between 11:50 and 18:10 h (380 min), this last hour being where the critical moisture (Mc) of the process was established, 2.01 g of water/g of dry solid (Figure 3), which is defined by the end of the constant drying rate period and the start of the falling drying rate period [50], showing the fitted straight line MR = −0.0019 Dt + 1.1071 (R 2 = 0.9941) of the constant rate period (Figures 2 and 3). rate period began, and the drying rate progressively decreased until equilibrium, reaching mean values of 0.00045 g water/g dry solid*min and final moisture, called equilibrium moisture (Meq), with mean values of 0.129 g of water/g of dry solid (Figure 3).
In other studies performed with carrots in the form of slices and/or cubes using convective solar drying, similar results have been obtained both in the behavior of the drying curves and in the drying rate curves [27,53,54].  Table 2 shows the results of the statistical parameters, R 2 , χ 2 , SSE and RMSE, which were used to determine the goodness of fit of the drying curves of the fifteen models studied in carrot strips [33]. The values closer to 1 for R 2 and closest to 0 for χ 2 , SSE and RMSE were considered optimal criteria to assess the quality of fit of the models [55]. Thus, the best fit model, according to the thermal history displayed in the process drying curve for the group of the three models derived from Newton's law of cooling, Newton, Page and modified Page, was the Page model with R 2 , χ 2 , SSE and RMSE values of 0.9971, 4.0 × 10 −4 , 3.7 × 10 −4 and 1.936 × 10 −2 , respectively. The values of the kinetic constant (k) and the empirical parameter (n) were: 1.836 × 10 −5 min −1 and 1.83288, respectively.

Modeling of Drying Curves Statistical Parameters
Amongst the tested models for the group of the nine models derived from Fick's second law of diffusion, Henderson and Pabis, modified Henderson and Pabis, logarithm, approach of diffusion, Midilli-Kucuk, two-term, two-term exponential, Aghbashlo and Critical moisture (M c ) is not a property of the food but depends on the ease of movement of moisture through the solid, the porosity of the food structure [51] and the drying process conditions. In a forced-air convective drying study using cubed carrots, M c values were obtained between 2.4 and 3.9 g water/g dry solid [52].
In Figure 3, the induction period began with an initial moisture content value (M o ) of 9.07 g of water/g of dry solid until reaching a moisture value of 8.32 g of water/g of dry solid at 80 min. From there, the constant drying rate period began where a drying rate of 0.010 to 0.023 g of water/g of dry solid * min is reached for a mean value of 0.0165 g of water/g of dry solid. This variation is due to fluctuations in air temperatures in the drying chambers, which is logical for this type of indirect solar dryer that depends on the air-water heat exchange that occurs when the constant temperature in the heat carrier liquid is not maintained. At this stage, the air temperature in drying chambers 1 and 2 reached values of 30.4 • C and 22.4 • C for the beginning and 55.5 • C and 40.3 • C for the end, respectively ( Figure S1 in Supplementary Materials). Once the M c was reached, the falling drying rate period began, and the drying rate progressively decreased until equilibrium, reaching mean values of 0.00045 g water/g dry solid*min and final moisture, called equilibrium moisture (M eq ), with mean values of 0.129 g of water/g of dry solid (Figure 3).
In other studies performed with carrots in the form of slices and/or cubes using convective solar drying, similar results have been obtained both in the behavior of the drying curves and in the drying rate curves [27,53,54]. Table 2 shows the results of the statistical parameters, R 2 , χ 2 , SSE and RMSE, which were used to determine the goodness of fit of the drying curves of the fifteen models studied in carrot strips [33]. The values closer to 1 for R 2 and closest to 0 for χ 2 , SSE and RMSE were considered optimal criteria to assess the quality of fit of the models [55]. Thus, the best fit model, according to the thermal history displayed in the process drying curve for the group of the three models derived from Newton's law of cooling, Newton, Page and modified Page, was the Page model with R 2 , χ 2 , SSE and RMSE values of 0.9971, 4.0 × 10 −4 , 3.7 × 10 −4 and 1.936 × 10 −2 , respectively. The values of the kinetic constant (k) and the empirical parameter (n) were: 1.836 × 10 −5 min −1 and 1.83288, respectively. Amongst the tested models for the group of the nine models derived from Fick's second law of diffusion, Henderson and Pabis, modified Henderson and Pabis, logarithm, approach of diffusion, Midilli-Kucuk, two-term, two-term exponential, Aghbashlo and Verma models (Table 2), the one with the best fit was the Midilli-Kucuk model with R 2 , χ 2 , SSE and RMSE values of 0.9989, 2.1 × 10 −4 , 1.9 × 10 −4 and 1.376 × 10 −2 , respectively; with empirical constants a, b and n of 0.967416, 1.3 × 10 −5 , 2.047330 and the kinetic constant (k) of 5.0 × 10 −6 min −1 , respectively.

Modeling of Drying Curves Statistical Parameters
Of the three tested models that belong to the group of so-called empirical models, the one that resulted in the best fit was the Weibull distribution model with values to R 2 , χ 2 , SSE and RMSE of 0.9988, 1.7 × 10 −4 , 1.6 × 10 −4 and 1.252 × 10 −2 , respectively; being the empirical constants a, b and n of 1.812 × 10 −2 , −0.95691, 1.9768 and the kinetic constant (k) of 7.87 × 10 −6 min −1 , respectively ( Table 2). Figure 4 shows the process curves made with the experimental values versus the predicted values from the three models that presented R 2 ≥ 0.997 ( Table 2). The moisture ratio (MR) values for the predicted drying time followed a trend in a straight line at 45 • , which showed the correlation of the data fitting to the best models that could predict the drying kinetics of carrot strips in the indirect solar drying system assisted with U-pipe solar collectors and thermal energy storage in the form of sensible heat, with a maximum air velocity of 1 m/s through in the drying chamber.  Figure 4 shows the process curves made with the experimental values versus the predicted values from the three models that presented R 2 ≥ 0.997 ( Table 2). The moisture ratio (MR) values for the predicted drying time followed a trend in a straight line at 45°, which showed the correlation of the data fitting to the best models that could predict the drying kinetics of carrot strips in the indirect solar drying system assisted with U-pipe solar collectors and thermal energy storage in the form of sensible heat, with a maximum air velocity of 1 m/s through in the drying chamber.   Table 2, those proposed by Page, Midilli-Kucuk and the Weibull distribution presented the highest values of determination coefficients (R 2 ), 0.9971, 0.9989 and 0.9988, representing the first two, the category of so-called semi-theoretical models and, finally, empirical models [29]. However, once the other statistical parameters (SSE, RMSE and χ 2 ) were analyzed, a good quality of fit was observed based on this evaluation since values of SSE = 1.6 × 10 −4 , RMSE = 0.01252 and χ 2 = 1.7 × 10 −4 were obtained through the Weibull distribution model closer to zero than the rest of the models according to a proposal by Erbay & Icier [29]. Therefore, the following Weibull Equation (17)

Model Validation
Although this model can be used to predict changes in moisture content over time, it does not have the combined effect of the other drying parameters, such as airspeed and drying temperature, since the present study is limited only to air velocity of 1 m/s and variable drying temperatures, which did not exceed 56.1 • C.
The Midilli-Kucuk model presented the best fit in investigations of thin-layer drying of carrot slices using the infrared dryer [56], by electric resistance convective dryer [17] and with ultrasound pretreatment [57]. On the other hand, for drying carrots in cubes and/or slices [6,42] or carrot pomace [28] employing a convective dryer, the best fit model corresponded to Page's model.

Chromatic Coordinates
The color of food represents a very important quality and acceptability parameter and is within the optical properties to be evaluated in dehydrated products [58]. In addition, it is a nutritional index valuation of carrots by significantly relating this attribute to the content of β-carotene determined [59]. Color affects consumers' preferences when choosing food products. Due to this, it is suggested that it is interrelated with other properties of food and can represent them. Additionally, color is controlled by reactions that occur throughout the supply chain. For these reasons, food industries and many researchers use this attribute to measure product quality [40]. Table 3 shows the values of the chromatic coordinates L*, a* and b* of the carrot strips, before and after the drying process, indicating that there were statistically significant differences (p < 0.05) analyzed by Duncan's multiple range test. The value of L* indicates that. once the drying process was completed, the carrot strips acquired a darker luminosity, with values that were 9.4% lower than the value of the state before drying. This condition can be attributed to a loss of brightness due to the effect of the measurement method since, as Nowacka & Wedzic [60] indicate, the instrumental concept is based on measuring the radiation reflected from the surface, and. in this case, the dry sample compared to the fresh sample has less water and stops reflecting the light that falls on the sample. This result can be compared with those obtained by Huang et al. [61] and Kowalski et al. [62], who observed a slight decrease in this chromatic coordinate when dehydrating carrots.
The results for the a* values showed a decrease in this chromatic coordinate at the end of the drying process, which indicated a red color change of 14.4%. This could be explained due to the degradation of the carotenoids in carrot tissue, which affects the color of the product [63]. A more pronounced decrease showed the b* values, which indicated that the fresh sample is 27.6% more yellow than the dry one. This behavior agrees with the results obtained by Demiray & Tulek [64] when drying sliced carrots in a cabin dryer with hot air, who indicated that the decrease in the value of the chromatic coordinate b* can be referred to as the degradation of carotenoid pigments [3,65], non-enzymatic Maillard browning and brown pigment formation [59]. Table 3. Chromatic coordinates of the carrot strips in their fresh state (before drying) (n = 10) and after the drying process (n = 10). Color difference (∆E), ratio of redness over yellowness (R), browning index (BI) and whiteness index (WI) (n = 3). The color differences ∆E = 14.11 ± 0.14 shown in Table 3 represent a perceptible color change to human vision since ∆E ≥ 5.0. Obón et al. [66] report that a difference in color from 0.0 ≤ ∆E ≤ 1.5 can be considered small and almost identical for visual observation; in the range from 1.5 ≤ ∆E ≤ 5.0, the color difference can be distinguished, while the color difference is evident for ∆E > 5. On the other hand, Hong et al. [67] state that it is generally known that values of ∆E < 3.0 cannot be easily detected by the human eye. This color change can be attributed to the degradation of compounds, long drying times and/or exposure to high drying temperatures [68]; however, this value is considered to be expected when compared with the results obtained for hot air-drying carrot slices of ∆E = 14.43 reported by Kowalski et al. [62]. Similarly, Liu et al. [63] reported values of 12.5 ≤ ∆E ≤ 22.5 when dehydrating carrots by convective drying with oxygen levels of 20.9% in the circulating air of the drying chamber. The other color parameters given by h • and C suffered losses of 9.1% and 23.0%, respectively. When the values of dehydrated carrots were compared with the samples before the drying process, with the decrease in hue, the saturation also decreased. For its part, R, the redness/yellowness ratio (a*/b* value), which is an index of great applicability when studying possible incidences of browning in food products, increased its ratio [69], going from the initial mean values for the fresh sample from R = 0.75 ± 0.01 to the dry state of R = 0.89 ± 0.01.

Chromatic Coordinates and Other Parameters
Non-enzymatic browning is a quality criterion for dried foods. It is caused by reactions between sugars and proteins in the form of either free amino acids or amino groups of the lysine residues as a result of high-temperature drying. The browning occurrence can lead to dark pigments and decreased nutritive value of amino acids and proteins and other undesirable changes in dried product's flavors and texture [70]. That is why non-enzymatic browning is another quality indicator in the drying of carrots. The extent of browning is mainly attributed to the color changes resulting in Maillard reactions in the carrot. The application of solar drying reduces the BI values and increases the WI values of the carrot strips at the end of the drying process, BI = 148.38 ± 0.26 and WI = 31.8 ± 0.17, compared to their initial values, BI = 209.82 ± 0.62 and WI = 24.5 ± 0.11, indicating that there was a degradation in compounds, such as sugars and organic acids. The reduction in BI values and increase in WI values samples show that the rate of brown pigment formation reduced with drying application [71]. Moreover, Wiktor & Witrowa-Rajchert [72] suggest that, in the case of carrots, nonenzymatic browning may occur in the drying process, for instance, due to the reactions between carbonylic molecules, which are the products of carotenoids oxidation, with other food constituents.

Estimation of the Effective Moisture Diffusivity (D eff )
Deff was calculated using the slope K from the representation of the natural logarithm of the moisture ratio concerning the time of drying process (( Figure 5), Equations (11) and (12)), where D eff is effective moisture diffusivity (m 2 /s) and L is the thickness dimension of the thin layer in which the carrot was placed in strips on the drying trays in (m). This value of Deff = 1.82 × 10 −8 m 2 /s was obtained for low temperatures of the drying chambers in the 26.6 h of the drying process and did not exceed 56.1 °C ± 0.4 as the maximum value of temperature reached by air in drying chamber 1, which agrees with Mahapatra & Tripathy [27], who studied the dried process of carrot slices of 0.03 m diameter and 0.005 m thickness in direct, indirect and mixed-mode type natural convection solar dryers without thermal energy storage and obtained results for Deff ranging from 2.59 × 10 −8 to 6.36 × 10 −8 m 2 /s, respectively, with the maximum values obtained for mixed-mode solar-dried samples followed by direct and indirect solar-dried samples. In other investigations, Doymaz [7] reported Deff values from 0.776 to 9.335 × 10 −9 m 2 /s for temperatures from 50 to 70 °C of carrots cut into cubes of 1 cm × 1 cm × 1 cm and 2 cm × 2 cm × 2 cm. In general, the time required to reduce the moisture ratio to any given level was dependent on the drying condition, being the highest at 50 °C and lowest 70 °C. With the drying process, the time needed to reduce the moisture content of the carrot cubes (d = 2 cm) from the initial 87.5% to the final 6% was 21, 16, 14 and 11 hours at 50, 60, 65 and 70 °C, respectively. The drying times of the carrot cubes (d = 1 cm) at 50 and 60 °C were 13.5 and 9.5 h, respectively. On the other hand, Aghbashlo et al. [73] carried out the drying process in a semi-industrial continuous belt dryer for carrot slices with a thickness of 5 mm, using, for the drying experiments, three air temperatures, 50, 60 and 70 °C, three speeds of the air, 0.5, 1.0 and 1.5 m/s, and Deff varied from 3.21 × 10 -7 to 8.98 × 10 -7 m 2 /s. Through another type of dryer, Hu et al. [74] studied the dehydration kinetics of previously treated carrot slices to prevent color loss using microwave-assisted fluidized bed drying at drying air temperatures of 40, 50 and 60 °C and evaluated the microwave power density between 0.44 and 1.21 W/g, reaching Deff values in the range of 0.720 × 10 -8 to 1.890 × 10 −8 m 2 /s.

Conclusions
The results of the present study led to the conclusion that the dryer that was used, of an indirect type, with storage of thermal energy in the form of sensible heat and with partial recovery of the exhaust air, allowed a continuous drying process of the carrot The mean value of D eff obtained for the present indirect solar dryer with a photovoltaic cell and thermal energy storage as a function of the dehydration of the carrot strips and an airspeed of 1 m/s was 1.82 × 10 −8 m 2 /s. This value agrees with the general ranges indicated for thin-layer drying characteristics, where the D eff for carrots, using methods of drying by ultrasound and vacuum at 65 • C, reached values of D eff 5.44 × 10 −8 and 2.85 × 10 −8 m 2 /s, respectively [68].
This value of D eff = 1.82 × 10 −8 m 2 /s was obtained for low temperatures of the drying chambers in the 26.6 h of the drying process and did not exceed 56.1 • C ± 0.4 as the maximum value of temperature reached by air in drying chamber 1, which agrees with Mahapatra & Tripathy [27], who studied the dried process of carrot slices of 0.03 m diameter and 0.005 m thickness in direct, indirect and mixed-mode type natural convection solar dryers without thermal energy storage and obtained results for D eff ranging from 2.59 × 10 −8 to 6.36 × 10 −8 m 2 /s, respectively, with the maximum values obtained for mixed-mode solar-dried samples followed by direct and indirect solar-dried samples. In other investigations, Doymaz [7] reported D eff values from 0.776 to 9.335 × 10 −9 m 2 /s for temperatures from 50 to 70 • C of carrots cut into cubes of 1 cm × 1 cm × 1 cm and 2 cm × 2 cm × 2 cm. In general, the time required to reduce the moisture ratio to any given level was dependent on the drying condition, being the highest at 50 • C and lowest 70 • C. With the drying process, the time needed to reduce the moisture content of the carrot cubes (d = 2 cm) from the initial 87.5% to the final 6% was 21, 16, 14 and 11 hours at 50, 60, 65 and 70 • C, respectively. The drying times of the carrot cubes (d = 1 cm) at 50 and 60 • C were 13.5 and 9.5 h, respectively. On the other hand, Aghbashlo et al. [73] carried out the drying process in a semi-industrial continuous belt dryer for carrot slices with a thickness of 5 mm, using, for the drying experiments, three air temperatures, 50, 60 and 70 • C, three speeds of the air, 0.5, 1.0 and 1.5 m/s, and D eff varied from 3.21 × 10 -7 to 8.98 × 10 -7 m 2 /s. Through another type of dryer, Hu et al. [74] studied the dehydration kinetics of previously treated carrot slices to prevent color loss using microwave-assisted fluidized bed drying at drying air temperatures of 40, 50 and 60 • C and evaluated the microwave power density between 0.44 and 1.21 W/g, reaching D eff values in the range of 0.720 × 10 -8 to 1.890 × 10 −8 m 2 /s.

Conclusions
The results of the present study led to the conclusion that the dryer that was used, of an indirect type, with storage of thermal energy in the form of sensible heat and with partial recovery of the exhaust air, allowed a continuous drying process of the carrot strips, verifying good operation in the hours without incidence of solar radiation, ending the night period with a difference of 9.8 • C concerning the ambient temperature after 26.6 h of processing. This total process time is relatively shorter considering that, in other types of indirect convective solar dryers, without thermal energy storage, the times would be longer. Of the 15 thin-film drying models studied, only the Weibull distribution model had the best fit according to the established statistical parameters. Moreover, the Page and Midilli-Kucuk models cannot be discarded since they had high adjustments and can be considered satisfactory in the investigation. The dried carrot strips' quality for this type of process was as expected, although they presented color differences concerning the fresh product, as well as different evaluated colorimetric parameters, such as tone angle, chroma, BI, WI and R. In the present investigation, the energy savings due to the non-use of fossil fuels in the drying installation were not quantified.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/su14042147/s1/ Solar irradiance and thermal history of the drying vehicles, water and air in the solar installation [75,76]. Figure S1. Variation of solar irradiance and temperatures in the stages of the solar drying process; Table S1. Average values of temperature and solar irradiance determined in different parts of the solar drying installation.