Estudo das isotermas de equilíbrio do farelo de soja
http://dx.doi.org/10.1590/S0101-20612006000200025
Study of the equilibrium isotherms of soybean meal
Gianini R. LuzII; Luiza H.C.D. SousaII; Luiz M.M. JorgeII; Paulo R. ParaísoII, *
IUniversidade Estadual de Maringá/DEQ, Avenida Colombo, 5790 – CEP 87020-900 – Maringá – Paraná – Brasil – (44) 3261-4752; E-mail: paulo@deq.uem.br
IIDEQ/UEM, e-mail: gianiniluz@yahoo.com.br
IIIDEQ/UEM, e-mail: lmmj@deq.uem.br
IVDEQ/UEM, luizah@deq.uem.br
Este trabalho tem como objetivo estudar o comportamento das isotermas de equilíbrio do farelo de soja. Para tanto, foram realizados experimentos utilizando-se o método estático com soluções salinas saturadas a 50, 60 e 70ºC. Os dados obtidos foram ajustados às equações existentes na literatura, utilizando-se o software Statistica 6.0®. Os resultados indicam que a temperatura não exerce influência significativa sobre os valores da umidade de equilíbrio do farelo de soja na faixa de condições experimentais exploradas e que os modelos de HALSEY e de LUIKOV podem ser utilizados para prever a umidade de equilíbrio do farelo soja.
Palavras-chave: farelo de soja, umidade de equilíbrio, isoterma, secagem.
This work has as aim to study the behavior of the isotherms of the soybean meal. Therefore, experiments had been carried out by using the static method with saturated salt solution at 50, 60 and 70ºC. The data had been adjusted to the equations of the literature, using Statistica software 6.0®. The results indicated that the temperature does not affect significantly the moisture equilibrium values of the soybean meal in the range of temperatures studied and that the LUIKOV and the HALSEY models can be used to foresee the moisture equilibrium of the soybean meal.
Keywords: soybean meal, moisture of the equilibrium, isotherms, drying.
O farelo de soja, resultante da extração do óleo de soja, possui grande importância no mercado agropecuário brasileiro e mundial como ração animal e complemento alimentar humano, conforme salientado por TEIXEIRA [12].
O farelo de soja, conforme PARAÍSO [7], representa 70% da massa sólida do grão de soja, sendo que 51% correspondem às proteínas, 43% aos carboidratos e 6% às cinzas.
A soja passa por várias fases até a obtenção do farelo de soja descritas por SARTORI [10] e PARAÍSO [7], que são: silagem, limpeza, secagem dos grãos, quebra dos grãos, cozimento, laminação, expansão, extração, dessolventização-tostagem, secagem do farelo de soja, peletização, moagem e expedição.
A secagem do farelo de soja é uma etapa importante deste processo, pois aumenta o tempo de validade do farelo e reduz custos de transporte. Esta operação, normalmente apresenta elevado consumo de energia, que pode ser de até 80 kg de vapor por tonelada de soja processada, aumentando os custos operacionais e prejudicando o meio ambiente em termos de aquecimento e emissões poluentes, segundo JONGENEELEN [5]. Sendo assim, na otimização do processo de secagem do farelo de soja, torna-se desejável estudar o processo, a fim de reduzir o consumo de energia e os danos ao meio ambiente. Para tanto, normalmente desenvolve-se um modelo matemático do processo de secagem, cuja simulação e otimização realistas dependem do conhecimento de dados de equilíbrio confiáveis do farelo de soja.
A umidade de equilíbrio é a umidade que o produto irá apresentar a uma determinada umidade relativa do ar e em uma determinada pressão e temperatura. Os dados de umidade de equilíbrio variam com o tipo de sólido e podem ser determinados, experimentalmente, por métodos estáticos ou dinâmicos, conforme BARROZO [1].
Estes dados de umidade de equilíbrio raramente são encontrados na literatura, uma destas exceções é o trabalho de PIXTON e WARBURTON [9], que apresentaram dados para farelo de soja a 15, 25 e 35ºC, uma faixa de temperatura adequada para o armazenamento. Podendo notar-se, a partir daí, a necessidade do conhecimento das umidades de equilíbrio para o farelo a temperaturas superiores a 35ºC com a finalidade de aplicação em processos industriais. Esses dados experimentais podem ser ajustados usando equações da literatura adequadas para cereais.
As correlações encontradas na literatura, normalmente, são empíricas e semi-empíricas, uma vez que nenhum modelo teórico foi capaz de predizer com precisão o teor de umidade de equilíbrio em todas as faixas de temperatura, conforme apresentado por BROOKER et al. [2].
As equações utilizadas para o ajuste neste trabalho foram as de LUIKOV (1978), de HALSEY (1976) e de HENDERSON-THOMPSON (1968), obtidas em SOUSA et al. [11] para estimar umidade de equilíbrio e podem ser vistas na Tabela 1.
Neste contexto, o objetivo deste trabalho foi determinar, experimentalmente, os dados de umidade de equilíbrio do farelo de soja, utilizando soluções salinas saturadas a 50, 60 e 70ºC, segundo a metodologia descrita por BARROZO [1], MUJUNDAR [6] e GREENSPAN [4], a qual é baseada em um procedimento estático. Posteriormente, os cálculos foram realizados no software Statistica 6.0® através do método de iteração de Rosenbrock e Quasi-Newton.
A dessorção do farelo foi realizada com amostras de 0,15 cm de tamanho médio, retiradas na saída do equipamento denominado de Dessolventizador/Tostador, fornecidas pela Cocamar (Cooperativa Agrícola de Maringá), localizada na Região Noroeste do Paraná, empregando-se soluções salinas de LiCl, CH3CO2K, K2CO3, NaNO2, NaCl e KCl que foram armazenadas juntamente com a amostra em vidros cilíndricos. Os recipientes apresentavam diâmetro de 7,6 cm, altura de 12 cm e tela localizada dentro deste com a altura 6 cm em relação ao fundo do vidro e sobre a solução para sustentar o farelo durante o processo de desidratação, conforme a Figura 1.
2.2 - Metodologia experimental
A desidratação das amostras de farelo de soja foi realizada nos laboratórios da UEM, segundo o procedimento descrito abaixo. Cada solução salina foi colocada cuidadosamente no fundo dos recipientes e as amostras de 0,025 kg ficaram dispostas sobre uma tela. Em seguida, fecharam-se os recipientes hermeticamente.
Posteriormente, as amostras foram colocadas em estufa a 50ºC e pesadas até atingir peso constante em três medidas sucessivas. Considerou-se que, normalmente este processo leva em torno de 21 dias para ser concluído. Na etapa seguinte, coloca-se a amostra, de forma isolada, para secar na estufa sob temperatura fixa de 105±2ºC, e então, é pesada após 24 h ou até atingir peso constante (o que indica que a secagem foi concluída). Também é medida a umidade relativa, UR, e a temperatura do ar, T. O processo foi repetido para as temperaturas de 60 e 70ºC. Todas as medidas foram realizadas em triplicata para dar confiabilidade aos dados obtidos. Os sais foram escolhidos a partir dos dados apresentados por MUJUNDAR [6].
2.3 - Metodologia de ajuste das isotermas de equilíbrio
Com os dados da umidade de equilíbrio, Xe, da umidade relativa do ar, UR, e as equações de LUIKOV, HALSEY e HENDERSON-THOMPSON, ajustaram-se os parâmetros de cada uma dessas equações pela técnica da regressão não-linear, utilizando o software Statistica 6.0®.
Posteriormente, conhecendo-se os parâmetros ajustados, An, construíram-se as isotermas de equilíbrio e obtiveram-se os coeficientes de correlação, R2, e o somatório dos desvios quadráticos, f, sendo f=S(Xe|modelo - Xe|experimental)2 , a fim de verificar qual delas melhor representa os dados experimentais.
Na Tabela 2, são apresentados os valores médios experimentais de umidade de equilíbrio, Xe, e umidade relativa, UR, para cada solução salina nas temperaturas 50, 60 e 70ºC. As concentrações das soluções salinas foram especificadas conforme PERRY e CHILTON [8], a fim de obterem-se os valores de UR apresentados na Tabela 2.
Observou-se que, para 50ºC, a umidade de equilíbrio do farelo variou de 0,034 a 0,240 kg H2O/kg sól. seco, em uma faixa de umidade relativa entre 11,1 e 79,9%. Já para 60ºC, a umidade de equilíbrio do farelo variou de 0,029 a 0,233 kg H2O/kg sól. seco e em uma faixa de umidade relativa entre 11 e 78,2% e, para 70ºC, a umidade de equilíbrio do farelo variou de 0,023 a 0,221kg H2O/kg sól. seco e em uma faixa de umidade relativa entre 10,9 a 76,5%.
Com isso, verificou-se que à medida que a temperatura do ar de secagem aumenta, a umidade de equilíbrio do farelo sofre uma pequena variação para uma determinada umidade relativa do ar, sendo que este comportamento foi observado por GEANKOPLIS [3] e BARROZO [1] em estudos de umidade de equilíbrio para cereais.
Na Tabela 3, são apresentados os resultados dos parâmetros obtidos para cada uma das equações quando ajustadas aos valores experimentais nas condições de 50, 60 e 70ºC. A avaliação do melhor modelo de ajuste se dá pelo maior coeficiente de correlação, R2, e o menor somatório dos desvios quadráticos.
A análise da qualidade de ajuste das isotermas de equilíbrio do farelo de soja apresentados na Tabela 3 indica que as equações de HALSEY e de LUIKOV se destacam em relação ao modelo de HENDERSON-THOMPSON, por apresentarem melhor coeficiente de correlação, R2. Ou seja, para as equações de HALSEY e LUIKOV os R2 obtidos foram de 0,97, enquanto que para HENDERSON-THOMPSON, foi de 0,95. Assim, para diferenciar os modelos de HALSEY e LUIKOV, calculou–se o somatório dos desvios quadráticos, f, de cada um e o resultado foi 0,006 para ambos, indicando que os dois modelos são adequados para a representação dos dados de umidade de equilíbrio.
Na Figura 5, são apresentados os dados experimentais da desidratação do farelo de soja classificados de acordo com a temperatura.
A análise do comportamento dos pontos experimentais apresentados na Figura 5 indica que a temperatura exerce pouca influência sobre a umidade de equilíbrio para o farelo de soja na faixa de condições experimentais explorada.
As Figuras 3, 4 e 5 foram construídas com o objetivo de apresentar as isotermas de equilíbrio, obtidas a partir das equações de HENDERSON-THOMPSON, HALSEY, LUIKOV, comparando-as com os dados experimentais nas temperaturas de 50, 60 e 70ºC de forma isolada.
Para tornar o estudo mais completo, somaram-se aos valores experimentais os dados apresentados por PIXTON e WARBURTON [9] para 15, 25 e 35ºC, os quais são apresentados na Tabela 4 e conjuntamente na Figura 6.
A Tabela 5 apresenta os parâmetros obtidos para cada uma das equações ajustadas conjuntamente aos dados deste trabalho e aos dados de PIXTON e WARBURTON [9].
Os resultados da Tabela 5 mostram que o melhor modelo matemático é o de LUIKOV, pois apresenta o maior coeficiente de correlação de 0,99 e menor somatório dos desvios quadráticos de 0,019.
A Figura 6 tem a finalidade de mostrar graficamente o grau de ajuste das equações de HALSEY, LUIKOV, HENDERSON-THOMPSON para os todos dados da Tabela 5.
Conforme previamente comentado e, apresentado na Figura 2, a temperatura não exerce influência significativa sobre os dados experimentais de umidade de equilíbrio do farelo de soja na faixa de 50 a 70ºC. Esta observação foi confirmada acrescentando os dados de PIXTON E WARBURTON [9] obtidos na faixa de 15 a 35ºC e foram apresentados na Figura 6. Desta forma, testaram-se as correlações tradicionais da literatura sem a variável temperatura, conforme a Tabela 6. Adicionalmente, propôs-se um modelo alternativo que relaciona a umidade de equilíbrio com a umidade relativa de forma exponencial.
A Tabela 7 apresenta os parâmetros obtidos para cada uma das equações modificadas, ajustadas conjuntamente aos dados deste trabalho e aos dados de PIXTON e WARBURTON [9].
Ao comparar as Tabelas 5 e 7, pode-se averiguar que o modelo de LUIKOV continua a ser o mais adequado de acordo com os resultados estatísticos. A diferença pode ser constatada apenas nos somatórios dos desvios quadráticos na terceira casa decimal, comprovando a interferência não significativa da temperatura na estimativa da umidade de equilíbrio do farelo de soja, o que também pode ser comprovado pelas Figuras 6 e 7. Apesar de ser proposto um modelo exponencial, ele não mereceu discussão por apresentar resultados insatisfatórios.
É interessante notar que, para uma grande faixa de umidade relativa do ar (0 a 60%), a umidade de equilíbrio do farelo de soja varia muito pouco e permanece em níveis inferiores a 0,15 kg H2O/kg sól. seco, representando uma faixa operacional atraente para processos industriais de secagem do farelo de soja. Este comportamento também é apresentado por GEANKOPLIS [3] para produtos alimentícios.
Os dados experimentais de 50, 60 e 70ºC também foram utilizados para o ajuste das equações modificadas, cujos parâmetros, coeficientes de correlação e somatório dos desvios médios são mostrados na Tabela 8.
Na Tabela 8, observou–se que os modelos de HALSEY e LUIKOV modificados continuam os mais adequados para estimar a umidade de equilíbrio do farelo de soja.
As isotermas de equilíbrio do farelo de soja, de acordo com os parâmetros apresentados na Tabela 8, são mostradas na Figura 8.
Ao comparar os modelos originais e modificados para os dados nas condições de 15, 25, 35, 50, 60 e 70ºC, averiguou–se que LUIKOV apresentou os melhores coeficientes de ajuste e para tanto, ambos podem ser observados na Figura 9.
Pela Figura 9, nota-se que a diferença entre os resultados previstos para os modelos aparece quando a umidade relativa está na faixa de 0,8 a 1.
Os dados experimentais para o farelo de soja indicam que as equações de HALSEY e de LUIKOV são as mais adequadas para predizer a umidade de equilíbrio do sólido.
Ao comparar as isotermas de equilíbrio para as diferentes condições de temperatura, verificou-se que esta exerce pouca influência sobre a umidade de equilíbrio do farelo de soja, sendo comprovada gráfica e estatisticamente pelos modelos na forma modificada ou original.
Pode-se observar que, de acordo com isotermas construídas, para uma grande faixa de umidade relativa do ar (0 a 60%), a umidade de equilíbrio do farelo de soja varia muito pouco e permanece em níveis inferiores a 0,15 kg H2O/kg sól. seco.
5 - NOMENCLATURA
A1, A2, A3 – Parâmetros de ajuste
T – Temperatura absoluta do ar (K)
UR – Umidade relativa do ar (%/100)
Xe – Umidade de equilíbrio do farelo kg H2O/kg sól. seco
F – Somatório dos desvios quadráticos
[1] BARROZO, M.A.S.; OLIVEIRA, D.T., SANCINETI, G.P., RODRIGUES, M.V. A study of the desorption isotherms of lentils. Brazilian Journal of Chemical Engineering, v. 17, p. 105-109, 2000. [ Links ]
[2] BROOKER, D.B.; BAKKER-ARKEMA, F.W.; HALL, C.W. Drying and storage of grains and oilseeds. The AVI Publishing Company, New York, p. 450, 1992. [ Links ]
[3] GEANKOPLIS, C. Transport processes and unit operations. USA: Prentice Hall, 1993. [ Links ]
[4] GREENSPAN, L. Humidity fixed – points of binary satured aqueous solutions. Journal Res. Natl. Bureau of Standards, 81(A), p. 89-93, 1977. [ Links ]
[5] JONGENEELEN, H.P.J. Energy conservation in solvent extration plants, JAOCS Journal of the American Oil Chemists Society, v .53, p. 291, June, 1976. [ Links ]
[6] MUJUNDAR, A.S. Handbook of industrial drying. New York: Marcel Dekker, Inc., 1987. [ Links ]
[7] PARAÍSO, P.R. Modelagem e análise do processo de obtenção do óleo de soja. Campinas, 2001, 200 p. Tese de doutorado, Faculdade de Engenharia Química, Universidade de Campinas (Unicamp), 2001. [ Links ]
[8] PERRY, R.H.; CHILTON, C.H. Manual de Engenharia Química, 5ª ed. Rio de Janeiro, Guanabara Dois S. A, 1980. [ Links ]
[9] PIXTON, S.W.; WARBURTON. The moisture content/equilibrium relative humidity relationship of soya meal. Journal of Stored Products Res. Report, n. 1.109, USDA, 1975. [ Links ]
[10] SARTORI, A.E.R. Modelagem, simulação e controle de um secador rotativo industrial. Maringá, 2001, 92 p. Dissertação (mestrado em Desenvolvimento de Processos), Faculdade de Engenharia Química, Universidade Estadual de Maringá (UEM), agosto, 2001. [ Links ]
[11] SOUSA, L.H.C.D.; PERRI, V.R.S.; PEREIRA, N.C.; MOTTA LIMA, O.C.; MENDES, E.S. Isotermas de equilíbrio para tecidos de algodão – adsorção e dessorção. In: XV Congresso Brasileiro de Engenharia Química, 2004, Curitiba, Paraná. [ Links ]
[12] TEIXEIRA, C.G. A soja na alimentação humana. Disponível em: http://www.vegetarianismo.com.br/artigos/soja1.html . Acesso em 18 abr. 2005. [ Links ]