Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Métodos computacionais em engenharia mecânica

Este estudo pretende mostrar algumas aplicações dos métodos computacionais na actividade de projecto em Engenharia Mecânica. Apresentam-se problemas concretos de engenharia que foram abordados durante um estágio realizado no CERN – Centre Européen pour la Recherche Nucléaire, e onde foram utilizados: a) o método dos elementos finitos para cálculo de temperaturas e fluxos de calor e a sua influência sobre os deslocamentos, tensões e deformações que ocorrem numa peça; b) o método híbrido dos elementos finitos/volumes finitos para a discretização das equações de Navier-Stokes e a análise do escoamento de fluidos; c) um algoritmo genético para a obtenção da solução óptima de um problema estrutural. O projecto em engenharia é uma actividade cada vez mais complexa, que requer o uso de ferramentas computacionais sofisticadas tais como os programas ANSYS e MATLAB que foram utilizados no estudo. A criação de modelos numéricos e a análise do seu comportamento com estas ferramentas requer simultaneamente um bom conhecimento dos princípios que estão na base do seu desenvolvimento e uma boa perícia na sua manipulação. Com elas é possível obter soluções quando os constrangimentos do projecto são exigentes e análises detalhadas do comportamento estrutural são necessárias. Neste estudo pretende-se também demonstrar que uma combinação inovadora destas ferramentas pode contribuir para obter aplicações úteis para a actividade de projecto em engenharia.

Estudo numérico de uma técnica de redução de esforços numa estrutura cilíndrica circular

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para obtenção do Grau de Mestre em Engenharia Mecânica

Simulação numérica do escoamento em turbinas colocadas numa falésia

Dissertação para obtenção do Grau de Mestre em Engenharia Mecânica

Estudo numérico do escoamento numa chaminé como ventilador estático

Dissertação para obtenção do Grau de Mestre em Engenharia Mecânica

ESTUDO TEÓRICO DO RESFRIAMENTO COM AR FORÇADO DE FRUTAS DE GEOMETRIAS DIFERENTES

Neste trabalho apresenta-se um estudo numérico para a obtenção das curvas de resfriamento de laranja Valência, (Citrus sinensis O.) e banana prata (Mussa balbisiana Colla), em posições diferentes ao longo do leito e acondicionadas em embalagens com 40% de área efetiva de abertura. Os frutos foram resfriados num sistema de circulação com ar forçado (1.933m³/h), à temperatura de 1°C e 7°C, umidade relativa de 88,4 ± 2,0% e velocidade do ar em torno de 1m/s. Um modelo matemático bidimensional, em coordenadas esferoidais prolato, foi aplicado para predizer a condução de calor transiente dentro das frutas, assumindo-se condição de contorno convectiva na superfície do corpo. As equações geradas foram resolvidas numericamente pelo método de volumes finitos e o coeficiente convectivo de transferência de calor foi obtido aplicando o método de ajuste do erro quadrático mínimo, entre os dados experimentais e numéricos. A análise comparativa das curvas teóricas e experimentais mostrou uma concordância satisfatória, com valores de erro entre 5% e 7%. Notou-se que o coeficiente de transferência de calor varia com a posição dos frutos no leito e que o tempo de resfriamento apresenta uma variação de aproximadamente 38% entre os diferentes pontos. A distribuição espacial de temperatura no interior dos frutos, para três instantes de tempo, demonstrou a existência de um diferencial de temperatura entre o centro e a superfície de aproximadamente 30%. O modelo matemático prediz que a taxa de resfriamento nas bananas é muito mais alta nas pontas que nas outras regiões do fruto. O número de Bi e de Fo mostraram-se adequados e representativos do processo estudado, caracterizando adequadamente o resfriamento de corpos em qualquer relação de área/volume diferente.

Simulação numérica do processo de fusão e solidificação da água no interior de cilindros

Durante a realização de projetos de sistemas de refrigeração assistidos por sistemas de armazenagem de energia, na forma latente, é necessário que se saiba de antemão as taxas de transferência de energia dos mesmos. Na literatura soluções para este problema são escassos e não abordam todas as geometrias utilizadas para o armazenamento de energia latente, inexistindo correlaçõesmatemáticas que permitam aos projetistas prever o comportamento da fusão e da solidificação do material de mudança de fase (MMF) e as taxas de transferência de energia térmica. Na solução destes problemas são usualmente adotadas hipóteses simplificativas para amenizar o grau de complexidade da descrição matemática do problema Uma destas hipóteses constitui-se em desconsiderar o processo advectivo na transferência de calor no MMF, em sua fase líquida. Outra hipótese adotada, quando o MMF é a água pura, consiste em utilizar um comportamento linear para a massa específica ou um comportamento médio, consistindo em um afastamento da realidade. Este trabalho descreve esforços no sentido de aproximar estas soluções cada vez mais da realidade fisica do problema. São apresentadas simulações para a mudança de fase da água pura levando em conta a inversão da massa específica, em geometria polar, na forma transiente, utilizando o método de Volumes Finitos para solução das equações governantes. As soluções apresentadas, otimizadas quanto às malhas espacial e temporal, envolvem condições de contorno de primeira e terceira espécies Como resultado das soluções são apresentados o Nusselt local, ao longo da parede do cilindro além do Nusselt médio. Correlações matemáticas para o número de Nusselt médio durante a fusão e para o volume de MMF fundido, são propostos e apresentados. A evolução geométrica da fronteira de mudança de fase também é avaliada para as duas condições de contorno e diversos números de Rayleigh. Comparações entre o volume fundido considerando apenas o processo difusivo e o volume fundido utilizando a convecção são realizadas. Para a verificação da independ~encia de malha aplica-se o índice de convergência de malha (GCI) . Algumas comparações qualitativas com soluções experimentais disponíveis na literatura são apresentadas.

Determinação dos campos de velocidade e temperatura em ambientes ventilados

No presente trabalho apresenta-se um estudo numérico do escoamento de ar em um ambiente interno tridimensional típico, sujeito a um sistema de ventilação por deslocamento. O estudo numérico consiste na obtenção dos campos de velocidade e temperatura presentes no interior deste ambiente. A determinação dos campos de velocidade e temperatura é necessária na avaliação precisa das condições de conforto térmico, produzidas por um determinado sistema de climatização. Este tipo de estudo não serve somente para determinar a viabilidade da aplicação de um determinado sistema de distribuição de ar, mas através dele também se pode determinar as condições ótimas de uso deste sistema. De forma crescente observa-se o interesse dos projetistas de ar condicionado em aperfeiçoar as suas técnicas de cálculo, de modo que a simulação numérica poderá tornar-se uma técnica usual em pouco tempo. Atualmente duas técnicas estão disponíveis para o estudo do escoamento de ar em ambientes internos: investigação experimental e simulação numérica. Em princípio, medidas diretas dão informações mais realísticas relativas ao escoamento do ar interno, contudo, medições experimentais são muito caras e consomem muito tempo, inviabilizando a sua aplicação caso a caso. Devido às limitações da investigação experimental e à disponibilidade de computadores cada vez mais rápidos, o uso da simulação numérica surge como uma opção considerável nos dias de hoje O método de solução numérica, utilizado no presente trabalho, baseia-se na resolução das equações de conservação de massa, de quantidade de movimento e de energia, em um domínio de cálculo discretizado com o uso do Método dos Volumes Finitos. Os efeitos da turbulência são levados em conta através da utilização do modelo de turbulência K- ε e de um modelo de turbulência simplificado. Resultados bem convergidos para o problema abordado somente são obtidos com a utilização do modelo de turbulência simplificado. As simulações numéricas apresentadas no presente trabalho foram realizadas com o uso do código Fluent. Com base na comparação dos resultados numéricos obtidos no presente trabalho e os resultados experimentais de Yuan et al., 1999a, pode-se concluir que o modelo simplificado é adequado, ao menos provisoriamente, para determinar a distribuição de velocidade e de temperatura no escoamento de ar em ambientes internos.

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.

Otimização de geometrias aerodinâmicas utilizando métodos inversos

O objetivo deste trabalho é a obtenção de uma técnica para a modelagem otimizada de corpos submetidos a fluxos de alta velocidade, como aerofólios em escoamentos transônicos e outras geometrias aerodinâmicas. A técnica é desenvolvida através de expansões em séries de Fourier para um conjunto de equações diferenciais com interrelação com as condições de contorno, sendo uma equação para a parte superior e outra para a parte inferior do aerofólio. O método de integração temporal empregado baseia-se no esquema explícito de Runge-Kutta de 5 estágios para as equações da quantidade de movimento e na relação de estado para a pressão. Para a aproximação espacial adota-se um esquema em volumes finitos no arranjo co-localizado em diferenças centrais. Utiliza-se dissipação artificial para amortecer as frequências de alta ordem do erro na solução das equações linearizadas. A obra apresenta a solução de escoamentos bi e tridimensionais de fluidos compressíveis transônicos em torno de perfis aerodinâmicos. Os testes num´ericos são realizados para as geometrias do NACA 0012 e 0009 e asas tridimensionais usando as equações de Euler, para número de Mach igual a 0.8 e ® = 0o. Os resultados encontrados comparam favoravelmente com os dados experimentais e numéricos disponíveis na literatura.

Modelo matemático de um agitador eletromagnético

O objetivo deste trabalho é estudar os efeitos eletromagnéticos e fluido-dinâmicos induzidos no aço, decorrentes do uso de um agitador eletromagnético. Para tal, foi proposta a construção de um modelo numérico que resolva, de forma acoplada, os problemas de eletromagnetismo e fluido-dinâmica. O modelo numérico do problema eletromagnético, em elementos finitos, foi construído utilizando-se o software Opera-3d/Elektra da Vector Fields. O mesmo foi validado com medidas experimentais de densidade de fluxo magnético feitas na usina. O escoamento decorrente da agitação eletromagnética foi resolvido fazendo-se o acoplamento das forças de Lorentz com as equações de Navier-Stokes. Essas últimas foram resolvidas pelo método de volumes finitos, usando-se o software CFX-4 da AEA Technology. O modelo eletromagnético mostrou que existe um torque máximo dependente da freqüência do campo magnético. Também foi observado que a força magnética aumenta em quatro vezes seu valor, quando a corrente é duplicada. O perfil de escoamento produzido no molde, sob agitação eletromagnética, indica, que as situações de lingotamento testadas, não propiciam o arraste da escória. A velocidade crítica de arraste, determinada via modelo físico, não foi atingida para nenhum caso testado. O modelo fluido-dinâmico e térmico apresentou um aumento do fluxo de calor cedido pelo fluido para a casca solidificada com o uso do agitador eletromagnético. Como conseqüência, observou-se uma queda na temperatura do banho. Também foi observado, que o uso do agitador propicia a remoção de inclusões das camadas mais externas do tarugo. Ao mesmo tempo, notou-se que o uso do agitador aumenta o índice de remoção de inclusões para as duas seções de molde analisadas.

Modelagem dos processos envolvidos nos sistemas de secagem e armazenamento de grãos

A produção de soja é uma das principais atividades econômicas na Região Noroeste do Estado do Rio Grande do Sul. As perdas de produto em condições de comercialização ocasionadas nas atividades de secagem e armazenamento são significativas, justificando a pesquisa e aprimoramento destes processos. Nesta tese foram pesquisados dois problemas: 1. Modelamento matemático dos processos de secagem, utilizando parâmetros conhecidos de soja e 2. Modelamento matemático do problema de aeração para o cálculo da distribuição da pressão e da velocidade do ar na massa de grãos em unidades de armazenamento de soja. No problema de secagem foi desenvolvido um sistema composto de quatro equações diferenciais parciais hiperbólicas acopladas não-lineares, que descreve o comportamento da temperatura e do teor de umidade do ar e dos grãos em função do tempo. Para resolver o sistema foram utilizados os métodos das diferenças finitas (p. ex., métodos de MacCormack e Crank- Nicolson.) e o método dos volumes finitos. A análise dos resultados permitiu recomendar o método mais adequado para cada tipo do problema. Para determinação da intensidade do fluxo de massa e de calor foram utilizados os dados experimentais de camada fina obtidos da literatura e complementados com dados experimentais desta tese. Foi desenvolvido um equipamento para obtenção das curvas de secagem de grãos em secador de leito fixo, a fim de identificar o modelo para secagem em camada espessa. A comparação entre os resultados experimentais e das simulações numéricas mostrou que o modelo descreve razoavelmente a dinâmica de secagem No problema de aeração foi desenvolvido um modelo matemático que descreve o escoamento do ar em sistemas de armazenamento de grãos, baseado em relações experimentais entre velocidade e gradiente de pressão. Para resolver o problema de aeração foi utilizado o método dos elementos finitos e desenvolvido um programa computacional. Um teste realizado com o programa mostrou que os resultados da solução numérica convergem para uma solução analítica conhecida. As simulações realizadas mostraram que o programa computacional pode ser usado como instrumento auxiliar para o projeto de silos, possibilitando o cálculo e a visualização gráfica da distribuição das pressões e das linhas de corrente em diferentes seções do armazém.

Simulação numérica da combustão turbulenta de gás natural em câmara cilíndrica

O presente trabalho apresenta uma modelagem detalhada de processos de combustao turbulentos para um jato concentrico de combustıvel e ar. A modelagem é fundamentada nas equacões de conservacão de massa, de quantidade de movimento, de energia e de espécies quımicas. A turbulencia é resolvida pela utilizacão do modelo k- padrão. Dois modelos de reacões quımicas são apresentados. O modelo SCRS – Simple Chemically-Reacting Systems, que assume taxas instantâneas de reacões quımicas. Também é abordado o modelo E-A – Eddy Breakup - Arrhenius, que assume taxas finitas de reacões quımicas. A radiacão térmica, fenômeno de grande importância devido as altas temperaturas alcancadas em processos de combustão,é modelada através do Método das Zonas. O modelo da soma ponderada de gases cinzas – WSGGM,é usado para determinar o espectro de emissão e absorcão dos gases no processo. Para a solucão destas equacões diferenciais, juntamente com os modelos de turbulência, de reaçõoes químicas e radiação térmica, faz-se o uso do Método dos Volumes Finitos. Para validar a modelagem apresentada resolve-se o processo de combustão em uma câmara cilíndrica. A câmara de combustão usada áa mesma abordada no First Workshop on Aerodynamics of Steady State Combustion Chambers and Furnaces, organizado pela ERCORTAC - European Research Community On Flow Turbulence And Combustion, em outubro de 1994, que apresenta dados experimentais de temperatura e concentração das espécies químicas para várias posições de interesse no interior da câmara. Utiliza-se o gás natural como combustível e o ar atmosférico como oxidante. O processo de combustão sem pré-mistura é resolvido para a condição de excesso de combustível de 5 % para ambos os modelos, onde o gás natural é injetado por um duto circular central, e o ar atmosférico por um orifício anular externo a esse duto, no mesmo plano Uma reação química não estagiada é assumida para o modelo SCRS. Para o modelo E-A duas situações são resolvidas: combustão não estagiada, com uma etapa global de reação química; e reação quımica estagiada, com duas etapas globais. Os resultados obtidos com o modelo SCRS para a distribuição de temperaturas, em termos de tendências gerais, são razoáveis. Já as concentrações de espécies químicas não apresentam dados satisfatórios para este modelo. Para o modelo E-A os resultados apresentam boa concordância com os dados experimentais, principalmente para a situação em que o processo de combustão é assumido em duas etapas globais. ´E analisado em detalhe o papel desempenhado pela transferencia de calor por radiacao, com meio participante. Para melhor verificar as trocas de calor, assume-se uma camara de combustao cilındrica com paredes d’agua. A injecao do combustıvel e do oxidante e feita atraves de um queimador central, semelhante ao usado para validar a modelagem, porem com dois orifıcios concentricos para injecao de combustıvel. Nesta situação o efeito do turbilhonamento (swril), assumido como 20 % da velocidade axial de entrada, sobre a injecao de ar e computado atraves da condicao contorno da equacao de conservacao da quantidade de movimento angular. Nesta fase apenas o modelo E-A, com duas etapas globais de reacoes quımicas, e considerado, ja que o mesmo apresenta os melhores resultados. O processo de combustao e simulado com e sem a presenca da radiacao termica. Verifica-se que a presenca da radiacao termica homogeneiza a temperatura dos gases no interior da camara. Com isso verifica-se tambem alterações nas taxas de reacoes quımicas, modificando a magnitude das fracoes das especies quımicas Quando a radiacao termica e considerada efeitos de extinção local da chama sao verificados nas regioes de temperaturas mais altas, diminuindo o consumo de oxigenio e aumentando a producao de monoxido de carbono, caracterizando assim uma combustao incompleta. Em algumas situacoes tem-se uma variacao de temperatura de ate 500 K, a montante da chama. A radiacao termica tambem aumenta a taxa de transferencia de calor dos gases quentes para as paredes da camara, e desta para o seu exterior. Com os resultados obtidos a partir desta modelagem e possıvel determinar o perfil da zona de combustao, a distribuicao de concentracoes de especies quımicas, o campo de velocidades e as taxas de transferencia de calor para as paredes da camara de combustao, total, por conveccao superficial e por radiacao. Estes resultados sao de extrema importancia para prever a performance de camaras de combustao, assim como auxiliar na sua otimizacao.

Análise numérica de uma torre de resfriamento de grande porte

Este trabalho tem como objetivo principal determinar a troca de calor em uma torre de resfriamento de grande porte, do tipo seca. A força motriz deste sistema é o empuxo resultante da diferença de temperatura entre o ar interior e exterior da torre. A partir de uma certa velocidade o vento externo também desempenha um papel importante. A convecção natural ocorre em regime turbulento, devido às grandes dimensões. A modelagem numérica é feita pelo Método dos Volumes Finitos com malhas criadas sobre Elementos Finitos, através dos softwares comerciais Fluent e CFX. Os resultados são obtidos para diferentes potências dissipadas no trocador de calor. A partir dos resultados são realizados estudos sobre a viabilidade de mudanças na geometria da torre, a fim de se obter melhorias no escoamento, bem como uma troca de calor mais eficiente. Faz-se ainda uma regressão a fim de obter correlações entre a potência dissipada no trocador de calor e a vazão mássica de ar necessária para a troca. Estas correlaçõoes são de utilidade, para aplicações posteriores, na implantação de um simulador para o processo de troca de calor nesta torre. O trocador de calor também é enfocado a fim de se determinar o coeficiente global de troca térmica, a fim de disponibilizá-lo para aplicação em outros trabalhos

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.