Resolução numérica de escoamentos compressíveis empregando um método de partículas livre de malhas e o processamento em paralelo (CUDA)

Os métodos numéricos convencionais, baseados em malhas, têm sido amplamente aplicados na resolução de problemas da Dinâmica dos Fluidos Computacional. Entretanto, em problemas de escoamento de fluidos que envolvem superfícies livres, grandes explosões, grandes deformações, descontinuidades, ondas de choque etc., estes métodos podem apresentar algumas dificuldades práticas quando da resolução destes problemas. Como uma alternativa viável, existem os métodos de partículas livre de malhas. Neste trabalho é feita uma introdução ao método Lagrangeano de partículas, livre de malhas, Smoothed Particle Hydrodynamics (SPH) voltado para a simulação numérica de escoamentos de fluidos newtonianos compressíveis e quase-incompressíveis. Dois códigos numéricos foram desenvolvidos, uma versão serial e outra em paralelo, empregando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA), que possibilita o processamento em paralelo empregando os núcleos das Graphics Processing Units (GPUs) das placas de vídeo da NVIDIA Corporation. Os resultados numéricos foram validados e a eficiência computacional avaliada considerandose a resolução dos problemas unidimensionais Shock Tube e Blast Wave e bidimensional da Cavidade (Shear Driven Cavity Problem).

Modelagem computacional aplicada na estimação de parâmetros cinéticos em processos químicos e cromatográficos

A estimação de parâmetros cinéticos em processos químicos e cromatográficos utilizando técnicas numéricas assistidas por computadores tem conduzido para melhoria da eficiência e o favorecimento da compreensão das fenomenologias envolvidas nos mesmos. Na primeira parte deste trabalho será realizada a modelagem computacional do processo de produção de biodiesel via esterificação, sendo que, o método de otimização estocástica Random Restricted Window (R2W) será correlacionado com os dados experimentais da produção de biodiesel a partir da esterificação do ácido láurico com etanol anidro na presença do catalisador ácido nióbico (Nb2O5). Na segunda parte do mesmo será realizada a modelagem computacional do processo de cromatografia de adsorção (batch process) onde serão correlacionados os dados provenientes dos modelos cinéticos de HASHIM, CHASE e IKM2 com os dados experimentais da adsorção de amoxicilina com quitosana, e também serão correlacionados os dados experimentais da adsorção de Bovine Serum Albumin (BSA) com Streamline DEAE com os dados provenientes de uma nova aplicação do método R2W mediante a implementação de um modelo cinético reversível. Ademais, as constantes cinéticas para cada processo supracitado serão estimadas levando em consideração o valor mínimo da função resíduos quadrados.

Uma dedução dos critérios de multicriticalidade e o cálculo de pontos tricríticos via otimização global

Uma dedução dos critérios de multicriticalidade para o cálculo de pontos críticos de qualquer ordem representa a formalização de ideias utilizadas para calcular pontos críticos e tricríticos e ainda amplia tais ideias. De posse desta dedução pode-se compreender os critérios de tricriticalidade e, com isso, através de uma abordagem via problema de otimização global pode-se fazer o cálculo de pontos tricríticos utilizando um método numérico adequado de otimização global. Para evitar um excesso de custo computacional com rotinas numéricas utilizou-se aproximações na forma de diferenças finitas dos termos que compõem a função objetivo. Para simular a relação P v - T optou-se pela equação de estado cúbica de Peng-Robinson e pela regra clássica de fluidos de van der Vaals, para modelagem do problema também se calculou os tensores de ordem 2, 3, 4 e 5 da função do teste de estabilidade. Os resultados obtidos foram comparados com dados experimentais e por resultados obtidos com outros autores que utilizaram métodos numéricos, equação de estado ou abordagem diferente das utilizadas neste trabalho.

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.

Inversão de funções do plano no plano aplicada ao cálculo de azeótropos

Azeotropia é um fenômeno termodinâmico onde um líquido em ebulição produz um vapor com composição idêntica. Esta situação é um desafio para a Engenharia de Separação, já que os processos de destilação exploram as diferenças entre as volatilidades relativas e, portanto, um azeótropo pode ser uma barreira para a separação. Em misturas binárias, o cálculo da azeotropia é caracterizado por um sistema não-linear do tipo 2 × 2. Um interessante e raro caso é o denominado azeotropia dupla, que pode ser verificado quando este sistema não-linear tem duas soluções, correspondendo a dois azeótropos distintos. Diferentes métodos tem sido utilizados na resolução de problemas desta natureza, como métodos estocásticos de otimização e as técnicas intervalares (do tipo Newton intervalar/bisseção generalizada). Nesta tese apresentamos a formulação do problema de azeotropia dupla e uma nova e robusta abordagem para a resolução dos sistemas não-lineares do tipo 2 × 2, que é a inversão de funções do plano no plano (MALTA; SALDANHA; TOMEI, 1996). No método proposto, as soluções são obtidas através de um conjunto de ações: obtenção de curvas críticas e de pré-imagens de pontos arbritários, inversão da função e por fim, as soluções esperadas para o problema de azeotropia. Esta metodologia foi desenvolvida para resolver sistemas não-lineares do tipo 2 × 2, tendo como objetivo dar uma visão global da função que modela o fenômeno em questão, além, é claro, de gerar as soluções esperadas. Serão apresentados resultados numéricos para o cálculo dos azeótropos no sistema benzeno + hexafluorobenzeno a baixas pressões por este método de inversão. Como ferramentas auxiliares, serão também apresentados aspectos numéricos usando aproximações clássicas, tais como métodos de Newton com técnicas de globalização e o algorítmo de otimização não-linear C-GRASP, para efeito de comparação.

Simulação numérica de dutos enterrados, submetidos à perda de apoio e elevação localizada.

Os recentes desastres ocorridos no país, como o rompimento da adutora em Campo Grande e os desastres relacionados às enchentes urbanas, mostram a necessidade de desenvolvimento de pesquisas científicas que auxiliem na compreensão e no dimensionamento das estruturas projetadas para atender a demanda da população. Os métodos analíticos e experimentais mais utilizados possuem algumas limitações de ordem teórica ou prática. Por outro lado, os métodos numéricos, capazes de simular etapas construtivas e envolver materiais com diferentes modelos constitutivos numa mesma análise, buscam atender às necessidades práticas dos projetos de geotecnia e, ao mesmo tempo, complementam os modelos analíticos e experimentais. Nesse trabalho foram realizadas comparações entre resultados obtidos em ensaios experimentais e resultados extraídos do modelo computacional, buscando aumentar a compreensão sobre a interação solo-estrutura em relação à distribuição de tensões mobilizadas e aos deslocamentos e deformações provocados. A simulação numérica foi feita com a utilização do PLAXIS/3D, software de análise geotécnica baseado no método dos elementos finitos. Os ensaios foram confeccionados na Escola de Engenharia de São Carlos/USP por Costa (2005) e envolveram dutos enterrados submetidos à perda de apoio ou elevação localizada. O estudo experimental foi realizado através de modelos físicos compostos por um maciço de areia pura, contendo um tubo repousando sobre um alçapão no centro do vão. Os modelos físicos foram equipados com instrumental capaz de medir as deflexões e as deformações específicas ao longo do duto, além das tensões totais no maciço de solo circundante e na base do equipamento.

Experimentos numéricos para deposição de parafinas com modelo Multisólido

Este trabalho objetiva a construção de estruturas robustas e computacionalmente eficientes para a solução do problema de deposição de parafinas do ponto de vista do equilíbrio sólido-líquido. São avaliados diversos modelos termodinâmicos para a fase líquida: equação de estado de Peng-Robinson e os modelos de coeficiente de atividade de Solução Ideal, Wilson, UNIQUAC e UNIFAC. A fase sólida é caracterizada pelo modelo Multisólido. A previsão de formação de fase sólida é inicialmente prevista por um teste de estabilidade termodinâmica. Posteriormente, o sistema de equações não lineares que caracteriza o equilíbrio termodinâmico e as equações de balanço material é resolvido por três abordagens numéricas: método de Newton multivariável, método de Broyden e método Newton-Armijo. Diversos experimentos numéricos foram conduzidos de modo a avaliar os tempos de computação e a robustez frente a diversos cenários de estimativas iniciais dos métodos numéricos para os diferentes modelos e diferentes misturas. Os resultados indicam para a possibilidade de construção de arcabouços computacionais eficientes e robustos, que podem ser empregados acoplados a simuladores de escoamento em dutos, por exemplo.

Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.

Estrutura populacional e avaliação de métodos de controle da espécie exótica invasora Artocarpus heterophyllus Lamk. (Moraceae) no Parque Estadual da Ilha Grande, Angra do Reis, RJ Rio de Janeiro 2011

Atualmente, com o aumento no número de estudos sobre invasões biológicas, sabemos como ocorre o processo, suas causas e conseqüências. A prevenção ainda é a melhor maneira de limitar e diminuir o aumento no número de problemas associados às espécies exóticas e invasoras biológicas. A jaqueira, Artocarpus heterophyllus Lamk. (Moraceae), é uma espécie exótica invasora que foi introduzida no Brasil no período colonial. A principal estratégia de controle de A. heterophyllus tem sido o método mecânico conhecido como anelamento e o arranque de plântulas. Utilizando o Manual da TNC para Controle de Espécies Invasoras, este estudo objetivou propor um novo método de controle da jaqueira, e caracterizar a estrutura populacional da espécie. O estudo foi conduzido na Ilha Grande, localizada no município de Angra dos Reis, RJ, que é coberta por Mata Atlântica em diferentes estágios sucessionais. Testou-se um novo método químico que consistiu na injeção de herbicida Garlon diluído a 4% no tronco de árvores com DAP>15 cm. Ao todo 684 indivíduos distribuídos em 10 parcelas medindo 0,64 ha cada foram encontrados. A densidade media encontrada foi de 107 ind. / ha, com densidades variando entre 3340 ind. / ha, na classe Juvenis 1 (DAP < 5 cm) a 13 ind / ha na classe Adltos 2 (20.1 25.0). A área basal média encontrada foi de 3,692 m / ha. Os resultados mostraram que o método químico foi mais eficiente que o método mecânico. Após 60, 150 e 240 dias do tratamento inicial, os métodos diferiram significativamente. Os resultados demonstram que não há correlação entre a eficiência do método mecânico em relação ao DAP. Entretanto o método químico é dependente do DAP. Os resultados das taxas de mortalidade foram significantes para o tempo de resposta aos 60 dias (p = 0,009), 150 dias (p = 0,039) e 240 dias (p = 0.013), após teste estatístico Kruskal Wallis. As vantagens do método químico em relação ao mecânico são claras, onde menos dinheiro é gasto e mais resultados são gerados.

Avaliação de ligações soldadas de perfis tubulares em estruturas de aço através do método dos elementos finitos.

A necessidade de se obter soluções para os diversos problemas estruturais na engenharia, associado ao amplo crescimento e difusão da construção metálica fez surgir, por volta dos anos 60, os perfis tubulares. As ligações soldadas entre perfis tubulares são as ligações de maior difusão nas construções com este tipo de perfil, sendo necessário portanto, estudos mais detalhados que tornem a utilização de ligações soldadas mais seguras e otimizadas. O Eurocode 3 possui um item específico para as ligações entre perfis tubulares, em especial as ligações soldadas devido a suas particularidades. Infelizmente, a norma brasileira NBR8800 não aborda as ligações envolvendo perfis tubulares. Para a execução da análise numérica no presente trabalho, faz-se necessário uma revisão bibliográfica. Com estes resultados, foram feitas modelagens de ligações tubulares com o método dos elementos finitos de forma a otimizar os modelos a serem utilizados em uma análise paramétrica futura. Desta forma, foram desenvolvidos dois modelos numéricos, um considerando ligação tipo T entre perfis tubulares quadrados e outro para uma ligação tipo K entre perfis circulares. Estes modelos foram caracterizados através de elementos de casca com seis graus de liberdade por nó considerando-se adicionalmente o efeito de membrana. A análise não-linear realizada considerou a não-linearidade do material através do critério de plastificação de Von Mises através de uma lei constitutiva tensão versus deformação bi-linear e a não-linearidade geométrica foi atribuída através da formulação de Lagrange atualizado. Dentre as principais conclusões obtidas no presente trabalho, pode-se citar que os resultados para as ligações tipo T, o Eurocode 3, fornece resultados que precisam ser observados com cautela. Todavia, para as ligações do tipo K, os resultados numéricos mostraram-se sempre inferiores aos valores através do Eurocode 3, representando um dimensionamento a favor da segurança.

Avaliação de ligações entre vigas de perfil I ou H e colunas em perfis tubulares circulares.

Considerando uma nova realidade com o incremento do uso de perfis tubulares, este trabalho apresenta uma análise de ligações soldadas tipo T com perfis tubulares circulares (CHS) para a coluna e perfil I ou H para a viga com flexão no plano, efetuado com base nas normas EC3, CIDECT e a NBR 16239 comparando com o critério de deformação limite, proposto por Lu et al., através de um modelo em elementos finitos desenvolvido no programa Ansys versão 12.0. A não-linearidade geométrica foi introduzida no modelo através da Formulação de Lagrange Atualizado. Com base nos resultados numéricos avaliados foram traçadas curvas momento-rotação para cada modelo com o objetivo de obter o momento resistente de cada ligação bem como a classificação da ligação quanto a capacidade de rotação, a influência dos parâmetros geométricos em cada modelo e o modo de falha que controlará o dimensionamento da ligação.

Análise espectral das equações de transporte de nêutrons na formulação de ordenadas discretas em meios multiplicativos

É presentada nesta dissertação uma análise espectral das equações de transporte de nêutrons, independente do tempo, em geometria unidimensional e bidimensional, na formulação de ordenadas discretas (SN), utilizando o modelo de uma velocidade e multigrupo, considerando meios onde ocorrem o fenômeno da fissão nuclear. Esta análise espectral constitui-se na resolução de problemas de autovalores e respectivos autovetores, e reproduz a expressão para a solução geral analítica local das equações SN (para geometria unidimensional) ou das equações nodais integradas transversalmente (geometria retangular bidimensional) dentro de cada região homogeneizada do domínio espacial. Com a solução geral local determinada, métodos numéricos, tais como os métodos de matriz de resposta SN, podem ser derivados. Os resultados numéricos são gerados por programas de computadores implementados em MatLab, versão 2012, a fim de verificar a natureza dos autovalores e autovetores correspondentes no espaço real ou complexo.

O pré-aquecimento do universo: da ressonância paramétrica à turbulência

A inflação consegue dar conta de uma série de problemas do Modelo padrão da Cosmologia, preservando ainda o sucesso do modelo do Big Bang. Na sua versão mais simples, a inflação é controlada por um campo escalar, o ínflaton, que faz com que o universo se expanda exponencialmente. Após, o ínflaton decai e ocorre o reaquecimento do universo. Contudo, alguns autores apontam a existência de uma fase intermediária, chamada de pré-aquecimento. O decaimento do ínflaton possui uma rica dinâmica não-linear. No primeiro estágio, a ressonância paramétrica promove o crescimento exponencial de alguns modos do ínflaton. Isto altera a dinâmica do modo homogêneo do ínflaton, promovendo uma reestruturação das cartas de ressonâncias da equação de movimento dos modos perturbativos. Desta forma, ocorre a transferência de energia para estes modos, até que o universo termaliza. Esta transferência de energia é típica de um sistema turbulento. Por se tratar de uma evolução não-linear, torna-se conveniente a implementação computacional de métodos numéricos. Neste contexto, os métodos espectrais têm se mostrado uma excelente ferramenta para estudar este tipo de sistema. Esta dissertação apresenta os resultados do esquema numérico desenvolvido para o modelo com potencial quártico, que será a base para os demais estudos a serem desenvolvidos. Como mostrado, este esquema é extremamente preciso e eficiente.

Simulação do processo de transferência de calor não linear condução-radiação por meio de um esquema linear em diferenças finitas.

Neste trabalho o processo não linear de transmissão de calor condução-radiação é abordado num contexto bidimensional plano e simulado com o uso de um esquema linear em diferenças finitas. O problema original é tratado como o limite de uma sequencia de problemas lineares, do tipo condução-convecção. Este limite, cuja existência é comprovada, é facilmente obtido a partir de procedimentos básicos, accessíveis a qualquer estudante de engenharia, permitindo assim o emprego de hipóteses mais realistas, já que não se tem o limitante matemático para a abordagem numérica de uma equação diferencial parcial elíptica. Neste trabalho foi resolvido o problema de condução de calor em regime permanente em uma placa com condições de contorno convectivas e radioativas utilizando-se o software MatLab, vale ressaltar, que a mesma metodologia é aplicável para geometrias mais complexas.