Parboiled Paddy Drying with Different Dryers: Thermodynamic and Quality Properties, Mathematical Modeling Using ANNs Assessment

The effect of hybrid infrared-convective (IRC), microwave (MIC) and infrared-convective-microwave (IRCM) drying methods on thermodynamic (drying kinetics, effective moisture diffusivity coefficient (Deff), specific energy consumption (SEC)) and quality (head rice yield (HRY), color value and lightness) characteristics of parboiled rice samples were investigated in this study. Experimental data were fitted into empirical drying models to explain moisture ratio (MR) variations during drying. The Artificial Neural Network (ANN) method was applied to predict MR. The IRCM method provided shorter drying time (reduce percentage = 71%) than IRC (41%) and microwave (69%) methods. The Deff of MIC drying (6.85 × 10−11–4.32 × 10−10 m2/s) was found to be more than the observed in IRC (1.32 × 10−10–1.87 × 10−10 m2/s) and IRCM methods (1.58 × 10−11–2.31 × 10−11 m2/s). SEC decreased during drying. Microwave drying had the lowest SEC (0.457 MJ/kg) compared to other drying methods (with mean 28 MJ/kg). Aghbashlo’s model was found to be the best for MR prediction. According to the ANN results, the highest determination coefficient (R2) values for MR prediction in IRC, IRCM and MIC drying methods were 0.9993, 0.9995 and 0.9990, respectively. The HRY (from 60.2 to 74.07%) and the color value (from 18.08 to 19.63) increased with the drying process severity, thereby decreasing the lightness (from 57.74 to 62.17). The results of this research can be recommended for the selection of the best dryer for parboiled paddy. Best drying conditions in the study is related to the lowest dryer SEC and sample color value and the highest HRY and sample lightness.


Introduction
Rice (Oryza sativa L.) is the major food for over half of the world's population. Fajr is the most consumed and exported rice variety in Iran, although it has a low milling efficiency. The parboiling (A) Soaking: Samples were soaked for 180 min in a bath with well-stirred water at a temperature of 65 ± 0.5 °C. Parboiled paddy (var. Fajr) has best quality at a temperature of 65 °C with a time of 4 min for the parboiled rice [1].
(B) Steaming: Following soaking, the paddy samples (1 kg) were drained and cooled down for 2 h at ambient temperature. The samples were steamed for 10 min and placed on a pot of a metal mesh at 96 °C temperature with 10 L of water boiling [24].
(C) Drying: After steaming, the paddy samples were dried using experimental drying different methods (IRC drying, IRCM drying, and MIC drying) for MC reduction in samples to 15% (d.b.). MC of samples is computable during drying based on sample dried weight using common equations.

Drying Techniques
An experimental setup of an IRCM dryer was manufactured and installed for parboiled paddy drying at the Department of Biosystems Engineering, Mohaghegh Ardabili University, Ardabil, Iran (Supplementary Figure S1). Temperature, velocity, humidity and sample weight were measured by means of a thermometer (Lutron Company, Taipei, Taiwan), anemometer (Lutron-YK Company, Taipei, Taiwan), hygrometer (Testo Company, Lenzkirch, German), and digital balance (A&D, Tokyo, Japan, ±0.01 g), respectively. Each drying experiment was conducted in three replications. The mean relative humidity and temperature of ambient were 24 ± 3% and 23 ± 2 °C, respectively.

Hybrid Infrared-Convection (IRC)
The paddy sample was transferred to a hybrid IRC dryer after steaming. The experiment was performed at different radiation intensity values (0.32 and 0.49 W/cm 2 ), hot air temperatures (40, 50 and (A) Soaking: Samples were soaked for 180 min in a bath with well-stirred water at a temperature of 65 ± 0.5 • C. Parboiled paddy (var. Fajr) has best quality at a temperature of 65 • C with a time of 4 min for the parboiled rice [1].
(B) Steaming: Following soaking, the paddy samples (1 kg) were drained and cooled down for 2 h at ambient temperature. The samples were steamed for 10 min and placed on a pot of a metal mesh at 96 • C temperature with 10 L of water boiling [24].
(C) Drying: After steaming, the paddy samples were dried using experimental drying different methods (IRC drying, IRCM drying, and MIC drying) for MC reduction in samples to 15% (d.b.). MC of samples is computable during drying based on sample dried weight using common equations.

Drying Techniques
An experimental setup of an IRCM dryer was manufactured and installed for parboiled paddy drying at the Department of Biosystems Engineering, Mohaghegh Ardabili University, Ardabil, Iran (Supplementary Figure S1). Temperature, velocity, humidity and sample weight were measured by means of a thermometer (Lutron Company, Taipei, Taiwan), anemometer (Lutron-YK Company, Taipei, Taiwan), hygrometer (Testo Company, Lenzkirch, German), and digital balance (A&D, Tokyo, Japan, ±0.01 g), respectively. Each drying experiment was conducted in three replications. The mean relative humidity and temperature of ambient were 24 ± 3% and 23 ± 2 • C, respectively.

Hybrid Infrared-Convection (IRC)
The paddy sample was transferred to a hybrid IRC dryer after steaming. The experiment was performed at different radiation intensity values (0.32 and 0.49 W/cm 2 ), hot air temperatures (40,50 and 60 • C), and a constant airflow rate (1 m/s). To determine the infrared power, the distance between the emitter and the sample was changed. The distances between the infrared lamp and samples were 20 and 10 cm for 0.32 and 0.49 W/cm 2 , respectively. The samples were spread on a 10 × 10 cm surface.

Infrared Convective-Microwave (IRCM)
In general, tempering duration was eight times longer than the drying time at the ambient temperature [25]. The drying process was performed in two stages: (I) IRC drying was conducted at the first stage (at 40, 50 and 60 • C using 0.32 and 0.49 W/cm 2 and 1 m/s velocity) to reduce MC from 54 to 27% (d.b.); and (II) the second stage included MIC drying (100, 200 and 300 W) to drop sample MC from 27 to 15% (d.b.).

Microwave (MIC)
A programmable microwave oven (M945, Samsung Electronics Inc., Ridgefield Park, NJ, USA) with a 1000 W maximum output at 2450 MHz was applied for the drying experiment. The dimensions of the oven were 330 × 370 × 210 mm 3 . Three power levels (100, 200 and 300 W) were used to study MC reduction in samples from 54 to 15% (d.b.).

Drying Kinetics Modeling
MC values of paddy were calculated and obtained at several equal time intervals during drying: where MC is the moisture content (g water /g dry matter ), W is the amount of evaporated moisture (g), W t is the initial weight of sample (g), W e is the dry matter content of sample (g). The moisture ratio (MR) for parboiled paddy was found using Equation (2) [26]: where MR is the moisture ratio (decimal), M t is the moisture content (% d.b.), M e is the equilibrium moisture content (% d.b.), M i is the initial moisture content (% d.b.). The equilibrium MC of parboiled paddy (long-grain) was calculated using the method of Ondier et al. They suggested the modified equation of Chung-Pfost as [27]: where RH is the relative humidity (%), T c is the hot air temperature ( • C), and A, B, and C for parboiled paddy are 406.92, 23.6172 and 0.2303, respectively. The equilibrium moisture content (M e ) was found to be 7.9% (d.b.).

Mathematical Modeling
Drying curves of parboiled paddy samples were obtained for three different drying processes (IRC, IRCM and MIC drying) and fitted into five different models (Table 1). Where the coefficients a, b, k, and n are empirical constants and coefficients in drying models and t and MR are drying time and moisture ratio, respectively. The curve fitting tool "Curve Expert (Ver. 1.4)" and the technique of nonlinear regression were used to fit the best models into the MR data. This model, which is MR Foods 2020, 9, 86 5 of 17 versus drying time (t), is called kinetic drying. The statistical factors, which include the coefficient of determination (R 2 ), chi-square (χ 2 ) and root mean square error (RMSE), were used for evaluation and comparison of models. The model with the highest R 2 and the lowest χ 2 and RMSE was chosen as the best model for estimating the drying curves [28,29]. Table 1. Mathematical empirical models [30].

Models Equation
Aghbashlo The effective diffusivity of (liquid) moisture transport is a fundamental parameter. Drying rates are represented using empirical models of explicit time-dependent functions. It is a simple and accurate analysis which can be performed using one-dimensional diffusion and Fick's diffusion equation. Fick's second law of diffusion with spherical coordinates was used to calculate the D eff and thermal diffusion. The effective moisture diffusivity determined in the thickness steady-state and transient conditions. The analytical solution for single-layer objects and several geometries was explained by Amiri Chayjan (2011), according to Equation (4) [31]: where D eff is the effective moisture diffusivity (m 2 /s), t is the drying time (s), r is the radius of paddy (m), n is the number of training patterns. Equation (4) can be simplified to the first term of the series for a long drying time: Thew D eff as also obtained typically by means of the slope of Equation (5). The slope (k) was computed using a plotting of ln(MR) versus time.
2.9. Specific Energy Consumption

IRC Drying
This dryer uses different sources of energy including thermal energy (infrared lamps, heaters) and mechanical energy (blowers). Thermal energy was determined using Equation (7) [31]: where EU ther is the thermal energy consumption (kJ), K is the lamp power (W), ρ a is the density of air (kg/m 3 ), v is the air velocity (m/s), A is the cross sectional area of container in which sample was placed (m 2 ), C a is the specific heat of air (kJ/kg • K), ∆T is the temperature difference ( • C). The mechanical energy consumed by the blower was obtained by Equation (8) [32]: where EU mec is the mechanical energy consumption (kJ), P is the pressure difference (mbar), M air is the drying air flow rate (m 3 /s).

MIC Drying
The SEC value of the MIC dryer was computed using Equation (9) [33].
where P mic is the microwave power (kW).

IRCM Drying
The total SEC values in the infrared-convection-microwave dryer was calculated from the sum of SEC of infrared, convection and microwave dryer (thermal energy and mechanical energy) [8]: EU mec = (A·v·ρ a C a ·∆T + K·t + P mic ·t)·3600 The SEC for drying of 1 kg of paddy for these three dryers was calculated using Equation (11) [6]: where SEC is the specific energy consumption (kJ/kg), M W is the weight of lost water (kg).

Artificial Neural Networks (ANN) Modeling
Different multilayer feed-forward-back-propagation (FFBP) and cascade-forward-backpropagation (CFBP) networks were developed and tested, using various hidden layers (1 and 2) and neurons, to determine the best network. Several topologies were investigated using neuron changes in the hidden layer. Training in these networks is an iterative process. The training process is validated if the error between the predicted and desired values is minimum. The optimization method was applied in order to select layers and neurons to evaluate various topologies. An ANN was separately specified for each dryer.

ANN 1: IRC DRYER
Three input parameters were used in the experiments with the IRC dryer. The MR values were measured. Networks with one neuron in the output layer (MR) and three neurons in the input layer (air temperature ( • C), radiation intensity (W/cm 2 ), and drying time (min)) were designed. In this stage, 163 data points were used for ANN. In the first group, 30% (48) was selected for testing, and in the second group, 70% of the data (115) was randomly selected from the entire data pool for training.

ANN 2: IRCM Drying
IRCM drying had four input variables including air temperature ( • C), radiation intensity (w/cm 2 ), microwave power (W), and drying time (min). The output variable for evaluating drying performance was MR. Therefore, the input layer and output layer consisted of four and one neurons respectively. In this research, 70% (315) of the data was randomly applied for training and 30% (135) was used for testing of neural networks.

ANN 3: MIC Drying
For MIC drying, networks with two neurons (microwave power (W) and drying time (min)) in the input layer and one neuron in the output layer (MR) were designed. About 70% (49) of the data were used for training to obtain the best topology.
The ANNs were also trained with Levenberg-Marquardt (LM) and Bayesian Regulation (BR) learning algorithms. Three transfer operations (Pureline, Tansig and Logsig) were applied to reach the suitable topology [29].
where X j are the weighted inputs for each neuron in the jth layer and are calculated as follows: where W ij is the weight between i-th and j-th layers, X j is the j-th input neuron, Y i is the j-th output neuron, b j is the bias of j-th neuron for FFBP and CFBP. The mean squared error (MSE) was computed using Equation (18) [34]: where S k is the network output for kth pattern, T k is the target output for kth pattern, N 0 is the number of output neurons. The efficiency of the ANN models were quantified using the R 2 and the MAE between the real and predicted values that were calculated according to the following equations [28]: MAE-mean absolute error.

Head Rice Yield (HRY)
After four days of drying, paddy was de-husked using laboratory stake rubber roller rice husker instrument (Satake Company, Hiroshima, Japan). Then the polishing of samples was performed by a polisher (Satake Company, Hiroshima, Japan) for 90 S. The HRY was computed with dividing of whole grains weight to paddy; the average of triplicates was considered [35].

Color Value and Lightness
To measure the color and lightness, precision color reader (Model 4510, Reston Company, Fairfax, VA, USA) was used. Before measurement, the instrument was calibrated with white and black plates. The results were expressed in term of L*, a*, b* values. L* shows brightness from black to white, positive and negative "a*" values represent redness/greenness, and "b*" value represents yellowness/blueness. The color measurement were done in ten replication and the average values were subjected to further analysis. The parboiled rice color value (B) was calculated using Equation (19) [36].

Drying Time
The values of moisture ratio versus drying time for different drying conditions are shown in Supplementary Figure S2. As can be seen, MR decreases continuously and exponentially with increasing drying time.

Hybrid IRC Drying
Moisture removal from samples was accelerated when air temperature and radiation intensity increased. As a result, the drying time decreased from 68 to 40 min.

IRCM Drying
The drying time of parboiled paddy samples using 0.32 and 0.49 W/cm 2 , 40-60 • C and 100-300 W conditions varied between 24 and 82 min. The drying time reduced with increase in the radiation intensity, air temperature, and MIC power.

Microwave Drying
The results showed that drying time decreased (72-22 min) when MIC power increased. This result is in agreement with the results of other studies [29]. In addition, the increased radiation intensity or air temperature or microwave power also increased the product surface temperature and moisture loss, leading to faster drying [37].

Mathematical Modeling
Five various empirical models were applied to analyze the MR data of parboiled paddy in various dryers and results are presented in Supplementary Table S1. For all dryers, Aghbashlo's model had the highest R 2 (0.99) and the lowest RMSE (0.0194) and χ 2 (0.0004) values. Therefore, Aghbashlo's model was chosen as the best model for explaining the drying properties of parboiled paddy in IRMC dryer. These results were in agreement with other researcher's reported data [32,38].

Effective Moisture Diffusivity (Deff)
D eff values of the parboiled paddy samples dried by different drying techniques are shown in Figure 2. In hybrid IRC drying, the lowest (1.32 × 10 −10 m 2 /s) and highest (1.87 × 10 −10 m 2 /s) D eff values belonged to the samples dried at 40 • C using 0.32 W/cm 2 radiation and at 60 • C using 0.49 W/cm 2 radiation, respectively. In IRCM drying, the lowest (1.58 × 10 −11 m 2 /s) and highest (2.31 × 10 −11 m 2 /s) D eff values belonged to the samples dried at 40 • C using 0.32 W/cm 2 and 100 W and at 60 • C using 0.49 W/cm 2 and 300 W, respectively. At the same time, in MIC drying, the lowest (6.85 × 10 −11 m 2 /s) and highest (4.32 × 10 −10 m 2 /s) D eff were reported for samples dried at 100 and 300 W power levels, respectively. For all dryers, the D eff values increased with increase in the drying temperature, radiation intensity, and MIC power. High temperatures increased the molecular activity of water and caused higher diffusivity. According to Figure 2, D eff in IRCM drying was lower than in hybrid and microwave dryers. Aydogdu et al. [39] proved that effective moisture diffusivity for eggplant in HA drying (5.07 × 10 −10 m 2 /s) is lower than that observed in microwave-infrared drying (7.10 × 10 −9 -1.44 × 10 −8 m 2 /s) [39]. However, the highest D eff was observed in samples dried using microwave drying compared to other drying techniques, because the high MIC power can accelerate the evaporation of water molecules in paddy samples and provide a faster decrease in the rice MC corresponding to the higher value of D eff [40]. D eff results were comparable with the values reported by other authors [41]. The reported D eff values were within the general range of 10 −8 to 10 −12 m 2 /s for agricultural and food products [29].

Specific Energy Consumption (SEC)
The values of SEC for different dryers are shown in Figure 3. It was observed that SEC values decreased as radiation intensity, drying temperature or microwave power increased. This is because low radiation intensity, air temperature, and microwave power can cause a relative decrease in Deff, leading to an increase in SEC values. For IRC drying, the highest (29.06 MJ/kg) and lowest (19.55 MJ/kg) SEC values were obtained at 40 °C using 0.32 W/cm 2 and at 60 °C using 0.49 W/cm 2 , respectively. In IRCM drying, the highest (40.23 MJ/kg) and lowest (22.88 MJ/kg) SEC values were observed at 40 °C using 0.32 W/cm 2 and 100 W and at 60 °C using 0.49 W/cm 2 and 300 W, respectively. According to Figure 3, the lowest SEC in microwave drying was 0.457 kWh/kg, which occurred at 100 W of microwave power, whereas the SEC value of 0.575 MJ/kg was observed when using the microwave power of 300 W. In other words, the highest SEC was 1.25 times greater than the lowest energy required. The reduction of SEC was related to the effect of volumetric heating of MIC, which substantially decreased the duration of drying. These SEC values are comparable with the values reported by other researchers [42,43]. The SEC in the IRCM dryer was higher than that in the hybrid IRC and MIC dryers. Similar results were reported by other researchers [44].

Specific Energy Consumption (SEC)
The values of SEC for different dryers are shown in Figure 3. It was observed that SEC values decreased as radiation intensity, drying temperature or microwave power increased. This is because low radiation intensity, air temperature, and microwave power can cause a relative decrease in D eff , leading to an increase in SEC values. For IRC drying, the highest (29.06 MJ/kg) and lowest (19.55 MJ/kg) SEC values were obtained at 40 • C using 0.32 W/cm 2 and at 60 • C using 0.49 W/cm 2 , respectively. In IRCM drying, the highest (40.23 MJ/kg) and lowest (22.88 MJ/kg) SEC values were observed at 40 • C using 0.32 W/cm 2 and 100 W and at 60 • C using 0.49 W/cm 2 and 300 W, respectively. According to Figure 3, the lowest SEC in microwave drying was 0.457 kWh/kg, which occurred at 100 W of microwave power, whereas the SEC value of 0.575 MJ/kg was observed when using the microwave power of 300 W. In other words, the highest SEC was 1.25 times greater than the lowest energy required. The reduction of SEC was related to the effect of volumetric heating of MIC, which substantially decreased the duration of drying. These SEC values are comparable with the values reported by other researchers [42,43]. The SEC in the IRCM dryer was higher than that in the hybrid IRC and MIC dryers. Similar results were reported by other researchers [44]. Figure 3, the lowest SEC in microwave drying was 0.457 kWh/kg, which occurred at 100 W of microwave power, whereas the SEC value of 0.575 MJ/kg was observed when using the microwave power of 300 W. In other words, the highest SEC was 1.25 times greater than the lowest energy required. The reduction of SEC was related to the effect of volumetric heating of MIC, which substantially decreased the duration of drying. These SEC values are comparable with the values reported by other researchers [42,43]. The SEC in the IRCM dryer was higher than that in the hybrid IRC and MIC dryers. Similar results were reported by other researchers [44].  Table 2 shows the best ANN topologies, threshold functions, and different algorithms applied for predicting MR in the drying of parboiled paddy samples using different dryers. The highest R 2  Table 2 shows the best ANN topologies, threshold functions, and different algorithms applied for predicting MR in the drying of parboiled paddy samples using different dryers. The highest R 2 between ANN-predicted MR values of paddy samples in different dryers and their measured values for training and testing datasets was >0.99.

IRC Drying
As shown in Table 2, the neural network with the 3-8-8-1 topology, TAN, PUR, and TAN transfer functions, and the LM training algorithm (trainlm) had the lowest MAE values (0.0059 for training and 0.0048 for testing) and the highest R 2 (0.9992 for training and 0.9993 for testing).

IRCM Drying
The best model for prediction of MR belonged to the cascade-forward-back-propagation (CFBP) network, the 4-10-10-1 topology and the TANSIG threshold function with the BR algorithm in the first strategy.

Microwave Drying
The widely used topologies and threshold functions have a good training process. The feed-forward back-propagation (FFBP) structure with TAN-TAN-TAN threshold function and the BR algorithm with 2-10-10-1 topology had the lowest MSE (0.00065), MAE (0.0054 for test and 0.0073 for training) values and the highest R 2 (0.9990 for testing and 0.9989 for training). These findings are in agreement with those reported by another author [7]. As shown in Figure 4, HRY values were between 60.2 and 74.07%. The maximum HRY (74.08%) was obtained when rice was dried using IRCM drying method (at 50 • C and 300 W). Taghinezhad and Brenner [1] found the maximum HRY (67.37%) after steaming to be 65 • C with a soaking time of 6 min [1]. The minimum HRY was observed (60.21%), when paddy was dried using MIC drying method (100 W). According to Figure 4, increasing drying temperature (from 40 • C to 60 • C) and microwave power (from 100 to 300 W) caused an increase of HRY. The amount of HRY is likely to be related to the higher temperatures [45]. Similar reports related to the effect of the drying conditions on the HRY have also been presented by other researchers [46]. They reported that the HRY of long-grain SP 1 parboiled rice was relatively higher when drying temperature increased. The experimental results for HRY agreed with results reported by other researchers [47]. They reported the HRY values in the range 60-80% for parboiled rice (KDML 105 paddy). Different drying systems have significant influence on the extent of HRY for parboiled rice as shown in Figure 4. Therefore, it seems to be established that an increase in drying temperature and power at the range we investigated in the present research may lead to an increase in the degree of starch gelatinization (DSG). This is because high temperature and radiation penetrates the rice grain kernel and leads to greater DSG [46]. Hence the appropriate condition of drying can lead to increase in HRY [36]. Also, the increase in HRY could also be related to facilitating the separation of the gelatinized kernels from the husk following the drying. As a result, the milling becomes easier following this separation of the husk from the kernel [48]. As shown in Figure 4, HRY values were between 60.2 and 74.07%. The maximum HRY (74.08%) was obtained when rice was dried using IRCM drying method (at 50 °C and 300 W). Taghinezhad and Brenner [1] found the maximum HRY (67.37%) after steaming to be 65 °C with a soaking time of 6 min [1]. The minimum HRY was observed (60.21%), when paddy was dried using MIC drying method (100 W). According to Figure 4, increasing drying temperature (from 40 °C to 60 °C) and microwave power (from 100 to 300 W) caused an increase of HRY. The amount of HRY is likely to be related to the higher temperatures [45]. Similar reports related to the effect of the drying conditions on the HRY have also been presented by other researchers [46]. They reported that the HRY of long-grain SP 1 parboiled rice was relatively higher when drying temperature increased. The experimental results for HRY agreed with results reported by other researchers [47]. They reported the HRY values in the range 60-80% for parboiled rice (KDML 105 paddy). Different drying systems have significant influence on the extent of HRY for parboiled rice as shown in Figure 4. Therefore, it seems to be established that an increase in drying temperature and power at the range we investigated in the present research may lead to an increase in the degree of starch gelatinization (DSG). This is because high temperature and radiation penetrates the rice grain kernel and leads to greater DSG [46]. Hence the appropriate condition of drying can lead to increase in HRY [36]. Also, the increase in HRY could also be related to facilitating the separation of the gelatinized kernels from the husk following the drying. As a result, the milling becomes easier following this separation of the husk from the kernel [48].

Color Value
The color of rice is a very important indicator for measuring the physical properties of processed rice. Figure 5 signifies that color value is greatly affected by different drying methods. The lightness value ranges from 18.08 to 19.63. The darkest color value was observed for IRC dryer (60 °C). The discoloration of rice was caused by the reaction of released sugar with amino acid of grain treated under high temperature [49]. Similar results have been presented by other researchers [43] studying the effect of the different parboiling conditions on the color value. Also, many researchers evaluated the parboiled rice color [50,51], noting that final color changes during parboiling can be explained as nonenzymatic browning (Maillard reactions) and the drying variables influenced the intensity of color. In other words, the color parameters revealed that during parboiling, yellow and red bran pigments

Color Value
The color of rice is a very important indicator for measuring the physical properties of processed rice. Figure 5 signifies that color value is greatly affected by different drying methods. The lightness value ranges from 18.08 to 19.63. The darkest color value was observed for IRC dryer (60 • C). The discoloration of rice was caused by the reaction of released sugar with amino acid of grain treated under high temperature [49]. Similar results have been presented by other researchers [43] studying the effect of the different parboiling conditions on the color value. Also, many researchers evaluated the parboiled rice color [50,51], noting that final color changes during parboiling can be explained as non-enzymatic browning (Maillard reactions) and the drying variables influenced the intensity of color. In other words, the color parameters revealed that during parboiling, yellow and red bran pigments diffused from the bran into the endosperm. Bran pigments diffused into the endosperm affect the color of parboiled rice [52].
Foods 2020, 9, x FOR PEER REVIEW 13 of 18 diffused from the bran into the endosperm. Bran pigments diffused into the endosperm affect the color of parboiled rice [52].

Lightness
The lightness value ranges from 57.74 to 62.17; these values were relatively similar to those reported by others [24]. The darkest rice was observed when rice was treated at 60 °C drying temperature (under IRC dryer). Figure 6 represents the effect of different drying methods on the lightness of parboiled rice. The degree of lightness decreased with the severity of drying process (increase of drying temperature (from 40 to 60 °C) and microwave power (from 100 to 300 W)). These findings are in agreement with those of other researchers [50,51]. Consequently, it can be explained that rice should not be treated under severe drying conditions. Higher temperature application could affect lightness of rice, leading to poor quality of rice and less demand in the market [53].

Conclusions
Drying is one of the commonly used methods for reducing product moisture, especially for long storage times. In this research, the effects of drying different methods on the thermodynamic and quality characteristics of parboiled paddy were evaluated in this study. Five mathematical methods were applied to predict the moisture ratio of parboiled rice. The results proved that with increasing The same letter over column shows that the mean amount had no significant difference (p < 0.05) based on Duncan's test. The error bars represent standard error of the means.

Lightness
The lightness value ranges from 57.74 to 62.17; these values were relatively similar to those reported by others [24]. The darkest rice was observed when rice was treated at 60 • C drying temperature (under IRC dryer). Figure 6 represents the effect of different drying methods on the lightness of parboiled rice. The degree of lightness decreased with the severity of drying process (increase of drying temperature (from 40 to 60 • C) and microwave power (from 100 to 300 W)). These findings are in agreement with those of other researchers [50,51]. Consequently, it can be explained that rice should not be treated under severe drying conditions. Higher temperature application could affect lightness of rice, leading to poor quality of rice and less demand in the market [53].
Foods 2020, 9, x FOR PEER REVIEW 13 of 18 diffused from the bran into the endosperm. Bran pigments diffused into the endosperm affect the color of parboiled rice [52].

Lightness
The lightness value ranges from 57.74 to 62.17; these values were relatively similar to those reported by others [24]. The darkest rice was observed when rice was treated at 60 °C drying temperature (under IRC dryer). Figure 6 represents the effect of different drying methods on the lightness of parboiled rice. The degree of lightness decreased with the severity of drying process (increase of drying temperature (from 40 to 60 °C) and microwave power (from 100 to 300 W)). These findings are in agreement with those of other researchers [50,51]. Consequently, it can be explained that rice should not be treated under severe drying conditions. Higher temperature application could affect lightness of rice, leading to poor quality of rice and less demand in the market [53].

Conclusions
Drying is one of the commonly used methods for reducing product moisture, especially for long storage times. In this research, the effects of drying different methods on the thermodynamic and quality characteristics of parboiled paddy were evaluated in this study. Five mathematical methods were applied to predict the moisture ratio of parboiled rice. The results proved that with increasing

Conclusions
Drying is one of the commonly used methods for reducing product moisture, especially for long storage times. In this research, the effects of drying different methods on the thermodynamic and quality characteristics of parboiled paddy were evaluated in this study. Five mathematical methods were applied to predict the moisture ratio of parboiled rice. The results proved that with increasing of radiation intensity or air temperature or microwave power increased the product surface temperature and moisture loss. As a result, it led to faster drying. Aghbashlo's model was the best model for the prediction of parboiled paddy MR. It was able to describe the thin-layer thermodynamic characteristics of samples in the three dryers. The highest D eff was obtained under the microwave drying conditions. The lowest SEC was calculated in MIC drying. ANN be able to predict the MR with high accuracy. The physical properties were influenced by means of different drying methods. The obtained results indicate that the amount of HRY and color value hardness increased during drying, while the amount of lightness decreased. The results of this research can be used for the selection of the best dryer in the parboiling industry of paddy. The best drying conditions were related to the highest HRY and sample lightness, and the lowest dryer SEC and sample color value.