41 resultados para High Definition Finite Difference Time Domain
Resumo:
Redes Ópticas Passivas estão cada vez mais difundidas como rede de acesso de banda larga. Devido à sua importância, tais redes necessitam constantemente de supervisão. Diversas técnicas vêm sendo empregadas para supervisão de redes ópticas, dentre as quais se destacam as que utilizam reflectometria, pois agregam vantagens, tais como: baixo custo, em comparação com outros métodos de supervisão, e não exigência de intervenção na casa do usuário. Esta dissertação apresenta uma estratégia capaz de quantificar, por meio de simulações, o impacto que a técnica de reflectometria no domínio do tempo gera como interferência em transmissões típicas de uma rede WDM-PON.
Resumo:
Este trabalho apresenta um estudo teórico-experimental sobre os efeitos da convecção termossolutal sobre os espaçamentos dendríticos primários de ligas hipoeutéticas Sn-Pb. Assim, um dispositivo experimental de configuração horizontal refrigerado a água foi construído e amostras do referido sistema de ligas foram solidificadas direcionalmente sob condições transientes de extração de calor. Um método teórico-experimental é aplicado para determinar as velocidades de deslocamento da isoterma liquidus, os gradientes térmicos e as taxas de resfriamento. São também apresentados resultados para os respectivos coeficientes de transferência de calor na interface metal/molde os quais foram calculados a partir de uma análise comparativa realizada entre os perfis experimentais de temperatura e valores teóricos fornecidos por um método numérico baseado em volumes finitos. Alguns modelos teóricos para a previsão dos espaçamentos dendríticos primários são comparados com os resultados experimentais obtidos. Finalmente, um estudo comparativo é realizado entre resultados encontrados neste trabalho e aqueles apresentados na literatura para os espaçamentos dendríticos primários das ligas Sn-Pb investigadas quando solidificadas direcionalmente nos sistemas verticais ascendente e descendente, sob as mesmas condições. A análise das microestruturas indica que os espaçamentos dendríticos primários são bastante influenciados pela direção de crescimento do sólido.
Resumo:
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.
Resumo:
O filtro de Kalman é aplicado para filtragem inversa ou problema de deconvolução. Nesta dissertação aplicamos o método de Kalman, considerado como uma outra visão de processamento no domínio do tempo, para separar sinal-ruído em perfil sônico admitido como uma realização de um processo estocástico não estacionário. Em um trabalho futuro estudaremos o problema da deconvolução. A dedução do filtro de Kalman destaca a relação entre o filtro de Kalman e o de Wiener. Estas deduções são baseadas na representação do sistema por variáveis de estado e modelos de processos aleatórios, com a entrada do sistema linear acrescentado com ruído branco. Os resultados ilustrados indicam a aplicabilidade dessa técnica para uma variedade de problemas de processamento de dados geofísicos, por exemplo, ideal para well log. O filtro de Kalman oferece aos geofísicos de exploração informações adicionais para o processamento, problemas de modelamento e a sua solução.
Resumo:
O método eletroresistivo é um dos métodos de prospecção geofísica mais empregados na exploração de água subterrânea. Todavia, ele apresenta a desvantagem de ser de difícil uso em locais habitados. De fato, para realizar sondagens elétricas é necessário, muitas vezes, expandir cabos elétricos até um quilômetro ou mais, dependendo da profundidade que se deseja investigar, o que nem sempre é possível em áreas urbanas. O método transiente eletromagnético no domínio do tempo, é usado também para a exploração de água subterrânea e não apresenta as mesmas desvantagens, porque este método indutivo emprega uma bobina quadrada cujo tamanho não é um impedimento para ser utilizado em locais habitados. Por outro lado, o método eletromagnético é mais sensível a ruídos devido a instalações elétricas. Para comparar as vantagens e desvantagens dos dois métodos em áreas urbanas foram utilizados programas para a modelagem sintética das sondagens elétricas verticais e das sondagens transientes eletromagnéticas, ambos desenvolvidos por RIJO (1994ab, 1996). Para tal propósito foram considerados dois casos: no primeiro, foram definidos quatro tipos de modelos (H, K, A e Q) com base nos horizontes geoelétricos em Palestina do Pará. No segundo caso, se considerou os modelos obtidos pelas 21 sondagens elétricas verticais tipo Schlumberger da mesma área. Com base nestes modelos, se utilizou um programa de transientes electromagnéticas para gerar as curvas sintéticas das resistividades aparentes. Como resultado, os valores de resistividades aparentes dados em função do espaçamento de eletrodos (método eletroresistivo) e em função do tempo (método eletromagnético), permitiu obter nos dois casos mencionados uma relação linear em escala bilogaritmica entre o espaçamento (metros) e o tempo (milisegundos). Deste modo é possível inferir a profundidade de investigação desejada em áreas urbanas usando o método eletromagnético, quando a aplicação do método eletroresistivo ficar limitada pela acessibilidade.
Resumo:
A motivação geológica deste trabalho reside no imageamento de estruturas de bacias sedimentares da região Amazônica, onde a geração e o acúmulo de hidrocarboneto estão relacionados com a presença de soleiras de diabásio. A motivação sísmica reside no fato de que essas rochas intrusivas possuem grandes contrastes de impedância com a rocha encaixante, o que resulta em múltiplas, externas e internas, com amplitudes semelhantes as das primárias. O sinal sísmico das múltiplas podem predominar sobre o sinal das reflexões primárias oriundas de interfaces mais profundas, o que pode dificultar o processamento, a interpretação e o imageamento da seção sísmica temporal. Neste trabalho, estudamos a atenuação de múltiplas em seções sintéticas fonte-comum (FC) através da comparação de dois métodos. O primeiro método resulta da combinação das técnicas Wiener-Hopf-Levinson de predição (WHLP) e o de empilhamento superfície-de-reflexão-comum (CRS), e denominando WHLP-CRS, onde o operador é desenhado exclusivamente no domínio do tempo-espaço. O segundo método utilizado é o filtro de velocidade (ω-k) aplicado após o empilhamento superfície-de-reflexão (CRS), onde o operador é desenhado exclusivamente no domínio bidimensional de freqüência temporal-espacial. A identificação das múltiplas é feita na seção de afastamento-nulo (AN) simulada com o empilhamento CRS, e utiliza o critério da periodicidade entre primária e suas múltiplas. Os atributos da frente de onda, obtidos através do empilhamento CRS, são utilizados na definição de janelas móveis no domínio tempo-espaço, que são usadas para calcular o operador WHLP-CRS. O cálculo do filtroω-k é realizado no domínio da freqüência temporal-espacial, onde os eventos são selecionados para corte ou passagem. O filtro (ω-k) é classificado como filtro de corte, com alteração de amplitude, mas não de fase, e limites práticos são impostos pela amostragem tempo-espaço. Em termos práticos, concluímos que, para o caso de múltiplas, os eventos separados no domínio x-t não necessariamente se separam no domínio ω-k, o que dificulta o desenho de um operador ω-k semelhante em performance ao operador x-t.
Resumo:
O efeito da topografia de um vale, uma colina, um declive e um aclive, em sondagens eletromagnéticas bidimensionais no domínio da freqüência, causa uma variação nos valores da amplitude e da fase da componente Hz devido a um meio homogêneo. A amplitude é menos afetada que a fase. A parede do vale mais próxima da linha de corrente causa uma forte diminuição dos valores da fase, enquanto que a parede do vale mais distante da linha causa um forte aumento. Os efeitos de uma colina são opostos aos do vale. Os efeitos do declive e do aclive num meio homogêneo, são similares, respectivamente, aos observados pelas paredes do vale e da colina mais próximas da linha de corrente. A resposta de um corpo condutivo retangular num meio homogêneo próximo a uma linha de corrente sofre pequenas variações devido a presença de um vale ou de uma colina situada longe da linha de corrente. Porém, se essas feições topográficas estiverem sobre o corpo, elas afetam fortemente a fase e a amplitude da componente Hz e apenas a amplitude, no caso da componente Hx. A resposta transiente não é afetada pela topografia para tempos muito baixos, pois nesse caso se investiga uma finíssima camada da superfície, assim como para tempos muito altos porque a profundidade de investigação é muito grande comparada com a dimensão da topografia. Para os modelos aqui estudados, a maior influência se dá para tempos intermediários, ao redor de 7 ms, ocasionando um retardo do ponto de "cross over" nas curvas de sondagens.
Resumo:
A implementação convencional do método de migração por diferenças finitas 3D, usa a técnica de splitting inline e crossline para melhorar a eficiência computacional deste algoritmo. Esta abordagem torna o algoritmo eficiente computacionalmente, porém cria anisotropia numérica. Esta anisotropia numérica por sua vez, pode levar a falsos posicionamentos de refletores inclinados, especialmente refletores com grandes ângulos de mergulho. Neste trabalho, como objetivo de evitar o surgimento da anisotropia numérica, implementamos o operador de extrapolação do campo de onda para baixo sem usar a técnica splitting inline e crossline no domínio frequência-espaço via método de diferenças finitas implícito, usando a aproximação de Padé complexa. Comparamos a performance do algoritmo iterativo Bi-gradiente conjugado estabilizado (Bi-CGSTAB) com o multifrontal massively parallel solver (MUMPS) para resolver o sistema linear oriundo do método de migração por diferenças finitas. Verifica-se que usando a expansão de Padé complexa ao invés da expansão de Padé real, o algoritmo iterativo Bi-CGSTAB fica mais eficientes computacionalmente, ou seja, a expansão de Padé complexa atua como um precondicionador para este algoritmo iterativo. Como consequência, o algoritmo iterativo Bi-CGSTAB é bem mais eficiente computacionalmente que o MUMPS para resolver o sistema linear quando usado apenas um termo da expansão de Padé complexa. Para aproximações de grandes ângulos, métodos diretos são necessários. Para validar e avaliar as propriedades desses algoritmos de migração, usamos o modelo de sal SEG/EAGE para calcular a sua resposta ao impulso.
Resumo:
A modelagem do mCSEM é feita normalmente no domínio da frequência, desde sua formulação teórica até a análise dos resultados, devido às simplificações nas equações de Maxwell, possibilitadas quando trabalhamos em um regime de baixa frequência. No entanto, a abordagem através do domínio do tempo pode em princípio fornecer informação equivalente sobre a geofísica da subsuperfície aos dados no domínio da frequência. Neste trabalho, modelamos o mCSEM no domínio da frequência em modelos unidimensionais, e usamos a transformada discreta de Fourier para obter os dados no domínio do tempo. Simulamos ambientes geológicos marinhos com e sem uma camada resistiva, que representa um reservatório de hidrocarbonetos. Verificamos que os dados no domínio do tempo apresentam diferenças quando calculados para os modelos com e sem hidrocarbonetos em praticamente todas as configurações de modelo. Calculamos os resultados considerando variações na profundidade do mar, na posição dos receptores e na resistividade da camada de hidrocarbonetos. Observamos a influência da airwave, presente mesmo em profundidades oceânicas com mais de 1000m, e apesar de não ser possível uma simples separação dessa influência nos dados, o domínio do tempo nos permitiu fazer uma análise de seus efeitos sobre o levantamento. Como parte da preparação para a modelagem em ambientes 2D e 3D, fazemos também um estudo sobre o ganho de desempenho pelo uso do paralelismo computacional em nossa tarefa.
Resumo:
O reconhecimento da intenção do condutor a partir de sinais de eletroencefalografia (EEG) pode ser útil no desenvolvimento de interfaces cérebro computador (BCI) para serem usadas em sinergia com veículos inteligentes. Isso pode ser benéfico para melhorar a qualidade de interação entre o motorista e o carro, por exemplo, fornecendo uma resposta do carro inteligente alinhada com a intenção do motorista. Neste estudo, considera-se a antecipação como sendo o estado cognitivo que leva a ações especificas durante a condução de um automóvel. Portanto, propomos investigar a presença de padrões antecipatórios em sinais EEG durante a condução de veículos para determinar duas ações especifícas (1) virar à esquerda e (2) virar à direita, alguns milissegundos antes que tais ações aconteçam. Um protocolo experimental foi proposto para gravar sinais EEG de 5 indivíduos enquanto eles operam um simulador de realidade virtual não invasiva - que foi projetado para tal experimento - que simula a condução de um carro virtual. O protocolo experimental é uma variante do paradigma da variação negativa contingente (CNV) com condições Go e No-go no sistema de condução de realidade virtual. Os resultados apresentados neste estudo indicam a presença de padrões antecipatórios em potenciais corticais lentos observados no domínio do tempo (medias dos sinais EEG) e da frequência (Power Spectra e coerência de fase). Isso abre um leque de possibilidades no desenvolvimento de sistemas BCI - baseados em sinais antecipatórios - que conectem o motorista ao veiculo inteligente favorecendo uma tomada de decisão que analise as intenções dos condutores podendo eventualmente evitar acidentes durante a condução.
Resumo:
The main purpose of this paper is to investigate both the columnar to equiaxed transition and primary dendritic arm spacings of Al-3wt.%Si alloy during the horizontal directional solidification. The transient heat transfer coefficient at the metal-mold interface is calculated based on comparisons between the experimental thermal profiles in castings and the simulations provided by a finite difference heat flow program. Simulated curve of the interfacial heat transfer coefficient was used in another numerical solidification model to determine theoretical values of tip growth rates, cooling rates and thermal gradients that are associated with both columnar to equiaxed transition and primary dendritic arm spacings. A good agreement was observed between the experimental values of these thermal variables and those numerically simulated for the alloy examined. A comparative analysis is carried out between the experimental data of this work and theoretical models from the literature that have been proposed to predict the primary dendritic spacings. In this context, this study may contribute to the understanding of how to manage solidification operational parameters aiming at designing the microstructure of Al-Si alloys.