Um algoritmo em paralelo para solução de equações diferenciais evolutivas

Este trabalho que envolve matemática aplicada e processamento paralelo: seu objetivo é avaliar uma estratégia de implementação em paralelo para algoritmos de diferenças finitas que aproximam a solução de equações diferenciais de evolução. A alternativa proposta é a substituição dos produtos matriz-vetor efetuados sequencialmente por multiplicações matriz-matriz aceleradas pelo método de Strassen em paralelo. O trabalho desenvolve testes visando verificar o ganho computacional relacionado a essa estratégia de paralelização, pois as aplicacações computacionais, que empregam a estratégia sequencial, possuem como característica o longo período de computação causado pelo grande volume de cálculo. Inclusive como alternativa, nós usamos o algoritmo em paralelo convencional para solução de algoritmos explícitos para solução de equações diferenciais parciais evolutivas no tempo. Portanto, de acordo com os resultados obtidos, nós observamos as características de cada estratégia em paralelo, tendo como principal objetivo diminuir o esforço computacional despendido.

Análise de unidades de pré-purificação de ar por adsorção

Os principais constituintes do ar, nitrogênio, oxigênio e argônio, estão cada vez mais presentes nas indústrias, onde são empregados nos processos químicos, para o transporte de alimentos e processamento de resíduos. As duas principais tecnologias para a separação dos componentes do ar são a adsorção e a destilação criogênica. Entretanto, para ambos os processos é necessário que os contaminantes do ar, como o gás carbônico, o vapor dágua e hidrocarbonetos, sejam removidos para evitar problemas operacionais e de segurança. Desta forma, o presente trabalho trata do estudo do processo de pré-purificação de ar utilizando adsorção. Neste sistema a corrente de ar flui alternadamente entre dois leitos adsorvedores para produzir ar purificado continuamente. Mais especificamente, o foco da dissertação corresponde à investigação do comportamento de unidades de pré-purificação tipo PSA (pressure swing adsorption), onde a etapa de dessorção é realizada pela redução da pressão. A análise da unidade de pré-purificação parte da modelagem dos leitos de adsorção através de um sistema de equações diferenciais parciais de balanço de massa na corrente gasosa e no leito. Neste modelo, a relação de equilíbrio relativa à adsorção é descrita pela isoterma de Dubinin-Astakhov estendida para misturas multicomponentes. Para a simulação do modelo, as derivadas espaciais são discretizadas via diferenças finitas e o sistema de equações diferenciais ordinárias resultante é resolvido por um solver apropriado (método das linhas). Para a simulação da unidade em operação, este modelo é acoplado a um algoritmo de convergência relativo às quatro etapas do ciclo de operação: adsorção, despressurização, purga e dessorção. O algoritmo em questão deve garantir que as condições finais da última etapa são equivalentes às condições iniciais da primeira etapa (estado estacionário cíclico). Desta forma, a simulação foi implementada na forma de um código computacional baseado no ambiente de programação Scilab (Scilab 5.3.0, 2010), que é um programa de distribuição gratuita. Os algoritmos de simulação de cada etapa individual e do ciclo completo são finalmente utilizados para analisar o comportamento da unidade de pré-purificação, verificando como o seu desempenho é afetado por alterações nas variáveis de projeto ou operacionais. Por exemplo, foi investigado o sistema de carregamento do leito que mostrou que a configuração ideal do leito é de 50% de alumina seguido de 50% de zeólita. Variáveis do processo foram também analisadas, a pressão de adsorção, a vazão de alimentação e o tempo do ciclo de adsorção, mostrando que o aumento da vazão de alimentação leva a perda da especificação que pode ser retomada reduzindo-se o tempo do ciclo de adsorção. Mostrou-se também que uma pressão de adsorção maior leva a uma maior remoção de contaminantes.

Modelagem e simulação da propagação de ondas em barras não homogêneas envolvendo materiais elásticos não lineares.

O objetivo deste trabalho é tratar da simulação do fenômeno de propagação de ondas em uma haste heterogênea elástico, composta por dois materiais distintos (um linear e um não-linear), cada um deles com a sua própria velocidade de propagação da onda. Na interface entre estes materiais existe uma descontinuidade, um choque estacionário, devido ao salto das propriedades físicas. Empregando uma abordagem na configuração de referência, um sistema não-linear hiperbólico de equações diferenciais parciais, cujas incógnitas são a velocidade e a deformação, descrevendo a resposta dinâmica da haste heterogénea. A solução analítica completa do problema de Riemann associado são apresentados e discutidos.

Avaliação dos algoritmos de Picard-Krylov e Newton-Krylov na solução da equação de Richards

A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards.

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.

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.

Posicionamento aproximado do estado final para sistemas térmicos descritos pela equação do calor.

Neste trabalho, será considerado um problema de controle ótimo quadrático para a equação do calor em domínios retangulares com condição de fronteira do tipo Dirichlet é nos quais, a função de controle (dependente apenas no tempo) constitui um termo de fonte. Uma caracterização da solução ótima é obtida na forma de uma equação linear em um espaço de funções reais definidas no intervalo de tempo considerado. Em seguida, utiliza-se uma sequência de projeções em subespaços de dimensão finita para obter aproximações para o controle ótimo, o cada uma das quais pode ser gerada por um sistema linear de dimensão finita. A sequência de soluções aproximadas assim obtidas converge para a solução ótima do problema original. Finalmente, são apresentados resultados numéricos para domínios espaciais de dimensão 1.

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.

Abordagens do tipo livre de jacobiana na simulação do escoamento de fluidos compressíveis em meios porosos

O desenvolvimento de software livre de Jacobiana para a resolução de problemas formulados por equações diferenciais parciais não-lineares é de interesse crescente para simular processos práticos de engenharia. Este trabalho utiliza o chamado algoritmo espectral livre de derivada para equações não-lineares na simulação de fluxos em meios porosos. O modelo aqui considerado é aquele empregado para descrever o deslocamento do fluido compressível miscível em meios porosos com fontes e sumidouros, onde a densidade da mistura de fluidos varia exponencialmente com a pressão. O algoritmo espectral utilizado é um método moderno para a solução de sistemas não-lineares de grande porte, o que não resolve sistemas lineares, nem usa qualquer informação explícita associados com a matriz Jacobiana, sendo uma abordagem livre de Jacobiana. Problemas bidimensionais são apresentados, juntamente com os resultados numéricos comparando o algoritmo espectral com um método de Newton inexato livre de Jacobiana. Os resultados deste trabalho mostram que este algoritmo espectral moderno é um método confiável e eficiente para a simulação de escoamentos compressíveis em meios porosos.

Estimativa do erro de discretização analítico na solução de equações diferenciais utilizando o Método de Volumes Finitos

As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado

Distribuição de tensões em próteses parciais fixas posteriores livres de metal com retentores intracoronários: análise em elementos finitos

Para reabilitar a ausência de um elemento dentário posterior, as próteses parciais fixas (PPF) com retentores intracoronários são uma alternativa aos implantes osseointegrados. O objetivo deste estudo foi avaliar a distribuição de tensões nessas próteses com três combinações de materiais: cerâmica de zircônia parcialmente estabilizada por ítria (ZPEI) revestida por cerâmica de fluorapatita (α), cerâmica de dissilicato de lítio (β) ou compósito fibrorreforçado (γ). Na composição α, foram analisadas a presença ou ausência da cerâmica de revestimento na parede cervical das caixas proximais e três variações na área total da seção transversal dos conectores (4 mm de largura x 3,2, 4,2 ou 5,2 mm de altura). Em 8 modelos bidimensionais de elementos finitos, uma carga vertical de 500 N foi aplicada na fossa central do pôntico e as tensões principais máximas (tração) e mínimas (compressão) foram apontadas em MPa. Inicialmente foram avaliados os 6 modelos com PPF de ZPEI e suas variações. Os maiores valores das tensões de tração foram encontrados no terço cervical dos conectores. Quando presente nestas regiões, a cerâmica de revestimento recebeu tensões acima do limite de sua resistência à flexão. Na comparação entre os modelos sem cerâmica de revestimento na parede cervical das caixas proximais, mesmo aquele com conectores de 3,2 x 4 mm, cuja infraestrutura apresentava 2,5 x 3 mm, poderia ser recomendado para uso clínico. Altos valores de tensões de compressão foram registrados entre o terço oclusal e médio dos conectores, correspondente à união entre as cerâmicas, o que poderia ocasionar, devido à flexão, falhas adesivas. Posteriormente, o modelo de ZPEI com a cerâmica de fluorapatita ausente da parede cervical das caixas proximais e área total dos conectores de 4,2 x 4 mm foi comparado aos dois outros materiais com conectores de mesma área. Na PPF de dissilicato de lítio, os valores representaram uma provável violação do limite de sua resistência à flexão. A PPF de compósito fibrorreforçado apresentou tensões bem abaixo do limite de resistência à flexão de sua infraestrutura, mas, como no modelo de ZPEI, tensões compressivas se concentraram com alto valor entre o terço oclusal e médio dos conectores, local de união entre a resina composta e a infraestrutura de fibras. Os resultados mostraram que a cerâmica de dissilicato de lítio e a presença da cerâmica de fluorapatita na parede cervical das caixas proximais deveriam ser contraindicadas para a condição proposta. Parece viável uma área de conectores na infraestrutura de ZPEI com no mínimo 2,5 x 3 mm. A PPF de compósito fibrorreforçado apresenta resistência estrutural para a situação estudada, mas, como também aquelas compostas de ZPEI, aparenta ter como pontos fracos a adesão entre a infraestrutura e o material de cobertura e a própria resistência deste último.

Determinação de integrais primeiras liouvillianas em equações diferenciais ordinárias de segunda ordem

Nesta Tese desenvolvemos várias abordagens "Darbouxianas"para buscar integrais primeiras (elementares e Liouvillianas) de equações diferenciais ordinárias de segunda ordem (2EDOs) racionais. Os algoritmos (semi-algoritmos) que desenvolvemos seguem a linha do trabalho de Prelle e Singer. Basicamente, os métodos que buscam integrais primeiras elementares são uma extensão da técnica desenvolvida por Prelle e Singer para encontrar soluções elementares de equações diferenciais ordinárias de primeira ordem (1EDOs) racionais. O procedimento que lida com 2EDOs racionais que apresentam integrais primeiras Liouvillianas é baseado em uma extensão ao nosso método para encontrar soluções Liouvillianas de 1EDOs racionais. A ideia fundamental por tras do nosso trabalho consiste em que os fatores integrantes para 1-formas polinomiais geradas pela diferenciação de funções elementares e Liouvillianas são formados por certos polinômios denominados polinômios de Darboux. Vamos mostrar como combinar esses polinômios de Darboux para construir fatores integrantes e, de posse deles, determinar integrais primeiras. Vamos ainda discutir algumas implementações computacionais dos semi-algoritmos.

Associação entre baixa estatura e sobrepeso em mulheres brasileiras

A associação entre o déficit nutricional em períodos recentes da vida com as doenças crônico-degenerativas e o excesso de peso, tem sido investigada principalmente em países desenvolvidos. No Brasil a associação tem se mostrado significativa apenas no sexo feminino. A realização de um estudo com a população feminina em idade reprodutiva contribuiria para explicar as diferenças de gênero observadas nas altas prevalências de sobrepeso em mulheres no país. Avaliar a associação entre baixa estatura, sobrepeso e ganho excessivo de peso. Este estudo foi baseado na Pesquisa Nacional de Demografia e Saúde (DHS-96). Trata-se de um estudo seccional de representatividade nacional. Todas as análises foram realizadas agrupando as regiões em Rural; Norte e Nordeste; Sul; Sudeste e Centro-Oeste. A população elegível para esta análise foi de 2.292 mães com idade entre 20 e 45 anos. Foram excluídas da análise: mulheres grávidas no momento do inquérito; com último parto anterior a seis meses ou menos da entrevista; não pesadas e não medidas. Os resultados mostraram associação entre a baixa estatura e sobrepeso anterior a primeira gravidez e indicam que a baixa estatura pode estar associada com ganho excessivo de peso no período reprodutivo. Entretanto, os diferenciais observados em relação a região Centro-Sul, consideram a possibilidade da baixa estatura ser um marcador de diferentes contextos sócio-econômicos associados a evolução dos problemas nutricionais no país.

Educação ambiental na educação formal: do paradigma moderno ao paradigma da complexidade

A educação ambiental vem sendo disseminada por grande parte das escolas do país, conforme atestam pesquisas oficiais, mas uma lenta transformação a acompanha no que diz respeito à mudança dos princípios do modelo de desenvolvimento industrial, que enfatiza o consumismo como estratégia de reprodução. Ao mesmo tempo em que as pesquisas sugerem uma ampliação da consciência ambiental da sociedade como um todo, percebe-se um acirramento dos problemas socioambientais e soluções sendo alcançadas em escala inferior ao patamar considerado desejado. Onde estaria, então, o problema? Qual a razão desta possível defasagem educação ambiental x resultados? Foi com o intuito de responder a estas perguntas que essa tese foi realizada. Parte-se de uma percepção de que a Educação Ambiental praticada na escola reflete e acentua o paradigma moderno hegemônico, que se baseia em uma concepção dualista homem / natureza, em que uma está a serviço do outro, sem promover o questionamento sobre os desdobramentos da adoção dos valores da sociedade de consumo no desequilíbrio da vida no planeta. Essa dicotomia entre cultura e natureza, sociedade e ciência, sujeito e objeto se reflete na busca de soluções parciais, incompletas, visando-se apenas a uma parte do problema socioambiental, que não é visto como um sistema complexo. A educação ambiental realizada nessas bases dificulta a reunião das condições necessárias à mudança das estruturas da atual sociedade brasileira e à busca das soluções dos seus problemas socioambientais. Nessa perspectiva, esse trabalho objetiva criar subsídios para um caminho para a educação ambiental que possa contribuir para uma visão complexa da realidade e dos problemas socioambientais, na busca de soluções mais abrangentes. Para isto procura entender: a) como a teoria da complexidade poderia colaborar para esta mudança; b) em quais modelos práticos e teóricos a Educação Ambiental se dá no Brasil, ou seja, quais as diversas tendências da educação ambiental brasileira; e c) como esses modelos se expressam nas práticas dos professores analisando-se artigos publicados em anais de seminários, congressos e/ou encontros sobre o tema.

Modelando evolução por endossimbiose

Nesta dissertação é apresentada uma modelagem analítica para o processo evolucionário formulado pela Teoria da Evolução por Endossimbiose representado através de uma sucessão de estágios envolvendo diferentes interações ecológicas e metábolicas entre populações de bactérias considerando tanto a dinâmica populacional como os processos produtivos dessas populações. Para tal abordagem é feito uso do sistema de equações diferenciais conhecido como sistema de Volterra-Hamilton bem como de determinados conceitos geométricos envolvendo a Teoria KCC e a Geometria Projetiva. Os principais cálculos foram realizados pelo pacote de programação algébrica FINSLER, aplicado sobre o MAPLE.