Acessibilidade / Reportar erro

Adjustment of mathematical models and quality of soybean grains in the drying with high temperatures

Ajustes de modelos matemáticos e qualidade de grãos de soja na secagem com altas temperaturas

ABSTRACT

The aim of this study was to evaluate the influence of the initial moisture content of soybeans and the drying air temperatures on drying kinetics and grain quality, and find the best mathematical model that fit the experimental data of drying, effective diffusivity and isosteric heat of desorption. The experimental design was completely randomized (CRD), with a factorial scheme (4 x 2), four drying temperatures (75, 90, 105 and 120 ºC) and two initial moisture contents (25 and 19% d.b.), with three replicates. The initial moisture content of the product interferes with the drying time. The model of Wang and Singh proved to be more suitable to describe the drying of soybeans to temperature ranges of the drying air of 75, 90, 105 and 120 °C and initial moisture contents of 19 and 25% (d.b.). The effective diffusivity obtained from the drying of soybeans was higher (2.5 x 10-11 m2 s-1) for a temperature of 120 °C and water content of 25% (d.b.). Drying of soybeans at higher temperatures (above 105 °C) and higher initial water content (25% d.b.) also increases the amount of energy (3894.57 kJ kg-1), i.e., the isosteric heat of desorption necessary to perform the process. Drying air temperature and different initial moisture contents affected the quality of soybean along the drying time (electrical conductivity of 540.35 µS cm-1g-1); however, not affect the final yield of the oil extracted from soybean grains (15.69%).

Key words:
dimensioning; optimization; post-harvest.

RESUMO

O objetivo deste estudo foi avaliar a influência do teor de água inicial da soja e as temperaturas do ar de secagem sobre a cinética de secagem, qualidade do grão e encontrar o melhor modelo matemático que se ajustar aos dados de secagem experimental, a difusividade efetiva e o calor isostérico de dessorção. Foi adotado o delineamento experimental inteiramente casualisado (DIC) com esquema fatorial (4 x 2), quatro temperaturas de secagem (75, 90, 105 e 120 ºC) e dois teores de água inicial (25 e 19% b.s.) com três repetições. O teor de água inicial do produto interferiu no tempo de secagem. O modelo de Wang e Singh mostrou-se mais adequado para descrever o processo de secagem dos grãos de soja para as faixas de temperaturas do ar de secagem de 75, 90, 105 e 120 ºC e teores de água iniciais de 19 e 25% (b.u.). A difusividade efetiva obtida da secagem dos grãos de soja foi maior (2,5 x 10-11 m2s-1) para a temperatura de 120 ºC e teores de água de 25% (b.s.). A secagem dos grãos de soja com temperaturas mais elevadas (acima de 105 ºC) e com maior teor de água inicial (25% b.u.) aumentou também a quantidade de energia (3894,57 kJ kg-1), ou seja, o calor isostérico de dessorção necessário para realizar o processo. A temperatura do ar de secagem e dos diferentes teores de água inicial afetou a qualidade dos grãos de soja ao longo do tempo de secagem (condutividade elétrica de 540,35 µS cm-1g-1), no entanto, não interferiu no rendimento de óleo final extraído dos grãos de soja (15,69%).

Palavras-chave:
dimensionamento; otimização; pós-colheita

Introduction

Soybean is one of the main cereals in the world, and Brazil accounts for 40% of the total grain production. This cereal is responsible for 27% of the global market. Brazil is the second largest producer and the largest exporter (CONAB, 2014CONAB - Acompanhamento da safra brasileira: grão safra 2013/2014, segundo levantamento, novembro 2013. Brasília, 2013. <http://www.conab.gov.br>. 17 Mar. 2014.
http://www.conab.gov.br...
). Soybeans are often harvested with water content between 16 and 18%, making it necessary to subject them to drying, in order to prevent enzymatic reactions and the development of fungi and microorganisms that produce hydrolytic enzymes such as lipases, leading to an increase of free fatty acids, which affect the quality of the grains (Ullmann et al., 2010Ullmann, R.; Resende, O.; Sales, J. F.; Chaves, T. H. Qualidade das sementes de pinhão manso submetidas à secagem em diferentes condições de ar. Revista Ciência Agronômica, v.41, p.442-447, 2010. http://dx.doi.org/10.1590/S1806-66902010000300017
http://dx.doi.org/10.1590/S1806-66902010...
). In a mass of moisture grain, this water vapor tends to occupy all the intercellular spaces, generating pressures in all directions.

Moreover, this drying air in the form of water vapor also exerts a partial pressure, designated as partial pressure of water vapor in air (Coradi et al., 2014aCoradi, P. C.; Melo, E. C.; Rocha, R. P. Evaluation of electrical conductivity as a quality parameter of lemongrass leaves (Cymbopogon Citratus Stapf) submitted to drying process. Drying Technology, v.32, n.8, p.969-980, 2014a. http://dx.doi.org/10.1080/07373937.2013.879593
http://dx.doi.org/10.1080/07373937.2013....
). Thus, the drying process aims to partially remove water from the grains, through the simultaneous transfer of heat from the air to the kernels and mass flow through the water vapor, from the beans into the air (Goneli et al., 2011Goneli, A. L. D.; Corrêa, P. C.; Magalhães, F. E. A.; Baptestini, F. M. Contração volumétrica e forma dos frutos de mamona durante a secagem. Acta Scientiarum. Agronomy, v.33, p.1-8, 2011. http://dx.doi.org/10.4025/actasciagron.v33i1.4629
http://dx.doi.org/10.4025/actasciagron.v...
; Siqueira et al., 2012Siqueira, V. C.; Resende, O.; Chaves, T. H. Drying kinetics of Jatropha seeds. Revista Ceres, v.59, p.171-177, 2012. http://dx.doi.org/10.1590/S0034-737X2012000200004
http://dx.doi.org/10.1590/S0034-737X2012...
). According to Corrêa et al. (2010)Corrêa, P. C.; Goneli, A. L. D.; Afonso Júnior, P. C.; Oliveira, G. H. H.; Valente, D. S. M. Moisture sorption isotherms and isosteric heat of sorption of coffee in different processing levels. International Journal of Food Science and Technology, v.45, p.2016-2022, 2010. http://dx.doi.org/10.1111/j.1365-2621.2010.02373.x
http://dx.doi.org/10.1111/j.1365-2621.20...
, information about the behavior of the phenomenon of heat and mass transfer between the biological material and the element of drying is essential for design, operation simulations of drying, dryers and grain quality systems. From the drying curves of soybeans, moisture data versus time can be adjusted (Corrêa et al., 2010Corrêa, P. C.; Resende, O.; Ribeiro, D. M. Drying characteristics and kinetics of coffee berry. Revista Brasileira de Produtos Agroindustriais, v.8, p.1-10, 2006. http://dx.doi.org/10.15871/1517-8595/rbpa.v8n1p1-10
http://dx.doi.org/10.15871/1517-8595/rbp...
; Goneli et al., 2011Goneli, A. L. D.; Corrêa, P. C.; Magalhães, F. E. A.; Baptestini, F. M. Contração volumétrica e forma dos frutos de mamona durante a secagem. Acta Scientiarum. Agronomy, v.33, p.1-8, 2011. http://dx.doi.org/10.4025/actasciagron.v33i1.4629
http://dx.doi.org/10.4025/actasciagron.v...
; Coradi et al., 2014bCoradi, P. C.; Borém, F. M.; Reinato, C. H. Coffee cherries drying process and the influence of environment relative humidity in the mathematical modeling, moisture content, and enthalpy of vaporization. Energia na Agricultura, v.29, p.148-157, 2014b. http://dx.doi.org/10.17224/energagric.2014v29n2p148-157
http://dx.doi.org/10.17224/energagric.20...
).

The study on drying systems, their dimensioning, optimization and determining the feasibility of their commercial application can be made by mathematical simulation, whose principle is based on the successive thin layers of drying of the product, using a mathematical model. The use of mathematical models to simulate the drying process in dryers operating at high temperature is an important tool for engineers working in the field of drying and storing grains (Corrêa et al., 2010Corrêa, P. C.; Goneli, A. L. D.; Afonso Júnior, P. C.; Oliveira, G. H. H.; Valente, D. S. M. Moisture sorption isotherms and isosteric heat of sorption of coffee in different processing levels. International Journal of Food Science and Technology, v.45, p.2016-2022, 2010. http://dx.doi.org/10.1111/j.1365-2621.2010.02373.x
http://dx.doi.org/10.1111/j.1365-2621.20...
; Goneli et al., 2011Goneli, A. L. D.; Corrêa, P. C.; Magalhães, F. E. A.; Baptestini, F. M. Contração volumétrica e forma dos frutos de mamona durante a secagem. Acta Scientiarum. Agronomy, v.33, p.1-8, 2011. http://dx.doi.org/10.4025/actasciagron.v33i1.4629
http://dx.doi.org/10.4025/actasciagron.v...
; Coradi et al., 2014bCoradi, P. C.; Borém, F. M.; Reinato, C. H. Coffee cherries drying process and the influence of environment relative humidity in the mathematical modeling, moisture content, and enthalpy of vaporization. Energia na Agricultura, v.29, p.148-157, 2014b. http://dx.doi.org/10.17224/energagric.2014v29n2p148-157
http://dx.doi.org/10.17224/energagric.20...
). There are numerous kinetic models of drying.

Drying of hygroscopic products can be described by theory, empirical and semi-theoretical methods. Theoretical methods usually consider the external conditions under which the operation occurs, as well as the internal mechanisms of energy transfer and mass and their effects. Brooker et al. (1992)Brooker, D. B.; Bakker-Arkema, F. W.; Hall, C. W. Drying and storage of grains and oilseeds. Westport: The AVI Publishing Company, 1992. 450p. cite the example of theoretical models for the hydrodynamic flow and broadcasts: liquid, capillary surface, steam and heat. Semi-theoretical models are based on Newton's cooling law applied to mass transfer, assuming that the conditions are isothermal and that the water transfer resistance is restricted only to the surface of the product (Brooker et al., 1992Brooker, D. B.; Bakker-Arkema, F. W.; Hall, C. W. Drying and storage of grains and oilseeds. Westport: The AVI Publishing Company, 1992. 450p.). Among semi-theoretical models, the model of two terms, Henderson & Pabis and Page have been widely used (Ullmann et al., 2010Ullmann, R.; Resende, O.; Sales, J. F.; Chaves, T. H. Qualidade das sementes de pinhão manso submetidas à secagem em diferentes condições de ar. Revista Ciência Agronômica, v.41, p.442-447, 2010. http://dx.doi.org/10.1590/S1806-66902010000300017
http://dx.doi.org/10.1590/S1806-66902010...
; Coradi et al., 2014bCoradi, P. C.; Borém, F. M.; Reinato, C. H. Coffee cherries drying process and the influence of environment relative humidity in the mathematical modeling, moisture content, and enthalpy of vaporization. Energia na Agricultura, v.29, p.148-157, 2014b. http://dx.doi.org/10.17224/energagric.2014v29n2p148-157
http://dx.doi.org/10.17224/energagric.20...
). According to Ullmann et al. (2010)Ullmann, R.; Resende, O.; Sales, J. F.; Chaves, T. H. Qualidade das sementes de pinhão manso submetidas à secagem em diferentes condições de ar. Revista Ciência Agronômica, v.41, p.442-447, 2010. http://dx.doi.org/10.1590/S1806-66902010000300017
http://dx.doi.org/10.1590/S1806-66902010...
, the rule of thumb is to form physical dimensionless groups that can be easily investigated by laboratory experiments and is based on external conditions such as temperature, mixing ratio and speed of the drying air. Among these models, there are Thompson, Peart & Foster and Wang & Singh, traditionally used in studies of drying agricultural products and foodstuff (Martinazzo et al., 2010Martinazzo, A. P.; Melo, E. C.; Corrêa, P. C.; Santos, R. H. S. Modelagem matemática e parâmetros qualitativos da secagem de folhas de capim-limão Cymbopogon citratus D. C. StapfRevista Brasileira de Plantas Medicinais, v.12, p.488-498, 2010. http://dx.doi.org/10.1590/S1516-05722010000400013
http://dx.doi.org/10.1590/S1516-05722010...
; Oliveira et al., 2013Oliveira, D. E. C.; Resende, O.; Souza, S. T. A.; Siqueira, V. C.; Neto, C. A. J. Alterações morfométricas em grãos de soja durante o processo de secagem. Semina: Ciências Agrárias, v.3, p.975-984, 2013. http://dx.doi.org/10.5433/1679-0359.2013v34n3p975
http://dx.doi.org/10.5433/1679-0359.2013...
; Nascimento et al., 2015Nascimento, V. R. G.; Biagi, J. D.; Oliveira, R. A. de. Modelagem matemática da secagem convectiva com radiação infravermelha de grãos deMoringa oleifera Revista Brasileira Engenharia Agrícola Ambiental, v.19, p.686-692, 2015. http://dx.doi.org/10.1590/1807-1929/agriambi.v19n7p686-692
http://dx.doi.org/10.1590/1807-1929/agri...
).

Thus, the aim of this study was to evaluate the influence of the initial water content of the soybeans and the temperatures of the drying air on drying kinetics and grain quality, as well as to find the best mathematical model that fit the experimental data of drying, effective diffusivity and isosteric heat of desorption.

Material and Methods

The experiment was conducted at the Federal University of Mato Grosso do Sul (UFMS), Campus of Chapadão do Sul (CPCS), at the Postharvest of Agricultural Products Laboratory, in the year of 2013. The soybean grains were harvested at random with 25 and 19% (d.b.) of moisture content, and the impurities was taken damaged grains using manual separation. Then, the grains were dried in a convection oven with forced air ventilation at temperatures of 75, 90, 105 and 120 °C. Drying was carried out until the seeds reached the hygroscopic moisture balance. For each temperature of the drying air, three tests were performed and, for each test, 5 kg of soybean grains were used. The temperature and relative humidity were monitored throughout drying with the aid of a psychrometer.

The moisture content of grains (% d.b.) was determined by weighing 15 g of sample. Then, the samples were placed in an oven with air heating and ventilation regulated at 105 ± 1 °C for 24 h, according to the recommendations of Brasil (2009)Brasil. Ministério da Agricultura e Reforma Agrária. Regras para análise de sementes. Brasília: SNDA/DNDV/CLAV, 2009.. After that, the samples were removed and placed in desiccators for cooling. The moisture content (% d.b.) was determined by mass difference between the initial and the final sample. Tests were performed in three replicates.

The drying curves were fitted to the experimental data using thirteen different semi-empirical and empirical equations (Eqs. 1 to 11), discriminated below:

- Newton

- Page

- Page Modified

- Henderson & Pabis

- Logarithmic

- Two Terms

- Two Exponential Terms

- Wang & Singh

- Henderson & Pabis

- Midilli

- Diffusion Approximation

where:

RX - moisture ratio, dimensionless;

t - drying time, h;

k, ko, k1 - drying constant, h-1; and,

a, b, c, n - model coefficients.

For determining the ratios of moisture during drying under different conditions, the following expression (Eq. 12) was used:

where:

U - moisture content of the product, d.b.;

Ui - initial moisture content of the product, d.b.; and,

Ue - equilibrium moisture content of the product, d.b.

It is usual to consider the value of the diffusion coefficient constant or linearly dependent on the temperature of the air drying. This relationship has been expressed by the Arrhenius model (Mohapatra & Rao, 2005Mohapatra, D.; Rao, P. S. A. Thin layer drying model of parboiled wheat. Journal of Food Engineering, v.66, p.513-518, 2005. http://dx.doi.org/10.1016/j.jfoodeng.2004.04.023
http://dx.doi.org/10.1016/j.jfoodeng.200...
), Eq. 13:

where:

A - constant, m2 s-1;

E - activation energy, kJ kmol-1;

R - universal gas constant, 8,314 kJ kmol-1 K-1; and,

T - absolute temperature, K.

The coefficients of the Arrhenius expression were linearized by applying the logarithm of the form (Eq. 14):

The values of water activity, temperature and equilibrium moisture content were obtained from the soybean grains desorption isotherms, using the model that best fitted the experimental data. For the calculation of the integral isosteric heat of desorption, it was used the following expression, Eq. 15:

where:

Qst - integral isosteric heat of sorption, kJ kg-1;

L - latent heat of vaporization of free water, kJ kg-1; and,

a, b - coefficients of the model.

The free water latent heat of vaporization (Mohapatra & Rao, 2005Mohapatra, D.; Rao, P. S. A. Thin layer drying model of parboiled wheat. Journal of Food Engineering, v.66, p.513-518, 2005. http://dx.doi.org/10.1016/j.jfoodeng.2004.04.023
http://dx.doi.org/10.1016/j.jfoodeng.200...
) was obtained using the average temperature according to Eq. 16:

where:

L - latent heat of vaporization of free water, kJ kg-1; and,

T m - average temperature, °C.

To adjust the mathematical models of analysis of drying grains, nonlinear regression was performed, through the Quasi-Newton method, using the computer program Statistica 7.0®. To check the degree of fit of each model, the significance of the regression coefficient by t-test was considered, adopting the 0.05 level of probability, the magnitude of the coefficient of determination (R2), the mean relative error values ​​(P) (Eq. 17) and the average estimated error (SE) (Eq. 18) and verified the behavior of distribution of residuals. The relative average error and the average error estimated for each model were calculated according to the following expressions, respectively:

where:

Y - experimentally observed value;

Ŷ - value calculated by the model;

n - number of experimental observations; and,

GLR - degrees of freedom of the model (the number of observations minus the number of model parameters).

The electrical conductivity test was conducted in sunflower seeds, according to the methodology described by Vieira & Krzyzanowski (1999)Vieira, R. D.; Krzyzanowski, F. C. Teste de condutividade elétrica. (Electrical conductivity test). In: Vieira, R. D.; Krzyzanowski, F. C.; França Neto, J. B. Vigor de sementes: conceitos e testes. Londrina: ABRATES, v.4, p.1-26, 1999.. To analyse the yield of oil in soybean grains, the methodology described by Oliveira et al. (2013)Oliveira, D. E. C.; Resende, O.; Souza, S. T. A.; Siqueira, V. C.; Neto, C. A. J. Alterações morfométricas em grãos de soja durante o processo de secagem. Semina: Ciências Agrárias, v.3, p.975-984, 2013. http://dx.doi.org/10.5433/1679-0359.2013v34n3p975
http://dx.doi.org/10.5433/1679-0359.2013...
was used.

The experiment was arranged in a completely randomized design (CRD) in a factorial scheme (4 x 2), and the treatments consisted of four drying air temperatures (75, 90, 105 and 120 °C) and two initial moisture contents of soybean grains (19 and 25% d.b.) with three replicates. The data for oil yield and electrical conductivity were analysed by analysis of variance, using the Tukey test at 0.01 and 0.05 probabilities.

Results and Discussion

The drying process aimed at the partial removal of water from the grains, through the simultaneous transfer of heat and the air mass flow through the water vapor in the air to the grains (Goneli et al., 2011Goneli, A. L. D.; Corrêa, P. C.; Magalhães, F. E. A.; Baptestini, F. M. Contração volumétrica e forma dos frutos de mamona durante a secagem. Acta Scientiarum. Agronomy, v.33, p.1-8, 2011. http://dx.doi.org/10.4025/actasciagron.v33i1.4629
http://dx.doi.org/10.4025/actasciagron.v...
; Coradi et al., 2014aCoradi, P. C.; Melo, E. C.; Rocha, R. P. Evaluation of electrical conductivity as a quality parameter of lemongrass leaves (Cymbopogon Citratus Stapf) submitted to drying process. Drying Technology, v.32, n.8, p.969-980, 2014a. http://dx.doi.org/10.1080/07373937.2013.879593
http://dx.doi.org/10.1080/07373937.2013....
). Figure 1A shows the curves of drying of soybean under different conditions of temperature and initial moisture content. It was found that the drying time of the soybeans, considering the reduction of the moisture content from 0.19 to 0.11 (decimal, d.b.) was 2.33 h for the temperature of 75 °C. Figure 1 shows that the content for the temperature of 120 °C is greater than that observed at 75, 90, and 105 °C. However, as the drying process continues toward the equilibrium moisture content, the drying representative curves are similar, differing in the time of drying. After 1 hour of drying the product subjected to a temperature of 120 °C, the moisture content reached 11.8% (d.b.), while for the drying at 75, 90, and 105 °C, the moisture contents were 17.5, 16.5, and 12% (d.b.), respectively. The drying curves for the soybean grains adjusted to the Wang & Singh model at temperatures of 75, 90, 105 and 120 °C, and initial moisture content of 19% (d.b.). It was found that the satisfactory adjustment was predicted by Wang & Singh curves to the experimental data model, regardless of the conditions tested (Figure 1B).

Figure 1
Moisture ratio (A) adjusted, experimental data and estimated moisture ratio (B) by the Wang & Singh model, for 19% (d.b.) of initial moisture content

The coefficients of the adjusted models, analysed during drying of soybean grains at different drying air temperatures for initial moisture content of 19% (d.b.) are shown in Table 1. Among the models that gave good statistical results (Table 1), the Wang & Singh model was selected to represent the phenomenon of drying of soybean grains due to its simplicity compared to other models and also was selected to present a number of significant coefficients. From the Wang & Singh model, it was observed that the drying coefficients "a" and "b" increased with increasing temperature, except at the temperature of 90 ºC.

Table 1
Parameters obtained, coefficient of determination (R2), mean relative error (P), estimated values of average error (SE) and distribution residue (DR) from Wang and Singh model fitted to the data of drying of soybean grains for 19% (d.b.) of initial water content

The coefficients "a" and "b" represent the effects of external drying conditions, according to Babalis & Belessiotis (2004)Babalis, S. J.; Belessiotis, V. G. Influence of the drying conditions on the drying constants and moisture diffusivity during the thin-layer drying of figs. Journal of Food Engineering, v.65, p.449-458, 2004.http://dx.doi.org/10.1016/j.jfoodeng.2004.02.005
http://dx.doi.org/10.1016/j.jfoodeng.200...
, and can be used as an approach to characterize the effect of temperature related to the effective diffusivity in the drying process period and the descending liquid diffusion controlling the process. On the other hand, the constant "n" reflects the internal strength of the product drying. For Corrêa et al. (2006)Corrêa, P. C.; Resende, O.; Ribeiro, D. M. Drying characteristics and kinetics of coffee berry. Revista Brasileira de Produtos Agroindustriais, v.8, p.1-10, 2006. http://dx.doi.org/10.15871/1517-8595/rbpa.v8n1p1-10
http://dx.doi.org/10.15871/1517-8595/rbp...
, the Verma model satisfactorily represented theCoffea arabica fruits phenomenon of drying. However, Coradi et al. (2014bCoradi, P. C.; Borém, F. M.; Reinato, C. H. Coffee cherries drying process and the influence of environment relative humidity in the mathematical modeling, moisture content, and enthalpy of vaporization. Energia na Agricultura, v.29, p.148-157, 2014b. http://dx.doi.org/10.17224/energagric.2014v29n2p148-157
http://dx.doi.org/10.17224/energagric.20...
) evaluated coffee grains and verified that the Midilli model adjusted best to the kinetics of drying.

Babalis & Belessiotis (2004)Babalis, S. J.; Belessiotis, V. G. Influence of the drying conditions on the drying constants and moisture diffusivity during the thin-layer drying of figs. Journal of Food Engineering, v.65, p.449-458, 2004.http://dx.doi.org/10.1016/j.jfoodeng.2004.02.005
http://dx.doi.org/10.1016/j.jfoodeng.200...
noted that the model of Two Terms showed the best fit to the experimental data of fig and dried seaweeds (Gelidium sesquipedale), respectively. Martinazzo et al. (2010)Martinazzo, A. P.; Melo, E. C.; Corrêa, P. C.; Santos, R. H. S. Modelagem matemática e parâmetros qualitativos da secagem de folhas de capim-limão Cymbopogon citratus D. C. StapfRevista Brasileira de Plantas Medicinais, v.12, p.488-498, 2010. http://dx.doi.org/10.1590/S1516-05722010000400013
http://dx.doi.org/10.1590/S1516-05722010...
found that the best model for drying sheets of lemon grass is the Midilli model, which also considered flat plate and noticed that the drying coefficient "k" increased with increasing temperature, and the constant "n" decreased, confirming the data obtained in this study. A similar result was found by Radünz et al. (2010)Radünz, L. L.; Melo, E. C.; Rocha, P. P.; Berbert, P. A.; Gracia, L. M. N. Study of essential oil from guaco leaves submitted to different drying air temperature. Engenharia na Agricultura v.18, p.241-247, 2010. http://dx.doi.org/10.13083/1414-3984.v18n03a07
http://dx.doi.org/10.13083/1414-3984.v18...
, who studied sage leaves drying conditions under room temperature between 40 and 90 °C and concluded that the models Handerson, Pabis, and Midilli were the ones that the best represented the drying kinetics for sage leaves.

The coefficient of determination is a good parameter to check the fit of nonlinear models and typically has higher values for the models that best fit the experimental data. Thus, it was observed that for most treatments, the coefficients of determination R2 were above 0.97, indicating, according to Corrêa et al. (2006)Corrêa, P. C.; Resende, O.; Ribeiro, D. M. Drying characteristics and kinetics of coffee berry. Revista Brasileira de Produtos Agroindustriais, v.8, p.1-10, 2006. http://dx.doi.org/10.15871/1517-8595/rbpa.v8n1p1-10
http://dx.doi.org/10.15871/1517-8595/rbp...
, a satisfactory representation of the phenomenon under study. The Wang & Singh model showed the lowest values of the mean relative error (P), average estimated error (SE) and distribution of residues for the air drying temperatures of 75, 90, 105, and 120 °C, and initial moisture content of 19% (d.b.) (Table 1).

The successful adjustment obtained by the Wang & Singh model occurred because the drying of soybeans occurred in the falling rate period, wherein the greater resistance to heat and mass transfer occurs from the product to the surface. It appears that during the drying the effective diffusion coefficient increased significantly (p < 0.05), with the rise of temperature (Figure 2A). Almeida et al. (2009)Almeida, D. P.; Resende, O.; Costa, L. M.; Mendes, U. C.; Sales, J. F. Cinética de secagem do feijão adzuki (Vigna angularis). Global Science and Technology, v.2, p.72-83, 2009. found that during the drying of adzuki, effective diffusion coefficients had magnitudes between 0.51x10-10 and 2.23x10-10m2 s-1 for the temperature above at 30 °C.

Figure 2
Mean values for the effective diffusion coefficient (m2s-1) (A), representation of the Arrhenius relationship for the effective diffusivity (B), experimental values and estimated integral isosteric heat of desorption as a function of equilibrium moisture content (C), due to different air temperatures in the drying for 19% (d.b.) of initial moisture content of soybeans grains

The linear fit obtained indicates a uniform variation of diffusivity with temperature, the value being the variation of diffusivity coefficient obtained at 105 and 120 °C, slightly higher than the temperatures of 75 and 90 °C (Figure 2A) This may be explained by the molecular vibration of water, because, according to Goneli et al. (2011)Goneli, A. L. D.; Corrêa, P. C.; Magalhães, F. E. A.; Baptestini, F. M. Contração volumétrica e forma dos frutos de mamona durante a secagem. Acta Scientiarum. Agronomy, v.33, p.1-8, 2011. http://dx.doi.org/10.4025/actasciagron.v33i1.4629
http://dx.doi.org/10.4025/actasciagron.v...
, the effective diffusion coefficient of variation occurs with increasing temperature, which increases the molecular vibration of water molecules and contributes to a faster diffusion. It is observed that the values ​​of Ln (D) as a function of the inverse absolute temperature (1/Ta) show similar behavior for the ranges of temperature 75 to 90 °C and 105 to 120 °C (Figure 2B). Thus, it can be inferred that there was no interference of external conditions on the drying temperatures. The observed activation energy increased with the increase of air temperature of drying (Figure 2C). The activation energy is a barrier that must be overcome so that the diffusion process can be triggered in the product (Kashaninejad et al., 2007Kashaninejad, M.; Mortazavi, A.; Safekordi, A.; Tabil, L. G. Thin - Layer drying characteristics and modeling of pistachionuts. Journal of Food Engineering, v.78, p.98-108, 2007. http://dx.doi.org/10.1016/j.jfoodeng.2005.09.007
http://dx.doi.org/10.1016/j.jfoodeng.200...
; Reis et al., 2011Reis, R. C.; Barbosa, L. S.; Lima, M. L.; Reis, J. S.; Devilla, I. A.; Ascheri, D. P. R. Modelagem matemática da secagem da pimenta Cumari do Pará. Revista Brasileira de Engenharia Agrícola e Ambiental, v.15, p.347-353, 2011.http://dx.doi.org/10.1590/S1415-43662011000400003
http://dx.doi.org/10.1590/S1415-43662011...
).

Figure 3A shows that the drying time of the soybeans, considering the reduction of the moisture content from 0.25 to 0.11 (decimal, d.b.), was 3.7 h for the temperature of 75 °C; however, in the drying at 120 °C the time to complete the process was 0.7 h. Thus, it was observed that, with increasing air temperature, the drying time decreased, showing greater water removal rate of the grains. This happened because of the increased availability of energy for vaporization of the water and increase in the coefficient of mass transfer by increasing the temperature of the drying air. It is observed in Figure 3A that, in the first 0.3 h of the process, the drying rate at 120 and 105 °C is greater than that observed at 75 and 90 °C. However, as the drying process continues toward the equilibrium moisture content, the representative curves of drying at 75, 90, and 105 °C were above the curve at 105 °C, indicating a possible change in the mechanisms that govern the internal movement of moisture. In any event, after 0.3 h drying, the product subjected to a temperature of 120 °C reached water content of 13% (d.b.), while those dried at 75, 90, and 120 °C reached 31, 28, and 17% (d.b.), respectively.

Figure 3
Moisture ratio (A) adjusted, experimental data and estimated moisture ratio (B) to the Wang & Singh model, for 25% (d.b.) of initial moisture content

From 0.3 h of drying, the product being dried at 75 °C has certainly greater amount of free water to provide higher drying rate compared to the product being dried at higher temperatures. The Wang & Singh model showed better fit to the experimental data at drying temperatures of 75, 90, 105 and 120 ºC and initial water content of 25% (d.b.) of soybean, showing similarity to the results obtained for drying and initial water content of 19% (d.b.). In Figures 3A and 3B, it is observed the adjustment of the experimental data and predicted moisture content to the Wang & Singh model to describe the drying of the soybean grains with initial water content of 25% (d.b.). It was found that the coefficient of determination was above 0.99, the mean relative error (P) and estimated average values ​​of error (SE) were below 10% and residues were randomly distributed (Table 2).

Table 2
Parameters obtained, coefficient of determination (R2), mean relative error (P), estimated values of average error (SE) and distribution residue (DR) from Wang & Singh model fitted to the data of drying of soybean grains for 25% (d.b.) of initial moisture content

It can be seen in Figure 4A that the temperature of 120 °C showed the highest effective diffusivity (2.5 x 10-11 m² s-1), indicating that the drying process in these conditions was faster than in the others. This result is consistent, since the temperature of the drying air of 120 °C was the highest one. When comparing the initial moisture contents of 25 and 19% (d.b.), it was found that increasing the moisture in the grains also increases the diffusivity. On the other hand, the drying temperature of 75 °C showed the lowest value of effective diffusivity (1.7 x 10-11 m² s-1) (Figure 4A). It was verified that values ​​of ln (D) as a function of the inverse absolute temperature (1/Ta) showed similar behavior for the ranges of temperature 75 to 90 °C and 105 to 120 °C for initial moisture content of 25% (d.b.) (Figure 4B). According to Figure 4C, reducing the moisture content in the grains increased the energy necessary for the removal of water from the product, represented by integral values of the isosteric heat of desorption (Qst) (Figure 4C) as observed for several agricultural products.

Figure 4
Mean values for the effective diffusion coefficient (m2s-1) (A), representation of the Arrhenius relationship for the effective diffusivity (B), experimental values and estimated integral isosteric heat of desorption as a function of equilibrium moisture content (C), due to different air temperatures in the drying for 25% (d.b.) of initial moisture content of soybeans grains

To remove water from grains with low moisture content, more energy is necessary on average than that required for wet products. These results confirm that, in products with higher moisture content, the force between water molecules and dry material decreases significantly. According to Table 3, it was observed that the variations of oil content for different temperatures and initial moisture contents were not significant at 0.01 and 0.05 probabilities.

Table 3
Analysis of variance of oil yield (%) of soybeans after drying (12% d.b.)

The results obtained in this study (Table 4) confirm those reported in the literature, not detecting immediate effect of drying methods on the yield of soybean oil. Similar results were observed by Oliveira et al. (2013)Oliveira, D. E. C.; Resende, O.; Souza, S. T. A.; Siqueira, V. C.; Neto, C. A. J. Alterações morfométricas em grãos de soja durante o processo de secagem. Semina: Ciências Agrárias, v.3, p.975-984, 2013. http://dx.doi.org/10.5433/1679-0359.2013v34n3p975
http://dx.doi.org/10.5433/1679-0359.2013...
, who studied the effect of five drying temperatures of soybean (ambient air, 40, 60, 80 and 100 °C) and reported that the yield of soybean oil with drying temperature not above 40 °C was significant; however, this value was more visible from a period of more than eight months of storage. When evaluating the effect of temperature on grain drying (20-25, 35-40, 55-60 and 75-80 ºC), Zeni (2010)Zeni, D. B. Efeitos da temperatura na secagem e da umidade relativa no armazenamento sobre a qualidade de grãos e óleo de canola para biocombustível. Pelotas: FAEM/LCTQ/UFPEL, 2010. 82p. Dissertação Mestrado also found no statistically significant difference in an immediate effect of the drying process. Pohndorf (2012)Pohndorf, R. S. Efeitos da umidade e do resfriamento no armazenamento sobre a qualidade de grãos e do óleo de soja para fins comestíveis e de produção de biodiesel. Pelotas: FAEM/LCTQ/UFPEL, 2012. 71p. Dissertação Mestrado, studying the yield of soybean oil extracted from grain stored at different moisture contents (12 and 16% d.b.) and temperatures (8, 13, 18, 23 and 28 ºC) found no statistically significant difference between the treatments employed.

Table 4
Average results for oil yield (%) in soybean after drying (12% d.b.) at different temperatures

In Table 5, it was observed significant variations of the variance at 0.01 probability for the results of electrical conductivity between the drying treatments with different temperatures and with different initial moisture contents.

Table 5
Analysis of variance of the electrical conductivity test (μS cm-1 g-1) in soybeans after drying (12% d.b.)

Analysing the results of Table 6, it was observed that increasing the drying air temperature increased electrical conductivity values ​​of the soybeans.

Table 6
Average results of the electrical conductivity test (μS cm-1 g-1) in soybeans after drying (12% d.b.) at different temperatures and initial moisture content (% d.b.)

The increase in electrical conductivity may be implicated in the major damage caused by the drying air temperature on the cellular structure of the grains during drying, causing them to lose physiological and nutritional quality (Barbosa et al., 2014Barbosa, R. M.; Silva, C. B.; Medeiros, M. A.; Centurion, M. A. P. C.; Vieira, R. D. Condutividade elétrica em função do teor de água inicial de sementes de amendoim. Ciência Rural, v.42, p.45-51, 2014. http://dx.doi.org/10.1590/S0103-84782012000100008
http://dx.doi.org/10.1590/S0103-84782012...
). Likewise, it was found that higher initial contents of water increase the risk of deterioration of grains during drying, due to the longer exposure to the temperatures of the drying air. Ullmann et al. (2010)Ullmann, R.; Resende, O.; Sales, J. F.; Chaves, T. H. Qualidade das sementes de pinhão manso submetidas à secagem em diferentes condições de ar. Revista Ciência Agronômica, v.41, p.442-447, 2010. http://dx.doi.org/10.1590/S1806-66902010000300017
http://dx.doi.org/10.1590/S1806-66902010...
reported a reduction in the percentage of germination of Jatropha while there was an increase in electrical conductivity. According to the authors, the damage at the cellular level in seeds caused by the removal of water from their interior initially can cause physiological and biochemical changes that can be identified by the electrical conductivity test and not be clearly seen in the germination test.

Conclusions

  1. The initial moisture content of the product interferes with the drying time, and the higher the temperature of the drying air (120 °C) and lower the initial moisture content of the grains (19% d.b.), the faster the product reaches the water balance level.

  2. The model of Wang & Singh proved to be more suitable to describe the drying of soybeans for temperature ranges of the drying air.

  3. The effective diffusivity obtained from the drying of soybeans was higher (2.5 x 10-11 m2 s-1) at a temperature of 120 °C and water content of 25% (d.b.).

  4. Drying of soybeans at higher temperatures (above 105 °C) and higher initial moisture content (25% d.b.) also increases the amount of energy (3894.57 kJ kg-1), i.e., the isosteric heat of desorption necessary to carry out the process.

  5. High influence of air temperature and of different initial moisture contents affected the quality of soybean along the drying time (electrical conductivity of 540.35 µS cm-1 g-1), however, does not affect the final yield of the oil extracted from soybean grains (15.69%).

Acknowledgement

The authors thank the FUNDECT - MS for the financial support for the development of this research.

Literature Cited

  • Almeida, D. P.; Resende, O.; Costa, L. M.; Mendes, U. C.; Sales, J. F. Cinética de secagem do feijão adzuki (Vigna angularis). Global Science and Technology, v.2, p.72-83, 2009.
  • Babalis, S. J.; Belessiotis, V. G. Influence of the drying conditions on the drying constants and moisture diffusivity during the thin-layer drying of figs. Journal of Food Engineering, v.65, p.449-458, 2004.http://dx.doi.org/10.1016/j.jfoodeng.2004.02.005
    » http://dx.doi.org/10.1016/j.jfoodeng.2004.02.005
  • Barbosa, R. M.; Silva, C. B.; Medeiros, M. A.; Centurion, M. A. P. C.; Vieira, R. D. Condutividade elétrica em função do teor de água inicial de sementes de amendoim. Ciência Rural, v.42, p.45-51, 2014. http://dx.doi.org/10.1590/S0103-84782012000100008
    » http://dx.doi.org/10.1590/S0103-84782012000100008
  • Brasil. Ministério da Agricultura e Reforma Agrária. Regras para análise de sementes. Brasília: SNDA/DNDV/CLAV, 2009.
  • Brooker, D. B.; Bakker-Arkema, F. W.; Hall, C. W. Drying and storage of grains and oilseeds. Westport: The AVI Publishing Company, 1992. 450p.
  • CONAB - Acompanhamento da safra brasileira: grão safra 2013/2014, segundo levantamento, novembro 2013. Brasília, 2013. <http://www.conab.gov.br>. 17 Mar. 2014.
    » http://www.conab.gov.br
  • Coradi, P. C.; Borém, F. M.; Reinato, C. H. Coffee cherries drying process and the influence of environment relative humidity in the mathematical modeling, moisture content, and enthalpy of vaporization. Energia na Agricultura, v.29, p.148-157, 2014b. http://dx.doi.org/10.17224/energagric.2014v29n2p148-157
    » http://dx.doi.org/10.17224/energagric.2014v29n2p148-157
  • Coradi, P. C.; Melo, E. C.; Rocha, R. P. Evaluation of electrical conductivity as a quality parameter of lemongrass leaves (Cymbopogon Citratus Stapf) submitted to drying process. Drying Technology, v.32, n.8, p.969-980, 2014a. http://dx.doi.org/10.1080/07373937.2013.879593
    » http://dx.doi.org/10.1080/07373937.2013.879593
  • Corrêa, P. C.; Goneli, A. L. D.; Afonso Júnior, P. C.; Oliveira, G. H. H.; Valente, D. S. M. Moisture sorption isotherms and isosteric heat of sorption of coffee in different processing levels. International Journal of Food Science and Technology, v.45, p.2016-2022, 2010. http://dx.doi.org/10.1111/j.1365-2621.2010.02373.x
    » http://dx.doi.org/10.1111/j.1365-2621.2010.02373.x
  • Corrêa, P. C.; Resende, O.; Ribeiro, D. M. Drying characteristics and kinetics of coffee berry. Revista Brasileira de Produtos Agroindustriais, v.8, p.1-10, 2006. http://dx.doi.org/10.15871/1517-8595/rbpa.v8n1p1-10
    » http://dx.doi.org/10.15871/1517-8595/rbpa.v8n1p1-10
  • Goneli, A. L. D.; Corrêa, P. C.; Magalhães, F. E. A.; Baptestini, F. M. Contração volumétrica e forma dos frutos de mamona durante a secagem. Acta Scientiarum. Agronomy, v.33, p.1-8, 2011. http://dx.doi.org/10.4025/actasciagron.v33i1.4629
    » http://dx.doi.org/10.4025/actasciagron.v33i1.4629
  • Kashaninejad, M.; Mortazavi, A.; Safekordi, A.; Tabil, L. G. Thin - Layer drying characteristics and modeling of pistachionuts. Journal of Food Engineering, v.78, p.98-108, 2007. http://dx.doi.org/10.1016/j.jfoodeng.2005.09.007
    » http://dx.doi.org/10.1016/j.jfoodeng.2005.09.007
  • Martinazzo, A. P.; Melo, E. C.; Corrêa, P. C.; Santos, R. H. S. Modelagem matemática e parâmetros qualitativos da secagem de folhas de capim-limão Cymbopogon citratus D. C. StapfRevista Brasileira de Plantas Medicinais, v.12, p.488-498, 2010. http://dx.doi.org/10.1590/S1516-05722010000400013
    » http://dx.doi.org/10.1590/S1516-05722010000400013
  • Mohapatra, D.; Rao, P. S. A. Thin layer drying model of parboiled wheat. Journal of Food Engineering, v.66, p.513-518, 2005. http://dx.doi.org/10.1016/j.jfoodeng.2004.04.023
    » http://dx.doi.org/10.1016/j.jfoodeng.2004.04.023
  • Oliveira, D. E. C.; Resende, O.; Souza, S. T. A.; Siqueira, V. C.; Neto, C. A. J. Alterações morfométricas em grãos de soja durante o processo de secagem. Semina: Ciências Agrárias, v.3, p.975-984, 2013. http://dx.doi.org/10.5433/1679-0359.2013v34n3p975
    » http://dx.doi.org/10.5433/1679-0359.2013v34n3p975
  • Nascimento, V. R. G.; Biagi, J. D.; Oliveira, R. A. de. Modelagem matemática da secagem convectiva com radiação infravermelha de grãos deMoringa oleifera Revista Brasileira Engenharia Agrícola Ambiental, v.19, p.686-692, 2015. http://dx.doi.org/10.1590/1807-1929/agriambi.v19n7p686-692
    » http://dx.doi.org/10.1590/1807-1929/agriambi.v19n7p686-692
  • Pohndorf, R. S. Efeitos da umidade e do resfriamento no armazenamento sobre a qualidade de grãos e do óleo de soja para fins comestíveis e de produção de biodiesel. Pelotas: FAEM/LCTQ/UFPEL, 2012. 71p. Dissertação Mestrado
  • Radünz, L. L.; Melo, E. C.; Rocha, P. P.; Berbert, P. A.; Gracia, L. M. N. Study of essential oil from guaco leaves submitted to different drying air temperature. Engenharia na Agricultura v.18, p.241-247, 2010. http://dx.doi.org/10.13083/1414-3984.v18n03a07
    » http://dx.doi.org/10.13083/1414-3984.v18n03a07
  • Reis, R. C.; Barbosa, L. S.; Lima, M. L.; Reis, J. S.; Devilla, I. A.; Ascheri, D. P. R. Modelagem matemática da secagem da pimenta Cumari do Pará. Revista Brasileira de Engenharia Agrícola e Ambiental, v.15, p.347-353, 2011.http://dx.doi.org/10.1590/S1415-43662011000400003
    » http://dx.doi.org/10.1590/S1415-43662011000400003
  • Siqueira, V. C.; Resende, O.; Chaves, T. H. Drying kinetics of Jatropha seeds. Revista Ceres, v.59, p.171-177, 2012. http://dx.doi.org/10.1590/S0034-737X2012000200004
    » http://dx.doi.org/10.1590/S0034-737X2012000200004
  • Ullmann, R.; Resende, O.; Sales, J. F.; Chaves, T. H. Qualidade das sementes de pinhão manso submetidas à secagem em diferentes condições de ar. Revista Ciência Agronômica, v.41, p.442-447, 2010. http://dx.doi.org/10.1590/S1806-66902010000300017
    » http://dx.doi.org/10.1590/S1806-66902010000300017
  • Vieira, R. D.; Krzyzanowski, F. C. Teste de condutividade elétrica. (Electrical conductivity test). In: Vieira, R. D.; Krzyzanowski, F. C.; França Neto, J. B. Vigor de sementes: conceitos e testes. Londrina: ABRATES, v.4, p.1-26, 1999.
  • Zeni, D. B. Efeitos da temperatura na secagem e da umidade relativa no armazenamento sobre a qualidade de grãos e óleo de canola para biocombustível. Pelotas: FAEM/LCTQ/UFPEL, 2010. 82p. Dissertação Mestrado

Publication Dates

  • Publication in this collection
    Apr 2016

History

  • Received
    06 Dec 2014
  • Reviewed
    18 Sept 2015
  • Accepted
    03 Mar 2016
Unidade Acadêmica de Engenharia Agrícola Unidade Acadêmica de Engenharia Agrícola, UFCG, Av. Aprígio Veloso 882, Bodocongó, Bloco CM, 1º andar, CEP 58429-140, Campina Grande, PB, Brasil, Tel. +55 83 2101 1056 - Campina Grande - PB - Brazil
E-mail: revistagriambi@gmail.com