Modelagem numérica do terreno na avaliação do volume admissível em canais de cordões de contorno nivelados

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

Simulação da produtividade de cana-de-açúcar nas regiões nordeste e sudeste do Brasil, utilizando o modelo DSSAT

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

Simulação de transporte eletrônico em dispositivos unimoleculares baseados em indicadores de pH

Nas últimas décadas, diversos pesquisadores têm tentado empregar moléculas em dispositivos eletrônicos de nanoescala. Por este motivo, diferentes parâmetros eletro/ópticos, que regem o transporte eletrônico em moléculas orgânicas, precisam ser analisados. Neste trabalho foi desenvolvido um estudo de transporte de carga para o composto Vermelho de Propila, popularmente utilizado como indicador de pH. A motivação para estudá-lo resulta de sua estrutura constituída por subunidades doadora-aceitadora, acopladas via grupo azo (N=N), uma característica bem conhecida em retificadores moleculares. A metodologia utilizada para tratar o sistema em equilíbrio é baseada em métodos de Mecânica Molecular e Hartree-Fock. Sendo que, para simular o sistema em não-equilíbrio, foi empregado o formalismo de Landauer-Büttiker. Através desses métodos, as curvas características do sistema molecular foram traçadas e comparadas. O resultado da comparação permitiu explicar os fenômenos que regem o transporte eletrônico na nanoestrutura. Além disso, foram analisados os efeitos de contatos metálicos, ligados a molécula na presença de campo elétrico externo.

Simulação espacial da perda de solo e nutrientes por erosão em um argissolo sob cultivo de cana-de-açúcar

Pós-graduação em Agronomia (Ciência do Solo) - FCAV

Metodologia de discretização espacial para o método radial de interpolação por pontos (RPIM) aplicada para solução numérica das equações de Maxwell

Neste trabalho, foi desenvolvido e implementado um método de discretização espacial baseado na lei de Coulomb para geração de pontos que possam ser usados em métodos meshless para solução das equações de Maxwell. Tal método aplica a lei de Coulomb para gerar o equilíbrio espacial necessário para gerar alta qualidade de discretização espacial para um domínio de análise. Este método é denominado aqui de CLDM (Coulomb Law Discretization Method ) e é aplicado a problemas bidimensionais. Utiliza-se o método RPIM (Radial Point Interpolation Method) com truncagem por UPML (Uniaxial Perfectlly Matched Layers) para solução das equações de Maxwell no domínio do tempo (modo TMz).

Desenvolvimento de framework para análise e simulação dinâmica de sistemas elétricos de potência

Esta tese apresenta uma metodologia flexível orientada a objetos (OO) para a aplicação no projeto e implementação de sistemas de software utilizados na realização de estudos dinâmicos de sistemas elétricos de grande porte. A metodologia OO proposta objetiva tornar mais simples o desenvolvimento, a atualização e a manutenção de complexos sistemas de software para estudos de transitórios eletromecânicos em sistemas elétricos de potência. Os requisitos de usuário são mapeados para um conjunto de classes básicas, as quais são usadas para efetuar a modelagem de dispositivos dinâmicos tais como geradores elétricos. Para avaliação da metodologia foram realizados dois estudos de casos. No primeiro estudo caso o Framework foi aplicado na simulação das unidades geradoras da Usina Hidrelétrica de Tucuruí. Os resultados da simulação foram comparados com medições obtidos em ensaios no campo e mostrou a boa performance do Framework na reprodução dos fenômenos eletromecânicos desta usina de grande porte. No segundo estudo de caso, por outro lado, o Framework foi aplicado na modelagem de um sistema de geração fotovoltaico (PV) com seu sistema de Rastreamento da Potência Máxima (MPPT). O controle MPPT foi implementado usando técnicas digitais. Os resultados das simulações demonstram a performance do Framework na modelagem do sistema de controle de corrente, assim como no controle MPPT, dos sistemas de geração PV.

Simulação computacional no auxílio à tomada de decisão nos processos industriais: utilizando a ferramenta Tecnomatix Plant Simulation 9.0

Atualmente existem inúmeros processos produtivos automatizados, os quais vêm se tornando cada vez mais complexos em função das necessidades do mundo moderno e, portanto, demandam nas fases de projetos e de implementação ferramentas de engenharia cada vez mais poderosas para modelá-los e analisá-los da maneira mais eficiente possível. Nesse ambiente de crescente pressão por resultados positivos, racionalização e aprimoramento de recursos internos, que a ferramenta computacional Tecnomatix Plant Simulation 9.0, surge como um caminho para a obtenção de competitividade produtiva. Ressaltando ainda que, o estudo proposto é de grande relevância para os profissionais da gestão produtiva, os quais almejam resultados que minimizem custos e maximizem lucros.

Controladores fuzzy aplicados ao conversor de geradores de indução duplamente excitados em sistemas eólicos integrados a sistemas de potência

Este artigo apresenta novas estratégias de controle fuzzy aplicadas ao conversor estático interligado ao rotor dos geradores de indução duplamente excitados (DFIG), em esquemas eólicos de velocidade variável, integrados a redes elétricas reais. As estratégias de controle propostas são do tipo fuzzy "look-up-table" com supervisão. O desempenho de tais reguladores, baseados em inteligência computacional, é comparado com o obtido com controladores PI's a parâmetros fixos, na ocorrência de faltas no sistema elétrico de potência. Esses controladores fuzzy são destinados a garantir a continuidade da operação dos conversores, e melhorar o desempenho transitório do sistema, em relação aos controladores convencionais. Os resultados apresentados confirmam a eficácia das estratégias de controle adotadas. O modelo físico dos aerogeradores, consistindo de um grande número de turbinas eólicas, foi implementado através de uma ferramenta de simulação dinâmica, desenvolvida no ambiente computacional MATLABTM.

Identificação de correntes de inrush na proteção diferencial de transformadores de potência através do gradiente da corrente diferencial e de mapas auto-organizáveis

A principal dificuldade encontrada na proteção diferencial de transformadores de potência é a correta distinção entre as correntes de inrush e as correntes de faltas internas. Tradicionalmente os relés diferenciais executam esta tarefa utilizando a técnica de restrição por harmônicos baseada na premissa de que as correntes de inrush possuem alta concentração de componentes harmônicas de segunda ordem, contudo essa técnica nem sempre é eficaz. O presente trabalho tem como objetivo apresentar a proposta de duas novas metodologias capazes de realizar a identificação e distinção entre as correntes de inrush das correntes de faltas internas na proteção diferencial de transformadores de potência através de metodologias que não dependem do conteúdo de harmônicos do sinal da corrente diferencial. A primeira metodologia proposta, denominada de método do gradiente da corrente diferencial, é baseada no comportamento do vetor gradiente, obtido através da diferenciação numérica do sinal da corrente diferencial. O critério de distinção utilizado é baseado no desvio padrão do ângulo do vetor gradiente que apresenta comportamento diferenciado para correntes de inrush e correntes de curto-circuito. A segunda metodologia proposta é baseada na capacidade de reconhecimento e classificação de padrões das redes neurais de Mapeamento Auto-organizável de Kohonen. Como padrão de entrada e de treinamento da rede neural é utilizado um vetor contendo quatro níveis do espectro do desvio padrão do ângulo do vetor gradiente da corrente diferencial nas três fases do transformador de potência. A eficácia dos métodos foi testada através da simulação de diversas situações de faltas internas e correntes de inrush, incluindo situações de “Sympathetic Inrush”, em um transformador de potência usando o software EMTP/ATP e através da implementação do algoritmo em MATLAB®, apresentando resultados altamente promissores.

Modelos equivalentes de parques eólicos usando algoritmos genéticos

Esta dissertação apresenta uma metodologia baseada em algoritmo genético (AG) para determinar modelos dinâmicos equivalentes de parques eólicos com geradores de indução em gaiola de esquilo ( GIGE) e geradores de indução duplamente alimentados ( GIDA), apresentando parâmetros elétricos e mecânicos distintos. A técnica se baseia em uma formulação multiobjetiva solucionada por um AG para minimizar os erros quadráticos das potências ativa e reativa entre modelo de um único gerador equivalente e o modelo do parque eólico investigado. A influência do modelo equivalente do parque eólico no comportamento dinâmico dos geradores síncronos é também investigada por meio do método proposto. A abordagem é testada em um parque eólico de 10MW composto por quatro turbinas eólicas ( 2x2MW e 2x3MW), consistindo alternadamente de geradores GIGE e GIDA interligados a uma barra infinita e posteriormente a rede elétrica do IEEE 14 barras. Os resultados obtidos pelo uso do modelo dinâmico detalhado para a representação do parque eólico são comparados aos do modelo equivalente proposto para avaliar a precisão e o custo computacional do modelo proposto.

Análise da influência do monitoramento OTDR em WDM-PON: um estudo baseado em simulação

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.

Análise de sistemas de controle de vibração em máquinas rotativas utilizando atuadores formados por ligas com memória de forma

A aplicação das ligas com memória de forma (shape memory alloys – SMA) têm se mostrado como uma alternativa promissora no controle de vibração de máquinas e estruturas, devido principalmente aos fenômenos de memória de forma e pseudoelástico que elas apresentam. Do mesmo modo, tais ligas proporcionam grandes forças de recuperação e capacidade de amortecimento quando comparadas aos materiais tradicionais. Na literatura científica encontra-se um grande número de trabalhos que tratam da aplicação das SMA no controle de vibração em estruturas. Contudo, a aplicação desse tipo de material em máquinas rotativas ainda é um assunto pouco abordado. Nesse sentido, busca-se explorar numericamente o comportamento de atuadores baseados em ligas com memória de forma para o controle de vibração em máquinas rotativas. Na primeira análise deste trabalho um rotor tipo Jeffcott com luvas SMA em um dos mancais é utilizado. São empregadas diferentes espessuras de luvas nos estados martensítico e austenítico e as variações em termos de amplitude e frequência são então comparadas. Posteriormente, dois diferentes sistemas rotativos com dois discos e molas SMA aplicadas em um e dois mancais são estudados sob configurações variadas. As molas foram posicionadas externamente aos mancais e a temperatura de operação desses componentes é ajustada de acordo com a necessidade do controle de vibração. Além disso, foi utilizado um código computacional para a representação do comportamento termomecânico de molas SMA assim como um programa baseado no Método de Elementos Finitos (MEF) para a simulação do comportamento dinâmico dos rotores. Os resultados das análises numéricas demonstram que as SMA são eficientes no controle de vibração de sistemas rotativos devido obterem-se reduções significativas das amplitudes de deslocamento, modificações nas velocidades críticas, supressão de movimentos indesejáveis e controle das órbitas de precessão.

Implantação de CGH's na Amazônia: uma visão global com ênfase nos aspectos hidrológicos, tecnológicos, ambientais e financeiros

Um problema de grande relevância social na Amazônia é o fato das pequenas comunidades isoladas não serem contempladas pelos benefícios dos grandes empreendimentos hidrelétricos instalados na região. Uma solução alternativa a para este problema, é o aproveitamento da grande malha de pequenos rios e igarapés da região a partir da implantação de Centrais Geradoras Hidrelétricas – CGH’s (antigas mini e microcentrais).. Dentro deste contexto, o presente trabalho analisa e discute certos aspectos gerais inerentes à implantação das CGH’s, e principalmente, trazendo-os à realidade das pequenas bacias Amazônicas. Os aspectos a serem abordados neste trabalho dizem respeito às avaliações preliminares de terreno, aspectos hidrológicos, tecnológicos, ambientais e financeiros. Os aspectos são analisados de forma global para implantação de CGH’s e suas etapas, também analisados e aplicados a um estudo de caso – implantação da CGH irmã Dorothy na pequena bacia hidrográfica do Igarapé são João em Anapú-pa, onde foi possível a utilização de metodologia existente para um estudo aprimorado de implantação de CGH’s na Amazônia. No âmbito hidrológico, utilizou-se um modelo chuva-vazão desenvolvido por Blanco 2005 e aplicou a pequenas bacias da Amazônia que não possuem registros de vazão. Nos aspectos tecnológicos utilizou ferramentas computacionais para predição de desempenho de turbinas axiais de baixa queda adaptas ao relevo da região, e simulados para a turbina axial a ser implantada na CGH irmã Dorothy. No contexto ambiental atualmente há forte cobrança pelas autoridades competentes locais para realização do estudo ambiental inerente ao aproveitamento, e se tratando da região Amazônica, ambientalmente muita agredida pela ação do homem, os estudos devem ser bem definidos, apontando os possíveis impactos que podem ser causados pela CGH, descritos no RAS (Relatório Ambiental Simplificado) anexo deste trabalho. Desta forma foi desenvolvido um método para cálculo e simulação da área inundada para implantação de CGH’s na Amazônia, aplicado à CGH irmã Dorothy. No âmbito financeiro, são raras as informações referentes aos custos de implantação de CGH’s na Amazônia. Assim, foram levantados os custos referentes à implantação da CGH irmã Dorothy e comparados a custos de CGH’s e de geradores a diesel disponíveis na literatura.

Modelagem numérica de dados MCSEM 3D usando computação paralela

Desenvolvemos a modelagem numérica de dados sintéticos Marine Controlled Source Electromagnetic (MCSEM) usada na exploração de hidrocarbonetos para simples modelos tridimensionais usando computação paralela. Os modelos são constituidos de duas camadas estrati cadas: o mar e o sedimentos encaixantes de um delgado reservatório tridimensional, sobrepostas pelo semi-espaço correspondente ao ar. Neste Trabalho apresentamos uma abordagem tridimensional da técnica dos elementos nitos aplicada ao método MCSEM, usando a formulação da decomposição primária e secundária dos potenciais acoplados magnético e elétrico. Num pós-processamento, os campos eletromagnéticos são calculados a partir dos potenciais espalhados via diferenciação numérica. Exploramos o paralelismo dos dados MCSEM 3D em um levantamento multitransmissor, em que para cada posição do transmissor temos o mesmo processo de cálculos com dados diferentes. Para isso, usamos a biblioteca Message Passing Interface (MPI) e o modelo servidor cliente, onde o processador administrador envia os dados de entradas para os processadores clientes computar a modelagem. Os dados de entrada são formados pelos parâmetros da malha de elementos nitos, dos transmissores e do modelo geoelétrico do reservatório. Esse possui geometria prismática que representa lentes de reservatórios de hidrocarbonetos em águas profundas. Observamos que quando a largura e o comprimento horizontais desses reservatório têm a mesma ordem de grandeza, as resposta in-line são muito semelhantes e conseqüentemente o efeito tridimensional não é detectado. Por sua vez, quando a diferença nos tamanhos da largura e do comprimento do reservatório é signi cativa o efeito 3D é facilmente detectado em medidas in-line na maior dimensão horizontal do reservatório. Para medidas na menor dimensão esse efeito não é detectável, pois, nesse caso o modelo 3D se aproxima de um modelo bidimensional. O paralelismo dos dados é de rápida implementação e processamento. O tempo de execução para a modelagem multitransmissor em ambiente paralelo é equivalente ao tempo de processamento da modelagem para um único transmissor em uma máquina seqüêncial, com o acréscimo do tempo de latência na transmissão de dados entre os nós do cluster, o que justi ca o uso desta metodologia na modelagem e interpretação de dados MCSEM. Devido a reduzida memória (2 Gbytes) em cada processador do cluster do departamento de geofísica da UFPA, apenas modelos muito simples foram executados.

Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais

A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10^-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10^-1 Hz e para freqüências menores que 10^-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.