terça-feira, 18 de novembro de 2008

Transferência de calor transiente na agitação linear intermitente de latas

RESUMO

Foi estudada a transferência de calor transiente na agitação linear e intermitente (ALI) de embalagens metálicas contendo simulantes de alimentos, objetivando-se sua aplicação em processos de pasteurização ou esterilização e conseqüentes tratamentos térmicos mais eficientes, homogêneos e com produto de melhor qualidade. Foram utilizados quatro meios fluidos simulantes de alimentos de diferentes viscosidades e massas específicas: três óleos e água. Foram combinados efeitos de cinco tratamentos, sendo: meio simulante (4 níveis), espaço livre (3 níveis), freqüência de agitação (4 níveis), amplitude de agitação (2 níveis) e posição das latas (4 níveis). Os ensaios de aquecimento e resfriamento foram feitos em tanque com água à temperatura de 98 °C e 17-20 °C, respectivamente. Com os dados de penetração de calor em cada experimento, foram calculados os parâmetros de penetração de calor fh, jh, fc e jc. Os resultados foram modelados utilizando-se grupos de números adimensionais e expressos em termos de Nusselt, Prandtl, Reynolds e funções trigonométricas (com medidas de amplitude e freqüência de agitação, espaço livre e dimensões da embalagem). Foram estabelecidas as duas Equações gerais para as fases de aquecimento e resfriamento:

Nu = ReA –0,199.Pr –0,288.sen(xa/AM)0,406.cos(xf/FA)–1,039.cos((xf/FA).(EL/H).p)–4,556 Aquecimento

Nu = 0,1295.ReA–0,047.Pr –0,193.sen(xa/AM)0,114.cos(xf/FA)–0,641.cos((xf/FA).(EL/H).p)–2,476 Resfriamento

O processo de ALI pode ser aplicado em pasteurizadores ou autoclaves estáticas horizontais e verticais, com modificações simples. Concluiu–se que a ALI aumenta significativamente a taxa de transferência de calor, tanto no aquecimento como no resfriamento.

1 Introdução

O processamento térmico é um processo de preservação de alimentos baseado na aplicação do calor cujo sucesso depende de vários parâmetros, como a difusividade térmica do alimento, as características da embalagem, as taxas de transferência de calor, tempo e temperatura efetivamente aplicados no alimento, dentre outros.

Processos de esterilização e pasteurização sob agitação vêm sendo estudados desde 19504,14,15,19,21,13, porém a agitação linear e intermitente (ALI) praticamente não é citada. A vantagem das autoclaves e pasteurizadores agitados está na elevada taxa de transferência de calor do meio de aquecimento para a embalagem. Esse aumento nas taxas de transferência de calor acarreta melhora da qualidade do alimento, levando a uma otimização do processo térmico.

A transmissão de calor está associada a uma diferença de temperatura entre duas regiões e é governada pela combinação de várias leis da física6,17. São reconhecidos pela literatura três mecanismos fundamentais para transmissão de calor: condução, convecção e radiação4,6.

Para condução foi desenvolvida experimentalmente uma relação matemática por Fourier, em 1822, descrita na Equação 1, que estabelece que o calor transmitido (q) é diretamente proporcional à área (A) que ele flui e ao gradiente de temperatura (T) na direção (x) do fluxo de calor. O coeficiente de proporcionalidade é definido como a condutividade térmica (k) do material4.

Casos de transferência de calor transiente em alimentos são comuns e, neste caso, a equação diferencial é descrita nos três eixos do sólido (x, y, z), da seguinte forma:

onde o primeiro termo da Equação 2 é a difusividade térmica do material. Ela representa, fisicamente, a relação entre a condução e a capacidade de armazenamento de calor no material (rCp = densidade x calor específico) e indica a velocidade com que a frente de temperatura se propaga no sentido do fluxo de calor4.

Os alimentos aquecidos tipicamente por condução são carnes em conservas, vegetais e frutas altamente compactadas na embalagem, além de produtos amiláceos em estado gelatinizado.

O calor transmitido por convecção é regido por uma relação empírica conforme a Equação 3:

A Equação 3 é uma definição do coeficiente médio de transferência de calor local, . Esse coeficiente é de difícil avaliação por ser a convecção um fenômeno complexo, envolvendo a geometria da embalagem, a rugosidade, o espaço livre da embalagem, a viscosidade, a massa específica, o calor específico, a condutividade térmica e alterações estruturais do alimento em função da temperatura3.

A efetividade da transmissão de calor por convecção depende do movimento da mistura do fluido e, conseqüentemente, o estudo desse mecanismo fundamenta-se no conhecimento das características físicas de escoamento do fluido3.

Teoricamente, as soluções dos problemas matemáticos de transferência de calor para fluidos puros podem ser obtidas analiticamente utilizando-se as três equações de transporte de energia. Contudo, a solução analítica destas equações é muito difícil para casos de alimentos sob agitação, nos quais os perfis de velocidades são muito complexos e as condições de contorno difíceis de serem estabelecidas3.

Na convecção, o transporte de energia é feito pela ação combinada de condução de calor, armazenamento de energia e movimento da mistura. A convecção pode ser natural ou forçada, respectivamente se o agente da mistura for ocasionado simplesmente por diferença da massa específica das partículas ou por um agente externo, como a agitação mecânica3.

Os alimentos aquecidos simplesmente por convecção são, por exemplo, sucos de frutas, sopas e produtos lácteos. Já os aquecidos por condução e convecção são os produtos sólidos imersos em salmoura ou xarope e sopas contendo partículas vegetais ou carne1.

Existem produtos que quando aquecidos inicialmente se comportam como convectivos, mas que durante o aquecimento passam a ter comportamento tipicamente condutivo, devido à mudança na sua estrutura, como por exemplo, um aumento na viscosidade. Tipicamente, são alimentos que contém amido, o qual gelatiniza durante o aquecimento. Como exemplo existe o feijão em salmoura, "cream style corn" e molhos para massas1.

O balanço de energia para o aquecimento por convecção de um alimento enlatado é escrito conforme Equação 4, onde m indica a massa e T¥ a temperatura do meio de aquecimento:

Rearranjando para integração na forma de uma equação de reta e com o tempo como variável independente, temos:

A Equação 5 é válida para fluxo de calor transiente em sistemas com resistência interna desprezível. Por exemplo, alumínio, cobre. Contudo, essa equação pode ser aplicada para alimentos envasados.

BALL e OLSON1 definiram dois parâmetros para descrever a transferência de calor em alimentos enlatados: fator de tempo fh (minutos) e fator de atraso de temperatura jh, ou fator de curvatura. O termo fh é definido conforme a Equação 6 e pode ser obtido pelo tempo necessário para que a curva de log(T¥ – T) em função do tempo atravesse um ciclo logarítmico.

O termo jh é definido como:

Tpih = temperatura pseudo-inicial para aquecimento obtida pelo intercepto da projeção da reta log(T¥ - T) em função do tempo com o eixo das ordenadas.

Enquanto

Os métodos de cálculos de BALL e ONSON1 são baseados na curva de aquecimento e resfriamento semilogarítmica. No eixo das ordenadas a escala é logarítmica e na abscissa a escala é linear.

A relação de fh com o coeficiente de película, h, e com o coeficiente global de transferência de calor, U, está apresentada na Equação 10:

A Equação 10 aproxima o valor de U (coeficiente global de transferência de calor) para o coeficiente de película, h. Dessa forma, nos ensaios de penetração de calor nas embalagens pode ser obtido U admitido como uma estimativa do valor de h.

De modo análogo ao aquecimento, têm-se as mesmas relações para a fase de resfriamento do alimento.

O cálculo de h e Nusselt (Nu) a partir de fh pode ser descrito conforme a Equação 11, onde V é o volume e D o diâmetro da lata.

QUAST e SIOSAWA14 mostraram que a rotação axial de latas contendo solução de glicose aumenta a transferência de calor. O tempo de aquecimento com latas a 150 rpm foi 30 a 60% menor que o tempo de aquecimento sem rotação, quando utilizados pasteurizadores a vapor.

QUAST et al.15 estudaram os efeitos do meio de resfriamento, da velocidade de rotação aparente, da posição da lata (horizontal, vertical e inclinada) da consistência e do espaço livre sobre a velocidade de resfriamento de soluções de carboximetilcelulose e de polpa de frutas em latas de 20 kg. O coeficiente global de transferência de calor (U) aumentou à medida que se utilizou espaço livre maior. A posição vertical foi a que apresentou maior eficiência na transferência de calor, devido à maior convecção do fluido. A posição horizontal foi a mais vantajosa em termos do projeto do equipamento. A consistência do fluido influi bastante na transferência de calor: fluidos mais viscosos (polpa de goiaba 12 °Brix) apresentaram maior tempo de resfriamento que fluidos pouco viscosos (2,5% de carboximetilcelulose).

TEIXEIRA NETO et al.19 estudaram o resfriamento de polpas concentradas de tomate e goiaba, previamente aquecida a 85~93 °C, em tambores de 200 L. O envase a quente e resfriamento de polpas de tomate em tambores foi muito utilizado pelas indústrias tomateira para estocagem de seus produtos na entressafra. A operação de processo consistia em encher o tambor com polpa, mantendo espaço livre e promover o resfriamento da massa através da rotação do tambor (24~36 rpm) sob pulverização de água potável e fria. O aumento da rotação aplicada ao tambor melhorava a transferência de calor, sendo mais aconselhável a rotação de 36 a 43 rpm. O resfriamento de 190 kg de polpa (com 10 cm de espaço livre) apresentou uma redução de tempo de 25% em relação ao tempo gasto para o resfriamento de 200 kg (5 cm de espaço livre).

CLIFCORN et al.4 estudaram a redução do tempo de atraso térmico, em função da posição das latas numa autoclave. O estudo foi realizado numa autoclave adaptada para ALI. Foram analisadas várias posições das latas e movimentos com variação de amplitude e freqüência. Esse método de agitação mostrou-se superior quando comparado com método de agitação longitudinal axial, porém, inferior à agitação "end-over-end". Nesse experimento, a situação otimizada foi de 50 mm de amplitude e 120 cpm (ciclos por minuto) de freqüência de agitação.

O estudo de transferência de calor em alimentos líquidos e não sólidos depende, em primeira instância, do estudo da reologia do alimento. Nesse sentido, é comum utilizar fluidos simulantes de alimentos e associar suas características às de alimentos, conhecendo-se suas propriedades físicas através de modelos matemáticos.

A análise dimensional permite deduzir correlações, transcritas em equações constituídas por grupos adimensionais de variáveis elevados a diversas potências6. O método de análise dimensional mais adequada para estudos de processos que envolvem transferência de calor é expressa por:

onde K, a, b, c, d, e e f são constantes adimensionais determinadas experimentalmente. Contudo, qualquer outra expressão algébrica ou mesmo relação gráfica, entre esses grupos de números adimensionais, pode ser obtida de dados experimentais12.

O número de Nusselt, , representa o coeficiente de transmissão de calor no aquecimento por convecção. Indica a porção entre o calor fornecido por convecção pura e aquele fornecido por condução pela película fluida próxima à parede da embalagem6.

O número de Prandtl, , é a relação entre a difusividade da quantidade de movimento e a difusividade térmica do fluido. Esse adimensional é utilizado na análise de transferência de calor por convecção6.

A ALI estudada por CLIFCORN et al.4 em 1950, ao que se sabe, não teve nenhuma aplicação prática na indústria de alimentos. A princípio, na época deste estudo, as principais preocupações dos processadores eram grandes escalas de produção, sem preocupação com o elevado consumo de energia nos ciclos de processamento. Isso levou ao desenvolvimento de sistemas de agitação contínuos do tipo axial em batelada FMC Orbitort® e contínuo FMC Sterilmatic® ou "end over end" em batelada, como a STOCK Rotomat®, todas de grande porte. Atualmente, face ao desenvolvimento de novos mercados, processos em menor escala, a custos mais baixos e com qualidade, a ALI poderá suprir esta necessidade com vantagens operacionais, de manutenção e homogeneidade no lote agitado.

Objetivou–se estudar a transferência de calor transiente por ALI em embalagens metálicas, pois este processo é de construção e manutenção relativamente simples se comparado ao sistema de agitação tradicional "end-over-end" (cabeça-sobre-cabeça); e a agitação axial, além de promover a agitação homogênea de todo o lote, pois não leva em consideração o raio de agitação. Para tal foi desenvolvido um sistema de agitação com movimento linear alternado que pode ser adaptado em pasteurizadores ou autoclaves fixas, tornando-os agitados. Também foram desenvolvidas as equações de transferência de calor, baseadas nos parâmetros de processo, com o intuito de se projetar processos térmicos a partir destas informações.

2 Material e métodos

2.1 Suporte para agitação

O suporte para agitação foi construído em aço carbono podendo receber no máximo quatro latas. As dimensões do suporte foram as seguintes: base para as latas 340 x 240 mm, altura 400 mm e haste para movimentação de 400 mm. Para efetuar a agitação, a haste do suporte foi conectada a um dispositivo de agitação constituído por um motor 1/4 CV trifásico, polias e correias. A amplitude de agitação foi de 25 mm (mínimo) a 42 mm (máximo).

2.2 Meio de aquecimento e resfriamento

Foi composto de dois tanques paralelos, sendo um para aquecimento e outro para resfriamento. Cada tanque tinha 82 litros (500 x 470 x 350 mm altura) de capacidade nominal. O aquecimento foi feito em água em ebulição aquecida por vapor direto à pressão de 2 kgf.cm–2. O resfriamento foi feito por circulação de água (6 litros/minuto) à temperatura ambiente (25 ± 1 °C). Ao final do aquecimento o suporte com as latas era transferido manualmente para o tanque de resfriamento.

2.3 Meios simuladores

Foram utilizados 4 meios para simulação: água destilada e óleos lubrificantes para uso automotivo, das marcas, Neutral Oil 150, Neutral Oil 500 e Bright Stock Oil.

2.4 Características físicas da lata experimental

Lata construída em folha de flandres, com verniz anticorrosivo, com anéis de expansão, 81 mm interno de altura x 73 mm de diâmetro interno.

2.5 Enchimento e recravação

As latas foram preenchidas com fluido simulante, T = 25 °C, mantendo-se espaço livre (5, 10 ou 15% do volume da lata), conforme o processo e foram recravadas sob vácuo de 600 mm Hg (23,5 polegadas Hg) ± 20 mm em recravadeira marca John Heine, Sydney, model 71D, series 2, previamente ajustada conforme normas do FDA (INSTITUTO DE TECNOLOGIA DE ALIMENTOS, 2001). Foram utilizadas 4 latas (capacidade do suporte), sendo uma em pé, uma deitada longitudinalmente ao eixo de movimento, outra transversalmente ao eixo de movimento e finalmente uma a 45° com o eixo do movimento. A lata em pé recebeu um suporte de termopar de 52 mm, com sensor colocado a 1/3 da altura da lata, por se tratar de um fluido convectivo (simulador de alimento). As outras três latas, deitadas, foram montadas com suporte de termopar de 42 mm, posicionado no centro da lata. As latas e o termopar, para medida da temperatura do meio de aquecimento, foram fixados no suporte com braçadeiras, diminuindo as chances de ocorrerem vazamentos e quebra dos cabos.

2.6 Termopares e aquisitor de dados

Para monitoramento das temperaturas de processos, foram utilizados 5 cabos com termopares (04 nas latas e 01 no meio de aquecimento), tipo agulha, de junta cobre-constantã, tipo T, um aquisitor de dados e software E-ValTM Ver 2.00 ELLAB A/S Krondalvej 9, DK-2610 Roedovre, Denmark, modelo TM 9616. Os históricos de temperaturas foram levantados por meio do aquisitor Ellab A/S conectado a um microcomputador compatível com PC. Todo o sistema de medida de temperatura, original Ellab (termopares, fios de extensão, aquisitor de dados) foi devidamente aterrado e previamente calibrado contra padrão de confiança.

Medidas de freqüência de agitação e comprimento de curso

Foi utilizado um gerador de freqüência marca Siemens, modelo Micromaster Vector 6SE3221–0CC40, com capacidade para motor de 3,0 CV. Esse gerador permitiu gerar freqüências de 0,1 até 65,0 Hz, com variação decimal. O gerador foi utilizado para acionar o motor trifásico instalado no suporte das latas, gerando as freqüências de agitações estudadas. A freqüência de agitação foi medida toda vez que iniciava-se o experimento.

2.7 Propriedades de engenharia dos fluidos

Foram determinados os valores da massa específica, viscosidade, calor específico, condutividade térmica e difusividade térmica dos fluidos conforme metodologias descritas em GRATZEK e TOLEDO7, KERN9, KESSLER10 e RAO et al.16. Para a água, estes valores foram obtidos teoricamente10.

2.8 Cálculo dos parâmetros de transferência térmica

O cálculo dos parâmetros fh e jh para aquecimento e fc e jc para resfriamento foi baseado em STUMBO18. Inicialmente foi feita uma avaliação gráfica das curvas de penetração de calor de cada lata, nas posições estudadas, selecionando-se os pontos da fase de aquecimento, log (T¥ – T) em função do tempo (segundo), e da fase de resfriamento, log (T – T¥) em função do tempo (segundo). Foram selecionados no mínimo 10 pontos para cada caso, escolhendo-se a porção retilínea das curvas com o maior resultado do R2 (coeficiente de determinação múltipla). Em cada condição estudada foi realizada pelo menos uma repetição e, no máximo, n repetições conforme procedimento descrito por GUMERATO e SCHMIDT8.

O valor numérico de fh é o inverso do coeficiente angular da reta e jh foi calculado conforme Equação 13.

O valor numérico de fc é o inverso do coeficiente angular da reta e jc foi calculado conforme Equação 14.

2.9 Planejamento experimental e avaliação do modelo matemático de transferência de calor

Foi utilizado um planejamento fatorial completo com as variáveis independentes: fluido (FL), espaço livre (EL), amplitude (AM), freqüência de agitação (FA) e posição (PO). Foram escolhidos 4 fluidos, estáveis a vários processos de aquecimento e resfriamento, sem sofrerem deteriorações microbiológicas e físico-químicas, cobrindo ampla faixa de viscosidade. Foram definidos 3 espaços livres (EL) comuns em aplicações industriais20: 5, 10 e 15% do volume total. Foram aplicadas AM de 25 e 42 mm e FA de 0, 60, 100 e 132 cpm. As 4 PO das latas (vertical, transversal, inclinada e longitudinal) foram planejadas para cobrir a faixa de disposição espacial que poderia ser influenciada pelo movimento de ALI. As condições de agitação foram baseadas em testes preliminares (não apresentados neste relato) que apresentaram garantia mecânica do funcionamento do sistema de agitação.

Na construção do modelo, utilizaram-se as variáveis de categoria (FL e PO) e as variáveis contínuas (EL, AM e FA), tendo o valor de fh e fc como variáveis dependentes. Utilizou-se a distribuição normal, e o modelo de predição dos resultados foi descrito conforme a função:

Para análise dos dados resultantes dos testes foi utilizada a regressão linear múltipla. As comparações dos parâmetros de transferência de calor entre as posições das latas foram feitas utilizando-se modelo linear geral2.

3. Resultados e discussão

3.1 Propriedades de engenharia dos fluidos simulantes

O modelo que melhor descreveu a dependência da massa específica dos óleos simuladores com a temperatura foi o modelo linear. Na Tabela 1, estão listados os parâmetros da regressão linear dos valores de massas específicas em função da temperatura. Para os óleos orgânicos, o modelo aplicado foi r = a.t + ro, onde ro é a intersecção da reta quando a temperatura é igual a zero.

No caso da água, a relação entre a temperatura e a massa específica foi quadrática e foi aplicado o modelo matemático r = 0,569 + 0,002 x T – 5,87 x 10–6 x T2 com os dados listados por10.

Quanto aos valores de calor específico, foram linearmente crescentes para os óleos, na faixa de temperatura de trabalho de 20 a 100 °C, o que já foi descrito anteriormente9. Os resultados em (kJ.kg–1.K–1) foram: água destilada, 4,18410; neutral oil 150, 1,859 ± 4,3%; neutral oil 500, 1,850 ± 7,4% e bright stock oil, 2,015 ± 7,3%.

Em relação aos valores da condutividade térmica para a água, foi realizada uma regressão linear, conforme10, utilizando um nível de significância de 10%. Para os três óleos, os resultados foram iguais a 0,147 W.m–1.K–1.

As viscosidades dos óleos foram quantificadas conforme metodologia descrita, utilizando-se um nível de significância de 1%. No caso da água, os resultados foram coletados e fez–se regressão linear10. Os parâmetros para modelo de Arrhenius estão apresentados na Tabela 2. Esses modelos foram avaliados através de regressão linear entre os valores de viscosidades e o inverso das temperaturas absolutas experimentadas. O modelo descreveu muito bem o comportamento de óleos e água, na faixa de temperatura deste trabalho.

3.2 Medidas da freqüência de oscilação

A função que relaciona a agitação (cpm) com a freqüência da tensão (sempre para o mesmo conjunto de polias) foi linear, apresentou R2 = 99,9%, e foi utilizada para estimar a velocidade da agitação.

3.3 Resultados dos parâmetros fh, jh, fc e jc obtidos nos ensaios

Foram feitos 294 ensaios de penetração de calor obtendo–se o mesmo número de valores para cada um dos parâmetros de transferência de calor (fh, jh, fc e jc). Na Figura 1, está representada a história térmica de um dos experimentos realizados (amplitude 42 mm; freqüência de agitação 120 cpm; espaço-livre 8,1 mm). A história térmica são os perfis de temperatura de cada uma das 4 latas e mais um termopar de medida da temperatura do meio em função do tempo de processo.

Os valores apresentados nas Figuras 2 a 5 são resultados de desvios respectivamente de fh, jh, fc e jc em função das médias (em %) de no mínimo duas observações para cada condição de ensaio. Estas figuras foram construídas para avaliar a dispersão e, conseqüentemente, a qualidade dos dados coletados nos experimentos, pois não existe modelo matemático descrito para esse mecanismo de transferência de calor.

Na Figura 2, pode ser observado que a distribuição dos desvios experimentais (em porcentagem) em função dos valores de fh (minutos) foi de no máximo 22% e de no mínimo 0%. Observa–se que não há nenhuma tendência entre os desvios com os valores de fh. De forma análoga, na Figura 3, apresenta-se a distribuição dos desvios experimentais (em porcentagem) em função dos valores de fc (minutos). O desvio máximo foi de 26% e o mínimo de 0%.

Os fatores de curvatura, jh e jc, foram calculados considerando-se todos os fluidos estudados e o valor médio para jh foi de 1,4 e jc 1,0. BALL e OLSON1 recomendam jc, fator de curvatura, de 1,4 para alimentos condutivos e 1,0 para alimentos convectivos. O Bright stock oil proporcionou uma curvatura (jh) maior no aquecimento, quando sob alta agitação (132 cpm) e amplitude de 42 mm. Os demais fluidos simuladores apresentaram valores menores e próximos de 1,4.

Na Figura 4, pode ser observada a distribuição dos desvios experimentais (em porcentagem) em função dos valores de jh (minutos). A maioria dos resultados de jh estão concentrados em 1,4, com desvio máximo de 26% e mínimo de 0%. Os fatores de curvilinearidade são dependentes do raio da lata, segundo BALL e ONSON1. Dois valores de jh estão fora da faixa de 1,4 e não foram explicados (3 e 6,5).

A Figura 5 apresenta a distribuição dos desvios experimentais (em porcentagem) em função dos valores de jc (minutos). O desvio máximo foi de 24% e o mínimo de 0%. Observou-se também que a maioria dos resultados de jc está concentrada em 1,4. Esses resultados estão de acordo com1.

3.4 Modelamento matemático do processo de transferência de calor no aquecimento para latas na posição vertical.

Foram analisados vários modelos na tentativa de associar fh (minutos) com medidas reológicas e características do processo. Os cálculos das propriedades dos fluidos (µ–viscosidade; r–massa específica; Cp–calor específico; e k–condutividade térmica) para compor os adimensionais Re, Nu e Pr foram feitos utilizando-se a temperatura média de processo de 60 °C (medida entre o aquecimento e o resfriamento). Os fluidos utilizados neste trabalho não sofreram alterações em suas propriedades de engenharia, as quais se assemelham à de vários alimentos fluidos e também não foram observadas curvas quebradas durante a fase de aquecimento e resfriamento. Os cálculos de fh e fc foram determinados por regressão linear, conforme etapas descritas no item 4.8.

Avaliou-se o efeito da posição das latas sobre o fh e concluiu-se que não teve influência significativa de aproximadamente 41% (1 – 0,595686) nos resultados. A posição vertical é a posição de melhor manejo nas indústrias de alimentos, facilitando o processo de carga e descarga mecanizadas de latas nos cestos de produção, diminuindo os danos mecânicos na litografia e evitando amassamentos em geral. Por esse motivo, na discussão de resultados subseqüentes apenas a posição vertical será considerada.

Da mesma forma, a combinação posição e freqüência foi pouco significativa. Por outro lado, os efeitos espaço-livre, freqüência de agitação, amplitude e fluido foram significativos. A amplitude foi um fator relevante na transferência de calor (fh), sendo que quanto maior, maior o valor de fh.

O modelo matemático foi construído baseado na relação de números adimensionais e em funções trigonométricas. No cálculo do Re foi aplicada a equação a seguir:

onde: xa = amplitude de agitação (mínimo = 1 m); xf = freqüência de agitação (mínimo = 1 cpm); e EL = espaço livre na lata (mm); para adequar o modelo no limite mínimo foi necessário atribuir amplitude = 1 mm e freqüência = 1 cpm para processo estático.

O valor de Re foi modificado para se adequar ao processo de agitação e, dessa forma, foi incluído o diâmetro da lata e a velocidade de agitação nos termos xa e xf, que são a proporção de amplitude e a proporção da freqüência de agitação, respectivamente. O número de Reynolds modificado variou de 8,7 x 10–3 a 2,4 x 104 e quantificou o comportamento do escoamento (teórico) do fluido na lata na posição vertical.

O número de Prandtl, calculado pela Equação 18, foi mantido sem alterações e variou de no mínimo 3,6 a no máximo 3,4 x 103, para as latas na posição vertical.

A Equação adimensional 18 foi então desenvolvida para as latas na posição vertical, utilizando-se, além do número de Reynolds e Prandtl, as funções trigonométricas seno e co-seno. Essas funções foram utilizadas para transformar medidas discretas (freqüência e amplitude de agitação) em uma função contínua (Nu). A aplicação dessas transformações foi inspirada na física quântica (funções de níveis quânticos)5.

Onde: xa/AM = razão entre amplitude de agitação (mínimo = 1) e amplitude máxima de agitação e xf/FA = razão entre freqüência de agitação (mínimo = 1) e freqüência máxima de agitação. H indica a altura da lata.

Os valores do coeficiente de película, h, variaram de 9,7 x 101 a 2,3 × 103 W.K.m–2. Os valores encontrados foram respectivamente aqueles observados para latas mantidas estáticas e sob agitação, na posição vertical. O número de Nusselt foi calculado utilizando-se os valores de coeficiente de película, h, e variou de 5,0 x 10–2 a 2,6 x 10–1.

Os resultados das constantes do modelo estão descritos na Tabela 3. A primeira coluna refere–se ao coeficiente da Equação do modelo, a segunda coluna, seu valor numérico e a terceira à precisão do resultado. Foi considerado nível de significância de 5%.

Na Figura 6, estão representados os valores dos resíduos de fh (fhexperimental – fhcalculado)/fhexperimental em função de fhexperimental. O resíduo mínimo foi de –50% e o máximo foi de +45%. Os resultados de erros maiores estão na faixa de 2 a 5 minutos de fh. Nessa faixa, pode–se observar a distribuição dos pontos em forma de três retas inclinadas para a direita. Elas são efeitos residuais do modelamento matemático e ocorreram com os óleos Bright stock oil e Neutral oil 150 e 500.

Modelamento matemático do processo de transferência de calor no resfriamento para latas na posição vertical

Aplicando as Equações 16 e 17 foi desenvolvida para o caso do resfriamento, de modo análogo ao aquecimento, a Equação 19:

Os resultados das constantes do modelo estão descritos na Tabela 4. A primeira coluna refere–se ao coeficiente da equação do modelo, a segunda coluna, a seu valor numérico e a terceira, à precisão do resultado. Foi considerado nível de significância de 5%.

Na Figura 7, estão representados os valores dos resíduos de fc (fcexperimental – fccalculado)/fcexperimental em função de fcexperimental. O resíduo mínimo foi –35% e o máximo foi 58%. Esse gráfico mostra os desvios das medidas de fc quando aplicado o modelo conforme Equação 19. Observa–se que os resíduos são bem distribuídos e que ainda persiste uma "reta" de inclinação positiva que descreve o fluido água. Para o resfriamento foi seguido o mesmo raciocínio que no aquecimento. Os valores de fc menores que 5 minutos apresentaram maior desvio, quando comparado com fc maiores que 5. O ajuste foi melhor na faixa de fc maior que 5 minutos.

Observações quanto aos modelos matemáticos propostos, tanto para aquecimento quanto para resfriamento

Observou-se que a dispersão dos dados aumentou à medida que o valor numérico de Nusselt também cresceu. A maior dispersão ocorreu com os dados obtidos com água, ou seja, no resfriamento e a maior distorção (diferença entre o experimental e o calculado) ocorreu com o Bright stock oil.

Foram calculados os mesmos coeficientes em outras duas temperaturas: 0 e 100 °C e não foi observada diferença estatística a um nível de significância de 5% entre os coeficientes adimensionais. Nos resultados apresentados neste trabalho para validação do modelo escolheu-se a temperatura média entre o final do aquecimento e o resfriamento, que foi de aproximadamente 60 °C.

4 Conclusões

Não houve diferença significativa (1% de significância) entre os valores de fh entre as posições das latas, para o mesmo fluido. Isso resultou no processamento das latas na posição vertical, mais favorável no ambiente industrial, com separadores dentro dos cestos que forem mecanizados para agitação.

Esse processo de agitação é homogêneo, pois todas as latas passam pela mesma intensidade de agitação e garante o mesmo fh entre as latas de todas as posições distribuídas pelos cestos.

O parâmetro de transferência de calor fh em latas agitadas pode ser estimado com auxílio da Equação 18 com desvio máximo de 40% para valores de fh entre 3 e 5 minutos.

Quanto maior a freqüência de agitação das latas, menor foi o valor de fh e, portanto, maior a transferência de calor para o alimento. O limite para o aumento da freqüência depende da construção mecânica da autoclave e um valor de 120 cpm pode ser considerado satisfatório, e está de acordo com CLIFCORN4. Não foram investigadas freqüências maiores que 132 cpm.

As Equações 18 e 19 permitiram relacionar amplitude, freqüência de agitação, massa específica, calor específico, condutividade térmica e viscosidade com fh e fc para o processo de agitação com movimento linear alternado.


Link retirado de:

http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0101-20612007000300034&lng=pt&nrm=iso

Nenhum comentário: