Simulação computacional do comportamento de estruturas de aço sob incêndio

A análise do comportamento estrutural sob incêndio constitui uma parte importante da engenharia de proteção contra incêndio, especialmente no caso de estruturas de aço, por sua alta condutividade térmica e relativa esbeltez das seções. As dificuldades econômicas e práticas, associadas à avaliação do comportamento estrutural, por meio de ensaios em escala real, têm estimulado o desenvolvimento e uso de métodos de simulação numérica. Esta tese trata da simulação numérica do comportamento de estruturas de aço sob condições de incêndio e se divide em três partes. As duas primeiras partes foram desenvolvidas na Universidade de Liége, na Bélgica, usando-se o programa SAFIR como ferramenta numérica. A terceira parte foi desenvolvida de forma independente, na Universidade Federal do Rio Grande do Sul. Na primeira parte, é feito um estudo comparativo entre o uso de elementos finitos de viga e de casca, na modelagem de vigas simplesmente apoiadas, sujeitas a flambagem lateral por torção. Os esforços de torção, presentes no caso de flambagem lateral, podem levar a uma plastificação da seção transversal e à conseqüente redução da rigidez à torção da seção. Da mesma forma, a degradação das propriedades do material, com o aumento da temperatura, leva à redução da rigidez à torção Havia dúvidas se o modelo com elementos de viga, com uma rigidez à torção constante, poderia fornecer uma resposta aceitável. O estudo mostrou que uma resposta com boa precisão pode ser conseguida, usando-se elementos de viga, desde que o módulo de elasticidade transversal seja ajustado para refletir seu valor correspondente à temperatura de falha. Isso implica um processo iterativo, uma vez que a temperatura de falha não é previamente conhecida. Por outro lado, a degradação da rigidez à torção, por efeitos de plastificação, pode ser ignorada. Na segunda parte, é feita a comparação entre as modelagens bidimensional e tridimensional, de galpões industriais de um andar, sob incêndio. Comumente, a estrutura de galpões industriais é composta por pórticos-tipo, dispostos em paralelo. A análise desses galpões é comumente feita pela simulação no plano do pórtico de aço ou, simplesmente, da treliça da cobertura Na análise bidimensional, importantes efeitos fora do plano são ignorados, como a redistribuição de esforços devido à degradação do material ou à expansão térmica, ou instabilidade lateral dos elementos. A importância desses efeitos e a adequabilidade do modelo 2D para representar o comportamento real são discutidas. Na terceira parte, um modelo numérico para a simulação tridimensional do comportamento de estruturas de aço sob incêndio é apresentado. O modelo é baseado no conceito de rótulas plásticas generalizadas, com modificações para melhor representar a formação e expansão da plastificação no elemento. A descrição cinemática adotada permite obter bons resultados, mesmo com o uso de poucos elementos. A determinação do vetor de esforços internos no elemento, incluindo os efeitos da temperatura, é detalhada. O procedimento foi validado por comparação com ensaios e com modelos numéricos mais sofisticados, como o programa SAFIR. Os resultados demonstram que o procedimento proposto pode ser usado como uma forma alternativa de análise 3D de estruturas sob incêndio, com precisão razoável e baixo esforço computacional.

Simulador em regime transiente de alguns equipamentos de um ciclo de Rankine regenerativo

Este trabalho apresenta um modelamento computacional dos diversos componentes de um ciclo de Rankine regenerativo em operação transiente, utilizando as leis fundamentais de balanço de massa e energia, além de algumas equações constitutivas. O cálculo das propriedades termodinâmicas foi realizado através de rotinas desenvolvidas no Grupo de Estudos Térmicos e Energéticos (GESTE) da UFRGS. Os componentes com comportamento inercial são o condensador e o tanque de água de alimentação. Os demais respondem instantaneamente. O método de solução das equações algébricas é seqüencial e o avanço no tempo é totalmente implícito. Primeiramente o sistema é colocado em regime permanente, a partir do qual é introduzida uma mudança de carga na turbina em um determinado intervalo de tempo. Algumas variáveis prescritas foram obtidas a partir de dados reais de uma usina termelétrica. A partir desse modelamento foi possível investigar o comportamento dinâmico do condensador e do tanque de alimentação e de que maneira estes efeitos eram percebidos pelos demais componentes. Os resultados obtidos demonstram que ambos possuem um comportamento dinâmico, sendo que o condensador apresenta uma constante de tempo menor que o tanque de alimentação. O condensador, por rejeitar o calor necessário, acaba por perceber rapidamente as oscilações introduzidas no ciclo Além disto é apresentado um método para determinação de coeficientes globais de troca térmica em condensadores e regeneradores. Os resultados obtidos neste trabalho foram utilizados na modelagem desses componentes durante a construção de um simulador em regime permanente da UTE Presidente Médici - Fase B. Também foi realizada uma pesquisa sobre os diversos programas simuladores de plantas de potência, assim como uma discussão das vantagens introduzidas por simuladores dinâmicos e em regime permanente nas decisões de engenharia.

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.

Essays in recursive equilibrium and market selection hypothesis

Este trabalho consiste em estudar modelos incluindo agentes com informação completa e incompleta sobre o ambiente econômico. Prova-se a existência de equilíbrio em que esses dois agentes coexistem e sob, algumas condições, obtêm-se que esse equilíbrio é recursivo e contínuo, ou seja, pode ser implementado por uma função contínua de transição que relaciona as variáveis de equilíbrio entre dois períodos consecutivos. Mostra-se, sob algumas hipóteses, que em equilíbrios recursivos contínuos, os agentes que cometem erros persistentes nas antecipações dos preços de equilíbrio são eliminados do mercado. Finalmente, exibimos diversos exemplos numéricos, no caso de mercados incompletos e informação completa, em que os agentes com expectativas racionais são eliminados do mercado. Usam-se métodos numéricos alternativos que possibilitam computar um equilíbrio em modelos com agentes heterogêneos.

O uso da transformada Wavelet bi-dimensional no conceito do espaço de escala

O processamento de imagens tem sido amplamente utilizado para duas tarefas. Uma delas é o realce de imagens para a posterior visualização e a outra tarefa é a extração de informações para análise de imagens. Este trabalho apresenta um estudo sobre duas teorias multi-escalas chamadas de espaço de escala e transformada wavelet, que são utilizadas para a extração de informações de imagens. Um dos aspectos do espaço de escalas que tem sido amplamente discutido por diversos autores é a sua base (originalmente a gaussiana). Tem se buscado saber se a base gaussiana é a melhor, ou para quais casos ela é a melhor. Além disto, os autores têm procurado desenvolver novas bases, com características diferentes das pertencentes à gaussiana. De posse destas novas bases, pode-se compará-las com a base gaussiana e verificar onde cada base apresenta melhor desempenho. Neste trabalho, foi usada (i) a teoria do espaço de escalas, (ii) a teoria da transformada wavelet e (iii) as relações entre elas, a fim de gerar um método para criar novas bases para o espaço de escalas a partir de funções wavelets. O espaço de escala é um caso particular da transformada wavelet quando se usam as derivadas da gaussiana para gerar os operadores do espaço de escala. É com base nesta característica que se propôs o novo método apresentado. Além disto, o método proposto usa a resposta em freqüência das funções analisadas. As funções bases do espaço de escala possuem resposta em freqüência do tipo passa baixas. As funções wavelets, por sua vez, possuem resposta do tipo passa faixas Para obter as funções bases a partir das wavelets faz-se a integração numérica destas funções até que sua resposta em freqüência seja do tipo passa baixas. Algumas das funções wavelets estudadas não possuem definição para o caso bi-dimensional, por isso foram estudadas três formas de gerar funções bi-dimensionais a partir de funções unidimensionais. Com o uso deste método foi possível gerar dez novas bases para o espaço de escala. Algumas dessas novas bases apresentaram comportamento semelhante ao apresentado pela base gaussiana, outras não. Para as funções que não apresentaram o comportamento esperado, quando usadas com as definições originais dos operadores do espaço de escala, foram propostas novas definições para tais operadores (detectores de borda e bolha). Também foram geradas duas aplicações com o espaço de escala, sendo elas um algoritmo para a segmentação de cavidades cardíacas e um algoritmo para segmentação e contagem de células sanguíneas.

Metodismo e setores populares: análise da proposta de trabalho do Instituto Central do Povo (Rio de Janeiro, 1906-1930)

A prática educativa protestante no Brasi I se ex pressou principalmente através de colégios destinados as camadas dominantes da sociedade. No entanto, existiram também outras escolas, para atendimento a alunos provenientes das classes populares. No caso específico dos metodistas, foram cria - das escolas paroquiais e, no Rio de Janeiro, na zona portuária, surgiu o Instituto Central do Povo, com ativida - des ligadas a educação, saúde e higiene, trabalho, lazer, além da pregação religiosa. Este tipo de proposta de trabalho foi liderado por setores da igreja metodista que revelavam especial in teresse pelas questões sociais. De fato, desde sua origem na Inglaterra existem evidências de que os metodistas se preocupavam com as condições sociais a que o povo da epoca estava submetido, como se pode ver nos discursos de seus líderes e mesmo na sua prática. Posteriormente, ao transplantar-se o meto - dismo para a América, persistiram as discussões acerca des sas questoes, inclusive provocando divisões internas. Embora a orientação predominante, trazida para o Brasil pelos missionários, tenha sido a influenciada p~ 10 pietismo e pela fi losofia 1 iberal, em íntima relação com o processo de consolidação do capitalismo na sociedade norte-americana durante o Século XIX, é necessário reconhecer entre os metodistas a existência de grupos que defenderam o envolvimento dos cristãos e da igreja nas lu tas sociais, em defesa dos setores populares. Dentro do quadro do protestantismo norte-americano esses grupos se fil iam à corrente conhecida como Evangelho Social.A tentativa de atuação dos missionários metodis tas, 1 igados a esta tendência, que vieram ao Brasil teve, porem, alcance limitado na medida em que esbarrou em di - versos obstáculos: nos fundamentos teóricos sobre os quais as atividades foram concebidas, nas dificuldades de inser ção num contexto social diferente e na própria situação de inferioridade numérica dessa corrente dentro da Igreja.

Desenvolvimento da técnica para determinação dos parâmetros térmicos de contorno aplicados ao projeto de peças forjadas

O presente trabalho descreve os resultados da pesquisa que tem por objetivo a determinação de parâmetros térmicos de um corpo sólido: o valor médio da emissividade total e o coeficiente médio de transferência de calor por convecção. Foi desenvolvido o modelo matemático do processo de resfriamento de um corpo forjado. Foram examinados os métodos da determinação dos parâmetros mencionados, que estão descritos na literatura, e também foram analisadas as possibilidades de aproveitá-los no processo de forjamento. No meio de métodos analisados, o método do Instituto IBF (Institut für Bildsame Formgebung, IBF, RWTH, Aachen, Alemanha), foi considerado o mais apropriado tendo a metodologia deste trabalho seguindo a mesma sistemática. Aproveitando a metodologia do IBF para a determinação dos parâmetros térmicos mencionados e o modelo matemático do resfriamento de um corpo forjado, foi desenvolvida uma nova técnica para determinar o valor médio da emissividade total e o coeficiente médio de transferência de calor por convecção de um corpo sólido. Os materiais analisados foram o aço ABNT 1045 e a liga de alumínio AA 6061. Foram feitas as comparações dos resultados obtidos pela técnica desenvolvida com o método do IBF.

Alguns aspectos numéricos e teóricos das equações de Navier-Stokes na modelagem do escoamento em torno de um vórtice

Neste trabalho desenvolvemos uma metodologia numérica para a solução do escoamento em torno de um vórtice. Como a análise completa deste tipo de fluxo não é uma tarefa fácil, simplificações quanto ao escoamento e ao método numérico são necessárias. Também investigamos o comportamento das soluções das equações governantes (Navier-Stokes) quando o tempo tende ao infinito. Nesse sentido, dividimos este trabalho em duas partes: uma numérica e outra analítica. Com o intuito de resolver numericamente o problema, adotamos o método de diferenças finitas baseado na formulação incompressível das equações governantes. O método numérico para integrar essas equações é baseado no esquema de Runge- Kutta com três estágios. Os resultados numéricos são obtidos para cinco planos bidimensionais de um vórtice com números de Reynolds variando entre 1000 e 10000. Na parte analítica estudamos taxas de decaimento das soluções das equações de Navier-Stokes quando os dados iniciais são conhecidos. Também estimamos as taxas de decaimento para algumas derivadas das soluções na norma L2 e comparamos com as taxas correspondentes da solução da equação do calor.

Mecanismos técnico-institucionais para a sustentabilidade da drenagem urbana

A drenagem urbana, especialmente nos grandes centros, tem sido efetuada de forma não sustentável com contaminação e alterações no regime de escoamento superficial devido à impermeabilização de superfícies, e deposição de resíduos sólidos. Estas condições representam ameaça considerável ao homem e ao ecossistema do corpo receptor. O emprego de práticas sustentáveis aparece como caminho a ser perseguido, encontrando no Desenvolvimento Urbano de Baixo Impacto (Low Impact Development, LID) o conjunto de técnicas que mais se aproxima desta meta. LID objetiva atingir paisagens com funções hidrológicas, apresentando comportamento mais similar ao natural, por controlar não somente o pico de vazões, mas volume, freqüência/duração, além de qualidade dos escoamentos pluviais. Este tipo de desenvolvimento atua recuperando a capacidade de infiltração das superfícies urbanas, além de estimular o reuso de água, reduzindo os impactos, com ganhos econômicos e paisagísticos em comparação ao controle efetuado pelos métodos tradicionais de controle por condutos e mesmo detenções. A utilização desta nova tecnologia se aplica à implantação de novos desenvolvimentos e redesenvolvimentos, a priori, apresentando ainda vantagens para implantação destas em empreendimentos antigos com relação a métodos tradicionais Buscou-se, por intermédio deste estudo, avaliar os potenciais mecanismos técnicoinstitucionais que possam ser empregados na realidade brasileira, em especial em Porto Alegre, vislumbrando a implementação da sustentabilidade, assim como, avaliar as respostas obtidas da simulação numérica da implantação das técnicas de LID a um condomínio hipotético. Para tanto, fez-se necessário aplicar metodologia diferente da proposta pelos manuais existentes, em virtude da insuficiência destes em representar o comportamento dos dispositivos de controle do escoamento. Os resultados obtidos confirmaram a necessidade de revisão nos mecanismos para controle da drenagem vigentes, bem como, as vantagens de aplicação de técnicas que apresentam abordagem mais integradora e de escala menor.

Transferência de calor e controle de temperatura em tubos capilares utilizados em eletroforese capilar

Este trabalho tem por objetivo estudar a transferência de calor em tubos capilares cilíndricos utilizados na técnica de separação de moléculas denominada Eletroforese Capilar. Esta técnica é usada, por exemplo, na análise de biomoléculas e no sequenciamento de DNA, onde o controle da temperatura está diretamente ligado ao desempenho destes métodos e à qualidade dos resultados. Para empregar esta técnica, tensões elétricas da ordem de 20 kV são aplicadas entre as extremidades dos tubos capilares, que possuem normalmente 50 cm de comprimento, 350 µm de diâmetro externo e 50 µm de diâmetro interno, preenchidos por uma solução aquosa. Tais tensões geram uma corrente elétrica na solução, provocando aquecimento distribuído por Efeito Joule. Os tubos capilares são construídos em quartzo amorfo e protegidos por uma camada de material polimérico (poliimida). Para implementar o controle da temperatura, os tubos capilares são colocados em contato com um fluido de resfriamento. Num primeiro momento, os estudos são realizados por simulação numérica, empregando o Método dos Volumes Finitos em rotinas escritas em FORTRAN. São simulados casos onde os tubos são recobertos por camadas cilíndricas de materiais com uma condutividade térmica relativamente boa, com o objetivo de aumentar a superfície de troca de calor com o fluido de resfriamento. Como resultado, obtêm-se curvas da temperatura no centro dos tubos capilares em função do coeficiente de transferência de calor por convecção. Um caso de interesse é quando os tubos capilares são posicionados excentricamente ao recobrimento cilíndrico Num segundo momento, é utilizado o software de simulação numérica ANSYS CFX®, onde é simulado o resfriamento dos mesmos tubos capilares expostos a um escoamento transversal de ar a 15°C. Neste caso, também são aplicados os recobrimentos cilíndricos e, além disso, opta-se por simular o resfriamento de um arranjo de vários tubos (sistema multicapilar) dispostos entre placas de vidro, no formato de um sanduíche. Como resultados mais importantes salientam-se: a) o aumento do raio do recobrimento resulta no aumento da transferência de calor, fazendo com que a temperatura no núcleo do capilar fique estacionada em valores baixos que não comprometem as separações/análises; b) chegou-se a um valor de raio crítico da ordem de 10 mm para a condição de operação mais típicas na área da Eletroforese Capilar; c) as montagens com o tubo capilar concêntrico e excêntrico ao recobrimento não apresentam diferenças significativas no perfil de temperatura da solução tampão; e finalmente d) observa-se que o uso de duas placas de material dielétrico com os capilares posicionados em forma de sanduíche entre elas permite uma eficiente dissipação do calor gerado na solução tampão.

Solução da equação de difusão-advecção para uma CLP não-homogênea e não-estacionária pelo método GILTT

Neste trabalho é apresentada a solução da equação de difusão-advecção transiente para simular a dispersão de poluentes na Camada Limite Planetária. A solução é obtida através do método analítico GILTT (Generalized Integral Laplace Transform Technique) e da técnica de inversão numérica da quadratura de Gauss. A validação da solução é comprovada utilizando as concentraçãos obtidas a partir do modelo com as obtidas experimentalmente pelo Experimento de Copenhagen. Nesta comparação foram utilizados os perfis de vento potencial e logaritmo e os parâmetros de turbulência propostos por Degrazia et al (1997) [19] e (2002) [17]. Os melhores resultados foram obtidos utilizando o perfil de vento potencial e o coeficiente de difusão propostos por Degrazia et al (1997). A influência da velocidade vertical é mostrada através do comportamento das concentrações de poluentes na pluma. Além disso, as velocidades verticais e longitudinais geradas pelo Large Eddy Simulation (LES) foram colocadas no modelo para poder simular uma camada limite turbulenta mais realística, a qual apresentou resultados satisfatórios quando comparados com os disponíveis na literatura.

Leilões de valor quase comum e o regime de partilha de produção brasileiro: uma abordagem teórica

Neste trabalho, estudamos o regime de partilha de produção brasileiro, instituído pela Lei N$^o$ 12.351, para exploração de petróleo através de uma abordagem teórica. Alguns artigos da literatura de leilões de valor quase comum foram estudados aqui. Além disso, desenvolvemos um modelo de partilha de produção com informação assimétrica a fim de capturar algumas características do modelo de partilha brasileiro. Através de solução numérica, fazemos uma análise das estratégias dos participantes e dos ganhos esperados.

Contribuição ao ensino e à prática da engenharia econômica

Esta dissertação versa sobre um conjunto de técnicas da análise numérica, baseadas nos cálculos de juros compostos, normalmente conhecidas pelo nome de Engenharia Econômica. A engenharia econômica tem cerca de um século de existência e, por sua importância, constitui disciplina obrigatória dos cursos de graduação e especialização em administração desta escola, sendo ensinada também em escolas de engenharia e outras escolas de administração, tanto no Brasil como no exterior. O autor não pretendeu escrever um tratado sobre a matéria mas sim contribuir com sua experiência para atingir dois objetivos.

Análise de sensibilidade aplicada a uma planta térmica de geração de potência

Este trabalho tem o objetivo de realizar uma análise de sensibilidade dos parâmetros operacionais e de projeto de uma usina termelétrica. A análise de sensibilidade serve para indicar qual a importância dessas grandezas no resultado do sistema. A planta modelada é a da Usina Termelétrica Presidente Médici – Fase B, localizada em Candiota, RS, onde cada módulo é capaz de produzir até 160 MW de potência. Para esta planta é apresentado o equacionamento de cada componente do sistema. Destacam-se duas modelagens para a turbina a vapor, usadas para a simulação de seu comportamento em carga parcial, o da Elipse de Stodola e o de Schegliáiev. Ambos os modelos são comparados. Também são apresentados os modos de controle de carga por válvula de estrangulamento, por válvula de bocais e por pressão deslizante. A análise de sensibilidade é feita por meio de três métodos distintos, cujos resultados são comparados. O Método Diferencial utiliza as derivadas parciais para estimar a sensibilidade. O Método de Monte Carlo realiza uma série de avaliações do sistema, com os dados de entrada gerados randomicamente em cada avaliação e análise estatística dos resultados para avaliar a sensibilidade das respostas. É desenvolvido neste trabalho o Método da Transformada de Fourier, que gera os dados de forma senoidal e utiliza a transformada de Fourier para medir as sensibilidades. A medida da sensibilidade é feita através de dois índices. O índice de importância é a derivada parcial de um resultado em relação a um parâmetro, adimensionalizado pelas médias das variáveis, e o índice de sensibilidade é a composição fracional da variância de um resultado. Como conclusão, observa-se que os modelos para operação em carga parcial da turbina de Schegliáiev e de Stodola produzem praticamente os mesmos resultados. Os três métodos de análise de sensibilidade conduzem aos mesmos índices de sensibilidade no caso estudado. A ordenação da importância dos parâmetros estudados é idêntica quando se usa tanto o índice de importância quanto o de sensibilidade, mesmo seus valores sejam diferentes para cada parâmetros.

Análise de escoamentos incompressíveis utilizando simulação de grandes escalas e adaptação de malhas

No presente estudo, são apresentadas soluções numéricas de problemas de Engenharia, na área de Dinâmica dos Fluidos Computacional, envolvendo fluidos viscosos, em escoamentos incompressíveis, isotérmicos e não isotérmicos, em regime laminar e turbulento, podendo envolver transporte de massa. Os principais objetivos deste trabalho são a formulação e a aplicação de uma estratégia de adaptação automática de malhas e a inclusão de modelos de viscosidade turbulenta, integrados com um algoritmo utilizado para simular escoamentos de fluidos viscosos bi e tridimensionais, no contexto de malhas não estruturadas. O estudo é dirigido no sentido de aumentar o conhecimento a respeito das estruturas de escoamentos turbulentos e de estudar os efeitos físicos no transporte de quantidades escalares propiciando, através de técnicas de adaptação automática de malhas, a obtenção de soluções numéricas precisas a um custo computacional otimizado. As equações de conservação de massa, de balanço de quantidade de movimento e de quantidade escalar filtradas são utilizadas para simular as grandes escalas de escoamentos turbulentos e, para representar as escalas submalha, são utilizados dois modelos de viscosidade turbulenta: o modelo de Smagorinsky clássico e o modelo dinâmico. Para obter soluções numéricas com precisão, é desenvolvida e implementada uma estratégia de adaptação automática de malhas, a qual é realizada simultaneamente e interativamente com a obtenção da solução. O estudo do comportamento da solução numérica é fundamentado em indicadores de erro, com o propósito de mapear as regiões onde certos fenômenos físicos do escoamento ocorrem com maior intensidade e de aplicar nestas regiões um esquema de adaptação de malhas. A adaptação é constituída por processos de refinamento/desrefinamento e por um processo de suavização laplaciana. Os procedimentos para a implementação dos modelos de viscosidade turbulenta e a estratégia de adaptação automática de malhas são incorporados ao código computacional de elementos finitos tridimensionais, o qual utiliza elementos tetraédricos lineares. Aplicações de escoamentos de fluidos viscosos, incompressíveis, isotérmicos e não isotérmicos em regime laminar e turbulento são simuladas e os resultados são apresentados e comparados com os obtidos numérica ou experimentalmente por outros autores.