Deconvolução de Euler: passado, presente e futuro – um tutorial

Neste tutorial apresentamos uma revisão da deconvolução de Euler que consiste de três partes. Na primeira parte, recordamos o papel da clássica formulação da deconvolução de Euler 2D e 3D como um método para localizar automaticamente fontes de campos potenciais anômalas e apontamos as dificuldades desta formulação: a presença de uma indesejável nuvem de soluções, o critério empírico usado para determinar o índice estrutural (um parâmetro relacionado com a natureza da fonte anômala), a exeqüibilidade da aplicação da deconvolução de Euler a levantamentos magnéticos terrestres, e a determinação do mergulho e do contraste de susceptibilidade magnética de contatos geológicos (ou o produto do contraste de susceptibilidade e a espessura quando aplicado a dique fino). Na segunda parte, apresentamos as recentes melhorias objetivando minimizar algumas dificuldades apresentadas na primeira parte deste tutorial. Entre estas melhorias incluem-se: i) a seleção das soluções essencialmente associadas com observações apresentando alta razão sinal-ruído; ii) o uso da correlação entre a estimativa do nível de base da anomalia e a própria anomalia observada ou a combinação da deconvolução de Euler com o sinal analítico para determinação do índice estrutural; iii) a combinação dos resultados de (i) e (ii), permitindo estimar o índice estrutural independentemente do número de soluções; desta forma, um menor número de observações (tal como em levantamentos terrestres) pode ser usado; iv) a introdução de equações adicionais independentes da equação de Euler que permitem estimar o mergulho e o contraste de susceptibilidade das fontes magnéticas 2D. Na terceira parte apresentaremos um prognóstico sobre futuros desenvolvimentos a curto e médio prazo envolvendo a deconvolução de Euler. As principais perspectivas são: i) novos ataques aos problemas selecionados na segunda parte deste tutorial; ii) desenvolvimento de métodos que permitam considerar interferências de fontes localizadas ao lado ou acima da fonte principal, e iii) uso das estimativas de localização da fonte anômala produzidas pela deconvolução de Euler como vínculos em métodos de inversão para obter a delineação das fontes em um ambiente computacional amigável.

A Law of the Wall Formulation for Recirculating Flows

ABSTRACT: This work presents a new law of the wall formulation for recirculating turbulent flows. An alternative expression for the internal length which can be applied in the separated region is also presented. The formulation is implemented in a numerical code which solves the k-e model through a finite volume method. The theoretical results are compared with the experimental data of Vogel and Eaton (J. of Heat Transfer, Transactions of ASME, vol.107, pp. 922-929, 1985). The paper shows that the present formulation furnishes better results than the standard k-e formulation.

Projeto e construção de câmaras reverberantes em escala reduzida para o estudo das características de perda de transmissão de divisórias confeccionadas a partir de materiais regionais

Projetar e construir câmaras reverberantes, em modelo reduzido, para a determinação da perda de transmissão de materiais usados para o controle de ruído nas áreas industriais e da construção civil foi o foco principal deste trabalho. Por outro lado, buscou-se estabelecer uma metodologia de análise e interpretação dos parâmetros acústicos a partir da teoria dos modelos em escala reduzida, de tal modo que se possa inferir a perda de transmissão dos materiais estudados. A necessidade desta construção ocorre devido aos estudos dos desempenhos acústicos dos materiais serem realizados em câmaras reverberantes (em escala real), disponíveis apenas em algumas regiões do Brasil, dificultando e encarecendo os testes de desempenho dos materiais para pesquisadores de outras regiões. Assim, de modo a se ter um custo baixo, foi realizada a construção das câmaras reverberantes, em escala reduzida de 1:6, sendo o objetivo principal o estudo da perda de transmissão de divisórias desenvolvidas com material regional, estabelecendo comparação entre os seus desempenhos, e ao mesmo tempo permitir a sua utilização em aulas práticas de acústica. Apesar das câmaras terem sido construídas em escala reduzida, o fator de escala (1:6) não foi considerado e as medições foram realizadas nas freqüências reais das bandas de oitava delimitadas entre 125 Hz e 8 kHz. As divisórias foram ensaiadas para a determinação dos parâmetros acústicos, tendo por base as normas ISO 140 e 354. Finalmente, a partir dos resultados encontrados, foi possível comprovar a eficiência dos materiais utilizados para a obtenção de divisórias que podem ser usadas para o isolamento de ruído em ambientes construídos e/ou enclausuramentos industriais. Além dos ensaios experimentais, foram realizados cálculos analíticos segundo as equações da Análise Estatística Energética (SEA) para efeito de comprovação e validação dos resultados.

Formulação e implementação de um pacote para interpretar mapas magnéticos utilizando os polinômios de Walsh

O presente trabalho consiste na formulação de uma metodologia para interpretação automática de dados de campo magnético. Desta forma, a sua utilização tornará possível a determinação das fronteiras e magnetização de cada corpo. Na base desta metodologia foram utilizadas as características de variações abruptas de magnetização dos corpos. Estas variações laterais abruptas serão representadas por polinômios descontínuos conhecidos como polinômios de Walsh. Neste trabalho, muitos conceitos novos foram desenvolvidos na aplicação dos polinômios de Walsh para resolver problemas de inversão de dados aeromagnéticos. Dentre os novos aspectos considerados, podemos citar. (i) - O desenvolvimento de um algoritmo ótimo para gerar um jôgo dos polinômios "quase-ortogonais" baseados na distribuição de magnetização de Walsh. (ii) - O uso da metodologia damped least squares para estabilizar a solução inversa. (iii) - Uma investigação dos problemas da não-invariância, inerentes quando se usa os polinômios de Walsh. (iv) - Uma investigação da escolha da ordem dos polinômios, tomando-se em conta as limitações de resolução e o comportamento dos autovalores. Utilizando estas características dos corpos magnetizados é possível formular o problema direto, ou seja, a magnetização dos corpos obedece a distribuição de Walsh. É também possível formular o problema inverso, na qual a magnetização geradora do campo observado obedece a série de Walsh. Antes da utilização do método é necessária uma primeira estimativa da localização das fontes magnéticas. Foi escolhida uma metodologia desenvolvida por LOURES (1991), que tem como base a equação homogênea de Euler e cujas exigências necessárias à sua utilização é o conhecimento do campo magnético e suas derivadas. Para testar a metodologia com dados reais foi escolhida uma região localizada na bacia sedimentar do Alto Amazonas. Os dados foram obtidos a partir do levantamento aeromagnético realizado pela PETROBRÁS.

Interpretação aeromagnética automática com uso da equação homogênea de Euler e sua aplicação na Bacia do Solimões

ABSTRACT: We present here a methodology for the rapid interpretation of aeromagnetic data in three dimensions. An estimation of the x, y and z coordinates of prismatic elements is obtained through the application of "Euler's Homogeneous equation" to the data. In this application, it is necessary to have only the total magnetic field and its derivatives. These components can be measured or calculated from the total field data. In the use of Euler's Homogeneous equation, the structural index, the coordinates of the corners of the prism and the depth to the top of the prism are unknown vectors. Inversion of the data by classical least-squares methods renders the problem ill-conditioned. However, the inverse problem can be stabilized by the introduction of both a priori information within the parameter vector together with a weighting matrix. The algorithm was tested with synthetic and real data in a low magnetic latitude region and the results were satisfactory. The applicability of the theorem and its ambiguity caused by the lack of information about the direction of total magnetization, inherent in all automatic methods, is also discussed. As an application, an area within the Solimões basin was chosen to test the method. Since 1977, the Solimões basin has become a center of exploration activity, motivated by the first discovery of gas bearing sandstones within the Monte Alegre formation. Since then, seismic investigations and drilling have been carried on in the region. A knowledge of basement structures is of great importance in the location of oil traps and understanding the tectonic history of this region. Through the application of this method a preliminary estimate of the areal distribution and depth of interbasement and sedimentary magnetic sources was obtained.

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.