Experimental Study on the Ventilation Resistance Characteristics of Paddy Grain Layer Modelled with Response Surface Methodology

: The ventilation resistance of air ﬂow through rice grain layers is one of the key parameters affecting drying uniformity as well as the energy consumption of the drying process. In order to reveal the variation of characteristics of the ventilation resistance with paddy grain moisture content, the air velocity and the bed layer depth are needed. A second order model was ﬁtted to pressure drop using the response surface methodology and the results are compared with those of the Ergun model. The results showed that the pressure drop increases with the increase of paddy grain moisture content, air velocity and the bed layer depth, and the interactions between the air velocity and the bed layer depth have the most signiﬁcant inﬂuence on the pressure drop. Moreover, a second-order polynomial pressure drop model based on RSM was established and compared with the Ergun model. The results showed that the pressure drop model established by RSM is similar to that of the Ergun model.


Introduction
The pressure drop of the grain bed layer is the key parameter for designing and optimizing the structure of the drying chamber, and it is also one of the most important references for selecting the power of a fan [1,2]. During the drying process, owing to the bed layer porosity as well as the air viscosity, there must be a pressure drop when the air flows through the bed layer. A high-pressure drop may make it difficult for the air to flow through the bed layer and further lead to the increase of non-uniformity of the drying process, while a low-pressure drop may lead to wasting the power of the fan [3]. Accordingly, it is necessary to find out a reasonable pressure drop under the premise of ensuring product quality and maintaining a reasonable power demand of the fan.
In the last few decades, researchers have done many studies on the ventilation resistance of the grain bed layer by applying the traditional pressure drop models [4][5][6][7]. For example, Gunasekaran and Jackson [8] conducted an experiment on the flow resistance of sorghum bed layer by applying the Ergun model, the resistance of sorghum moisture content of 16.5, 18.5, and 23 %w.b. were determined under the grain bed depth of 150 to 1200 mm with the airflow ranges of 0.05 to 0.3 m/s. The results showed that the resistance of airflow increases with the increase of air velocity and bed layer depth, while resistance increases with the decrease of moisture content. Li et al. [9] applied the Shedd model to simulate the flow resistance of the hot air flow through the paddy grain bed layer; they found that the flow resistance of the hot air increases with the increase of hot air velocities. Zhang et al. [10] investigated the pressure drop of the air (with velocity ranges of 0.1 to 0.6 m/s) flow through the bed layer with depth ranges of 100 to 1000 mm, and found that the simulation accuracy of the pressure drop using Ergun model is reliable only when the air velocity is under 0.2 m/s. It can be summarized from the above analyses that the flow resistance of the grain bed layer are mainly determined by the traditional models and the characteristics of the air flow resistance are investigated by using single factor test. Though the works provided some valuable references for guiding the practical production, the air flow resistance prediction accuracy based on traditional models should be further improved and the interaction of the parameters on the flow resistance should be further investigated.
Respond surface methodology (RSM) is one of the most effective modeling tools to investigate the interactions of two or more variables on the responses, which is widely adopted to optimize the process parameters in the agriculture field [11,12]. Zeng et al. [13] conducted an experimental study on the vacuum drying extraction technology of konjac glucomannan using RSM, and an artificial neuron network (ANN) was also adopted and then compared with the RSM. The results showed that the model established by RSM is slightly better than that of an ANN. In order to improve the dynamic performance of orchard ditching fertilizer and to avoid resonance, Liu et al. [14] adopted the RSM to optimize the structure of the rack of orchard ditching-fertilizer machine. The results showed that the optimized structure has a better first-order modal frequency of 38.31 Hz than that of the traditional structure, and the optimized frequency is far away from the input frequency (35 Hz) of a tractor, indicating that the frequency of the structure optimized by RSM can avoid the resonance. Moreover, Kowalczyk et al. [15] adopted the RSM to optimize the production process of a micro-elemental feed additive for laying hens, the most favorable enrichment conditions of the process were determined by the model established by the RSM. It can be summarized that the RSM is widely used in the agricultural field, and it can be also adopted to model the ventilation resistance of paddy grain layer.
The present work adopted the RSM to model the ventilation resistance of paddy grain layer. In detail, an experiment followed by Box-Behnken Design was conducted to investigate the influences of the independent variable including paddy grain moisture content (Mc), air velocity (v air ), and paddy layer depth (L) on the pressure drop (∆p). The influences of the single factor as well as the interactions between any two independent variables on ∆p were investigated and a ternary quadratic model for ∆p was established and compared with the traditional Ergun model.

Materials and Methods
To fully understand the variation rules of ventilation resistance of the air flowing through paddy grain layer and further specify the optimization strategy for deep bed paddy drying process. A laboratory-scale experimental test was conducted based on the engineering application consideration; the detailed analysis steps are as follows.

Materials
The paddy grains sample with initial moisture content of 12.2 %w.b. are purchased from a local farmer at Leizhou Guangdong Province, China. The samples are, respectively manual humidified to 13.6 %w.b., 14.8 %w.b., 16.2 %w.b., 17.8 %w.b., 18.9 %w.b., 21.3 %w.b., 22.9 %w.b., 24.5 %w.b., and 25.6 %w.b. according to the methodology introduced in the literature [16]. Additionally, the prepared samples are sealed in plastic bags and stored in a refrigerator at approximately 4 • C for about 48 h prior to the experiments [13].

Experimental Apparatus
As can be seen from Figure 1, the apparatus mainly consists of 6 parts including the centrifugal fan, frequency converter, air supply duct, screen mesh, test pipeline, and the fixing bracket. There are 10 pressure taps with a diameter of 10 mm on the test pipeline, the height of the test pipeline is 1100 mm, and the distance between any two adjacent pressure taps is 100 mm. Moreover, there is one air velocity tap closed to connection between the test pipeline and air supply duct and the connection is sealed by a rubber ring and silicone [17]. During the tests, the three levels of air velocities (0.1, 0.35, and 0.6 m/s) are produced by the centrifugal fan connected with a frequency converter. The air flows through the air supply duct, then flows cross the screen mesh, further flows cross the grain layer, and finally flows into the ambient space. The v air was controlled by the frequency converter. The inlet v air was measured by an anemometer. The air pressure was measured by a digital pressure gauge. The details of the instruments used in the tests are tabulated in Table 1. Compared with the industrial dryer-although the physical structure of the laboratory-scale apparatus is different, the grain bed is dynamic (actually for a single grain bed layer, and the velocities of the grain bed are very low during the cyclic drying process-the ratio of air volume to the mass of the bed layer for the industrial dryer and the apparatus are the same [18]. This indicates that the main results obtained by the laboratory-scale apparatus can provide some references for the industrial dryer, especially for power selection of the fan and to ascertain of the bed layer depth in the design link of the industrial dryer. the height of the test pipeline is 1100 mm, and the distance between any two adjacent pressure taps is 100 mm. Moreover, there is one air velocity tap closed to connection between the test pipeline and air supply duct and the connection is sealed by a rubber ring and silicone [17]. During the tests, the three levels of air velocities (0.1, 0.35, and 0.6 m/s) are produced by the centrifugal fan connected with a frequency converter. The air flows through the air supply duct, then flows cross the screen mesh, further flows cross the grain layer, and finally flows into the ambient space. The vair was controlled by the frequency converter. The inlet vair was measured by an anemometer. The air pressure was measured by a digital pressure gauge. The details of the instruments used in the tests are tabulated in Table 1. Compared with the industrial dryer-although the physical structure of the laboratory-scale apparatus is different, the grain bed is dynamic (actually for a single grain bed layer, and the velocities of the grain bed are very low during the cyclic drying process-the ratio of air volume to the mass of the bed layer for the industrial dryer and the apparatus are the same [18]. This indicates that the main results obtained by the laboratory-scale apparatus can provide some references for the industrial dryer, especially for power selection of the fan and to ascertain of the bed layer depth in the design link of the industrial dryer.  In the present work, three different paddy grain Mc of 12.2 %w.b., 18.9 %w.b., and 25.6 %w.b.; three different vair of 0.10, 0.35, and 0.60 m/s; three different L of 0.10, 0.55, and 1.0 m, were considered to be the independent variables, while the Δp was regarded as the response. Owing to the fact that the Box-Behnken Design has the advantages of the less experimental operations comparing with Central Composite Design [19], and by  In the present work, three different paddy grain Mc of 12.2 %w.b., 18.9 %w.b., and 25.6 %w.b.; three different v air of 0.10, 0.35, and 0.60 m/s; three different L of 0.10, 0.55, and 1.0 m, were considered to be the independent variables, while the ∆p was regarded as the response. Owing to the fact that the Box-Behnken Design has the advantages of the less experimental operations comparing with Central Composite Design [19], and by comprehensive considering that the mass of each test is huge, the Box-Behnken Design method was adopted to analyze the effects of the interaction of independent variables (Mc, v air, and L) on the response (∆p) and further optimize the drying process. The experimental design representing the natural and coded values of independent variables are shown in Table 2. In order to develop the mathematical relations between the independent variables Appl. Sci. 2021, 11, 7826 4 of 15 and the response, the second-order polynomial model (Equation (1)) was adopted to establish the RSM model [20].
where Y is the response, X i and X j are the independent variables affecting the response, β 0 , β i , and β ij are the regression coefficients for intercept quadratic liner and interaction terms.

Model Validation
In statistics, the p-value is the probability of obtaining results as extreme as the results observed from a statistical hypothesis test, assuming that the null hypothesis is correct. The p-value is used as an alternative to rejection points to provide the smallest level of significance where the null hypothesis would be rejected [21,22]. A smaller p-value means stronger evidence for the alternative hypothesis [23].
Based on the analysis above, the p-value is used to validate hypothesis test of the model. The null hypothesis is defined as the response (∆p) that cannot be explained by the obtained model, while the alternative hypothesis is defined as the response (∆p) can be explained by the obtained model.

RSM Model Validation Metrics
In the present work, the statistical indexes including determination coefficient (R 2 ), mean absolute error (MAE), mean square error (MSE), and root mean square (RMSE) were adopted to evaluate the fitting performance of the ∆p model developed by the RSM and the Ergun model, which can be, respectively, calculated followed by the following Equations (2)-(5) [24]: where, Y pre,i is the value predicted by the models, Y exp,i is the experimental value, − Y exp is the mean of the experimental values, and n is the number of the data points. Moreover, the experimental results are tabulated in Table 3.

ANOVA Validation: Anderson-Darling (AD) Test for Normality
According to some of the literature, the distribution of residuals for a set of data can reflect the validity of a predicted model [25]. In detail, if the residuals follow a normal distribution (occurrences are random), the established model has a well-predicting capacity. The Anderson-Darling (AD) test is the statistical tool used to quantify the deviation for a set of residuals from a normal distribution. The validity of the distribution of residuals in the AD test at a 5% level of significance was used to confirm the model accuracy and to determine if a sample data set was taken from a population with a specific distribution [26]. In the present work, the AD test was conducted by the software Minitab 17.0, the null hypothesis is defined as the distribution of the residuals that do not follow the normal distribution, and the hypothesis can be rejected if the test statistic (AD value) is greater than a critical value [27].

Single Factor Experiments
In order to investigate influence of each independent variable on the response, the ∆p under ten levels of Mc, v air , and L were, respectively, measured followed by the 3 sets of experimental design, the single factor experiment design and the corresponding results are tabulated in Table 4.

Ergun Model
Ergun model is one of most effective tools for revealing the relationship between fluid properties and porous media properties, which has been widely used in Hydrodynamics, thermodynamics, and engineering [28]. In order to investigate the ventilation resistance characteristics of the air flow through the paddy bed layer as well as to verify the ∆p model developed by the RSM. The most widely used Ergun model (∆p Eur ) was applied in the present work, which can be expressed as following expressions [29][30][31]: where µ is the dynamic viscosity of the air, 1.81 × 10 −5 kg/(m·s) [32]; ρ is the density of the air, ρ = 1.205 kg/m 3 [33]; ε is the porosity of the bed layer; and d g is the equivalent diameter of a single paddy grain, m. The ε for different Mc levels of the paddy grains can be determined by a self-developed porosimeter; the details of the porosimeter can be found in our previous work [34].
On the other hand, the d g for a single paddy grain with different Mc can be calculated according to the length (L 0 ), width (W 0 ), and height (H 0 ) of the paddy grain. The calculation methodology introduced by Mohsenin et al. [35] are expressed as follows: Each parameter used for calculating the d g was measured 3 times and the average value was used for the calculation.

ANOVA Analysis
In the present work, the software Design Expert 8.0.6 was applied to analyze the experimental results and the ANOVA analysis was employed to analyze the significance of the Mc, vair, and L on the Δp. The results are tabulated in Table 5. It can be seen from the table that the sum of squares, mean square, F-value, and p-value for each of the terms of the model were calculated, of which the F-value and p-value are the criteria for evaluating the fitting performance of the specific term [36]. Specifically, the larger the F-value is, the more significant the term is and the better the fitting performance is; on the other hand, the p-value with the range of (0.01 < p < 0.05) indicates the term is significant (characterized by the symbol "*") while a p-value less than 0.01 indicates that the term is extremely significant (characterized by the symbol "**").

ANOVA Analysis
In the present work, the software Design Expert 8.0.6 was applied to analyze the experimental results and the ANOVA analysis was employed to analyze the significance of the Mc, v air , and L on the ∆p. The results are tabulated in Table 5. It can be seen from the table that the sum of squares, mean square, F-value, and p-value for each of the terms of the model were calculated, of which the F-value and p-value are the criteria for evaluating the fitting performance of the specific term [36]. Specifically, the larger the F-value is, the more significant the term is and the better the fitting performance is; on the other hand, the p-value with the range of (0.01 < p < 0.05) indicates the term is significant (characterized by the symbol "*") while a p-value less than 0.01 indicates that the term is extremely significant (characterized by the symbol "**"). Note: The symbol "**" indicates the term is extremely significant.
As can be seen from Table 5, the ∆p model F value of 148.56 corresponds to a low p-value (p < 0.0001), which implies that the model was significant as the test failed to accept the null hypothesis that variations in the response could not be explained by the model. This indicates that the ∆p model developed by the RSM is significant, and the experimental ∆p can be explained by the obtained model. Moreover, the values of R 2 , Adj R 2 , and Pred R 2 are similar and all close to 1 [37], which also indicates that the ∆p model is in a well-fitting performance and can be used to predicted the ∆p. The value of the C.V.% term is 10.98%, indicating that there are only 10.98% of the experimental data that cannot be explained by the ∆p model developed by the RSM, in statistics, if the value of the C.V.% is greater than 15%, the data may be abnormal and should be eliminated [19].

Anderson-Darling Normality Test Results
In the present work, the residuals are the difference between the experimental ∆p and the predicted ∆p by the established model (Equation (2)). The AD test was conducted by the software Minitab 17.0, and the probability distribution plot of the residuals is shown in the following Figure 3. As can be seen from the figure, the obtained AD value (0.484) is less than a critical value of 0.752, and the associated p-value (0.594) of the AD test is significant at a 5% level (p-value > 0.05) [25]. The results indicate that the residuals follow a normal distribution, and the deviation between the experimental ∆p and predicted ∆p are random [26]. Accordingly, it can be concluded that the null hypothesis can be rejected and the model predictions are correlated with the experimental data over the factor-space under evaluation in the present work.

Anderson-Darling Normality Test Results
In the present work, the residuals are the difference between the experimental Δp and the predicted Δp by the established model (Equation (2)). The AD test was conducted by the software Minitab 17.0, and the probability distribution plot of the residuals is shown in the following Figure 3. As can be seen from the figure, the obtained AD value (0.484) is less than a critical value of 0.752, and the associated p-value (0.594) of the AD test is significant at a 5% level (p-value > 0.05) [25]. The results indicate that the residuals follow a normal distribution, and the deviation between the experimental Δp and predicted Δp are random [26]. Accordingly, it can be concluded that the null hypothesis can be rejected and the model predictions are correlated with the experimental data over the factor-space under evaluation in the present work.

The Interactions of the Independent Variables on the Pressure Drop
Based on the ANOVA analysis in Section 3.1, the Mc·v air , Mc·L, and v air ·L, respectively, have significant influences on the ∆p. The response surface as well as its corresponding contour plots are depicted in the following Figures 4-6. As can be seen from Figure 4a, the ∆p increases with the increase of v air under constant Mc, which might be due to the fact that the inertia resistance of the air increases with the increase of v air , just as reported by Zhang et al. [10]. On the other hand, the ∆p increases with increase of Mc under constant v air , which might be caused by the increased ε in high Mc, just as reported by Li et al. that the ε of paddy layer increases with increase of paddy moisture content [38]. Though the suitable ∆p for the engineering application should be further investigated by comprehensively considering the power of the fan, the contours shown in Figure 4b can be the reference for the single factor analysis.    [38]. Though the suitable Δp for the engineering application should be further investigated by comprehensively considering the power of the fan, the contours shown in Figure 4b can be the reference for the single factor analysis.    The interaction of vair and L on the Δp are shown in the following Figure 6. As can be seen from Figure 6a

The Comparison of the Δp Model Developed by RSM and Ergun Model
In order to verify the feasibility of the Δp model developed by RSM, the comparison analysis focused on the RSM experiments for the Δp model developed by RSM and Ergun model was conducted, and the statistical indexes including R 2 , MAE, MSE, and RMSE were calculated. The detailed analysis is in Section 3.4.1.  The interaction of vair and L on the Δp are shown in the following Figure 6. As can be seen from Figure 6a that the Δp varies from the minimum 17.11 Pa to the maximum 689.56 Pa, the variation range is much larger than that of the interaction terms of Mc•L and Mc•vair, indicating that the influence of term L•vair on Δp is more significant than that of Mc•L and Mc•vair, just as analyzed in Section 3.1. Accordingly, in order to reach the objective of the efficient and energy-saving drying, efforts should be firstly made to optimize the drying parameters related to the L and vair in the engineering drying application. Especially when the air with high vair flow through the bed layer with large L, the Δp rapidly increases with the interaction term of L•vair (as shown in Figure 6b, in other words, the larger driven force for supplying the drying medium is needed when the L•vair is high.

The Comparison of the Δp Model Developed by RSM and Ergun Model
In order to verify the feasibility of the Δp model developed by RSM, the comparison analysis focused on the RSM experiments for the Δp model developed by RSM and Ergun model was conducted, and the statistical indexes including R 2 , MAE, MSE, and RMSE were calculated. The detailed analysis is in Section 3.4.1.  The interaction of v air and L on the ∆p are shown in the following Figure 6. As can be seen from Figure 6a that the ∆p varies from the minimum 17.11 Pa to the maximum 689.56 Pa, the variation range is much larger than that of the interaction terms of Mc·L and Mc·v air , indicating that the influence of term L·v air on ∆p is more significant than that of Mc·L and Mc·v air , just as analyzed in Section 3.1. Accordingly, in order to reach the objective of the efficient and energy-saving drying, efforts should be firstly made to optimize the drying parameters related to the L and v air in the engineering drying application. Especially when the air with high v air flow through the bed layer with large L, the ∆p rapidly increases with the interaction term of L·v air (as shown in Figure 6b, in other words, the larger driven force for supplying the drying medium is needed when the L·v air is high.

The Comparison of the ∆p Model Developed by RSM and Ergun Model
In order to verify the feasibility of the ∆p model developed by RSM, the comparison analysis focused on the RSM experiments for the ∆p model developed by RSM and Ergun model was conducted, and the statistical indexes including R 2 , MAE, MSE, and RMSE were calculated. The detailed analysis is in Section 3.4.1.

The Comparison of the ∆p in the RSM Experiments
Based on the ANOVA analysis results shown in Table 5, the second-order polynomial expression of the ∆p model developed by the RSM (∆p RSM ) are expressed as following Equation (8) The 17 sets of independent variables tabulated in Table 5 were substituted into the above ∆p RSM model and the predicted ∆p RSM values were obtained. On the other hand, the 17 sets of independent variables were also substituted into the ∆p Eur model shown in Equation (2) by transmitting the Mc into the corresponding ε (just as tabulated in Figure  2), and the predicted ∆p Eur values were obtained and compared with the corresponding predicted ∆p RSM . The comparison between the ∆p RSM and the ∆p Eur were evaluated by the validation metrics (Equations (2)-(5)), and the statistical results are tabulated in Table 6. As can be seen from Table 6, the index R 2 (0.9948) of the RSM model is higher than that of the Ergun model (0.9868), while the indexes MAE, MSE, and RMSE of the RSM model are, respectively, lower than the corresponding index of the Ergun model, which indicates that the pressure drop model developed by the RSM is better fitted than the Ergun model.

The Single Factor Analysis of the ∆p
In order to validate fitting performance of the ∆p RSM and ∆p Erg . The second-order polynomial RSM model was transmitted into three first-order quadratic models by substituting the independent variables shown in Table 3 into Equation (8) Based on the above Equations (6), (9)-(11), the variations of ∆p with Mc, v air , and L can be obtained and, respectively, depicted in the following Figures 7-9. As can be seen from Figure 7 that the ∆p is proportional to the Mc, and the variation range of the ∆p was ascertained to be 86. 38-197.57 Pa. It can be summarized from Figure 2 [41] and Ye et al. [10]. However, in this case, the developed ΔpRSM has a better performance than the Ergun model when the vair is above 0.2 m/s, wh be more suitable for guiding the actual process.   Mahmood et al. [41] and Ye et al. [10]. However, in this case, the developed has a better performance than the Ergun model when the vair is above 0.2 m be more suitable for guiding the actual process.   ΔpErg. In order to further investigate the fitting performance of the ΔpRSM and the ΔpErg in the above three single factor experiment, the statistical indexes mentioned in Section 2.3.3 were also calculated and the results are tabulated in the following Table 7.  As can be seen from Table 7 that the R 2 value of each single factor model based on the ΔpRSM is higher than that of the corresponding model based on the ΔpErg, while the other indexes (MSE, MAE, and RMSE) of each single factor model based on the ΔpRSM are lower than the corresponding index of the corresponding model based on the ΔpErg. Which indicates that the single factor model ΔpMc, Δpvair, and ΔpL based on the ΔpRSM has better fitting performance than that of the corresponding single factor model based on the ΔpErg. By combining the results obtained in Section 3.4.1, it can be concluded that the developed ventilation resistance model based on the RSM has a better predicting performance than the traditional Ergun model.

Conclusions
By considering the moisture content (Mc), air velocity (vair), and grain bed layer depth (L) as the key factors affecting the ventilation resistance of the air, one must regard the pressure drop (Δp) as the key index characterizing the ventilation resistance of the air. The present work investigated the influences of each factor and the interactions of the factors on the pressure drop; a second-order polynomial model was also obtained based on the RSM. The detailed conclusions depending on the main results are listed as follows: (1) The relations among the Mc, ε, and dg were revealed and the mathematical expressions were also given in the present work.   Figure 7, ∆p also increases with increase of v air , and the variation range of the ∆p was ascertained to be 38.44-352.92 Pa. Moreover, it can be seen from Figure 8 that the ∆p Erg with the v air of 0.1 m/s and 0.15 m/s are more closed to the corresponding experimental ∆p than the ∆p RSM , indicating that the Ergun model has a better prediction performance than the developed RSM model when the v air is under 0.2 m/s. The same conclusions have been made by Mahmood et al. [41] and Ye et al. [10]. However, in this case, the developed ∆p RSM model has a better performance than the Ergun model when the v air is above 0.2 m/s, which may be more suitable for guiding the actual process.
The variations of ∆p with L under the constant v air of 0.35 m/s and Mc of 18.9 %w.b. are depicted in Figure 9. It can be seen from the figure that the ∆p linearly increases with the increase of L, and the ∆p gets the minimum (27.68 Pa) and maximum value (276.83 Pa), respectively, at the L of 0.1 and 1m. On the other hand, it can be also seen from Figure 9 that, for the all levels of L, the fitting performance of the ∆p RSM is better than that of the ∆p Erg . In order to further investigate the fitting performance of the ∆p RSM and the ∆p Erg in the above three single factor experiment, the statistical indexes mentioned in Section 2.3.3 were also calculated and the results are tabulated in the following Table 7. As can be seen from Table 7 that the R 2 value of each single factor model based on the ∆p RSM is higher than that of the corresponding model based on the ∆p Erg , while the other indexes (MSE, MAE, and RMSE) of each single factor model based on the ∆p RSM are lower than the corresponding index of the corresponding model based on the ∆p Erg . Which indicates that the single factor model ∆p Mc , ∆p vair , and ∆p L based on the ∆p RSM has better fitting performance than that of the corresponding single factor model based on the ∆p Erg . By combining the results obtained in Section 3.4.1, it can be concluded that the developed ventilation resistance model based on the RSM has a better predicting performance than the traditional Ergun model.

Conclusions
By considering the moisture content (Mc), air velocity (v air ), and grain bed layer depth (L) as the key factors affecting the ventilation resistance of the air, one must regard the pressure drop (∆p) as the key index characterizing the ventilation resistance of the air. The present work investigated the influences of each factor and the interactions of the factors on the pressure drop; a second-order polynomial model was also obtained based on the RSM. The detailed conclusions depending on the main results are listed as follows: (1) The relations among the Mc, ε, and d g were revealed and the mathematical expressions were also given in the present work. The present work revealed the influences of the variables that can be directly measured in an engineering drying application on the pressure drop and established a prediction model of the pressure drop. Thus, the main results would be helpful for further optimizing the drying process. Further study is recommended to investigate the pressure drop variation characteristics by taking the air temperature as well as the paddy grain flow velocity into consideration.