Estado térmico de produtos armazenados em silos com sistema de aeração: estudo teórico e experimental

Thermal state of stored products in storage bins with aeration system: experimental-theoretical study

ARTIGOS CIENTÍFICOS

CIÊNCIA E TECNOLOGIA PÓS-COLHEITA

Estado térmico de produtos armazenados em silos com sistema de aeração: estudo teórico e experimental

Thermal state of stored products in storage bins with aeration system: experimental-theoretical study

Fabiane A. de Oliveira^I; Oleg A. Khatchatourian^II; Anderson Bihain^III

^IFísica, Profa. Dra., Departamento de Física, Estatística e Matemática, UNIJUI, Ijuí - RS, Fone: (0XX55) 3332.0204, [e-mail]

^IIEngº Mecânico, Prof. Doutor, Departamento de Física, Estatística e Matemática, UNIJUI, Ijuí - RS

^IIIBolsista de Iniciação Científica - CNPq, UNIJUI

RESUMO

O resfriamento de grãos por ventilação de ar ambiente aplica-se amplamente nas etapas finais de secagem e para controle posterior da temperatura de grãos armazenados. O objetivo deste trabalho é realizar estudo teórico-experimental sobre o estado térmico de massa de grãos de soja sujeita à aeração. Foram obtidos dados experimentais sobre a dinâmica de resfriamento de massa pré-aquecida de grãos de soja para diferentes alturas da coluna de grãos e velocidades do ar. A análise dos resultados mostrou que a taxa de resfriamento varia significativamente durante todo o processo e em todo o domínio e a difusividade térmica das camadas não poderia ser considerada constante. Para simular a dinâmica de resfriamento, foram apresentados dois modelos matemáticos. No primeiro modelo, o domínio de resfriamento foi dividido pela fronteira móvel em duas zonas, representadas por diferentes difusividades térmicas (análogo de problema de Stefan). No outro modelo, todo o domínio foi dividido hipoteticamente em pequenas camadas e foi considerado que o processo de equilíbrio térmico entre o ar e a massa de grãos para essas camadas é atingido instantaneamente ("reatores homogêneos"). Os resultados de simulações mostraram concordância satisfatória com os dados experimentais.

Palavras-chave: armazenamento, aeração, resfriamento.

ABSTRACT

The cooling of grains by ventilation with ambient air applies widely in the final stages of drying and for the subsequent temperature control of stored grains. The objective of this work is to carry out an experimental-theoretical study on the thermal state of a soybean mass, subject to the aeration. The experimental data on cooling dynamics of preheated soybean mass were obtained for different heights of the grain column and air velocities. The analysis of the results showed that the cooling rate varies significantly during all process in all domains and the thermal diffusivity of the layers cannot be considered as a constant. To simulate the cooling dynamics two mathematical models were introduced. In the first model the cooling domain was divided by means of the movable frontier in two zones represented by different coefficients of thermal diffusivity (analogous of Stefan's problem). In another model all the domain was divided hypothetically in small layers and was considered that the thermal balance between the air and the mass of grains for these layers is reached instantly ("homogeneous reactors"). The results of simulations showed a satisfactory agreement with the experimental data.

Keywords: storage, aeration, cooling.

INTRODUÇÃO

Durante o armazenamento de grãos em silos, esses ficam sujeitos às alterações da temperatura ambiente que, junto com o teor de água da massa de grãos, são responsáveis pela evolução da população de insetos e microorganismos que influenciam significativamente na qualidade do produto. Para preservar essa qualidade, utiliza-se o processo de aeração, que consiste no resfriamento dos grãos através da passagem forçada de ar ambiente pela massa de grãos.

Apesar do número considerável de recomendações sobre administração de armazéns, os projetistas e operadores de sistemas de armazenamento encontram várias dificuldades, desenvolvendo e operando os sistemas de aeração de grãos. A modelagem matemática adequada do funcionamento desses sistemas e a simulação da distribuição dos parâmetros em silos, como a pressão estática, a velocidade do ar e a temperatura poderiam facilitar significativamente o aperfeiçoamento do armazenamento.

Atualmente, existem modelos que descrevem satisfatoriamente o escoamento do ar em silos em condições isotérmicas (KHATCHATOURIAN & SAVICKI, 2004). Ao mesmo tempo, a simulação do estado térmico de um silo para condições não-isotérmicas apresenta algumas dificuldades relacionadas com a ausência de dados confiáveis sobre coeficientes de transferência de calor entre o grão e o ar, coeficientes de difusividade térmica do próprio grão e da camada de grãos.

Vários autores tentaram descrever a dinâmica de variação de temperatura da massa de grãos. Alguns deles, como, por exemplo, BAKKER-ARKEMA & BICKERT (1966) e KHATCHATOURIAN & OLIVEIRA (2005, 2006), admitiram que o fenômeno de transferência de massa de água pode ser desprezado. Outros, como BAKKER-ARKEMA et al. (1974) e MONTROSS & MAIER (2000), levaram em conta a transferência de calor e de massa em camadas de grãos. Esses últimos modelos seriam mais adequados para a simulação de secagem de grãos, apesar de, às vezes, estarem unidos pelo termo "Aeration models" (NAVARRO & NOYES, 2002).

Os objetivos deste artigo são: obter dados experimentais da temperatura da massa de grãos de soja em função do tempo e da posição na coluna para diferentes velocidades do ar de aeração; desenvolver modelo matemático que descreva a dinâmica de resfriamento da massa de soja; realizar simulações numéricas utilizando dois modelos matemáticos e compará-las com os dados experimentais obtidos, e realizar a escolha do modelo que melhor se adapte e descreva o processo de resfriamento/aquecimento de grãos.

MATERIAL E MÉTODOS

Estudo Experimental

Equipamento

Os dados experimentais do resfriamento foram obtidos de um pequeno "silo" composto por um tubo de PVC (diâmetro de 0,15 m e 1,0 m de altura) com isolamento térmico nas laterais. O esquema da bancada é mostrado na Figura 1.

Os grãos de soja com teor de água médio de 12 % (b.u) foram previamente selecionados e limpos e para a realização dos experimentos, e aquecidos em estufa até temperatura aproximada de 52 ºC (temperatura inicial - T₀), durante aproximadamente 8 h, suficiente para uniformizar a temperatura em toda a massa de grãos.

O resfriamento dos grãos aquecidos foi realizado por meio da ventilação de uma camada de grãos de soja de 0,60 m de altura pelo ar ambiente, por um ventilador centrífugo, acionado por motor elétrico com ¾ cv. Esse ventilador estava ligado a um inversor de freqüência que controlava as velocidades do ar do ventilador. O ar era impulsionado e passava por uma tubulação de 0,05 m de diâmetro onde uma placa de orifício, ligada a um manômetro, media a diferença de pressão do ar na placa. A partir dessa medida, foi possível calcular a vazão do ar que atravessava a placa, dada pela expressão abaixo:

em que,

C = 0,72 - coeficiente da placa de orifício calibrada utilizada nos experimentos;

_ra - densidade do ar;

Dp - diferença de pressão na placa medida pelo manômetro, e

A₀ = 0,001017 m² - área do orifício da placa.

Experimentos preliminares realizados na bancada para o intervalo de velocidades do ar estudadas neste trabalho (0,05 a 0,5 m s^-1) mostraram que a diferença da temperatura do ar e do grão durante o resfriamento não é significativa (excluindo a secção próxima da entrada do ar, nos instantes iniciais). Por isso, a temperatura dos grãos foi medida por termopares modelo MTK - tipo K, inseridos dentro da massa de grãos ao longo do tubo, nas seguintes secções da coluna de grãos: y = 0,15 m; y = 0,27 m; y = 0,40 m; y = 0,54 m e ligados a uma placa de aquisição de dados A/D do tipo LR 7018, que permitia realizar até dez leituras por segundo. Na entrada da câmara, foi medida a temperatura do ar de entrada utilizado para aeração (temperatura ambiente - T_a).

Em condições reais, a diferença entre as temperaturas no interior dos silos e a temperatura do ar ambiente (utilizado para a aeração) raramente excedem 10 ºC, porém, nas experiências realizadas com resfriamento de grãos de soja, essa diferença foi mantida entre 20 ºC e 25 ºC, com o objetivo de aumentar a sensibilidade das medidas e reduzir o erro relativo.

Resultados

Na Figura 2, os pontos apresentam os dados experimentais da temperatura da massa de grãos de soja obtidos em diferentes secções do tubo para uma velocidade do ar na entrada do tubo constante.

A análise dos resultados mostra que a taxa de resfriamento varia significativamente durante todo o processo, em todas as secções. Para secções mais próximas da entrada do ar fresco, o resfriamento ocorre mais rápido. A dinâmica de resfriamento com diferentes velocidades para a distância y = 0,27 m é mostrada na Figura 3 e possui o mesmo comportamento.

Para secções mais distantes da entrada, podem-se identificar três intervalos principais de variação de taxa de resfriamento: A) no primeiro intervalo (início do processo de resfriamento), a taxa de resfriamento não é muito alta (resfriamento lento). Isso ocorre porque a temperatura do ar que retira o calor dos grãos das secções anteriores, mais próximas à entrada de ar, cresce e reduz a capacidade do ar de resfriar as camadas mais distantes. A extensão desse intervalo depende da secção analisada e da velocidade do ar; B) no segundo intervalo, quando o resfriamento da massa de grãos nas secções iniciais já ocorreu, a taxa de resfriamento aumenta consideravelmente e, por conseqüência, ocorre aceleração do processo de resfriamento das camadas mais distantes. A parte principal do processo de resfriamento encontra-se neste intervalo, e C) no terceiro intervalo, novamente ocorre resfriamento lento, onde a temperatura de grãos assintoticamente tende para a temperatura inicial do ar.

No início, a camada considerada está sendo ventilada por ar aquecido pelas camadas anteriores, "esperando" a chegada do ar frio. À medida que ocorre o resfriamento das camadas anteriores, o ar que passa pela camada considerada possui a temperatura cada vez menor. A taxa de resfriamento dos grãos, que é proporcional à diferença entre as temperaturas do grão e do ar, aumenta. Posteriormente, essa diferença reduz-se devido à diminuição da temperatura do grão pelo resfriamento. Todos esses fatores implicam que a difusividade térmica do meio particulado (massa de grãos) varia com a temperatura, que, por sua vez, varia tanto com o tempo em cada secção quanto com a altura para cada instante de tempo considerado.

O aumento da velocidade do ar acelera o processo de resfriamento, fornecendo à massa de grãos mais ar frio por unidade do tempo (Figura 3). Os resultados obtidos permitiram avaliar dois modelos matemáticos do processo de resfriamento apresentados na próxima seção.

Estudo Teórico

Modelos matemáticos

Inicialmente, aplicando o modelo de estado térmico de uma barra com isolamento térmico lateral com temperatura fixa em uma das extremidades e com transferência de calor por convecção na outra, é possível mostrar que os resultados obtidos desse modelo não descrevem com boa concordância os dados experimentais, trabalhando com a difusividade térmica a e o coeficiente de transferência de calor h constantes durante todo o processo de resfriamento. Por outro lado, a variação desses coeficientes torna o problema não-linear e exige o conhecimento das dependências desses coeficientes com a temperatura.

Neste trabalho, estão apresentados dois modelos matemáticos que permitiram obter razoável concordância entre dados experimentais e simulados para a descrição da dinâmica de resfriamento de massa de grãos previamente aquecidos sujeitos à aeração à temperatura ambiente.

No primeiro modelo (modelo 1), para transformar o problema não-linear de propagação do calor numa seqüência de problemas lineares com contorno móvel, foi desenvolvido modelo matemático parecido com o problema de Stefan.

Seja o domínio dividido por dois subdomínios e em que S(t) é alguma fronteira móvel inicialmente desconhecida. Os subdomínios W1 e W2 possuem diferentes difusividades térmicas e as temperaturas T₁(y,t) e T₂(y,t) nesses subdomínios satisfazem às correspondentes equações de condução de calor:

em que,

T(y,t) - temperatura da massa de grãos;

y - posição na coluna de grãos, e

a - difusividade térmica da massa de grãos.

Um esquema do problema pode ser visualizado na Figura 4.

Conforme o modelo proposto e satisfazendo às condições de continuidade da solução, foi obtida a temperatura constante da fronteira móvel (T_S) de separação dos dois domínios, como no Problema de Stefan.

A posição da fronteira móvel S, em função do tempo, foi escolhida na forma linear, dada pela expressão:

em que,

k - constante de proporcionalidade e representa a velocidade da frente de resfriamento que depende da velocidade do ar no tubo.

Conforme a posição dessa fronteira, a camada de grãos foi dividida em dois intervalos separados pela posição da fronteira.

À direita da fronteira, relacionando a posição da fronteira com a nova origem, foi introduzida uma nova coordenada dada por:

Fazendo a mudança de variáveis, substituindo as expressões (4) e (5) na equação de condução de calor unidimensional (eq.2), tem-se :

As condições de contorno para esse caso são:

A temperatura inicial da massa de grãos (T₀) é constante para toda a massa de grãos.

A solução da eq.(6) foi obtida na seguinte forma:

A distribuição da temperatura no intervalo à esquerda da fronteira pode ser descrita pela equação de condução de calor para uma barra limitada com extremos mantidos a temperaturas constantes:

A solução da eq.(3) foi obtida na forma:

em que,

Para esse caso, j(y) é a distribuição inicial da temperatura.

Na fronteira que separa os subdomínios, estão satisfeitas as condições:

O problema considerado é intermediário entre o problema de estado térmico de um corpo constituído por dois materiais diferentes (com fronteira fixa) e o problema de Stefan, no qual a fronteira móvel separa os domínios de fases diferentes. As difusividades térmicas de cada domínio foram escolhidas minimizando a discrepância entre os dados experimentais e simulados.

No segundo modelo (modelo 2), o estado térmico da massa de grãos sujeita à aeração foi calculado considerando a conservação de energia do sistema composto pelos grãos de soja e o ar que entra em contato com esses grãos. Para realizar os cálculos, todo o domínio de integração foi dividido em elementos retangulares com nodos em pontos n= 0,1,...,N. Nesse caso, M representa o número de divisões espaciais e N o número de divisões temporais.

A função contínua desconhecida T = T(y,t) foi substituída por um conjunto de pontos isolados Tⁱ(n) @ T(yi,t(n)).

Admitindo que o equilíbrio térmico entre o ar e o grão se estabelece instantaneamente, isto é, o ar e os grãos em cada camada considerada sempre têm as temperaturas idênticas, a energia cedida pela massa de grãos em uma camada Dy_i por um tempo Dt seria igual à energia recebida pelo ar. Dessa forma, pode-se escrever:

Aplicando a regra de trapézios para o cálculo da integral:

que representa a temperatura média temporal do ar na entrada de camada "i" no intervalo e considerando , foi obtida a fórmula de recorrência para o cálculo da temperatura das camadas de grãos em correspondentes instantes de tempo (KHATCHATOURIAN & OLIVEIRA, 2006):

sendo i=1,2,..., M e n=0,1,2,..., N. Os índices a e g nas eqs.(13), (14) e (15) referem-se, respectivamente, ao ar e à camada de grãos de soja. Para as propriedades da camada de grãos (c_g, r_g e e), foram utilizados os dados de MOHSENIN (1980), citado por ANDRADE et al. (2001).

Análogo ao cálculo de parâmetros de reatores homogêneos, as camadas obtidas pela divisão espacial foram chamadas de "reatores homogêneos" (KHATCHATOURIAN & OLIVEIRA, 2006).

RESULTADOS E DISCUSSÃO

Aplicando os modelos apresentados, foram realizadas as simulações numéricas do processo de resfriamento da massa de grãos para diferentes velocidades do ar.

Para realizar as simulações utilizando o modelo 1, foi necessário determinar os valores das difusividades térmicas de cada domínio, a posição da fronteira móvel e os valores da temperatura nessa fronteira para cada velocidade do ar estudada. Esses valores foram obtidos a partir da análise dos dados experimentais.

As simulações da posição da fronteira móvel S(t) foram realizadas utilizando a eq.(4), na qual os valores de k foram calculados pelo Método de Mínimos Quadrados (M.M.Q), cuja variação com a velocidade do ar é mostrada na Figura 7. Pelos resultados da simulação apresentados na Figura 6, pode-se notar que, conforme a hipótese inicial [eq.(4)], a posição da fronteira varia linearmente com o tempo e depende da velocidade do ar.

A relação entre a velocidade da frente de resfriamento (velocidade da fronteira) e a velocidade do ar é mostrada na Figura 7.

Na Figura 8, são mostrados os valores da difusividade térmica do primeiro e do segundo domínios para diferentes velocidades do ar.

Pode-se observar que a difusividade térmica a₁ do primeiro domínio é maior que a₂ do segundo domínio. Com aumento de velocidade, os valores absolutos de coeficientes a₁ e a₂ crescem.

A comparação com dados experimentais (Figuras 2 e 3) mostrou que tanto o primeiro como o segundo modelo descrevem razoavelmente o processo de resfriamento.

Uma discrepância entre os dados experimentais e calculados, no caso do primeiro modelo, é causada, provavelmente, pela linearização do modelo, escolhendo as difusividades térmicas constantes. Na Figura 10, são mostradas as curvas obtidas pelos dois modelos teóricos junto com os dados experimentais para os valores da temperatura adimensional ao longo da coluna de grãos para um tempo fixo.

No segundo modelo, os passos temporal Dt = 2 s e espacial Dy_i = 0,01m foram escolhidos de maneira que a diminuição posterior deles não influa sobre o resultado dos cálculos. As simulações mostraram que esse modelo tem melhor desempenho em comparação com o outro para o intervalo estudado de variação de velocidade. Foi calculado o coeficiente de correlação (R²) para cada curva da Figura 9. Os valores obtidos para o modelo 1 estão no intervalo entre 0,9676 e 0,9897, e para o modelo 2, entre 0,9919 e 0,9957.

Uma boa coincidência entre a simulação utilizando o modelo 2 e os dados experimentais aprova a hipótese de equilíbrio térmico. Isso significa que o valor da difusividade térmica do próprio grão é bastante alto para que o tempo de permanência do ar na mesma camada (que depende da velocidade) seja compatível com o tempo característico de aquecimento do grão (que depende da difusividade térmica). Era esperado que o aumento da velocidade atrasasse o resfriamento da parte periférica do grão, o que levaria parte do ar a passar sem aquecimento completo, de acordo com o modelo de reator homogêneo. Nesse caso, a temperatura simulada de resfriamento deveria ser menor que a temperatura experimental da massa de grãos. O efeito contrário foi observado na Figura 9, indicando que, para a posição escolhida da junção do termopar no silo, na realidade, foi medida a temperatura da superfície do grão, sendo essa igual à temperatura do ar. Essa temperatura é sempre menor que a temperatura média do grão, por isso as temperaturas simuladas ficaram acima das obtidas nos experimentos. Para evitar (ou reduzir) essa discordância, o modelo pode ser complementado por um coeficiente de não-uniformidade (que varia entre 0 e 1). Uma grande vantagem desse modelo simples é a possibilidade de realizar as simulações em silos com secção transversal variável. Nesse caso, o campo de velocidades tem de ser calculado previamente (KHATCHATOURIAN & SAVICKI, 2004).

CONCLUSÕES

O coeficiente de difusividade térmica das camadas não pode ser considerado constante, variando com o tempo, durante o resfriamento.

Foi mostrado que a divisão de toda a câmara pela fronteira móvel em dois domínios com coeficientes de difusividade térmica diferentes, permite resolver o problema do estado térmico nessas condições. Foi mostrado que a posição da fronteira móvel depende da velocidade do ar e varia linearmente com o tempo. Foram encontradas as dependências da difusividade térmica do primeiro e segundo domínios para diferentes velocidades do ar.

Foi proposto um modelo ("reatores homogêneos") que permitiu evitar a utilização da difusividade térmica para a simulação do resfriamento no intervalo estudado de variação dos parâmetros.

A comparação das simulações com dados experimentais mostrou que ambos os modelos descrevem razoavelmente o processo de resfriamento.

A aplicação do segundo modelo para simulação do estado térmico de um silo é preferível pela melhor coincidência com os dados experimentais e pela facilidade de considerar as secções transversais variáveis.

Recebido pelo Conselho Editorial em: 7-6-2006

Aprovado pelo Conselho Editorial em: 22-3-2007

• ANDRADE, E.T.; MOLIN, L.; COUTO, S.M.; CARDOSO, G.E. Distribuição de temperatura em grãos armazenados em silos: modelagem por elementos finitos. Revista Brasileira de Armazenamento, Viçosa, n.26, p.28-35, 2001.
• BAKKER-ARKEMA, F.W.; BICKERT, W.G. A deep-bed computational cooling procedure for biological products. Transactions of the ASAE, St. Joseph, v.9, n.6, p.834-6, 1966.
• BAKKER-ARKEMA, F.W.; LEREW, L.E.; DEBOER, S.F.; ROTH, M.G. Grain dryer simulation. Agricultural Experiment Station Report East Lansing: Michigan State University, 1974. n.224. 80 p.
• KHATCHATOURIAN, O.A.; SAVICKI, D.L. Mathematical modelling of airflow in an aerated soya bean store under non-uniform conditions. Biosystems Engineering, Edinburg, v.88, n.2, p.201-11, 2004.
• KHATCHATOURIAN, O.A.; de OLIVEIRA, F.A. Mathematical modelling of air flow and thermal state in large aerated grain storage. In: INTER-AMERICAN DRYING CONFERENCE, 3., 2005, Montreal. Proceedings.. Montreal: McGill University, 2005. 1 CD-ROM.
• KHATCHATOURIAN, O.A.; de OLIVEIRA, F.A. Mathematical modelling of airflow and thermal state in large aerated grain storage. Biosystems Engineering, Edinburg, v.95, n.2, p.159-69, 2006.
• MOHSENIN, N.N. Thermal properties of foods and agricultural materials London: Gordon & Breach Publishers, 1980. 404 p.
• MONTROSS, M.D.; MAIER, D.E. Simulated performance of conventional high-temperature drying, dryaeration and combination drying of shelled corn with automatic conditioning. Transactions of the ASAE, St. Joseph, v.43, n.3, p.691-9, 2000.
• NAVARRO, S.; NOYES, R.T. (Eds.). The mechanics and physics of modern grain aeration management New York: CRC Press, 2002. 647 p.

Datas de Publicação

• Publicação nesta coleção
  15 Abr 2008
• Data do Fascículo
  Abr 2007

Histórico

• Recebido
  07 Jun 2006
• Aceito
  22 Mar 2007