Estudo experimental do escoamento bifásico líquido gás em golfadas com leve mudança de direção

Gas-liquid two-phase flow is very common in industrial applications, especially in the oil and gas, chemical, and nuclear industries. As operating conditions change such as the flow rates of the phases, the pipe diameter and physical properties of the fluids, different configurations called flow patterns take place. In the case of oil production, the most frequent pattern found is slug flow, in which continuous liquid plugs (liquid slugs) and gas-dominated regions (elongated bubbles) alternate. Offshore scenarios where the pipe lies onto the seabed with slight changes of direction are extremely common. With those scenarios and issues in mind, this work presents an experimental study of two-phase gas-liquid slug flows in a duct with a slight change of direction, represented by a horizontal section followed by a downward sloping pipe stretch. The experiments were carried out at NUEM (Núcleo de Escoamentos Multifásicos UTFPR). The flow initiated and developed under controlled conditions and their characteristic parameters were measured with resistive sensors installed at four pipe sections. Two high-speed cameras were also used. With the measured results, it was evaluated the influence of a slight direction change on the slug flow structures and on the transition between slug flow and stratified flow in the downward section.

Estudo experimental e numérico de juntas adesivas em T com aderentes de alumínio

Nos dias de hoje, a ligação adesiva de estruturas complexas que não poderiam ou não seriam tão fáceis de ser fabricadas numa só peça é cada vez mais usual. As juntas adesivas têm vindo a substituir muitos outros métodos de ligação, como por exemplo ligações aparafusadas, rebitas ou soldadas, devido às vantagens de facilidade na sua fabricação, resistência superior e capacidade de unir materiais diferentes. Por esta razão as juntas adesivas têm vindo a ser aplicadas cada vez mais em várias industrias como aeroespacial, aeronáutica, automóvel, naval e calçado. O tipo de adesivo a usar em determinada aplicação é principalmente escolhido consoante as suas características mecânicas e o tipo de resposta pretendida às solicitações impostas. Como exemplo de adesivo resistente e frágil existe o Araldite® AV138. Por outro lado, o adesivo Araldite® 2015 é menos resistente, mas apresenta maior ductilidade e flexibilidade. Além dos adesivos Araldite® comerciais, existem adesivos de poliuretano que combinam características de elevada resistência com características de grande ductilidade e flexibilidade, como por exemplo o Sikaforce® 7752. Esta dissertação tem como objetivo estudar experimentalmente e numericamente, através de modelos de dano coesivo (MDC), o comportamento de diferentes configurações de junta em T quando sujeitas a solicitações de arrancamento. Consideram-se os adesivos anteriormente mencionados para testar as juntas sob diferentes tipos de adesivos. A junta em T é constituída por 2 aderentes em L de alumínio e um aderente base também em alumínio, unidos por uma camada de adesivo. Experimentalmente é feito um estudo da resistência da junta com a variação da espessura dos aderentes em L (tP2). Com a análise numérica são estudadas as distribuições de tensões, evolução do dano, modos de rotura e resistência. Além disso, realizou-se um estudo numérico da existência ou não de adesivo de preenchimento na zona da curvatura dos aderentes em L nas tensões e na resistência da junta. Mostrouse que a variação da geometria nos aderentes em L, a presença de adesivo de preenchimento e o tipo de adesivo têm uma influência direta na resistência de junta. Os ensaios experimentais validaram os resultados numéricos e permitiram concluir que os MDC são uma técnica precisa para o estudo das geometrias das juntas em T.

Aplicação do método de Elementos Finitos eXtendido (MEFX) para a previsão da resistência de juntas adesivas de sobreposição dupla

As ligações adesivas têm sido utilizadas em diversas áreas de aplicação. A utilização das juntas adesivas em aplicações industriais tem vindo a aumentar nos últimos anos, por causa das vantagens significativas que apresentam comparativamente com os métodos tradicionais de ligação tais como soldadura, ligações aparafusadas e rebitadas. A redução de peso, redução de concentrações de tensões e facilidade de fabrico são algumas das principais vantagens das ligações adesivas. Devido à crescente utilização das ligações adesivas, torna-se necessário a existência de ferramentas que permitam prever a resistência das juntas com elevada precisão. Assim, para a análise de juntas adesivas, está a ser cada vez mais utilizado o método de Elementos Finitos. Neste âmbito o Método de Elementos Finitos eXtendido (MEFX) perfila-se como um método capaz de prever o comportamento da junta, embora este ainda não esteja convenientemente estudado no que diz respeito à aplicação a juntas adesivas. Neste trabalho é apresentado um estudo experimental e numérico pelo MEFX de juntas de sobreposição dupla, nas quais são aplicados adesivos que variam desde frágeis e rígidos, como o caso do Araldite® AV138, até adesivos mais dúcteis, como o Araldite® 2015 e o Sikaforce® 7888. Foram considerados substratos de alumínio (AW6082-T651) em juntas com diferentes comprimentos de sobreposição, sendo sujeitos a esforços de tração de forma a avaliar o seu desempenho. Na análise numérica foi realizada uma análise da distribuição de tensões na camada adesiva, a previsão da resistência das juntas pelo MEFX segundo critérios de iniciação de dano baseados em tensões e deformações, e ainda um estudo sobre o critério energético de propagação de dano. A análise por MEFX revelou que este método é bastante preciso quando usados os critérios de iniciação de dano MAXS e QUADS, e parâmetro com valor de 1 no critério energético de propagação de dano. Apesar de ser um método pouco estudado na literatura comparativamente com outros, o MEFX apresentou resultados muito satisfatórios.

Influência das tensões tectônicas na permeabilidade da formação Açu - Bacia Potiguar

The behavior of the fluid flux in oil fields is influenced by different factors and it has a big impact on the recovery of hydrocarbons. There is a need of evaluating and adapting the actual technology to the worldwide reservoirs reality, not only on the exploration (reservoir discovers) but also on the development of those that were already discovered, however not yet produced. The in situ combustion (ISC) is a suitable technique for these recovery of hydrocarbons, although it remains complex to be implemented. The main objective of this research was to study the application of the ISC as an advanced oil recovery technique through a parametric analysis of the process using vertical wells within a semi synthetic reservoir that had the characteristics from the brazilian northwest, in order to determine which of those parameters could influence the process, verifying the technical and economical viability of the method on the oil industry. For that analysis, a commercial reservoir simulation program for thermal processes was used, called steam thermal and advanced processes reservoir simulator (STARS) from the computer modeling group (CMG). This study aims, through the numerical analysis, find results that help improve mainly the interpretation and comprehension of the main problems related to the ISC method, which are not yet dominated. From the results obtained, it was proved that the mediation promoted by the thermal process ISC over the oil recovery is very important, with rates and cumulated production positively influenced by the method application. It was seen that the application of the method improves the oil mobility as a function of the heating when the combustion front forms inside the reservoir. Among all the analyzed parameters, the activation energy presented the bigger influence, it means, the lower the activation energy the bigger the fraction of recovered oil, as a function of the chemical reactions speed rise. It was also verified that the higher the enthalpy of the reaction, the bigger the fraction of recovered oil, due to a bigger amount of released energy inside the system, helping the ISC. The reservoir parameters: porosity and permeability showed to have lower influence on the ISC. Among the operational parameters that were analyzed, the injection rate was the one that showed a stronger influence on the ISC method, because, the higher the value of the injection rate, the higher was the result obtained, mainly due to maintaining the combustion front. In connection with the oxygen concentration, an increase of the percentage of this parameter translates into a higher fraction of recovered oil, because the quantity of fuel, helping the advance and the maintenance of the combustion front for a longer period of time. About the economic analysis, the ISC method showed to be economically feasible when evaluated through the net present value (NPV), considering the injection rates: the higher the injection rate, the higher the financial incomes of the final project

Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems

Análise experimental e numérica do comportamento térmico de um coletor solar acumulador

Esta tese apresenta um estudo do comportamento térmico de um coletor solar acumulador e desenvolve uma metodologia para medir a sua eficiência diária. O coletor solar acumulador está instalado na face norte do prédio de Energia Solar da UFRGS e possui cerca de 26 m2. É constituído de uma massa espessa de concreto com uma superfície absorvente feita de tijolos, possuindo uma cobertura dupla de vidros colocada de modo a deixar um espaço para a circulação de ar. Os raios solares atravessam a cobertura de vidro e aquecem a massa absorvente de tijolo, a qual aquece o ar que é introduzido no interior da construção por efeito de termossifão. Uma das principais características do coletor solar acumulador consiste no fato de que a resposta do coletor é defasada no tempo. Este fenômeno permite que o coletor entregue calor ao ambiente mesmo após o término da radiação solar. Essa defasagem dos picos de energia térmica ocorre devido ao baixo valor da difusividade térmica do concreto. Este trabalho foi dividido em duas etapas. A primeira etapa consistiu na montagem de um calorímetro para controle e monitoração das variáveis envolvidas. No interior do calorímetro foram instaladas 36 garrafas com água. As temperaturas dos conteúdos das garrafas, do coletor solar e as radiações envolvidas foram monitoradas através de 26 sensores de temperatura de CI, 8 sensores resistivos PT100 e dois sensores de radiação fotovoltaicos. Para obter as medidas dos sensores instalados foi feita a montagem de um sistema de aquisição de dados interfaceado a um microcomputador A segunda etapa consistiu na produção de um programa computacional, escrito em linguagem Fortran 90, para simular o comportamento térmico dos diversos elementos constituintes do coletor, determinar a potência térmica do coletor solar e sua eficiência diária. Para a simulação numérica do coletor solar acumulador, adotou-se um modelo simplificado bidimensional do mesmo. Foi integrada, através do Método dos Volumes Finitos, a equação de difusão de calor transiente em 2 dimensões. Na formulação das equações lineares optou-se pelo emprego das diferenças centrais no espaço e formulação explícita no tempo. Ao todo foram produzidas 4 malhas computacionais, com distintos refinamentos e foi realizado o estudo da estabilidade numérica das diversas malhas. Através da montagem experimental obtiveram-se várias características térmicas do comportamento do sistema, entre as quais, a transmitância da cobertura, curvas de temperatura do ar fornecido ao calorímetro e curva da eficiência diária do coletor solar . Através da simulação numérica foi possível determinar a potência térmica que o coletor entrega para o laboratório, a eficiência do coletor, os campos de temperatura e a vazão mássica nos diversos canais interiores do coletor solar.

Fundamentação numérica da Análise em Portugal em Anastácio da Cunha, Gomes Teixeira e Vicente Gonçalves

Neste trabalho estudamos a fundamentação numérica da Análise em Portugal, centrando particularmente este estudo nos trabalhos de José Anastácio da Cunha, Francisco Gomes Teixeira e Vicente Gonçalves. Num capítulo introdutório apresentamos uma perspectiva cronológica da procura de uma fundamentação rigorosa para a matemática, com o intuito de enquadrar historicamente as obras destes matemáticos Portugueses e reconhecer possíveis influências prestadas por trabalhos de outros autores. Relacionado com Anastácio da Cunha, analisamos os aspectos fundamentais da sua obra Principios Mathematicos, procurando evidenciar os resultados mais importantes avançados pelo autor, bem como as suas preocupações axiomáticas que não eram usuais no século XVIII, em que se insere a sua obra. Neste trabalho foi igualmente efectuada uma análise às quatro edições do Curso de Analyse Infinitesimal — Calculo Integral de Francisco Gomes Teixeira, particularmente centrada na definição do conceito de número irracional. Finalmente, analisamos o Curso de Álgebra Superior de Vicente Gonçalves, particularmente as duas últimas edições. A 2a edição do referido Curso foi objecto de duras críticas por parte de Neves Real e um dos objectivos deste trabalho foi o de procurar analisar essas críticas e verificar até que ponto influenciaram a reformulação de alguns aspectos da 3a edição.

Análise experimental e numérica do comportamento de estruturas de solo reforçado com geossintéticos : solo granular versus fino

As estruturas de solo reforçado com geossintéticos são normalmente constituídas por solos granulares com boas propriedades físicas e mecânicas. O uso de apenas este tipo de solos pode proporcionar o aumento, por vezes insustentável, do custo da execução das estruturas e o aumento do seu impacto ambiental. Deste modo, as estruturas de solo reforçado perdem a sua vantagem competitiva em relação a outros tipos de estruturas (muros de betão, muros de gravidade, muros de gabiões, etc.). Para resolver este problema podem ser utilizados outros tipos de solos (solos locais, finos, com propriedades físicas e mecânicas piores mas, no entanto, mais baratos) para a execução deste tipo de estruturas. De forma geral, com este estudo pretendeu-se contribuir para o incremento do conhecimento sobre a utilização de solos finos para a construção de estruturas de solo reforçado (muros e taludes). Para tal avaliaram-se as diferenças no comportamento mecânico dos materiais compósitos (solo granular reforçado versus solo fino reforçado) e das estruturas de solo reforçado constituídas com os dois tipos de solos. Assim, os objetivos deste estudo foram avaliar: a influência de vários parâmetros nas propriedades mecânicas e na capacidade de carga dos solos reforçados com geossintéticos; a influência de vários parâmetros no dimensionamento das estruturas de solo reforçado; e o comportamento das estruturas dimensionadas (incluindo a estabilidade global e a influência do processo construtivo) recorrendo a uma ferramenta numérica (PLAXIS). Para cumprir os objetivos propostos foram realizadas análises experimentais em laboratório (análise do comportamento do solo reforçado através de ensaios triaxiais e de California Bearing Ratio) e análises numéricas (dimensionamento de estruturas de solo reforçado; modelação numérica do comportamento através de uma ferramenta numérica comercial com o método dos elementos finitos). Os resultados dos ensaios experimentais mostraram que o comportamento mecânico e a capacidade de carga do solo foram incrementados com a inclusão das camadas de geossintético. Este efeito variou com os diversos parâmetros analisados mas, de forma geral, foi mais importante no solo fino (solo com propriedades mecânicas piores). As análises numéricas mostraram que as estruturas de solo fino precisaram de maior densidade de reforços para serem estáveis. Além disso, as estruturas de solo fino foram mais deformáveis e o efeito do seu processo construtivo foi mais importante (principalmente para estruturas de solo fino saturado).

Análise experimental e numérica de quedas de pressão em canais estenosados

O interesse nos estudos relacionados com o escoamento de fluidos em microcanais tem assumido um lugar de destaque na área de biomédica, uma vez que que estes estudos permitem analisar escoamentos sanguíneos evitando questões de ética que se colocam quando se trabalha fisicamente com o sangue. Numa primeira fase deste trabalho foram fabricados modelos anatómicos tridimensionais, mais concretamente canais circulares, um sem patologia e dois com ateromas de alturas diferentes, onde foram efetuadas medições de quedas de pressão. Com os resultados obtidos foi possível concluir que a presença de ateromas influencia o escoamento sanguíneo, verificando-se que as quedas de pressão aumentam com o aumento do grau da estenose. Numa segunda fase os escoamentos estudados experimentalmente foram alvo de um estudo numérico com o objetivo de comparar os resultados da simulação com os obtidos experimentalmente. As simulações numéricas foram efetuadas recorrendo a um software comercial que implementa o método dos volumes finitos e os resultados obtidos revelaram-se concordantes com os obtidos experimentalmente, o que permitiu concluir que a dinâmica de fluidos computacional é de facto uma boa estratégia a ser utilizada. Na fase final do trabalho, procedeu-se à identificação da melhor modelação matemática para caracterizar a diferença de valores obtidos nos estudos experimental e numérico. Após diversas estratégias conclui-se que a diferença de valores pode ser caracterizada por uma função polinomial.

Avaliação numérica de estabilidade lateral de vigas casteladas

Restrições de espaço e altura são frequentemente impostas às edificações residenciais, comerciais, industriais, depósitos e galpões com um ou diversos pavimentos em função de aspectos de regulamentos regionais, técnicos, econômicos ou ainda de natureza estética. A fim de proporcionar a passagem de tubulações e dutos de grande diâmetro sob vigas de aço, grandes alturas são normalmente requeridas, demandando por vezes, magnitudes de altura inviáveis entre pavimentos de edificações. Diversas soluções estruturais podem ser utilizadas para equacionar tais obstáculos, onde dentre outras, pode-se citar as vigas com inércia variável, stub-girders, treliças mistas, vigas misuladas e vigas com uma ou múltiplas aberturas na alma com geometrias variadas. No que tange às vigas casteladas, solução estrutural pautada neste estudo, a estabilidade é sempre um motivo de preocupação tipicamente durante a construção quando os contraventamentos laterais ainda não estão instalados. De qualquer forma, o comprimento destravado em geral alcançado pelos vãos destas vigas, são longos o suficiente para que a instabilidade ocorra. Todavia, o acréscimo substancial da resistência à flexão de tais membros devido ao aumento da altura oriundo de seu processo fabril em relação ao perfil matriz, aliada a economia de material e utilidade fim de serviço, garante a atratividade no aproveitamento destas, para grandes vãos junto aos projetistas. Não obstante, este aumento proporcional no comprimento dos vãos faz com que a instabilidade lateral ganhe importância especial. Neste contexto, o presente trabalho tem por objetivo desenvolver um modelo numérico que permita a realização de uma avaliação paramétrica a partir da calibração do modelo com resultados experimentais, efetuar a análise do comportamento de vigas casteladas e verificar seus mecanismos de falha, considerando comportamento elasto-plástico, além das não-linearidades geométricas. Também é objetivo deste trabalho, avaliar, quantificar e determinar a influência das diferenças geométricas características das vigas casteladas em relação às vigas maciças com as mesmas dimensões, analisando e descrevendo o comportamento estrutural destas vigas de aço para diversos comprimentos de vãos. A metodologia empregada para tal estudo baseou-se em uma análise paramétrica com o auxílio do método numérico dos elementos finitos.

Análise de erros da equação de advecção unidimensional no Método de Volumes Finitos

Uma análise utilizando a série de Taylor é apresentada para se estimar a priori os erros envolvidos na solução numérica da equação de advecção unidimensional com termo fonte, através do Método dos Volumes Finitos em uma malha do tipo uniforme e uma malha não uniforme. Também faz-se um estudo a posteriori para verificar a magnitude do erro de discretização e corroborar os resultados obtidos através da análise a priori. Por meio da técnica de solução manufaturada tem-se uma solução analítica para o problema, a qual facilita a análise dos resultados numéricos encontrados, e estuda-se ainda a influência das funções de interpolação UDS e CDS e do parâmetro u na solução numérica.

Modelagem numérica de elementos tracionados em aço inoxidável com parafusos defasados

Atualmente, a utilização do aço inoxidável em elementos estruturais é considerada uma solução cara para os problemas da engenharia estrutural. Todavia, mudanças de atitudes dentro da construção civil, uma transição global para um desenvolvimento sustentável e redução em impactos ambientais têm seguramente provocado um aumento na utilização do aço inoxidável. As normas de projeto de aço inoxidável atuais são, em grande parte, baseadas em analogias assumidas com o comportamento de estruturas desenvolvidas com aço carbono. Todavia, o aço inoxidável apresenta quatro curvas não-lineares tensão versus deformação (tensão e compressão, paralela e perpendicular a laminação do material), sem patamar de escoamento e região de encruamento claramente definidos, modificando assim, o comportamento global das estruturas que o utilizam. Em elementos estruturais submetidos a forças axiais de tração, a ruptura da seção líquida representa um dos estados limites últimos a serem verificados. Com o objetivo de se avaliar a resistência a tração de elementos estruturais aparafusados em aço inoxidável S304, este trabalho apresenta um modelo numérico baseado no método dos elementos finitos através do programa Ansys (versão 11). A não-linearidade do material foi considerada através do critério de plastificação de Von Mises e curvas tensão versus deformação verdadeira. A não-linearidade geométrica foi introduzida no modelo através da Formulação de Lagrange atualizado. O modelo numérico foi calibrado com resultados experimentais obtidos em ensaios de laboratório, a partir de ligações aparafusadas alternadas rígidas, onde não se ocorre nenhuma rotação entre os membros, transferindo nenhum momento fletor, apenas esforço normal e cisalhante.

Análise da resposta dinâmica experimental de uma passarela tubular mista, aço-concreto, submetida ao caminhar humano.

Este trabalho de pesquisa apresenta como objetivo principal o desenvolvimento de investigação experimental dinâmica sobre estrutura real de uma passarela tubular mista aço-concreto. O sistema estrutural objeto deste trabalho corresponde a uma passarela composta por três vãos (32,5m, 17,5m e 20,0m, respectivamente) e dois balanços (7,50m e 5,0m, respectivamente), com comprimento total de 82,5m. A passarela com estrutura contínua de aço com as ligações soldadas se apoia em quatro pórticos também de aço. Estruturalmente está constituída por duas treliças planas que se interligam através de contraventamentos horizontais fixados na corda superior e inferior da treliça e lajes de concreto, formando um sistema misto com interação completa. A estrutura está submetida correntemente à travessia de pedestres e ciclistas. Testes experimentais foram realizados sobre o sistema estrutural e confrontados com resultados numéricos. Para a modelagem numérica do sistema são empregadas técnicas usuais de discretização, via método dos elementos finitos (MEF), por meio do programa ANSYS. Os resultados experimentais são analisados de acordo com a metodologia desenvolvida, sendo realizada análise modal experimental para a determinação das propriedades dinâmicas: freqüências, modos e taxa de amortecimento, enquanto que os resultados da estrutura, em termos de aceleração de pico, são comparados com os valores limites propostos por diversos autores, normas e recomendações de projeto, para uma avaliação do desempenho da estrutura em relação a vibração quando solicitada pelo caminhar dos pedestres no que diz respeito a critério para conforto humano.

Modelagem numérica de ligações viga-coluna em aço sob momento fletor e força normal.

As ligações desempenham um papel fundamental no comportamento global das estruturas de aço. Inúmeros trabalhos de pesquisa têm sido desenvolvidos para entender o comportamento real de uma ligação e sua influência na resistência global dos pórticos. Atualmente, a Norma Brasileira de estruturas de aço de edificações, NBR 8800, considera o comportamento das ligações entre duas situações extremas: rígidas, onde não ocorre nenhuma rotação entre os membros conectados, transferindo momento fletor, força cortante e força normal; ou flexíveis, caracterizadas pela liberdade de rotação entre os membros conectados, impedindo a transmissão de momento fletor. Outras normas de projeto de estruturas de aço, consideram que as ligações apresentam um comportamento intermediário, ou seja, semi-rigídas, que podem estar submetidas a uma combinação de momento fletor e esforço normal. Porém, mesmo com a combinação, estas normas não consideram a presença de esforço normal (tração e/ou compressão). Uma limitação empírica de 5% da resistência plástica da viga é a única condição imposta no Eurocode 3. Para o estudo da ligação semi-rigída será utilizada a filosofia do Método das Componentes, que verifica a resistência da ligação e sua classificação quanto à rigidez rotacional, desenvolvida através de modelos mecânicos (modelos de molas). O objetivo deste trabalho é descrever alguns resultados de caracterização de ligações viga-coluna com placa de extremidade ajustada à altura da viga obtidos através de um modelo de elementos finitos. Para tal, será realizada uma análise não-linear geométrica e de material. Esta análise possibilitará avaliar os principais parâmetros que influenciam no comportamento deste componente no que diz respeito a sua avaliação em termos de distribuição de tensões e deformações no modelo de forma global.

Simulação numérica do processo de alteamento de áreas de deposição de resíduos pelo método a montante.

A previsão do comportamento de resíduos constituiu-se em um desafio geotécnico, uma vez que estes materiais apresentam uma resposta distinta dos materiais usualmente encontrados em depósitos naturais. A análise dos recalques da fundação, decorrentes da sobrecarga imposta pelo alteamento, é complexa, tendo em vista que o adensamento de resíduos pressupõe grandes deformações, invalidando o uso de teorias clássicas de adensamento. Atualmente, no Brasil, a técnica de disposição de resíduos de bauxita prevê uma operação inicial de lançamento no interior de lagos artificiais, em forma de polpa. Após o esgotamento do lago e ressecamento do resíduo, inicia-se o lançamento pelo método a montante. Neste método, a polpa é lançada sobre o resíduo pré-existente, que se encontra em processo de adensamento. O presente trabalho tem como objetivo reproduzir numericamente o comportamento de áreas de resíduos durante a etapa de alteamento a montante. A pesquisa tem como enfoque 2 áreas de resíduos de bauxita. Uma delas encontra-se em fase de reabilitação e dispõe de instrumentação de campo (recalques e deslocamentos horizontais). A outra se encontra em fase de operação do alteamento e dispõe de dados experimentais. Desta forma, a metodologia consistiu na reprodução numérica do processo de alteamento da área instrumentada e comparação dos resultados com a instrumentação de campo, com objetivo de avaliar o modelo numérico e os parâmetros do resíduo. Posteriormente, realizou-se a previsão do comportamento do resíduo de fundação da área em fase de alteamento. Os parâmetros geotécnicos foram definidos a partir de um extenso programa de ensaios de campo e laboratório, executado no local em estudo, fazendo-se uso de um tratamento estatístico dos dados experimentais. Os resultados numéricos mostraram a potencialidade do programa na previsão do comportamento de áreas de resíduos durante o alteamento a montante, com previsões de recalques e deslocamentos horizontais coerentes com a instrumentação de campo.