CFD-Simulink Modeling of the Inﬂatable Solar Dryer for Drying Paddy Rice

: Small-scale farmers in developing Asian countries have minimal agricultural mechanisms available to them. In the Philippines, postharvest losses in rice production can reach about 36% in the drying process alone. Thus, the inﬂatable solar dryer (ISD) was developed through the collaboration of the University of Hohenheim, the International Rice Research Institute, and GrainPro Philippines Inc. Although the ISD was successfully tested with different agricultural products, further characterization of the ISD design is required for predicting the drying performance. To this end, the airﬂow behavior in the ISD was simulated using computational ﬂuid dynamics (CFD) via ANSYS Fluent. Moreover, a thermal model was developed in MATLAB/Simulink by taking into account heat transfer in the heating area and coupled heat and mass transfer within the drying area. Three batches of drying experiments were performed and airﬂow measurements were taken inside the dryer to validate the models. The MATLAB/Simulink model was further used to predict the drying performance under various weather conditions spanning 10 years. The simulated temperatures and moisture content in the ISD showed high accuracy (mean absolute percentage error (MAPE) < 10%) with the experimental data. The proposed dynamic model provides an efﬁcient computational tool that can be applied to predict the drying performance and to optimize the ISD design.


Introduction
Rice is one of the most produced and consumed staple cereals globally. In the Asia-Pacific region [1], 144 million rice farms, most of which are less than one hectare in area [2], produce over 90% of the world's production. Small-scale farmers in developing Asian countries have minimal agricultural mechanisms available to them. Almost all agricultural activities, including postharvest operations such as drying, are carried out manually. As a result, losses estimated during postharvest and processing stages across the rice value chain can range from 20% to 30% [3].
In the Republic of the Philippines, postharvest losses found in paddy rice can reach about 36% in the drying process alone [4,5]. Nowadays, most farmers still perform traditional sun drying practices. However, the farmers face several challenges when drying paddy rice due to high ambient relative humidity. This can be the case during unfavorable weather conditions or the rainy season, when there is lower solar radiation, or when sudden rainfall occurs. Under these conditions, the safe moisture content of 14% wet basis (w.b.) cannot be reached [6], and grains with high moisture contents are susceptible to attacks by microorganisms, insects, and pests [7]. Thus, decreasing postharvest losses is a realistic solution to increase the rice supply. Organizations and governments are supporting the development and the adaption of different technologies for small-and medium-scale farmers. Some examples are low-cost grain dryers [8], solar tunnel dryers [9], two-stage grain dryers [10,11], and flat-bed dryers [12,13].
The dissemination of flat-bed dryers for paddy rice has been promoted since the 1970s in the Republic of the Philippines, but this technology was not well received by farmers. Some of the reasons behind the lack of adoption were the higher operating and maintenance costs compared to traditional sun drying practices [14]. As an alternative, in the 1980s, the solar tunnel dryer was tested for paddy rice at the International Rice Research Institute (IRRI) in the Republic at the Philippines headquarters. Despite the drying performance being promising, mixing the grains along the tunnel dryer was challenging [9,15]. Several interventions followed that offered alternative drying practices. However, in most instances, energy usage and loading capability restricted the acceptance of these technologies. Therefore, under this scenario, the inflatable solar dryer (ISD) was developed through a collaboration of the University of Hohenheim, the IRRI, and GrainPro Philippines Inc., who manufactured the dryer (trade name "Solar Bubble Dryer™") [16].
The ISD is a further development of the solar tunnel dryer, which was developed by the University of Hohenheim in the 1980s. The ISD design characteristics are: ease of transport and installation; simple operation compared to mechanical dryers; low maintenance since it is made of plastic material. Currently, the ISD has been successfully tested on paddy rice [16,17], amaranth leaves [18], mushrooms [19], and maize [20]. However, information on airflow distribution or heat and mass transfer of the ISD is still lacking.
According to Ghaffari and Mehdipour [21], four methods are widely used to model and analyze solar dryer systems. The first method is thermodynamic modeling using semiempirical equations. The second method is to establish differential equations that are solved by using computational fluid dynamics (CFD). The third method is based on statistical methods that use a series of experiments to develop empirical equations. The fourth method is the reacting engineering approach based on the chemical reaction subject to the drying process.
Researchers such as Hossain, et al. [22] have mostly used the first method for modeling solar tunnel dryer designs. Their study provided an approach for the simulation of a solar tunnel dryer, considering heat transfer and mass transfer equations to predict temperatures, product temperature, and moisture content. The numerical solution was programmed in BASIC language. A similar method was reported by Janjai, et al. [23]. The study generates a comprehensive numerical solution programmed in Compaq Visual Fortran. The results showed a good agreement between the simulated and the experimental values of temperatures and moisture content.
Esper [24] combines the first and second modeling methods where the program ESATAN was applied to solve the thermal model around the solar tunnel dryer. The simulated crop temperature and moisture content were in good agreement with the experimental data.
Based on recent developments in computer simulation developments, this study applied the first and the second methods for modeling the ISD. Therefore, this study aims: (i) to graphically portray the spatial distribution of airflow across the ISD by applying a computational fluid dynamic (CFD) approach and (ii) to ascertain the heat transfer in the heating area and coupled heat and mass transfer within the drying area using transient thermal modeling. The collaborative models were validated with field experiments and were used to simulate the ISD's performance when drying paddy rice drying under various weather conditions. Furthermore, this computational approach will allow specific development of the ISD design to dry different agricultural commodities with a reduced number of experiments.

Dryer Description
The design and function of the ISD have already been described in a previous study [16]. Since then, however, the materials used to build the ISD have been reduced to a black reinforced polyvinyl chloride film (PVC, 0.52 mm) as the bottom layer, and a UV-stabilized transparent polyethylene film (PE, 150 µm) as the cover layer. A heavy-duty zipper connects the plastic films sewed along the edges of each plastic film, inflated by two axial flow ventilators of 72.6 W each, operated at 220 V (Runda Electronics CO., LTD, Shenzhen, China) ( Figure 1).

Dryer Description
The design and function of the ISD have already been described in a previous study [16]. Since then, however, the materials used to build the ISD have been reduced to a black reinforced polyvinyl chloride film (PVC, 0.52 mm) as the bottom layer, and a UV-stabilized transparent polyethylene film (PE, 150 µm) as the cover layer. A heavy-duty zipper connects the plastic films sewed along the edges of each plastic film, inflated by two axial flow ventilators of 72.6 W each, operated at 220 V (Runda Electronics CO., LTD, Shenzhen, China) (Figure 1).

Field Experiments
The experiments were performed at the International Rice Research Institute (IRRI) experimental station in Los Baños, Laguna, Republic of the Philippines (14°11'N, 121°15'E, 21 m a.s.l.). Mixed rice varieties were used for the three batches of experiments during the rainy season between October and November 2013. The dryer was assembled on the ground. The two plastic films were unzipped to load the paddy rice. Then, the zipper of the dryer was closed, and the blowers were turned on to provide pressure to inflate the unit. The special roller provided for mixing the grains was pulled underneath the PVC film every hour between 8:00 and 17:00. After mixing, the moisture content was monitored. Ten grams of paddy rice was collected at three positions, as shown in Figure 2a. The samples were stored in plastic bags and transported to the IRRI laboratories. The moisture content was determined by the gravimetric method at 103 ± 2 °C for 24 h [25]. The experiments were terminated when the rice reached a target moisture content of 14% wet basis (MCwb).

Instrumentation for the Field Experiments and Velocity Measurements
Solar radiation was measured using a pyranometer (CMP6, Kipp and Zonen, Delftechpark, Netherlands). Temperature and relative humidity capacity sensors (OM-EL-USB-2, Omega, Stamford, CT, USA) were placed inside the dryer at 1.5, 4, 14, 24 m from the inlet (Figure 2a). A further sensor was set to record the ambient temperature and relative humidity. The sensors were connected to a data logger (34970A, Agilent Technologies Inc., Loveland, CO, USA). The frequency of measurements was at five minutes intervals. A hot-wire anemometer (HHF42, Omega, Taipei, Taiwan) was used to measure the drying air velocity at the outlet position at one-hour intervals during the daytime.

Field Experiments
The experiments were performed at the International Rice Research Institute (IRRI) experimental station in Los Baños, Laguna, Republic of the Philippines (14 • 11 N, 121 • 15 E, 21 m a.s.l.). Mixed rice varieties were used for the three batches of experiments during the rainy season between October and November 2013. The dryer was assembled on the ground. The two plastic films were unzipped to load the paddy rice. Then, the zipper of the dryer was closed, and the blowers were turned on to provide pressure to inflate the unit. The special roller provided for mixing the grains was pulled underneath the PVC film every hour between 8:00 and 17:00. After mixing, the moisture content was monitored. Ten grams of paddy rice was collected at three positions, as shown in Figure 2a. The samples were stored in plastic bags and transported to the IRRI laboratories. The moisture content was determined by the gravimetric method at 103 ± 2 • C for 24 h [25]. The experiments were terminated when the rice reached a target moisture content of 14% wet basis (MC wb ).

Instrumentation for the Field Experiments and Velocity Measurements
Solar radiation was measured using a pyranometer (CMP6, Kipp and Zonen, Delftechpark, Netherlands). Temperature and relative humidity capacity sensors (OM-EL-USB-2, Omega, Stamford, CT, USA) were placed inside the dryer at 1.5, 4, 14, 24 m from the inlet (Figure 2a). A further sensor was set to record the ambient temperature and relative humidity. The sensors were connected to a data logger (34970A, Agilent Technologies Inc., Loveland, CO, USA). The frequency of measurements was at five minutes intervals. A hot-wire anemometer (HHF42, Omega, Taipei, Taiwan) was used to measure the drying air velocity at the outlet position at one-hour intervals during the daytime.
To monitor the air velocity profile over the paddy rice, additional grid measurements were taken with the hot wire anemometer at nine positions along the dryer length and the six positions over the dryer width. The anemometer was fixed at 0.08 m above the paddy layer, facing the main airflow direction (Figure 2b). Contour plots were generated by kriging (Surfer ® 19.1.189, Golden software, Golden, CO, USA).
To monitor the air velocity profile over the paddy rice, additional grid measurements were taken with the hot wire anemometer at nine positions along the dryer length and the six positions over the dryer width. The anemometer was fixed at 0.08 m above the paddy layer, facing the main airflow direction (Figure 2b). Contour plots were generated by kriging (Surfer ® 19.1.189, Golden software, Golden, CO, USA).
(a) (b) Figure 2. Position of measuring devices for (a) solar drying experiments (solar radiation, temperature, relative humidity, moisture content, and velocity) and (b) velocity measuring points across the dryer.

Governing Equations
The differential form for a three-dimensional ISD design was simulated from the mass conservation presented according to [26]: where is the air velocity vector and ϱair is the density of air. The Navier-Stokes equation is represented by the momentum conservation and it is written as [26]: where p is the static pressure, ̿ is the Reynolds stress tensor, g is the gravitational acceleration, and Sm is the source term for momentum defined as [27]: where vair is the velocity of the air, μair is the dynamic viscosity of air, ∈ is the porosity of paddy rice bulk and dk is the kernel diameter. The paddy rice porosity and kernel diameter were 0.57 [28] and 1.99 mm [29]. According to ANSYS [30], the first term represents the viscous resistance, and the second term represents the inertial resistance. The corresponding values for each term were 1.634 × 10 7 m −2 and 2047 m −1 . Symbols, terms, and units for each notation in all equations are described in the nomenclature list.

Governing Equations
The differential form for a three-dimensional ISD design was simulated from the mass conservation presented according to [26]: where → v is the air velocity vector and ρ air is the density of air. The Navier-Stokes equation is represented by the momentum conservation and it is written as [26]: where p is the static pressure, = τ is the Reynolds stress tensor, g is the gravitational acceleration, and S m is the source term for momentum defined as [27]: where v air is the velocity of the air, µ air is the dynamic viscosity of air, ∈ is the porosity of paddy rice bulk and d k is the kernel diameter. The paddy rice porosity and kernel diameter were 0.57 [28] and 1.99 mm [29]. According to ANSYS [30], the first term represents the viscous resistance, and the second term represents the inertial resistance. The corresponding values for each term were 1.634 × 10 7 m −2 and 2047 m −1 . Symbols, terms, and units for each notation in all equations are described in the nomenclature list.

Domain Description
The 3D geometry of the ISD used for the simulations has a total length of 25 m, where 3 m corresponds to the heating area and 22 m to the drying area ( Figure 3). The inlet has a diameter of 0.25 m for each blower. Likewise, the outlets have diameters of 0.25 m. The dryer's maximum height reached 1 m, and the paddy rice layer was set to a thickness of 0.04 m.

Domain Description
The 3D geometry of the ISD used for the simulations has a total length of 25 m, where 3 m corresponds to the heating area and 22 m to the drying area ( Figure 3). The inlet has a diameter of 0.25 m for each blower. Likewise, the outlets have diameters of 0.25 m. The dryer's maximum height reached 1 m, and the paddy rice layer was set to a thickness of 0.04 m. The fluid volume of the dryer was meshed using the ANSYS Workbench meshing tool. The mesh quality was assessed by ensuring that all elements' skewness and element quality fell within the optimal ranges [31]. A grid independence study was performed with four mesh refinements (mesh I-813490, mesh II-1072594, mesh III-1572592, and mesh IV-2498880). Finally, mesh III was selected as a good compromise between accuracy and computation time [32].

CFD Model Simulation
The simulation of the airflow distribution in the ISD was carried out using ANSYS Fluent 19.1 (ANSYS Workbench, ANSYS, Inc., Canonsburg, PA, USA) installed on a HP Z420 Workstation (Hewlett-Packard, Palo Alto, CA, USA) with an Intel Xeon CPU E5-1660 3.30 GHz processor and RAM 32 GB. The simulations were performed in a steady state and the k-є realizable turbulence model was used with enhanced wall treatment. The governing equations were computed by applying a pressure-based segregated solver and the default under-relaxation factors set as solution controls. The materials used were subdivided into the fluid (drying air across the ISD) and the solid (paddy rice). The properties of incompressible gas were applied. Based on the experimental data, the calculated inlet velocity was 5 m s −1 . The assumptions considered in the ISD simulation were nonadiabatic and no-slip walls, as well as no shrinkage or deformation of paddy rice kernels during drying. The paddy rice bulk was considered as a porous medium.

Mathematical Modeling of the ISD
The ISD is divided into segments of one m in the direction of the airflow. The math- The fluid volume of the dryer was meshed using the ANSYS Workbench meshing tool. The mesh quality was assessed by ensuring that all elements' skewness and element quality fell within the optimal ranges [31]. A grid independence study was performed with four mesh refinements (mesh I-813490, mesh II-1072594, mesh III-1572592, and mesh IV-2498880). Finally, mesh III was selected as a good compromise between accuracy and computation time [32].

CFD Model Simulation
The simulation of the airflow distribution in the ISD was carried out using ANSYS Fluent 19.1 (ANSYS Workbench, ANSYS, Inc., Canonsburg, PA, USA) installed on a HP Z420 Workstation (Hewlett-Packard, Palo Alto, CA, USA) with an Intel Xeon CPU E5-1660 3.30 GHz processor and RAM 32 GB. The simulations were performed in a steady state and the k-єrealizable turbulence model was used with enhanced wall treatment. The governing equations were computed by applying a pressure-based segregated solver and the default under-relaxation factors set as solution controls. The materials used were subdivided into the fluid (drying air across the ISD) and the solid (paddy rice). The properties of incompressible gas were applied. Based on the experimental data, the calculated inlet velocity was 5 m s −1 . The assumptions considered in the ISD simulation were nonadiabatic and no-slip walls, as well as no shrinkage or deformation of paddy rice kernels during drying. The paddy rice bulk was considered as a porous medium.

Mathematical Modeling of the ISD
The ISD is divided into segments of one m in the direction of the airflow. The mathematical model for the ISD was developed, taking into account heat transfer in the heating area and a coupled heat and mass transfer within the drying area. The drying of paddy Appl. Sci. 2021, 11, 3118 6 of 22 rice in the ISD was modeled following the approach of Hossain, et al. [22]. A schematic diagram of energy transfers inside the inflatable solar dryer is shown in Figure 4, and the following heat and mass balances are formulated. area and a coupled heat and mass transfer within the drying area. The drying of paddy rice in the ISD was modeled following the approach of Hossain, et al. [22]. A schematic diagram of energy transfers inside the inflatable solar dryer is shown in Figure 4, and the following heat and mass balances are formulated.
(a) (b)  Figure 4 shows the energy equilibration of a control volume in the heating area. The energy balance of the cover is given as: where Tcov is the temperature of the cover, Tamb is the temperature of the ambient conditions, Tair is the temperature of the air, Tsky is the temperature of the sky, Tabs is the temperature of the absorber, hw, hc,cov-air, hc,abs-air are the convection heat transfer coefficients between the cover and the air outside of the heating area due to wind, the cover and the air inside the heating area, and the surface of the heating area and the air inside, respectively. Additionally, hr,abs-air and hr,cov-air are the radiative transfer coefficients between the cover and the air outside of the heating area, and the absorber and the cover, respectively. αcov is the absorbance of the cover, τcov is the transmittance of the cover, ρabs is the reflectance of the absorber and It is the instantaneous solar radiation. Complementary information on the heat transfer coefficients and how to determine them are provided in Appendix A.

Energy Balance of the Absorber for the Heating area
The energy balance corresponding to the heating area relates to the change over time of the energy stored in the absorber: where Tasp is the temperature of the asphalt and Kabs-asp is the thermal conductance from the absorber to the asphalt.

Energy Balance of the Airflow for the Heating Area
The energy balance in the air stream in the collector is given as: where δair is the layer thickness of the air, vair is the air velocity, cair is the air specific heat of air, and ∂Tair/∂x is the temperature gradient in the direction of the airflow.

Mathematical Modeling of the ISD
The ISD is divided into segments of one m in the direction of the airflow. Th mathematical model for the ISD was developed, taking into account heat transfer in th heating area and a coupled heat and mass transfer within the drying area. The drying o paddy rice in the ISD was modeled following the approach of Hossain, et al. [22]. schematic diagram of energy transfers inside the inflatable solar dryer is shown in Figur 4, and the following heat and mass balances are formulated.  Figure 4 shows the energy equilibration of a control volume in the heating area. Th energy balance of the cover is given as: where Tcov is the temperature of the cover, Tamb is the temperature of the ambien conditions, Tair is the temperature of the air, Tsky is the temperature of the sky, Tabs is th temperature of the absorber, hw, hc,cov-air, hc,abs-air are the convection heat transfer coefficien between the cover and the air outside of the heating area due to wind, the cover and th air inside the heating area, and the surface of the heating area and the air insid respectively. Additionally, hr,abs-air and hr,cov-air are the radiative transfer coefficien between the cover and the air outside of the heating area, and the absorber and the cove respectively. αcov is the absorbance of the cover, τcov is the transmittance of the cover, ρa is the reflectance of the absorber and It is the instantaneous solar radiation Complementary information on the heat transfer coefficients and how to determine them are provided in Appendix A.

Energy Balance of the Absorber for the Heating area
The energy balance corresponding to the heating area relates to the change over tim of the energy stored in the absorber: where Tasp is the temperature of the asphalt and Kabs-asp is the thermal conductance from the absorber to the asphalt.

Energy Balance of the Airflow for the Heating Area
The energy balance in the air stream in the collector is given as: on of the airflow. The unt heat transfer in the ing area. The drying of Hossain, et al. [22]. A ryer is shown in Figure   - rature of the ambient re of the sky, Tabs is the eat transfer coefficients ind, the cover and the ea and the air inside, e transfer coefficients absorber and the cover, ttance of the cover, ρabs eous solar radiation. how to determine them to the change over time mal conductance from s: -convection, and →-conduction). Figure 4 shows the energy equilibration of a control volume in the heating area. The energy balance of the cover is given as:

Energy Balance of the Cover for the Heating Area
where T cov is the temperature of the cover, T amb is the temperature of the ambient conditions, T air is the temperature of the air, T sky is the temperature of the sky, T abs is the temperature of the absorber, h w , h c,cov-air , h c,abs-air are the convection heat transfer coefficients between the cover and the air outside of the heating area due to wind, the cover and the air inside the heating area, and the surface of the heating area and the air inside, respectively. Additionally, h r,abs-air and h r,cov-air are the radiative transfer coefficients between the cover and the air outside of the heating area, and the absorber and the cover, respectively. α cov is the absorbance of the cover, τ cov is the transmittance of the cover, ρ abs is the reflectance of the absorber and I t is the instantaneous solar radiation. Complementary information on the heat transfer coefficients and how to determine them are provided in Appendix A.

Energy Balance of the Absorber for the Heating area
The energy balance corresponding to the heating area relates to the change over time of the energy stored in the absorber: where T asp is the temperature of the asphalt and K abs-asp is the thermal conductance from the absorber to the asphalt.

Energy Balance of the Airflow for the Heating Area
The energy balance in the air stream in the collector is given as: ρ air ·δ air ·v air ·c air · ∂T air ∂x = h c,cov−air ·(T cov − T air ) + h c,abs−air ·(T abs − T air ) where δ air is the layer thickness of the air, v air is the air velocity, c air is the air specific heat of air, and ∂T air /∂x is the temperature gradient in the direction of the airflow.
The radiative and convective heat transfer coefficients were taken from Hossain, et al.'s work [22]. The convective heat transfer coefficient was computed as: where λ air is thermal conductivity of the air, D h is the hydraulic diameter, and Nu is the Nusselt number. The Reynolds number for a turbulent regime was calculated according to Cengel [33]. Nu = 0.03808·Re 0.8 ·Pr 1/3 (8) where Pr is the Prandtl number and Re is the Reynolds number, which is given according to [33]: where L c is the characteristic length of the geometry and υ air is the kinematic viscosity of air.

Energy Balance of the Cover in the Drying Area
According to Figure 4, the energy balance for the cover in the drying area results from the convective heat exchange between the cover and the heat transport by the drying air, which flows between the cover and the paddy rice in the drying area: where T pad is the temperature of paddy rice and ρ pad is the reflectance of paddy rice.
As found in a previous study [34], the reflectance is dependent on the zenith angle (θ):

Energy Balance of the Paddy Rice in the Drying Area
The change over time of the energy stored in the product results from the radiant energy on the product, the convective heat exchange between the product and the drying air, the radiant heat exchange with the cover, the environment, and the thermal losses to the ground as well as the energy required in the unit of time to evaporate the moisture. The energy balance is presented as: where ρ pad is the bulk density of paddy rice, δ pad is the thickness of the paddy rice bulk, c vap is the specific heat of water vapor, c liq is the specific heat of liquid, c pad is the specific heat of paddy rice, K pad-abs is the thermal conductance from paddy rice to absorber, M is the moisture content dry basis (d.b.), L pad is the latent heat of paddy rice [35] and α pad is the absorbance of paddy rice.

Energy Balance of the Airflow in the Drying Area
The energy balance for the moist air in the drying area is given as: ρ air ·δ air ·v air · c pad + c vap ·H · ∂T air ∂x = h c,cov−air ·(T cov − T air ) + h c,pad−air · T pad − T air (13) where H is the humidity ratio of air inside the dryer. The heat transfer coefficients are provided in Appendix A. As the ISD does not have an insulation material, the influence of the ground layers i (asphalt, soil, and absorber) is considered in the model following our previous research [34], in which the sun drying practice was modeled. The energy balance of the bottom layers under the ISD is driven by heat conduction. Figure 4 shows the thermal network across the bottom layers. The rate of thermal energy flow into the bottom layers is given as: where T i is the material temperature of the respective layer i, C i is the thermal capacity of the material, calculated by the specific heat capacity c i , and the thickness δ i of the layer i: The thermal conductance K i is given as: where λ i is the thermal conductivity of the material i (asphalt, soil, and absorber). The temperatures T top and T bottom of the layers are known boundary conditions between the layers and are given as:

Mass Balance
The moisture evaporating from the crop in the drying area is absorbed by the air. This causes a change in the enthalpy of the air, which can be described by the energy and mass balance. The rate of moisture transfer between the product and the surrounding air is expressed as: ρ air ·δ air ·v air · ∂H ∂x ·dx·dt = −ρ pad ·dx· ∂M ∂t ·dt (19) where ∂H/∂x is the humidity ratio and ∂M/∂t is the drying rate gradient in the direction of the airflow. The drying rate is calculated as: where M e is the equilibrium moisture content and k is the drying constant. The calculation of M e was obtained from the Chung-Pfost equation using the coefficients given by Iguaz and Vírseda [36]: where ϕ is the relative humidity of the drying air. The drying constant was obtained from paddy rice drying experiments conducted under controlled conditions. This experiment was necessary as no other paper on the subject had been reported. The drying experiments were carried out using the overflow chamber of a high-precision laboratory dryer from the Institute of Agricultural Engineering, University of Hohenheim in Stuttgart (Germany). The experimental system has been previously described in detail by Udomkun, et al. [37]. The drying experiments were conducted at air temperatures of 30, 40, 50, and 60 • C and humidity ratio of 0.020 kg kg −1 .
The paddy rice was filled in a tray of 0.30 × 0.30 × 0.04 m having a layer thickness of 0.04 m. The paddy rice was dried to the equilibrium moisture content. The experiments were performed per triplicate for each drying condition.
The drying equation was built following three steps. The first step corresponds to the computation of the moisture ratio from the drying experiments; the second step was to obtain the drying coefficient (k) per experiment attributed to the Newton model. The third step was to fit the drying coefficient of the different experimental conditions by using an Arrhenius type equation: The goodness of fit for Equation (22) resulted in a coefficient of determination (R 2 ) of 0.87 and a root mean square error (RMSE) of 0.10 h −1 .

Solution Procedure
The partial differential equations generated for the heating area and the drying area describe heat and heat-mass transfer over small time increments. The numerical integration with respect to time and position can be solved using the finite difference method [22]. The transformation from the physical domain to the computational domain was obtained by creating a numerical simulation model of algebraic equations in one dimension. The finite volume differences of formulated heat and mass balances were replaced by the derivatives in the differential equations.
The drying air temperature given in Equation (6) was calculated over the finite distance (∆x) to give the change in air temperature in the heating area (∆T air ) as: ∂T air ∂x = h c,cov−air ·(T cov − T air ) + h c,abs−air ·(T abs − T air ) ρ air ·δ air ·v air ·c air Introducing auxiliary terms A 1 and A 2 and noting h c for h c,cov-air and h c,abs-air : ρ air ·δ air ·v air ·c air A 2 = T cov + T air 2 (24) leads to: The solution for the energy balance of the air in the drying area given in Equation (13) is solved as: where auxiliary terms A 3 and A 4 are: ρ air ·δ air ·v air · c air + c vap ·H A 4 = T cov + T air ρ air ·δ air ·v air · c air + c vap ·H (27) leading to: The temperature of paddy rice generated in Equation (12) can be solved as: with auxiliary terms A 5 and A 6 : +ρ pad ·δ pad ·L pad · ∂M ∂t ρ pad ·δ pad ·(c pad +c liq ·M) + h c,pad−air ·T air +h r,pad−sky ·T sky +h r,pad−cov ·T c +K pad−abs ·T abs ρ pad ·δ pad ·(c pad +c liq ·M) ∂t +h c,pad−air +h r,pad−sky +h r,pad−cov +K pad−abs ρ pad ·δ pad ·(c pad +c liq ·M) (30) leading to: The rate of change of moisture content inside the dryer is represented in the finite difference form according to Bala [35]: The change of humidity solution can be calculated as: Table 1 presents the material properties of the ISD, the paddy rice, and the ground layers used for the MATLAB/Simulink ® model. MATLAB/Simulink ® 2019a (MathWorks Inc., Natick, MA, USA) was used to compute the temperatures and the moisture content behavior of the paddy rice drying inside the ISD. The lumped equations are a set of partial differential equations that are converted to coupled ordinary equations to solve using Simulink graphical interface [41]. The flowchart of the solution is illustrated in Figure 5. The coefficient of determination (R 2 ), the root mean square error (RMSE), and mean absolute percentage error (MAPE) were used to assess the accuracy of prediction models.
Specific heat csoi 870.0 J kg −1 K −1 Thermal conductivity λsoi 2.50 W m −1 K −1 MATLAB/Simulink ® 2019a (MathWorks Inc., Natick, MA, USA) was used to compute the temperatures and the moisture content behavior of the paddy rice drying inside the ISD. The lumped equations are a set of partial differential equations that are converted to coupled ordinary equations to solve using Simulink graphical interface [41]. The flowchart of the solution is illustrated in Figure 5. The coefficient of determination (R 2 ), the root mean square error (RMSE), and mean absolute percentage error (MAPE) were used to assess the accuracy of prediction models.

Model Implementation
The weather data were collected from the weather station at the experimental station of the International Rice Research Institute in Los Baños, Laguna, Republic of the Philippines at 14°11'N, 121°15'E. The annual weather records from 2008 to 2017 were collected from the weather station every 15 min. Ambient temperature, relative humidity, and daily solar radiation data were averaged (see Figure 6). The weather data were used for the model implementation corresponding to the peak harvest season for the dry season (April and May) and the rainy season (October and November). Table 2 shows the initial values for the ISD model in MATLAB/Simulink.

Model Implementation
The weather data were collected from the weather station at the experimental station of the International Rice Research Institute in Los Baños, Laguna, Republic of the Philippines at 14 • 11 N, 121 • 15 E. The annual weather records from 2008 to 2017 were collected from the weather station every 15 min. Ambient temperature, relative humidity, and daily solar radiation data were averaged (see Figure 6). The weather data were used for the model implementation corresponding to the peak harvest season for the dry season (April and May) and the rainy season (October and November). Table 2 shows the initial values for the ISD model in MATLAB/Simulink.

Simulation of the Airflow Distribution
The simulation of the airflow pattern in the inflatable solar dryer (ISD) is shown in Figure 7, where the airflow direction travels from left to right. Figure 7a shows the vortex caused by the fans located at the inlet. The vortex covers the first meters of the ISD, where the air velocity is the highest (Figure 7b). Furthermore, the airflow in the drying area forms a transition zone, where the airflow presents an even distribution over the paddy rice. However, backflow was observed one meter before exiting the dryer. The turbulence model k-є predicted the presence of the air turbulence in the area located at 3 to 4 m from the fans well, creating a recirculation pattern.

Simulation of the Airflow Distribution
The simulation of the airflow pattern in the inflatable solar dryer (ISD) is shown in Figure 7, where the airflow direction travels from left to right. Figure 7a shows the vortex caused by the fans located at the inlet. The vortex covers the first meters of the ISD, where the air velocity is the highest (Figure 7b). Furthermore, the airflow in the drying area forms a transition zone, where the airflow presents an even distribution over the paddy rice. However, backflow was observed one meter before exiting the dryer. The turbulence model k-є predicted the presence of the air turbulence in the area located at 3 to 4 m from the fans well, creating a recirculation pattern.

Validation of the CFD Model
A contour plot was used to compare the CFD simulation and the experimental values and to analyze the airflow distribution over the paddy rice ( Figure 8). The simulated results were consistent with Figure 7. Thus, the airflow velocity was higher at the left side. It can be seen that the air velocity profile at 0.08 m decreased in the heating area from 5.0 to 0.5 m s −1 , and the velocity in the drying area was around 0.1 m s −1 . Figure 9 shows the computational and experimental velocities over the ISD's length. The results indicate that air velocity decreased along the dryer length and that air velocity magnitude in the x-direction at 0.08 m above the paddy rice verifies the computational results with the experimental data. The comparison with the computed and experimental data revealed a good correlation, with a R 2 of 0.78 and a RMSE of 0.72 m s −1 .

Validation of the CFD Model
A contour plot was used to compare the CFD simulation and the experimental values and to analyze the airflow distribution over the paddy rice (Figure 8). The simulated results were consistent with Figure 7. Thus, the airflow velocity was higher at the left side. It can be seen that the air velocity profile at 0.08 m decreased in the heating area from 5.0 to 0.5 m s −1 , and the velocity in the drying area was around 0.1 m s −1 .        Figure 10 presents the solar radiation, the ambient temperature, and the relative humidity during the three experimental batches of solar drying of paddy rice in the ISD. The solar radiation during the rainy season shows higher variations due to clouds and sudden rainfall. When the sun rises, the ambient relative humidity declines gradually and keeps to a minimum value of 55% at around noon. This corresponds to the maximum temperature of 36 • C. The ambient temperature during nighttime reaches about 22 • C, while the relative humidity increases to 88%. midity during the three experimental batches of solar drying of paddy rice in the ISD solar radiation during the rainy season shows higher variations due to clouds and sud rainfall. When the sun rises, the ambient relative humidity declines gradually and k to a minimum value of 55% at around noon. This corresponds to the maximum temp ture of 36 °C. The ambient temperature during nighttime reaches about 22 °C, while relative humidity increases to 88%. Figure 10. The variation of ambient temperature, relative humidity, and solar radiation for the three experimental batches during rainy season.
The simulated temperature was validated against the experimental paddy rice d as shown in Figure 11. It can be observed that the temperature increase in the drying was well reproduced throughout the drying process. During the daytime, the ave Figure 10. The variation of ambient temperature, relative humidity, and solar radiation for the three experimental batches during rainy season.
The simulated temperature was validated against the experimental paddy rice data, as shown in Figure 11. It can be observed that the temperature increase in the drying area was well reproduced throughout the drying process. During the daytime, the average drying air temperatures varied from 5 to 25 • C above the ambient temperature. The highest experimental and simulated drying temperatures were 63.7 and 60.2 • C at 13:05 and 12:00, respectively. The high fluctuation of weather conditions such as cloudiness or sudden showers influenced the temperature in the ISD. This phenomenon can be also observed from the simulated temperature. Therefore, the simulation was in good agreement with experimental data. drying air temperatures varied from 5 to 25 °C above the ambient temperature. The highest experimental and simulated drying temperatures were 63.7 and 60.2 °C at 13:05 and 12:00, respectively. The high fluctuation of weather conditions such as cloudiness or sudden showers influenced the temperature in the ISD. This phenomenon can be also observed from the simulated temperature. Therefore, the simulation was in good agreement with experimental data.
(a) (b) Figure 11. Simulated (a) and experimental (b) temperature profile over time and along the length of the ISD of the 3 experimental batches. Figure 12 shows the comparison between the experimental and predicted moisture contents during solar drying of paddy rice in the ISD. Rewetting during night or rainy hours was not observed as the fans were operating continuously. The target moisture content was reached with the ISD between 27 and 29.5 h according to the experimental values. From the experimental data, it can be seen that increased fluctuation of weather conditions strongly influences the drying behavior of paddy rice. Nevertheless, the simulation model could effectively predict the moisture content changes during drying, particularly under good weather conditions, as shown in batch 3.  Figure 12 shows the comparison between the experimental and predicted moisture contents during solar drying of paddy rice in the ISD. Rewetting during night or rainy hours was not observed as the fans were operating continuously. The target moisture content was reached with the ISD between 27 and 29.5 h according to the experimental values. From the experimental data, it can be seen that increased fluctuation of weather conditions strongly influences the drying behavior of paddy rice. Nevertheless, the simulation model could effectively predict the moisture content changes during drying, particularly under good weather conditions, as shown in batch 3.

Accuracy of the Model
The goodness of fit of the simulation model is presented in Table 3. High accuracy can be observed from the temperature in the heating area and the outlet, as well as the moisture content of the paddy rice during drying. The highest R 2 can be observed from batch 3 where the weather conditions were stable, especially the first day compared to the other batches. For the three batches, the RMSE for temperature was less than 6.6 • C and 1.4% for moisture content. In general, the MAPE for temperatures is below 6.9%, while the MAPE of moisture content is below 7.4%. As the MAPE of the simulation model was still below 10%, high prediction accuracy was therefore proven [42].

Accuracy of the Model
The goodness of fit of the simulation model is presented in Table 3. High accuracy can be observed from the temperature in the heating area and the outlet, as well as the moisture content of the paddy rice during drying. The highest R 2 can be observed from batch 3 where the weather conditions were stable, especially the first day compared to the other batches. For the three batches, the RMSE for temperature was less than 6.6 °C and 1.4% for moisture content. In general, the MAPE for temperatures is below 6.9%, while the MAPE of moisture content is below 7.4%. As the MAPE of the simulation model was still below 10%, high prediction accuracy was therefore proven [42].    Figure 13a shows the simulated moisture content during the relevant months in the dry season and the rainy season for the 10-year period (2008-2017). Figure 13b summarizes the average drying time for the dry season and the rainy season, where 14% MC wb was achieved. The prediction shows that during the dry season, the target moisture content can be reached within 10 h, while during the rainy season, the drying time increases, and the drying time may be higher than 30 h. Figure 13a shows the simulated moisture content during the relevant months in the dry season and the rainy season for the 10-year period (2008-2017). Figure 13b summarizes the average drying time for the dry season and the rainy season, where 14% MCwb was achieved. The prediction shows that during the dry season, the target moisture content can be reached within 10 h, while during the rainy season, the drying time increases, and the drying time may be higher than 30 h.

Discussion
The current study investigated the airflow patterns in the ISD. The results obtained with the CFD simulation model helped to visualize the airflow distribution (Figure 7). A velocity reduction is expected due to an abrupt expansion from the fan to the dryer. These results are likely to be related to a turbulent regime [43] and generated vortex in the heating area, as seen in Figure 7a. A similar flow structure has been shown in another study [44].
Previous studies evaluating the airflow patterns in solar dryers have been conducted using a CFD approach [21,45,46], where the study of the turbulence effect was proportionally related to the rate of heat transfer. This effect was also observed in the ISD where the presence of higher velocity in the heating area and lower velocities over the paddy rice spread along the drying area shows a temperature rise across the dryer (Figure 11b).
The simulated results from the MATLAB/Simulink model are consistent with the experimental data ( Figure 11). A similar behavior was also observed from drying various agricultural products using the Hohenheim solar tunnel dryer type [23,[47][48][49][50].
Several studies have reported different simulation models of forced convection solar drying systems, where a thin layer of drying materials was commonly assumed and described [22,23,51]. In this study, we have used the same approach and have further developed the drying kinetics of paddy rice. For example, we show that the prediction of moisture content during drying is more accurate compared to the prediction of temperatures (Table 3). Based on the results, the drying process in the ISD can be computationally evaluated with a good agreement to the empirical measurements.
The results of the application of the model over the dry and the rainy seasons further support the results of the previous ISD experiments [16], where drying during the rainy season required a long period of time, while in the dry season it was accomplished within 10 h of the drying process. Further research should be undertaken with ANSYS fluent and

Discussion
The current study investigated the airflow patterns in the ISD. The results obtained with the CFD simulation model helped to visualize the airflow distribution (Figure 7). A velocity reduction is expected due to an abrupt expansion from the fan to the dryer. These results are likely to be related to a turbulent regime [43] and generated vortex in the heating area, as seen in Figure 7a. A similar flow structure has been shown in another study [44].
Previous studies evaluating the airflow patterns in solar dryers have been conducted using a CFD approach [21,45,46], where the study of the turbulence effect was proportionally related to the rate of heat transfer. This effect was also observed in the ISD where the presence of higher velocity in the heating area and lower velocities over the paddy rice spread along the drying area shows a temperature rise across the dryer (Figure 11b).
The simulated results from the MATLAB/Simulink model are consistent with the experimental data ( Figure 11). A similar behavior was also observed from drying various agricultural products using the Hohenheim solar tunnel dryer type [23,[47][48][49][50].
Several studies have reported different simulation models of forced convection solar drying systems, where a thin layer of drying materials was commonly assumed and described [22,23,51]. In this study, we have used the same approach and have further developed the drying kinetics of paddy rice. For example, we show that the prediction of moisture content during drying is more accurate compared to the prediction of temperatures (Table 3). Based on the results, the drying process in the ISD can be computationally evaluated with a good agreement to the empirical measurements.
The results of the application of the model over the dry and the rainy seasons further support the results of the previous ISD experiments [16], where drying during the rainy season required a long period of time, while in the dry season it was accomplished within 10 h of the drying process. Further research should be undertaken with ANSYS fluent and MATLAB/Simulink to investigate the application of the collaborative model approach with different agricultural commodities and weather conditions.

Conclusions
This study presents CFD-Simulink models used to simulate and predict the ISD's performance when drying paddy rice. The airflow patterns were successfully investigated across the ISD using a CFD approach via ANSYS Fluent. A good agreement was revealed between the CFD simulation and experimental velocity measurements with an R 2 of 0.78 and RMSE of 0.72 m s −1 . A variation of air velocity from 5 to 0.5 m s −1 and 0.5 to 0.1 m s −1 was obtained at the heating and drying area, respectively. A vortex was observed from 3 to 4 m in the heating area caused by the air conveyance of the axial fans. Moreover, we successfully modeled the drying process of paddy rice using the MATLAB/Simulink. A good agreement between the experimental and the simulated moisture contents of paddy rice (R 2 > 0.82, RMSE < 1.4 • C, MAPE < 7.4%) and temperatures (R 2 > 0.73, RMSE < 6.6%, MAPE < 6.9%) was achieved. The validated model allowed an accurate assessment of the drying behavior of paddy rice over the rainy and dry seasons to be made.
The CFD-Simulink models proved to be an efficient and effective means of simulating and assessing the ISD's performance under different weather conditions. The computational approach used in this study shows great potential as a decision support tool to improve the ISD design. Further investigations should include the mixing effect of paddy rice during drying and the ratio of heating and drying area.