Desenvolvimento de Interface gráfica como suporte para soluções numéricas das equações de Maxwell em coordenadas gerais – 3D

Neste trabalho, é implementada uma interface gráfica de usuários (GUI) usando a ferramenta Qt da Nokia (versão 3.0). A interface visa simplificar a criação de cenários para a realização de simulações paralelas usando a técnica numérica Local Nonorthogonal Finite Difference Time-Domain (LN-FDTD), aplicada para solucionar as equações de Maxwell. O simulador foi desenvolvido usando a linguagem de programação C e paralelizado utilizando threads. Para isto, a biblioteca pthread foi empregada. A visualização 3D do cenário a ser simulado (e da malha) é realizada por um programa especialmente desenvolvido que utiliza a biblioteca OpenGL. Para melhorar o desenvolvimento e alcançar os objetivos do projeto computacional, foram utilizados conceitos da Engenharia de Software, tais como o modelo de processo de software por prototipagem. Ao privar o usuário de interagir diretamente com o código-fonte da simulação, a probabilidade de ocorrência de erros humanos durante o processo de construção de cenários é minimizada. Para demonstrar o funcionamento da ferramenta desenvolvida, foi realizado um estudo relativo ao efeito de flechas em linhas de baixa tensão nas tensões transitórias induzidas nas mesmas por descargas atmosféricas. As tensões induzidas nas tomadas da edificação também são estudadas.

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.

Solução híbrida da equação de advecção-dispersão em meios porosos

Este trabalho consiste na solução híbrida da Equação de Advecção-dispersão de solutos unidimensional em meios porosos homogêneos ou heterogêneos, para um único componente, com coeficientes de retardo, dispersão, velocidade média, decaimento e produção dependentes da distância percorrida pelo soluto. Serão estudados os casos de dispersão-advecção em que o retardamento, dispersão, velocidade do fluxo, decaimento e produção variem de forma linear enquanto a dispersividade assuma os modelos linear, parabólico ou exponencial. Para a solução da equação foi aplicada a Técnica da Transformada Integral Generalizada. Os resultados obtidos nesta dissertação demonstram boa concordância entre os problemas-exemplo e suas soluções numéricas ou analíticas contidas na literatura e apontam uma melhor adequação no uso de modelos parabólico no estudo da advecção-dispersão em curto intervalo de tempo, enquanto que o modelo linear converge mais rapidamente em tempos prolongados de simulação. A convergência da série mostrou-se ter dependência direta quanto ao comprimento do domínio, ao modelo de dispersão e da dispersividade adotada, convergindo com até 60 termos, podendo chegar a NT = 170, para os casos heterogêneos, utilizando o modelo de dispersão exponencial, respeitando o critério adotado de 10^-4.

Simulação de perfis nucleares de poço em formações complexas

A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P₁, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm³ e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e^-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas.

Modelagens direta e inversa de dados mCSEM 3D

Marine Controlled Source Electromagnetic - mCSEM é um método geofísico eletromagnético que nos últimos dez anos vem sendo usado na prospecção de hidrocarbonetos com bastante êxito. Este método consiste em um dipolo elétrico horizontal (DEH) localizado um pouco acima do assoalho marinho, operando em baixa frequência (0,1-1,0 Hz) e receptores regularmente distribuídos no fundo do mar que captam os campos eletromagnéticos provenientes da difusão de energia gerada pelo dipolo transmissor. Neste trabalho vamos apresentar o problema direto do método mCSEM 3D, propondo soluções numéricas, através do método dos elementos finitos tridimensionais, para modelos geoelétricos mCSEM 3D. Para fins de análise de coerência, os resultados obtidos são comparados com soluções disponíveis na literatura. Em seguida, apresentaremos a inversão de um de seus modelos segundo uma proposta de metodologia de inversão juntamente com a proposta de solução direta para o mCSEM 3D, acima mencionada, realizando assim a inversão de um modelo geoelétrico do mCSEM 3D para duas frequências.

A conversão de registros de representações semióticas no estudo de sistemas de equações algébricas lineares

Este trabalho consiste na proposta de uma sequencia didática para o ensino de Sistemas de Equações Algébricas Lineares na qual estabelecemos uma conexão entre o Método da Substituição e o buscando a conversão de registros de representação. O objetivo da proposta foi verificar se os alunos conseguem realizar a conexão entre os dois métodos desenvolvendo a conversão do método da substituição no Método do escalonamento caracterizando assim, o aprendizado do objeto matemático estudado, segundo a teoria de registros de representação semiótica de Raimund Duval. A pesquisa foi realizada com alunos do ensino médio em uma escola da rede pública estadual da cidade de Belém e os resultados apontaram para o estabelecimento de uma conexão entre os dois métodos empregados no processo de resolução de sistemas.

Sólitons latentes, cordas negras, membranas negras e equações de estado em modelos de Kaluza-Klein

Neste trabalho investigamos soluções solitônicas em modelos de Kaluza-Klein com um número arbitrário de espaços internos toroidais, que descrevem o campo gravitacional de um objeto massivo compacto. Cada toro d_i-dimensional possui um fator de escala independente C_i, i = 1, ..., N, que é caracterizado pelo parâmetro ᵞi. Destacamos a solução fisicamente interessante correspondente à massa puntual. Para a solução geral obtemos equações de estado nos espaços externo e interno. Estas equações demonstram que a massa pontual solitônica possui equações de estado tipo poeira em todos os espaços. Obtemos também os parâmetros pósnewtonianos que nos possibilitam encontrar as fórmulas da precessão do periélio, do desvio da luz e do atraso no tempo de ecos de radar. Além disso, os experimentos gravitacionais levam a uma forte limitação nos parâmetros do modelo: T = ƩNi=1 diYi = −(2, 1±2, 3)×10⁻⁵. A solução para massa pontual com Y1 = . . . = YN = (1+ƩNi=1 di)⁻¹ contradiz esta restrição. A imposição T = 0 satisfaz essa limitação experimental e define uma nova classe de soluções que são indistinguíveis para a relatividade geral. Chamamos estas soluções de sólitons latentes. Cordas negras e membranas negras com Yi = 0 pertencem a esta classe. Além disso, a condição de estabilidade dos espaços internos destaca cordas/membranas negras de sólitons latentes, conduzindo exclusivamente para as equações de estado de corda/membrana negra p_i = −ε/2, i = 1, . . . ,N, nos espaços internos e ao número de dimensões externas d₀ = 3. As investigações do fluido perfeito multidimensional estático e esfericamente simétrico com equação de estado tipo poeira no espaço externo confirmam os resultados acima.

Interpretação de anomalias magnéticas usando integrais finitas dos momentos das componentes

O objetivo do trabalho consiste em desenvolver uma metodologia para determinar os parâmetros de um corpo causador de uma anomalia magnética. Essa metodologia baseia-se na utilização de expressões analíticas, deduzidas para as integrais finitas dos momentos da anomalia de intensidade total e das componentes H_x, H_y e H_z ao longo dos três eixos de um sistema de coordenadas cartesianas. Por meio do ajuste entre essas expressões analíticas e as integrais numéricas finitas dos momentos das componentes magnéticas obtidas a partir de um levantamento geofísico, pode-se computar, através de um processo iterativo, os parâmetros de magnetização, posição, profundidade e dimensões de um corpo anômalo. No caso em que são conhecidas as medidas de somente uma componente magnética, ainda é possível aplicar o método pois as demais componentes podem ser obtidas por um esquema de filtragem matemática. A metodologia foi testada com sucesso para modelos dipolo pontual e linha de dipolos os quais são muito utilizados em Geofísica para interpretar anomalias magnéticas produzidas por corpos geológicos que possuem uma ou duas dimensões horizontais muito menores do que a profundidade.

Inversão gravimétrica rápida do relevo do embasamento aplicando o funcional da variação total

Apresentamos um novo método de inversão linear bidimensional de dados gravimétricos produzidos por bacias sedimentares com relevo do embasamento descontínuo. O método desenvolvido utiliza um modelo interpretativo formado por um conjunto de fitas horizontais bidimensionais justapostas cujas espessuras são os parâmetros a serem estimados. O contraste de densidade entre o embasamento e os sedimentos é presumido constante e conhecido. As estimativas das espessuras foram estabilizadas com o funcional da Variação Total (VT) que permite soluções apresentando descontinuidades locais no relevo do embasamento. As estimativas do relevo são obtidas através da resolução de um sistema de equações lineares, resolvido na norma L₁. Como métodos lineares subestimam as estimativas de profundidade do embasamento de bacias maiores que cerca de 500 m, amplificamos as estimativas de profundidade através da modificação da matriz associada ao modelo interpretativo de fitas. As estimativas obtidas através deste procedimento são em geral ligeiramente superestimadas. Desse modo, elas são corrigidas através de uma correção definida pela expressão da placa Bouguer. Testes em dados sintéticos e reais produziram resultados comparáveis aos produzidos pelo método não linear, mas exigiram menor tempo computacional. A razão R entre os tempos exigidos pelo método não linear e o método proposto cresce com o número de observações e parâmetros. Por exemplo, para 60 observações e 60 parâmetros, R é igual a 4, enquanto para 2500 observações e 2500 parâmetros R cresce para 16,8. O método proposto e o método de inversão não linear foram aplicados também em dados reais do Steptoe Valley, Nevada, Estados Unidos, e da ponte do POEMA, no Campus do Guamá em Belém, produzindo soluções similares às obtidas com o método não linear exigindo menor tempo computacional.

Mapeamento das fundações mais usadas na cidade de Belém-PA: aspectos gerais e proposta preliminar de mapeamento de soluções utilizadas em casos recentes

Realiza o levantamento dos vários tipos de fundações mais comumente utilizadas na construção civil, nas últimas décadas na cidade de Belém PA, definindo suas propriedades, peculiaridades, aspectos técnicos de seu dimensionamento, capacidade de carga e processo executivo. Identifica o tipo de solo, em que estão assentes, com respectivos perfil e características geotécnicas e efetua o mapeamento da cidade por regiões, de acordo com perfil geológico; tipo de fundação; profundidade de assentamento; processo executivo; e custo relativo. O trabalho é dividido em quatro etapas: a primeira correspondente à pesquisa bibliográfica a respeito do estudo das fundações, abrangendo desenvolvimento histórico das fundações e mecânica dos solos, no mundo, Brasil e Belém; pesquisa geotécnica e características do solo de Belém e; estudo dos Sistemas de Informações Geográficas SIG,; a segunda envolvendo levantamento técnico junto as empresas construtoras e firmas de projetos locais, para catalogação e formação de um banco de dados acerca das fundações praticadas na cidade; a terceira relativa a ordenação e análise técnica dos dados coletados e a quarta, consiste na redação do texto final e na elaboração de mapas da cidade, baseados em Sistema de Informação Geográfica - SIG, em escala 1:10.000, concernentes ao tipo de fundação; o perfil geotécnico do terreno; o tipo de solo e processo executivo; profundidade provável para a fundação e, o custo relativo da fundação em relação ao total da obra.

Probabilidade de emergência de classes ordinais após o ensino de relações numéricas

O paradigma de equivalência de estímulos tem se mostrado útil na explicação de processos comportamentais complexos, como aqueles envolvidos em comportamentos conceituais numéricos. Vários estudos têm buscado a compreensão de como desempenhos sob controle da função de ordem são estabelecidos e mantidos. O objetivo do presente trabalho foi verificar se classes ordinais poderiam emergir após o ensino por emparelhamento arbitrário e de produção de seqüência. Participaram do estudo três alunos com atraso no desenvolvimento. Os estímulos visuais foram formas abstratas indicando numerosidade (A), numerais cardinais (B)e nomes escritos em letras maiúsculas de numerais (C). As sessões experimentais foram conduzidas em uma sala da APAE-BELÉM e um software controlou e registrou os dados comportamentais. As relações AB/AC foram ensinadas e testou-se a emergência de três classes de equivalência. Em seguida, houve um ensino por encadeamento de respostas com estímulos de um dos conjuntos (A1A2A3) e uma sonda de seqüenciação. Então, foi avaliada a emergência de novas seqüências (B1B2B3 e C1C2C3). Posteriormente, testes de substitutabilidade foram aplicados para verificar a formação de classes ordinais (por exemplo: A1B2C3). Testes de generalização também foram apresentados para verificar se um responder envolvendo numerosidade ocorreria com novos estímulos (por exemplo: E1E2E3). Os resultados demonstraram que os participantes responderam a novas seqüências prontamente ou com emergência gradual. A análise de topografias de controle de estímulos envolvidas nesse tipo de tarefa mostrou-se útil para a compreensão da ordinalidade. Todos os participantes responderam a seqüências com novos estímulos (generalização). O procedimento mostrou-se também eficiente na transferência de funções ordinais em pessoas com atraso no desenvolvimento.

Emergência de relações numéricas em crianças surdas

Os estudos sobre a formação de classes ordinais têm apresentado diversos achados experimentais relevantes para a compreensão das relações entre estímulos em seqüências. O objetivo deste trabalho foi o de replicar sistematicamente os resultados obtidos com crianças e adultos normais (Estudo 1) usando um procedimento de ensino por encadeamento e verificar se esses resultados se manteriam sob controle condicional (Estudo 2) com crianças surdas, usando apenas estímulos visuais, através de um ensino de encadeamento de respostas motoras. Participaram do Estudo 1 sete crianças com grau profundo de surdez, na fase de alfabetização em língua de sinais. Foi programado um procedimento de encadeamento de respostas com três conjuntos de estímulos: A nomes impressos dos números, B numerais em língua brasileira de sinais (LIBRAS), C formas abstratas em quantidades diferentes. Para os três tipos de estímulos os valores eram de 1 a 6. No Estudo 2 cinco outras crianças surdas foram submetidas ao mesmo procedimento de ensino, porém, o participante deveria responder na presença da cor vermelha na seqüência A1A2A3 A4A5A6 e na presença da cor verde na seqüência 6A5A4A3A2A1. Nos dois estudos, após responder corretamente a cada seqüência, uma animação gráfica era apresentada na tela, a experimentadora elogiava com as palavras muito bem, certo, legal em linguagem de sinais. Caso o participante emitisse qualquer outra resposta, a tela se embranquecia por um segundo e uma nova configuração dos estímulos era apresentada aleatoriamente. As sessões experimentais eram realizadas na sala de aula de uma instituição pública especializada. Foi utilizado um Laptop e um software especialmente desenvolvido para coleta de dados comportamentais. Após treino e revisão da linha de base com as três seqüências de respostas, testes eram aplicados para avaliar transitividade e conectividade na emergência de classes ordinais através da disposição de pares de estímulos não 14 adjacentes e de substitutabilidade de estímulos. Os resultados mostraram que os participantes responderam prontamente em ambos os estudos. Conclui-se que o procedimento é eficiente na aquisição de comportamentos conceituais numéricos e que a emergência de classes de estímulos equivalentes também ocorre fora do formato de emparelhamento de acordo com o modelo, mesmo em crianças surdas.

Análise de tensões induzidas em linhas de distribuição de baixa tensão frente a uma descarga atmosférica

Neste trabalho são apresentadas simulações computacionais inéditas para o cálculo de tensões induzidas em linhas de baixa tensão provenientes de descargas atmosféricas em estações rádio-base de telefonia celular (ERBs). Foram construídas estruturas representativas que denotam um grau de complexidade bastante avançado e semelhante ao encontrado em campo, visando assim a obtenção o de resultados bem próximos aos da realidade. Para tal, desenvolveu-se um software, no qual as equações de Maxwell são resolvidas numericamente utilizando o Método das Diferenças Finitas no Domínio do Tempo (FDTD), associado à truncagem de domínio de análise pela técnica da UPML e representação de condutores elétricos pela formulação de fio fino para meios condutivos, gerando soluções de onda completa para o problema.

Provendo soluções de comunicação em ambientes industriais inóspitos via PLC - um estudo de caso aplicado a uma empilhadeira de bauxita da refinaria da Alunorte

Este trabalho tem como principal objetivo fazer um estudo sobre um sistema comunicação de dados, utilizando-se a rede elétrica existente entre equipamentos, como meio físico em um ambiente industrial inóspito. Nesse caso, foi implementado um piloto em uma refinaria de alumina, localizada em Barcarena no estado do Pará, em uma máquina móvel de empilhamento de bauxita, que possuía um sistema de comunicação de dados via rádio, no qual ocorriam muitas falhas de comunicação devido às interferências eletromagnéticas, RFI, etc. prejudicando o processo de produção da refinaria. Como alternativa, apresentar-se-á uma tecnologia menos susceptível a problema da interferência que permite melhorar o desempenho do sistema de comunicação entre equipamentos. Trata-se da tecnologia PLC - Power Line Communication- que é full-duplex e utiliza modulação OFDM, o que reduz o efeito indesejável de interferências. Na conclusão, foram feitos comentários pertinentes aos resultados obtidos e os esperados acerca da confiabilidade da tecnologia, bem como sua utilização para futuras aplicações, nas diversas áreas onde a mesma pode ser utilizada.

Análise da área desmatada municipal na Amazônia brasileira no período 2000-2004: uma abordagem com modelos não lineares

Existe uma certa assimetria na distribuição do desmatamento da Amazônia Legal brasileira, concentrando-se particularmente nos estados do Mato Grosso, Pará e Rondônia, o que estaria, em princípio, associado a certas características comuns de suas atividades produtivas predominantes, como a pecuária, a exploração madeireira e, mais recentemente, a expansão da produção de grãos. Todavia, pode-se indagar se é possível a formação de grupos de economias municipais com características semelhantes, quanto às causas ligadas ao desmatamento. Com este objetivo em mente, procura-se neste artigo testar a hipótese de convergência em clubes a partir das áreas desmatadas dos municípios da Amazônia Legal no período de 2000 a 2004. A fim de alcançar o objetivo proposto foi utilizado um modelo com efeito limiar (threshold) baseado em Hansen (2000). Os resultados empíricos corroboram a hipótese levantada, no que concerne a formação de quatro clubes de convergência.