A Case Study of Turbulent Free Jet Flows Issuing from Rectangular Slots on Process Performances and Quality of Hot-Air-Dried Apple

: This study deals with the improvement in drying process performances and the quality of the ﬁnal product for industrial equipment in the food industry. Designers need to optimize the design parameters of devices to create synergies between the greater energy efﬁciency of the process and high-quality dried products. Air impingement drying was carried out on apple cylinders at 323 K and with air velocities ranging between 30 and 60 m s − 1 . The studied drying process presents a particular setup of jets as they are multiple rectangular slot jets issued from triangular nozzles. The effect of four design jet parameters (slot width, nozzle-to-surface height, nozzle-to-nozzle spacing, and airﬂow) on the drying process performances and the quality of the ﬁnal product was analyzed and optimized using response surface methodology (RSM). A minimal inﬂuence of design jet parameters on the process performances was shown, while an important impact was observed on the quality of dried apple. The slot width and the nozzle-to-nozzle spacing had a signiﬁcant effect on the textural and functional properties. Predictive models were established and good agreements were found between predictive and observed values. Sorption isotherms were properly modeled by the Guggenheim–Anderson–de Boer (GAB) model.


Introduction
The drying process is largely used to reduce the moisture content of agricultural products to a certain level; this is achieved by inhibiting the growth of the microorganisms and minimizing deterioration reactions to extend the shelf life of products. However, the high temperature and long drying time needed in conventional hot-air dryers negatively impact the organoleptic (color and texture), functional (rehydration), and nutritional qualities of food products [1][2][3][4].
Among the different drying methods available, the drying process that uses impinging jets seems to be promising for food thermal treatment [5]. During the air impingement drying process, the air reaches the surface of the product at high velocities which reduces the air layer that is directly in contact with the product (boundary layer), increases heat and mass transfers, diminishes the drying time, and increases the global quality of the product [6,7]. Air impingement drying has become popular over the past 25 years in the food industry and has been applied to vegetables and fruits (seedless grapes, Hami melon flakes, apricots, carrots, line peppers), herbs (American ginseng, yam), foods (tortilla chips, potato chips, fish cakes, fish particles, shrimp), seeds, and soy residue or okara [8]. Among the parameters in impingement flows, the most important parameters are the nozzle-to-product spacing, the nozzle cross-section (circular or rectangular), the nozzle exit conditions (length, smoothness, shape of the nozzle, and Reynolds number (calculated from the airflow)) [9]. Most of the quoted studies have been published by the same research group and have been carried out on the same dryer composed of a series of round nozzles in lines. The distance between the impingement nozzle and the surface of the product has been approximately 8 × 10 −2 m, and other parameters characterizing the airflow have not been reported. These studies have shown that air impingement drying reduces the drying time and increases the quality of dried products (color, texture, nutrient content, rehydration ratio, and microstructure by scanning electron microscope). The following conditions have previously been investigated: 323-353 K and 3-20 m s −1 [8]. At present, the available information concerning the parameters involved in impingement flows during the drying of food and agricultural products is limited.
The optimization of different parameters involved in impingement flows has been carried out for cooling processes. Hardisty and Can [10] studied the effects of changes in the geometry of a slot nozzle on the heat characteristics of an impinging air jet. They showed that the heat transfer coefficient is at maximum value for a ratio H/B (ratio between the nozzle-to-product height and the slot width) in the region of 6-8. Can et al. [11] experimentally determined heat transfer coefficients under arrays of impinging air jets from slot nozzles and circular holes. The results were discussed in terms of free area (ratio of the area covered by jets on the total surface area swept by jets). The authors characterized a maximum heat transfer coefficient of 254 W m −2 K −1 for a free area of 4.25%, a slot width of 2 × 10 −3 m, a height of 1.2 × 10 −2 m (ratio H/B = 6), and an air jet velocity of 40 m s −1 . More recently, Etemoglu et al. [12] conducted a numerical study to optimize the thermal processes of impinging air jets using single and multiple nozzle systems. They found that an increase in the slot width decreases the heat transfer coefficient. For example, when the slot width increases from 3 to 6 × 10 −3 m, and then to 1.2 × 10 −2 m, this results in a decrease in heat transfer coefficients from 133 to 95.6, and then to 68.6 W m −2 K −1 . The reduction factors are retained (1.4 and 2, respectively) when the air velocity increases (from 20 to 40, and then to 60 m s −1 ). For multi-nozzle systems, the heat transfer coefficient depends on the free area. A maximum coefficient of 315 W m −2 K −1 was determined for a 3% free area, 3 × 10 −3 m slot width, 1.8 × 10 −2 m height (ratio H/B = 6), and 60 m s −1 air velocity. It is worth mentioning that these studies have been carried out within the same research group.
Optimization is an essential tool in food engineering for the efficiency of processing operations and the high acceptability of the processes yield. Response surface methodology (RSM) is a statistical technique that includes optimization. This methodology has been widely used in many food processing applications, such as the optimization of osmotic dehydration, spray drying, and fluid bed drying conditions, since the early 2000s [13]. In this work, we are primarily interested in the nonlinear effect of the parameters and their interactions. For example, quadratic effects have been identified for air temperature and air velocity in the convective drying of apple slice [14].
Parthasarathi and Anandharamakrishnan [15] reviewed different modeling techniques to understand food microstructures. Among physical, chemical, and structural changes during the process, they focused on three properties that play an important role in determining food characteristics: skrinkage (the change in volume of food material due to water loss during dehydration), rehydration (the restoration of properties from the dehydrated food material), and texture (the enjoyment and acceptance of food sold to customers). Hardness is a textural property of fruits that provides important information on the storability and resistance to injury of the food products during handling and processing. The loss of texture of rehydrated apples may be highly significant, even for a gentle process such as vacuum freeze-drying [16]. Values of hardness, defined as puncture rupture force, have been used to define the textural impacts of different operations, such as air drying, texturing vacuum treatment and freezing/thawing [17]. Vega-Gálvez et al. [18] studied the effect of temperature and air velocity on the drying kinetics and quality attributes of apples (Granny Smith variety). Depending on the treatment parameters, hardness was very sensitive to drying and decreased by up to 67% with respect to the fresh sample. The sorption isotherms of fresh and processed apple (Golden Delicious variety) were measured by Sá et al. [19]. They successfully compared the Guggenheim-Anderson-de Boer (GAB) model of their data with similar curves obtained from five similar works. Moraes et al. [20] applied GAB, Brunauer-Emmet-Teller (BET), Henderson-Thompson and Oswin models to moisture sorption isotherms of apple at 323 K. They mentioned that the parameters of the GAB and BET models have a physical meaning, and identified the best fitting for the GAB model. These results confirmed that the GAB model presented the best fit for more than 50% of the fruits, meats, and vegetables, as stated by Lomauro et al. [21].
This study has two major objectives. The first is to determine the influence of four design jet parameters, namely the slot width, the nozzle-to-surface height, the nozzle-tonozzle spacing, and the airflow, on the drying process performances and the quality of the final products. The second is to propose empirical predictive models according to the most relevant design jet parameters. Cylinder-shaped apples were chosen as a model food.

Materials and Methods
This research was carried out on apples because they are a product largely studied in the literature, and are easily available and cheap. Moreover, the apple is one of the most popular and widespread fruits around the world, with a production of 87 million tons per year [22]. The Gala variety was selected as it is one of the most appreciated varieties worldwide.

Raw Material and Sample Preparation
Apples (Malus domestica Gala) were harvested in September/October and kept in storage rooms at 277-279 K for six months. They were issued from a local production near to the supermarket (6 km distance). Before they reached the store, the apples were calibrated, and the overall quality of the product was regularly controlled. The apples were purchased at a supermarket and stored in our lab at 277 K until use. Cylindrical samples were cut through the central axis with a cutting shape (2 × 10 −2 m in diameter and 4 × 10 −2 m in height). Cylinder-shaped samples were chosen as the geometry most extensively investigated is the slab. It was specified that samples were placed vertically. To maintain the samples in a vertical position during the drying, polytetrafluoroethylene (PTFE) sample holders were designed (Figure 1a). The initial mass of the samples was 13.4 ± 0.3 × 10 −3 kg.

Drying Processes
Two types of drying, air impingement and oven drying, at 323 K were compared, particularly for the sorption capacity of the dried apples, as oven drying at a moderate temperature and with natural convection is a slow drying process. In such conditions, tissue shrinkage and collapse should be limited, and the formation of the crust layer on the surface of the dried materials should be avoided due to low thermal gradients.

Air Impingement Drying
Boy et al. [23] and Boy et al. [24] previously described the air impingement dryer equipment. The dryer mostly consists of a dehumidifier and a centrifugal fan to supply the air in two separate cylindrical ducts (0.25 m diameter), each with an electrical heater and joined to each other in the drying chamber. At the exit of the drying chamber, the air circulates back to the fan with a degree of recirculation of 100% (Figure 1b). Inside the drying chamber, the air temperature was controlled. Both the hot air heating and the airflow were in continuous operation, meaning that the airflow continuously circulated in a closed circuit surrounding the samples.
The setup of the air jets was original as they were created by triangular nozzles pierced with rectangular slots (Figure 1c). Two nozzle-to-nozzle spacings (∆ = 5 × 10 −2 m and 10 × 10 2 m), two slot widths (B = 2 × 10 −3 m and 4 × 10 −3 m), and two heights between the nozzle outlet and the sample surface for the upper part were tested (H upper = 1.2 × 10 −2 m and 2.4 × 10 −2 m). The height for the lower part was fixed (H lower = 2.5 × 10 −2 m). The airflow was circulated perpendicular to the sample surface and the sample was placed Processes 2021, 9,1900 4 of 21 on a grid subjected to a back-and-forth motion which allowed for the entire surface to be covered by the airflow. storage rooms at 277-279 K for six months. They were issued from a local production near to the supermarket (6 km distance). Before they reached the store, the apples were calibrated, and the overall quality of the product was regularly controlled. The apples were purchased at a supermarket and stored in our lab at 277 K until use. Cylindrical samples were cut through the central axis with a cutting shape (2 × 10 -2 m in diameter and 4 × 10 −2 m in height). Cylinder-shaped samples were chosen as the geometry most extensively investigated is the slab. It was specified that samples were placed vertically. To maintain the samples in a vertical position during the drying, polytetrafluoroethylene (PTFE) sample holders were designed (Figure 1a). The initial mass of the samples was 13.4 ± 0.3 × 10 −3 kg.

Drying Processes
Two types of drying, air impingement and oven drying, at 323 K were compared, particularly for the sorption capacity of the dried apples, as oven drying at a moderate temperature and with natural convection is a slow drying process. In such conditions, tissue shrinkage and collapse should be limited, and the formation of the crust layer on the surface of the dried materials should be avoided due to low thermal gradients.

Air Impingement Drying
Boy et al. [23] and Boy et al. [24] previously described the air impingement dryer equipment. The dryer mostly consists of a dehumidifier and a centrifugal fan to supply the air in two separate cylindrical ducts (0.25 m diameter), each with an electrical heater and joined to each other in the drying chamber. At the exit of the drying chamber, the air circulates back to the fan with a degree of recirculation of 100% (Figure 1b). Inside the drying chamber, the air temperature was controlled. Both the hot air heating and the airflow were in continuous operation, meaning that the airflow continuously circulated in a closed circuit surrounding the samples.
The setup of the air jets was original as they were created by triangular nozzles pierced with rectangular slots (Figure 1c). Two nozzle-to-nozzle spacings (Δ = 5 × 10 −2 m and 10 × 10 2 m), two slot widths (B = 2 × 10 −3 m and 4 × 10 −3 m), and two heights between the nozzle outlet and the sample surface for the upper part were tested (Hupper = 1.2 × 10 −2 m and 2.4 × 10 −2 m). The height for the lower part was fixed (Hlower = 2.5 × 10 −2 m). The airflow was circulated perpendicular to the sample surface and the sample was placed on a grid subjected to a back-and-forth motion which allowed for the entire surface to be covered by the airflow.
Experiments were performed at a moderate air temperature of 323 K for two airflows (3460 m 3 h −1 and 4282 m 3 h −1 ), which corresponded to 80% and 99% of the nominal frequency of the fan motor (T2A90S-2, TECHTOP © Motor, Shanghai, China). The corresponding air velocities at the slot outlets were measured using an anemometer with a pitot tube (TA400, range 1-80 m s −1 , accuracy ± 2.5% at 10 m s −1 , Trotec, Heinsberg, Germany). Air velocities ranged from 30 to 60 m s −1 . During the experiments, a T-type thermocouple (accuracy ± 0.1 K) was used to measure the air temperature under the nozzles, and the data were collected every 10 s by an acquisition unit (SR Mini System, TC Experiments were performed at a moderate air temperature of 323 K for two airflows (3460 m 3 h −1 and 4282 m 3 h −1 ), which corresponded to 80% and 99% of the nominal frequency of the fan motor (T2A90S-2, TECHTOP © Motor, Shanghai, China). The corresponding air velocities at the slot outlets were measured using an anemometer with a pitot tube (TA400, range 1-80 m s −1 , accuracy ± 2.5% at 10 m s −1 , Trotec, Heinsberg, Germany). Air velocities ranged from 30 to 60 m s −1 . During the experiments, a T-type thermocouple (accuracy ± 0.1 K) was used to measure the air temperature under the nozzles, and the data were collected every 10 s by an acquisition unit (SR Mini System, TC S.A., Dardilly, France). The air humidity under the nozzles was measured using a hygrometer (HygroPalm 22-A, 0-100%, accuracy ± 1.2%, Rotronic, Saint-Priest, France). The mean values of air temperature and humidity were as follows: 323.0 ± 2.4 K and 11.8 ± 1.6%.
Four samples were dried at the same time and the weightings were made using a precision scale (Acculab, ACL-110.4, accuracy 1 × 10 −7 kg, Sartorius, Göttingen, Germany); this was firstly completed at 10 min intervals during the first two hours, then the time interval was changed to 30 min for the next four hours, to finally be expanded to 1 h after 6 h of drying. All experiments lasted from 7 to 10 h. The experiments were continued until a moisture ratio ω of 0.1 was obtained. The moisture ratio (ω) was defined according to the following expression: where M t is the moisture content at time (t) (kg kg −1 , dry basis (d.b.)), M 0 is the initial moisture content (kg kg −1 , d.b.), and M e is the equilibrium moisture content (kg kg −1 , d.b.) M t and M 0 were determined by the ratio of the mass of water contained in the Processes 2021, 9, 1900 5 of 21 sample at time (t) and at t = 0, respectively, on the dry solids content (see Section 2.2.3 for further details). Dried samples were cooled in ambient conditions, packed into heat-sealed polyethylene bags, kept at 277 K, and protected from light until further analysis.

Oven Drying
Triplicate samples were dried in an oven (BM 600, Memmert, Schwabach, Germany) at 323 ± 2 K. Samples were withdrawn and weighed every 30 min until they reached a moisture ratio of 0.1. After the experiment, samples were preserved in the same manner as those from the air impingement drying.

Equilibrium Moisture Content and Dry Solids Content
The equilibrium moisture content was determined by placing dried samples in a drying cabinet (Nalgene 5317-0120, Sigma-Aldrich, Saint-Louis, MO, USA), all samples were left in an oven at 323 K to a constant weight (deviation less than 10 −6 kg). Then, the samples were transferred to an oven at 378 ± 2 K for at least 24 h to obtain the dry solids content.

Sorption Isotherms
Before performing the sorption isotherms, dried samples were stabilized by a gravimetric method [20]. Samples were transferred to a drying cabinet (Nalgene ® , Thermo Fisher Scientific, Waltham, MA, USA) and placed in an oven in which the relative humidity was controlled by a saturated salt solution. To maintain the samples at a constant water activity a w of 0.8 (corresponding to ω = 0.1), a saturated salt solution of (NH 4 ) 2 SO 4 was used at 323 K. Samples were then weighed until equilibrium with the surrounding air was reached (constant weight). Once stabilized, sorption isotherms were measured using a Dynamic Vapor Sorption (DVS) equipment (IGAsorp HT, Hiden Instruments Ltd., Warrington, UK). About 2.1 × 10 −5 kg of dried apples were placed under an air flow (4.2 × 10 −6 m 3 s −1 ) until a constant weight was achieved. The mass uptake and sorption kinetics were recorded through a fine balance (10 −7 kg resolution) and a control of both temperature and humidity. Curves were registered for increasing humidity from 80 to 96%, in 2 increments, and the temperature was kept constant at 323 K. The equilibrium mass was determined by the exponential increase in mass over time for each increment. Samples were then exposed to zero humidity to estimate the dry solids content. Sorption isotherms were determined through triplicate samples.
Experimental data were fitted using the Guggenheim-Anderson-de Boer (GAB) model because this model has been found to be the most fitting model to the sorption isotherm of apple (Fuji variety) at 323 K [20,25]. The equation was defined as: where M m is the monolayer moisture content (kg kg −1 , d.b.), C g and K are constants, and a w is the water activity.

Heat and Mass Transfer Parameters
The determination of heat and mass transfer coefficients h and k m , respectively, requires the knowledge of some physical properties which were estimated from the literature [26,27].

Heat Transfer Coefficient
The heat transfer coefficient was determined from the thermal history of an aluminum plate, which is known as the lumped capacitance method [28,29]. This method assumes that the internal temperature gradient is negligible because the aluminum thermal conductivity is high. Thus, the resistance to the heat transfer is located on the surface. The plate was Processes 2021, 9,1900 6 of 21 placed in the refrigerator to maintain a maximum temperature difference with the hot airflow. The increase in temperature of the aluminum plate was measured until it reached air temperature (323 K). The heat transfer coefficient (h) (W m −2 K −1 ) was calculated by solving the heat balance that occurred when the plate was surrounded by air: where ρ Al is the density of aluminum (kg m −3 ), C pAl is the specific heat of aluminum (J kg −1 K −1 ), V is the volume of the aluminum plate (m 3 ), T is the temperature (K), t is the time (s), A is the surface area of the aluminum plate (m 2 ), and T a is the hot air temperature (K). After integration, Equation (3) becomes: where T 0 is the initial temperature of the aluminum plate (K). The value of h was deduced from the slope of the line in the graph representing ln T* vs. time.

Mass Transfer Parameters
The mass transfer coefficient (k m ) can be estimated by an analogy with the heat transfer coefficient according to the correlation developed by Chilton and Colburn [30]: where ρ air is the density of air (kg m −3 ), C p air is the specific heat of air (J kg −1 K −1 ), Pr is the Prandtl number, and Sc is the Schmidt number. The effective moisture diffusivity (D eff ) (m 2 s −1 ) was calculated from Fick's second law of diffusion (Equation (6)) when the moisture diffusion is the mass transport mechanism that limits the drying process: where M db is the moisture content (kg kg −1 , d.b.) and x is the characteristic dimension of the sample (m). The mathematical solution of Equation (6), established by Crank [31] for an infinite slab, was expressed by: where i is the number of terms considered in the series of the exponential functions and L is the sample half thickness (m). It should be noted that Equation (7) involves numerous assumptions, including a one-dimensional transport of moisture within the product, a negligible shrinkage, and a constant diffusion coefficient (for further details, see Boy et al. [23]). Among the different geometries, we chose the solution developed for an infinite slab. This assumes an isotropic material, which is the case for apple. However, the solution developed for an infinite cylinder implies that the diffusion is radial. This assumption does not hold in our case. To take into account the sample's shape, we introduced the equivalent thickness (Le) in Equation (7) instead of the sample thickness (2 L). Le was calculated by dividing the volume of the cylinder by the surface area (Le = 4 × 10 −3 m). As stated previously by Boy et al. [24], two exponential terms (i = 1) were sufficient to identify the effective moisture diffusivity (D eff ) of apple. Accordingly, Equation (7) was simplified to the first two terms.
Although the velocities studied were high (between 30 and 60 m s −1 ), it is more appropriate to call the diffusivity characterized in Equation (7) as apparent moisture diffusivity (D app ), as it may be a function of the air velocity. Indeed, the diffusivity introduced in Equation (7) should be an intrinsic material property and independent of the external drying conditions. The evolution of D app versus the air velocity was not studied as it goes beyond the scope of the present work.

Goodness of Fit
Non-linear regressions were used to identify the parameters of the selected models, namely, M m , C g , K, and D eff (Matlab, The Mathworks, Inc., Natick, MA, USA). The coefficient of determination (R 2 ), the adjusted coefficient of determination (adjusted R 2 ), the root-mean-square error (RMSE), the chi-square (χ 2 ), and the residual sum of squares (RSS) were determined to evaluate the fitting of Equations (2) and (7) to the experimental data. These parameters were calculated as mentioned below: where ω exp,j and ω pre,j are the experimental and predicted moisture ratios, respectively, ω exp,j is the average of the experimental moisture ratio, N is the number of observations, and p is the number of the parameters involved in the selected models (three parameters in Equation (2) and a single parameter in Equation (7), respectively).

Texture Analysis
The penetration force required to puncture the apple tissue was measured using a Texture Analyzer (TA HD Pluc C, Stable Micro Systems Ldt, Godalming, UK). A probe (2 × 10 −3 m diameter convex tip) was used to puncture the sample through the central axis at a constant speed of 1.5 × 10 −3 m s −1 over a distance of 3 × 10 −2 m. The force versus distance curves were plotted and two parameters were analyzed ( Figure 2): hardness (F) (calculated in N), and the work associated with hardness (W) (calculated in N m). Hardness (F) is defined as the maximum force on the curve and represents the force required for the rupture of the flesh. The work associated with hardness (W) is calculated as the area under the curve up to F, and reflects the mechanical work needed for the rupture of the flesh. Both definitions and methods of calculation were based on earlier works [32]. Ten and five measurements were made in duplicate for fresh and dried apples, respectively. curves were plotted and two parameters were analyzed ( Figure 2): hardness (F) (calculated in N), and the work associated with hardness (W) (calculated in N m). Hardness (F) is defined as the maximum force on the curve and represents the force required for the rupture of the flesh. The work associated with hardness (W) is calculated as the area under the curve up to F, and reflects the mechanical work needed for the rupture of the flesh. Both definitions and methods of calculation were based on earlier works [32]. Ten and five measurements were made in duplicate for fresh and dried apples, respectively.
(a) (b) Figure 2. Force versus distance curve obtained during a penetration test on: (a) fresh apple; (b) dried apple by air impinging jets (trial n°13, Table 2).

Rehydration Ratio
The rehydration ratio (RR) was carried out to measure the capacity of dried samples to rehydrate. First, the rehydration time was assessed. For this purpose, the samples were immersed in 5 × 10 −5 m 3 distilled water at ambient temperature and were weighed regularly with a precision scale (ENTRIS 124i-1, accuracy 10 -7 kg, Sartorius, Göttingen, Germany) until the mass was stable. Before weighing, the samples were drained on absorbent paper to remove the surface water [33]. The time needed for the rehydration process was 15,800 s (equal to 4 h and 23 min). Then, the rehydration ratio was calculated by weighing the samples before and after soaking in distilled water for 15,800 s (4 h and 38 min) [34]: where mτ = 15,800 s is the mass after rehydration (kg) and m0 is the initial mass before drying (kg). Rehydration measurements were run in duplicate.  Table 2).

Rehydration Ratio
The rehydration ratio (RR) was carried out to measure the capacity of dried samples to rehydrate. First, the rehydration time was assessed. For this purpose, the samples were immersed in 5 × 10 −5 m 3 distilled water at ambient temperature and were weighed regularly with a precision scale (ENTRIS 124i-1, accuracy 10 −7 kg, Sartorius, Göttingen, Germany) until the mass was stable. Before weighing, the samples were drained on absorbent paper to remove the surface water [33]. The time needed for the rehydration process was 15,800 s (equal to 4 h and 23 min). Then, the rehydration ratio was calculated by weighing the samples before and after soaking in distilled water for 15,800 s (4 h and 38 min) [34]: RR = m τ =15,800 s m 0 (13) where m τ = 15,800 s is the mass after rehydration (kg) and m 0 is the initial mass before drying (kg). Rehydration measurements were run in duplicate.

Design of Experiments
The effect of design jet parameters on both the drying process performances and the quality attributes of apples were studied using the response surface methodology (RSM). The design jet parameters (independent variables) were slot width (B), nozzle-to-surface height (H), nozzle-to-nozzle spacing (∆), and airflow (υ). Two levels were selected, and the variation ranges are presented in Table 1.
Among the response parameters (dependent variables), the drying process performances were the drying time to reach a moisture ratio ω of 0.5 t, the heat and mass transfer coefficients (h and k m ), hardness (F), the work associated with hardness (W), and the rehydration ratio (RR). A full factorial design of experiments (DoE) was completed and included 16 trials (2 k = 2 4 trials), as shown in Table 2. Each trial was completed in duplicate.  The response variables were analyzed using Statgraphics Plus 5.1 software (Statgraphics Technologies Inc., Warrenton, VA, USA). A statistical model (Equation (14)) was established to predict the response as a function of the studied parameters: where Y is the response variable, β are the regression coefficients, x i is the independent variable, and k is the number of studied parameters. The statistical model was developed by deleting the least significant parameters. The analysis of variance (ANOVA) was carried out to evaluate the significant differences between the studied parameters for each response at p-value ≤ 0.1 (Table 3).  The lower the p-value, the more significant the parameter. The relevance of the developed empirical model was estimated using the coefficient of determination (R 2 ) and adjusted R 2 , the mean absolute error, and the standard error.
The optimization of the jet impingement process was achieved to find the best possible combination of the design jet parameters to ensure, if possible, the highest process performances and the best quality of the dried apples. The optimum parameters were identified to minimize the drying time to reach a moisture ratio ω of 0.5, hardness, and the work associated with hardness, and to maximize the heat and mass transfer coefficients and the rehydration ratio. Low textural properties (hardness and the work associated with hardness) were selected to prevent damage to the product, such as the case hardening phenomenon.

Results and Discussion
The initial moisture content of apple was 87.7 ± 0.7% in wet basis (w.b.), which corresponds to 7.2 ± 0.4 kg kg −1 in d.b.
The statistical study will be presented by distinguishing, on the one hand, the effect of the four operating parameters on the drying process performances (i.e., drying time, heat and mass transfer coefficients) and, on the other hand, the quality of the final product (i.e., apparent diffusivity and both textural and functional properties).

Effect of Design Jet Parameters on Drying Time
The drying time to reach a moisture ratio of 0.5 (corresponding to 78.6 ± 1.2% (w.b.) and 3.7 ± 0.3 kg kg −1 (d.b.)) ranges from 4860 s to 6060 s (1 h and 35 min and 2 h and 10 min, respectively) for all the configurations. From Table 3, both interactions (slot width × nozzleto-surface height) and (nozzle-to-nozzle spacing × airflow) have a significant positive effect (p-value ≤ 0.1) on the drying time. In contrast, the main factors (the slot width, the nozzleto-surface height, the nozzle-to-nozzle spacing, and the airflow) have no significant effect (p-value ≥ 0.1) on the drying time.
The RSM study established a statistical model, Equation (15) in Table 4, by adjusting the drying time to the studied jet setup parameters. Table 4. Statistical models expressing the drying process performances and the final quality of the apple as a function of the studied parameters and their interactions. The coefficient of determination (R 2 ) and the adjusted R 2 calculated by the model are 52.1% and 20.2%, respectively, which indicates a large disparity between experimental and modeled values. Therefore, the drying time cannot be accurately calculated using the parameters examined.
The moisture ratio versus time is plotted on Figure 3 for two extreme drying times (trials n • 8 and n • 9, Table 2). The coefficient of determination (R 2 ) and the adjusted R 2 calculated by the model are 52.1% and 20.2%, respectively, which indicates a large disparity between experimental and modeled values. Therefore, the drying time cannot be accurately calculated using the parameters examined.
The moisture ratio versus time is plotted on Figure 3 for two extreme drying times (trials n°8 and n°9, Table 2). The moisture ratio decreases continually over time. Both curves show a rapid moisture elimination from the product at the initial stage of the process (up to 2000 s), to slow down with time. Initially, the surface moisture evaporates very rapidly due to high heat and mass transfer coefficients. The drying rate is exceedingly high during the first moments, and then declines exponentially when all the surface moisture is evaporated and the evaporation front is moved within the product. The curves have the same shape, regardless of the configuration of jets. The time to reach a moisture ratio of 0.1 depends on the configurations studied; it is 25,260 s for trial n°8, and 28,860 s for trial n°9 ( Figure  3).
The drying time obtained at 50 °C for two types of drying (air impingement and oven drying) can be compared. The time required to reduce the moisture ratio to 0.5 is 23,460 s (6 h and 52 min, data not presented). This means that the drying time determined using impinging jets is around 75 to 80% lower than the corresponding oven drying time. The drying time can be discussed in relation to other results from the literature for convective drying. Contreras et al. [35] carried out research on apples (Granny Smith variety) in a cylindrical shape (2 × 10 −2 m in diameter and 7 × 10 −3 m in height) and subjected to air drying (323 K and 2.5 m s −1 ). Samples were positioned perpendicular to the airflow. The results revealed that 5,400 s (1 h and 30 min) is required to reach a moisture ratio of 0.5.

Effect of Design Jet Parameters on Heat and Mass Transfer Coefficients
Heat and mass transfer coefficients under the different conditions were in the range of 157-273 W m −2 K −1 and 1.46 10 −4 -2.54 10 −4 m s −1 , respectively ( Table 2). The nozzle-to-nozzle spacing and the interaction (nozzle-to-surface height × nozzle-to-nozzle spacing) have a significant effect (p-value ≤ 0.1) on the transfer coefficients (Table 3). In The moisture ratio decreases continually over time. Both curves show a rapid moisture elimination from the product at the initial stage of the process (up to 2000 s), to slow down with time. Initially, the surface moisture evaporates very rapidly due to high heat and mass transfer coefficients. The drying rate is exceedingly high during the first moments, and then declines exponentially when all the surface moisture is evaporated and the evaporation front is moved within the product. The curves have the same shape, regardless of the configuration of jets. The time to reach a moisture ratio of 0.1 depends on the configurations studied; it is 25,260 s for trial n • 8, and 28,860 s for trial n • 9 ( Figure 3).
The drying time obtained at 50 • C for two types of drying (air impingement and oven drying) can be compared. The time required to reduce the moisture ratio to 0.5 is 23,460 s (6 h and 52 min, data not presented). This means that the drying time determined using impinging jets is around 75 to 80% lower than the corresponding oven drying time. The drying time can be discussed in relation to other results from the literature for convective drying. Contreras et al. [35] carried out research on apples (Granny Smith variety) in a cylindrical shape (2 × 10 −2 m in diameter and 7 × 10 −3 m in height) and subjected to air drying (323 K and 2.5 m s −1 ). Samples were positioned perpendicular to the airflow. The results revealed that 5400 s (1 h and 30 min) is required to reach a moisture ratio of 0.5.

Effect of Design Jet Parameters on Heat and Mass Transfer Coefficients
Heat and mass transfer coefficients under the different conditions were in the range of 157-273 W m −2 K −1 and 1.46 10 −4 -2.54 10 −4 m s −1 , respectively ( Table 2). The nozzle-tonozzle spacing and the interaction (nozzle-to-surface height × nozzle-to-nozzle spacing) have a significant effect (p-value ≤ 0.1) on the transfer coefficients (Table 3). In either case, the regression coefficient is negative for the nozzle-to-nozzle spacing, which suggests a facilitated transfer when the nozzle-to-nozzle spacing is low. This may be explained by the fact that the lower the nozzle-to-nozzle, the more often the product is in contact with the air, and the more the heat and mass transfers increase.
Equations (16) and (17) (Table 4) were used to calculate heat and mass transfer coefficients for the levels of parameters studied, respectively. As for the drying time, R 2 of 54.1% and 54.4%, and the adjusted R 2 of 31.2% and 31.5%, were found. These values are insufficient to properly predict heat and mass transfer coefficients. The statistical analysis shows that all parameters studied do not have a significant impact on both heat and mass transfer coefficients.
Etemoglu et al. [12] reported that the heat transfer coefficient rises to attain a maximum value for a given free area (ratio of the area covered by slots on the total surface area swept by jets). The coefficient h reaches 315 W m −2 K −1 (maximal value) at a free area of 3%, then decreases to 270 W m −2 K −1 at 6.7% of opening (multiple nozzles, B = 3 × 10 −3 m and 60 m s −1 ). The highest coefficient calculated in this investigation was 273 W m −2 K −1 at 4% of the free area (trial n • 9: B = 2 × 10 −3 m and 35 m s −1 , Table 2). For information purposes, the free area in this study ranged from 1 and 6% for all configurations. Can et al. [11] revealed the dependence of the free area on the heat transfer coefficient on multiple nozzle systems. For free area values less than 5% and H/B values (ratio between the nozzle-tosurface height and the slot width) below 6, the heat transfer coefficient under each jet is independent, i.e., not influenced by the presence of the surrounding jets, and the coefficient (h) reaches a maximum. In our study, these two conditions are satisfied for eight trials (trials n • 1, 2, 5, 6, 7, 8, 9, and 10; Table 2). In these conditions, the heat transfer coefficient increases linearly with the free area and the H/B ratio. Moreover, Can et al. [11] observed that the heat transfer coefficient goes up when decreasing the slot width; it increases from 241 W m −2 K −1 for B = 2.5 × 10 −3 m to 254 W m −2 K −1 for B = 2 × 10 −3 m (an increase of 5%). Only half of our results (trials n • 1-8, 3-6, 9-16, and 12-13) are in agreement with this observation. This difference is most likely due to the uncertainty in the determination of the heat transfer coefficient. Zhu et al. [36] studied both effects of nozzle-to-surface height and nozzle-to-nozzle spacing on the heat transfer coefficient. The authors established that, for 5 × 10 −2 m and 10 × 10 −2 m of height, the heat transfer coefficient is at its maximum value. When the height is fixed, the coefficient (h) is reduced as the nozzle-to-nozzle spacing increases. For instance, h diminishes from 314 W m −2 K −1 to 243 W m −2 K −1 and then 207 W m −2 K −1 for a nozzle spacing of 3 × 10 −3 m, 5 × 10 −3 m, and 7 × 10 −3 m, respectively. This result confirms the findings of our study since we have noticed that a lower nozzle spacing improves transfer coefficients.
Despite the variability in results due to difficulties in designing and measuring near jets, our results are perfectly in line with previous works.
Finally, the statistical study carried out to optimize the performances of industrial dryers makes it possible to conclude that all tested factors have a minimal impact on the analyzed responses. The reason that may explain these results is that the levels of the selected factors are too close to the optimum operating conditions [11,12]. The levels of each factor have been selected according to previous studies, in particular those carried out by Can et al. [11] and Etemoglu et al. [12].

Effect of Design Jet Parameters on the Apparent Moisture Diffusivity
The influence of the design jet parameters on the apparent moisture diffusivity of apple was investigated (Figure 4).
The average value obtained in the 16 trials of the DoE was 1.39 ± 0.06 × 10 −10 m 2 s −1 . A one-way ANOVA on D app values shows few significant differences between the tests (p-value ≤ 0.1). From Figure 4, the highest value of D app (1.56 ± 0.19 × 10 −10 m 2 s −1 ) is obtained in trial n • 10 (B = 2 × 10 −3 m, H = 1.2 × 10 −2 m, ∆ = 5 × 10 −2 m, and υ = 99%), while the lowest value (1.13 ± 0.01 × 10 −10 m 2 s −1 ) is achieved in trial n • 2 (B = 2 × 10 −3 m, H = 1.2 × 10 −2 m, ∆ = 10 × 10 −2 m, and υ = 99%). The difference between the two tests comes from the nozzle-to-nozzle spacing value. This is an interesting result because it states that a low nozzle-to-nozzle spacing helps to increase the moisture diffusivity within the product during the drying. Furthermore, this result supports the ANOVA analysis on heat and mass transfer coefficients (previous section) where it has been noticed that a low nozzle-to-nozzle spacing tends to improve the transfer phenomena that occur during the drying process.

Effect of Design Jet Parameters on the Apparent Moisture Diffusivity
The influence of the design jet parameters on the apparent moisture diffusivity of apple was investigated (Figure 3). The difference between the two tests comes from the nozzle-to-nozzle spacing value. This is an interesting result because it states that a low nozzle-to-nozzle spacing helps to increase the moisture diffusivity within the product during the drying. Furthermore, this result supports the ANOVA analysis on heat and mass transfer coefficients (previous section) where it has been noticed that a low nozzleto-nozzle spacing tends to improve the transfer phenomena that occur during the drying process.
The literature presents many values of Deff for different cultivars of apple and numerous scenarios of drying. Kaya et al. [37] reported effective moisture diffusivities within a range from 0.48 × 10 −10 m 2 s −1 to 2.02 × 10 −10 m 2 s −1 during the air drying of Red Delicious apples (6 × 10 −3 m thick) at 308 K, 318 K, and 328 K, velocities at 0.2 m s −1 , 0.4 m s -1 , and 0.6 m s −1 , and humidities of 40%, 55%, and 70%. In the same way, Doymaz [38] obtained diffusivities between 3.00 × 10 −10 m 2 s −1 , and 4.13 × 10 −10 m 2 s −1 for Amasya red apples (5 × 10 −3 m thick) dried at 328 K, 338 K, and 348 K at 2 m s −1 . More recently, ElGamal et al. [39] calculated effective diffusivities of Red Delicious apples (6 × 10 −3 m thick) dried in a solar dryer equipped with light sensors to track the position of the sun at a different time of day. These results have been compared with a sample batch dried in a fixed solar dryer and another dried under the temperature and velocity of ambient air. Diffusivities have ranged from 2.42 × 10 −10 m 2 s −1 and 5.43 × 10 −10 m 2 s −1 . Not surprisingly, the highest value has been obtained in the adjustable solar dryer at the highest airflow rate (44 m 3 h −1 ), while the lowest value has been recorded for apples subjected to ambient drying. Overall, reported values of Deff in our study are consistent with the literature's data. Not surprisingly, the highest value has been obtained in the adjustable solar dryer at the highest airflow rate (44 m 3 h −1 ), while the lowest value has been recorded for apples subjected to ambient drying. Overall, reported values of D eff in our study are consistent with the literature's data.

Effect of Design Jet Parameters on Textural Properties
The hardness of dried apple varies from 1.52 ± 0.33 N to 2.67 ± 0.48 N, and the work associated with hardness ranges between 21.88 ± 9.18 10 −3 N m and 38.39 ± 7.86 10 −3 N m ( Table 2). The results of the ANOVA show that the slot width has a very significant effect (p-value ≤ 0.01) on hardness ( Table 3). The nozzle-to-nozzle spacing and the interaction (nozzle-to-nozzle spacing × airflow) have a significant effect (p-value ≤ 0.1) on the work associated with hardness. In addition, the slot width has a very significant effect (p-value ≤ 0.01) on the work associated with hardness (Table 3). Moreover, the latter increases with the rise of the nozzle-to-nozzle spacing, while it decreases when the slot width diminishes as illustrated in Figure 5a,b.
(p-value ≤ 0.01) on hardness ( Table 3). The nozzle-to-nozzle spacing and the interaction (nozzle-to-nozzle spacing × airflow) have a significant effect (p-value ≤ 0.1) on the work associated with hardness. In addition, the slot width has a very significant effect (p-value ≤ 0.01) on the work associated with hardness (Table 3). Moreover, the latter increases with the rise of the nozzle-to-nozzle spacing, while it decreases when the slot width diminishes as illustrated in Figure 5a,b. The mathematical models that allow adjusting the hardness and the work associated with hardness according to the studied parameters through the RSM study are presented The mathematical models that allow adjusting the hardness and the work associated with hardness according to the studied parameters through the RSM study are presented by Equations (18) and (19) (Table 4), respectively. To realize the fitting level of the models, predicted versus observed values are plotted for the work associated with hardness. The adjustment of the model is reviewed by a good distribution of the values around the bisector (Figure 5c). The coefficients of determination (R 2 ) of 60.5% and 69.3% and adjusted R 2 of 46.2% and 54.0% indicate the discrepancies between the response values and models for hardness and the work associated with hardness, respectively (Table 3). Consequently, the two models (Equations (18) and (19)) are satisfying to predict, in a relevant manner, hardness and the work associated with hardness of the dried apple in the considered conditions. The optimal conditions, identical for the two responses, are cited in Table 3 (B = 4 × 10 −3 m, H = 2.4 × 10 −2 m, ∆ = 5 × 10 −2 m, and υ = 80%).
The hardness of fresh apple is 5.74 ± 1.17 N, which is 2 to 4 times higher than the hardness measured during the tests (Figure 1 and Table 2). A similar value has been found in the case of Granny Smith apples [18]. A significant reduction in hardness is observed in the course of the drying process; this result is supported by earlier studies [40]. Textural properties of dried products rely on the behavior of the cellular matrix, the soluble solid phase within the tissues, all combined with their interaction with water. A decline in the moisture content of the product caused by the drying process affects the structure and the textural properties. Changes in structure and textural properties of dried apple may be from the amount and composition of pectin in apple tissues [41].

Effect of Design Jet Parameters on Functional Properties
Rehydration consists of restoring the properties of the raw material to a dry product by bringing it into contact with water. However, drying and rehydration processes induce changes in the structure and composition of the product, leading to damaged reconstitution properties [42]. Therefore, rehydration is a quality criterion that informs on chemical and structural injuries that occur during the drying process.
Values of the rehydration ratio of apple samples dried in different operating conditions are between 39.68 ± 5.07% and 49.32 ± 1.95% (trials n • 5 and 15, Table 2). The interaction (slot width × airflow) has a significant effect (p-value ≤ 0.1), whereas the interaction (slot width × nozzle-to-nozzle spacing) and the nozzle-to-nozzle spacing have a very significant effect (p-value ≤ 0.01) on the rehydration ratio (Table 3). Besides, the regression coefficient is negative for the nozzle-to-nozzle spacing, which means that improved rehydration is obtained when the nozzle spacing is low. On another note, the effect of nozzle-to-nozzle spacing on the rehydration ratio is shown in Figure 5a. The ratio declines when nozzle-tonozzle spacing increases.
The statistical model established by Equation (20) (Table 4) can be used to predict the rehydration ratio, even if all the parameters are not significant. The coefficients of determination R 2 of 84.5% and the adjusted R 2 of 74.2% reveal the quite acceptable quality of the model to predict the rehydration ratio of dried apple. Figure 5d depicts the distribution of predicted versus observed values around the bisector. The developed model may be used to properly predict the rehydration ratio of dried apples. As indicated earlier, the optimal conditions to get the best rehydration are reported in Table 3 (B = 4 × 10 −3 m, H = 1.2 × 10 −2 m, ∆ = 5 × 10 −2 m, and υ = 99%).
Given that most of the dried products need to be rebuilt before use, it is fundamental to know their rehydration behavior. The highest values of rehydration ratio reveal stiffer and more open structures [34]. Following Vadivambal and Jayas [43], the highest values are characteristic of high quality dried products. It is widely accepted that tissue rehydration consists of three simultaneous processes: the imbibition of water within the dried product, the swelling of the product, and the leaching of soluble solids derived from the dried product [44]. Lewicki and Wiczkowska [45] compared the features of the rehydration of apple (Idared variety) dried by convection (343 K and 1.5 m s −1 ) and freeze-dried (frozen at 252 K for 20 h and then freeze-dried at 37 Pa and 293 K for 24 h). They showed that the collapsed structure of dried apple by convection rebuilds during rehydration, unlike those freeze-dried. Based on these results, the authors established a potential mechanism during the rehydration process. The shrunken and compact structure of the product absorbs the water by capillary suction and the wetted dry matrix starts to swell. Then, the swelling causes an increase in volume and the matrix continuity is the source of the mechanical strength of the product. These results suggest that the cell wall network is, at least, partly continuous after convective drying. The lowest values of rehydration ratio indicate extensive cellular and structural tissue damage occurring during the drying process [46]. This leads to denser structures with collapsed capillaries and shrunken pores that lock the water reabsorption of the product [34].
From the sorption curves ( Figure 6), apple dried by air impingement and oven drying adsorbs the same quantity of water vapor over the water activity range investigated (0.8 to 0.96).
is the source of the mechanical strength of the product. These results suggest that the cell wall network is, at least, partly continuous after convective drying. The lowest values of rehydration ratio indicate extensive cellular and structural tissue damage occurring during the drying process [46]. This leads to denser structures with collapsed capillaries and shrunken pores that lock the water reabsorption of the product [34].
From the sorption curves ( Figure 6), apple dried by air impingement and oven drying adsorbs the same quantity of water vapor over the water activity range investigated (0.8 to 0.96). This result suggests that the porous structure between apples dried by impinging jets (fast drying) and by oven drying (slow drying) is the same. This is in accordance with Prawiranto et al. [47] who established, by simulation, that the apple tissue sorption isotherm and its water permeability have a high dependency on the way in which the cellular structure changes during drying. In an earlier review, Lewicki [48] underlined that changes in sorption processes are in relationship with the effective diffusion coefficient. He also clarified that the diffusion coefficients account for the way the molecules are transported within the solid and includes porosity and tortuosity.
It should be specified that the sorption curve for the air impingement drying was determined in optimal conditions to attain the best rehydration (B = 4 × 10 −3 m, H = 1.2 × 10 -2 m, Δ = 5 × 10 −2 m and υ = 99%, trial n°15, Table 2). Nonetheless, this result can be confounded by the fact that the adsorptive capacity appears to be impacted by the porous structure only in low water activity areas (<0.6), as indicated by Saravacos [49].
The water sorption isotherms were modelled using the GAB (Guggenheim-Anderson-de Boer) model according to Equation (2). Values of the corresponding parameters for the two drying processes and the statistical parameters that allow for evaluating the quality of the fit are presented in Table 5. This result suggests that the porous structure between apples dried by impinging jets (fast drying) and by oven drying (slow drying) is the same. This is in accordance with Prawiranto et al. [47] who established, by simulation, that the apple tissue sorption isotherm and its water permeability have a high dependency on the way in which the cellular structure changes during drying. In an earlier review, Lewicki [48] underlined that changes in sorption processes are in relationship with the effective diffusion coefficient. He also clarified that the diffusion coefficients account for the way the molecules are transported within the solid and includes porosity and tortuosity.
It should be specified that the sorption curve for the air impingement drying was determined in optimal conditions to attain the best rehydration (B = 4 × 10 −3 m, H = 1.2 × 10 −2 m, ∆ = 5 × 10 −2 m and υ = 99%, trial n • 15, Table 2). Nonetheless, this result can be confounded by the fact that the adsorptive capacity appears to be impacted by the porous structure only in low water activity areas (<0.6), as indicated by Saravacos [49].
The water sorption isotherms were modelled using the GAB (Guggenheim-Andersonde Boer) model according to Equation (2). Values of the corresponding parameters for the two drying processes and the statistical parameters that allow for evaluating the quality of the fit are presented in Table 5. Statistical parameters reveal a satisfactory modeling of the experimental points (R 2 > 99.9%, RMSE < 7.7 × 10 −3 , χ 2 < 9.2 × 10 −5 and RSS < 5.5 × 10 −4 , Table 5).

Conclusions
This research has focused on studying the impact of design jet parameters (slot width, nozzle-to-surface height, nozzle-to-nozzle spacing, and airflow) on the drying process performances (drying time, and heat and mass transfer coefficients) and the quality attributes of apples (apparent diffusivity, textural properties, and rehydration) using RSM.
A greater significant influence of the examined parameters was shown through the quality of the product compared with the drying process performances. The slot width had a very significant effect on textural properties (hardness and work associated with hardness). The nozzle-to-nozzle spacing and the interaction (slot width × nozzle-to-nozzle spacing) had a very significant effect on the rehydration ratio. The nozzle-to-nozzle spacing and the interaction (nozzle-to-nozzle spacing × airflow) had a significant effect on the work associated with hardness. The interaction (slot width × airflow) had a significant effect on the rehydration ratio. From the ANOVA, the conditions for which the best quality of products, commonly known as 'optimal conditions', were obtained. In this study, optimal conditions were estimated for low textural properties and high rehydration ratios. Low textural properties were preferred to avoid the case hardening phenomenon which occurs during the drying process. The optimal predicted conditions were B = 4 × 10 −3 m, H = 2.4 × 10 −2 m, ∆ = 5 × 10 −2 m, and υ = 80% to obtain the lowest textural properties, and B = 4 × 10 −3 m, H = 1.2 × 10 −2 m, ∆ = 5 × 10 −2 m, and υ = 99% to obtain the highest rehydration ratio. Empirical models between design jet variables and parameters characterizing the quality of the product presented adequate levels of R 2 and adjusted R 2 .
The sorption curve at 323 K of apple dried by air impingement drying in optimal conditions to get the highest rehydration was overlapped by that obtained after oven drying. This result indicates that, after a fast drying (impinging jets), the structure of apple was no further deteriorated than after a slow drying (oven). The GAB model provided a good fit to the experimental data. Funding: This research was funded by CIMS Company (Sablé-sur-Sarthe, France) during the end-ofstudies internship of S. Mlayah (no funding number was assigned to this work).

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The datasets generated during and/or analyzed during this investigation are available from the corresponding author on reasonable request.