Preparação e caracterização de nanopartículas de ródio(0) imobilizadas em líquidos iônicos : novo sistema para hidrogenação de aromáticos

O presente trabalho investigou como a natureza de um líquido iônico pode influenciar a formação e o crescimento de nanopartículas de [Rh(0)], bem como estabeleceu o desempenho catalítico das nanopartículas obtidas em reações de hidrogenação de arenos. A simples redução de compostos de ródio RhCl3.3H2O, [Rh(cod)Cl]2 e [Rh(cod)2]X (X = BF4 -, CF3SO3 - e cod = 1,5-ciclooctadieno), dispersas em líquidos iônicos tetrafluoroborato (BF4 -), hexafluorofosfato (PF6 -) e trifluorometanosulfonato (CF3SO3 -) de 1- n-butil-3-metilimidazólio (BMI+), com hidrogênio molecular, rendeu nanopartículas de [Rh(0)], isoladas como pó escuro. As nanopartículas de [Rh(0)] formadas nos três líquidos iônicos foram caracterizadas por Microscopia Eletrônica de Transmissão (MET), Difração de Raios X (DRX) e Espectroscopia de Fotoelétrons de Raios X (XPS). O tamanho médio e a dispersão de tamanho destas nanopartículas depende pouco do precursor de ródio e mais do líquido iônico empregado. O tamanho médio das nanopartículas geradas em BMI.BF4 centram-se em 2,8 nm e em contraste, tamanhos médios maiores foram observados para nanopartículas de [Rh(0)] obtidas em BMI.PF6 (4,7 nm) e em BMI.CF3SO3 (5,0 nm). A atividade catalítica destas nanopartículas em reações de hidrogenação de arenos depende tanto da natureza do precursor quanto do líquido iônico empregado para prepará-las. Independentemente da natureza do líquido iônico, nanopartículas progressivamente mais ativas foram obtidas, primeiro, mudando o precursor na seqüência [Rh(cod)Cl]2 < [Rh(cod)2]CF3SO3 < RhCl3.3H2O < [Rh(cod)2]BF4 e, segundo, mudando o líquido iônico na sequência BMI.CF3SO3 < BMI.PF6 < BMI.BF4. É digno de nota, que as nanopartículas de [Rh(0)] obtidas pela combinação de [Rh(cod)2]BF4 / BMI.BF4 mostraram, na hidrogenação do benzeno, atividade catalítica similar à amostra de Rh/C.

Aplicação de técnicas de gerenciamento por objetivos à área industrial de uma empresa eletro-eletrônica

Este trabalho descreve a implantação da Gestão por objetivos na Área Industrial da Empresa PARKS Comunicações Digitais, que é uma empresa do segmento de telecomunicações que desenvolve, fabrica e distribui equipamentos destinados à comunicação de dados denominados de "última milha", ou seja, dentro de um sistema de tranmissão de dados e voz, seus equipamentos estão instalados próximos ao usuário final do serviço. São apresentados a geração dos indicadores de desempenho, o treinamento dos gestores da empresa e os resultados auferidos no decorrer dos primeiros anos. A implementação desta metodologia de gestão seguiu o processo de pesquisa e ação, usando técnicas de elaboração de objetivos e de seu acompanhamento. Os resultados foram avaliados através de dados financeiros da empresa, usando um indicador geral, o custo do produto vendido gerencial(CPV), e seus componentes, o custo do produto comprado, o de industrialização e o da não qualidade, cada um com seus próprios indicadores. Verificou-se que, tomando-se como base o CPV, ao longo do período, houve uma redução do mesmo que alcançou cerca de 16% do faturamento bruto obtido. Outro resultado positivo foi a formação de um grupo de gestores já treinado em levantamento, análise e determinação de metas, aguardando as diretrizes para desenvolver o Plano Estratégico da Empresa.

Solução LTSn para problemas de transferência radiativa com polarização em geometria plana

O método LTSN tem sido utilizado na resolução de uma classe abrangente de problemas de transporte de partículas neutras que são reduzidos a um sistema linear algébrico depois da aplicação da transformada de Laplace. Na maioria dos casos estudados os autovalores associados são reais e simétricos. Para o problema de criticalidade os autovalores associados são reais ou imaginários puros e simétricos, e para o o problema de multigrupo podem aparecer autovalores complexos. O objetivo deste trabalho consiste na generalização da formulação LTSN para problemas de transporte com autovalores complexos. Por esse motivo é focada a solução de um problema radiativo de transporte com polarização em uma placa plana. A solução apresentada fundamenta-se na aplicação da transformada de Laplace ao conjunto de equações SN dos problemas resultantes da decomposição da equação de transferência radiativa com polarização em série de Fourier, seguindo o procedimento de Chandrasekhar. Esse procedimento gera 2L + 2 sistemas lineares de ordem 4N dependentes do parâmetro complexo "s". Aqui, L é o grau de anisotropia e N a ordem de quadratura. A solução desse sistema simbólico é obtida através da aplicação da transformada inversa de Laplace depois da inversão da matriz simbólica pelo método da diagonalização. Para a obtenção das constantes de integração é assumido que os componentes do vetor de Stokes são reais e as matrizes dos autovalores e autovetores são separadas em suas partes real e imaginária. A solução LTSN para autovalores complexos é validada através da comparação da solução para uma placa com espessura unitária, grau de anisotropia L = 13, albedo de espalhamento simples $ = 0:99, coe ciente de re exão de Lambert ¸0 = 0:1 e N = 150, segundo dados da literatura consultada.

Avanços na determinação da distribuição em profundidade de elementos leves com resolução subnanométrica utilizando reações nucleares ressonantes

A utilização de reações nucleares ressonantes estreitas em baixas energias é uma ferramenta importante para a determinação de distribuições de elementos leves em filmes finos com resolução em profundidade subnanométrica. O objetivo do trabalho descrito ao longo da presente dissertação é aprimorar os métodos utilizados para a aquisição e interpretação das curvas experimentais. A obtenção das curvas experimentais consiste na detecção dos produtos das reações nucleares em função da energia das partículas incidentes, fazendo necessária variar a energia das partículas do feixe em passos discretos no intervalo desejado. Neste trabalho implementou-se um sistema automático para o controle e incremento da energia do feixe e monitoramento dos produtos das reações nucleares. Esse sistema de varredura automático de energia, além de aumentar consideravelmente a velocidade da medida, aumenta a qualidade das curvas experimentais obtidas. Para a interpretação das curvas de excitação experimentais, foi implementado um programa em linguagem de programação C, baseado na teoria estocástica, que permite simular numericamente as curvas de excitação. Para demonstrar a sua funcionalidade o método desenvolvido foi aplicado para verificar a estabilidade termodinâmica de filmes dielétricos ultrafinos depositados sobre silício, que foram posteriormente nitretados por plasma.

Adição de óxido de zinco (ZnO) em pastas de cimento visando viabilizar a reciclagem de pós de aciaria elétrica (PAE) na construção civil

O Pó de Aciaria Elétrica (PAE) é um resíduo sólido perigoso gerado durante o processo de fabricação do aço. Devido ao grande volume gerado, de 12 a 14 kg por tonelada de aço fabricado, faz-se necessário o estudo de alternativas de reciclagem do mesmo. A caracterização química, física e de fases é etapa fundamental para avaliar a viabilidade de reciclagem deste resíduo. Uma das alternativas de reciclagem do PAE é a adição deste resíduo na fabricação de artefatos para a construção civil. Entretanto, esta adição retarda o início das reações de hidratação dos artefatos, dificultando seu uso. Esse retardo é atribuído aos compostos de zinco (Zn) presentes no resíduo. Este trabalho foi desenvolvido em duas etapas. Na primeira etapa foi realizada a caracterização química, física e de fases. Os resultados analíticos indicam que o PAE possui um tamanho médio de partícula que possibilita sua adição na produção de artefatos de cimento. Em relação às fases cristalinas as amostras de PAE estudadas contêm: ZnFe2O4, Fe3O4, FeCr2O4, Ca0,15Fe2,85O4, SiO2 e ZnO. Na segunda etapa deste trabalho, o ZnO, uma das formas cristalinas do Zn detectadas no PAE, foi adicionado em pastas de cimento. O objetivo desta etapa foi estudar os fenômenos relacionados ao início das reações de hidratação destas pastas. Essas amostras foram analisadas via Difração de Raios-X (DRX), Microscopia Eletrônica de Varredura (MEV) e mapeamento por raios-x característicos em diferentes idades (1, 2, 4, 7 e 14 dias), a fim de rastrear os fenômenos ocorridos durante a hidratação das pastas de cimento Os resultados obtidos nesta segunda etapa indicam que no primeiro momento antes da hidratação (1 dia) o Zn apresentava-se na forma de ZnO. A partir do quarto dia observa-se a formação de uma nova fase contendo Zn: CaZn2(OH)6 . 2H2O e, sugerindo a transição de Zn na forma de óxido para hidróxido. Isto se torna mais evidente no sétimo dia, onde é possível constatar que as reações de hidratação aconteceram (pasta totalmente endurecida), pois neste caso, o Zn é encontrado somente na forma de hidróxido. A formação do composto CaZn2(OH)6 . 2H2O permite que as reações de hidratação da pasta de cimento contendo Zn sejam efetivadas. O trabalho finaliza com a apresentação de um modelo onde se busca explicar o fenômeno observado durante a execução do mesmo.

Síntese de pós de alumina nanocristalina por combustão em solução

Esta dissertação teve como objetivo o estudo e desenvolvimento da técnica de síntese por combustão em solução para obtenção de nanocristais de alumina, e a caracterização desse material por diferentes técnicas. O estudo da obtenção incluiu o modelamento termodinâmico das reações investigadas e a influência da razão combustíveloxidante no tamanho de partícula. A síntese da alumina por combustão foi realizada a partir do par precursor químico glicina-nitrato de alumínio em água, para diferentes razões combustível-oxidante, temperaturas de chama e número total de moles gasosos, relacionando esses parâmetros com tamanho de cristalito no pó, deformação de rede, área superficial e porosidade total. As reações químicas de combustão em solução resultaram na formação in situ de fases metaestáveis de Al2O3 para o pó como-sintetizado. A fase Al2O3-α, cuja formação não ocorreu nem na presença de um extra-oxidante, somente foi obtida após um tratamento térmico a 1100°C do pó como-sintetizado. O modelamento termodinâmico da reação de combustão em solução mostrou que quando a razão combustível-oxidante aumenta, obtém-se uma elevação da temperatura de chama adiabática e da quantidade de gás produzida, definindo características do particulado como morfologia, tamanho de cristalito, área superficial e nível de aglomeração. A caracterização do pó utilizou técnicas como análise termodiferencial (ATD) e termogravimétrica (ATG), granulometria por dispersão de laser (GDL), BET para análise de área superficial e porosidade total, difração de raios X (DRX), microscopia eletrônica de varredura (MEV), e microscopia eletrônica de transmissão (MET). Os pós de alumina alfa obtidos apresentaram-se como constituídos de cristalitos nanométricos, dispostos na forma de agregados de tamanho micrométrico. O tamanho de cristalito médio, calculado pelo método de Williamson-Hall, a partir de dados de análises por DRX, foi de 97 nm, e deformação média calculada de 5x10-3 % A análise granulométrica dos pós como-sintetizados verificou um tamanho de aglomerado de 40 µm, com distribuição monomodal. As análises por MEV e MET confirmaram a aglomeração das partículas do pó como-sintetizado, a partir de partículas com morfologia irregular.

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.

Reações de acoplamento para a obtenção de olefinas

Este trabalho teve como objetivo desenvolver novas metodologias para as reações de acoplamento C-C e C-O para a obtenção de olefinas. Na primeira etapa do trabalho foi desenvolvido um procedimento simples e eficiente para a preparação de estirenos funcionalizados, a partir de reações de acoplamento catalisadas por paládio entre ácidos arilborônicos e o brometo de vinila, gerado in situ a partir do 1,2-dibromoetano. A reação utiliza condições brandas e uma variedade de grupos funcionais é tolerada. Os ácidos arilborônicos contendo grupos eletrodoadores mostraram-se mais reativos na reação de acoplamento. O sistema é muito simples quando comparado com outros sistemas propostos para a obtenção de vinilarenos, não tem o inconveniente de usar etileno sob pressão ou reagentes tóxicos de estanho. Os ésteres α-arilacrílicos são precursores para os ácidos α-arilpropiônicos, importantes não esteroidais, e podem ser obtidos pela reação de acoplamento catalisada por paládio do 2-bromoacrilato de butila com ácidos arilborônicos. Nesse trabalho foram realizados vários experimentos para determinar as melhores condições reacionais. Os melhores resultados foram obtidos com a utilização de Pd(OAc)2 como catalisador, a PCy3 como fosfina ligante, o Cs2CO3 como base e DMA como solvente. Utilizando essas condições, a uma temperatura de 120 oC durante 1.5 h, o produto de acoplamento, o 2-fenilacrilato de butila, foi obtido com 83% de rendimento Na última parte foi investigada a formação de ligações C-O catalisadas por metais. A reação de acoplamento entre o 2-bromoacrilato de butila com fenóis possibilita a obtenção de precursores insaturados para os ácidos α-ariloxipropanóicos, uma importante classe de herbicidas.O sistema catalítico baseado em cloreto de cobre mostrou-se mais ativo do que os sistemas com paládio. O melhor sistema catalítico foi obtido com cloreto de cobre como catalisador, na presença de Cs2CO3 como base e tolueno como solvente. Utilizando essas condições para o acoplamento entre o 2-bromoacrilato de butila e t-butilfenol foi possível obter o 2-(p-tbutilfenoxi)acrilato de butila com 37% de rendimento.

Síntese, caracterização e aplicação catalítica de nanopartículas de rutênio suportadas em líquidos iônicos

A decomposição controlada do precursor organometálico [Ru(cod)(cot)] (cod = ŋ4-1,5 ciclooctadieno e cot = ŋ6-1,3,5 ciclooctatrieno) com hidrogênio molecular disperso em líquidos iônicos, derivados da associação de cátions 1-n-butil-3-metilimidazólio com ânions PF6-, BF4- e CF3SO3-, leva a formação de nanopartículas de rutênio. O diâmetro médio destas partículas determinado por Microscopia Eletrônica de Transmissão (MET) e Difração de Raios X (DRX) foi de 1,9 – 2,4 e 2,7 nm, respectivamente. A análise por Espectroscopia de Fotoelétrons Induzidos por Raios X (XPS) mostra claramente uma contribuição importante da componente Ru-O indicando a formação de uma camada de passivação nas nanopartículas de rutênio, que desaparece após bombardeamento com Ar+. Estudos de Espalhamento de Raios X em Baixo Ângulo (SAXS) mostram a formação de uma camada semicristalina de líquido iônico em torno das nanopartículas. Na hidrogenação do 1-hexeno com nanopartículas de [Ru(0)]n.PF6 foi possível realizar de 8 e 9 reciclos nos líquidos iônicos BMI.PF6 e BMI.BF4, respectivamente. Na hidrogenação do benzeno pela ação de nanopartículas de rutênio verificou-se seletividade para cicloexeno de até 39% a baixas conversões. Uma série de constantes de velocidades relativas das reações competitivas de hidrogenação de compostos aromáticos monosubstituídos por grupos alquila foram correlacionadas com os parâmetros estéreos de Taft, obtendo-se valores satisfatórios para os coeficientes de correlação. Os resultados mostraram que a velocidade relativa de hidrogenação diminui com o aumento do substituinte alquila.

Obtenção de alfa-olefinas em reações de oligomerização do eteno utilizando complexos de níquel em meio homogêneo e bifásico

Neste trabalho são descritos estudos comparativos de reações de oligomerização do etileno utilizando complexos do tipo: Ni-[P^O] e Ni-[Diimina] (na presença de MAO como cocatalisador). A pesquisa foi realizada com o objetivo de verificar e comparar o comportamento dos dois complexos em estudo, quando estes passam do meio homogêneo para o meio bifásico (na presença de um líquido iônico). O líquido iônico utilizado é o tetrafluorborato de 1-butil-3-metilimidazol. Foram avaliados parâmetros reacionais como pressão, temperatura, tempo de reação, velocidade de agitação e volume de líquido iônico para o catalisador Ni-[P^O] nos meios homogêneo e bifásico. Para o precursor Ni-[Diimina], os parâmetros avaliados foram pressão e temperatura nos dois meios. Os produtos obtidos são olefinas lineares e ramificadas, principalmente α-olefinas; com faixa de distribuição de C4 a C20 (Distribuição Shulz-Flowry). As atividades obtidas nos ciclos de reações variam entre 21 h-1 a 32.772 h-1, e os resultados são interpretados envolvendo os mecanismos de reação. Devido as diferenças existentes entre nos complexos em estudo, verifica-se diferentes atividades e diferentes faixas de produtos oligoméricos. O comportamento obtido para cada complexo depende diretamente dos parâmetros reacionais utilizados, além do meio ao qual foram submetidos durante as reações.

Contraste iodado em tomografia computadorizada : reações adversas imediatas em pacientes hospitalizados

Pesquisa quantitativa do tipo exploratório-descritivo, prospectivo, de caráter não experimental. O objetivo consistiu em conhecer as reações adversas imediatas ao contraste iodado intravenoso em pacientes internados, submetidos a tomografia computadorizada num hospital escola, sendo campo do estudo a Unidade de Tomografia do Hospital de Clínicas de Porto Alegre. A coleta de dados ocorreu através de planilha de registros preenchida pela equipe de enfermagem atuante na Unidade, durante as 24 horas de funcionamento do setor, no período entre outubro e dezembro de 2004, totalizando 351 pacientes observados. Os registros incluíram dados relativos à dinâmica do exame, características da clientela e ocorrência dos eventos em estudo. Para o tratamento dos dados, recorreu-se à estatística descritiva e descritiva, com emprego dos softwares SPSS, EPI INFO e PEPI. Os eventos adversos foram considerados imediatos quando ocorreram até 30 minutos após a injeção do contraste. Durante o período do estudo, todas as reações adversas apresentaram intensidade leve, manifestando-se com freqüência de 12,5% entre os 160 pacientes que receberam contraste iodado iônico e 1,0% entre os 191 pacientes que receberam contraste não iônico (p=0,000). O emprego do meio não iônico mostrou-se eficaz na prevenção de reações adversas ao contraste iodado, mesmo na presença de condições clínicas que aumentam o risco para ocorrência desses eventos. A administração intravenosa do contraste através de bomba injetora aumentou significativamente o percentual de reações adversas, em comparação com a injeção manual (p=0,013). O extravasamento de contraste, considerado um evento adverso local, ocorreu em 2,2% das 317 injeções em veia periférica, e os volumes extravasados oscilaram entre 1 e 15 mililitros, com média equivalente a 4,4 mililitros, não ocasionando complicações em nenhum dos casos. A administração intravenosa de contraste através de cateter plástico ocasionou freqüência significativamente menor de extravasamentos do que o emprego de agulha metálica (p=0,041). Os índices evidenciados no presente estudo encontram-se dentro dos limites que constam na revisão de literatura, apesar das pesquisas internacionais apresentarem diferenças entre si na seleção das amostras ou critérios para definição das reações adversas imediatas. Entre as recomendações, sugere-se que os serviços de tomografia conheçam as próprias taxas de reações adversas ao contraste iodado e as condições em que elas ocorrem, a fim de obter evidências para avaliação dos respectivos processos assistenciais.

Reações adversas a medicamentos: uma coorte em hospital universitário

Resumo não disponível.

Projeto e implementação de um sistema de instrumentação eletro-eletrônica para caracterização de escoamentos através de processamento digital de imagens

Historicamente, devido à incompreensibilidade dos fenômenos físicos envolvidos, diversas técnicas de análise experimental de escoamentos foram desenvolvidas. Inicialmente, estas técnicas consistiam em uma simples visualização do escoamento em questão. Visava-se com isto produzir a intuição física a respeito do sistema físico para que futuras teorias pudessem ser formuladas. No entanto, nos anos posteriores, tornou-se necessária a criação de técnicas quantitativas. Neste segmento, a evolução tecnológica da eletrônica propiciou que técnicas intrusivas e não intrusivas fossem criadas. Tal processo culminou com o surgimento da microeletrônica, quando sensores eletrônicos de imagem puderam ser fabricados. Assim, técnicas não intrusivas de velocimetria através de processamento digital de imagens tornaram-se referência, possibilitando a medição de campos de velocidade do escoamento em análise. Tais técnicas são atualmente utilizadas não apenas no meio acadêmico, na investigação da turbulência, como na indústria, em projetos de engenharia que envolvam alta tecnologia. Neste trabalho, apresenta-se uma metodologia para a concepção de um sistema de velocimetria de 2 dimensões por processamento digital de imagens. Para tanto, o sistema de velocimetria foi dividido em duas partes, sendo os equipamentos bem como os conceitos físicos envolvidos em cada uma analisados e descritos Por fim, projetou-se e implementou-se um sistema de velocimetria de 2 dimensões através de processamento digital de imagens, modo PIV. A sua validação foi efetuada sobre um escoamento altamente turbulento interno a um duto circular, cuja solução é conhecida e exaustivamente estudada no decorrer dos anos. O resultado obtido se apresentou totalmente satisfatório e demonstrou a grande importância atual de tal técnica de medição.

Um protocolo de perfil lite para colaboração visual que utiliza transferência de arquivos com t.127

Atualmente, quando falamos em Internet, não podemos somente pensar em serviços que utilizam correio eletrônico ou acessam informações de sites através de um navegador. Dentre alguns serviços utilizados na Internet, podemos destacar a videoconferência. Videoconferência contempla, além do intercâmbio de áudio e vídeo entre duas pessoas ou grupos, também o compartilhamento de dados. Os serviços e protocolos definidos nas recomendações H.323 (para videoconferência) e T.120 (para a colaboração de dados) do ITU são bastante complexos e ocorrem muitos problemas na sua utilização em redes de pacotes funcionando segundo o princípio de best effort da Internet. Nesta dissertação de mestrado são apresentados alguns resultados do estudo realizado sobre o contexto de videoconferência, suas soluções e protocolos, com ênfase nos protocolos padronizados para colaboração de dados dentro de uma videoconferência e a sua estrutura. Esta dissertação também apresenta uma proposta de solução para uma aplicação que utilize transferência de arquivos nos moldes do padrão ITU T.127, mas que atenda aos requisitos de menor complexidade (tráfego e processamento). A proposta utiliza as estratégias de simplificação dos protocolos usados no contexto de colaboração de dados em ambiente que ainda mantenha compatibilidade com o ambiente T.120 e, sobretudo utilizando a Internet.

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.