3.1. Drying Characteristics
The initial moisture content of the apple slices was 4.25 ± 0.124 (d.b). The drying rate exhibited higher values during the middle phase of the drying process across all microwave power settings (
Figure 3a–c). Moreover, there was a notable enhancement in the drying rate as the microwave power was increased. The elevated drying rate during the middle phase resulted from the rapid mass transfer rate and the evaporation of a substantial volume of free moisture from the apple’s surface. Furthermore, the drying rates were high at higher microwave power due to accelerated evaporation, promoting moisture diffusion from the interior to the exterior [
7].
Also, PEF pretreatment had a significant effect on increasing the drying rate for all drying conditions. The application of PEF treatment led to tissue damage and enhanced the permeability of cell membranes. This led us to propose the hypothesis that following PEF treatment, there would be an acceleration in the drying rate during the middle drying phase due to the expedited moisture supply to the surface. The influence of the PEF treatment on the drying rate of the current research was in agreement with [
16,
19,
44].
The impact of pretreatment with PEF on the drying rate and changes in spinach quality during hot air drying was investigated [
16]. They reported the drying times for the PEF samples were shortened to 1.5 h by air drying. They supposed that the moisture supply to the surface was accelerated after PEF treatment, so the drying rates were increased in the early stage of drying. Also, the influence of PEF pretreatment on the convective drying kinetics of onions was assessed by [
19]. The drying time and the diffusion coefficient for the onion in this study were optimal at a medium PEF treatment intensity (4–6 kJ/kg), giving a Z value in the range of 0.53–0.60. They reported that applying PEF treatment resulted in an enhancement of the drying rate throughout the drying process. The sharp incline observed in the middle phase of the drying curve also indicated that liquid moisture diffusion was the primary mechanism for eliminating water from the sample’s interior.
Figure 3d–f illustrates the correlation between the moisture ratio of the samples as MR versus drying time progresses under different pretreatment and drying conditions. It is evident that as the drying time increased, the moisture ratio of the apple consistently decreased. The moisture ratio dropped rapidly at the beginning of the drying process and a gradual decrease was noted during the concluding phases of the drying procedure. The difference in water vapor pressure between the interior and surface of the apple mass accelerates the internal movement of moisture toward the surface. This phenomenon explains the decrease in moisture content.
As the microwave power increased, the amount of moisture decreased considerably (
Figure 3d–f). For example, the required time to reach MR of 0.1 with powers of 100, 200 and 300 W were 54, 36 and 30 min, respectively. As the microwave power increased, there were improvements in the temperature gradients and surface evaporation rates, leading to an acceleration in the diffusion of moisture from the interior to the surface. These outcomes align with the results of other studies that linked higher microwave power to reduced drying times [
42,
45].
Also, PEF treatment reduced the rate of MR significantly in all experimental conditions. However, increasing the frequency from 100 to 300 Hz (which enhanced the specific energy) had no significant effect. For instance, the MR of control, PEF-100, PEF-200 and PEF-300 were 0.29 ± 0.02, 0.12 ± 0, 0.14 ± 0.01, 0.12 ± 0, respectively, after 21 min (
Figure 3d–f). Based on the findings, implementing PEF reduced the drying time from 4.2 to 31.4% compared to the untreated sample. The results could be elucidated by releasing a significant amount of free water from the samples subjected to PEF electroporation. The released water from the electroporated apple cells can move unrestrictedly and be efficiently removed during the drying procedure [
13,
30]. Similar behavior of PEF treatment on the MR trends over the drying time was reported for onion [
19], carrot [
3], spinach [
16] and potato [
18]. The effects of PEF on the vacuum drying and quality characteristics of dried carrot were evaluated by [
3] under different drying temperatures (25 °C, 50 °C, 75 °C and 90 °C). When the temperature increased from 25 to 90 °C, the total drying time decreased from 21,600 s to 5400 s and from 9720 s to 3600 s for untreated and PEF-pretreated samples, respectively. The application of PEF treatment caused a noticeable decrease in the drying time (by 33–55%) and acceleration drying kinetics even at a higher temperature (90 °C); it evidently reflected the effects of electroporation by the PEF treatment.
3.2. Mathematical Modeling of Drying
The alteration in moisture content throughout the MV drying of apple slices was assessed by fitting the Page and Weibull models and, subsequently, the best model was employed to depict the drying behavior of the sample (
Table 2). The applied mathematical model demonstrated an acceptable fit with the experimental data and exhibited R, RMSE and MAE values in the range of 0.917–0.977, 0.011–0.075 and 0.006–0.046, respectively. Based on the statistical parameters, it was recommended that the Weibull model yielded an average R
2 of 0.963 to effectively characterize the moisture loss process in apple slices during the MV drying. The recent literature has confirmed the robustness of the Weibull model in accurately fitting the kinetic drying profiles of diverse products exposed to varying drying conditions [
14,
25].
Although the models successfully created the correlation between average moisture ratio and drying time, they lack consideration for the underlying principles of the drying process, such as PEF, resulting in their parameters lacking physical significance. Consequently, understanding the kinetics of drying enables the comprehension and prediction of drying times, which ultimately lead to enhanced process efficiency through optimization [
14]. The kinetic constants of each model under different experimental conditions are presented in
Table 2. The kinetic constant of the Page model revealed a notable increase as microwave power intensified the PEF treatment. The lowest and highest kinetic constants of the Page model were observed in the control-100 W (1.241) and PEF3-300 W (3.931) cases, respectively. Conversely, the kinetic constant of the Weibull model decreased with the augmentation of microwave power and PEF treatment. The influence of PEF treatment on the upward and downward trends of Page’s and Weibull’s constant values in this study was similar to the findings of [
14,
46].
To assess the relationship between the kinetic parameter and microwave power, the kinetic constants were graphed against the power and a linear regression was computed (
Figure 4). A linear correlation between the kinetic parameter and microwave power was observed. The linear regression models offer the potential to predict the drying rate at a power ranging from 100 to 300 W. The maximum calculated correlation coefficients by Page’s and Weibull’s models were 0.962 and 0.986, respectively. The pretreatments led to a more pronounced decrease in drying time at 300 W, whereas the escalation in the specific energy of PEF did not facilitate any further reduction. This pattern might be associated with the alteration of the initial water content due to pretreatments, which could potentially augment the mass transfer rate through the application of PEF.
Although the overall accuracy and errors of the Page and Weibull models were acceptable, it appears that all the models assumed a form of drying kinetics that may not capture all the complexities of the drying process. Thus, it should be combined with other additional empirical models. Further, the models have the assumption that the drying characteristics are homogeneous, which may not be the case for sliced apple samples; thus, this model could prove to be useful for materials with relatively uniform properties or perhaps be integrated together with other models accounting for heterogeneity.
3.3. Chemical Analysis and Chromatographic Method
The results obtained from HPLC analyses show a wide range of chemical compounds constituting the extract (
Figure 5a). It is noticeable that the HPLC method was performed for using its data as an input dataset to increase the accuracy of the ML model. The antioxidant activity, assessed using DPPH radicals, ranged from 0.98 to 1.24 mg TE/100 g dry weight for the obtained dried apple slices (
Figure 5b). There was a notable decrease in antioxidant activity with increasing microwave energy. However, pretreatment with PEF showed a positive impact on reducing antioxidant levels. For instance, with PEF1, the antioxidant content at an energy level of 200 W was ten units higher and, intriguingly, as the energy level of PEF increased, the antioxidant levels remained consistent. Nevertheless, no significant difference was observed with PEF2 and PEF3. During treatment with PEF, it can potentially disrupt cell membranes and eventually lead to the release of intracellular components. The disruption may result in the degradation or oxidation of certain sensitive antioxidants during the process when otherwise protected in cell structures [
47]. Similarly, the total polyphenol content exhibited a notable decrease with increasing microwave energy, with a significant distinction observed between PEF1 and PEF2 at 100 W (
Figure 5c). However, there was no significant contrast between the PEF treatments when the device operated at 200 and 300 W. Polyphenols are very sensitive to heat and drying is carried out at high temperatures, so the degradation of polyphenols can easily occur by oxidation or simple thermal degradation. The applied PEF may enhance the cell membrane permeability, facilitating easy access to polyphenols and, hence, their bioavailability.
Also, both the level of microwave energy and the application of PEF treatment influenced the total flavonoid (
Figure 5d) and vitamin C contents (
Figure 5e). As the microwave energy escalated from 100 to 300 W, the flavonoid content decreased from 265.5 to 97.3 TFmg QE/100 g dry weight. Nevertheless, the findings indicated that PEF treatment mitigated this decline rate. Based on the literature, PEF treatment could induce reversible breakdowns in the cell membrane at the moment of application, called electroporation. As a result, there was an increase in the permeability of the cell membrane and intracellular flavonoids were released into the surrounding medium for preservation [
47]. While PEF treatment led to enhanced preservation of vitamin C at 300 W microwave energy, no significant contrast was noted between the control and PEF-treated samples. It seems PEF decreased the generation of some reactive oxygen species, inducing oxidative stress, and thus oxidized and could preserve the vitamin C content. Ascorbic acid is easily oxidized and high oxidative stress due to PEF treatment could maintain the vitamin C content [
24,
44].
The influence of microwave energy and the application of PEF treatment on bioactive compounds in this study mirrored findings by [
29,
33,
48]. These studies revealed that PEF pretreatment could expedite the drying process without compromising the nutritional integrity of foods. Samples subjected to PEF before drying retained higher levels of polyphenols, anthocyanins and flavonoids, exhibiting improved color and flavor. Moreover, this method mitigated the risk of localized high temperatures, which could otherwise induce undesirable changes in color, flavor, nutrition and texture. However, it was observed that after intense PEF treatment, the antioxidant capacity of dried apples diminished.
3.4. Machine Learning Approach
The outcome of training using various training functions is depicted in
Figure 6a, enabling the identification of the optimal number of hidden neurons that yield the highest R
2 value and the lowest RMSE value. The statistical parameters revealed that the Levenberg–Marquardt (R
2 = 0.998, RMSE = 0.041) and RPROP (R
2 = 0.86, RMSE = 0.129) were the best and worst training functions, respectively. The optimum topology of neuron number for LM, BR, SCG, BN, GDM and RP were 20–20, 16–16, 18–18, 13–13, 15–15 and 17–17, respectively. Notably, the most favorable outcomes for each training function were achieved within the neuron range of 13–20. This suggests that augmenting the number of neurons to a certain threshold can enhance accuracy [
39].
Afterward, the performance of different transfer functions was evaluated (
Figure 6b). The
Logsig function, with the corresponding train and test of correlation coefficient and root mean squared error with the neuron topology of 18–18 and epoch of 156 had the best performance. It was followed by
tansig,
radbas,
elliotosig,
poslin and
tribas. It is noticeable that the transfer function of the first and second hidden layers was the same in
Figure 5b. Thus, to achieve improved outcomes through distinct transfer functions for the first and second layers, the network was subjected to retraining using the Levenberg–Marquardt algorithm and the top 10 models are presented in
Table 3.
For each architecture of the developed ANN model trained using the complete dataset, the cross-validation procedure was implemented across the whole data set, excluding one group from the dataset at each iteration. Statistical parameters indicated that employing diverse transfer functions in the first and second hidden layers can enhance model performance. Notably, the most proper outcomes were achieved when the network utilized
Logsig and
Tansig functions in the first and second layers, respectively. Furthermore, the findings revealed that all of the top-performing networks reached convergence within fewer than 170 training epochs and an increase in the number of neurons corresponded to an enhancement in network accuracy. The investigation of different transfer functions along with varying network structures in this study is consistent with the observations made by [
36,
39]. The findings from their studies were largely in agreement, with the primary difference between their research lying in the arrangement of hidden layers. In their work, [
36] asserted that one hidden layer was sufficient to attain high performance.
Upon establishing the optimal topology (4–18–18–1), the performance of the developed network to predict the moisture content of the apple slices was evaluated. The statistical parameters presented favorable outcomes in both the test and cross-validation sets (
Table 4). The maximum error rates for test and cross-validation were 0.038 and 0.047, respectively. Although the network’s performance was satisfactory, a decline in accuracy was observed with higher microwave energy and PEF treatment levels. This outcome suggests that alterations in temperature and cellular structure can introduce complexity into prediction conditions, thereby influencing network performance (similar to the mathematical modeling findings).
The predictive ability of artificial neural networks (ANN) in determining the moisture ratio kinetic has been corroborated in previous studies by [
26,
27]. Ref. [
26] developed an intelligent interface to assess various drying techniques for pomelo fruit (
Citrus maxima) peel, reporting a minimum correlation coefficient of 0.993. Similarly, ref. [
27] employed the ANN approach to predict moisture content in coated pineapple cubes. Their findings demonstrated a highly effective model (R
2 = 0.999, RMSE = 0.001 and MAE = 0.0007) with a topology of 3–14–14–1. However, they underscored that the appropriateness of the chosen topology depends on factors such as the type of sample, the drying methodology, the inputs, the number of samples, the number of neurons and the configuration of hidden layers employed in the analysis. Furthermore, their findings affirmed that the selection of the neuron number is contingent upon the available number of samples. In instances where there are numerous samples yet few neurons, the ANN learning process might prove inadequate, leading to an underfitted estimation. Conversely, when confronted with few samples and an abundance of neurons, the process could become time consuming, resulting in an overfitted estimation.
Also, to comprehensively assess the robustness of the ML approach, the performance of the developed SVR model was examined. Various regression parameters were employed to fine-tune different kernels and the most favorable outcomes are presented (
Table 5). Although all the optimized kernel functions exhibited satisfactory performance, the radial basis (RB) model was identified as the superior choice to predict the moisture ratio across diverse experimental conditions.
Figure 7 depicts the prediction models developed using the RBF network under various experimental conditions. The most optimal and least favorable outcomes were exhibited in the laboratory conditions of 200 W-control (R
2 = 0.994, RMSE = 0.012) and 300 W-PEF3 (R
2 = 0.847, RMSE = 0.127), respectively. Similar to the findings from mathematical and ANN modeling, the prediction accuracy was notably influenced by the quantity of microwave energy and the application of PEF treatment. The potential of SVR to predict the moisture ratio in the drying process has been confirmed by previous studies [
36,
42]. Ref. [
36] optimized the SVR model by investigating the Gaussian kernel regression parameters for the online prediction of the moisture ratio of lentil seeds in a microwave-fluidized bed dryer. Further, Ref. [
42] presented a novel modeling approach employing support vector regression methods enhanced by the dragonfly algorithm techniques to predict the drying kinetics of pea pods. They found that the optimized hyperparameters derived from the dragonfly algorithm effectively revealed the nonlinear characteristics of pea pod drying (R
2 = 0.9983 and RMSE = 0.0132).
The comparison of the employed models in this study aimed to ascertain their efficacy in predicting the drying kinetics of apple slices using the applied MV (
Table 6). The findings underscored the ANN’s robust performance as the highest among the considered models. Overall, the accuracy achieved through machine learning approaches surpassed that of the mathematical models, a notion supported by [
9,
25,
31].
Ref. [
31] employed diverse mathematical models and machine learning techniques, including ANN, KNN and SVR, to model the drying process of chanterelle mushrooms using a heat pump dryer. Due to the elevated accuracy of the ML model, they suggested its applicability in the online monitoring and control of the drying processes of chanterelle mushrooms. Additionally, Ref. [
25] developed an ANN model to predict enzyme inactivation kinetics of Irish potatoes using infrared and microwave as dry blanching tools. They reported that ML accuracy surpassed that of mathematical methods. They highlighted the superior performance of ANNs, as evidenced by an elevated coefficient of determination (R
2 = 0.963–1) compared to the best-fitting mathematical model, Weibull (R
2 = 0.862–0.969), underscoring the advantage of ANNs over mathematical modeling. Moreover, their study revealed that ANNs provide a more comprehensive comprehension of the modeling and data prediction process compared to mathematical models. This is achieved by elucidating the relationship between input and output factors through a set of rules governing data management.
Based on the results of this paper and the literature [
6,
25,
26,
27,
28,
39,
40], it can be claimed that the ANN and SVR models are prominent MLs that can be used for evidence predicting and optimizing the drying process in industrial settings. ANN models, as complex, nonlinear relationship learners, are characterized by their layered architecture approach. Their major advantage is high computational performance and adaptability to real-time applications, which are extremely necessary for under-equipped modern hardware with optimized algorithms. They can be used to process volumes of large datasets from drying operations quickly to provide correct predictions and adjustments in real time [
26,
28]. On the other side, whereas, in general, less computationally intensive and quite robust with smaller datasets, SVR models may not scale to achieve real-time performance because the computational complexity in dealing with large feature spaces and nonlinear mappings may be an issue [
40]. In this sense, in the industrial field, where dryers work against variable conditions and quick, adaptive responses to their actions are required, always an important choice will be the ANN model, due to better handling of dynamic and complex data, while the SVR models can be very useful in cases with constrained dimensions of data and simpler relationships. While both models help improve operational efficiency, ANN’s computational efficiency closely relates to it as being the instant model of choice for large industrial models.
3.5. Characterization of the Dried Samples
Variations in the color attributes of dried apples were explored across different experimental conditions. The initial color parameters for fresh fruit were measured as follows: L* = 77.25 ± 3.19, a* = 3.48 ± 0.27 and ∆E = 12.09 ± 0.89. As the microwave energy increased, the L* and ∆E values exhibited a significant decrease. For instance, in the case of the control dried apples, the L* values under 100 W, 200 W and 300 W were recorded as 75.68 ± 3.82, 69.81 ± 2.17 and 55.42 ± 2.86, respectively. Conversely, the a* values significantly decreased, indicating that the microwave’s lower average energy could inhibit the browning of apple slices. An additional factor contributing to the notable decline in the L* index during drying at 300 W power was the occurrence of darkened areas resulting from overheating. This effect can be attributed to the elevated energy levels (which result in higher temperatures) during the drying process, leading to the deactivation of enzymes such as polyphenol oxidase and peroxidase, while simultaneously augmenting the incidence of enzymatic browning reactions. Nevertheless, the vacuum part effectively eliminated the oxygen within the chamber and the inclusion of a condensation unit expedited the drying process. This combination of vacuum drying assisted by condensation is a highly effective method for preserving the color of dried apples. Ref. [
45] evaluated the color change of peaches under a combination of convective and microwave methods and they claimed that increasing the temperature considerably affected the color of the samples. Similarly Ref. [
32] used microwave drying and combined microwave-convective drying for drying blueberry fruit and assessed the effect of increasing the power on the color and some chemical aspects of blueberry fruit. They revealed that the results could be attributed to the fact that microwave-assisted hot-air-dried blueberries had no difference compared to convective drying due to pigment loss related to the higher temperature attained by samples during the process, which led to surface damage [
32]. Therefore, we can conclude that the utilization of PEF treatment can significantly improve the drying process of perishable fruit.
The PEF treatment influenced the maintenance of the color changing during the drying process; however, the value was not statistically significant. The samples subjected to PEF treatment and dried at 100 W exhibited the least noticeable color differences. This favorable impact on the sample’s color could be attributed to the alteration in cell membrane permeability caused by electroporation. This phenomenon led to reduced enzyme release and the release of substrates involved in enzymatic browning reactions as well as a decrease in pigment oxidation through thermal decomposition. Similarly, Ref. [
17] evaluated the effect of PEF on the quality parameters of osmotically dehydrated tomatoes and they revealed that PEF treatment could improve the color change of the samples during the heat treatment. On the other hand, Ref. [
29] assessed the effect of air humidity and temperature on the convective drying of apples with PEF and reported that the color of the untreated and the PEF-pretreated dried samples was similar.
The firmness and Young’s modulus values for both fresh and dried samples were assessed (
Table 7). The initial firmness and Young’s modulus of the fresh apple slices were 16.47 ± 1.57 N and 0.14 ± 0.02 MPa, respectively. The application of different levels of microwave energy and PEF treatment exerted an influence on the mechanical attributes. The findings indicated a significant decrease in the mechanical properties with escalating energy levels. Specifically, based on Young’s modulus and firmness values of the fruit dried under 300 W (F = 6.76 ± 0.71 N and Y = 0.05 ± 0.00 MPa), it became evident that the dried sample was susceptible to early breakage. This occurrence was due to the uneven dispersion of water molecules within the fruit matrix, rendering their reorganization within the structure unattainable [
30].
Although the firmness of the treated samples decreased, this decline was not statistically significant at 100 W and the firmness was able to be moderately maintained during the drying process. However, the impact of PEF treatment on the mechanical properties did not exhibit significance when subjected to drying conditions of 200 and 300 W. This mechanical response indicates an antiplasticizing effect induced by the presence of slight amounts of adsorbed water. At elevated microwave levels, a reduction in the maximum force, along with a significant refinement in the force-deformation relationship, has been observed. This corresponds to the plasticizing influence of water on the structure, resulting in an enhanced cohesion and toughness of the material. The effect of PEF on similar mechanical behavior was studied by [
14,
49]. Ref. [
14] reported that PEF improved the firmness parameter of kiwifruit when a low temperature of drying (50 and 60 °C) was used, while an increasing temperature resulted in less firmness (apart from the samples dried at 60 °C) compared to other pretreated samples. Similarly, ref. [
49] revealed that the PEF maintained the texture during the storage period. The 25, 50 and 100 kV/m of firmness loss in atemoya were 78.56%, 72.23% and 88.41%, respectively. They reported that when PEF was administered to partially dehydrated tissue, it led to increased cell disintegration.
A sensory evaluation was conducted by experts and the results have been detailed in
Table 7. The most favorable appearance quality in terms of color was observed when the sample was treated by PEF dried at 100 W. As the microwave energy increased, the visual color quality of the dried samples deteriorated due to the appearance of scorch marks caused by excessive heat. Consequently, the PEF treatment exhibited a beneficial effect on color preservation (only for the 100 W treatment). Similarly, experts expressed satisfaction with the taste quality of samples treated under 100 W, both in the controlled and treated groups, receiving a high score of 9. In contrast, the sweetness score for samples dried at 200 W and 300 W was below the acceptable threshold (i.e., 5). Finally, the total admissibility scores for dried samples subjected to 100 W, 200 W and 300 W conditions were high (9), near the threshold of acceptability (5) and under the acceptability threshold (2), respectively. Notably, the PEF treatment displayed no discernible impact on the toral acceptability score. In a similar way, ref. [
14] reported the sensory parameters of kiwifruit samples that were treated by PEF, when dried at a lower temperature (50 and 60 °C) showed an intermediate value of all parameters while, when dried at 70 °C, in addition to having the lowest score for the overall acceptability, showed the minimum texture and a balance between the sweetness and acidity level. Overall, a comparison of the current research and literature demonstrates that PEF could improve the physical, chemical and physical properties of the fruit under the heating process.
While the model developed in this study focused on pulsed electric fields and apple slices, its optimization renders it applicable to various devices and products. Initially, artificial intelligence models are employed to determine the optimal network topology. Subsequently, network inputs are tailored according to specified parameters. There is even flexibility to adjust input parameters and tailor the network based on datasets encompassing dryer specifications, PEF system characteristics and product attributes.