Integração energética da rede de trocadores de calor em extração por solvente para a produção de farelo branco de soja

A síntese de redes de trocadores de calor foi muito bem estudada pela comunidade científica nos últimos trinta anos, porém, o desenvolvimento de métodos e a melhoria dos já existentes sempre foram o foco principal da grande maioria dos autores. São poucos os trabalhos que aplicam qualquer uma das metodologias a situações reais e discorrem sobre as dificuldades encontradas - provavelmente porque as empresas não têm interesse em tornar público os resultados obtidos. Este trabalho aplica a metodologia do ponto de estrangulamento na nova fábrica de farelo branco do Parque Industrial de Esteio e, portanto, pode servir como um guia prático aos profissionais que desejem fazer o mesmo em outras plantas industriais. Esta dissertação sintetiza duas redes de trocadores de calor para a fábrica de farelo branco de Esteio, levando em consideração as condições econômicas e os custos de equipamento e de montagem no Brasil. A primeira rede estabelecida é baseada em aproveitamentos térmicos utilizados pelos fabricantes e pelas empresas que dominam a tecnologia de extração de óleos vegetais, já a segunda rede sintetizada utiliza a metodologia do ponto de estrangulamento. Os dois resultados obtidos são bastante viáveis economicamente, sendo que a rede sintetizada pelo ponto de estrangulamento é mais econômica. O Valor Presente Líquido é utilizado como critério de cálculo de viabilidade das duas redes. Foi demonstrado que para a realidade brasileira este critério é mais adequado que o Custo Total Anualizado. A flexibilidade da rede produzida pelo método do ponto de estrangulamento foi avaliada de maneira preliminar por simulação e tabelas de sensibilidade. O resultado obtido é robusto na posta em marcha e quando a planta é submetida a variações no extrator - situações que produzem distúrbios importantes e freqüentes. A aplicação da metodologia do ponto de estrangulamento na integração energética da fábrica de farelo branco confirmou a simplicidade do método e a sua forte interação com o engenheiro de processos. O Rio Grande do Sul tem muito a ganhar se esta metodologia for disseminada nas universidades e nas empresas que compõe o parque industrial do estado.

Características dinâmicas do escoamento turbulento em bancos de tubos inclinados

Bancos de tubos ou de barras são encontrados nas indústrias nucleares e de processos, constituindo-se nos dispositivos mais comumente utilizados no desenvolvimento de trocadores de calor. O procedimento usual na análise de bancos de tubos considera o escoamento incidindo uniformemente no banco ao longo de toda a extensão dos tubos. Em aplicações reais de trocadores de calor este fato não é observado, pois o escoamento não é uniforme devido a presença de defletores e condições de entrada e saída. A proposta deste trabalho é apresentar um estudo experimental das flutuações de pressão e velocidade e sua interdependência do escoamento turbulento incidindo em arranjos de cilindros inclinados, com diferentes ângulos de incidência (30º, 45º, 60º e 90º) simulando bancos de tubos inclinados com arranjo quadrangular e razão de espaçamento P/D=1,26. A seção de testes é um canal retangular com 146 mm de altura e largura de 193 mm. O fluido de trabalho é o ar, impulsionado por um ventilador centrífugo, passando por uma câmara de tranquilização e um conjunto de grades e telas, que tem a finalidade de uniformizar o escoamento e reduzir a intensidade de turbulência para a ordem de 2%, antes de alcançar a seção onde são realizadas as medições. Medições da distribuição de pressão na parede e flutuações de velocidade e pressão em torno de um dos tubos no interior dos bancos foram realizadas para diferentes ângulos de incidência, utilizando a técnica de anemometria de fio quente e um transdutor de pressão. Os resultados experimentais são apresentados na forma adimensional A distribuição média de pressão em uma das paredes laterais do canal, com a presença do banco de tubos, mostrou que as linhas isobáricas não são paralelas ao eixo dos tubos, indicando que o gradiente de pressão não é perpendicular aos tubos, nem paralelo a direção principal do escoamento. As médias quadráticas das flutuações de pressão indicam uma distribuição mais uniforme das cargas dinâmicas em torno dos tubos do bancos com inclinações de 30º e 90º. A análise dos resultados dos espectros das flutuações de pressão mostra que conforme o ângulo de incidência do escoamento aumenta, picos aparecem para freqüência adimensionais cada vez menores. O comportamento dos espectros de flutuação de velocidade mostram um decaimento suave das curvas quando o ângulo de incidência do escoamento aumenta. A correlação cruzada entre as flutuações de velocidade e pressão indicam que, quando o ângulo de incidência do escoamento diminui, importantes efeitos tridimensionais aparecem. Esses efeitos, caracterizados pelas oscilações na correlação cruzada para o ângulo de incidência de 45º, parecem desaparecer se o ângulo de incidência continua-se a ser reduzido.

Estudo experimental dos campos de pressão e velocidade em bancos de tubos com a utilização de defletores

A utilização de placas defletoras em banco de tubos é muito difundida em trocadores de calor, aumentando os níveis de turbulência do escoamento e o tempo de contato do fluido refrigerante e a superfície aquecida dos tubos. Os vórtices gerados pela turbulência tem a finalidade de imprimir a camada limite térmica dos tubos, onde o calor é trocado basicamente por condução, tornando a troca de calor mais eficiente. No entanto a maioria dos estudos em bancos de tubos considera um escoamento incidindo uniformemente sobre os tubos, sem qualquer tipo de análise dos efeitos causados pela utilização de placas defletoras, transformando-se assim em idealizações pouco reais. A proposta deste trabalho é analisar experimentalmente os efeitos causados no escoamento pelo uso de defletores, tomando-se como referência a região entre-tubos (linhas), onde para tanto foram estudadas duas geometrias distintas. A primeira geometria consiste uma placa defletora voltada para o escoamento e nela encontram-se montados dez tubos, dispostos horizontalmente em duas linhas, e a partir de uma certa distância a jusante ocorre um estreitamento súbito do canal prolongando-se até a saída. Com este tipo de geometria, pretende-se estudar os efeitos no desenvolvimento do escoamento para uma entrada assimétrica. A segunda geometria estudada é uma simplificação de bancos de tubos utilizados em trocadores de calor, porém, agora, com o uso de defletores. Nesta seção o escoamento é obrigado a contornar três placas defletoras, duas na parede superior do canal e uma na parede inferior, montadas de forma eqüidistantes O número de Reynolds, para ambas as seções, foi calculado com o diâmetro externo do tubo (32.1 mm) e a velocidade média do escoamento abaixo da primeira placa defletora, Ue. Para a primeira seção o número de Reynolds foi de 3.32 x 104 e para segunda 2.99 x 104. Os testes foram realizados em um canal retangular de 146 mm de altura e 193 mm de largura, com 900 mm de comprimento. Os bancos de tubos possuem uma relação de espaçamento, passo-diâmetro de 1,26, com dez tubos dispostos em duas linhas. As velocidades médias, suas direções e flutuações de velocidade foram investigadas ao longo das fendas, linhas, permitindo um conhecimento do desenvolvimento das componentes de velocidades, bem as intensidades de turbulência desde a região de entrada. Este estudo revelou que o desenvolvimento do escoamento , após a passagem por defletores se dá para distâncias superiores aquelas aqui estudadas.

Avaliação da eficiência energética em um sítio industrial utilizando análise pinch com heurísticas e técnicas de projeto

Este trabalho envolve a avaliação da eficiência energética em um sítio industrial, que contém três plantas de produção de polietilenos, utilizando conceitos da análise pinch associados com heurísticas e técnicas de projetos. Na primeira etapa do trabalho, foi realizada uma análise preliminar, na qual foram verificadas as possibilidades de melhorias nas condições de operação das unidades. Com esta avaliação, obteve-se uma economia de 500kg/h de vapor de baixa pressão com a simples alteração na sua pressão de fornecimento para uma das unidades. Na seqüência, com base no conhecimento dos processos e na avaliação das curvas compostas das unidades, verificou-se o potencial de integração energética das plantas individualmente, selecionando-se uma unidade (Planta 3) para ser reprojetada. Para esta unidade foram eliminadas do problema as correntes utilizadas em condições esporádicas, as correntes fora da faixa de integração e, por fim, aquelas cuja integração seria demasiadamente complexa. Após estas simplificações, verificou-se o Ponto de Estrangulamento Energético da unidade, os trocadores que empregam utilidades de uma forma inapropriada, o consumo atual e o mínimo de vapor. Com base nestas informações foram verificadas as modificações estruturais na rede de trocadores de calor da Unidade 3 que resultariam em uma unidade energeticamente mais eficiente. Para isto, utilizou-se o programa Aspen Pinch (da Aspen Tech) que gerou uma lista com vinte possibilidades de adição de novos trocadores. Considerando-se somente as modificações com economia superior a 200kW e as restrições de segurança, layout e controle, apenas duas oportunidades foram detalhadas. A primeira delas, foi o aproveitamento do vapor de flash resultante da expansão do condensado de média pressão da unidade. A segunda oportunidade envolveu a adição de um novo trocador, que foi projetado utilizando-se o programa Bjac conjuntamente com o Aspen Plus. Ambas as oportunidades foram avaliadas economicamente, apresentando Taxas Interna de Retorno muito superiores à Taxa mínima de atratividade de 15% considerada neste trabalho.

Análise teórico experimental dos evaporadores de uma bomba de calor

Faz-se um estudo da construção de uma bomba de calor residencial de baixo custo utilizando um condicionador de ar do tipo janela. Com tal equipamento operando fora das condições de projeto é preciso fazer uma avaliação das novas condições de operação, onde os evaporadores passam a ser objetivo de estudo do presente trabalho. A análise foi conduzida de duas formas. Uma tomada de dados experimental, feita em duas épocas distintas do ano, nas condições mais severas de consumo, no inverno, e nas condições mais severas de temperatura, no verão. A outra, uma análise teórica, gerando dados que podem ser comparados com a literatura. Inicialmente, é feita uma avaliação do coeficiente de performance visando a competitividade do equipamento no que tange o custo operacional. Como objetivo principal é avaliada a troca de calor nos evaporadores e os parâmetros que a influenciam, bem como o comportamento do fluido refrigerante na entra e saída. Foi preciso desenvolver uma metodologia para trocadores aletados que apresentam mudança de fase já que esta não é vista com facilidade na literatura. Faz-se uso do método NUT e- , que é adaptado de forma a ser função somente dos dados psicrométricos do ar. Para a efetividade média externa e o coeficiente de convecção externo é utilizada a metodologia desenvolvida por Schmidt em 1949 Como resultados, observou-se que o coeficiente de performance fica na faixa dos equipamentos comercias (2,6 na média anual) e o custo operacional é competitivo com o uso do gás liquefeito de petróleo, em aquecedores de passagem. A efetividade do trocador e o coeficiente de convecção externo apresentam um valor baixo, atribuído a grande área de troca de calor externa. O coeficiente de contato, apresenta um valor alto se comparado ao coeficiente de convecção externo. Quanto ao fluido refrigerante, este apresenta condições normais nas saídas dos evaporadores e título baixo na entrada de um dos evaporadores. Também é observado que um dos evaporadores apresenta baixo uso efetivo, sendo o problema atribuído a baixa vazão de refrigerante.

Estudo do calor específico de compostos Heusler paramagnéticos, da série do níquel Ni 2 Tal, onde T=Ti, Zr, Hf, V, Nb e Ta

Neste trabalho apresentamos os resultados experimentais da medida do calor específico dos compostos Heusler da série Ni2TA1, onde T =Ti, Zr,Hf,V,Nb, Ta. Estas medidas foram feitas utilizando-se um calorímetro adiabático que atua na faixa de 1,84 K a 10,3 K e visaram o estudo da estrutura elêtrônica dos compostos em termos da densidade de estados eletrônicos à nível de Fermi, bem como da dinâmica da rede cristalina, com o uso dos modelos de Debye e de Einstein para o calor específico. Apresentamos, também, uma descrição do calorímetro, seu funcionamento e as alterações ocorridas em função deste trabalho.

Simulação numérica de um sistema de armazenamento de calor em estufas de plasticultura

Apresentar um modelo para simular um sistema de armazenamento de calor no solo em estufas para plasticultura é o objetivo do presente trabalho. O sistema consiste num feixe de tubos enterrados no solo. A convecção forçada de ar no seu interior realiza a troca térmica necessária para manter as estufas sob faixas desejadas de temperatura. O objetivo do modelo é investigar os efeitos no calor armazenado e a influência das variáveis, tais como diâmetro, comprimento, espaçamento entre os tubos e a velocidade de ar no canal provocam no sistema. O solo é tratado como um meio difusivo e avalia-se a contribuição do termo de condensação e evaporação da água contida no ar em escoamento nos tubos. A equação da energia é resolvida para o solo e para o ar. Os tubos de seção transversal circular são modelados como tubos de seção transversal quadrada com o objetivo de que as simulações possam ser processadas em coordenadas cartesianas. O programa resolve situações tridimensionais, transientes e emprega o Método dos Volumes Finitos para integrar as equações diferenciais governantes. O modelo original é baseado no modelo de Gauthier et al., 1997, tendo sido os resultados do mesmo foram usados para a validação do presente estudo. Um circuito de água quente é também projetado e apresentado para o aquecimento das estufas. A água circula através de mangueiras sobre o solo e é aquecida por um sistema de queimadores a gás liqüefeito de petróleo ou óleo combustível, transferindo assim calor para o interior da mesma. O projeto de aquecimento foi realizado através de um programa de parceria entre a Ufrgs, Sebrae, Fapergs e a Agropecuária Clarice.

Uso da função de transferência em problemas de condução do calor com a lei de Fourier modificada

Este trabalho visa o uso da função de transferência, a qual relaciona distribuição de temperatura e fluxo de calor, na comparação no domínio freqüência, entre o modelo de difusão usual (parabólico) e um modelo ondulatório (hiperbólico) que inclue o efeito de propagação do calor não instantâneo, sendo avaliados os casos de meio semi-infinito e finito. Para o caso de meio semi-infinito, são determinadas as expressões para as características de amplitude e de fase, considerando tanto a abordagem parabólica quanto hiperbólica. É observada a relação entre estas duas abordagens, mostrando que a abordagem parabólica é uma caso particular da abordagem hiperbólica, podendo ser obtida através de um processo de limite envolvendo o tempo de relaxação r . Para o caso de meio finito, são determinadas as expressões para as caracteríısticas de amplitude de ambas as faces da placa unidimensional, considerando tanto a abordagem parabólica quanto hiperbólica. Estas expressões são transformadas para a forma adimensional quando então são deduzidas as expressões correspondentes das características de amplitude. Mais uma vez, todos os resultados para o caso parabólico podem ser determinados a partir dos resultados do caso hiperbólico, através de um processo de limite envolvendo o tempo de relaxação r São apresentados resultados numéricos referentes às características de amplitude, onde é apontada a existência de uma freqüência limite, acima da qual a diferença entre os dois modelos, do tipo parabólico ou hiperbólico, aumenta rapidamente. Também é apresentada uma forma alternativa de cálculo da distribuição de temperatura transiente que faz uso da função de transferência do sistema.

Solução da equação de condução de calor bidimensional, em meios multicompostos, pelos métodos nodal, com parâmetros concentrados, e a técnica da transformada de Laplace

Neste trabalho, desenvolvemos uma metodologia semi-analítica para solução de problemas de condução de calor bidimensional, não-estacionária em meios multicompostos. Esta metodologia combina os métodos nodal, com parâmetros concentrados, e a técnica da transformada de Laplace. Inicialmente, aplicamos o método nodal. Nele, a equação diferencial parcial que descreve o problema é integrada, transversalmente, em relação a uma das variáveis espaciais. Em seguida, é utilizado o método de parâmetros concentrados, onde a distribuição de temperatura nos contornos superior e inferior é substituída pelo seu valor médio. Os problemas diferenciais unidimensionais resultantes são então resolvidos com o uso da técnica da transformada de Laplace, cuja inversão é avaliada numericamente. O método proposto é usado na solução do problema de condução de calor, em paredes de edificações. A implementação computacional é feita, utilizando-se a linguagem FORTRAN e os resultados numéricos obtidos são comparados com os disponíveis na literatura.

Modelos cinéticos da equação linearizada de Boltzmann e um problema de transferência de calor em microescala

Neste trabalho, um problema de transferência de calor da dinâmica de gases rarefeitos, causado pela diferença de temperaturas nas superfícies de um canal, é abordado. O problema é formulado através dos modelos cinéticos BGK, S e Gross-Jackson da equação linearizada de Boltzmann e resolvido, de forma unificada, pelo método analítico de ordenadas discretas (método ADO). Resultados numéricos para as perturbações de densidade e temperatura e também para o fluxo de calor são apresentados e comparados, mostrando que não se pode dizer que algum dos três modelos seja uma melhor aproximação da solução aos resultados da equação linearizada de Boltzmann.

Influência dos fotoativadores PAC e QTH na polimerização de resina composta e na geração de calor

O objetivo deste estudo foi avaliar e comparar in vitro o desempenho dos fotoativadores Apollo® 95E /DenMed Technologies, com luz xênon (arco de plasma - PAC), e o Curing Light® XL 3000 /3M-Unitek, com luz halógena (convencional - QTH), comparando-se, desta forma, um fotoativador de alta intensidade de luz (PAC) com um convencional (QTH). A efetividade dos aparelhos foi verificada pelo grau de polimerização apresentado pela resina composta (Transbond™XT/3MUnitek), avaliado em espectrofotometria por infravermelho e pela análise de microdureza Knoop. Os tempos de ativação para a fonte PAC foram um, dois e três segundos e, para a fonte QTH, foram 10, 20 e 30 segundos, buscando-se avaliar a equivalência, ou não, dos diferentes tempos das diferentes unidades fotoativadoras. Foram usadas matrizes metálicas circulares, que receberam a resina para posterior polimerização. Para cada tempo de ativação, para cada aparelho e para cada teste, foi fotoativada a resina composta de cinco corpos-de-prova. Foi avaliada, também, a temperatura no extremo distal do condutor de luz de cada aparelho, nos seus respectivos tempos predeterminados, colocando-se justaposto, à ponta do fotoativador, o sensor do dispositivo de aferição Sekundenthermometer/GTH 1160 – Phywe. Os resultados obtidos foram submetidos ao teste não-paramétrico Kolmogorov-Smirnov, à Análise de Variância (ANOVA) e ao teste de comparações múltiplas de Tukey. Os resultados da espectrofotometria apontaram diferença estatística significativa entre os grupos estudados. O grupo QTH alcançou as maiores médias e diferiu do grupo PAC (p=0,01). O ensaio de microdureza Knoop também apresentou resultados que apontaram a existência de diferença estatística significativa entre os grupos estudados, tendo o grupo QTH alcançado as maiores médias (p=0,01). A avaliação da temperatura igualmente apresentou diferença estatística significativa entre os grupos estudados (p=0,01), sendo que o grupo QTH atingiu as maiores médias. Com base na metodologia empregada e na análise dos dados obtidos no presente estudo, pode-se concluir que o grau de polimerização da resina composta e a temperatura no extremo distal do condutor de luz apresentaram diferença estatística significativa com o uso das unidades fotoativadoras Apollo® 95E (xênon) e Curing Light® XL 3000 (halógena) nos respectivos tempos predeterminados (p=0,01), ficando sempre, com as maiores médias, o aparelho convencional com luz halógena.

Métodos multigrid paralelos em malhas não estruturadas aplicados à simulação de problemas de dinâmica de fluidos computacional e transferência de calor

Fenômenos naturais, tecnológicos e industriais podem, em geral, ser modelados de modo acurado através de equações diferenciais parciais, definidas sobre domínios contínuos que necessitam ser discretizados para serem resolvidos. Dependendo do esquema de discretização utilizado, pode-se gerar sistemas de equações lineares. Esses sistemas são, de modo geral, esparsos e de grande porte, onde as incógnitas podem ser da ordem de milhares, ou até mesmo de milhões. Levando em consideração essas características, o emprego de métodos iterativos é o mais apropriado para a resolução dos sistemas gerados, devido principalmente a sua potencialidade quanto à otimização de armazenamento e eficiência computacional. Uma forma de incrementar o desempenho dos métodos iterativos é empregar uma técnica multigrid. Multigrid são uma classe de métodos que resolvem eficientemente um grande conjunto de equações algébricas através da aceleração da convergência de métodos iterativos. Considerando que a resolução de sistemas de equações de problemas realísticos pode requerer grande capacidade de processamento e de armazenamento, torna-se imprescindível o uso de ambientes computacionais de alto desempenho. Uma das abordagens encontradas na literatura técnica para a resolução de sistemas de equações em paralelo é aquela que emprega métodos de decomposição de domínio (MDDs). Os MDDs são baseados no particionamento do domínio computacional em subdomínios, de modo que a solução global do problema é obtida pela combinação apropriada das soluções obtidas em cada um dos subdomínios Assim, neste trabalho são disponibilizados diferentes métodos de resolução paralela baseado em decomposição de domínio, utilizando técnicas multigrid para a aceleração da solução de sistemas de equações lineares. Para cada método, são apresentados dois estudos de caso visando a validação das implementações. Os estudos de caso abordados são o problema da difusão de calor e o modelo de hidrodinâmica do modelo UnHIDRA. Os métodos implementados mostraram-se altamente paralelizáveis, apresentando bons ganhos de desempenho. Os métodos multigrid mostraram-se eficiente na aceleração dos métodos iterativos, já que métodos que utilizaram esta técnica apresentaram desempenho superior aos métodos que não utilizaram nenhum método de aceleração.

Tratamento com calor no controle de Botryosphaeria dothidea (Moug.) Ces. & De Not.) em maçãs cv. Fuji

A necessidade de se prolongar o período pós-colheita fez com que a utilização de tratamentos em pós-colheita com o uso de defensivos agrícolas aumentasse nas últimas décadas. No entanto, a preocupação de consumidores quanto a resíduos as deficiências de controle observadas incentivaram a procura por tratamentos alternativos, entre os quais os tratamentos com calor vem sendo testados em frutos de várias espécies. No presente trabalho estudou-se o efeito do tratamento com calor no crescimento do micélio e germinação dos conídios de Botryosphaeria dothidea e, avaliou-se a viabilidade da termoterapia com calor, isolada ou em combinação com carbonato de sódio no controle da podridão branca e na manutenção da qualidade das maçãs ‘Fuji’ em diferentes estádios de maturação, mantidas ou não sob armazenamento refrigerado. Foi objetivo também verificar a expressão de genes que são regulados durante a interação maçã-Botryosphaeria dothidea, após tratamento com calor. Os experimentos foram conduzidos testando-se combinações de diferentes temperaturas e períodos de exposição ao calor sobre o patógeno B. dothidea e sobre frutos inoculados e não inoculados artificialmente com o patógeno. Os resultados obtidos sobre o patógeno mostraram que, independente da concentração de caldo batata-dextrose (BD) utilizada na suspensão, 45% dos conídios de B. dothidea germinaram nos frutos após 24h e 91% após 48h de incubação a 26°C. Constatou-se um aumento significativo da germinação dos conídios quando acrescentado 2,4% de caldo BD a suspensão. Sobre o efeito da temperatura no patógeno, as temperaturas inferiores a 62°C, em qualquer período de exposição não foram capazes de inibir o micélio do patógeno. Os conídios submetidos à temperatura de 58°C somente foram inativados quando o tratamento foi estabelecido por 60s. A termoterapia com calor associada ou não ao carbonato de sódio reduziu a incidência e severidade da doença podridão branca. A severidade da doença foi influenciada pelo estádio de maturação dos frutos verificando pelo aumento no número de lesões/fruto quanto mais avançada a maturação das maçãs. Os resultados obtidos no RT-PCR evidenciaram que os protocolos utilizados não produziram RNA de alta qualidade e desta forma inviabilizaram a análise da expressão diferencial dos genes na interação maçã-B.dothidea.

Simulação numérica da transferência simultânea de energia e umidade através do solo em um sistema trocador-armazenador de calor

O presente trabalho pesquisa o fenômeno da transferência simultânea de calor e umidade em solos insaturados que envolvem dutos enterrados. Estes dutos e o solo envolvente compõem o chamado sistema trocador-armazenador de calor no solo, que é muitas vezes utilizado como fonte complementar de aquecimento em estufas solares agrícolas. O funcionamento do sistema trocador-armazenador de calor é baseado na energia armazenada no solo durante os horários de máxima insolação, sendo que durante a noite parte desta energia é recuperada. Durante o dia o ar quente do meio interno da estufa solar é bombeado para dentro do feixe de tubos enterrados, que devolvem o ar mais frio na outra extremidade. Por outro lado, durante a noite o ar mais frio do meio interno da estufa é bombeado para dentro dos dutos, os quais efetuam troca térmica com o solo envolvente, que neste horário possui a temperatura mais elevada do sistema. O objetivo geral deste trabalho é resolver o problema transiente periódico e tridimensional da transferência simultânea de calor e umidade em solos não-saturados, que compõem o sistema trocador-armazenador de calor, utilizando a simulação numérica. Como objetivos secundários têm-se: o melhoramento do sistema de troca térmica, a quantificação da parcela de calor transportado pela difusão da umidade no solo e a análise dos campos de temperatura e de conteúdo de umidade, sendo que a análise dos campos de umidade permite verificar se existe formação de frentes de secagem significativas na vizinhança dos dutos durante os sucessivos períodos de aquecimento e resfriamento a que o sistema é submetido. Na resolução do problema em questão é empregado o modelo clássico de Philip e De Vries para a transferência simultânea de calor e massa em meios porosos insaturados Neste modelo, as equações de conservação de energia e massa obtidas trazem explicitamente as influências combinadas dos gradientes de temperatura e de conteúdo de umidade nos processos de transporte de calor e umidade. O sistema de equações diferenciais governantes do problema em questão é resolvido numericamente utilizando o Método dos Volumes Finitos e na discretização destas equações é usada uma integração temporal totalmente implícita. Todas as propriedades difusivas e termofísicas empregadas são consideradas variáveis com a temperatura e o conteúdo de umidade. Os dutos de seção transversal circular do sistema trocadorarmazenador de calor no solo, são modelados como dutos de seção transversal quadrada de área equivalente para que coordenadas cartesianas possam ser utilizadas nos modelos analisados. Neste trabalho são simulados quatro modelos computacionais associados ao sistema trocador-armazenador de calor no solo. Estes modelos são compostos por: um duto isolado, um duto com convecção, dois dutos isolados e dois dutos com convecção. A variação da temperatura do ar na entrada do(s) escoamento(s), assim como a temperatura do meio ambiente, para os modelos com convecção, é dada por uma senóide com uma amplitude de 14 ºC. No modelo de um duto isolado, são realizadas simulações utilizando várias combinações dos parâmetros do modelo em questão e os resultados, assim obtidos, são comparados com aqueles encontrados na literatura Visando melhorar o sistema de troca térmica dos modelos computacionais investigados, são selecionados valores e intervalos de valores recomendados para os parâmetros do modelo de um duto isolado. Para este modelo, com um diâmetro de 0,1 m, são escolhidos valores (ou intervalos de valores) recomendados: de 4 m/s para a velocidade do escoamento interno dentro do duto, de 0,25 para o conteúdo de umidade do solo, de 5 até 20 metros para o comprimento do duto e de 0,20 até 0,30 m para a distância entre centros do dutos. As simulações dos quatro modelos computacionais realizadas utilizando as várias combinações dos valores recomendados para os parâmetros destes modelos, mostrou que não há diferença significativa entre os valores de calor volumétrico armazenado no solo empregando a resolução acoplada das equações de energia e de massa e a resolução da equação da temperatura. Mesmo para os modelos de um e de dois dutos com convecção a diferença percentual encontrada foi insignificante. Finalmente, são apresentados e analisados os campos de temperatura e de conteúdo de umidade para os quatro modelos computacionais avaliados. Os perfis de temperatura e de conteúdo de umidade em diferentes horários mostraram que, durante o dia, o solo absorve calor dos escoamentos internos de ar e, uma vez que, junto à superfície dos dutos tem-se regiões de maior temperatura, há, conseqüentemente, uma migração da umidade nestas regiões. Durante a noite, ocorre o contrário, o solo fornece calor aos escoamentos dentro dos dutos, e, desta forma, as regiões próximas aos dutos apresentam níveis de conteúdo de umidade superiores ao inicial. Ainda, os perfis de conteúdo de umidade para todas as situações analisadas mostraram que, não há formação de frentes de secagem significativas nas proximidades dos dutos que compõem os quatro modelos computacionais avaliados.