O presente trabalho apresenta um estudo numérico e experimental sobre o desempenho aerodinâmico de turbinas eólicas de eixo vertical envolvendo rotores Savonius convencional de duplo-estágio e helicoidal. O estudo experimental é realizado no Túnel Aerodinâmico Professor Debi Pada Sadhu, do Laboratório de Mecânica dos Fluidos da UFRGS. As simulações numéricas são realizadas com o software Fluent/ANSYS utilizando o Método dos Volumes Finitos. São comparados os coeficientes de torque estático e dinâmico, o coeficiente de potência, além de uma análise aerodinâmica das duas turbinas. As medições são realizadas empregando Tubos de Pitot, um torquímetro estático digital e um torquímetro simples construído para a medição do torque dinâmico. As turbinas são fabricadas através da técnica de prototipagem 3D, com uma semelhança de dimensões e parâmetros. As soluções numéricas são resolvidas através da equação da continuidade, das equações de Navier-Stokes com médias de Reynolds (RANS) e pelo modelo de turbulência k-ω SST. A qualidade da malha utilizada é avaliada através do método de Índice de Convergência de Malha (GCI), para três diferentes tamanhos de malha. São feitas análises dos rotores na forma estática para diferentes ângulos de incidência e com a turbina em rotação são feitas análises para diferentes razões de velocidades de ponta de pá (λ). Resultados demonstram que a turbina helicoidal apresenta um coeficiente de torque positivo para todos os ângulos do rotor, assim como a turbina convencional de dois estágios. O coeficiente de torque dinâmico da turbina helicoidal é superior ao da turbina de duplo-estágio para a maioria dos casos, e também apresenta menor oscilação de torque ao longo de cada rotação. Por consequência, o coeficiente de potência do rotor helicoidal também se tornou superior, com um valor máximo encontrado na ordem de 11,8% para um λ de 0,65 no caso experimental, e de 8,4% para o mesmo λ no caso numérico, quando comparado com o rotor de duplo-estágio. Os erros relativos entre as simulações numéricas e os resultados experimentais estão entre 2,16% e 13,4%. Uma estimativa de potência gerada é feita para ambos os casos, para uma razão de velocidade de ponta de 0,65, onde a turbina helicoidal apresenta melhores resultados em relação ao rotor de duplo-estágio, na ordem de 13,6% para uma velocidade de 10,4 m/s. / This paper presents a numerical and experimental study of vertical axis wind turbine performance comparison involving two-stage and helical Savonius rotors. The experimental study is conducted in the Aerodynamic Tunnel Professor Debi Pada Sadhu at the Fluid Mechanics Laboratory of the UFRGS. The numerical simulations are performed with the Fluent/ANSYS software using the Finite Volumes Method. The static and dynamic torque coefficients, the power coefficients, and an aerodynamic analysis of the two turbines are compared. Measurements are made using Pitot tubes, a digital static torque wrench and a simple wrench constructed for the dynamic torque measurement. The aerodynamics rotors are manufactured by 3D prototyping technique with similar dimensions and parameters. Numerical solutions are solved by the continuity equation, the Reynolds Averaged Navier-Stokes (RANS) equations and the turbulence model k-ω SST. The quality of the mesh used is evaluated used the Grid Convergence Index (GCI) method, for three different mesh sizes. The rotors analyzes are made in static form for different angles of incidence and for the rotating turbine analyzes are made for differents tip speed ratio (λ). Results show that the helical turbine has a positive static torque coefficient for any rotor angles, as well as conventional two-stage turbine. The dynamic torque coefficient of the helical turbine is higher than the two-stage turbine for most cases and also shows less torque variation along each rotation. Consequently, the power coefficient of the helical rotor also become higher, with a maximum value found on the order of 11.8% for a λ of 0.65 in the experimental case, and 8.4% for the same λ number when compared with the two-stage rotor. The relative errors between the numerical simulations and the experimental results are between 2.16% and 13.4%. A generated power estimate is made for both cases, for a tip speed ratio of 0.65, where the helical turbine provides better results compared to two-stage rotor in order of 13.6% for a velocity of 10.4 m/s.

Os métodos de controle baseados em dados têm como diferencial utilizar os dados coletados de experimentos para realizar a sintonia do controlador de uma forma direta e sem utilizar nenhum modelo do processo. Um dos métodos já consagrados é o conhecido Iterative Feedback Tuning, método este caracterizado como um algoritmo iterativo de otimização que busca a minimização do critério de desempenho norma H2. O presente trabalho propõe uma extensão baseada no clássico algoritmo Iterative Feedback Tuning para sistemas em formato cascata com duas malhas de controle. O formato cascata com duas malhas de controle é uma estrutura comum e muito utilizada em aplicações de cunho industrial. A sintonia deste tipo de estrutura em metodologias clássicas de controle envolve geralmente duas etapas. Primeiramente a malha interna é sintonizada para então posteriormente a malha externa. A vantagem do método de controle baseado em dados proposto está em realizar a sintonia de ambas as malhas de forma simultânea e concomitantemente. É escolhida como aplicação um veículo aéreo não tripulado com quatro rotores, mais conhecido como quadricóptero. Esta aplicação tem uma dinâmica complexa, uma vez que o sistema é classificado como sendo não-linear e multivariável, algo demonstrado através do equacionamento dinâmico do quadricóptero que é apresentado. Uma série de simulações e experimentos práticos são apresentadas com o intuído de demonstrar a dinâmica descrita no modelo e também demonstrar a aplicação do algoritmo Cascade Iterative Feedback Tuning para os mesmos. / The methods of data-driven control have the distinction to use the data collected from experiments to perform the tuning of the controller directly and without using any process model. One method already established is the Iterative Feedback Tuning, this method is characterized as an iterative optimization algorithm that aims to minimize the H2 norm performance criterion. The present work proposes an extension based on the classical Iterative Feedback Tuning algorithm for cascade systems with two control loops. The cascade format with two loops is a common structure and widely used in industrial applications. In classical control methods, the tuning of this type of structure usually involves two steps. First, the inner loop is tuned and subsequently the outer loop. The advantage of the proposed data-driven control method is to perform this tuning of both controller simultaneously and concurrently. We choose as an application a unmanned aerial vehicle (UAV) with four rotors, better known as quadcopter. This application has a complex dynamic, since the system is classified as being non-linear and multivariable, demonstrated by the quadcopter dynamic equations. A series of simulations and practical experiments are presented to demonstrate the dynamicmodel described and also to validate the method CIFT algorithms.

A simulação numérica de componentes estmturais necessita de dados do material que são obtidos por ensaios mecânicos. Materiais não-lineares, como termoplásticos, podem apresentar em ensaios de tração estricção quando submetidos à deformações firútas. Este fenômeno normalmente ainda está associado a campos de deformações heterogêneos que possuem uma cinemática particular, se propagando ao longo do corpo. A formação e propagação da estricção podem mascarar o real comportamento mecârúco do material quando obtido por um ensaio de tração, levando o processo de caracterização da curva tensão-deformação real do material a um procedimento não-trivial. Através de um ensaio de tração de um termoplástico, dados experimentais de força e do campo de deslocamentos da região de estricção, obtido através de métodos ópticos, são utilizados em um procedimento numérico-experimental de otimização para a deternúnação dos parâmetros (FEMU) de alguns modelos constitutivos capazes de representar a tendência do comportamento de termoplásticos. De modo a estudar a capacidade representativa da resposta mecânica do ensaio, um modelo constitutivo multi linear e um modelo constitutivo variacional, são estudados. Uma função objetivo conveniente que utiliza dados experimentais e numéricos é usada para considerar a cinemática da estricção, responsável pela mudança geométrica que diferencia as respostas nonúnais e reais da curva tensão-deformação. Por meio dos resultados obtidos neste trabalho verificou-se que a caracterização dos modelos constitutivos utilizando somente a resposta de força, não garante uma caracterização constitutiva adequada, sendo necessária a resposta de deslocamento para garantir a representação da cinemática da região de estricção. / Numerical simulation of structural components requires material data obtained from mechanical testing. Nonlinear materials, like thermoplastics, submitted to tensile testing presents necking undergoing finite strain. This phenomenon is still typically associated to heterogeneous strain fields which possess a particular kinematic, propagating through the body. Necking and colddrawing may mask the actual mechanical behavior measured by a tensile testing machine, turning the real stress-strain curve characterization process into a non-trivial procedure. Through a thermoplastic tensile testing, force experimental data and displacement experimental data from the necking region, obtained by optical methods, are used in a numerical-experimental optimization procedure in order to determine the constitutive parameters (FEMU) from some constitutive models which are able to represent the thermoplastic behavior tendency. In order to study the representative capacity of the mechanical testing response, a multilinear constitutive model and a variational constitutive model, are studied. A suitable objective function which uses experimental and numerical data is used to consider the necking kinematics, responsible for geometric change that distinguishes the nominal and real responses of the stress-strain curve. Through the obtained results, it can be seen that the model parameters determination using, in addition to force, the displacement field in the objective function is necessary to represent the kinematic behavior in the necking region.

A radiação térmica é um mecanismo de transferência de calor muito importante em processos que envolvem gases participantes, como CO2 e H2O, em temperaturas elevadas. A dependência altamente irregular do coeficiente de absorção em relação ao número de onda torna a integração linha-por-linha (LBL) da equação da transferência radiativa (RTE) impraticável, sobretudo em situações onde a radiação pode ser apenas parte de um problema mais complexo. Modelos espectrais globais, como a soma-ponderada-de-gases- cinza (WSGG) e a soma-ponderada-de-gases-cinza baseada em linhas espectrais (SLW), representam alternativas à integração LBL. Entretanto, alguns modelos requerem o uso da aproximação de escala (que assume que as dependências espacial e espectral da seção transversal de absorção são separáveis). Essa aproximação pode produzir erros consideráveis, principalmente quando existem gradientes elevados de temperatura e de concentração. Neste trabalho, os resultados de alguns modelos espectrais globais foram comparados com os da solução LBL. E, a partir de algumas proposições dos modelos WSGG e SLW (integração espectral através de gases cinza e intervalos espectrais fixos por meio da definição de uma temperatura de referência), foi desenvolvida a soma-ponderada-de-gases-cinza com coeficientes não constantes (NCC-WSGG). Nesse modelo, os coeficientes de absorção e de emissão para cada gás cinza são funções polinomiais da temperatura. O modelo NCC-WSGG foi aplicado em problemas unidimensionais e bidimensionais, envolvendo CO2, H2O e mistura dessas espécies químicas. Comparações com a solução LBL mostraram que o NCC-WSGG pode fornecer resultados muito satisfatórios para o fluxo de calor e para o termo fonte radiativos em problemas não isotérmicos e não homogêneos. / The thermal radiation is a very important mechanism of heat transfer in processes that embody participating gases, as CO2 and H2O, at high temperature. The highly irregular dependence of the absorption coefficient with respect to the wavenumber makes the application of line-by-line (LBL) integration of the radiative transfer equation (RTE) prohibitive, principally in situations where the radiation can be only part of a more complex problem. Global spectral models, like the weighted-sum-of-gray-gases (WSGG) and the spectral line based weighted-sum-of-gray-gases (SLW), are alternatives to the LBL integration. However, some models demand the application of the scaling approximation (which assumes that the spatial and spectral dependences of the absorption cross-section are separable). This approximation can lead to pronounced errors, mainly under high gradients of temperature and concentration. In this work, some results obtained from global models were compared with the LBL solution. And, applying some proposals from WSGG and SLW models (spectral integration through gray gases and fixed spectral intervals by the definition of a reference condition), it was developed the nonconstant coefficient weighted-sum-of-gray-gases (NCC-WSGG). In this model, the absorption and emission coefficients of each gray gas are polynomials functions of the temperature. The model was applied to solve one and two dimensional problems, which were comprise of CO2, H2O and mixtures of these chemical species. Comparisons with the LBL solution showed that the NCC-WSGG can provide very good results for the heat flux and for the radiative heat source under nonisothermal and nonhomogeneous conditions.

Análises numérico-experimentais são fundamentais na simulação numérica do comportamento de alguns tecidos biológicos ou mesmo de sua interação com outras estruturas, como o que acontece entre artérias e stents. Simulações numéricas podem auxiliar no projeto da estrutura tubular ou na escolha do modelo de stent mais adequado para uma dada indicação clínica, bem como no aprimoramento da técnica cirúrgica de posicionamento. A obtenção de dados experimentais do comportamento mecânico de artérias é um desafio devido a fatores biológicos e experimentais, como a retirada da amostra, controle de temperatura e de umidade no armazenamento e manuseio dos corpos de prova. Outras dificuldades se apresentam na fixação da amostra e na realização de medições de deformação, além da mudança que pode ocorrer no tecido para diferentes indivíduos ou mesmo no mesmo indivíduo. Do ponto de vista numérico, é necessário um modelo de material capaz de reproduzir o comportamento experimental observado. O objetivo do presente trabalho é estudar problemas envolvidos na caracterização numérico-experimental de paredes arteriais para então propor uma metodologia apropriada para ser aplicada em uma amostra de aorta de um suíno da raça Landrace. Modelos numéricos adequados para simular o comportamento do tecido em sua região fisiológica de deformação foram estudados, e através dos dados obtidos de procedimentos experimentais foi realizada a identificação de parâmetros do modelo escolhido. / Numerical and experimental analyses are essential to understand the numerical simulation of the behavior of some biological tissues or even their interaction with other structures, such as what happens between arteries and stents. Numerical simulations can help in the tubular structure project or in the choice of the most suitable stent for a given clinical indication, as well as in the improvement of the surgical techniques in the positioning. The experimental data acquisition of arteries’ mechanical behavior is a challenge due to biological and experimental factors such as the removal of the sample, temperature and humidity control in the storage and handling of specimens. Other difficulties are the sample fixation and the strain measurements as well as the change that can occur in the tissue for different individuals or even in the same one. From the numerical point of view is necessary a material model capable to reproduce the observed experimental behavior. The objective of this work is to study the issues involved in the numerical and experimental characterization of arterial walls and then propose a suitable methodology to be applied in a pig's aorta sample of Landrace's breed. Numerical models suitable for simulating the behavior of physiological tissue deformation were studied and, through the data from experimental procedures, parameters were identified to represent it.

A microestrutura resultante do processo de forjamento a quente é influenciada por variáveis como temperatura, deformação, tamanho de grão austenítico inicial e taxa de resfriamento. A proposta deste trabalho é estudar o comportamento termomecânico do aço ABNT 4140 mantido nas temperaturas de 900, 1000, 1100 e 1200ºC por 22, 66 e 200 minutos e, posteriormente, forjado e tratado termicamente em patamares semelhantes aos de recozimento direto e revenimento e, ainda, analisar a influência da temperatura sobre o fenômeno de crescimento de grão. Esta pesquisa correspondeu a uma das etapas do projeto BRAGECRIM (Iniciativa Brasil-Alemanha para Pesquisa Colaborativa em Tecnologia de Manufatura), realizado em parceria entre o Laboratório de Transformação Mecânica (LdTM) da Universidade Federal do Rio Grande do Sul e o Instituto de Conformação Mecânica (IBF) da Universidade Técnica de Aachen (RWTH - Aachen). O projeto objetiva a realização de um estudo para o desenvolvimento do processo de forjamento incremental de eixos vazados para aerogeradores, sendo que este trabalho analisa, de forma independente, as propriedades do material utilizado no forjamento referido. Foram elaboradas sequências e rotas de aquecimento-forjamento-resfriamento, sendo observado o fenômeno metalúrgico de recristalização completa na amostra aquecida a 1200ºC por 200min. / The microstructure resulting from the forging process is studied as a function of variables such as temperature, deformation, starting austenite grain size and cooling rate. The purpose of this work is to study the thermomechanical behavior of steel ABNT 4140 maintained at temperatures of 900, 1000, 1100 and 1200°C for 22, 66 and 200 minutes and subsequently forged and post heat treated at similar levels to those of direct annealing and tempering. Furthermore, the influence of temperature on the phenomenon of grain growth is also analyzed. This work is a part of a research project named BRAGECRIM (Brazil-German Initiative for Collaborative Research in Manufacturing Technology), in collaboration with the Laboratory of Mechanical Manufacturing (LDTM) of the Federal University of Rio Grande do Sul and the Institute of Mechanical Forming (IBF) of the Technical University of Aachen. The aim of the project is to study the development of incremental forging process of hollow shafts for wind turbines. However, the work presented in this study analyzes independently, the properties of the material in relation to the forging process. Sequences and routes of heating-forging-cooling were elaborated, being observed metallurgical phenomenon of complete recrystallization in the specimen heated at 1200°C for 200min.

A intransitividade é uma propriedade dos grafos conectados e orientados que representam interações entre espécies. Essa propriedade está associada `a persistência da coexistência mesmo na presença de competição, sendo o jogo Pedra-Papel-Tesoura o exemplo padrão. Neste trabalho, consideramos uma generalização com quatro espécies, que é o número mínimo necessário para que haja outras interações além do ciclo único (um predador, uma presa). Além disso, introduzimos sítios vazios dinâmicos ao sistema, os quais afetam o nível de intransitividade nas interações. Verificamos então que, ao contrário do que o campo médio prevê, na rede quadrada o modelo apresenta duas transições, cujas localizações dependem das probabilidades de predação e reprodução. A primeira transição ocorre entre um estado absorvente, com apenas uma espécie, para um estado em que todas as espécies coexistem. A transição seguinte ocorre quando uma das quatro estratégias é extinta. Essa dependência com as probabilidades de predação e reprodução demonstra que a estrutura do grafo de interação sozinha não é suficiente para prever os resultados finais em modelos desse tipo. Adicionalmente, probabilidades diferentes de predação permitem ajustar o nível de transitividade do sistema, indicando que é necessária uma quantidade mínima de intransitividade para que a coexistência entre todas as espécies persista. / Intransitivity is a property of connected, oriented graphs representing species interactions that may drive their coexistence even in the presence of competition, the standard example being the three species Rock-Paper-Scissors game. We consider here a generalization with four species, the minimum number of species allowing other interaction beyond the single loop (one predator, one prey). We also consider dynamic empty sites, which affect the amount of intransitivity in the system. We show that contrary to the mean field prediction, in the square lattice the model presents two transitions as the parameters setting the predation and reproduction rates change, from an absorbing state with one species only, to a coexistent state first, and then to a state in which one species gets extinct. Such a dependence on the predation and reproduction rates shows that the interaction graph structure alone is not enough to predict the outcome of such models. In addition, different predation rates permit to tune the level of transitiveness, indicating that for the coexistence of all species to persist, there must be a minimum amount of intransitivity.

Em túneis, a deformação do maciço e a pressão do solo no revestimento dependem das tensões e características do maciço, das tensões do maciço, da geometria, da rigidez, do instante da colocação do revestimento. A variação das tensões no revestimento e as deformações são causadas pelo avanço da escavação e as propriedades dependentes do tempo do maciço rochoso e revestimento. O objetivo desta dissertação é apresentar um modelo numérico para a análise tridimensional de túneis em maciços rochosos no programa ANSYS. Ele constitui uma ferramenta para calcular um túnel, revestido de concreto em diferentes situações de rigidez do revestimento, processo construtivo e características do maciço. O processo de escavação e colocação do revestimento é simulado pelo método da ativação/desativação de elementos em etapas sequenciais. O maciço apresenta comportamento elástico na primeira parte, e posteriormente uma lei viscoplástica. O concreto projetado e o concreto pré-moldado do revestimento são modelados como elásticos na primeira parte, e posteriormente, como materiais viscoelásticos utilizando o modelo da cadeia de Maxwell e de Kelvin, propostos por Bazant. Ambos os modelos foram calibrados por comparação com o Código Modelo fib 2010. A simulação do túnel é realizada através da incorporação de sub-rotinas para o concreto viscoelástico no programa ANSYS. Para validar o modelo são feitas comparações dos resultados obtidos no programa ANSYS com a solução analítica, com o programa GEOMEC91 e com o Novo Método Implícito, dependendo do caso. Um estudo paramétrico em elasticidade é realizado com o código ANSYS variando alguns parâmetros relevantes. Na segunda parte, o mesmo túnel revestido é modelado com dois diferentes tipos de revestimento e os dois modelos de cadeia citados. Finalmente, é comparado o comportamento com o túnel experimental Kielder, do qual se tem dados disponíveis da instrumentação. A comparação destes dados com os resultados do programa ANSYS mostra uma aproximação muito boa. / Stabilizing underground openings such as tunnels excavated in rock mass remain a major concern of geotechnical engineers dealing with this kind of structures. In tunnels, the rock mass strain and the ground pressure on lining depend on the stress and characteristics of the rock mass as well as of the geometry, the stiffness and the moment of the lining installation. Pressure variation on lining and strain are caused by the advance of excavation and the time-dependent properties of the rock mass and lining. This dissertation describes the numerical implementation of constitutive laws in the numerical code ANSYS. It is performed a numerical simulation with 3D finite elements of a tunnel. The rock mass presents elastic behavior in the first part and then a viscoplastic law. The shotcrete and precast concrete are modeled as elastic in the first part and then viscoelastic material with the constitutive laws of the Maxwell and the Kelvin chain model, proposed by Bazant. So the characteristics of the viscoelastic concrete changes with time until the stabilization. Both chain models had to be calibrated by comparing with de fib Model Code 2010. The finite elements simulation is performed by incorporating subroutines for the viscoelastic concrete model in the ANSYS code. The method to simulate the tunnel excavations is by activating and deactivating elements in sequential steps. In the first part of the dissertation two validations are performed. The analytical solution and the deformation achieved on the stabilization in the ANSYS code are compared for the validation of an unlined tunnel. For a lined tunnel, with an elastic and viscoplastic rock mass and elastic lining, validation is performed by comparing the results of the GEOMEC91 code with the ANSYS code. The results show a very good approximation in viscoplasticity and agree perfectly in elasticity. Also a parametric study is undertaken with the ANSYS code varying some relevant parameters. In the second part, it is compared the same tunnel with the two different concrete lining and the two chain models. Finally, it is compared with the Kielder experimental tunnel, which in situ measured data is available. The comparison of this data with the results from the ANSYS code shows a very good approximation.

O interesse no desenvolvimento de sistemas eficientes de combustão para reduzir a poluição ambiental e aumentar a eficiência de queima tem chamado a atenção para a combustão em meios porosos inertes. A recirculação de calor, induzida pela matriz sólida a partir dos produtos quentes para os reagentes frios, aumenta a temperatura da chama melhorando a sua estabilidade e permitindo a utilização de combustíveis com baixo poder calorífico. Um estudo teórico recente mostra que uma chama estabilizada por um plano de estagnação imersa em um meio poroso pode, sob certas condições, estender os limites de inflamabilidade de uma mistura de ar e combustível. Por outro lado, o plano de estagnação é um problema que simula o efeito da taxa de deformação do escoamento sobre a estabilidade da chama, o que é relevante para várias configurações de queimador poroso. Portanto, o foco deste trabalho é o estudo da combustão em um queimador poroso com um plano de estagnação. Um experimento é conduzido com empacotamento de esferas, onde uma chama pode ser estabilizada por plano de estagnação devido a um anteparo. A razão de equivalência e a taxa de deformação são controladas pelos fluxos de ar e de combustível e da distância entre injetor e anteparo. A posição da chama é aproximadamente determinada pelo campo de temperaturas medidas por termopares. Complementarmente é realizada uma análise numérica simplificada do problema na qual se pode verificar o efeito da taxa de deformação sobre a estabilidade de chamas em queimadores porosos. Os resultados mostram que é possível estabilizar chamas no interior do meio poroso com plano de estagnação, porém, não foi possível atribuir um aumento de temperatura devido ao aumento da taxa de deformação. / The interest in developing efficient combustion systems to reduce environmental pollution and increase the burning efficiency has called attention to the combustion in inert porous media. The heat recirculation, induced by the solid matrix, from the hot products to the incoming cold reactants, increases the flame temperature and improves its stability, allowing for the use of fuels with low heat content. A recent study shows theoretically that a flame stabilized by a stagnation plane immersed in a porous medium may, under certain conditions, to extend the flammability limits of a mixture of fuel and air. On the other hand, the stagnation plane imposes a certain strain rate on the flow field, which is relevant to various porous burner configurations. Therefore, the focus of this work is the study of combustion in a porous burner with a stagnation plane. An experiment is conducted with packing bed of spheres where a flame can be stabilized against a stagnation plane. The equivalence ratio and the strain rate are controlled by the flows of air and fuel and the distance between the injector and the stagnation plane. The flame position is approximately determined by the temperature field measured by thermocouples. In addition, it is performed a simplified numerical analysis of the problem in which one can see the effect of the strain rate on the stability of flames in porous burners. The results show that it is possible to stabilize flames within the porous medium with stagnation plane, however, it has not been possible to assign a temperature increase due to the increased strain rate.

O estudo do comportamento reológico de fluidos não Newtonianos tem grande importância em diversas áreas da engenharia. O aumento na demanda destes fluidos por exemplo, no uso doméstico, pessoal e em processos químicos acarreta em dificuldades que vão desde o processo de sua mistura ao seu manuseio. Dentre os fluidos não Newtonianos estão os viscoelásticos, os quais exibem deformação aparente quando os níveis de tensões são inferiores ao limite de escoamento do material e, dentro desta classe, alguns ainda podem apresentam comportamento elástico quando submetidos à baixa taxa de cisalhamento. Juntamente com os efeitos elasto-viscoplásticos, os materiais podem apresentar comportamento tixotrópico, onde devido as tensões sua reestruturação não é instantânea. Na presente Tese, fez-se um estudo numérico a fim de analisar o problema especifico de escoamentos de fluidos elasto-viscoplásticos com comportamento tixotrópico em uma expansão planar abrupta na razão de aspecto 4 : 1. O modelo mecânico aplicado consiste de uma equação viscoelástica para o campo de tensões, uma evolutiva para o parâmetro de estrutura do material, bem como as equações de conservação de massa e momento. O modelo mecânico aplicado mostrou-se capaz de prever o comportamento tixotrópico. A aproximação numérica do modelo aplicado foi feita através do método estabilizado de elementos finitos, especificamente o método Galerkin Mínimos-Quadrados (GLS), o qual foi implementado no código de elementos finitos para fluidos não Newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Os fenômenos reológicos presentes no problema foram analisados a partir da influência da cinemática do escoamento, da elasticidade e da tixotropia, no nível de estruturação do material, na posição e tamanho das regiões não escoadas e na deformação elástica do material. Os resultados mostraram-se satisfatórios, pois condizem com os apresentados na literatura, apontando a boa predição do modelo mecânico aplicado bem como a robustez de sua implementação computacional. / The rheological behavior of non-Newtonian uids study is of great then importance in many areas of engineering. The increase in demand of these uids - for example, domestic use, personal and processes chemical - causes di culties ranging from the process of mixing it to handling. Among the non-Newtonian uids are viscoelastic, which exhibit apparent deformation when stress levels are lower than the yield limit of the material and, within this class, some still have elastic behavior when subjected to low shear rates. Together with the elastic-viscoplastic e ects, materials may exhibit thixotropic behavior, ie, due the restructuring strain is not instantaneous. In this thesis was made a numerical study to simulate the speci c problem of ow of elastic-viscoplastic uids with thixotropic behavior in abrupt planar expansion { common geometry in industrial systems associated with elastic-viscoplastic uids, whose ratio the aspect is 4:1. The mechanical behavior of most of these structured materials, are highly non-Newtonian, with this, there is need to make them more pseudoplastic, causing undesirable behaviors such as thixotropic which is a phenomenon somewhat characterized and modeled in the literature. The mechanical model applied is able to predict thixotropic behavior and is composed of a viscoelastic equation for the stress eld and an evolutionary to the material structure parameter in addition to the mass and momentum conservation equations. This mechanical model is approximated by a stabilized nite element model, called the Galerkin method of least squares (GLS). In order to study the rheological phenomena present in the problem is analyzed the in uence of the ow kinematics, elasticity and thixotropy in the level of structure of the material, in position and size of unyielded regions and the elastic deformation of the material. The results were satisfactory, since the the study of intensity the U the results agree with those reported in the literature pointing to good prediction of the mechanical model applied well as the hardiness of their computational implementation.The results of the elasticity showed quite spectacular behavior of unyielded regions of the material since for high values of relaxation time, the unyielded region of channel larger has the form of " ngers"because of the high exibility of the material along the line of symmetry of the channel. The results of the analysis of the elastic deformation show, that model correctly dosing the elasticity in the unyielded regions. The thixotropic e ects reported a slower structuring of the material with increasing characteristic time, in response to strain change caused by the expansion, results indicate that for high values of relaxation time, higher the distance between the unyielded regions of smaller channel and larger.