Next Article in Journal
Microstructural Changes Influencing the Magnetoresistive Behavior of Bulk Nanocrystalline Materials
Next Article in Special Issue
GPU-Enabled Shadow Casting for Solar Potential Estimation in Large Urban Areas. Application to the Solar Cadaster of Greater Geneva
Previous Article in Journal
Influence of Medium Frequency Light/Dark Cycles on the Cultivation of Auxenochlorella pyrenoidosa
Previous Article in Special Issue
The Buying Time Argument within the Solar Radiation Management Discourse
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Surface Reflection (Albedo) in Simulating the Sun Drying of Paddy Rice

by
Ana Salvatierra-Rojas
*,
Victor Torres-Toledo
and
Joachim Müller
Institute of Agricultural Engineering, Tropics and Subtropics Group (440e), University of Hohenheim, 70599 Stuttgart, Germany
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(15), 5092; https://doi.org/10.3390/app10155092
Submission received: 30 June 2020 / Revised: 17 July 2020 / Accepted: 22 July 2020 / Published: 24 July 2020
(This article belongs to the Special Issue Solar Radiation: Measurements and Modelling, Effects and Applications)

Abstract

:
The sun drying of agricultural products is a complicated process involving heat transfer, mass transfer, and variable weather conditions. Surface reflection (albedo), a crop’s radiative property, plays an essential role in energy balance, and understanding its contribution can improve the thermal analysis. In this study, field experiments were conducted in the Philippines to explore the influence of surface albedo on the sun drying of paddy rice. First, we implemented energy and mass balance equations in a transient model with the surroundings using a graphical programming language in Matlab/Simulink®. Second, we identified the influence of albedo on the sun drying model by using a sensitivity analysis. Third, we investigated the relationship of paddy rice albedo and the solar zenith angle. Lastly, we integrated the albedo function into the sun drying model. The simulation outputs were validated with field experiments. A better estimation of the measured exit temperature and instantaneous mass were obtained when a variable albedo was applied. This study makes clear that introducing a variable albedo has a positive impact on model improvement. This information is important for application in solar drying technologies, so that the drying process can be better assessed.

1. Introduction

The sun drying of paddy rice is a common practice in the Philippines and many other countries in the tropics and subtropics. Immediately after harvesting/threshing, paddy rice is spread over mats positioned alongside roads or on other paved ground when weather conditions are favorable. To ensure and optimum milling process, it is recommended to dry the paddy rice to a moisture content (MC) of around 14% wet basis [1,2]. During the dry season, the grain temperature can reach up to 55 °C due to the direct exposure of the paddy rice to high solar radiation. Under these conditions, the paddy rice can reach a very low MC (<10% wet basis) [3,4], which may result in fissured grains after milling, thus decreasing its market value. However, during the rainy season, the lower radiation and extended rainy periods may lead to a prolonged drying process. This delay in drying increases the risk of paddy rice deterioration. This deterioration of paddy rice can be observed in the milled rice.
The sun drying of paddy rice is a complex process, which involves the removal of water through direct exposure of the crop to solar radiation. Heat is transferred from the surroundings and from the sun to the exposed crop surface. Raising the crop temperature helps to diffuse water vapor from the interior of the kernel to the surface. Hence, mathematical models have been developed to understand and predict the sun drying process of different agricultural products. Jain and Tiwari [5] proposed the calculation of the convective heat transfer coefficient (hc), followed by the generation of a mathematical model for the estimation of the crop temperature (Tcrop) and moisture evaporated (mev). To simplify the model’s solution, a steady state condition is commonly used, and indeed several authors follow this methodology [6,7,8,9]. Even though meteorological conditions have a strong influence on the solar radiation incident on the crop surface, constant radiative properties, such as the reflectivity and absorptivity of the crop, are commonly assumed during the modeling of the sun drying of medicinal plants, cereals, vegetables, and fruits.
A fraction of the incident solar radiation is partly absorbed and scattered by the atmosphere and the earth’s different surfaces. Here, the surface reflection (albedo) is the ratio of the reflected solar radiation to the incoming solar radiation [10,11]. Albedo plays a crucial role in the fraction of solar energy absorbed and is a fundamental parameter applied in different research fields, including agriculture, climatology, forestry, building energy, and photovoltaics [12,13,14,15]. Previous studies have explored the parametrization of albedo, and most of them reported a close relationship with the solar zenith angle [16,17,18,19]. The parametrization of albedo has commonly been used for the improvement of the accuracy of models such as land, weather, and climate forecast models [20].
In rice production, several studies have assessed albedo measurements during the whole cropping period, and some even monitored it up until the harvesting stage. The values ranged from 0.09 to 0.25 [10,21,22,23,24,25,26,27,28,29,30]. This information is commonly used for the estimation of surface energy, water budget, and potential cropping yield. However, limited research on albedo in the postharvest processes was found, particularly in the sun drying of paddy rice, where solar energy plays a critical role in the success of drying. Belessiotis and Delyannis [31] suggest albedo values depending on the surface color of the product: dark-colored materials have an albedo value of 0.10, while gray, red, and green materials have a value about 0.30, and light-colored materials have a surface value of about 0.50. For brown rice, Arnize [21] and Bala [32] proposed an albedo value of 0.16.
Three general procedures for measuring albedo can be found in the literature: Laboratory measurements, remotely sensed measurements, and field measurements. First, the laboratory measurements comprise the use of spectral radiometers and solar simulators to measure the spectral reflectivity of small samples; more extensive measurements can be carried out remotely by fixing the spectral radiometer on aircraft or satellites [33]. The second measurement is the remote sensing technique, which provides spectral data that can generate essential information. However, this technique might not be practical for field applications as it requires a careful estimation of atmospheric attenuation and has limited accuracy for small test surfaces in heterogeneous field environments [34].
The third procedure, taking field measurements of albedo, typically involves installing pyranometers over the test area at a height of less than half a meter, according to the ASTM E1918-06 standard [35]. This standard requires a target area of a minimum of 4 m in diameter, and the measurement should be performed in an open area [36]. Sailor et al. [37] proposed the use of a circular shield limiting the target area to 1 m and with a pyranometer installed at a height of 0.18 m. The shadows produced by the sensor and the surroundings were estimated by calculating the view factors of the mentioned surfaces. This last approach was applied in the present study.
In this article, an optimized model approach is presented for the sun drying of paddy rice to predict crop temperature and mass evaporation. The objectives of this study were to (i) develop a dynamic model to simulate the sun drying of paddy rice, (ii) determine the relationship of albedo and the solar zenith angle, and (iii) evaluate the accuracy of the model at constant and variable albedos in different environmental conditions in the Philippines.

2. Materials and Methods

2.1. Model Description

The sun drying of paddy rice was modeled based on the methodology proposed by Jain and Tiwari [5]. Figure 1 schematically shows the sun drying principle as a one-dimensional model of the thermal network with seven nodes. The layers from top to bottom are paddy rice, polyvinyl chloride (PVC) tarpaulin, asphalt, and soil. The energy and mass balance equations were formulated as ordinary differential equations (ODEs).

2.1.1. External Node

According to the thermal network shown in Figure 1, energy balance at the external node can be expressed as:
m t · c p a · d T p a d t = ( 1 α p a ) · G u p · A h r c . p a a e · ( T p a T a e ) · A Q e v · A K p a · ( T p a T b o t t o m ) · A
where mt (kg) is the instantaneous mass of the paddy, αpa (-) is the albedo of paddy rice, Gup (W m−2) is the total solar radiation, A (m−2) is the area of the crop exposed to the sun, hrc,pa-ae (W m−2 K−1) is the combined radiative and convective heat transfer coefficient from the paddy rice surface to the air, Tae (°C) is the temperature of the air exiting the paddy rice surface—henceforth called exit temperature—Tpa (°C) is the temperature of the paddy rice, and Qev (kJ m−2 s−1) is the rate of heat utilized to evaporate moisture. Qev is calculated as [38]:
Q e v = 0.016 · h c , p a a e   · [ P T p a r h a m · P T a e ]
where hc,pa-ae (W m−2 K−1) is the convective heat transfer coefficient from the paddy rice surface to the air, rham (-) is the ambient relative humidity, and PTpa and PTae (N m−2) are the saturated vapor pressure at the respective temperatures. The values of PT were extracted from the temperature range from 25 to 55 °C [39], which corresponds to the operating condition of the sun drying experiments. A 1D lookup table block built in Simulink was used to retrieve the saturated vapor pressure at the respective temperatures [40].
The energy balance of moist air above the paddy rice is calculated as:
h r , p a a e · ( T p a T a e ) · A + Q e v · A = h w · ( T a e T a m ) · A
where Tam (°C) is the ambient temperature and hw (W m−2 K1) is the convective heat transfer coefficient for wind. hw is computed according to Duffie and Beckmann [41]:
h w = 5.7 + 3.8 · V w
where Vw (m s−1) is the wind velocity.
The radiative heat transfer coefficient hr,pa-ae is written as:
h r , p a a e = ε · σ · ( T p a 2 + T a e 2   ) · ( T p a + T a e )
where the emissivity ε of a solid body is equal to the absorbance (1-α) measured at the same wavelengths according to Kirchhoff’s law [32,42] and σ is the Boltzmann constant (5.6696 × 10−8 W m−2 K−4).
The calculation of the mass of paddy rice mt during drying is done by subtracting the evaporated water from the initial mass m0. The evaporated mass of water mev is obtained as:
m e v = Q e v L T p a     · A · t
where t (h) is the time of the sampling interval and L (kJ kg−1) is the latent heat of vaporization, which is calculated as [43]:
L = 2502.535 2.386 · T p a

2.1.2. Intermediate Nodes

The energy balance of the intermediate nodes is driven by heat conduction. Figure 1 shows the thermal network across the bottom layers, and the rate of thermal energy flow into the bottom layers is given as:
C i · d T i d t = K i · ( T t o p T i ) + K i · ( T b o t t o m T i )
where Ti corresponds to the material temperature of the respective node and Ci is the lumped thermal capacity of the material around the node, which can be calculated by knowing the specific heat capacity cpi (kJ m−3 K−1), and the thickness zi (m) of the material (9).
C i = c p i · z i
The thermal conductance Ki is given as:
K i = λ i 1 2 · z i
where λi (kJ m−2 K−1) is the thermal conductivity of the material i (soil, asphalt, and PVC tarpaulin).
The temperatures Ttop and Tbottom of the layers are known boundary conditions between the layers and are given in (11) and (12), respectively:
T t o p = K i · T i + K i + 1 · T i + 1 K i + K i + 1
T b o t t o m = K i · T i + K i 1 · T i 1 K i + K i 1
Table 1 presents the default input parameter values for the model that were found in the literature and average values that were used in the model.
Equations (1)–(12) were solved numerically using a variable order solver based on the numerical differentiation formulas in the ordinary differential equations (ODE15s solver Matlab/Simulink® Version 9.6 [56]). Figure 2 illustrates the flow chart procedure of the simulation.

2.2. Measurement of Albedo

To analyze the actual values for the albedo of paddy rice, experiments were performed at the International Rice Research Institute (IRRI), located in Los Baños, Laguna province, Republic of the Philippines (14°9′55″ N, 121°15′3″ E, 21 m a.s.l.). The albedo of paddy rice was measured according to the protocol developed by Sailor et al. [37]. In total, three sets of experiments were conducted, namely on 29 and 30 November 2012, 1 December 2012, and on 14, 15, and 25 April 2013, mostly under clear sky conditions. According to Sailor et al. [37], the measurements should be performed around noon. The measurements carried out in this study were between 10:30 and 15:00. Figure 3 shows the setting of the sensors for the measurement of albedo.
Albedo was determined by measuring total solar radiation (direct and diffuse solar radiation) with an upward-facing pyranometer, while measuring reflected solar radiation was done with a downward-facing pyranometer (CMP11 and CMP6, Kipp and Zonen, Delftechpark, Delft, the Netherlands) from a vertical distance of 0.18 m above the paddy surface. The properties of the two pyranometers are presented in Table 2. Separate outputs for each sensor were automatically recorded using a data logger (34970A, Agilent Technologies Inc., Loveland, CO, USA). Data were logged at 1-min intervals.
The advantages of using the circular shield (Figure 3) are the reduction of the area subjected to the measurement from the minimum 4 m (standard [35]) to 1 m, the performance of the experiment without the restriction of a larger testing space, and the avoidance of the solar reflection of the surrounding surfaces on the downward pyranometer. Therefore, this setup allowed the reflected radiation to be measured over the course of the day by excluding the areas shaded by the pyranometer, its support, and the circular shield. Albedo was calculated as:
α p a = G d o w n G u p ( F s e n s o r b + γ d i f f · F s e n s o r f )   · α c s F s e n s o r b + γ d i f f · ( F s e n s o r b + F s e n s o r b + F s e n s o r b )
where αpa (-) is the albedo of paddy rice, αcs (-) is the albedo of the circular shield, Ddown (W m−2) is the radiation measured by the inverted pyranometer, Gup (W m−2) is the radiation measured by the upright pyranometer, γdiff (-) is the diffuse fraction, and the view factors F between the sensor and the six surfaces are: (a) unshaded paddy sample F(sensor-a), (b) unshaded circular shield F(sensor-b), (c) pyranometer shadow F(sensor-c), (d) support rod shadow F(sensor-d), (e) crescent shade of the circular shield F(sensor-e), (f) and shaded part of the circular shield F(sensor-f), these surfaces are shown in Figure 4.
The clearness index k t was used for the calculation of γdiff, according to Berrizbeitia et al. [57] and is valid for latitudes between 8.5° and 28° N.
γ d i f f =                   { 0.98   ·   k t 1.024 + 0.47 ·   k t 3.62 ·   k t 2 + 2 ·   k t 3   0.16               k t < 0.2 0.2 < k t < 0.77 k t > 0.77
where kt is the ratio of total solar radiation Gup and extraterrestrial radiation Gext. Gext is calculated as:
G e x t = 1353 · [ 1 + 0.033 ·   c o s ( 360 ·   N 365 ) ] ·   cos θ
where N is the Julian day of the year and θ is the solar zenith angle, which is represented in (16):
θ = arccos ( sin ϕ · sin δ + cos δ · cos ϕ · cos ω )
where δ is the declination angle, ϕ is the latitude, and ω is the hour angle.
For modeling the daytime-dependent course of albedo, an approach of Zheng et al. [58] was chosen:
α = a + b · e x p c · θ
where a, b, and c denote statistically estimated coefficients.
Finally, the albedo was integrated as a variable into the model for the sun drying of paddy rice.

2.3. Drying Experiments

To parametrize the model, sun drying experiments were performed on the premises of IRRI. Three batches of experiments were performed during the rainy season (October 2011 and November 2012) and dry season (May 2012 and April 2013). The paddy rice was spread out on a black PVC tarpaulin with a bulk height (zpa) of 40 mm. The grains were mixed manually by using a wooden rake at 1-h intervals during daylight.
To investigate the drying performance, omega sensors (OM-EL-USB-2, Omega, Stamford, CT, United States of America) were used to monitor temperature and relative humidity over the paddy rice surface, and a separate omega sensor was used to measure ambient conditions. The precision of the temperature and relative humidity measurements of the omega logger was ±2 and ±5%, respectively. These sensors were placed at different positions, as shown in Figure 5. Total solar radiation was measured by a pyranometer (CMP11, Kipp & Zonen, Delftechpark, Delft, Netherlands). Pyranometers were connected to a data logger (34970A, Agilent Technologies, Inc., Loveland, CO, United States of America). Data were logged at 5-min intervals and merged with the data from the temperature/humidity loggers after each drying experiment.
Immediately after the hourly mixing operation, a sample of about 10 g of paddy rice was taken from each of the positions shown in Figure 5. Due to the preceding mixing, a sufficient mixture of the rice layers was assumed. The samples were weighed by an electronic balance with a precision of ±0.001 (PS-20, Voltcraft, Hirschau, Germany). Dry matter mass mDM (g) was determined by drying the sample in a ventilated oven (ESP-400 Series, Blue M, White Deer, PA, United States of America) at 130 °C for 16 h [59]. The moisture content MC was calculated on dry basis.
M C = m i n i m D M
where mini (g) is the initial mass of the sample.
From the initial mass of paddy rice m0 (kg) used per experimental batch (about 500 kg), the instantaneous mass mt was calculated for one-hour intervals:
m t = m 0 ( 1 + M C 0 )   · ( 1 + M C t )
where MC0 is the initial moisture content and MCt is the instantaneous moisture content at time t.

2.4. Statistical Analysis

2.4.1. Sensitivity Analysis

The parameters with a significant influence in the model were explored with the sensitivity analysis (SA) tool of the Matlab/Simulink® software. The SA was used to test how the variation of the input parameters affects the output behavior of the model [60,61]. Random samples of the individual parameters were generated by using a uniform distribution of its values collected from the literature. Experiments to determine the parameters′ influence on the thermal model were performed using 100 random samples of the exit temperature Tae and the instantaneous mass mt. Then, a Monte Carlo evaluation was carried out on Tae and mt. Finally, the standardized regression coefficient ( R ) was calculated between all the parameters under the partial derivatives by using (20).
R = b x · σ x σ y
where b x is the regression coefficient, σ x is the standard deviation of the corresponding sample, and σ y is the standard deviation of the partial derivatives. A tornado plot shows the results of the SA analysis [61].

2.4.2. Albedo Model

Matlab Version 9.6 software [62] was used for fitting the experimental values of paddy rice albedo by applying the non-linear least squares method to the two proposed models. The model with the highest coefficient of determination (R2) and the lowest root mean square error (RMSE) was selected.
Code developed in Matlab Version 9.6 software was used to assess model performance at constant and variable albedos. R2, RMSE, and mean absolute percentage error (MAPE) were used to indicate the accuracy of the prediction.

3. Results

3.1. Measurements of Albedo

Figure 6 exemplarily shows the course of solar radiation, the solar zenith angle, and the albedo of paddy rice for the rainy season (19 November 2012) and for the dry season (15 April 2013). Fluctuations in albedo were encountered when clouds passed through the experimental area. A higher albedo of paddy rice was observed in the rainy season than in the dry season, with values of 0.29 and 0.25, respectively, at noon. Nevertheless, the obtained values were higher than the mean value calculated from literature data (Table 1).
Based on Equation (17), the albedo can be modeled for the course of the day as a function of the solar zenith angle:
α = 0.2704 + 1.816 · 10 5 · e 0.1761 · θ
With an R2 of 0.7880 and RSME of 0.0299.

3.2. Parameters′ Influence

Figure 7 shows the tornado plot in which parameters are ranked by influence. An example of these results is shown for one batch during the rainy season (4 June 2012). According to the sensitivity analysis, the albedo has the most significant influence, and the other parameters have little effect. This means that at a higher albedo, lower air temperatures are expected. In the case of the instantaneous mass of paddy rice during drying, albedo is also the main parameter of influence, and the other parameters have less impact. Higher albedo values result in lower energy. Therefore, using a constant albedo value in the model may result in the over/underestimation of paddy rice temperature, and consequently, drying rate. Similar results were found for all the experiments and they showed a significant influence of albedo, with a correlation R ranging from −0.81 to −0.67 for the estimation of the exit temperature and from 0.92 to 0.98 for the instantaneous mass during drying.

3.3. Meteorological Conditions during the Drying Experiments

Figure 8 shows an example of the weather conditions during the sun drying of paddy rice during the rainy season (19 October 2011) and dry season (22 May 2012). Weather data collected for the other experiments were similar to the corresponding season. The sun drying conditions in the dry season were characterized by an ambient relative humidity of about 40%; the ambient temperature ranged from 30 to 38 °C, and solar irradiation ranged from 4 to 5 kWh m−2 d−1. Compared with conditions in the rainy season, higher values of relative humidity of about 60%, lower temperatures of about 30 °C, and lower solar radiation values of about 3 kWh m−2 d−1 were expected, primarily due to cloudiness. In both seasons, relative humidity reached values of over 80%, with temperatures of around 25 °C during the night. Furthermore, a higher variability of solar radiation is typically experienced during the rainy season.

3.4. Estimation of Temperatures

Figure 9 shows the simulated and the experimental temperatures over the paddy rice during the dry and rainy seasons. Due to lower cloudiness, a higher temperature was expected during the dry season compared to the rainy season. Nevertheless, the simulated exit temperature in the rainy season showed higher variation at night compared with the experimental temperature over the paddy rice. This variation was not observed during the dry season. During the course of the day, the simulated temperature over the paddy rice for the rainy season showed higher values compared to the dry season. It could be as much as 5 °C higher than the experimental value. Therefore, we can assume that the dynamic sun drying model provides an accurate estimation of the crop temperature.
Figure 10 provides the simulated temperatures of the soil Tso, asphalt Tas, PVC tarpaulin Tta, paddy rice Tpa, and air exiting the paddy rice surface Tae. It can be seen that the paddy rice temperature reached temperatures higher than 50 °C during the dry season. While temperatures reached almost 50 °C or lower during the rainy season.

3.5. Estimation of Moisture Content during Drying

Figure 11 shows a comparison of the simulated and the experimental moisture content of paddy rice during sun drying. The drying curves for the rainy season and the dry season are shown exemplarily. It was observed that the moisture content decreased continuously with drying time. The absence of connection lines in Figure 11 indicates night periods. The MC of the paddy rice for all the experiments decreased from the range of 0.18 to 0.38 to about 0.16 within 3.5 to 76.7 h (Table 3). Although there were slight discrepancies between the simulated and the experimental moisture content, the results obtained by the sun drying model are in good agreement with the experimental data, as the MAPE values were below 2%.

3.6. Accuracy of the Model of Variable vs. Constant Albedo

To assess the model with constant or variable albedo values, RMSE, MAPE, and R2 were used for Tae and mt. Figure 12 shows that for Tae, the RMSE and MAPE are lower in most cases for variable albedo in both seasons, with values for RMSE ranging from 3 to 7.5 and below 10% for MAPE, while R2 shows a minor increment for variable albedo compared to the constant values. Interestingly for the instantaneous mass, RMSE and MAPE show a more significant improvement of the model outputs when a variable albedo is used. Furthermore, the R2 shows slightly higher values of the model with variable albedo compared to the constant albedo (Table 3). Overall, these results indicate that there is a strong influence of variable albedo in the developed sun drying model.

4. Discussion

In previous studies on simulating sun drying, constant albedo values have been used, and the energy exchange within the layers under the paddy rice has been neglected. Moreover, steady state rather than transient state conditions are commonly assumed for solving equations in thermal analysis. Thus, temporary heat and mass transfer analyses were explored and presented in this paper.
In reviewing the literature, Sailor et al.’s [37] approach resulted in a practical procedure to assess the measurement of albedo. We obtained an albedo value of around 0.25–0.28, corresponding to the period from 11:00 to 13:00 (Figure 6). This albedo value is higher than the values reported by other studies such as Arninze et al. [21] and Bala [32], who proposed an albedo of 0.16 for brown rice. On the other hand, the presented range found in this study lies within the range provided by Arya [17], who proposed values for wheat and rice in a range from 0.18 to 0.25. The albedo of paddy rice was shown to correlate with solar zenith angle, showing an R2 of 0.7880 and RSME of 0.0299. Hence, the solar zenith angle should be used for the parametrization of the albedo of paddy rice in modeling schemes in solar drying applications [10,58].
The results of the sun drying model show that when a variable, daytime-dependent albedo was applied, most of the RMSE and MAPE values were smaller for the estimation of the exit temperature and for the instantaneous mass than when using a constant albedo. Another important finding is that RMSE and MAPE with variable albedos are much lower for the instantaneous mass than for the exit temperature, showing that the model provides a better estimation of instantaneous mass. This was already observed from the sensitivity analysis, where albedo showed higher standardized regression values for the instantaneous mass than for the exit temperature.

5. Conclusions

A clear understanding of the incident solar radiation is essential for the effective use of the solar energy.
Albedo, i.e., the solar energy reflected from a surface, is an important component of the energy balance in the solar drying process. This study results suggest that the albedo of the drying product should not be regarded as a constant value, but that the day-dependent variation should be taken into account when estimating energy balances, and thus the course of the moisture content during drying. Furthermore, the development of a mathematical model serves to improve understanding of the dynamic behavior of paddy rice during sun drying. Thus, the modeling of the sun drying of paddy rice can support management decisions for the proper handling of paddy rice.

Author Contributions

Conceptualization, A.S.-R., V.T.-T., and J.M.; data curation, A.S.-R.; formal analysis, A.S.-R., V.T.-T., and J.M.; funding acquisition, A.S.-R. and J.M.; investigation, A.S.-R.; methodology, A.S.-R. and J.M.; resources, A.S.-R. and J.M.; software, A.S.-R. and V.T.-T.; supervision, J.M.; validation, A.S.-R., V.T.-T., and J.M.; visualization, A.S.-R.; writing—original draft preparation, A.S.-R.; writing—review and editing, A.S.-R. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the scholarships awarded by the Carl Zeiss Foundation (Stuttgart, Germany), Fiat Panis Foundation (Ulm, Germany), and Global Rice Science Scholarships—IRRI (Los Baños, Republic of Philippines).

Acknowledgments

We wish to express our appreciation to Marcus Nagle for this academic supervision and Sabine Nugent for language editing, and to the kind and supportive staff members at IRRI, specifically Martin Gummert, Eduardo Secretario, and Elenita Suñaz.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bhattacharya, K.R. 3-Milling Quality of Rice. In Rice Quality; Bhattacharya, K.R., Ed.; Woodhead Publishing: Cambridge, UK, 2013; pp. 61–99. [Google Scholar] [CrossRef]
  2. Verma, D.K.; Srivastav, P.P.; Nadaf, A.B. Nutritional Quality Evaluation of Different Rice Cultivars. In Agronomic Rice Practices and Postharvest Processing: Production and Quality Improvement; Apple Academic Press Inc.: Oakville, Canada, 2018; pp. 299–348. [Google Scholar]
  3. Salvatierra-Rojas, A.; Nagle, M.; Gummert, M.; de Bruin, T.; Müller, J. Development of an inflatable solar dryer for improved postharvest handling of paddy rice in humid climates. Int. J. Agric. Biol. Eng. 2017, 10, 269–282. [Google Scholar] [CrossRef]
  4. Meas, P.; Paterson, A.H.J.; Cleland, D.J.; Bronlund, J.E.; Mawson, A.J.; Hardacre, A.; Rickman, J.F. A mathematical model of solar drying of rice. Int. J. Food Eng. 2012, 8. [Google Scholar] [CrossRef]
  5. Jain, D.; Tiwari, G.N. Thermal aspects of open sun drying of various crops. Energy 2003, 28, 37–54. [Google Scholar] [CrossRef]
  6. Hande, A.R.; Swami, S.B.; Thakor, N.J. Open-Air Sun Drying of Kokum (Garcinia indica) Rind and Its Quality Evaluation. Agric. Res. 2016, 5, 373–383. [Google Scholar] [CrossRef]
  7. Anwar, S.I.; Tiwari, G.N. Evaluation of convective heat transfer coefficient in crop drying under open sun drying conditions. Energy Convers. Manag. 2001, 42, 627–637. [Google Scholar] [CrossRef]
  8. Kumar, M.; Khatak, P.; Sahdev, R.K.; Prakash, O. The effect of open sun and indoor forced convection on heat transfer coefficients for the drying of papad. J. Energy S. Afr. 2011, 22, 40–46. [Google Scholar] [CrossRef]
  9. Togrul, I.T. Determination of convective heat transfer coefficient of various crops under open sun drying conditions. Int. Commun. Heat Mass Transf. 2003, 30, 285–294. [Google Scholar] [CrossRef]
  10. Tsai, J.L.; Tsuang, B.J.; Lu, P.S.; Yao, M.H.; Shen, Y. Surface energy components and land characteristics of a rice paddy. J. Appl. Meteorol. Climatol. 2007, 46, 1879–1900. [Google Scholar] [CrossRef]
  11. Chu, S. Albedo. In Solar and Infrared Radiation Measurements, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2012; pp. 193–204. [Google Scholar]
  12. Chappell, A.; Webb, N.P. Using albedo to reform wind erosion modelling, mapping and monitoring. Aeolian Res. 2016, 23, 63–78. [Google Scholar] [CrossRef] [Green Version]
  13. Cierniewski, J.; Ceglarek, J.; Kaźmierowski, C. Estimating the diurnal blue-sky albedo of soils with given roughness using their laboratory reflectance spectra. J. Quant. Spectrosc. Radiat. Transf. 2018, 217, 213–223. [Google Scholar] [CrossRef]
  14. Jacobs, A.F.G.; Van Pul, W.A.J. Seasonal changes in the albedo of a maize crop during two seasons. Agric. For. Meteorol. 1990, 49, 351–360. [Google Scholar] [CrossRef]
  15. Favero, A.; Sohngen, B.; Huang, Y.; Jin, Y. Global cost estimates of forest climate mitigation with albedo: A new integrative policy approach. Environ. Res. Lett. 2018, 13. [Google Scholar] [CrossRef] [Green Version]
  16. Nkemdirim, L.C. A Note on the Albedo of Surfaces. J. Appl. Meteorol. 1972, 11, 867–874. [Google Scholar] [CrossRef] [Green Version]
  17. Arya, P.S. Radiation Balance Near the Surface. In Introduction to Micrometeorology, 2nd ed.; Holton, J.R., Ed.; Elsevier Science Publishing Co Inc.: San Diego, CA, USA, 2001; Volume 79, pp. 28–45. [Google Scholar]
  18. Ineichen, P.; Guisan, O.; Perez, R. Ground-reflected radiation and albedo. Sol. Energy 1990, 44, 207–214. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, Z.; Barlage, M.; Zeng, X.; Dickinson, R.E.; Schaaf, C.B. The solar zenith angle dependence of desert albedo. Geophys. Res. Lett. 2005, 32, 1–4. [Google Scholar] [CrossRef]
  20. Yang, F.; Mitchell, K.; Hou, Y.-T.; Dai, Y.; Zeng, X.; Wang, Z.; Liang, X.-Z. Dependence of Land Surface Albedo on Solar Zenith Angle: Observations and Model Parameterization. J. Appl. Meteorol. Climatol. 2008, 47, 2963–2982. [Google Scholar] [CrossRef] [Green Version]
  21. Arinze, E.A.; Schoenau, G.J.; Bigsby, F.W. Determination of soar energy absoption and thermal radiative properties of some agricultural products. Trans. Am. Soc. Agric. Eng. 1987, 30, 259–265. [Google Scholar] [CrossRef]
  22. Murata, Y.; Miyasaka, A.; Munakata, K.; Akita, S. On the Solar Energy Balance of Rice Population in Relation to the Growth Stage. Jpn. J. Crop Sci. 1968, 37, 685–691. [Google Scholar] [CrossRef] [Green Version]
  23. Iqbal, M. Ground Albedo. In An Introduction To Solar Radiation; Academic Press: Toronto, Canada, 1983; pp. 281–293. [Google Scholar]
  24. Higuchi, A.; Kondoh, A.; Kishi, S. Relationship among the surface albedo, spectral reflectance of canopy, and evaporative fraction at grassland and paddy field. Adv. Space Res. 2000, 26, 1043–1046. [Google Scholar] [CrossRef]
  25. Iwamoto, A.; Urano, S.I.; Aragaki, M. Water Budget and Estimation of Net Water Requirement for Paddy Field in the Zambezi River Flood Plain, Zambia. J. Agric. Meteorol. 1998, 54, 125–131. [Google Scholar] [CrossRef]
  26. Nakagawa, K.; Ooi, Y. The surface albedo distribution and its seasonal change over the Nagaoka area, Niigata Prefecture, central Japan, estimated with Landsat/MSS data. Geogr. Rev. Jpn. Ser. A 1992, 65, 769–790. [Google Scholar] [CrossRef]
  27. Susaki, J.; Yasuoka, Y.; Kajiwara, K.; Honda, Y.; Hara, K. Validation of MODIS albedo products of paddy fields in Japan. IEEE Trans. Geosci. Remote Sens. 2007, 45, 206–217. [Google Scholar] [CrossRef]
  28. Linacre, E. Estimating Climate Data. In Climate Data and Resources: A Reference and Guide; Routledge: London, UK, 1992; pp. 54–91. [Google Scholar]
  29. Maruyama, A.; Kuwagata, T.; Ohba, K.; Maki, T. Dependence of solar radiation transport in rice canopies on developmental stage. Jpn. Agric. Res. Q. 2007, 41, 39–45. [Google Scholar] [CrossRef] [Green Version]
  30. Gao, Z.; Bian, L.; Zhou, X. Measurements of turbulent transfer in the near-surface layer over a rice paddy in China. J. Geophys. Res. Atmos. 2003, 108, ACL 6-1–ACL 6-13. [Google Scholar] [CrossRef]
  31. Belessiotis, V.; Delyannis, E. Solar drying. Solar Energy 2011, 85, 1665–1691. [Google Scholar] [CrossRef]
  32. Bala, B.K. Fluid Mechanics and Heat Transfer in Solar Drying. In Solar Drying Systems: Simulations and Optimization; Agrotech Publishing Academy: Delhi, India, 1998; pp. 49–86. [Google Scholar]
  33. Berdahl, P.; Bretz, S.E. Preliminary survey of the solar reflectance of cool roofing materials. Energy Build. 1997, 25, 149–158. [Google Scholar] [CrossRef] [Green Version]
  34. Cao, C.; Lee, X.; Muhlhausen, J.; Bonneau, L.; Xu, J. Measuring landscape albedo using unmanned aerial vehicles. Remote Sens. 2018, 10, 1812. [Google Scholar] [CrossRef] [Green Version]
  35. ASTM. E1918-06-Standard Test Method for Measuring Solar Reflectance of Horizontal and Low-Sloped Surfaces in the Field; ASTM International: West Conshohocken, PA, USA, 2006. [Google Scholar]
  36. Qin, Y.; He, H. A new simplified method for measuring the albedo of limited extent targets. Solar Energy 2017, 157, 1047–1055. [Google Scholar] [CrossRef]
  37. Sailor, D.J.; Resh, K.; Segura, D. Field measurement of albedo for limited extent test surfaces. Solar Energy 2006, 80, 589–599. [Google Scholar] [CrossRef]
  38. Das, T.; Bora, G. Greenhouse Solar Thermal Application. In Handbook of Research on Solar Energy Systems and Technologies; Anwar, S., Ed.; IGI Global: Hershey, PA, USA, 2013; pp. 462–479. [Google Scholar]
  39. Brooker, D.B.; Bakker-Arkema, F.W.; Hall, C.W. Grain structure, composition and properties. In Drying and Storage Of Grains and Oilseeds; Springer Science & Business Media: New York, NY, USA, 1992; pp. 19–26. [Google Scholar]
  40. MathWorks. PS Lookup Table (1D). Available online: https://de.mathworks.com/help/physmod/simscape/ref/pslookuptable1d.html (accessed on 28 June 2020).
  41. Duffie, J.A.; Beckman, W.A. Selected Heat Transfer Topics. In Solar Engineering of Thermal Processes, 4th ed.; John Wiley & Sons Inc.: Hoboken, NJ, USA, 2013; pp. 138–172. [Google Scholar]
  42. Palmer, J. The Measurement of Transmission, Absoption, Emission and Reflection. In Handbook of Optics: Devices, Measurements, and Properties, 2nd ed.; Bass, M., Ed.; McGraw-Hill Inc.: New York, NY, USA, 1995; pp. 25.1–25.25. [Google Scholar]
  43. ASABE. Psychrometric Data. ASAE D271.2 APR1979, R2014. Available online: http://elibrary.asabe.org/abstract.asp?aid=32006&t=2 (accessed on 3 June 2020).
  44. Campbell, G.S. Soil Temperature and Heat Flow. In Soil Physics with BASIC: Transport Models for Soil-Plant Systems; Elsevier Science: Amsterdam, The Netherlands, 1985; Volume 14, pp. 26–39. [Google Scholar]
  45. Takakura, T.; Fang, W. Heat Balance of Bare Ground. In Climate Under Cover, 2nd ed.; Springer Science & Business Media: New York, NY, USA, 2002; pp. 45–64. [Google Scholar]
  46. Campbell, G.S.; Norman, J. Radiation Fluxes in Natural Environments. In An Introduction to Environmental Biophysics, 2nd ed.; Springer: New York, NY, USA, 2012; pp. 167–184. [Google Scholar]
  47. Sodha, M.S.; Bansal, P.K.; Dang, A.; Sharma, S.B. Open Sun Drying: An Analytical Study. Dry. Technol. 1985, 3, 517–527. [Google Scholar] [CrossRef]
  48. Bai, B.C.; Park, D.W.; Vo, H.V.; Dessouky, S.; Im, J.S. Thermal Properties of Asphalt Mixtures Modified with Conductive Fillers. J. Nanomater. 2015, 2015. [Google Scholar] [CrossRef] [Green Version]
  49. Sreedhar, S.; Biligiri, K.P. Development of pavement temperature predictive models using thermophysical properties to assess urban climates in the built environment. Sustain. Cities Soc. 2016, 22, 78–85. [Google Scholar] [CrossRef]
  50. Hassn, A.; Chiarelli, A.; Dawson, A.; Garcia, A. Thermal properties of asphalt pavements under dry and wet conditions. Mater. Des. 2016, 91, 432–439. [Google Scholar] [CrossRef]
  51. Subramanian, M.N. Polymers. In Polymer Blends and Composites: Chemistry and Technology; John Wiley & Sons Inc.: Beverly, MA, USA, 2017; pp. 7–55. [Google Scholar]
  52. Crawford, R.J. General Properties of Plastics. In Plastics Engineering, 2nd ed.; Pergamon Press: Oxford, UK, 1987; pp. 1–40. [Google Scholar]
  53. Meas, P. Mathematical Modelling and Improvement of Operating Practices of Sun Drying of Rice: A Thesis Presented in Partial Fulfilment of the Requirements for the Degree of Doctor of Phylosophy [i.e. Philosophy] at Massey University. Ph.D. Thesis, Massey University, Wellington, New Zealand, 2008. [Google Scholar]
  54. Iguaz, A.; San Martín, M.B.; Arroqui, C.; Fernández, T.; Maté, J.I.; Vírseda, P. Thermophysical properties of medium grain rough rice (LIDO cultivar) at medium and low temperatures. Eur. Food Res. Technol. 2003, 217, 224–229. [Google Scholar] [CrossRef]
  55. Garratt, J.R. Energy Fluxes at the Land Surface. In The Atmospheric Boundary Layer; Houghton, J., Rycroft, M., Dessler, A., Eds.; Cambridge University Press: Cambridge, UK, 1994; pp. 115–144. [Google Scholar]
  56. MathWorks. Ordinary Differential Equations. Available online: https://de.mathworks.com/help/matlab/ordinary-differential-equations.html?s_tid=CRUX_lftnav (accessed on 2 March 2020).
  57. Berrizbeitia, S.E.; Gago, E.J.; Muneer, T. Empirical models for the estimation of solar sky-diffuse radiation. A review and experimental analysis. Energies 2020, 13, 701. [Google Scholar] [CrossRef] [Green Version]
  58. Zheng, Z.; Wei, Z.; Wen, Z.; Dong, W.; Li, Z.; Wen, X.; Zhu, X.; Ji, D.; Chen, C.; Yan, D. Inclusion of Solar Elevation Angle in Land Surface Albedo Parameterization Over Bare Soil Surface. J. Adv. Model. Earth Syst. 2017, 9, 3069–3081. [Google Scholar] [CrossRef] [Green Version]
  59. RKB. Measuring rice moisture content. Available online: http://www.knowledgebank.irri.org/training/fact-sheets/postharvest-management/rice-quality-fact-sheet-category/item/measuring-rice-moisture-content-fact-sheet (accessed on 1 February 2020).
  60. Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. Global Sensitivity Analysis. The Primer; John Wiley & Sons: Hoboken, NJ, USA, 2008; pp. 1–292. [Google Scholar] [CrossRef] [Green Version]
  61. MathWorks. Sensitivity Analysis. Available online: https://de.mathworks.com/help/sldo/sensitivity-analysis.html (accessed on 2 March 2020).
  62. The MathWorks Inc. MATLAB, version 9.6.0 (R2019a); The MathWorks Inc.: Natick, MA, USA, 2019. [Google Scholar]
Figure 1. Thermal network of a one-dimensional model for the sun drying of paddy rice with layer thickness z and temperatures of ambient air, Tam, air exiting from the surface, Tae, paddy rice, Tpa, polyvinyl chloride (PVC) tarpaulin, Tta, asphalt, Tas, and soil, Tso.
Figure 1. Thermal network of a one-dimensional model for the sun drying of paddy rice with layer thickness z and temperatures of ambient air, Tam, air exiting from the surface, Tae, paddy rice, Tpa, polyvinyl chloride (PVC) tarpaulin, Tta, asphalt, Tas, and soil, Tso.
Applsci 10 05092 g001
Figure 2. Flow chart of the simulation procedure.
Figure 2. Flow chart of the simulation procedure.
Applsci 10 05092 g002
Figure 3. Setup of the pyranometers for measuring albedo at International Rice Research Institute (IRRI) in Los Baños, Republic of the Philippines (14°9′55″ N, 121°15′37″ E).
Figure 3. Setup of the pyranometers for measuring albedo at International Rice Research Institute (IRRI) in Los Baños, Republic of the Philippines (14°9′55″ N, 121°15′37″ E).
Applsci 10 05092 g003
Figure 4. Top view of shadows created by the circular shield, the pyranometers and the support rod with (a) unshaded paddy rice, (b) unshaded circular shield, (c) pyranometer shadow, (d) support rod shadow, (e) crescent shade of the circular shield, and (f) shaded part of the circular shield.
Figure 4. Top view of shadows created by the circular shield, the pyranometers and the support rod with (a) unshaded paddy rice, (b) unshaded circular shield, (c) pyranometer shadow, (d) support rod shadow, (e) crescent shade of the circular shield, and (f) shaded part of the circular shield.
Applsci 10 05092 g004
Figure 5. Manual mixing during the drying and position of measuring devices, top view and cross-section view, at IRRI in Los Baños, Republic of the Philippines (14°9′55″ N, 121°15′37″ E).
Figure 5. Manual mixing during the drying and position of measuring devices, top view and cross-section view, at IRRI in Los Baños, Republic of the Philippines (14°9′55″ N, 121°15′37″ E).
Applsci 10 05092 g005
Figure 6. Total solar radiation Gup, zenith angle θ, and albedo of paddy rice α during the rainy season (19 November 2012) and the dry season (15 April 2013) at IRRI in Los Baños, Republic of the Philippines (14°9′55″ N, 121°15′37″ E).
Figure 6. Total solar radiation Gup, zenith angle θ, and albedo of paddy rice α during the rainy season (19 November 2012) and the dry season (15 April 2013) at IRRI in Los Baños, Republic of the Philippines (14°9′55″ N, 121°15′37″ E).
Applsci 10 05092 g006
Figure 7. Standardized regression coefficient R from the sensitivity analysis of exit temperature (Tae) and instantaneous mass (mt) of the thermal sun drying model on 4 June 2012 (paddy rice: αpa—albedo, dpa—density, cppa—volumetric heat capacity, λpa—thermal conductivity; PVC tarpaulin: cpta—volumetric heat capacity, λpvc—thermal conductivity; Asphalt: cpas—volumetric heat capacity, λas—thermal conductivity; Soil: cpso—volumetric heat capacity, λso—thermal conductivity, hc heat transfer coefficient from the paddy rice to the air).
Figure 7. Standardized regression coefficient R from the sensitivity analysis of exit temperature (Tae) and instantaneous mass (mt) of the thermal sun drying model on 4 June 2012 (paddy rice: αpa—albedo, dpa—density, cppa—volumetric heat capacity, λpa—thermal conductivity; PVC tarpaulin: cpta—volumetric heat capacity, λpvc—thermal conductivity; Asphalt: cpas—volumetric heat capacity, λas—thermal conductivity; Soil: cpso—volumetric heat capacity, λso—thermal conductivity, hc heat transfer coefficient from the paddy rice to the air).
Applsci 10 05092 g007
Figure 8. Total solar radiation Gup, ambient temperature Tam, and relative humidity rham during an exemplary drying batch in the dry season (22 May 2012) and the rainy season (19 October 2012) at IRRI in Los Baños, Republic of the Philippines (14°9′55″ N, 121°15′37″ E).
Figure 8. Total solar radiation Gup, ambient temperature Tam, and relative humidity rham during an exemplary drying batch in the dry season (22 May 2012) and the rainy season (19 October 2012) at IRRI in Los Baños, Republic of the Philippines (14°9′55″ N, 121°15′37″ E).
Applsci 10 05092 g008
Figure 9. The experimental and simulated exit temperature Tae surface during the rainy season (22 November 2012) (a) and the dry season (4 June 2012) (b).
Figure 9. The experimental and simulated exit temperature Tae surface during the rainy season (22 November 2012) (a) and the dry season (4 June 2012) (b).
Applsci 10 05092 g009
Figure 10. The simulated temperatures during the rainy season (22 November 2011) (a) and the dry season (4 June 2012) (b); temperature: air exiting the paddy rice surface Tae, paddy rice Tpa, PVC tarpaulin Tta, asphalt Tas, and soil Tso.
Figure 10. The simulated temperatures during the rainy season (22 November 2011) (a) and the dry season (4 June 2012) (b); temperature: air exiting the paddy rice surface Tae, paddy rice Tpa, PVC tarpaulin Tta, asphalt Tas, and soil Tso.
Applsci 10 05092 g010
Figure 11. Experimental and simulated moisture content MC during the rainy season (22 November 2012) (a) and dry season (4 June 2012) (b).
Figure 11. Experimental and simulated moisture content MC during the rainy season (22 November 2012) (a) and dry season (4 June 2012) (b).
Applsci 10 05092 g011
Figure 12. The lowest root mean square error RMSE and the mean absolute percentage error MAPE for constant and variable albedo αpa for exit temperature Tae (a) and for predicted instantaneous mass mt (b) during the rainy season and dry season.
Figure 12. The lowest root mean square error RMSE and the mean absolute percentage error MAPE for constant and variable albedo αpa for exit temperature Tae (a) and for predicted instantaneous mass mt (b) during the rainy season and dry season.
Applsci 10 05092 g012
Table 1. Model input parameters from the literature.
Table 1. Model input parameters from the literature.
MaterialParameterDefinitionValues
MinMaxMean
Soil [4,44,45,46,47]λsoThermal conductivity, kJ m−1 K−1 h−10.51905.50002.2598
cpsoVolumetric heat capacity, kJ m−3 K−1200038352780
Asphalt [48,49,50]λasThermal conductivity, kJ m−1 K−1 h−10.88001.86601.3020
cpasVolumetric heat capacity, kJ m−3 K−1218822872247
PVC tarpaulin [51,52]λtaThermal conductivity, kJ m−1 K−1 h−10.14000.15000.1450
cptaVolumetric heat capacity, kJ m−3 K−1143919501694
Paddy rice [39,53,54]λpaThermal conductivity, kJ m−1 K−1 h−10.08950.12500.1068
cppaVolumetric heat capacity, kJ m−3 K−1669891754
ρpaBulk density, kg m−3585600592
[10,17,21,24,30,31,32,55]αpaAlbedo of paddy rice, -0.09000.50000.1957
Table 2. Pyranometer sensor type, properties, and specifications.
Table 2. Pyranometer sensor type, properties, and specifications.
SpecificationsPyranometer Type
CMP6 DownwardCMP11 Upward
Field of view180°180°
Response time18 s5 s
Sensitivity9.05 × 10−6 V W−1 m24.96 × 10−6 V W−1 m2
Temperature dependence of sensitivity<4%<1%
Spectral range285–2800 nm285–2800 nm
Operating temperature−40 to + 80°C−40 to + 80°C
Table 3. Coefficient of determination R2 of air temperatures Tae and instantaneous mass mt at a constant albedo (C) of αpa = 0.1957 and variable (V) albedo αpa.
Table 3. Coefficient of determination R2 of air temperatures Tae and instantaneous mass mt at a constant albedo (C) of αpa = 0.1957 and variable (V) albedo αpa.
SeasonDateMCiniDrying TimeAlbedoTaemt
--hR2R2
----
Rainy12.10.20110.336352.0C0.63960.2029
V0.61140.2182
19.10.20110.249052.3C0.83630.8952
V0.81610.9049
25.10.20110.290976.7C0.65620.8524
V0.67530.8682
14.11.20120.381043.0C0.81990.8611
V0.87290.8160
19.11.20120.353054.4C0.83350.9021
V0.74270.9012
22.11.20120.257849.5C0.87770.8491
V0.83280.8352
Dry22.05.20120.243126.0C0.80120.6700
V0.72340.6722
29.05.20120.244130.0C0.80120.7000
V0.81350.7340
04.06.20120.242733.0C0.87250.8426
V0.88180.8490
09.04.20130.189729.5C0.80620.9328
V0.74710.9405
11.04.20130.16803.5C0.92460.9847
V0.93040.9859
16.04.20130.258929.5C0.76820.9692
V0.76110.9571

Share and Cite

MDPI and ACS Style

Salvatierra-Rojas, A.; Torres-Toledo, V.; Müller, J. Influence of Surface Reflection (Albedo) in Simulating the Sun Drying of Paddy Rice. Appl. Sci. 2020, 10, 5092. https://doi.org/10.3390/app10155092

AMA Style

Salvatierra-Rojas A, Torres-Toledo V, Müller J. Influence of Surface Reflection (Albedo) in Simulating the Sun Drying of Paddy Rice. Applied Sciences. 2020; 10(15):5092. https://doi.org/10.3390/app10155092

Chicago/Turabian Style

Salvatierra-Rojas, Ana, Victor Torres-Toledo, and Joachim Müller. 2020. "Influence of Surface Reflection (Albedo) in Simulating the Sun Drying of Paddy Rice" Applied Sciences 10, no. 15: 5092. https://doi.org/10.3390/app10155092

APA Style

Salvatierra-Rojas, A., Torres-Toledo, V., & Müller, J. (2020). Influence of Surface Reflection (Albedo) in Simulating the Sun Drying of Paddy Rice. Applied Sciences, 10(15), 5092. https://doi.org/10.3390/app10155092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop