Modelagem numérica de dados MCSEM 3D usando computação paralela

Desenvolvemos a modelagem numérica de dados sintéticos Marine Controlled Source Electromagnetic (MCSEM) usada na exploração de hidrocarbonetos para simples modelos tridimensionais usando computação paralela. Os modelos são constituidos de duas camadas estrati cadas: o mar e o sedimentos encaixantes de um delgado reservatório tridimensional, sobrepostas pelo semi-espaço correspondente ao ar. Neste Trabalho apresentamos uma abordagem tridimensional da técnica dos elementos nitos aplicada ao método MCSEM, usando a formulação da decomposição primária e secundária dos potenciais acoplados magnético e elétrico. Num pós-processamento, os campos eletromagnéticos são calculados a partir dos potenciais espalhados via diferenciação numérica. Exploramos o paralelismo dos dados MCSEM 3D em um levantamento multitransmissor, em que para cada posição do transmissor temos o mesmo processo de cálculos com dados diferentes. Para isso, usamos a biblioteca Message Passing Interface (MPI) e o modelo servidor cliente, onde o processador administrador envia os dados de entradas para os processadores clientes computar a modelagem. Os dados de entrada são formados pelos parâmetros da malha de elementos nitos, dos transmissores e do modelo geoelétrico do reservatório. Esse possui geometria prismática que representa lentes de reservatórios de hidrocarbonetos em águas profundas. Observamos que quando a largura e o comprimento horizontais desses reservatório têm a mesma ordem de grandeza, as resposta in-line são muito semelhantes e conseqüentemente o efeito tridimensional não é detectado. Por sua vez, quando a diferença nos tamanhos da largura e do comprimento do reservatório é signi cativa o efeito 3D é facilmente detectado em medidas in-line na maior dimensão horizontal do reservatório. Para medidas na menor dimensão esse efeito não é detectável, pois, nesse caso o modelo 3D se aproxima de um modelo bidimensional. O paralelismo dos dados é de rápida implementação e processamento. O tempo de execução para a modelagem multitransmissor em ambiente paralelo é equivalente ao tempo de processamento da modelagem para um único transmissor em uma máquina seqüêncial, com o acréscimo do tempo de latência na transmissão de dados entre os nós do cluster, o que justi ca o uso desta metodologia na modelagem e interpretação de dados MCSEM. Devido a reduzida memória (2 Gbytes) em cada processador do cluster do departamento de geofísica da UFPA, apenas modelos muito simples foram executados.

Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais

A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10^-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10^-1 Hz e para freqüências menores que 10^-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.

Estudos dos efeitos do manto de intemperismo sobre anomalias de VLF geradas por corpos inclinados: modelagem numérica através de programa eGs

As referências apresentadas e utilizadas através de métodos indutivos de propagação eletromagnética E.M., são utilizadas na prospecção de minerais, entre outros, haja vista que em regiões tropicais existem camadas superficiais condutivas (manto intemperismo) que levam a influenciar, evidentemente, na interpretação do resultado real. Em nossa região Amazônica, o manto apresenta-se bastante desenvolvido e condutivo. O objetivo deste trabalho é estudar os efeitos dos mantos sobre anomalias VLF, devidos a corpos tabulares inclinados em contato ou não com o manto (manto galvânico e manto indutivo). Deste modo, foram realizadas umas séries de experimentos em modelagens numéricas através da eGs, considerando-se diferentes parâmetros de resposta para o manto e para o corpo. Os trabalhos de modelagem numérica, bem como seus resultados foram apresentados por diversas situações, como: variação da profundidade do topo do condutor, variação do mergulho do condutor, variação da condutância do manto, entre outros. Para interpretação dos resultados, foi gerado um conjunto de curvas considerando-se valores pico-a-pico das anomalias de Tilt Angle e de Elipsidade. Os resultados foram sintetizados em um outro conjunto de curvas reunidas em diagramas de Argand. A medida em que se aumenta a profundidade, o pico-a-pico e a forma do pico da anomalia do Tilt Angle e da Elipsidade diminuem na presença do manto, ou na ausência do mesmo, fazendo com que este se afaste do ponto de “Cross-Over”. Para um condutor da mesma espessura, o aumento da condutância causa um aumento nas anomalias de Tilt Angle e da Elipsidade. O efeito geral na variação do mergulho do condutor é o de causar uma assimetria nos perfis de Tilt Angle e Elipsidade. O aumento da condutância do manto de intemperismo causa uma diminuição nas anomalias de Tilt Angle e de Elipsidade. Há rotação de fase no sentido anti-horário, sendo mais intensa para grandes valores de condutância do corpo. Na presença do manto, o corpo parece estar a uma profundidade inferior à verdadeira e, tendendo, a ser menos condutivo.

Migração em profundidade usando a solução numérica da equação da eiconal

Nos últimos anos tem-se verificado um interesse crescente no desenvolvimento de algoritmos de imageamento sísmico com a finalidade de obter uma imagem da subsuperfície da terra. A migração pelo método de Kirchhoff, por exemplo, é um método de imageamento muito eficiente empregado na busca da localização de refletores na subsuperficie, quando dispomos do cálculo dos tempos de trânsito necessários para a etapa de empilhamento, sendo estes obtidos neste trabalho através da solução da equação eiconal. Primeiramente, é apresentada a teoria da migração de Kirchhoff em profundidade baseada na teoria do raio, sendo em seguida introduzida a equação eiconal, através da qual são obtidos os tempos de trânsitos empregados no empilhamento das curvas de difrações. Em seguida é desenvolvido um algoritmo de migração em profundidade fazendo uso dos tempos de trânsito obtidos através da equação eiconal. Finalmente, aplicamos este algoritmo a dados sintéticos contendo ruído aditivo e múltiplas e obtemos como resultado uma seção sísmica na profundidade. Através dos experimentos feitos neste trabalho observou-se que o algoritmo de migração desenvolvido mostrou-se bastante eficiente e eficaz na reconstrução da imagem dos refletores.

Análise numérica de seções reforçadas de concreto armado considerando o estado inicial de tensões

O presente trabalho trata da formulação e da implementação computacional, em MATLAB^®, para a análise numérica de seções reforçadas de concreto armado, submetidas à flexão composta, considerando o estado de tensões anterior ao reforço. A referida análise se dá com a geração de diagramas de interação momento fletor esforço normal por dois métodos, quais sejam: a) varredura dos domínios de deformação da NBR6118/2003; b) determinação dos picos de diagramas momento fletor – esforço normal – curvatura. Em ambos os procedimentos utiliza-se o método numérico do ponto médio na integração do cálculo dos esforços resistentes, e o método de Newton Raphson para a obtenção de raízes é usado na determinação da deformação no eixo de referência da seção, durante a determinação dos diagramas momento fletor -esforço normal - curvatura. Preliminarmente, concluiu-se que a primeira das duas metodologias aplicadas é inviável. Posteriormente, com a confirmação da eficácia da segunda metodologia, foi possível expandir o escopo do trabalho de modo a permitir a análise de seções de formatos quaisquer executadas em várias etapas, considerando o estado de tensões inicial em cada uma das etapas. A implementação computacional referente a este trabalho se baseou no programa para análise numérica de seções SECLAB, desenvolvido pelo professor Remo Magalhães de Souza.

Simulação numérica do processo de secagem de madeiras

Pós-graduação em Engenharia Mecânica - FEG

Desenvolvimento de uma metodologia numérica para escoamentos viscoelásticos não-isotérmicos

Pós-graduação em Matematica Aplicada e Computacional - FCT

Desenvolvimento experimental e simulação numérica de uma réplica térmica de placa de circuito impresso

Pós-graduação em Engenharia Mecânica - FEIS

AVALIAÇÃO ROBUSTA DA MODELAGEM NEUTROSFÉRICA SOBRE O TERRITÓRIO BRASILEIRO BASEADA EM MODELOS DE PREVISÃO NUMÉRICA DE TEMPO DA AMÉRICA DO SUL

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Simulação numérica através do método dos elementos finitos e comparação experimental para obtenção de tensões e deformações específicas em fêmur com prótese total de quadril

Pós-graduação em Engenharia Mecânica - FEG

Modelagem e simulação numérica de um sistema de refrigeração por absorção utilizando o par H2O-LiBr

The term refrigeration solar refers to any air conditioning system that uses solar energy as a primary energy source. The use of solar radiation for cooling purposes is divided according to their technological possibilities which are distinguished from one another as the way that energy is involved in the cycle, work or heat. The first case is related to vapor compression cycles, in which the work input is provided by the photovoltaic conversion of solar energy into electrical energy. In the second case, an absorption refrigeration cycle is used and the thermal energy collected from the solar radiation is provided at the generator of this cycle.. In this work a system with an absorption cycle using the pair BrLi-water, using solar energy as input is modeled. It is considered a simple refrigeration cycle whose the equations of mass and energy conservation in each component are developed in order to obtain an algebraic equation set and a simulation routine using the EES software. Although the simulation operates under certain specified thermal load it is possible to estimate the necessary areas of heat exchangers and solar collectors

Simulação numérica de fluxo de águas subterrâneas na bacia do Rio Grande (BA)

The Urucuia Aquifer System represents a strategic water source in western Bahia. Its baseflow is responsible for the flow rate of the main tributaries of São Francisco river left bank in the dry season, including the Rio Grande, its main tributary in Bahia state. This river has a hydrological regime heavily affected by groundwater and is located in a region with conflicts about water resources. The aquifers geology is constituted by neocretacious sandstones of Urucuia Group subdivided in Posse Formation and Serra das Araras Formation. The embasement is formed by neoproterozoic rocks of Bambuí Group. This work focuses on an important tool application, the mathematical model, whose function is represent approximately and suitably the reality so that can assist in different scenarios simulations and make predictions. Many studies developed in this basin provided the conceptual model basis including a full free aquifer, lithological and hydraulical homogeneity in entire basin, null flux at plateau borders and aquifer base. The finite element method is the numerical method used and FEFLOW the computational algorithm. The simulated area was discretized in a single layer with 27.357,6 km² (314.432 elements and 320.452 nodes) totaling a 4.249,89 km³ volume. Were utilized 21 observation wells from CERB to calibrate the model. The terrain topography was obtained by SRTM data and the impermeable base was generated by interpolation of descriptive profiles from wells and electric vertical drilling from previous studies. Works in this area obtained mean recharge rates varying approximately from 20% to 25% of average precipitation, thus the values of model recharge zones varying in this range. Were distributed 4 hydraulic conductivity zones: (K1) west zone with K=6x10-5 m/s; (K2) center-east zone with K=3x10-4 m/s; (K3) far east zone with K=5x10-4 m/s; e (K4) east-north zone with K=1x10-5 m/s. Thereby was incorporated to the final conceptual model...

Estudo da fluência dos aços inoxidáveis ferríticos estabilizados ao titânio e nióbio e simulação numérica do processo no abaqus

The purpose of this research was tested a finite element model (FEM) that represented the creep of a slab during the reheating process of hot rolling. The aim is to prevent creep phenomenon changing the reheating profile with hot tensile test in Gleeble 3500, and, also, understand the former defect crisis. The goal of this work is to have a predictive tool to optimize the reheating process changing parameters (length and thickness). Then, use input parameters obtained from the tests to approximate the solution of the problem aided by Abaqus CAE. The results have showed that the ferritic stainless steel AISI 409 has a lower sensitivity to creep comparing to the stainless steel AISI 409, AISI 430Ti, AISI 441 and AISI 444

Simulação numérica de fluxo de águas subterrâneas na fazenda Cara Preta - Santa Rita do Passa Quatro - SP

The term model refers to any representation of a real system. The use of models in Hydrogeology can be valuable predictive tools for management of groundwater resources. The numeric models of groundwater flow, object of this study, consist on a set of differential equations that describe the water flow in the porous medium. In this context, numeric simulations were made for a sub-basin located at Cara Preta farm – Santa Rita do Passa Quatro – SP. The aquifer at the local is composed by rocks of Pirambóia Formation, which is part of Guarani Aquifer System. It was developed a conceptual model from previous studies in the area, and from that, simulations were made through the software Visual Modflow®. The conceptual model established previously was considered consistent through the results of simulation.

Simulação geoestatística e classificação numérica para o delineamento de unidade de mapeamento de atributos e erosão do solo

Pós-graduação em Agronomia (Ciência do Solo) - FCAV