Acessibilidade / Reportar erro

DRYING KINETICS OF THE SLICED PULP OF BIOFORTIFIED SWEET POTATO (Ipomoea batatas L.)

ABSTRACT

Biofortified sweet potato (Ipomoea batatas) is one of the foods with the highest contributions of carotenoids in the diet, especially provitamin A carotenoids. Thus, this study aimed to analyze the drying kinetics of the biofortified sweet potato pulp using the Akaike (AIC) and Schwarz's Bayesian (BIC) information criteria for model selection, as well as determine the effective diffusion coefficient and activation energy under different drying conditions. The biofortified sweet potatoes were sliced into chips and submitted to drying in an air circulation oven at 1.0 m s−1 at temperatures of 45, 55, 65, and 75 °C until constant mass. The mathematical models Wang and Singh, Verma, Thompson, Page, Newton, Midilli et al., logarithmic, Henderson and Pabis, two-term exponential, two-term, diffusion approach, frequently used to predict the drying of vegetal products, were adjusted to the data. The Wang and Singh model was selected to represent the drying of the biofortified sweet potato pulp by exhibiting the best adjustment for most conditions. The AIC and BIC criteria were suitable for selecting the Wang and Singh model. The effective diffusion coefficient increased as drying air temperature increase and the activation energy for liquid diffusion was 29.18 kJ mol−1.

KEYWORDS
liquid diffusion; activation energy; Wang and Singh

INTRODUCTION

Sweet potato (Ipomoea batatas Lam.) has a prominent role, as it is one of the most important food crops in the world (Kim et al., 2012Kim HW, Kim JB, Cho SM, Chung MN, Lee YM, Shu SM, Che JH, Kim SN, Kim SY, Cho YS, Kim JH, Park HJ, Lee DJ (2012) Anthocyanin changes in the Korean purple-fleshed sweet potato, Shinzami, as affected by steaming and baking. Food Chemistry 130(4):966-972.). It also has a significant β-carotene content, whose regular intake can prevent and combat blindness and infant mortality caused by vitamin A deficiency, especially in poorer populations that do not have access to other vitamin A sources (Nascimento et al., 2013Nascimento KO, Rocha DGCM, Silva EB, Barbosa Júnior JL, Nascimento MIMJB (2013) Caracterização química e informação nutricional de fécula de batata-doce (Ipomoea batatas L.) orgânica e biofortificada. Revista Verde de Agroecologia e Desenvolvimento Sustentável 8(1): 132-138.).

Drying is a process used to preserve food quality by reducing the availability of water to deterioration reactions, increasing food stability, and reducing the volume and mass of the product (Casarin et al., 2016Casarin F, Mendes CE, Lopes TJ, Moura NF (2016) Planejamento experimental do processo de secagem da amora-preta (Rubus sp.) para a produção de farinha enriquecida com compostos bioativos. Brazilian Journal of Food Technology 19(2) : 116-125.). Reducing the amount of water in agricultural products provides longer shelf life and avoids the degradation of benign substances to the consumer (Tontul & Topuz, 2017Tontul I, Topuz A (2017) Spray-drying of fruit and vegetable juices: Effect of drying conditions on the product yield and physical properties. Trends in Food Science & Technology 63:91-102.). Drying is also the commercial process most used to preserve food and, when compared to other long-term preservative methods, it has a low cost and simpler operation (Alexandre et al., 2013Alexandre HV, Silva FLH, Gomes JP, Silva OS, Carvalho JPD, Lima EE (2013) Cinética de secagem do resíduo de abacaxi enriquecido. Revista Brasileira de Engenharia Agrícola e Ambiental 17(6):640-646.).

The application of reliable mathematical models allows predicting the behavior of the various phenomena that occur during the drying process, which implies the reduction of the operational cost (Dionello et al., 2009Dionello RG, Berbert PA, Molina MB, Pereira RC, Viana AP, Carlesso VO (2009) Secagem de fatias de abacaxi in natura e pré-desidratadas por imersão-impregnação: cinética e avaliação de modelos. Ciência e Tecnologia de Alimentos 29(1):232-240.). Several criteria can be used to verify the adjustment of mathematical models of drying in plant products. However, some parameters have limitations and it is necessary to adopt additional criteria to reinforce and endorse decision-making, such as the Akaike (AIC) and Schwarz's Bayesian (BIC) information criteria (Gomes et al. 2018Gomes FP, Resende O, Sousa EP, Oliveira DEC, Araújo Neto FR (2018) Drying kinetics of crushed mass of ‘jambu': Effective diffusivity and activation energy. Revista Brasileira de Engenharia Agrícola e Ambiental 22(7):499-505.).

Chemical composition and physical structure differ from one product to another, making water outlet specific to each material. In drying studies involving water diffusion there are variations in the values of the effective diffusivity coefficient due to the complexity of plant products, such as different prediction methods, type of material, moisture content, drying process, as well as the methodology used to obtain it (Goneli et al., 2007Goneli ALD, Corrêa PC, Resende O, Reis Neto SA (2007) Estudo da difusão de umidade em grãos de trigo durante a secagem. Ciência e Tecnologia de Alimentos 27(1): 135-140.).

Thus, the aim of this study was to analyze the drying kinetics of the biofortified sweet potato pulp using the AIC and BIC criteria for model selection, as well as determine the effective diffusion coefficient and activation energy for the process.

MATERIAL AND METHODS

The experiment was carried out at the Laboratory of Post-Harvest of Plant Products of the Federal Institute of Education, Science, and Technology of Goiás, campus of Rio Verde. An amount of 3 kg of CNPH1210 biofortified sweet potatoes, with the initial moisture content of 3.796±0.018 (decimal dry basis), were manually peeled using metal knives. Subsequently, the pulp was sliced into chips of approximately 4.6 × 4.0 × 0.2 cm (length, width, and thickness), adapted according to Borges et al. (2008)Borges SV, Mancini MC, Corrêa JLG, Nacimento DA (2008) Secagem de fatias de abóboras (Cucurbita moschata, L.) por convecção natural e forçada. Revista Ciências e Tecnologia de Alimentos 28:245-251., with a domestic grater. The slices of potato pulps were dried in an air circulation oven at 1.0 m s−1 and temperatures of 45, 55, 65, and 75 °C. Samples of 150 g were uniformly distributed in a 0.54 cm layer in rectangular aluminum trays (25 × 10 cm) without drilling, with four replications for each drying temperature.

The reduction of moisture content during the drying was performed by the gravimetric method (loss of mass). The initial moisture content of the product is known and the drying process occurs until reaching a constant mass. The mass reduction was monitored with a 0.01 g resolution scale.

To determine the equilibrium moisture content, samples of biofortified sweet potato were weighed every 24 hours until reaching three weighing measures until constant mass. The moisture contents of the material were determined in an oven regulated at 105 °C for 24 hours.

For determining the moisture content ratios of biofortified sweet potato during drying, [eq. (1)] was used.

(1) RX = X X e X i X e

Where,

  • RX is the moisture content ratio of the product (dimensionless);

  • X is the moisture content of the product (db);

  • Xi is the initial moisture content of the product (db), and

  • Xe is the equilibrium moisture content of the product (db).

Mathematical models frequently used to represent the drying of plant products (Table 1) were adjusted to the experimental data of biofortified sweet potato drying.

TABLE 1
Mathematical models used to predict the drying of plant products.

Mathematical models were adjusted from nonlinear regression analysis by the Gauss-Newton method. The significance of model parameters was evaluated by the t-test adopting a 5% significance level. The degree of adjustment of each model was verified according to the magnitudes of the coefficient of determination (R2), relative mean error (P), and estimated mean error (SE). The relative and estimated mean errors were calculated according to eqs (13) and (14), respectively.

(13) P = 100 n | Y - Y ^ | Y
(14) SE = ( Y Y ^ ) 2 DF

Where,

  • n is the number of experimental observations;

  • Y is the value of the moisture content ratio experimentally observed;

  • Ŷ is the value of the moisture content ratio estimated by model, and

  • DF is the degree of freedom of the model (difference between the number of observations and the number of model parameters). Ŷ

To select a single model for describing the drying process at each condition, the models with the best adjustments were submitted to the Akaike (AIC) and Schwarz's Bayesian (BIC) information criteria.

The AIC and BIC criteria were used as a counterpart in choosing the best mathematical model to predict the phenomenon. AIC allows using the principle of parsimony in choosing the best model, i.e. according to this criterion, the most parameterized model is not always the best (Burnham & Anderson, 2004Burnham KP, Anderson DR (2004) Multimodel Inference: Understanding AIC and BIC in Model Selection. Sociological methods & research 33(2):261-304.).

AIC (Equation 15) is used to compare non-nested models or when three or more models are being compared. Lower AIC values reflect a better adjustment (Akaike, 1974Akaike H (1974) A new look at the statistical model identification. IEEE Transaction on Automatic Control 19:716-723.).

(15) AIC = 2 log l i k e + 2 p

Where,

  • p is the number of parameters, and

  • loglike is the value of the logarithm of the likelihood function considering the estimates of parameters.

BIC (Equation 16) also considers the degree of parameterization of the model and, similarly, the lower the BIC value is, the better the model adjustment. It is an asymptotic criterion whose adequacy is strongly related to the magnitude of the sample size. In relation to the penalty applied in the number of parameters, it will be more rigorous than AIC for small samples.

(16) BIC = 2 log l i k e + p ln ( n )

Where,

  • n is the number of observations used to adjust the curve.?

The effective diffusivity was determined based on the Fick's second law by assuming the geometric shape of an infinite flat plate, which was calculated by [eq. (17)].

(17) RX= 8 π 2 i = 0 1 ( 2 i + 1 ) 2 exp [ ( 2 i+1 ) 2 π 2 D .t 4 . ( S V ) 2 ]

Where,

  • RX is the moisture content ratio (dimensionless),

  • D is the effective diffusion coefficient (m2 s−1), and i is the number of terms.

Volume and surface area of sweet potato slices were determined by measuring with a digital caliper the three orthogonal axes (length, width, and thickness) in fifteen slices of sweet potatoes before drying them. Volume was calculated by [eq. (18)] and surface area by [eq. (19)].

(18) V = S C

Where,

  • V is the volume (m3);

  • S is the surface area (m2), and

  • C is the thickness (m).

(19) S = π ( A . B ) 4

Where,

  • S is the surface area (m2);

  • A is the largest axis (m), and

  • B is the average axis (m).

The relationship between the diffusion coefficient and the drying air temperature was analyzed by the Arrhenius [eq. (20)].

(20) D ef = D 0 exp ( E a R . T a )

Where,

  • Def is the effective diffusion coefficient (m2 s−1),

  • Do is the pre-exponential factor (m2 s−1),

  • Ea is the activation energy (J mol−1),

  • R is the universal gas constant (8.134 kJ kmol−1 K−1), and

  • Ta is the absolute temperature (K−1).

RESULTS AND DISCUSSION

The eleven tested models had coefficients of determination (R2) above 0.96 decimal (Table 2), except for the Verma model (3) for a temperature of 65 °C. According to Sozzi & Ramos (2015)Sozzi G, Ramos DS (2015) Avaliação do ágio no preço da energia convencional no mercado brasileiro de contratos de curto prazo: metodologia e aplicação. Revista Espaço Ciência 22:24-32., the closer to 1 the R2 value is, the more elucidative the model and the better it will be adjusted to the experimental data.

TABLE 2
Statistical parameters: coefficient of determination (R2), Chi-square test (χ2), relative mean error (P), and estimated mean error (SE) obtained for the models (Mod.) adjusted to the drying of slices of biofortified sweet potato pulp.

All drying temperatures presented reduced values (< 0.015) when comparing the SE values with each other and only the Verma model (3) presented a discrepant value, with a high magnitude for the temperature of 55 °C (0.340). For temperatures of 45 and 55 °C, the Wang and Singh (2), Verma (3), Page (5), Midilli et al. (7), logarithmic (8), and two-term (11) models presented values lower than 10% for P, while the Wang and Singh, logarithmic, and Midilli et al. models obtained low values of P for 65 and 75 °C, standing out in the representation of the drying phenomenon.

The models considered satisfactory under all the evaluated temperature conditions were Wang and Singh, logarithmic, and Midilli et al. These three models presented an R2 higher than 99.15% and P values lower than 10%. Thus, through the joint analysis of statistical parameters (R2, P, and SE), these models presented a better adjustment to the drying process of the potato pulp.

In order to select the best model (Wang and Singh, logarithmic, and Midilli et al.), the Akaike (AIC) and Schwarz's Bayesian (BIC) information criteria were adopted as complementary precepts for selection, as shown in Table 3.

TABLE 3
Akaike (AIC) and Schwarz's Bayesian (BIC) selection criteria for the Wang and Singh, logarithmic, and Midilli models.

The joint analysis of the AIC and BIC criteria, in which the lower values indicate better the model adjustment, showed that the values estimated by AIC had the best adjustment for the Midilli et al. at temperatures of 45 and 75 °C and for the Wang and Singh model at temperatures of 55 and 65 °C. In addition, the values estimated by BIC presented better experimental adjustments for the Midilli et al. model at 45 °C and for the Wang and Singh model at temperatures of 55, 65, and 75 °C.

Thus, the Wang and Singh model was selected to represent the drying kinetics of biofortified sweet potato pulp because it presented the best adjustment for most of the drying conditions. Ribeiro et al. (2014)Ribeiro TD, Mattos RWP, Muniz JÁ, Borges SV (2014) Descrição da cinética de secagem de frutos de banana Prata e D'Água por modelos de regressão não linear. Revista da Estatística UFOP 3:78-85. carried out similar studies with sliced banana and Gomes et al. (2018)Gomes FP, Resende O, Sousa EP, Oliveira DEC, Araújo Neto FR (2018) Drying kinetics of crushed mass of ‘jambu': Effective diffusivity and activation energy. Revista Brasileira de Engenharia Agrícola e Ambiental 22(7):499-505. with ground jambu mass and also used the AIC and BIC selection criteria to indicate the most suitable model to represent their drying kinetics.

Figure 1 shows that different temperatures influenced water loss from the material and drying was faster for higher temperatures. According to Fiorentin et al. (2010)Fiorentin LD, Menon BT, Alves JA, Barros STD, Pereira NC, Módenes N (2010) Determinação da cinética e das isotermas de secagem do bagaço da laranja. Acta Scientiarum Technology 32:147-152., an increased temperature accelerates the drying process, promoting moisture content reduction more quickly at the beginning of the drying process, which consequently results in shorter drying time.

FIGURE 1
A) Moisture content ratio during the drying time to different temperatures and B) experimental values estimated by the Wang and Singh model for the drying of slices of biofortified sweet potato pulp at temperatures of 45, 55, 65, and 75 °C.

This behavior was expected and was also recorded by other authors when drying pumpkin slices (Borges et al., 2008Borges SV, Mancini MC, Corrêa JLG, Nacimento DA (2008) Secagem de fatias de abóboras (Cucurbita moschata, L.) por convecção natural e forçada. Revista Ciências e Tecnologia de Alimentos 28:245-251.), fibrous cassava mass (Castiglioni et al., 2013Castiglioni GL, Silva FA, Caliari M, Soares Júnior MS (2013) Modelagem Matemática do Processo de Secagem da Massa Fibrosa de Mandioca. Revista Brasileira de Engenharia Agrícola e Ambiental 17(9):987-994.), lemon slices (Wang et al., 2018Wang J, Law C, Nema PK, Zhao J, Liu Z, Deng L, Gao Z, Xiao H (2018) Pulsed vacuum drying enhances drying kinetics and quality of lemon slices. Journal of Food Engineering 224:129-138.), and green banana peel and pulp (Gonçalves et al., 2017Gonçalves JQ, Silva MAP, Placido GR, Caliari M, Silva RM, Souza DG, Moura LC (2017) Cinética de secagem da casca e polpa da banana verde (Musa acuminata): Propriedades físicas e funcionais da farinha. Global Science and Technology 9:62-72.). In these cases, the highest rates of water vaporization in a shorter drying time was obtained by raising the temperature.

The Wang and Singh mathematical model was the best adjustment for the experimental data of drying of biofortified sweet potato pulp at temperatures of 45 to 75 °C. This result is in accordance with Doymaz (2004aDoymaz I (2004a) Effect of dipping treatment on air drying of plums, Journal of Food Engineering 64:465-470.), who studied the drying of plums, and with Sousa et al. (2016)Sousa FT, Silva MAS, Oliveira DEC, Placido GR, Cagnin C, Moura LC, Souza DG, Caliari M, Lima MS (2016) Modelagem matemática da secagem e propriedades físicas e funcionais do bagaço de malte. Global Science and Technology 9(3):51-61., who studied the drying of malt bagasse.

Table 4 shows the coefficients of the Wang and Singh model. The values found for the parameter b increased as temperature increased, while the parameter a decreased as water removal increased.

TABLE 4
Coefficients a and b obtained from the Wang and Singh model for the drying kinetics of slices of biofortified sweet potato pulp at temperatures of 45, 55, 65, and 75 °C.

The diffusion coefficient increased linearly as the temperature of the drying air increased and its influence can be described by means of the Arrhenius representation, as shown in Figure 2B. During the drying, the diffusion coefficients showed magnitudes of 7.55 × 10−11, 8.9 × 10−11, 13.01 × 10−11, and 19.24 × 10−11 m2 s−1 for temperatures of 45, 55, 65, and 75 °C, respectively. These results are in accordance with those obtained by Fernando et al. (2011)Fernando WJN, Low HC, Ahmad AL (2011) The Effect of Infrared on Diffusion Coefficients and Activation Energies in Convective Drying: A Case Study for Banana, Cassava and Pumpkin. Journal of Applied Sciences 103(4):310-316., who dried slices of banana, cassava, and pumpkin and found values ranging from 0.69 × 10−10 to 0.52 × 10−10 m2 s−1 for banana, 0.80 × 10−10 to 1.12 × 10−10 m2 s−1 for cassava, and 5.76 × 10−10 to 5.60 × 10−10 m2 s−1 for pumpkin.

FIGURE 2
Effective diffusion coefficient (A) and Arrhenius representation for the effective diffusion coefficient (B) obtained for the drying of biofortified sweet potato pulp at temperatures of 45, 55, 65, and 75 °C.

The activation energy for the liquid diffusion of biofortified sweet potato pulp was 29.18 kJ mol−1 for the temperature range between 45 and 75 °C, which is in accordance with the results obtained by Doymaz (2004bDoymaz I (2004b) Convective air drying characteristics of thin layer carrots. Journal of Food Engineering 61(3):359-364.), who studied the drying of sliced carrots with an activation energy of 28.36 kJ mol−1. According to Corrêa et al. (2010)Corrêa PC, Oliveira GHH, Botelho FM, Goneli ALD, Carvalho FM (2010) Modelagem matemática e determinação das propriedades termodinâmicas do café (Coffeaarabica L.) durante o processo de secagem. Revista Ceres 57(5):595-601., the lower the activation energy in the drying processes is, the higher the water diffusivity in the product, i.e. the lower the energy required for the physical transformation to occur and for the liquid water pass to steam (drying of the product).

CONCLUSIONS

Among the studied models, the Wang and Singh model was selected to represent the drying kinetics of biofortified potato pulp since it exhibits the best adjustment for most conditions. The AIC and BIC criteria were suitable for selecting a single model. The effective diffusion coefficient increased as the temperature of the drying air increased and the activation energy for the liquid diffusion was 29.18 kJ mol−1.

ACKNOWLEDGEMENTS

The authors extend thanks to IF Goiano, CAPES, FAPEG, FINEP and CNPq for their financial support, which was indispensable to the execution of this study.

REFERENCES

  • Akaike H (1974) A new look at the statistical model identification. IEEE Transaction on Automatic Control 19:716-723.
  • Alexandre HV, Silva FLH, Gomes JP, Silva OS, Carvalho JPD, Lima EE (2013) Cinética de secagem do resíduo de abacaxi enriquecido. Revista Brasileira de Engenharia Agrícola e Ambiental 17(6):640-646.
  • Borges SV, Mancini MC, Corrêa JLG, Nacimento DA (2008) Secagem de fatias de abóboras (Cucurbita moschata, L.) por convecção natural e forçada. Revista Ciências e Tecnologia de Alimentos 28:245-251.
  • Burnham KP, Anderson DR (2004) Multimodel Inference: Understanding AIC and BIC in Model Selection. Sociological methods & research 33(2):261-304.
  • Castiglioni GL, Silva FA, Caliari M, Soares Júnior MS (2013) Modelagem Matemática do Processo de Secagem da Massa Fibrosa de Mandioca. Revista Brasileira de Engenharia Agrícola e Ambiental 17(9):987-994.
  • Casarin F, Mendes CE, Lopes TJ, Moura NF (2016) Planejamento experimental do processo de secagem da amora-preta (Rubus sp.) para a produção de farinha enriquecida com compostos bioativos. Brazilian Journal of Food Technology 19(2) : 116-125.
  • Corrêa PC, Oliveira GHH, Botelho FM, Goneli ALD, Carvalho FM (2010) Modelagem matemática e determinação das propriedades termodinâmicas do café (Coffeaarabica L.) durante o processo de secagem. Revista Ceres 57(5):595-601.
  • Dionello RG, Berbert PA, Molina MB, Pereira RC, Viana AP, Carlesso VO (2009) Secagem de fatias de abacaxi in natura e pré-desidratadas por imersão-impregnação: cinética e avaliação de modelos. Ciência e Tecnologia de Alimentos 29(1):232-240.
  • Doymaz I (2004a) Effect of dipping treatment on air drying of plums, Journal of Food Engineering 64:465-470.
  • Doymaz I (2004b) Convective air drying characteristics of thin layer carrots. Journal of Food Engineering 61(3):359-364.
  • Fernando WJN, Low HC, Ahmad AL (2011) The Effect of Infrared on Diffusion Coefficients and Activation Energies in Convective Drying: A Case Study for Banana, Cassava and Pumpkin. Journal of Applied Sciences 103(4):310-316.
  • Fiorentin LD, Menon BT, Alves JA, Barros STD, Pereira NC, Módenes N (2010) Determinação da cinética e das isotermas de secagem do bagaço da laranja. Acta Scientiarum Technology 32:147-152.
  • Gomes FP, Resende O, Sousa EP, Oliveira DEC, Araújo Neto FR (2018) Drying kinetics of crushed mass of ‘jambu': Effective diffusivity and activation energy. Revista Brasileira de Engenharia Agrícola e Ambiental 22(7):499-505.
  • Gonçalves JQ, Silva MAP, Placido GR, Caliari M, Silva RM, Souza DG, Moura LC (2017) Cinética de secagem da casca e polpa da banana verde (Musa acuminata): Propriedades físicas e funcionais da farinha. Global Science and Technology 9:62-72.
  • Goneli ALD, Corrêa PC, Resende O, Reis Neto SA (2007) Estudo da difusão de umidade em grãos de trigo durante a secagem. Ciência e Tecnologia de Alimentos 27(1): 135-140.
  • Kim HW, Kim JB, Cho SM, Chung MN, Lee YM, Shu SM, Che JH, Kim SN, Kim SY, Cho YS, Kim JH, Park HJ, Lee DJ (2012) Anthocyanin changes in the Korean purple-fleshed sweet potato, Shinzami, as affected by steaming and baking. Food Chemistry 130(4):966-972.
  • Nascimento KO, Rocha DGCM, Silva EB, Barbosa Júnior JL, Nascimento MIMJB (2013) Caracterização química e informação nutricional de fécula de batata-doce (Ipomoea batatas L.) orgânica e biofortificada. Revista Verde de Agroecologia e Desenvolvimento Sustentável 8(1): 132-138.
  • Ribeiro TD, Mattos RWP, Muniz JÁ, Borges SV (2014) Descrição da cinética de secagem de frutos de banana Prata e D'Água por modelos de regressão não linear. Revista da Estatística UFOP 3:78-85.
  • Sousa FT, Silva MAS, Oliveira DEC, Placido GR, Cagnin C, Moura LC, Souza DG, Caliari M, Lima MS (2016) Modelagem matemática da secagem e propriedades físicas e funcionais do bagaço de malte. Global Science and Technology 9(3):51-61.
  • Sozzi G, Ramos DS (2015) Avaliação do ágio no preço da energia convencional no mercado brasileiro de contratos de curto prazo: metodologia e aplicação. Revista Espaço Ciência 22:24-32.
  • Wang J, Law C, Nema PK, Zhao J, Liu Z, Deng L, Gao Z, Xiao H (2018) Pulsed vacuum drying enhances drying kinetics and quality of lemon slices. Journal of Food Engineering 224:129-138.
  • Tontul I, Topuz A (2017) Spray-drying of fruit and vegetable juices: Effect of drying conditions on the product yield and physical properties. Trends in Food Science & Technology 63:91-102.

Publication Dates

  • Publication in this collection
    25 Apr 2019
  • Date of issue
    Mar-Apr 2019

History

  • Received
    25 July 2018
  • Accepted
    20 Dec 2018
Associação Brasileira de Engenharia Agrícola SBEA - Associação Brasileira de Engenharia Agrícola, Departamento de Engenharia e Ciências Exatas FCAV/UNESP, Prof. Paulo Donato Castellane, km 5, 14884.900 | Jaboticabal - SP, Tel./Fax: +55 16 3209 7619 - Jaboticabal - SP - Brazil
E-mail: revistasbea@sbea.org.br