Estudo de uma equação de onda não-linear

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Controle exato para a equação de onda: o método HUM

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Migração por equação de onda em meios anisotrópicos com correção de amplitude

Este trabalho discute dois aspectos da migração em profundidade através da continuação para baixo dos campos de onda: o tratamento de modos evanescentes e a correção da amplitude dos eventos migrados. Estes dois aspectos são discutidos em meios isotrópicos e para uma classe de meios anisotrópicos. Migrações por diferenças finitas (FD) e por diferenças finitas e Fourier (FFD) podem ser instáveis em meios com forte variação lateral de velocidade. Estes métodos utilizam aproximações de Padé reais para representar o operador que descreve a propagação de ondas descendentes. Estas abordagens não são capazes de tratar corretamente os modos evanescentes, o que pode levar à instabilidades numéricas em meios com forte variação lateral de velocidade. Uma solução possível para esse problema é utilizar aproximação de Padé complexa, que consegue melhor representar os modos evanescentes associados às reflexões pós-críticas, e neste trabalho esta aproximação é utilizada para obter algoritmos FD e híbrido FD/FFD estáveis para migração em meios transversalmente isotrópicos com eixo de simetria vertical (VTI), mesmo na presença de forte variação nas propriedades elásticas do meio. A estabilidade dos algoritmos propostos para meios VTI foi validada através da resposta ao impulso do operador de migração e pela sua aplicação na migração de dados sintéticos, em meios fortemente heterogêneos. Métodos de migração por equação de onda em meios heterogêneos não tratam corretamente a amplitude dos eventos durante a propagação. As equações de onda unidirecionais tradicionais descrevem corretamente apenas a parte cinemática da propagação do campo de onda. Assim, para uma descrição correta das amplitudes deve-se usar as equações de onda unidirecionais de amplitude verdadeira. Em meios verticalmente heterogêneos, as equações de onda unidirecionais de amplitude verdadeira podem ser resolvidas analiticamente. Em meios lateralmente heterogêneos, essas equações não possuem uma solução analítica. Mesmo soluções numéricas tendem a ser instáveis. Para melhorar a compensação de amplitude na migração, em meios com variação lateral de velocidade, é proposto uma aproximação estável para solução da equação de onda unidirecional de amplitude verdadeira. Esta nova aproximação é implementada nas migrações split-step e diferenças finitas e Fourier (FFD). O algoritmo split-step com correção de amplitude foi estendido para meios VTI. A migração pré e pós-empilhamento de dados sintéticos, em meios isotrópicos e anisotrópicos, confirmam o melhor tratamento das amplitudes e estabilidade dos algoritmos propostos.

Migração usando equação da onda unidirecional de abertura arbitrária

Neste trabalho avaliamos uma classe de operadores de continuação de campos de onda, baseados em equações unidirecionais e com aplicação direta à migração sísmica. O método de representação de equações de onda unidirecionais, desenvolvido neste trabalho, é válido para abertura angular arbitrária, baseia-se no conceito de rigidez de um semiespaço, na transformação Dirichlet-Neumann e em sua discretização por elementos finitos. O método de construção dos operadores de continuação requer a introdução de variáveis auxiliares cujo número cresce em função da maior abertura angular desejada para o operador. Efetuamos a implementação no domínio do espaço e da frequência o que permite sua imediata paralelização. Baseados em experimentos numéricos, que avaliam a relação de dispersão e a resposta ao impulso do operador, propomos prescrições que permitem especificar o número de variáveis auxiliares e o passo de continuação para o operador de migração. A aplicação do algoritmo nos dados do modelo de domo salino da SEG-EAGE demonstra a capacidade do algoritmo em migrar refletores com forte mergulho em meios com forte variação lateral de velocidade.

Vibrações naturais em um sistema de interação viga-água sem o efeito de onda de superfície livre

Neste trabalho é apresentado um estudo sobre o comportamento dinâmico de sistemas flexíveis com interação viga-água. Especificamente, é escolhida uma barragem freqüentemente encontrada em problemas de engenharia de irrigação e represas. O modelo matemático utilizado para descrever tal fenômeno tem como variáveis principais a pressão hidrodinâmica e os deslocamentos horizontais da barragem. Uma vez formuladas as equações governantes de tal sistema, que resultam em uma equação de onda para a pressão e uma equação de Euler-Bernoulli para os deslocamentos, é utilizada a técnica de separação de variáveis para proceder à sua solução. Por simplicidade, assumiu-se que no modelo não havia o efeito de onda da superfície livre. Primeiro, foi resolvida a equação para a pressão, sendo calculadas as suas componentes temporais e espaciais. Depois, foi resolvida a equação de Euler-Bernoulli junto com a condição de interface. Foram simulações para o material específico, e apresentados os resultados conseguidos. No final do trabalho são enunciadas as conclusões pertinentes.

Comportamento assintótico e controlabilidade exata para a equação de Klein-Gordon

Pós-graduação em Matemática - IBILCE

Soluções analíticas da equação de Duffin-Kemmer-Petiau

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Equações de onda generalizadas e quantização funtorial para teorias de campo escalar livre

Nesta dissertação apresentamos um método de quantização matemática e conceitualmente rigoroso para o campo escalar livre de interações. Trazemos de início alguns aspéctos importantes da Teoria de Distribuições e colocamos alguns pontos de geometria Lorentziana. O restante do trabalho é dividido em duas partes: na primeira, estudamos equações de onda em variedades Lorentzianas globalmente hiperbólicas e apresentamos o conceito de soluções fundamentais no contexto de equações locais. Em seguida, progressivamente construímos soluções fundamentais para o operador de onda a partir da distribuição de Riesz. Uma vez estabelecida uma solução para a equação de onda em uma vizinhança de um ponto da variedade, tratamos de construir uma solução global a partir da extensão do problema de Cauchy a toda a variedade, donde as soluções fundamentais dão lugar aos operadores de Green a partir da introdução de uma condição de contorno. Na última parte do trabalho, apresentamos um mínimo da Teoria de Categorias e Funtores para utilizar esse formalismo na contrução de um funtor de segunda quantização entre a categoria de variedades Lorentzianas globalmente hiperbólicas e a categoria de redes de álgebras C* satisfazendo os axiomas de Haag-Kastler. Ao fim, retomamos o caso particular do campo escalar quântico livre.

Análise de guias de ondas pelos métodos vetorial magnético e dos elementos-finitos

Neste trabalho, é apresentada uma formulação apropriada à análise de guias de ondas eletromagnéticos, cobrindo do espectro de microondas até o da óptica. Nas regiões a partir do ultravioleta, os comprimentos de onda são equivalentes às dimensões atômicas e a formulação necessita de uma abordagem quântica, que não é considerada neste estudo. A formulação é fundamentada nos métodos vetorial magnético e dos elementos finitos (MEF), em meios não homogêneos, anisotrópicos e não dissipativos, embora a dissipação possa ser facilmente introduzida na análise. Deu-se preferência à formulação com o campo magnético em vez do elétrico, pelo fato do campo magnético ignorar descontinuidades elétricas. Ele é contínuo em regiões de permeabilidade homogênea, propriedade dos meios dielétricos em geral ( = 0), independente da permissividade dos respectivos meios, conquanto os campos elétricos sejam descontínuos entre regiões de permissividades diferentes.

Análise de transição entre cabos coaxiais constituída por corrugações e desníveis angulares.

Esta dissertação apresenta um formalismo baseado no Método dos Elementos Finitos (MEF), adequado a análise de descontinuidades coaxiais com simetria axial, entre duas linhas coaxiais quaisquer, incluindo corrugações nos tubos internos e externos. Pelo método de Galerkin-Budnov deduz-se um operador integral bilinear, aplicado ao campo magnético, expandido em todo o domínio da estrutura, coaxial-descontinuidade-coaxial. As portas de entrada e saída da estrutura são posicionadas distantes da descontinuidade, de forma que nelas só haja o modo Transversal Eletromagnético (TEM). O campo magnético procurado e obtido pelo MEF. Os resultados encontrados; perdas de retorno, comportamento do campo magnético e as equi-fases nas portas de entrada e saída da estrutura, foram calculados e confrontados com as do Método de Casamento de Modos (MCM), com um alto grau de concordância.

Um estudo de caso envolvendo a retroanálise de tensões residuais de cravação.

As tensões residuais de cravação ocorrem quando a estaca e o solo não são totalmente descomprimidos após dissipação completa da energia transferida à estaca durante sua instalação. Este assunto tem merecido a atenção de vários pesquisadores, uma vez que a presença das tensões residuais na estaca causa uma alteração no seu comportamento quando carregada. A presente dissertação reuniu estudos anteriores que estabeleceram os principais fatores que influenciam as tensões residuais de cravação. Procurou-se executar uma série de simulações objetivando melhor entender o desenvolvimento das tensões residuais com a porcentagem de resistência de ponta, em continuidade aos estudos de Costa (1994). Costa (1994) observou que a porcentagem de carga na ponta é um fator relevante na avaliação das tensões residuais de cravação. Uma análise paramétrica efetuada confirmou estudos anteriores, verificando que a razão entre a carga residual na ponta em relação à capacidade de carga global cresce, a medida que a porcentagem de carga na ponta aumenta, chegando a um valor máximo para, em seguida, diminuir. Este comportamento é similar ao da curva de compactação do solo, quando, ao aumentar a umidade, o peso específico seco aumenta, até um valor máximo, correspondente à umidade ótima, para em seguida reduzir. Verificouse, ainda, semelhantemente ao aumento da energia de compactação, que vai transladando a curva para cima e para à esquerda do gráfico umidade x peso específico seco, que o aumento do comprimento da estaca apresenta um comportamento similar. O aumento do comprimento leva a curva para cima e para a esquerda. Finalmente, selecionou-se um caso de obra com condições propícias para o desenvolvimento de altas tensões residuais de cravação. Trata-se de um caso de estacas metálicas longas, embutidas em solo residual jovem, de elevada resistência, bem documentado com extensa instrumentação, por ocasião da instalação. A retroanálise de cinco estacas do banco de dados mostrou que as cargas residuais previstas, num programa de simulação de cravação, se aproximaram muito dos valores experimentais. Foi observado também que a profundidade do ponto neutro previsto e medido apresentou uma excelente concordância. O resultado mais relevante desta pesquisa foi quando os valores de porcentagem de ponta foram introduzidos no eixo das abscissas e os valores da razão entre a carga residual na ponta e a capacidade de carga global, no eixo das ordenadas, e se observou o aspecto da curva semelhante à de compactação. A dissertação ilustra ser possível e simples a previsão das tensões residuais de cravação através de uma análise pela equação da onda.

Utilização de redes TIN em modelos distribuídos de precipitação/escoamento superficial

A comunicação apresenta as directivas de desenvolvimento de um modelo distribuído e a sua aplicação à Ribeira da Pradiela (distrito de Évora, afluente do Dejebe). A utilização de modelos distribuídos de precipitação/escoamento superficial para modelação integrada dos fenómenos hidrodinâmicos, transporte de sedimentos e qualidade da água à escala da bacia hidrográfica, é fundamental para uma correcta percepção dos binómios causas/efeitos e consequente gestão optimizada dos recursos hídricos. Os modelos distribuídos existentes assentam em malhas estruturadas de células regulares adjacentes (DEM) ou em malhas de triângulos irregulares adjacentes (TIN). Os primeiros têm um processamento relativamente simples e são compatíveis com a estrutura matricial de uma imagem raster, contudo obrigam a uma densidade de informação uniforme sobre a área em estudo, apresentam problemas de escala na representação das linhas de água e das variáveis distribuídas como a altitude, classes de solo, classes de uso do solo, dotações de rega, fertilizantes e pesticidas. As redes TIN são mais versáteis e no caso particular de os pontos levantados sobre o terreno formarem uma métrica regular, a rede TIN representa uma malha estruturada. Os modelos existentes que utilizam uma malha TIN definem a rede hidrográfica pelas arestas comuns dos triângulos que formam um ângulo côncavo entre si. Esta implementação pode forçar a existência de linhas de água onde não existem traços morfológicos da sua existência. A metodologia proposta para gerar as linhas de água passa por definir uma área de influência associada a cada nó (polígono de Voronoi) que drena por uma linha de água definida entre o próprio nó e o nó vizinho com o qual forma o maior declive. Esta metodologia é relativamente simples, não cria ambiguidades e permite criar uma rede hidrográfica cuja estrutura é perfeitamente compatível com a implementação de um algoritmo de diferenças finitas. As secções transversais das linhas de água são definidas como função da área a montante da respectiva secção. A precipitação efectiva é calculada à escala da bacia pelo método da curva número (CN) do Soil Conservation Service (SCS) e à escala de uma parcela pela equação de Richardson. O escoamento é resolvido pela equação da onda cinemática nas linhas de água de ordem superior, pela equação de conservação da massa nas depressões e pelas equações de Saint-Venant na linha de água principal. É considerado o destacamento nas encostas e o destacamento/ transporte/ deposição nas linhas de água. Também é considerado o transporte de poluentes adsorvidos aos sedimentos, bem como fenómenos de wash-off .

Implementação do algoritmo (RTM) para processamento sísmico em arquiteturas não convencionais

With the growth of energy consumption worldwide, conventional reservoirs, the reservoirs called "easy exploration and production" are not meeting the global energy demand. This has led many researchers to develop projects that will address these needs, companies in the oil sector has invested in techniques that helping in locating and drilling wells. One of the techniques employed in oil exploration process is the reverse time migration (RTM), in English, Reverse Time Migration, which is a method of seismic imaging that produces excellent image of the subsurface. It is algorithm based in calculation on the wave equation. RTM is considered one of the most advanced seismic imaging techniques. The economic value of the oil reserves that require RTM to be localized is very high, this means that the development of these algorithms becomes a competitive differentiator for companies seismic processing. But, it requires great computational power, that it still somehow harms its practical success. The objective of this work is to explore the implementation of this algorithm in unconventional architectures, specifically GPUs using the CUDA by making an analysis of the difficulties in developing the same, as well as the performance of the algorithm in the sequential and parallel version

Controle na fronteira para um sistema de equações de ondas

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Um estudo do teorema de unicidade de Holmgren

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)