Inversão de momentos de fonte em métodos potenciais

A inversão de momentos de fonte gravimétrica tridimensional é analisada em duas situações. Na primeira se admite conhecer apenas a anomalia. Na segunda se admite conhecer, além da anomalia, informação a priori sobre o corpo anômalo. Sem usar informação a priori, mostramos que é possível determinar univocamente todo momento, ou combinação linear de momentos, cujo núcleo polinomial seja função apenas das coordenadas Cartesianas que definem o plano de medida e que tenha Laplaciano nulo. Além disso, mostramos que nenhum momento cujo núcleo polinomial tenha Laplaciano não nulo pode ser determinado. Por outro lado, informação a priori é implicitamente introduzida se o método de inversão de momentos se baseia na aproximação da anomalia pela série truncada obtida de sua expansão em multipolos. Dado um centro de expansão qualquer, o truncamento da série impõe uma condição de regularização sobre as superfícies equipotenciais do corpo anômalo, que permite estimar univocamente os momentos e combinações lineares de momentos que são os coeficientes das funções-bases da expansão em multipolos. Assim, uma distribuição de massa equivalente à real é postulada, sendo o critério de equivalência especificado pela condição de ajuste entre os campos observado e calculado com a série truncada em momentos de uma ordem máxima pré-estabelecida. Os momentos da distribuição equivalente de massa foram identificados como a solução estacionária de um sistema de equações diferenciais lineares de 1a. ordem, para a qual se asseguram unicidade e estabilidade assintótica. Para a série retendo momentos até 2a. ordem, é implicitamente admitido que o corpo anômalo seja convexo e tenha volume finito, que ele esteja suficientemente distante do plano de medida e que a sua distribuição espacial de massa apresente três planos ortogonais de simetria. O método de inversão de momentos baseado na série truncada (IMT) é adaptado para o caso magnético. Para este caso, mostramos que, para assegurar unicidade e estabilidade assintótica, é suficiente pressupor, além da condição de regularização, a condição de que a magnetização total tenha direção e sentido constantes, embora desconhecidos. O método IMT baseado na série de 2a. ordem (IMT2) é aplicado a anomalias gravimétricas e magnéticas tridimensionais sintéticas. Mostramos que se a fonte satisfaz as condições exigidas, boas estimativas da sua massa ou vetor momento de dipolo anômalo total, da posição de seu centro de massa ou de momento de dipolo e das direções de seus três eixos principais são obtidas de maneira estável. O método IMT2 pode falhar parcialmente quando a fonte está próxima do plano de medida ou quando a anomalia tem efeitos localizados e fortes de um corpo pequeno e raso e se tenta estimar os parâmetros de um corpo grande e profundo. Definimos por falha parcial a situação em que algumas das estimativas obtidas podem não ser boas aproximações dos valores verdadeiros. Nas duas situações acima descritas, a profundidade do centro da fonte (maior) e as direções de seus eixos principais podem ser erroneamente estimadas, embora que a massa ou vetor momento de dipolo anômalo total e a projeção do centro desta fonte no plano de medida ainda sejam bem estimados. Se a direção de magnetização total não for constante, o método IMT2 pode fornecer estimativas erradas das direções dos eixos principais (mesmo se a fonte estiver distante do plano de medida), embora que os demais parâmetros sejam bem estimados. O método IMT2 pode falhar completamente se a fonte não tiver volume finito. Definimos por falha completa a situação em que qualquer estimativa obtida pode não ser boa aproximação do valor verdadeiro. O método IMT2 é aplicado a dados reais gravimétricos e magnéticos. No caso gravimétrico, utilizamos uma anomalia situada no estado da Bahia, que se supõe ser causada por um batólito de granito. Com base nos resultados, sugerimos que as massas graníticas geradoras desta anomalia tenham sido estiradas na direção NNW e adelgaçadas na direção vertical durante o evento compressivo que causou a orogênese do Sistema de Dobramentos do Espinhaço. Além disso, estimamos que a profundidade do centro de massa da fonte geradora é cerca de 20 km. No caso magnético, utilizamos a anomalia de um monte submarino situado no Golfo da Guiné. Com base nos resultados, estimamos que o paleopolo magnético do monte submarino tem latitude 50°48'S e longitude 74°54'E e sugerimos que não exista contraste de magnetização expressivo abaixo da base do monte submarino.

Modelagem e simulação do processo de separação por membrana difusa polimérica para a extração de metais

O processo de extração de metais em membrana polimérica é analisado empregando o método das Linhas-Gear nas equações diferenciais resultantes da modelagem matemática do problema físico. O estudo tem a finalidade de averiguar o comportamento das curvas de extração de metal extraído e a influência dos parâmetros físicos da membrana no processo de extração por membrana, e também analisar, a partir de um estudo paramétrico, a influência dos parâmetros mais importantes do processo, verificando-se assim se o processo de extração se encontra sob controle cinético ou difusivo. A modelagem do problema é realizada assumindo-se a extração dos metais cádmio, Cd (II), ouro, Au (lII), e Paládio, Pd (11) sobre a membrana difusiva polimérica Aliquat 336/PVC. Os resultados são então comparados com dados obtidos em trabalhos previamente reportados na literatura, concluindo-se que os mesmos são satisfatórios.

Modelagem matemática nas aulas de cálculo: uma estratégia que pode contribuir com a aprendizagem dos alunos de engenharia

Este trabalho objetiva analisar os possíveis efeitos que o uso da Modelagem Matemática, enquanto estratégia de ensino, provoca no processo de aprendizagem dos alunos da disciplina Cálculo III – EDO (Equações Diferenciais Ordinárias). A pesquisa foi desenvolvida em uma turma de alunos do 2° ano do curso de Engenharia da Computação, na Universidade Federal do Pará. O trabalho é de cunho qualitativo onde foram levados em consideração os aspectos sociais que permeiam uma sala de aula universitária. Importante destacar que houve a participação direta da professora-pesquisadora de Matemática. Para que eu pudesse fazer a coleta dos dados, utilizei alguns instrumentos que considerei essenciais, tais como: observações, gravações em áudio, questionários semiestruturados e registros escritos dos alunos. De posse de alguns resultados preliminares, me foi possível observar o quanto a Modelagem Matemática desempenha um papel relevante na aprendizagem dos conteúdos matemáticos por parte dos alunos, pois foi possível eles interagirem com outras áreas do conhecimento sendo, desta forma, estimulados a realizarem pesquisa e, simultaneamente, serem parte do processo de ensino e aprendizagem que foi gerado no ambiente de sala de aula. Observei, também, que a utilização da Modelagem Matemática, enquanto estratégia de ensino e aprendizagem, conduziu os alunos a despertarem para os aspectos reflexivos e críticos até então adormecidos, uma vez que são necessários para uma aprendizagem com qualidade para, assim, construírem seus conhecimentos acadêmicos e profissionais.

Avaliação da influência do eletrojato equatorial na interpretação de dados magnetotelúrico bidimensionais

O campo magnetotelúrico em regiões equatoriais viola a condição de ondas planas por causa de uma fonte fortemente concentrada na direção E-W na ionosfera, denominada eletrojato equatorial. No presente trabalho, procurou-se analisar a resposta magnetotelúrica de fontes que simulam o efeito do eletrojato equatorial. Foram considerados dois tipos de fontes para simular o eletrojato: uma linha infinita de corrente e uma distribuição gaussiana de densidade de corrente em relação a uma das coordenadas horizontais. A resistividade aparente foi obtida através da relação de Cagniard e comparada com os resultados de ondas planas. É mostrada também a comparação entre a fase da impedância na superfície, para os três tipos de fontes (ondas planas, eletrojato gaussiano e linha de corrente). O problema de meios com heterogeneidades laterais foi resolvido em termos de campos secundários, sendo as equações diferenciais solucionadas através da técnica de elementos finitos bidimensionais. Os resultados mostram que o eletrojato tem pouca influência nas respostas (resistividade aparente e fase) de estruturas geológicas rasas. Entretanto, a influência pode ser considerável nas estruturas profundas (maior que 5000 m), principalmente se suas resistividades são altas (maior que 100 Ω.m). Portanto, a influência do eletrojato equatorial deve ser considerada na interpretação de dados magnetotelúricos de bacias sedimentares profundas ou no estudo da crosta terrestre.

A transmission problem for the Timoshenko system

In this work we study a transmission problem for the model of beams developed by S.P. Timoshenko [10]. We consider the case of mixed material, that is, a part of the beam has friction and the other is purely elastic. We show that for this type of material, the dissipation produced by the frictional part is strong enough to produce exponential decay of the solution, no matter how small is its size. We use the method of energy to prove exponential decay for the solution.

Modelagem computacional de dados magnetotelúricos marinhos 2-D

O presente trabalho tem como objetivo investigar o espalhamento das ondas planas causados pela variação lateral das propriedades físicas das rochas. Essa é uma das questões mais importante para o sucesso da exploração geofísica. O método geofísico usado nesta dissertação foi método magnetotelúrico (MT) e método Magnetotelúrico marinho (mMT). Umas das ferramentas utilizadas foram o método dos elementos finitos, este método é notável pela flexibilidade em resolver numericamente as equações diferenciais para campos eletromagnéticos secundários, gerados por estruturas geológicas com geometrias complexas, sendo utilizado para obter a solução numérica do espalhamento do campo eletromagnético da onda plana por um dique vertical (modelos bidimensionais). Os algoritmos desenvolvidos foram aplicados em diversos modelos de dique vertical, com diferentes parâmetros geoelétricos.

Detecção e localização de descargas parciais em transformadores de potência utilizando sinais acústicoas transitórios e o método GSPS

Neste trabalho, é desenvolvido um método de localização de descargas parciais, em transformadores de potência, baseado no algoritmo GPS (Global Positioning System). Para a análise da estrutura, foi desenvolvido um solftware, no qual as equações diferenciais que representam a propagação de ondas acústicas são resolvidas numericamente através do método Acoustic Finite Difference Time Domain (AFDTD), cujo domínio computacional é truncado através da técnica CPML (Convolutional Perfectly Matched Layer). Os resultados obtidos são comparados a estimativas produzidas utilizando-se sinais elétricos relativos às descargas.

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.

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.

Estabilizador neural não-linear para sistemas de potência treinado por rede de controladores lineares

A utilização de Estabilizadores de Sistemas de Potência (ESP), para amortecer oscilações eletromecânicas de pequena magnitude e baixa freqüência, é cada vez mais importante na operação dos modernos sistemas elétricos. Estabilizadores convencionais, com estrutura e parâmetros fixos, têm sido utilizados com essa finalidade há algumas décadas, porém existem regiões de operação do sistema nas quais esses estabilizadores lineares não são tão eficientes, especialmente quando comparados com estabilizadores projetados através de modernas técnicas de controle. Um ESP Neural, treinado a partir de um conjunto de controladores lineares locais, é utilizado para investigar em quais regiões de operação do sistema elétrico o desempenho do estabilizador a parâmetros fixos é deteriorada. O melhor desempenho do ESP Neural nessas regiões de operação, quando comparado com o ESP convencional, é demonstrado através de simulações digitais não-lineares de um sistema do tipo máquina síncrona conectada a um barramento infinito e de um sistema com quatro geradores.

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.

Simulação de sistemas de simples estágios de refrigeração por compressão de vapor

Na atualidade, o estudo do desempenho térmico de um sistema de refrigeração por compressão de vapor representa uma ferramenta importante no auxílio do desenvolvimento de novos produtos ou melhoria dos já existentes. Um modelo de simulação em regime permanente foi elaborado para avaliar o desempenho do sistema frigorífico. O sistema estudado inclui uma Central de Ar Condicionado, modelo PA HILTON, constituída de um compressor alternativo do tipo semi-hermético, evaporador e condensador compacto de tubos e aletas e uma válvula de expansão termostática. O modelo do condensador considera três regiões distintas de troca de calor as quais são respectivamente a região de dessuperaquecimento, condensação e subresfriamento. Para a modelagem do evaporador, foram consideradas as regiões de evaporação e superaquecimento. No modelo de simulação foram utilizadas correlações adequadas para a estimativa dos coeficientes de transferência de calor e perda de pressão para cada região do evaporador e condensador. Não foram consideradas a transferência de calor e queda de pressão nas linhas de conexão entre os componentes. A solução do sistema de equações não lineares resultantes da modelagem matemática dos componentes do sistema simulado foi obtida utilizando-se o método das substituições sucessivas com o emprego do software Engineenng Equation Solver . Os resultados obtidos pelo modelo de simulação apresentaram erros inferiores a 9% em relação aos valores experimentais.

Predição dos campos elétricos na faixa de ondas médias em ambiente rural utilizando equações parabólicas

O Brasil está realizando testes para selecionar o padrão de transmissão digital a ser adotado. Sistemas como Digital Radio Mondiale (DRM) e Rádio de Alta Definição (HD Radio), desenvolvido apenas para frequências abaixo de 30 MHz permitem a operação com largura de banda compatível com a utilizada no país, abrindo a possibilidade de coexistência de radiodifusão analógica e digital. Para qualquer sistema a ser adotado são necessários estudos que permitem uma melhor gestão do espectro eletromagnético, que exige conhecimento real do alcance do sinal. A propagação de ondas eletromagnéticas na faixa de ondas médias (MW) é caracterizada pela dependência do campo em relação aos parâmetros elétricos do solo. Para contribuir para o planejamento e adaptação do Plano Básico de Radiodifusão em Onda Média – PMWB, a inclusão de novas estações que operam principalmente em simulcast-canal, torna-se necessário desenvolver ferramentas que permitem a avaliação das características do solo onde não fornece dados precisos, permitindo uma revisão de modelos teóricos para prever a adoção, contribuindo para a implantação do rádio digital em nosso país. Este trabalho apresenta os resultados dos ensaios de campo realizados pela ANATEL (Agência Nacional de Telecomunicações) e Radiobrás para analisar um sistema de DRM (Digital Radio Mondiale) na faixa de ondas médias. Foi gerada uma potência de 50kW, com antenas omni-direcionais operando na frequência de 980kHz e utilizando um veículo de medição para recepção fixa e móvel. Várias vias radiais foram percorridas a partir do transmissor localizado em uma área urbana e rural nos entorno da capital do Brasil, Brasília. A partir desses dados é proposto um modelo para avaliação das características elétricas do solo correspondente ao campo elétrico, através da aplicação do método de Equações Parabólicas e comprovação da eficácia do modelo proposto.

Otimização de antenas lineares e planares para aplicações em sistemas de banda larga

O objetivo deste trabalho é a otimização da largura de banda de antenas linear e planar para aplicações em sistemas de banda larga. Nesse sentido, foi feito um estudo das técnicas de análise, aumento da largura de banda e otimização adequadas para o problema em questão. Como técnica de análise, foi utilizado o método dos momentos, o qual está apresentado no capítulo II. Para aumentar a largura de banda, foram utilizadas as técnicas de colocação de elementos parasitas e construção de fendas no radiador, descritos sucintamente no capítulo III. Como algoritmo de otimização, foi utilizado o algoritmo genético, descrito sucintamente no capítulo II. Neste trabalho, são apresentadas duas propostas de antenas, uma antena dipolo linear combinada com quatros espiras parasitas, capítulo IV, e uma antena planar do tipo espira, capítulo V. No primeiro caso, foram utilizados elementos parasitas e o algoritmo genético para aumentar a largura de banda e, no segundo, foram empregadas fendas no radiador e a otimização paramétrica para este objetivo.

Redução de ordem de modelos matemáticos lineares usando algoritmos genéticos

Nesta dissertação apresenta-se o problema de redução de ordem de modelos dinâmicos lineares, sob o ponto de vista de otimização via Algoritmos Genéticos. Uma função custo, obtida a partir da norma dos coeficientes do numerador da função de transferência do erro entre o modelo original e o reduzido, e minimizada por meio de um algoritmo genético, com consequente calculo dos parâmetros do modelo reduzido. O procedimento e aplicado em alguns exemplos que demonstram a validade da abordagem.