Uso do modelo de simulação 'ABC 6' para avaliação de áreas sujeitas à ocorrência de ondas de cheia e corridas de massa/detritos: estudo de caso na bacia do Ribeirão da Prata (SP)

Pipelines are linear engineering works, designed mostly for transporting oil and its derivatives for long distances, furnishing even the farthermost zones of the country. Due to oil sector needs to ensure for the safety and conservation of its properties, several geotechnical studies are being held at the pipelines field, in order to preserve this important transportation, and also to prevent accidents, which might seriously compromise the environment and the population who lives around it. The OSBRA pipeline, who connects the city of Paulínia to the capital Brasília, is one of these engineering works that deserves to be pointed out. This research, performed at the Ribeirão da Prata Basin, was a pilot study conducted with the main objective of testing the current methodology efficiency, for future applications in the closest watersheds to the OSBRA pipeline. The objective of this research is to analyze flood wave and debris flow processes in a non-fictional watershed, by comparing two different kinds of methods: the first one based on simulation models (software ‘ABC 6’), and the other one by flood wave and debris flow susceptibility mapping. The results from the hydrological modeling were both hydrographs and ietographs that estimated values of outputs and infiltration. To construct the susceptibility maps were necessary three other maps: ground use and occupation maps, divided according to the different protection degrees that were offered to the ground; maps of dam locations in the area and physiographic compartimentation maps, divided according to the local geology. To complete the methodology, the results were collected from both methods for comparison. The obtained product for this methodology was series of data whose different susceptibility degrees to flood wave and debris flow could define the safest route for a pipeline crossing in this watershed...

Estudo de materiais supercondutores em forma de SQUID com uma constrição usando métodos de paralelização computacional

Pós-graduação em Ciência e Tecnologia de Materiais - FC

Simulação do transporte de poeiras radioactivas resultantes de um acidente no complexo de Angra dos Reis, Brasil

O objetivo da avaliação de impactos ambientais (AIA) é permitir uma análise integrada de possíveis impactos diretos ou indiretos ao meio ambiente decorrentes da implantação e operação de empreendimentos, de forma a propor de medidas ou programas que visem evitar, mitigar ou compensar tais impactos. Para tanto é necessário conhecer as diversas características das áreas direta e indiretamente afetadas pela instalação de um projeto, tais como as condições meteorológicas e climatológicas. Estas também são relevantes no estudo das emissões em cenários de operação regular ou acidental de empreendimentos, dada sua influência nas condições de transporte e de dispersão de poluentes na atmosfera. Neste trabalho é realizado um estudo das condições de dispersão de poluentes na atmosfera para a região da Central Nuclear Almirante Álvaro Alberto (CNAAA) em Angra dos Reis, no Estado do Rio de Janeiro, utilizando o modelo WRF, considerando um cenário acidental com liberações por 48 horas. Os dois episódios simulados representam os regimes de tempo predominantes na região obtidos a partir da análise pelo o método k-means sobre as EOFs para o campo de pressões ao nível médio do mar entre os anos de 1985 e 2014. A aplicação da metodologia dos regimes de tempo permite observar os fenômenos meteorológicos de grande escala persistentes e recorrentes sobre uma dada região, servindo como uma ferramenta para a elaboração de estudos e documentos técnicos que fundamentem a decisão dos órgãos reguladores.

Efecto de diferentes métodos de almacenamiento sobre la viabilidad de la semilla de soya (Glycine Max (L) Merr), variedad "cristalina"

Se sometió a prueba de almacenamiento semilla de soya variedad “cristalina” en bodegas con distinto ambiente, embaladas en diferentes empaques con los objetivos de determinar las posibles causas que deterioran la semilla de soya en condiciones de almacenamiento natural y artificial existentes en nuestro país, comparar tres ambientes de bodegas y tres empaques en la conservación de la viabilidad de la semilla y determinar el periodo de almacenamiento que permita mantener la semilla agronómicamente aceptable durante el periodo de almacenamiento se registraron los factores ambientales y se realizaron análisis de semilla en el laboratorio. Se determinó que las condiciones ambientales de alta temperatura en la bodega Nº3 /EMPROSEM) con promedio de 32.5ºC y alta humedad relativa en bodega Nº1 (ESTELI) con promedio de 77% y la lata humedad inicial de la semilla en todas las bodegas , (13.9%) fueron los factores que provocaron el rápido deterioro de la semilla por ataque principalmente de hongos de almacén además se logró determinar Nº2 (EMPROSEN) con ambiente artificial que la semilla se conserva por encima del 80% de germinación por más tiempo (pero no lo que se necesita de un ciclo a otro). También se observó que el empaque de papel kraff asfaltado conserva mejor la calidad de la semilla por otra parte se comprobó que esta semilla es susceptible al deterioro en almacenamiento natural, de aquí que se sugiere continuar los estudios sobre este tema.

euskadiMoving: sistema de seguimiento de transporte público

El tema abordado por este proyecto es la información disponible acerca de los horarios del transporte público en las ciudades, concretamente centrado en la ciudad de Donostia - San Sebastián. En esta ciudad, la información ofrecida a los usuarios está centrada en: Tiempo que falta para que lleguen los autobuses. Únicamente accesible en las marque- sinas. Horarios con la primera hora de salida y la frecuencia de la línea. Estos métodos de información se pueden mejorar mediante un sistema de información en tiempo real, el cual informe de la posición actual de cada autobús de cada línea a petición del usuario. Asimismo, siendo esta información accesible a través de un servicio web, se evita el tener que desplazarse a las marquesinas para informarse. Por tanto, estos son los objetivos que se consideran para el proyecto. La realización de un servicio web que muestre a petición del usuario las líneas de transporte público disponibles y este pueda seleccionar cuales visualizar, de manera que al seleccionar una línea se muestren todos los vehículos de dicha línea y su posición actual.

Estudo dos fenômenos de transporte em cromatografia através da aplicação de modelos de rede e estocásticos

O estudo dos diferentes fenômenos de separação tem sido cada vez mais importante para os diferentes ramos da indústria e ciência. Devido à grande capacidade computacional atual, é possível modelar e analisar os fenômenos cromatográficos a nível microscópico. Os modelos de rede vêm sendo cada vez mais utilizados, para representar processos de separação por cromatografia, pois através destes pode-se representar os aspectos topológicos e morfológicos dos diferentes materiais adsorventes disponíveis no mercado. Neste trabalho visamos o desenvolvimento de um modelo de rede tridimensional para representação de uma coluna cromatográfica, a nível microscópico, onde serão modelados os fenômenos de adsorção, dessorção e dispersão axial através de um método estocástico. Também foram utilizadas diferentes abordagens com relação ao impedimento estérico Os resultados obtidos foram comparados a resultados experimentais. Depois é utilizado um modelo de rede bidimensional para representar um sistema de adsorção do tipo batelada, mantendo-se a modelagem dos fenômenos de adsorção e dessorção, e comparados a sistemas reais posteriormente. Em ambos os sistemas modelados foram analisada as constantes de equilíbrio, parâmetro fundamental nos sistemas de adsorção, e por fim foram obtidas e analisadas isotermas de adsorção. Foi possível concluir que, para os modelos de rede, os fenômenos de adsorção e dessorção bastam para obter perfis de saída similares aos vistos experimentalmente, e que o fenômeno da dispersão axial influência menos que os fenômenos cinéticos em questão

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 método espectronodal livre de erros de truncamento espacial para problemas adjuntos de transporte de partículas neutras monoenergéticas na formulação de ordenadas discretas em geometria unidimensional

Um método numérico nodal livre de erros de truncamento espacial é desenvolvido para problemas adjuntos de transporte de partículas neutras monoenergéticas em geometria unidimensional com fonte fixa na formulação de ordenadas discretas (SN). As incógnitas no método são os fluxos angulares adjuntos médios nos nodos e os fluxos angulares adjuntos nas fronteiras dos nodos, e os valores numéricos gerados para essas quantidades são os obtidos a partir da solução analítica das equações SN adjuntas. O método é fundamentado no uso da convencional equação adjunta SN discretizada de balanço espacial, que é válida para cada nodo de discretização espacial e para cada direção discreta da quadratura angular, e de uma equação auxiliar adjunta não convencional, que contém uma função de Green para os fluxos angulares adjuntos médios nos nodos em termos dos fluxos angulares adjuntos emergentes das fronteiras dos nodos e da fonte adjunta interior. Resultados numéricos são fornecidos para ilustrarem a precisão do método proposto.

Problemas inversos aplicados à identificação de parâmetros hidrodinâmicos de um modelo do estuário do rio Macaé

Esta tese tem por objetivo propor uma estratégia de obtenção automática de parâmetros hidrodinâmicos e de transporte através da solução de problemas inversos. A obtenção dos parâmetros de um modelo físico representa um dos principais problemas em sua calibração, e isso se deve em grande parte à dificuldade na medição em campo desses parâmetros. Em particular na modelagem de rios e estuários, a altura da rugosidade e o coeficiente de difusão turbulenta representam dois dos parâmetros com maior dificuldade de medição. Nesta tese é apresentada uma técnica automatizada de estimação desses parâmetros através deum problema inverso aplicado a um modelo do estuário do rio Macaé, localizado no norte do Rio de Janeiro. Para este estudo foi utilizada a plataforma MOHID, desenvolvida na Universidade Técnica de Lisboa, e que tem tido ampla aplicação na simulação de corpos hídricos. Foi realizada uma análise de sensibilidade das respostas do modelo com relação aos parâmetros de interesse. Verificou-se que a salinidade é uma variável sensível a ambos parâmetros. O problema inverso foi então resolvido utilizando vários métodos de otimização através do acoplamento da plataforma MOHID a códigos de otimização implementados em Fortran. O acoplamento foi realizado de forma a não alterar o código fonte do MOHID, possibilitando a utilização da ferramenta computacional aqui desenvolvida em qualquer versão dessa plataforma, bem como a sua alteração para o uso com outros simuladores. Os testes realizados confirmam a eficiência da técnica e apontam as melhores abordagens para uma rápida e precisa estimação dos parâmetros.

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.

Identificação de areias siliciclásticas para a recuperação de praias em erosão através de métodos geofísicos acústicos

Devido a sua grande extensão latitudinal, a costa brasileira é influenciada por diferentes regimes climáticos e oceanográficos. Adicionalmente, a distribuição da população brasileira é caracterizada por uma alta concentração nas capitais litorâneas. Todos esses fatores levam à construção de inúmeras estruturas de engenharia que podem de alguma forma impactar o transporte de sedimento e consequentemente o balanço sedimentar de algumas praias. Uma das formas mais eficientes de recuperar esse balanço sedimentar é a alimentação artificial de praias, com sedimentos provenientes da plataforma continental com características semelhantes. Os métodos geofísicos acústicos permitem mapear de forma eficiente o fundo e o subfundo marinho para a busca de áreas fontes para a lavra de forma eficiente. O objetivo deste estudo é a identificação de padrões geoacústicos a partir de analises quantitativas e qualitativas, para a caracterização de áreas fontes de areias siliciclásticas compatíveis com sedimentos de praias em erosão na cidade do Rio de Janeiro. O mapeamento foi realizado na plataforma continental interna do Rio de Janeiro, em área adjacente à praia de Itaipuaçu (Maricá, RJ). As análises evidenciaram depósitos, localizados entre 19 e 30 metros, de areias com granulometria, textura e selecionamento em condições ideias para a recomposição de praias do Rio de Janeiro.

Método numérico de Matriz Resposta acoplado a um esquema de reconstrução espacial analítica para cálculos unidimensionais de transporte de nêutrons na formulação de ordenadas discretas multigrupo de energia com fonte fixa

Um método de Matriz Resposta (MR) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas multigrupo de transporte de nêutrons com fonte fixa e em geometria unidimensional na formulação de ordenadas discretas (SN). Portanto, o método multigrupo MR com esquema iterativo de inversão nodal parcial (NBI) converge valores numéricos para os fluxos angulares nas fronteiras das regiões que coincidem com os valores da solução analítica das equações multigrupo SN, afora os erros de arredondamento da aritmética finita computacional. É também desenvolvido um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em cada grupo de energia em um intervalo qualquer do domínio definido pelo usuário, com um passo de avanço também escolhido pelo usuário. Resultados numéricos são apresentados para ilustrar a precisão do presente método em cálculos de malha grossa.

Uma abordagem matricial para modelagem e simulação de redes de trocadores de calor com aplicações para o gerenciamento da deposição

Uma rede de trocadores de calor pode ser definida como um grupo de trocadores de calor interligados com o objetivo de reduzir a necessidade de energia de um sistema, sendo largamente usada nas indústrias de processos. Entretanto, uma rede está sujeita à deposição, a qual causa um decréscimo na efetividade térmica dos trocadores. Este fenômeno é provocado pelo acúmulo de materiais indesejáveis sobre a superfície de troca térmica. Para compensar a redução de efetividade térmica causada pela deposição, torna-se necessário um aumento no consumo de utilidades. Isto eleva os custos de operação, assim como os custos de manutenção. Estima-se que os custos associados à deposição atinjam bilhões de dólares anualmente. Em face a este problema, vários trabalhos de pesquisa têm investigado métodos para prevenir a deposição e/ou gerenciar as operações em uma rede. Estudos envolvem desde a otimização de trocadores de calor individuais, simulação e monitoramento de redes, até a otimização da programação das paradas para limpeza de trocadores de calor em uma rede. O presente trabalho apresenta a proposição de um modelo para simulação de redes de trocadores de calor com aplicações no gerenciamento da deposição. Como conseqüência, foi desenvolvido um conjunto de códigos computacionais integrados, envolvendo a simulação estacionária de redes, a simulação pseudo-estacionária do comportamento de redes em relação à evolução da deposição, a estimação de parâmetros para diagnóstico do problema da deposição e a otimização operacional deste tipo de sistema. Com relação ao simulador estacionário, o modelo da rede foi formulado matricialmente e os balanços de massa e energia são resolvidos como sistemas de equações lineares. Do ponto de vista da otimização, o procedimento proposto redistribui as vazões, visando um melhor aproveitamento térmico dos trocadores da rede, como, por exemplo, buscando as vazões da rede que maximizem a temperatura da corrente de entrada no forno em unidades de destilação atmosférica de óleo cru. Os algoritmos foram implementados em alguns exemplos da literatura e em um problema de uma refinaria real. Os resultados foram promissores, o que sugere que a proposta deste trabalho pode vir a ser uma abordagem interessante para operações envolvendo redes de trocadores de calor

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.

Implementação paralela de um código de elementos finitos em 2D para as Equações de Navier-Stokes para fluidos incompressíveis com transporte de escalares.

O estudo do fluxo de água e do transporte escalar em reservatórios hidrelétricos é importante para a determinação da qualidade da água durante as fases iniciais do enchimento e durante a vida útil do reservatório. Neste contexto, um código de elementos finitos paralelo 2D foi implementado para resolver as equações de Navier-Stokes para fluido incompressível acopladas a transporte escalar, utilizando o modelo de programação de troca de mensagens, a fim de realizar simulações em um ambiente de cluster de computadores. A discretização espacial é baseada no elemento MINI, que satisfaz as condições de Babuska-Brezzi (BB), que permite uma formulação mista estável. Todas as estruturas de dados distribuídos necessárias nas diferentes fases do código, como pré-processamento, solução e pós-processamento, foram implementadas usando a biblioteca PETSc. Os sistemas lineares resultantes foram resolvidos usando o método da projeção discreto com fatoração LU por blocos. Para aumentar o desempenho paralelo na solução dos sistemas lineares, foi empregado o método de condensação estática para resolver a velocidade intermediária nos vértices e no centróide do elemento MINI separadamente. Os resultados de desempenho do método de condensação estática com a abordagem da solução do sistema completo foram comparados. Os testes mostraram que o método de condensação estática apresenta melhor desempenho para grandes problemas, às custas de maior uso de memória. O desempenho de outras partes do código também são apresentados.