Simulação de escoamentos incompressíveis não newtonianos em dutos com expansão brusca

O objetivo deste trabalho é a simulação numérica de escoamentos incompressíveis bidimensionais em dutos com expansão brusca, considerando o raio de expansão de 3 : 1. As equações governantes implementadas são as de Navier, que junto com relações constitutivas para a tensão visam representar comportamentos não newtonianos. A integração temporal é feita usando o esquema explícito de Runge-Kutta com três estágios e de segunda ordem; as derivadas espaciais são aproximadas pelo método de diferenças finitas centrais. Escoamentos em expansões bruscas para fluidos newtonianos apresentam um número de Reynolds crítico, dependente do raio de expansão, na qual três soluções passam a ser encontradas: uma solução sim étrica instável e duas soluções assimétricas rebatidas estáveis. Aumentando o número de Reynolds, a solução passa a ser tridimensional e dependente do tempo. Dessa forma, o objetivo é encontrar as diferenças que ocorrem no comportamento do fluxo quando o fluido utilizado possui características não newtonianas. As relações constitutivas empregadas pertencem à classe de fluidos newtonianos generalizados: power-law, Bingham e Herschel-Bulkley. Esses modelos prevêem comportamentos pseudoplásticos e dilatantes, plásticos e viscoplásticos, respectivamente. Os resultados numéricos mostram diferenças entre as soluções newtonianas e não newtonianas para Reynolds variando de 30 a 300. Os valores de Reynolds críticos para o modelo power-law não apresentaram grandes diferenças em comparação com os da solução newtoniana. Algumas variações foram percebidas nos perfis de velocidade. Entretanto, os resultados obtidos com os modelos de Bingham e Herschel-Bulkley apresentaram diferenças significativas quando comparados com os newtonianos com o aumento do parâmetro adimensional Bingham; à medida que Bingham é aumentado, o tamanho dos vórtices diminui. Além disso, os perfis de velocidade apresentam diferenças relevantes, uma vez que o fluxo possui regiões onde o fluido se comporta como sólido.

Transferência de calor combinando radiação e convecção no interior de dutos de geradores de vapor fumotubulares

Esta dissertação de mestrado considera a transferência de calor combinando convecção e radiação térmica no escoamento de gases participantes em dutos de seção circular. Partindo de uma metodologia geral, o trabalho enfoca principalmente os casos típicos de aplicação em geradores de vapor fumotubulares de pequeno e médio porte, em que gases em alta temperatura escoam através de um tubo mantido em temperatura uniforme. O escoamento é turbulento e o perfil de velocidade é plenamente desenvolvido desde a entrada do duto. A temperatura do gás, contudo, é uniforme na entrada, considerando-se a região de desenvolvimento térmico. Duas misturas de gases são tratadas, ambas constituídas por dióxido de carbono, vapor d’água e nitrogênio, correspondendo a produtos típicos da combustão estequiométrica de óleo combustível e metano. As propriedades físicas dos gases são admitidas uniformes em todo o duto e calculadas na temperatura de mistura média, enquanto que as propriedades radiantes são modeladas pela soma-ponderada-de-gases-cinzas. O campo de temperatura do gás é obtido a partir da solução da equação bidimensional da conservação da energia, sendo os termos advectivos discretizados através do método de volumes de controle com a função de interpolação Flux-Spline; as trocas de energia radiantes são avaliadas por meio do método das zonas, onde cada zona de radiação corresponde a um volume de controle. Em um primeiro passo, a metodologia é verificada pela comparação com resultados apresentados na literatura para a transferência de calor envolvendo apenas convecção e combinando convecção com radiação. Em seguida, discutem-se alguns efeitos da inclusão da radiação térmica, por exemplo, no número de Nusselt convectivo e na temperatura de mistura do gás. Finalmente, são propostas correlações para o número de Nusselt total, que leva em conta tanto a radiação quanto a convecção. Essa etapa exige inicialmente uma análise dos grupos adimensionais que governam o processo radiante para redução do número elevado de parâmetros independentes. As correlações, aplicáveis a situações encontradas em geradores de vapor fumotubulares de pequeno e médio porte, são validadas estatisticamente pela comparação com os resultados obtidos pela solução numérica.

A climatologia aplicada ao ambiente construído : uma contribuição à caracterização climática de Porto Alegre/RS

Com caráter informativo, inicia-se apresentando a relação entre a climatologia e o ambiente construído, em seus campos especificos de aplicação (urbanismo, projeto de edificações, execução de obras, e manutenção e uso das construções), considerando-se os aspectos históricos, econômico-ecológicos e os ligados à súde e ao conforto do usuário, abrangidos por tal relação. A importância da disponibilidade de informações climatológicas é então comentada, bem como diversas metodalagias de análise e representação de dados climáticos. Com caráter aplicativo, apresenta-se uma contribuigão à caracterização climática da cidade de Porto Alegre, RS. São analisadas os valores médios horários-mensais doa principais parâmetros climáticos disponíveis (temperatura, umidade, velocidade e direção do vento, e nebulosidade), propondo-se uma representação gráfica sintética alternativa para os mesmos. A análise é complementada pala execução de dois procedimentos técnicos: avaliação do conforto térmico, com identificação do período critico de verão, e projeto de dispositivo de sombreamento de aberturas retangulares. Embora esteja em foco o caso específico de Porto Alegre, tais metodologias encontram-se descritas de modo a serem aplicadas a qualquer outra localidade, procurando-se ampliar a utilidade prática do presente estudo.

Contribuição ao estudo da viabilidade da cura do concreto por energia solar

O presente trabalho visou a utilização da energia solar na cura de pré-fabricados de concreto. Para tanto, foi utilizado dois modelos. O primeiro modelo constou da utilização de coletores solares planos, ligados até uma câmara de cura por um sistema de dutos; A radiação direta do sol incide sobre os coletores, que acumulam energia, que é removida pela água que circula pelos mesmos, levando essa energia através do sistema de distribuição (dutos) até a câmara de cura; Essa câmara, possui uma pista metálica em sua base, onde a energia é cedida para os pré-fabricados de concreto; Esse ciclo de captação distribuição e consumo é repetido a cada hora do dia, onde há radiação solar; Um segundo modelo constou de uma heliocâmara, princípio de um destilador solar convencional, onde a radiação incide diretamente sobre a câmara de cura; A radiação, penetra pela cobertura e forma-se um efeito estufa no interior da câmara, aquecendo seu interior e consequentemente os pré-fabricados de concreto; Avaliou-se os dois sistemas para as estações de verão e inverno na cidade de Porto Alegre (dias típicos) e concluiu-se que para o verão, tanto o sistema com coletores solares quanto o sistema com heliocâmaras poderão ser utilizados. No período de inverno, nenhum destes sistemas serão satisfatórios, porém suas utilizações serão viáveis em conjunto com uma fonte de energia alternativa; A energia solar constitui-se assim numa importante fonte de energia a ser utilizada em pré-fabricados de concreto.

Simulação numérica de esteiras em transição utilizando o método dos contornos virtuais

Simulações Numéricas são executadas em um código numérico de alta precisão resolvendo as equações de Navier-Stokes e da continuidade para regimes de escoamento incompressíveis num contexto da turbulência bidimensional. Este código utiliza um esquema compacto de diferenças finitas de sexta ordem na aproximação das derivadas espaciais. As derivadas temporais são calculadas usando o esquema de Runge-Kuta de terceeira ordem com baixo armazenamento. Tal código numérico fornece uma representação melhorada para uma grande faixa de escalas de comprimento e de tempo. As técnicas dos contornos imersos acopladas ao método dos contornos virtuais permitem modelar escoamentos não-estacionários sobre geometrrias complexas, usando simplesmente uma malha Cartesiana uniforme. Por meio de procedimentos de aproximação/interpolação, as técnicas dos contornos imersos (aproximação Gaussiana, interpolação bilinear e redistribuição Gaussiana), permitem a representação do corpo sólido no interior do campo de escoamento, com a superfície não coincidindo com a malha computacional. O método dos contornos virtuais, proposto originalmente por Peskin, consiste, basicamente, na imposição na superfície e/ou no interior do corpo, de um termo de força temporal acrescentando às equações do momento. A aplicação deste campo de força local leva o fluido ao repouso na superfície do corpo, permitindo obter as condições de contorno de não-deslizamento e de não penetração de fluido na parede. A análise das oscilações induzidas no escoamento-contorno pelo processo de desprendimento de vórtices na esteira do cilindro circular e de geometria retangulares na incidência, para números de Reybolds variando de 40 a 400, confirma a eficiência computacional e a aplicabilidade das técncias implementadas.

Análise comparativa de desempenho de redes IP e ATM com tráfego multimídia interativo

A década de 80 é um marco para a área de comunicação de dados. Muitos estudos, pesquisas e especificações foram iniciados para permitir a integração de sistemas de comunicação para voz, vídeo e dados. Desde essa década, os cientistas procuram investir na Internet, que é um dos mais importantes e populares meios de transmissão de informações no mundo, para acomodar aplicações que permitam enviar imagens e sons por esta imensa rede de comunicação. Também nessa década, no final, o International Telecommunications Union – Telecommunication (ITU-T), especificou a tecnologia ATM, cujas pesquisas vinham ocorrendo desde o seu início. O serviço de rede tradicional na Internet é transmissão de datagramas ´besteffort ´, conforme será visto neste trabalho. Neste serviço, os pacotes da origem são enviados para o destino, sem garantia de entrega. Para aquelas aplicações que requerem garantia de entrega, o protocolo TCP (equivalente à camada 4 do RM-OSI) se encarregará da retransmissão dos pacotes que falharam ao atingir o destino para então se conseguir uma recepção correta. Para aplicações de comunicação tradicionais, como FTP e Telnet, em que uma entrega correta é mais importante que a perda de tempo da retransmissão, este serviço é satisfatório. Entretanto, uma nova classe de aplicações, as quais usam mídias múltiplas (voz, som e dados), começam a aparecer na Internet. Exemplos desta classe de aplicação são: vídeo teleconferência, vídeo sob demanda e simulação distribuída. Operações de modo tradicional para estas aplicações resultam em redução da qualidade da informação recebida e, potencialmente, ineficiência do uso da largura de banda. Para remediar este problema é desenvolvido um ambiente de serviço para tempo real, no qual múltiplas classes de serviços são oferecidas. Este ambiente estende o modelo de serviços existentes para ir ao encontro das necessidades das aplicações multimídia com obrigatoriedade de tempo real, porém esta não é uma meta muito fácil. Contudo, a comunidade pesquisadora tem conseguido desenvolver alguns mecanismos que vão pouco a pouco permitindo que este objetivo se transforme em realidade. O ATM é uma tecnologia que provê dutos de alta velocidade para a transferência de qualquer tipo de informação em pequenos pacotes de tamanho fixo, chamados células. A idéia básica é estabelecer entre dois pontos que desejam se comunicar um circuito virtual que é mantido pelos comutadores de células para levar a informação de um lado a outro. A característica marcante do ATM é a Qualidade de Servico – QoS, que garante o desempenho predefinido que determinado serviço necessita. Isso permite suportar aplicações de tempo real que são sensíveis ao atraso. O atendimento à diversidade de características de tráfego e exigências distintas de grandes quantidades de serviços é possível pelo ATM graças ao controle de tráfego reunido à capacidade de multiplexação estatística do meio em altas taxas de transmissão. O objetivo principal desta dissertação é elaborar uma comparação quantitativa e qualitativa do funcionamento de aplicações multimídia sobre IP com RSVP - Protocolo desenvolvido para realizar reservas de recursos integrante da arquitetura de Serviços Integrados (IntServ) proposta pelo IETF para oferecer qualidade de serviço para aplicações tais como aplicações multimídia. Essa proposta também inclui duas classes de serviços, sendo que nessa dissertação, o serviço de carga controlada é que está sendo utilizado. Isso deve-se a implementação dos módulos apresentados em [GRE 2001] e que foram utilizados na realização desse trabalho - e sobre ATM. A proposta final é a elaboração de uma técnica de análise baseado nas principais métricas de desempenho em redes que deve permitir uma melhor visualização do comportamento das tecnologias sob determinadas cargas de tráfego, permitindo assim uma tomada de decisão entre qual das tecnologias que deverá ser adotada em um dado momento, em uma dada situação, ou seja, a indicação do ponto de quebra tecnológica na situação modelada. Para que fosse possível fazer esta comparação, foi necessário dividir-se este trabalho em 3 grandes etapas, que são: • Estudo e desenvolvimento da técnica para análise do elemento carga de tráfego na tecnologia ATM; • Estudo e desenvolvimento da técnica para análise do elemento carga de tráfego na tecnologia IP com RSVP; • Comparativo quantitativo e qualitativo dos estudos acima.

Verificação do estado limite de serviço de abertura das fissuras em seções de concreto armado submetidas à flexão simples

Este trabalho apresenta um programa que tem por finalidade automatizar os procedimentos de verificação do estado limite de serviço de abertura das fissuras em seções retangulares e T de concreto armado, submetidas à flexão simples, conforme as recomendações da norma NBR 6118/2003 – Projeto de Estruturas de Concreto. O programa foi implementado na linguagem de programação Microsoft Visual Basic, aproveitando os recursos fornecidos por esta linguagem, que propiciam o desenvolvimento de uma entrada de dados amigável, típica do ambiente Windows. A fissuração, em elementos de concreto armado, ocorre quando as tensões de tração no concreto excedem a sua capacidade resistente. Nos elementos fissurados, as aberturas das fissuras devem ser controladas. Segundo as recomendações da norma, o valor característico da abertura das fissuras não deve ultrapassar determinados valores admissíveis. No procedimento implementado, são determinadas as aberturas das fissuras considerando-se o conjunto de barras da armadura tracionada e a camada mais tracionada da armadura, com o propósito de estabelecer a situação mais desfavorável. O programa tem também a pretensão de auxiliar estudantes de engenharia na compreensão do fenômeno da fissuração, fornecendo valores de elementos que caracterizam o funcionamento da peça fletida.

Estudo de cargas impulsivas com ênfase em explosões : estratégias de análises utilizando métodos numéricos

A avaliação da solicitação produzida por explosões, assim como da resposta de estruturas, são temas de muito interesse na engenharia atualmente, tanto pela quantidade crescente de acidentes relacionados com explosões quanto pelas ações terroristas muitas vezes vinculadas a estes tipos de ações. Neste contexto, o presente trabalho tem por objetivo explorar técnicas de análise tanto na modelagem da excitação quanto na resposta de estruturas consideradas como alvos. Para isto, são utilizadas metodologias de diferentes tipos: implementações baseadas em sistema de elementos finitos comerciais como Ansys [2000] e LS-Dyna [2003] e técnicas simplificativas que permitem realizar uma avaliação preliminar. As aplicações consideradas são indicadas a seguir: • Análise da Resposta de Estruturas Laminares Submetidas à Ação de Cargas Explosivas: determina-se a pressão produzida por explosivos sólidos a certa distância do epicentro, através de métodos simplificados, determinando a resposta esperada em placas retangulares; • Efeito da Pressão Interna em Vasos de Pressão (Extintores de Incêndio): comparando resultados numéricos e experimentais verifica-se a influência da pressão interna nas propriedades dinâmicas do sistema; • Estudo de Um Vaso Esférico de GLP Sob Ação de Uma Carga Explosiva: aplica-se a ação de uma onda explosiva produzida por um gás inflamável pesado sobre uma estrutura de vaso de pressão esférico com fluido e gás em seu interior, determinando sua resposta, avaliando também a influência de diferentes quantidades de líquido e pressão interna na resposta da estrutura; • Modelamento de uma Camada de Solo / Propagação das Ondas: verifica-se o comportamento da propagação de ondas em um meio elástico, comparando valores encontrados em testes experimentais. Analisa-se numericamente o efeito da inserção de uma valeta na atenuação de tais ondas; • Simulação Numérica Completa de uma Explosão: modela-se um semi-espaço submetido à ação de um explosivo sólido sobre sua superfície, avaliando os campos de pressão gerados. Ao final de cada aplicação são apresentadas conclusões parciais obtidas e as possibilidades de trabalhos futuros vislumbrados. Finalmente, conclui-se que as técnicas empregadas para as simulações são extremamente eficientes, considerando conhecidos todos os parâmetros envolvidos em cada modelo. Neste ponto é fundamental o trabalho do engenheiro, utilizando-se de seus conhecimentos técnicos para transformação do evento real em um modelo numérico, considerando e selecionando as simplificações necessárias.

Biologia reprodutiva do lagostim Parastacus varicosus Faxon, 1898 (Decapoda : Parastacidae) da Bacia do Rio Gravataí, Rio Grande do Sul

Considerando à lacuna sobre as informações da biologia reprodutiva dos crustáceos de água doce no Brasil, principalmente dos lagostins, a presente pesquisa teve como objetivo caracterizar aspectos da dinâmica da reprodução, como a descrição e desenvolvimento das gônadas, e analisar as variações bioquímicas do metabolismo intermediário de Parastacus varicosus Faxon, 1898 para relacioná-las ao período reprodutivo. Os animais foram amostrados de junho de 2004 a maio de 2005 em um arroio pertencente à Bacia do Rio Gravataí, localidade Cova do Touro, Gravataí, Rio Grande do Sul. As fêmeas foram separadas em três estágios gonadais de acordo com a maturação. Todos animais apresentaram os dois pares de aberturas e dutos genitais. Foram identificados histologicamente 40 machos, 32 fêmeas e 6 espécimes transicionais. Os índices gonadossomático (IG) e hepatossomático (IH) foram determinados. Ocorreu um aumento do IG e uma diminuição no IH, no período de passagem da primavera para o verão, indicando uma maior atividade reprodutiva no verão e uma possível transferência das reservas do hepatopâncreas para a gônada nesse período. Nas análises bioquímicas da hemolinfa, hepatopâncreas, músculo e gônadas os níveis de glicose, glicogênio, proteínas totais, colesterol e lipídeos totais foram quantificados. As variações sazonais observadas mostraram estar relacionada com o período reprodutivo, mostrando uma transferência de parte das reservas energéticas do hepatopâncreas e do músculo para o desenvolvimento das gônadas durante período reprodutivo. Nessa pesquisa portanto foi observada uma correlação entre os resultados obtidos das investigações das gônadas e das análises bioquímicas indicando assim uma conexão entre essas pesquisas e contribuindo para um melhor entendimento da biologia reprodutiva de P. varicosus.

Sedimentação carbonática de algas vermelhas coralináceas da plataforma continental da bacia de campos um modelo carbonático análogo para o terciário

O presente trabalho faz uma análise da plataforma continental média e externa da porção central da Bacia de Campos, entre as batimetrias de 65 e 300 metros. Nesta porção da plataforma são encontrados depósitos sedimentares formados por carbonatos e, secundariamente, por siliciclásticos com influência carbonática. Os carbonatos do Membro Grussaí da Formação Emborê se configuram como os mais importantes sedimentos presentes na área de estudo, sendo compostos majoritariamente por algas vermelhas coralináceas (lithothamniun). Os depósitos algálicos ocorrem sob duas formas: a) formas livres formando concreções nodulares centimétricas (rodolitos) e grainstones. b) Incrustações, formando substratos crustosos e lineamentos. Para a análise destes depósitos foram utilizados dados sísmicos tridimensionais relativos à seção terciária da Bacia de Campos e principalmente do fundo marinho. A associação destes dados com amostras de fundo e filmagens do fundo marinho ao longo do traçado de dutos de escoamento de hidrocarbonetos, permitiu a definição de um modelo de distribuição faciológica para os carbonatos formados por algas vermelhas coralináceas e feições fisiográficas associadas a eles. A morfologia do fundo do mar foi dividida em províncias fisiográficas diretamente relacionadas com o padrão deposicional das algas vermelhas e, secundariamente, por sedimentos siliciclásticos. Os padrões morfológicos refletem a presença de bancos carbonáticos de morfologia rugosa, lineamentos carbonáticos em forma de cristas com direção paralela à linha de costa atual, lineamentos em forma de barras obliquas à linha de costa atual, e platôs ou patamares, caracterizados por relevo suave Os sedimentos de fundo foram agrupados em seis fácies sedimentares, as quais estão relacionadas com a morfologia do fundo marinho: (a) Fácies Areias Siliciclásticas com Bioclastos - depositada nos patamares e formada por sedimentos mistos. (b) Fácies Grainstones com Rodolitos – desenvolve-se nos patamares e é composta, basicamente, por grainstones dominados por algas vermelhas, rodolitos e, secundariamente, por areias siliciclásticas. (c) Fácies Boundstone Ridges – relaciona-se com a formação dos lineamentos carbonáticos paralelos à linha de costa; apresentam grainstones e rodolitos na base dos lineamentos e um predomínio de boundstones nas partes altas; marcam a passagem entre patamares. (d) Fácies Boundstone Mounds - relaciona-se fisiograficamente aos lineamentos de barras carbonáticas com orientação obliqua à linha de costa, e se localizam no interior dos patamares. (e) Fácies Boundstone Banks e Grainstones - corresponde a bancos carbonáticos de relevo irregular e que resultam da bioconstrução algálica sob forma de boundstones, com deposição de grainstones nos interbancos. (f) Fácies Mista - posicionada no talude superior, é resultado da mescla entre os sedimentos siliciclásticos de granulometria silte até areia muito fina com lamas e areias carbonáticas. A distribuição faciológica das algas vermelhas coralináceas presentes na área de estudo se desenvolve em batimetrias muito superiores ao que normalmente é esperado para os carbonatos. Sua atividade se desenvolve desde a batimetria de 45 m até 150 m, sendo que o intervalo mais produtivo ocorre entre 55 e 90 m Esses valores decorrem da análise de incrustações presentes em dutos situados na plataforma e talude continental. Outra indicação da existência de deposição em batimetrias mais elevadas é a presença de relictos de boundstones entre as atuais batimetrias de 170 a 230 m. Estes relictos indicam que, durante o rebaixamento do nível do mar de 120 m que ocorreu no final do Pleistoceno, a fábrica carbonática responsável por estes depósitos permaneceu ativa em batimetrias próximas a 50 e 110 m. Por estarem se desenvolvendo em batimetrias elevadas, os depósitos sedimentares estão submetidos aos efeitos hidrodinâmicos de correntes marinhas, principalmente a Corrente do Brasil. Devido a isto, a distribuição faciológica dos carbonatos algálicos apresenta diferença quando comparada com o ambiente convencional de deposição carbonática. Os substratos dos principais lineamentos carbonáticos, posicionados nas atuais cotas batimétricas de 110 m, 100 m, 85 m e 70 m, são interpretados como decorrentes das estabilizações episódicas de linha de praia que ocorreram durante o último processo transgressivo. A erosão desenvolvida pelas ondas nestas paleolinhas de praia esculpiu pequenas escarpas que serviram de substrato para a instalação posterior da Fácies Boundstone Ridges a partir de batimetrias entre 45 a 50 m. O substrato gerado nos patamares não desenvolveu rigidez suficiente para a instalação de formas incrustantes, exceto onde a remobilização do fundo marinho, durante o último evento transgressivo, desenvolveu formas de leitos tipo dunas transversas que serviram de substrato para a deposição da Fácies Boundstone Mounds, as quais atualmente estão estabilizadas pelas incrustrações de boundstones O modelo de deposição de algas coralináceas de “águas profundas” proposto neste trabalho pode ser o mais apropriado como análogo para a sedimentação dominada por algas vermelhas do Terciário das bacias meridionais Brasileiras, quando não houver uma clara associação com outros constituintes carbonáticos exclusivos de águas rasas. Uma evidência que sinaliza a viabilidade deste modelo decorre da interpretação e análise de atributos sísmicos da seção Terciária da Bacia de Campos, onde as feições carbonáticas, semelhantes às encontradas no recente, tendem a se repetir com freqüência.

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.