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

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).

Desenvolvimento de um simulador numérico empregando o método Smoothed Particle Hydrodynamics para a resolução de escoamentos incompressíveis. Implementação computacional em paralelo (CUDA)

Neste trabalho, foi desenvolvido um simulador numérico baseado no método livre de malhas Smoothed Particle Hydrodynamics (SPH) para a resolução de escoamentos de fluidos newtonianos incompressíveis. Diferentemente da maioria das versões existentes deste método, o código numérico faz uso de uma técnica iterativa na determinação do campo de pressões. Este procedimento emprega a forma diferencial de uma equação de estado para um fluido compressível e a equação da continuidade a fim de que a correção da pressão seja determinada. Uma versão paralelizada do simulador numérico foi implementada usando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA) da NVIDIA Corporation. Foram simulados três problemas, o problema unidimensional do escoamento de Couette e os problemas bidimensionais do escoamento no interior de uma Cavidade (Shear Driven Cavity Problem) e da Quebra de Barragem (Dambreak).

Simulação de escoamentos incompressíveis empregando o método Smoothed Particle Hydrodynamics utilizando algoritmos iterativos na determinação do campo de pressões

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas.

Simulação numérica do escoamento de água em áreas alagáveis da floresta amazônica com a utilização da estrutura de dados Autonomous Leaves Graph

A Amazônia exibe uma variedade de cenários que se complementam. Parte desse ecossistema sofre anualmente severas alterações em seu ciclo hidrológico, fazendo com que vastos trechos de floresta sejam inundados. Esse fenômeno, entretanto, é extremamente importante para a manutenção de ciclos naturais. Neste contexto, compreender a dinâmica das áreas alagáveis amazônicas é importante para antecipar o efeito de ações não sustentáveis. Sob esta motivação, este trabalho estuda um modelo de escoamento em áreas alagáveis amazônicas, baseado nas equações de Navier-Stokes, além de ferramentas que possam ser aplicadas ao modelo, favorecendo uma nova abordagem do problema. Para a discretização das equações é utilizado o Método dos Volumes Finitos, sendo o Método do Gradiente Conjugado a técnica escolhida para resolver os sistemas lineares associados. Como técnica de resolução numérica das equações, empregou-se o Método Marker and Cell, procedimento explícito para solução das equações de Navier-Stokes. Por fim, as técnicas são aplicadas a simulações preliminares utilizando a estrutura de dados Autonomous Leaves Graph, que tem recursos adaptativos para manipulação da malha que representa o domínio do problema

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.

Uma formulação implícita para o método Smoothed Particle Hydrodynamics

Em uma grande gama de problemas físicos, governados por equações diferenciais, muitas vezes é de interesse obter-se soluções para o regime transiente e, portanto, deve-se empregar técnicas de integração temporal. Uma primeira possibilidade seria a de aplicar-se métodos explícitos, devido à sua simplicidade e eficiência computacional. Entretanto, esses métodos frequentemente são somente condicionalmente estáveis e estão sujeitos a severas restrições na escolha do passo no tempo. Para problemas advectivos, governados por equações hiperbólicas, esta restrição é conhecida como a condição de Courant-Friedrichs-Lewy (CFL). Quando temse a necessidade de obter soluções numéricas para grandes períodos de tempo, ou quando o custo computacional a cada passo é elevado, esta condição torna-se um empecilho. A fim de contornar esta restrição, métodos implícitos, que são geralmente incondicionalmente estáveis, são utilizados. Neste trabalho, foram aplicadas algumas formulações implícitas para a integração temporal no método Smoothed Particle Hydrodynamics (SPH) de modo a possibilitar o uso de maiores incrementos de tempo e uma forte estabilidade no processo de marcha temporal. Devido ao alto custo computacional exigido pela busca das partículas a cada passo no tempo, esta implementação só será viável se forem aplicados algoritmos eficientes para o tipo de estrutura matricial considerada, tais como os métodos do subespaço de Krylov. Portanto, fez-se um estudo para a escolha apropriada dos métodos que mais se adequavam a este problema, sendo os escolhidos os métodos Bi-Conjugate Gradient (BiCG), o Bi-Conjugate Gradient Stabilized (BiCGSTAB) e o Quasi-Minimal Residual (QMR). Alguns problemas testes foram utilizados a fim de validar as soluções numéricas obtidas com a versão implícita do método SPH.

Simulação numérica da cavitação em turbomáquinas usando uma formulação Euler-Lagrange.

Turbomáquinas são máquinas operacionais que transferem energia mecânica entre um rotor e um fluido. Estas máquinas têm muitas aplicações industriais. Um dos componentes de uma turbomáquina responsável pela transferência da energia, ou receber a rotação do eixo e transformar em energia de fluido em caso de bomba ou transferir a energia do fluido para o eixo em caso de uma turbina, é o impelidor ou rotor. O fenómeno da cavitação envolve escoamento bifásico: o líquido a ser bombeado e as bolhas de vapor que são formadas durante o processo de bombeamento. O processo de formação dessas bolhas é complexo, mas ocorre principalmente devido a presença de regiões de pressões muito baixas. O colapso dessas bolhas pode muitas vezes levar a deterioração do material, dependendo da intensidade ou da velocidade de colapso das bolhas. O principal objetivo deste trabalho foi estudar o comportamento hidrodinâmico do escoamento nos canais do impelidor de uma turbomáquina do tipo radial usando recursos de fluidodinâmica computacional (CFD). Uma abordagem Euler-Lagrange acoplada com o modelo da equação de Langevin foi empregada para estimar a trajetória das bolhas. Resultados das simulações mostram as particularidades de um escoamento líquido-bolha de vapor passando em um canal de geometria curva, fornecendo assim informações que podem nos ajudar na prevenção da cavitação nessas máquinas.

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.

Estudo de fluxos numéricos para modelos magneto-hidrodinâmicos utilizando o código Flash

A teoria magneto-hidrodinâmicos permite a estruturação de modelos computacionais, designados modelos MHDs, que são uma extensão da dinâmica dos fluidos para lidar com fluidos eletricamente carregados, tais como os plasmas, em que se precisa considerar os efeitos de forças eletromagnéticas. Tais modelos são especialmente úteis quando o movimento exato de uma partícula não é de interesse, sendo que as equações descrevem as evoluções de quantidades macroscópicas. Várias formas de modelos MHD têm sido amplamente utilizadas na Física Espacial para descrever muitos tipos diferentes de fenômenos de plasma, tais como reconexão magnética e interações de ventos estelares com diferentes objetos celestiais. Neste trabalho, o objetivo é analisar o comportamento de diversos fluxos numéricos em uma discretização de volumes finitos de um modelo numérico de MHD usando um esquema de malha entrelaçada sem separação direcional considerando alguns casos testes. Para as simulações, utiliza-se o código Flash, desenvolvido pela Universidade de Chicago, por ser um código de amplo interesse nas simulações astrofísicas e de fenômenos no espaço próximo à Terra. A metodologia consiste na inclusão de um fluxo numérico, permitindo melhoria com respeito ao esquema HLL.

Dinâmica populacional de Halodule wrightii Ascherson e sua fauna e flora associada em Cabo Frio RJ

Fanerógamas marinhas (gramas marinhas) são plantas com flores adaptadas ao ambiente marinho costeiro da maioria dos continentes do mundo. As gramas marinhas formam extensos bancos e proveem valiosos recursos em águas costeiras rasas em todo o mundo, servindo de alimento e berçário para espécies importantes de pescados comerciais e recreacionais. Nesse estudo foi realizada uma revisão sobre o estado de conhecimento das fanerógamas marinhas no Brasil até o presente momento; avaliou-se a importância do monitoramento em longo prazo e a influência de fatores ambientais, como o número de manchas solares; pesquisou-se também a distribuição espacial da grama marinha, bem como a fauna e flora associada; e o crescimento de Halodule wrightii em duas condições ambientais extremas (exposta no ciclo de maré baixa e permanentemente submersa). A revisão bibliográfica sobre as gramas marinhas foi abrangente e verificou a existência de algumas lacunas no conhecimento. Através do monitoramento a longo prazo pôde ser observado que o número de manchas solares tem forte relação negativa sobre a altura do dossel das gramas marinhas de região entre marés. A variação de marés na região de mediolitoral está relacionada diretamente com a distribuição espacial de Halodule wrightii e, consequentemente na distribuição da fauna e flora associada. A diferença de crescimento nos eixos de Halodule wrightii em condições ambientais diferentes é compensada pelas variações nas características de distribuição da planta no ambiente, tais como a altura do dossel, a densidade e biomassa de eixos. O monitoramento a longo prazo pode permitir a tomada de ações que auxiliem no manejo e na recuperação desses importantes habitats costeiros.

Dinâmica de crescimento de Cedrela odorata L. (Meliaceae) na Floresta Atlântica do Estado do Rio de Janeiro: fenologia, atividade cambial e dendrocronologia

Ao longo do século XX, poucos estudos de dendrocronologia foram desenvolvidos com espécies de ambientes tropicais, em função da crença de que as condições climáticas nessas regiões não apresentavam variações suficientemente marcantes e regulares para induzir um ritmo anual de crescimento radial. A realização de trabalhos sobre esse tema nas últimas décadas revelou que a formação de anéis de crescimento anuais nos trópicos pode estar associada a fatores diversos, como: existência de estação seca bem definida, ocorrência de inundações sazonais, respostas ao comportamento fenológico, respostas ao fotoperíodo e a ritmos endógenos. O presente estudo tem por objetivo compreender a dinâmica de crescimento radial de uma espécie da Mata Atlântica se desenvolvendo em ambiente natural. Para tanto, propôs-se: i) investigar a periodicidade da atividade cambial e dos fatores que a influenciam; ii) estimar a idade e taxa de crescimento diamétrico e iii) correlacionar os fatores ambientais com os anéis de crescimento, em indivíduos de Cedrela odorata L. Para o estudo da atividade cambial, foram obtidas amostras de caule a 1,30 m do solo, contendo periderme, faixa cambial e xilema e floema secundários, por métodos não destrutivos. A fenologia vegetativa e a frutificação dos indivíduos amostrados foram acompanhadas durante todo o período do experimento. O material coletado foi processado segundo técnicas usuais em Anatomia Vegetal e analisado sob microscopia óptica e de fluorescência. Os dados de fotoperíodo, precipitação, temperatura e fenologia vegetativa foram correlacionados à atividade cambial. Para o estudo dos anéis de crescimento, as coletas também foram realizadas a 1,30 m do solo, por meio de sonda de Pressler. As amostras obtidas foram polidas e analisadas sob microscópio estereoscópio, para demarcação e aferição do número de anéis de crescimento, e a largura dos anéis foi mensurada para a determinação das taxas de crescimento radial. A série histórica de temperatura e precipitação foi correlacionada à cronologia dos anéis de crescimento. Os resultados indicaram que a atividade cambial segue um ritmo anual de crescimento, correlacionado à sazonalidade do fotoperíodo, da precipitação e da fenologia vegetativa. A análise dos anéis de crescimento permitiu estimar a idade dos indivíduos e determinar a taxa média de incremento e as taxas de incremento diamétrico acumulado e incremento médio anual para a espécie no sítio de estudo. Os dados de incremento radial evidenciaram a ausência de relação entre a idade e o diâmetro das árvores. A análise da variação na largura dos anéis não apresentou correlações significativas com os fatores climáticos analisados.

Análise dinâmica de plataformas de aço para produção de petróleo com base na consideração do efeito da interação solo-estrutura

Atualmente, as tendências competitivas do mercado mundial, têm forçado os engenheiros estruturais a desenvolver soluções de projeto que acarretem em menor peso e custo de execução. Uma consequência direta desta nova tendência de projeto é o aumento considerável de problemas relacionados a vibrações de piso indesejadas. Por esta razão, os sistemas estruturais de pisos podem tornar-se vulneráveis a vibrações excessivas, como por exemplo, aquelas induzidas por equipamentos mecânicos (máquinas rotativas). Deste modo, este trabalho objetiva investigar o comportamento dinâmico de uma plataforma de aço para produção de petróleo, localizada na bacia de Santos (campo de Merluza), São Paulo, Brasil. Para tal, investiga-se a influência das ações dinâmicas oriundas dos equipamentos mecânicos localizados sobre os decks metálicos da plataforma. A resposta dinâmica do modelo estrutural foi determinada através de um extenso estudo numérico, a partir da análise de suas frequências naturais, deslocamentos, velocidades e acelerações de pico. Nesta investigação, as cargas dinâmicas provenientes dos equipamentos mecânicos (máquinas rotativas) foram aplicadas sobre o piso metálico do sistema estrutural. Com base obtenção da resposta dinâmica da estrutura (deslocamentos, velocidades e acelerações), foi possível avaliar a performance do modelo estrutural em termos de critérios de conforto humano e das tolerâncias máximas referentes aos equipamentos mecânicos, de acordo com normas e recomendações de projeto.

A complexidade da dinâmica socioambiental em território insular: desafio ao planejamento - Ilha Grande, Angra dos Reis (RJ)

A Ilha Grande situa-se na chamada Costa Verde do Estado do Rio de Janeiro, tem a totalidade de seu território protegido legalmente e vem sendo alvo de especuladores e investidores cujos interesses são, principalmente, explorar o potencial natural do bioma de Mata Atlântica. Observa-se um processo de transformação local em diversos aspectos que englobam o meio ambiente e a sociedade, como a degradação ambiental e o deslocamento da comunidade tradicional. O Poder Público manteve-se distante ou mesmo ausente por décadas, o que possibilitou a instalação de problemas ambientais e a configuração de conflitos na esfera do planejamento e gestão. Considerando-se a relevância socioambiental deste território insular e, em especial, com a implantação do Centro de Estudos Ambientais e Desenvolvimento Sustentável (CEADS) da Universidade do Estado do Rio de Janeiro (UERJ) na Enseada de Dois Rios, o presente estudo tem como objetivo identificar a dinâmica ambiental estabelecida na ilha. Os aspectos relacionados ao uso e ocupação do solo e à fragilidade ambiental são fundamentais para a análise integrada, para a compreensão da configuração atual do território e para a formulação de diretrizes de desenvolvimento adequadas à complexa realidade da Ilha Grande, cuja atividade econômica está pautada no turismo. Identificou-se como as principais causas dos problemas ambientais e dos conflitos de gestão: o caótico ordenamento do território, a falta de infraestrutura de serviços, a dificuldade de mobilização social, a falta de integração entre as políticas públicas e competências administrativas, a precariedade da fiscalização, a impunidade aos que não cumpriram a normatização em vigor e, a não observância da fragilidade ambiental do geossistema Ilha Grande.

A dinâmica dos espaços de cultura e lazer em Botafogo

O processo de terciarização da economia, o qual se faz presente claramente desde os primórdios da cidade do Rio de Janeiro, promoveu profundas transformações na organização interna da cidade, o que implica em novos usos do espaço urbano e gera novos processos na cidade, sobretudo em nossa área de estudo, o bairro de Botafogo. O processo de globalização provoca mudanças significativas na economia mundial. A cidade, portanto se insere nos espaços de globalização e, assim, sofre mudanças que vão se refletir em sua estrutura urbana. Desta forma, o setor de serviços é dinamizado, reforçando, então, o processo de terciarização. Neste caso, se torna evidente a relação entre urbanização as atividades terciárias. Percebe-se, assim, o surgimento de um novo modelo de urbanização, atrelado à presença destas atividades terciárias, as quais serão de suma importância para a constituição de novas centralidades no espaço interno da cidade, como é o caso de Botafogo. Neste contexto, os espaços de cultura e lazer têm fundamental importância, visto que atraem fluxos de consumidores, os quais se aproveitam da infraestrutura proporcionada pelo Estado, sobretudo de transportes, através do metrô e das inúmeras linhas de ônibus que servem o bairro. Os espaços de cultura e lazer possuem uma grande importância no contexto de renovação de Botafogo, uma vez que o bairro apresenta uma grande concentração de equipamentos desta natureza, os quais contribuem para uma maior circulação de pessoas em seu espaço interno. Observam-se duas áreas de concentração destes espaços em Botafogo, as quais se localizam nas extremidades e adjacências do principal eixo viário do bairro, onde podemos encontrar também onde se concentram serviços especializados. Então, podemos dizer que, na esfera do consumo, estes espaços de cultura e lazer surgem como resultado de ações públicas e privadas. Os espaços de cultura e lazer promovem uma nova dinâmica a Botafogo e gera impactos na economia urbana do Rio de Janeiro. São responsáveis pelo surgimento de uma nova centralidade no âmbito da cidade e contribuem, portanto, para o processo de reestruturação urbana da cidade do Rio de Janeiro