El objetivo del siguiente trabajo es establecer el comportamiento de un perfil aerodinámico cuando éste se acerca al suelo, además de buscar otros fenómenos que se producen en dicho acercamiento (como en la estabilidad por ejemplo). Para alcanzar los objetivos y luego de una revisión de las bases teóricas del tema, se preparó un programa de simulación basado en el método de vórtices discretos en fluidos bidimensionales cuyos resultados permitieran observar como variaban estos fenómenos en relación a la cercanía al suelo. Para comprobar la fidelidad de los resultados obtenidos se compararon algunos de éstos con conceptos teóricos ideales y con resultados experimentales de túneles de viento, además de fuentes diversas sobre experimentos realizados sobre el fenómeno WIG (Wing in Ground). Luego de poner en práctica conjunta todos estos conceptos se pudo llegar a diversas conclusiones, las más importantes obviamente relacionadas con los objetivos de trabajo:  La sustentación de un perfil aumenta cuando éste se acerca al suelo (a una distancia de 50% de la cuerda del perfil puede haber un aumento de aprox 30% en la sustentación).  Este perfil pierde estabilidad cuando se acerca al suelo (aumenta el momento de giro generado por las fuerzas de sustentación) Al tener gran cantidad de datos se pudo llegar a otras conclusiones como por ejemplo:  El efecto suelo es proporcional a la longitud de cuerda.  El efecto suelo tiene una variación geométrica respecto de la cercanía al suelo.  Un aeroplano puede sacar partido de este fenómeno para lograr performances mayores en alcance o capacidad de carga. A manera de sugerencia anexo en este trabajo algunos posibles temas que pueden ser desarrollados tomando éste como partida en cualquier parte, sea por aerodinámica, estudio del efecto suelo o la programación de la simulación en si. / Tesis

Since the turn of the century, the talk about the limited reserve of fossil fuels and the effects of their burning on our climate has become a major topic of the media. The evidence is staggering and as a consequence most of our world’s countries have started an energy transition. The main goal is to get away from fossil fuels and use “renewable energies”, so called because the resources are constantly renewed and compared to fossil fuels seem infinite. Energies such as solar, wind, biomass and geothermal are examples of such sources of renewables. Denmark is the current leader in wind generated electricity, California and Spain are showing how to harness the power of the Sun, and all those efforts to generate more with renewables has to be matched with the effort to make those solutions more efficient and more attractive to other countries who still view fossil fuels as the easy solution and keep on using them. With the knowledge available now renewables it feels for some that burning fossil fuels is a primitive solution. Nonetheless it should be the duty of engineers to enlarge and better that knowledge for everyone to use. Now more specifically about wind energy. Humans have harvested the energy in the wind for more than 2000 years (the Persians used windmills around 200 B.C.) and with time our technology has improved. They have designed incredible new machines such as the Savonius and Darreus type turbines and their knowledge of fluid dynamics has permitted the implementation of new streamlined blades that harvest more energy from the wind. Albert Betz has shown that a maximum of 59% was the limit for the efficiency of a wind-turbine. They have been getting closer to this number with the years but there is still room for improvement on certain types of turbines. Vertical axis wind turbines (or VAWTs) have always been considered not as suitable for energy production as horizontal type wind turbines. It is true because not all blades are exposed to the wind at all times (like in a horizontal axis wind turbine), but new studies have proved that streamlining the blades a certain way and adding a wing like thickness to them improved the overall efficiency of the turbine. Knowing that and considering that HAWTs are significantly cheaper to produce and maintain than HAWTS, it makes them a more viable solution notably for local decentralized production in isolated areas of the world.That leads to Peru. Peru is a fast growing still yet a 3rd world country. Its potential for renewable energies production (especially wind energy) is tremendous, yet the great amount of gas and oil available in the underground and them coming at a cheap price does not encourage the government to subsidize renewables. It leaves Peru dependent on foreign investments to develop this sector and takes away a great opportunity to forego its energetical transition and get ahead of competition in South America. Some projects have surfaced notably in northern Peru, in the Trujillo region, but there are few compared to the mega industry of oil and gas. This thesis paper has for goal to further the knowledge of wind turbines in the context of hoping to change Peru’s view on their use and also for the world to use as a database for further research and other works. / Tesis

O presente trabalho simula processos de conformação mecânica em domínios bidimensionais. Para tanto, considera-se o material como sendo um fluido não newtoniana e utiliza-se a formulação da mecânica dos fluidos com função de penalidade. O método empregado para realizar a discretização espacial é o método dos elementos finitos com integração reduzida do termo da função de penalidade. Para discretizar o tempo usa-se diferenças finitas com o esquema de Crank-Nicolson. E apresentado também o método Lagrangiano-Euleriano arbitrário para descrever o movimento da malha de elementos finitos. / Simulations o f metal forming processes in two dimensional domains are presented in this work. The material is considered as a non newtonian flu id and the flow formulation with the penalty function approach is used. Spatial discretization is accomplished using the finite element method with reduced integration for the pressure term. Finite differences is used for time integration with a Crank-Nicolson scheme. The arbitrary Lagrangian-Eulerian scheme is employed to describe the motion of the finite element mesh.

Este trabalho tem por objetivo apresentar um modelo numérico para simular as vibrações induzidas por escoamentos isotérmicos em estruturas imersas em fluidos quase­ incompressiveis. Utiliza-se o método dos elementos finitos e considera-se o acoplamento do fluido com a estrutura imersa. Admite-se que as estruturas sao corpos rígidos com restrições elásticas para giros e deslocamentos. Para a análise do escoamento do fluido é utilizado um método explicito de dois passos de Taylor-Galerkin. Para a análise da estrutura imersa é utilizado o método de Newmark. O fluido e a estrutura são acoplados através da imposição de certas condições de compatibilidade e de equilíbrio, e da consideraçao da contribuição do fluido na massa e no amortecimento da estrutura. sao apresentados exemplos de aplicaçao à análise do escoamento em torno de um cilindro fixo e do escoamento em torno de um cilindro oscilante. / A numerical model to simulate isothermal flow induced vibrations of immersed structures in quasi-incompressible fluids, is presented in this work. Finito cle ent mcthod is uscd and fluid-structure interaction is considered. It is assumed that immersed structures are rigid bodies with elastic restrains for rotations and displacements. For the fluid flow analysis is used an explicit two steps Taylor-Galerkin method. For the immersed structure is used Newmark's method. The fluid and the structure are coupled through the imposition of compatibility and equilibrium conditions, and the consideration of the fluid contribution to the mass and the darnping of the structure. Examples of aplication to the analysis of flow around one fixed cilinder and of flow around one oscillating cilinder, are presented.

In questo lavoro vengono presentati dei modelli di dispersione analitici, per la qualità dell'aria, basati sulle equazioni di diffusione ed avvezione multidimensionali, in geometria Cartesiana. Questi modelli sono validi per un contaminante passivo e per condizioni di turbolenza omogenea e uniforme e velocità dei vento uniforme, che equivale ad emissioni a quote elevate che si verificano in condizioni di stabilità intermedia e in assenza di di una forte forza di "buoyancy". Le soluzioni sono ottenute con l'uso dei metodo della Trasformata Integralle Generalizzata e sono analitiche nel senso che nessuna approssimazione e fatta nella derivazione. L'idea di base dei metodo si sviluppa nelle fasi seguenti: costruzione di un paio di trasformate attraverso la soluzione di un problema di Sturm-Liouville; applica.zione della trasformata e della soluzione dei sistema di equazioni differenziali ordinarie risultanti; applicazione della trasformata inversa. Sono presentati le simulazioni numeriche e i confronti con risultati disponibili in letteratura. In queste simulazioni i coefficienti di dispersione sono stati considerati dipendenti sia dali a distanza dali a sorgente (caso bidimensionale e tridimensionale) sia dai loro valore asintotico (caso tridimensionale). / Neste trabalho são apresentados modelos de dispersão analíticos de qualidade do ar, baseados nas equações de difusão e advecção estacionárias bidimensional e tridimensional, na geometria cartesiana. Estes modelos são válidos para um contaminante passivo, sujeito a situações de turbulência homogênea e com velocidades do vento médio uniformes, ou seja, para abandonos elevados que ocorrem em condições de estabilidade intermediária e na ausência de fortes empuxos. As soluções são alcançadas com o uso do método da Transformada Integral Generalizada e são analíticas no sentido de que nenhuma aproximação é feita ao longo de sua derivação. A idéia básica do método consiste: elaboração de um par de transformadas pela solução do problema de Sturm-Liouville; aplicação da transformada integral e solução do sistema de equações diferenciais ordinárias resultantes e aplicação da transformada inversa. Simulações numéricas e comparações com resultados disponíveis na literatura são apresentadas. Nestas simulações foram considerados coeficientes de dispersão dependentes da distância da fonte (caso bidimensional e tridimensional) e o valor assintótico do coeficiente de dispersão (caso bidimensional). / In this work, analytical dispersion models for air quality are reportedo The models are based upon the multidimensional steady-state diffusion-advection equation, in the cartesian geometryo They are valid for a passive contaminant and for physical conditions of homogeneous turbulence and uniform mean wind speed, that means, for elevated contaminant releases that occur in intermediate stability conditions without strong buoyancyo The solutions are obtained using the Generalized Integral Transform Technique; they are analytical in the sense that no approximation is made along their derivationo The main idea consists of the steps: construction of a pair of transformation through the solution of a Sturm-Liouville problem, application of the integral transform and solution of result system of ordinary differential equations followed by the application of the inverse transformo Numerical simulations and comparisons with available results in the literature are presentedo In these simulations the dispersion coefficients were assumed to be dependent on the distance o f the source (two and three dimensional cases) as well its asymptotic bebaviour (two dimensional case).

Orientador: Edson Del Rio Vieira / Coorientador: Sergio Said Mansur / Banca: João Batista Aparecido / Banca: Francisco Augusto Aparecido Gomes / Resumo: Este trabalho experimental realizado em um túnel hidrodinâmico vertical estudou os padrões do escoamento ao redor de cilindros elípticos com diferentes razões entre eixo. A investigação foi executada devido à importância destes escoamentos na engenharia e também devido à relativa escassez de pesquisas que aborde as características pertinentes a este trabalho. Os ensaios foram realizados com cilindros elípticos com razões de eixos de 0,33; 0,4; 0,5 e 0,67 e com um cilindro circular, para números de Reynolds variando até cerca de 2000. Os cilindros elípticos foram posicionados de forma que o maior eixo da elipse fosse paralelo à direção do escoamento. Os resultados obtidos mostram imagens dos escoamentos visualizados em torno dos cilindros para diferentes regimes de escoamento e a variação da frequência de emissão de vórtices (número de Roshko e número de Strouhal) com o número de Reynolds. Foram comparados os resultados referentes aos escoamentos em torno do cilindro circular com trabalhos precedentes e observou-se boa concordância, validando os demais resultados. A comparação entre as visualizações dos escoamentos para os diferentes cilindros elípticos revelaram que os cilindros com menores razões entre eixos (mais delgados) se aproximam mais da definição de corpos aerodinâmicos. Assim como observado nas visualizações de escoamento, os resultados obtidos com a anemometria de filme quente mostraram que, para números de Reynolds inferiores a 500, os escoamentos em torno dos cilindros ensaiados apresentaram comportamento semelhante em relação à frequência de emissão de vórtices. Quando os valores de Re ultrapassam a faixa de 500, a variação da frequência de emissão referente a cada cilindro difere entre os casos ensaiados. Os resultados referentes à variação do número adimensional... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: This work is an experimental study about the flow around elliptical cylinders with different axis ratios carried out in a hydrodynamic vertical tunnel. The investigation was performed due to the relevance of these flows in an engineering view point and because the lack of research in this area. Tests were performed with elliptical cylinders with axis ratios of 0.33, 0.4, 0.5 and 0.67 and with a circular cylinder, for Reynolds numbers up to about 2,000. The elliptical cylinders were positioned transversally to the non-perturbed flow with null attack angle. Qualitative results were presented as the images of the visualized flow and quantitative results shown the vortex shedding frequency, the Roshko number and the Strouhal number in function of Reynolds number. The results for the flow around the circular cylinder were compared with previous work and showed a good agreement, validating the other results. A comparison of the flow visualizations for the different cylinders revealed that the cylinders with smaller axis ratios (thinner) are closer to the definition of aerodynamic bodies. The results obtained with hot film anemometry showed that for Reynolds numbers less than 500, the flows around the tested cylinders showed very similar behavior to the vortex shedding frequency. For Reynolds numbers more than 500, the variation of the vortex shedding frequency differs significantly in each cylinder tested. The results concerning of the Roshko number in function of Reynolds number exhibited similar behavior to those observed in the variation of frequency. As regards the variation of the Strouhal, for Reynolds numbers up to about 500, there is an irregularity in the behavior of the Strouhal number. However, for Reynolds... (Complete abstract click electronic access below) / Mestre

Devido à grande magnitude com que as correntes de turbidez ocorrem no ambiente natural, a simulação física em escala reduzida tem sido a ferramenta usada pelos pesquisadores para expandir a compreensão do desenvolvimento de tal fenômeno na natureza. O presente trabalho apresenta a caracterização reológica de misturas de água e sedimento que são utilizadas em laboratório na simulação física de correntes de turbidez. O objetivo principal deste trabalho é comparar os dados obtidos via ensaio reológico com as equações reológicas frequentemente usadas para o cálculo da viscosidade. O estudo contemplou o uso de um equipamento específico, o reômetro, para realização da medição das amostras e obtenção de reogramas (relação tensão/taxa de deformação) e o uso de 3 spindles diferentes (cone e prato, vane e cilindro concêntrico). Os ensaios foram realizados com uma mistura de água, caulim e/ou carvão em diferentes proporções de sedimentos nas concentrações volumétricas de 5, 10, 15, 20, 25 e 30%. Os resultados experimentais da caracterização reológica ajustaram-se bem aos modelos de Newton e Herschel-Bulkley. Com os resultados também foi possível determinar a máxima fração de empacotamento por meio do método gráfico de Liu (2000), cujos valores encontrados estão consonantes com os mencionados na literatura. Os resultados demonstraram que se deve ter cautela ao se generalizar a utilização de um modelo para prever a viscosidade relativa de suspensões. Por fim, foram propostas duas equações para estimativa da viscosidade relativa em função da concentração volumétrica de cada sedimento, da máxima fração de empacotamento e da porcentagem de argila presente na mistura. / Due to the large magnitude that the turbidity currents occur on the natural environment, physical simulation on a small scale has been the tool used by researchers to expand the understanding of the development of such a phenomenon in nature. This paper presents the rheological characterization of mixtures of water and sediment that are used in the laboratory in the physical simulation of turbidity currents. The main objective of this work is to compare the data obtained via rheological test with the rheological equations often used to calculate the viscosity. The study comprises the use of a specific tool, the rheometer, for realization of the measurement of the samples and obtaining of rheograms (shear stress/shear rate) and the use of three different spindles (cone and plate, concentric cylinder and vane). Assays were performed with a mixture of water, kaolin and/or coal in different proportions of sediment in the volumetric concentrations of 5, 10, 15, 20, 25 and 30%.With the results it was also possible to determine the maximum packaging fraction through the graphical method of Liu (2000), whose found values are correspondents with those mentioned in the literature. The results demonstrate that we must be careful to generalize the use of a model to predict the relative viscosity of suspensions. Finally, have been proposed two equations to estimate the relative viscosity due of volumetric concentration of each sediment, the maximum packing fraction and the percentage of clay present in the mixture.

Made available in DSpace on 2016-12-23T14:08:12Z (GMT). No. of bitstreams: 1 Estudo de um Novo Modelo Constitutivo de Fluidos.pdf: 1977718 bytes, checksum: 7a35032e9d735176c39ffb1da6e70cef (MD5) Previous issue date: 2008-11-28 / Flow of viscoplastic material is extremely common in many industrial processes. A practical application includes the drilling of petroleum wells where the drilling mud must carry the drill chips with a minimum pumping power. It is obtained using a drilling mud as a highly viscoplastic material. The material used on cementation of drilling wells also exhibit a highly viscoplastic behavior. Other interesting phenomenon involving viscoplastic materials is the mucus displacement in pulmonary airways. The good understanding of the physical mechanism involved in all processes mentioned before is dependent of a good rheological characterization of the viscoplastic material involved. The viscoplastic behavior of the liquid is generally modelled by the generalized Newtonian liquid model with a typical viscosity function equation. One of the most common equation employed is the Herschel-Bulkley viscosity function. This equation predict an infinite viscosity in the limit of zero-shear-rate and this behavior is not compatible with continuity equation for many complex flows. In the present work a viscosity function recently proposed by Souza Mendes and Dutra (2004) is used to perform numerical simulations via finite element method of the flow of yield-stress materials. This viscosity function encompasses most of the viscoplastic and pseudoplastic models such as Herschel-Bulkley, Papanastasiou, Bingham and Carreau as special cases. In order to compare with the results available in the literature, the geometry chosen is the classic 4:1 abrupt contraction. One important dimensionless number analyzed is the jump number J that gives a relative measure of the shear rate jump that occurs when the extra-stress-tensor reaches the yield-stress plateau on the flow field / O escoamento de materiais viscoplásticos é extremamente importante devido sua vasta aplicação industrial. Uma aplicação prática do escoamento de materiais viscopláticos ocorre na perfuração de poços de petróleo. O carreamento do cascalho produzido na perfuração de poços de petróleo é obtido com uso da lama de perfuração. A utilização de lamas de perfuração com propriedades viscoplásticas diminui a potência necessária ao bombeamento desse material. O material utilizado na cimentação dos poços de petróleo também possui propriedades altamente viscoplásticas. Um outro fenômeno que está associado a um escoamento de material tipicamente plástico é o escoamento do muco nas vias áereas pulmonares. O bom entendimento dos mecanismos físicos envolvidos nos processos mencionados depende de uma boa caracterização reológica dos materiais viscoplásticos envolvidos. O comportamento viscoplástico é geralmente modelado pelo modelo de Fluido Newtoniano Generalizado com uso de alguma função típica para descrever a viscosidade. Uma equação para a viscosidade muito utilizada é a de Herschel-Bulkley. Essa equação prevê viscosidade infinita para valores de tensão abaixo da tensão limite de escoamento. Esse comportamento não é compatível com o comportamento previsto pela equação da continuidade para fluidos complexos. No presente trabalho utiliza-se a função viscosidade recentemente proposta por Souza Mendes e Dutra (2004) para uma análise numérica, utilizando-se o método dos elementos finitos, do escoamento de materiais viscoplásticos. Essa função engloba a maioria dos modelos viscoplásticos e pseudoplásticos tais como Herschel-Bulkley, Papanastasiou, Bingham e Carreau. A fim de comparar os resultados previstos pelo presente trabalho com os disponíveis na literatura, a geometria escolhida é a clássica contração abrupta 4:1. Um importante número adimensional analisado é o número J que dá uma medida relativa do salto na taxa de deformação quando atingida a tensão limite de escoamento.

A Província Carajás é notável por sua diversidade metalogenética. Além de conter depósitos de ferro de classe mundial, destaca-se por apresentar a maior quantidade de depósitos de óxido de ferro-cobre-ouro (IOCG) conhecida mundialmente. Na porção centro-sul da Província Carajás ao longo da Zona de Cisalhamento de Carajás encontram-se importantes depósitos IOGC como Sossego (corpos Sequeirinho e Pista), Bacaba, Castanha, Bacuri e Cristalino. O depósito de Borrachudos localiza-se a aproximadamente 6 km a sudoeste do depósito Cristalino (500 Mt @ 1,0% Cu e 0,3 ppm Au), próximo à Serra do Rabo, no limite entre a Bacia Carajás e o seu embasamento mesoarqueano. As rochas hospedeiras do Borrachudo compreendem metavulcânicas ácidas (riolitos) e intermediárias (andesitos) atribuídas à Formação Parauapebas do Grupo Grão-Pará, de idade neoarqueana e a diques de microgabro. O sistema hidrotermal desenvolveu-se em dois eventos. O primeiro relaciona-se a estruturas dúcteis, com desenvolvimento de estágios de alterações hidrotermais semelhantes à IOCG, alteração sódica, sódica-calcica, ferro-(P) e potássica, com o qual vincula-se fraca e subordinada mineralização de cobre. No segundo evento, há recorrência de alguns estágios hidrotermais como silicificação I, alteração sódica II e alteração potássica II, além de cloritização e carbonatização, todos vinculados a estruturas rúpteis, bem como a mineralização principal de cobre, representada por calcopirita, magnetita, pirita, siegenita, cobaltita, galena, anglesita, uraninita e thorita. Os fluidos hidrotermais aquosos e carbônicos do primeiro evento, no qual se formou quartzo I pré-mineralização, possuem temperatura (300 a > 600º C) e salinidade mais elevados (40 a 60 % equiv. NaCl), e possivelmente foram originados a partir da imiscibilidade de um fluido magmático. Os fluidos hidrotermais aquosos, mais tardios, relacionados ao segundo evento, que formou quartzo II sin-mineralização, possuem menores valores de temperatura (150 a 350° C) e salinidade (30 a 40 % equival. NaCl) e indicam um processo mais avançado de misturas de fluidos externos diluídos, que teve efetiva importância para a precipitação do cobre na forma de sulfetos. O primeiro evento hidrotermal não possui ainda idade conhecida, entretanto devido à estruturação dúctil e a correlação com outros depósitos da Província Carajás, possivelmente são de idade neoarqueana, o que deve ser confirmada em trabalhos futuros. O segundo evento hidrotermal possui idade de 2.011 ± 6,8 Ma (MSWD = 3,9), obtida em titanita hidrotermal associada com a mineralização principal de cobre, e está relacionada a eventos tectônicos riacianos. Esse estudo aponta para a recorrência de múltiplos eventos hidrotermais responsáveis pela mealogênese de cobre, intrinsicamente controlados pela evolução magmática e tectônica da Província Carajás. / The Carajás Province is notable for its metallogenic diversity. It contains world class iron deposits and stands out for presenting the highest amount of iron-oxide copper-gold deposits (IOCG) known worldwide. In the south-central portion of the Carajás Province, within the Carajás Shear Zone, occur important IOGC deposits, such as Sossego, Bacaba, Castanha, Bacuri, and Cristalino. The Borrachudo copper deposit is located about 6 km southwest of the Cristalino deposit (500 Mt @ 1.0% Cu and 0.3 ppm Au), near Serra do Rabo region, at the boundary between the Carajás Basin and its Mesoarchean basement. Host rocks at Borrachudo comprise acid (rhyolites) and intermediate (andesite) metavolcanic rocks assigned to the Neoarchean Parauapebas Formation, Grão Pará Group, Itacaiúnas Supergroup and possibly coeval microgabbro dikes. The hydrothermal system developed in two events. The first was related to ductile structures, with development of hydrothermal alteration types similar to those of IOCG deposits, including sodic, sodic-calcic, iron-phospate, and potassium alteration. Weak and subordinate copper mineralization was associated with this early hydrothermal event. In the latter event, there was recurrence of hydrothermal alteration stages, sucs as silicification, sodic alteration e potassic alteration, besides chlorite and carbonate alteration. These alterations were linked to the brittle structures. The main copper mineralization represented by chalcopyrite, magnetite, pyrite, siegenite, cobaltite, galena, anglesite, uraninite and thorite was developed during this late event. The aqueous and carbonic hydrothermal fluid associated with the first hydrothermal event, in which quartz I pre-mineralization was formed, have higher temperature (300 to> 600 °) higher salinity (40 to 60% equivalence. NaCl). These possibly were originated from the immiscibility of a magmatic fluid exsolved during granite crystallization. The aqueous hydrothermal fluids, related to the second event which formed sin-ore quartz II, have lower values of temperature (150 to 350 °C) and salinity (30 to 40% equivalence. NaCl), indicating a more advanced mixing process involving diluted external fluids. Those fluids had an effective importance for the copper precipitation. The timing of the first hydrothermal event is yet unkown, however its correlation with the reactivation of ductile shear zones suggest a Neoarchean age, similar to that recorded at Cristalino and Sequeirinho (Sossego mine). The second hydrothermal event has an U-Pb age of 2011 ± 6.8 Ma (MSWD = 3.9) obtained in hydrothermal titanite associated with the main copper mineralization of copper, is related to Orosirian tectonic events. This study points to the recurrence of multiple hydrothermal events responsible for the copper metallogenesis, which were intrinsically controlled by the magmatic and tectonic evolution of the Carajás Province.

O presente trabalho apresenta o estudo e implementação de um algoritmo numérico para análise de escoamentos turbulentos, tridimensionais, transientes, incompressíveis e isotérmicos, através da Simulação de Grande Escalas, empregando o Método de Elementos Finitos. A modelagem matemática do problema baseia-se nas equações de conservação de massa e quantidade de movimento de um fluido quase-incompressível. Adota-se um esquema de Taylor-Galerkin, com integração reduzida e fórmulas analíticas das funções de interpolação, para o elemento hexaédrico de oito nós, com funções lineares para as componentes de velocidade e constante no elemento para a pressão. Para abordar o problema da turbulência, emprega-se a Simulação de Grandes Escalas, com modelo para escalas inferiores à resolução da malha. Foram implementados o modelo clássico de Smagorinsky e o modelo dinâmico de viscosidade turbulenta, inicialmente proposto por Germano et al, 1991. Uma nova metodologia, denominada filtragem por elementos finitos independentes, é proposta e empregada, para o processo de segunda filtragem do modelo dinâmico. O esquema, que utiliza elementos finitos independentes envolvendo cada nó da malha original, apresentou bons resultados com um baixo custo computacional adicional. São apresentados resultados para problemas clássicos, que demonstram a validade do sistema desenvolvido. A aplicabilidade do esquema utilizado, para análise de escoamentos caracterizados por elevados números de Reynolds, é discutida no capítulo final. São apresentadas sugestões para aprimorar o esquema, visando superar as dificuldades encontradas com respeito ao tempo total de processamento, para análise de escoamentos tridimensionais, turbulentos e transientes .