Um metodo numerico com paralelismo no tempo para aproximar solucoes de EDPs

Este trabalho de pesquisa tem por objetivo apresentar e investigar a viabilidade de um método numérico que contempla o paralelismo no tempo. Este método numérico está associado a problemas de condição inicial e de contorno para equações diferenciais parciais (evolutivas). Diferentemente do método proposto neste trabalho, a maioria dos métodos numéricos associados a equações diferencias parciais evolutivas e tradicionalmente encontrados, contemplam apenas o paralelismo no espaço. Daí, a motivação em realizar o presente trabalho de pesquisa, buscando não somente um método com paralelismo no tempo mas, sobretudo, um método viável do ponto de vista computacional. Para isso, a implementação do esquema numérico proposto está por conta de um algoritmo paralelo escrito na linguagem C e que utiliza a biblioteca MPI. A análise dos resultados obtidos com os testes de desempenho revelam um método numérico escalável e que exige pouco nível de comunicação entre processadores.

Previsão do comportamento hidráulico de um rio com base na estimativa de coeficientes que controlam seu escoamento. Estudo de caso: Rio Bengalas, Nova Friburgo-RJ

Nesta tese é realizada a modelagem do comportamento hidráulico dos principais rios que compõem a bacia hidrográfica do Rio Bengalas, localizada no município de Nova Friburgo-RJ, a qual abrange a área mais urbanizada da referida cidade. Para a realização das simulações foi utilizado o Sistema de Modelagem de Águas MOHID, ferramenta MOHID Land. Já para a calibração do modelo foram adotados alguns métodos de otimização, mais precisamente, os algoritmos de Luus- Jaakola (LJ) e Colisão de Partículas (PCA), acoplados ao referido sistema, com o intuito de determinar os principais parâmetros necessários à modelagem de corpos hídricos, bem como suas bacias hidrográficas. Foram utilizados dados topográficos do IBGE disponibilizados pela prefeitura após a elaboração do Plano de Águas Pluviais da região de interesse. Com o modelo devidamente calibrado por meio de dados experimentais, foi realizada a validação do mesmo através da simulação de inundações nesta região. Apesar de técnicas de otimização acopladas à plataforma MOHID terem sido utilizadas pela primeira vez em um rio de montanha, os resultados apresentaram-se importantes e qualitativamente satisfatórios do ponto de vista de auxílio à tomada de decisões, tendo como base a prevenção de danos causados pelas elevações da lâmina dágua que ocorrem frequentemente em Nova Friburgo, como por exemplo, a recente tragédia de janeiro de 2011 ocorrida na Região Serrana do Estado do Rio de Janeiro.

Análise numérica do problema de difusão anômala unidimensional

A presente dissertação tem como objetivo analisar o comportamento da solução numérica da equação de difusão anômala com distribuição de fluxo bimodal, no regime estacionário, através de dois métodos numéricos. Foram desenvolvidos modelos utilizando o Método de Elementos Finitos e o Método de Volumes Finitos para a solução numérica desta equação. No modelo do Método de Elementos Finitos utilizou-se polinômios cúbicos de Hermite como funções de interpolação. No modelo de Volumes Finitos foi utilizada uma discretização de ordem superior para a avaliação das derivadas da equação em estudo. Em ambos os métodos, os modelos desenvolvidos consideram a utilização de diferentes tipos de condições de contorno para a solução do problema. Foram analisadas as influências de parâmetros da equação, das condições de contorno e do refinamento da malha na solução numérica. Os resultados apresentam a análise de erros da solução numérica através da comparação desta com a solução analítica.

Um modelo físico-matemático para escoamentos em meios porosos com transição insaturado-saturado.

Neste trabalho é apresentada uma nova modelagem matemática para a descrição do escoamento de um líquido incompressível através de um meio poroso rígido homogêneo e isotrópico, a partir do ponto de vista da Teoria Contínua de Misturas. O fenômeno é tratado como o movimento de uma mistura composta por três constituintes contínuos: o primeiro representando a matriz porosa, o segundo representando o líquido e o terceiro representando um gás de baixíssima densidade. O modelo proposto possibilita uma descrição matemática realista do fenômeno de transição insaturado/saturado a partir de uma combinação entre um sistema de equações diferenciais parciais e uma desigualdade. A desigualdade representa uma limitação geométrica oriunda da incompressibilidade do líquido e da rigidez do meio poroso. Alguns casos particulares são simulados e os resultados comparados com resultados clássicos, mostrando as consequências de não levar em conta as restrições inerentes ao problema.

Redes neurais artificiais aplicadas à determinação do tamanho ótimo da malha para o cálculo da intensidade útil

Neste trabalho é apresentado um estudo para a determinação do tamanho ótimo da malha de elementos, utilizando redes neurais artificiais, para o cálculo da intensidade útil. A ideia principal é treinar as redes de modo a possibilitar a aprendizagem e o reconhecimento do melhor tamanho para diversas áreas superficiais em fontes sonoras com geometria plana. A vantagem de se utilizar redes neurais artificiais deve-se ao fato de apresentarem um único tamanho para a obtenção da intensidade útil, consequentemente, uma redução significativa de tempo computacional quando comparado com o tempo de cálculo de uma malha bem refinada. Ensaios numéricos com placas planas - geometria separável que permite uma solução analítica - são utilizados para se realizar comparações. É apresentado um estudo comparativo entre o tempo computacional gasto para a obtenção da intensidade útil e o mesmo com a malha otimizada via redes neurais artificiais. Também é apresentada uma comparação do nível de potência sonora mediante solução numérica, a fim de validar os resultados apresentados pelas redes neurais.

Uma nova descrição para a transferência de massa em meios porosos com transição saturado-insaturado.

Esse texto trata do problema de um fluido contaminado escoando por um meio poroso, tratando os componentes na mistura como meios contínuos. Na primeira parte, desenvolvemos a teoria de misturas de meios contínuos e discutimos equações da continuidade, momento linear e momento angular. A seguir, descrevemos o problema em detalhe e fazemos hipóteses para simplificar o escoamento. Aplicamos as equações encontradas anteriormente para encontrarmos um sistema de equações diferenciais parciais. Desse ponto em diante, o problema se torna quase puramente matemático. Discutimos o caso insaturado, e depois a saturação do meio poroso. Finalmente, adicionamos um contaminante à mistura e, em seguida, N contaminantes.

Construção de quadrados mágicos pelo método do passo uniforme

Lehmer (1929) analisa matematicamente o método do passo uniforme para construção de quadrados mágicos de ordem impar. Ele divide sua análise em várias etapas. Na primeira delas, envolvendo a discussão de condições necessárias e suficientes para o preenchimento do quadrado pelo método, o autor afirma que se dois números guardarem entre si uma certa relação, eles serão designados a ocupar a mesma célula do quadrado causando seu não preenchimento. A análise do preenchimento pelo método do passo uniforme envolve a resolução de um sistema linear módulo n. Nesse trabalho, discutimos o comportamento das soluções desse sistema quando o método falha no preenchimento. Como consequência, concluímos que números que guardam a relação mencionada nunca ocupam a mesma célula. A análise das condições necessárias e suficientes para obter quadrados mágicos segundo a definição de Lehmer (1929) envolve a resolução de equações de congruências lineares a duas variáveis. Nesse trabalho, detalhamos os resultados de Lehmer (1929). A análise das condições necessárias e suficientes para obtenção de quadrados mágicos, como são reconhecidos usualmente, também envolve a resolução de equações de congruências lineares a duas variáveis. Discutimos o comportamento das soluções dessas equações para obter diagonais principais mágicas. Como consequência, mostramos que diagonais principais mágicas são obtidas se e somente se as coordenadas iniciais guardarem certas relações

Localização em poço quântico de Ga_(1-x)Al_(x)AsGaAs com campo magnético inclinado intenso

Nesta dissertação, apresentamos resultados da análise de um gás de elétrons bidimensional, confinado em um poço quântico parabólico de Ga_(1-x)Al_(x)As com campo magnético aplicado inclinado e intenso. Mostramos uma solução analítica exata para o caso de um gás de elétrons não interagente e provamos que o fator de Landé efetivo g*, que neste caso varia com a posição devido à própria variação da concentração de Al na liga GaAlAs, também contribui com um termo parabólico nas frequências normais, tornando-as dependente do spin. O termo normal do splitting Zeeman aparece, no entanto, com um fator g dado por aquele do GaAs. Indo além desse tratamento, incluímos um campo elétrico externo, aplicado paralelamente a direção de crescimento, de modo a modular a distribuição de cargas dentro do poço, e com isso controlar externamente os efeitos dos cruzamentos dos níveis. Calculamos numericamente os níveis de energia no poço parabólico quântico, variando o campo magnético aplicado e a inclinação do campo.

Estudo da precisão e da exatidão da técnica da birrefringência acústica na avaliação de tensões residuais e aplicadas.

Este trabalho apresenta um procedimento numérico para o estudo da precisão e da exatidão da técnica de birrefringência acústica como usada no Instituto de Engenharia Nuclear (IEN) para avaliação de tensões residuais e aplicadas em estruturas metálicas. Esse procedimento deverá ser incorporado posteriormente ao módulo de processamento de sinal do sistema ultrassônico utilizado no Laboratório de Ultrassom do IEN para levar em conta, de forma automática e sistemática, as incertezas nos dados de entrada e as suas propagações ao longo dos cálculos efetuados. A birrefringência acústica é geralmente definida a partir das velocidades de duas ondas ultrassônicas volumétricas de incidência normal e ortogonais entre si. A birrefringência pode ser definida diretamente a partir dos tempos de percurso dessas duas ondas, uma vez que elas percorrem o mesmo espaço físico. Os tempos de percurso das ondas podem ser assim considerados como as variáveis primárias de interesse. Por meio da teoria da acustoelasticidade é possível relacionar a birrefringência acústica com as tensões atuantes no material explorando o fato da velocidade da onda ultrassônica ser afetada pela presença de um campo de tensões. Nesta dissertação elaborou-se um conjunto de planilhas eletrônicas no aplicativo Microsoft Excel para efetuar de forma automática todos os cálculos necessários ao objetivo proposto, levando em conta as incertezas nos dados, sua propagação ao longo dos cálculos e o número de dígitos significativos nos resultados. Como exemplo do procedimento desenvolvido, o estudo da precisão e da exatidão na determinação das tensões em uma viga de aço sob flexão pela técnica da birrefringência acústica é apresentado. Como valor de referência, para fins do cálculo da exatidão alcançada, utilizou-se uma solução analítica derivada da teoria clássica da resistência de materiais. Os resultados encontrados indicaram uma razoável precisão, com um erro relativo menor do que 8%, e uma baixa exatidão, menor do que 57%, nos pontos da viga sujeitos às maiores tensões principais de flexão.

Um modelo determinístico para avaliação do controle biológico de praga de citros.

O objetivo desse trabalho é apresentar um modelo matemático compartimental e determinístico para analisar a dinâmica populacional do inseto-praga Phyllocnistis citrella em interação com seus inimigos naturais Galeopsomyia fausta (parasito nativo) e Ageniaspis citricola (parasito exótico). Com o estudo analítico do sistema de equações diferenciais que descreve o modelo, avalia-se a eficiência do controle biológico da P.citrella através da possibilidade de coexistência entre as espécies e da competição entre os dois inimigos naturais.

Optimal control and numerical optimization applied to epidemiological models

A relação entre a epidemiologia, a modelação matemática e as ferramentas computacionais permite construir e testar teorias sobre o desenvolvimento e combate de uma doença. Esta tese tem como motivação o estudo de modelos epidemiológicos aplicados a doenças infeciosas numa perspetiva de Controlo Ótimo, dando particular relevância ao Dengue. Sendo uma doença tropical e subtropical transmitida por mosquitos, afecta cerca de 100 milhões de pessoas por ano, e é considerada pela Organização Mundial de Saúde como uma grande preocupação para a saúde pública. Os modelos matemáticos desenvolvidos e testados neste trabalho, baseiam-se em equações diferenciais ordinárias que descrevem a dinâmica subjacente à doença nomeadamente a interação entre humanos e mosquitos. É feito um estudo analítico dos mesmos relativamente aos pontos de equilíbrio, sua estabilidade e número básico de reprodução. A propagação do Dengue pode ser atenuada através de medidas de controlo do vetor transmissor, tais como o uso de inseticidas específicos e campanhas educacionais. Como o desenvolvimento de uma potencial vacina tem sido uma aposta mundial recente, são propostos modelos baseados na simulação de um hipotético processo de vacinação numa população. Tendo por base a teoria de Controlo Ótimo, são analisadas as estratégias ótimas para o uso destes controlos e respetivas repercussões na redução/erradicação da doença aquando de um surto na população, considerando uma abordagem bioeconómica. Os problemas formulados são resolvidos numericamente usando métodos diretos e indiretos. Os primeiros discretizam o problema reformulando-o num problema de optimização não linear. Os métodos indiretos usam o Princípio do Máximo de Pontryagin como condição necessária para encontrar a curva ótima para o respetivo controlo. Nestas duas estratégias utilizam-se vários pacotes de software numérico. Ao longo deste trabalho, houve sempre um compromisso entre o realismo dos modelos epidemiológicos e a sua tratabilidade em termos matemáticos.

Descarga costeira de águas subterrâneas com elevado conteúdo em nitratos para a Ria Formosa

Dissertação de Mestrado, Estudos Marinhos e Costeiros, Faculdade de Ciências do Mar e Ambiente, Universidade do Algarve, 2007

Modelos de avaliação de hipotecas comerciais: soluções analíticas

Tese de dout., Gestão, Faculdade de Economia, Universidade do Algarve, 2005

Modelação do comportamento térmico do betão em massa : aplicação a uma barragem de BCC

A presente dissertação tem como objectivo efectuar o estudo do comportamento térmico de estruturas com betão em massa através do Método de Schmidt e compará-lo com outros métodos. A análise do comportamento do betão reveste-se de extrema importância em diversos domínios, de entre os quais se pode referir o que visa evitar a fissuração de origem térmica decorrente do calor de hidratação. Pelo motivo enunciado torna-se necessária dispor de metodologias de análise que requerem a mobilização de meios de cálculo. Existem diferentes métodos para o cálculo das temperaturas: uns mais simples, outros mais complexos, uns mais rápidos, outros mais demorados, uns mais onerosos e outros mais económicos. O conhecimento dos vários métodos permite adequá-los à especificidade de cada obra. Para realizar os vários procedimentos pode recorrer-se a programas informáticos de cálculo, como por o Excel ou utilizar Software Específico. Nesta dissertação apresenta-se uma metodologia de cálculo baseada no Método de Schmidt e faz-se a comparação com o Método de Solução Analítica. Utiliza-se também este Método no cálculo das temperaturas das várias camadas de uma barragem. O Método de Schmidt é ainda aplicado a um caso real de uma barragem da Tailândia construída em Betão Compactado com Cilindros e os resultados obtidos são comparados com outros Métodos de Análise Numérica de Barragens.

Inversão numérica da transformada de Laplace por polinômios trigonométricos e de Laguerre

Neste trabalho são desenvolvidos métodos numéricos para inversão da transformada de Laplace, fazendo-se uso de polinômios trigonométricos e de Laguerre. Sua utilização é ilustrada num problema de fronteira móvel da área de engenharia nuclear, através do algoritmo computacional ALG-619. Uma revisão dos aspectos analíticos básicos da transformada de Laplace e sua utilização na resolução de equações diferenciais parciais é apresentada de maneira suscinta.