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.

Modelagem e simulação do escoamento imiscível em meios porosos fractais descritos pela equação de Kozeny-Carman Generalizada

O presente trabalho trata do escoamento bifásico em meios porosos heterogêneos de natureza fractal, onde os fluidos são considerados imiscíveis. Os meios porosos são modelados pela equação de Kozeny-Carman Generalizada (KCG), a qual relaciona a porosidade com a permeabilidade do meio através de uma nova lei de potência. Esta equação proposta por nós é capaz de generalizar diferentes modelos existentes na literatura e, portanto, é de uso mais geral. O simulador numérico desenvolvido aqui emprega métodos de diferenças finitas. A evolução temporal é baseada em um esquema de separação de operadores que segue a estratégia clássica chamada de IMPES. Assim, o campo de pressão é calculado implicitamente, enquanto que a equação da saturação da fase molhante é resolvida explicitamente em cada nível de tempo. O método de otimização denominado de DFSANE é utilizado para resolver a equação da pressão. Enfatizamos que o DFSANE nunca foi usado antes no contexto de simulação de reservatórios. Portanto, o seu uso aqui é sem precedentes. Para minimizar difusões numéricas, a equação da saturação é discretizada por um esquema do tipo "upwind", comumente empregado em simuladores numéricos para a recuperação de petróleo, o qual é resolvido explicitamente pelo método Runge-Kutta de quarta ordem. Os resultados das simulações são bastante satisfatórios. De fato, tais resultados mostram que o modelo KCG é capaz de gerar meios porosos heterogêneos, cujas características permitem a captura de fenômenos físicos que, geralmente, são de difícil acesso para muitos simuladores em diferenças finitas clássicas, como o chamado fenômeno de dedilhamento, que ocorre quando a razão de mobilidade (entre as fases fluidas) assume valores adversos. Em todas as simulações apresentadas aqui, consideramos que o problema imiscível é bidimensional, sendo, portanto, o meio poroso caracterizado por campos de permeabilidade e de porosidade definidos em regiões Euclideanas. No entanto, a teoria abordada neste trabalho não impõe restrições para sua aplicação aos problemas tridimensionais.

Jato transversal de gotas: simulações por ALE/FEM e efeitos interfaciais.

Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos.

A bacia hidrográfica Lagos-São João pós década de 1960: um estudo das transformações espaciais da região e suas influências sobre a qualidade ambiental

A Bacia Hidrográfica Lagos São João, localizada no sudoeste do Estado do Rio de Janeiro, abrange 13 municípios que abrigam cerca de 520 mil habitantes. Na temporada de férias esse número sobe para mais de 1 milhão de pessoas. A pastagem constitui o principal tipo de uso do solo, em seguida vem as áreas urbanas e as salinas. A partir da década de 1960 essa região passou a receber maior contingente populacional, tanto de veranistas quanto de moradores fixos, beneficiados pela implantação de novas vias de acesso, como a Ponte Rio-Niteroi e pela construção da represa de Juturnaíba, que ampliou o abastecimento de água dos municípios. Surge neste contexto a especulação imobiliária, que acelera a ocupação das terras próximas a Lagoa de Araruama. Rapidamente essas terras foram loteadas e o setor da construção civil foi ganhando força. Entretanto, a região não contou com adequado planejamento, e os investimentos em saneamento básico e outras infraestruturas urbanas não acompanharam o ritmo da construção civil, que cada vez mais investia em casas, prédios e condomínios, que ampliaram consideravelmente a área urbana e a ocupação da zona costeira. Sendo assim, ficou visível o aumento da malha urbana e a ocupação de áreas impróprias, como as margens dos corpos hídricos, os manguezais, dunas e restingas, além da redução da cobertura vegetal. Dessa forma, foi substancial a perda de qualidade ambiental na região, sobretudo, com relação a água da lagoa e dos rios, que passaram a receber maior volume de efluentes sem tratamento. O potencial turístico da região tem sido explorado e provocado altos investimentos dos agentes de especulação imobiliária, entretanto além de promover a ocupação em áreas irregulares, leva a privatização de espaços públicos e incentiva o fenômeno da segunda residência. A chegada de novos turistas iniciou o processo de desenvolvimento do turismo e, consequentemente, a redução da produção salineira. Com isso, o espaço local ganhou novos significados, inseridos pela lógica da urbanização turística. Foi essa nova lógica transformadora que, gradativamente, valorizou a paisagem local, ampliando e encarecendo o seu consumo. Além de ampliar as transformações espaciais, tendo em vista a expansão da malha urbana verificada nas imagens de satélites, atuais, que foram comparadas com fotografias aéreas de décadas anteriores. Todas as transformações ocorridas na região apresentam alguma relação com o desencadeamento de novos problemas ambientais identificados nos seus ecossistemas, sobretudo a Lagoa de Araruama, ou a ampliação de problemas anteriormente existentes.

Comportamento de uma área de resíduos de bauxita durante a operação de enchimento.

A técnica de deposição de resíduos de processamento de bauxita prevê uma operação inicial de lançamento de resíduos, por via úmida, no interior de reservatórios artificiais. Nesta fase de deposição, os resíduos são lançados em forma de polpa, sofrendo processos de sedimentação e adensamento por peso próprio. A vida útil dos reservatórios depende fundamentalmente do comportamento geotécnico do resíduo, cujas propriedades variam ao longo do tempo e da profundidade. A presente dissertação tem como objetivo a análise do comportamento de uma área de resíduos de bauxita durante a operação de enchimento do reservatório. Para tanto, foi utilizado o programa Plaxis, que executa análises tensão-deformação pelo método dos elementos finitos. Os parâmetros geotécnicos do resíduo foram definidos com base em ensaios de adensamento e piezocone, executados em amostras representativas da polpa de resíduos. Para calibração do modelo, a batimetria executada no reservatório foi confrontada com as cotas de enchimento previstas numericamente. Finalmente, a distribuição dos parâmetros geotécnicos (peso específico, índice de vazios, permeabilidade) com a profundidade prevista numericamente foi confrontada com os valores medidos no campo, ao final da etapa de enchimento. O conhecimento da variação dos parâmetros geotécnicos com a profundidade ao final do enchimento é importante para avaliar o comportamento do reservatório na etapa de operação de alteamento, e para a estimativa do ganho de vida útil do reservatório. Os resultados permitiram avaliar o funcionamento da drenagem de fundo, a distribuição complexa dos excessos de poropressão ao final do enchimento, os perfis dos índices físicos com a profundidade e o ganho de vida útil. As análises mostraram que o programa foi capaz de reproduzir os processos simultâneos de sedimentação e adensamento por peso próprio, em análises de enchimento de reservatórios de resíduos.