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.

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.

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

Modelagem e simulação da formação de hidrocarbonetos na combustão do gás natural

A exaustão de um veículo de motor a diesel é uma importante fonte de poluentes atmosféricos, pois forma uma matriz complexa composta de poluentes regulados e não regulados pelos órgãos governamentais. Dentre os poluentes regulados podemos citar óxidos de nitrogênio (NOx) e material particulado. Os poluentes não regulados são pouco estudados até hoje e dentre estes encontra-se a classe dos hidrocarbonetos policíclicos aromáticos e seus derivados nitrados (nitro-HPA). Estes são encontrados na exaustão do diesel na forma gasosa ou agregados ao material particulado. Hoje, o interesse em estudos destes compostos vem aumentando, devido às suas atividades carcinogênicas e mutagênicas às quais estão sujeitas as populações dos centros urbanos. O impacto causado pelos nitro-HPA emitidos por motores a ciclo diesel ao ambiente não está ainda completamente estabelecido. Este estudo consiste na modelagem e simulação do processo de combustão de hidrocarbonetos na faixa de C1 a C4 com o objetivo de descrever a formação de compostos aromáticos, principalmente HPA, e óxidos de nitrogênio a partir de modelos cinéticos de combustão propostos na literatura como referência e fazendo uso do software de simulação Kintecus. Este projeto tem como objetivo em longo prazo propor um modelo cinético para combustão do óleo Diesel. Foi iniciada a construção de um modelo cinético de combustão a partir de modelos de hidrocarbonetos simples de C1 a C4, com formação de aromáticos, HPA e óxidos de nitrogênio. Os modelos originais foram avaliados e modificados a fim de estudar como parâmetros do modelo afetam a concentração das espécies de interesse. Foi observado a tendência de formação de benzeno e fulveno em baixas temperaturas e a tendência de formação de antraceno, pireno, fenantreno a temperaturas mais altas. Foi avaliado que a conversão NO-NO2 ocorre em maiores proporções em reações iniciadas a baixas temperaturas, 600 K. Os resultados indicam que propano é o maior responsável por esta conversão. O modelo final obtido resultou da união dos modelos de combustão Hori e Marinov mais inclusão do GRI-Mech 3.0 e reações adicionais de NOx retiradas da base de dados NIST

Modelagem e simulação da decomposição térmica do óleo mineral isolante aplicadas à classificação de defeitos em transformadores de potência

A análise de gases dissolvidos tem sido aplicada há décadas como a principal técnica de manutenção preditiva para diagnosticar defeitos incipientes em transformadores de potência, tendo em vista que a decomposição do óleo mineral isolante produz gases que permanecem dissolvidos na fase líquida. Entretanto, apesar da importância desta técnica, os métodos de diagnóstico mais conhecidos são baseados em constatações de modelos termodinâmicos e composicionais simplificados para a decomposição térmica do óleo mineral isolante, em conjunto com dados empíricos. Os resultados de simulação obtidos a partir desses modelos não reproduzem satisfatoriamente os dados empíricos. Este trabalho propõe um modelo termodinâmico flexível aprimorado para mimetizar o efeito da cinética de formação de sólidos como restrição ao equilíbrio e seleciona, entre quatro modelos composicionais, aquele que apresenta o melhor desempenho na simulação da decomposição térmica do óleo mineral isolante. Os resultados de simulação obtidos a partir do modelo proposto apresentaram uma melhor adequação a dados empíricos do que aqueles obtidos a partir dos modelos clássicos. O modelo propostofoi, ainda, aplicado ao desenvolvimento de um método de diagnóstico com base fenomenológica.Os desempenhos desta nova proposta fenomenológica e de métodos clássicos de diagnóstico por análise de gases dissolvidos foram comparados e discutidos; o método proposto alcançou desempenho superior a vários métodos usualmente empregados nessa área do conhecimento. E, ainda, um procedimento geral para a aplicação do novo método de diagnóstico é descrito

Modelagem e simulação de padrões arquiteturais de plantas.

2009

Modelagem e simulação da dinâmica do arrastamento por ar

Os compostos orgânicos voláteis constituem uma fonte vulgar de contaminação da água subterrânea, a qual pode ser eliminada pela tecnologia do arrastamento por ar (air stripping) em colunas com enchimento desordenado e utilizando fluxos das fases em contra-corrente. Propõe-se neste trabalho uma nova metodologia de dimensionamento destas colunas, para qualquer tipo de enchimento e de contaminante, onde não há necessidade de se arbitrar nenhum diâmetro, onde se evita o recurso a ábacos experimentais e onde o regime hidráulico conveniente é seleccionado à partida. O procedimento proposto foi algoritmizado e convertido num programa em linguagem C++. Para verificar e testar não só o dimensionamento mas também o comportamento teórico estacionário e dinâmico construiu-se de raiz uma coluna experimental. Seleccionou-se como contaminante uma solução de clorofórmio em água destilada. A experimentação permite, ainda, corrigir o coeficiente de transferência de massa global teórico estimado pelas correlações de Onda e que depende de inúmeros parâmetros nem sempre controláveis experimentalmente. Apresenta-se, em seguida, um modelo original de simulação dinâmica do comportamento da coluna e que é constituído por um sistema de equações diferenciais não lineares (parâmetros distribuidos). No entanto, se os débitos forem arbitrados como constantes, o sistema passa a ser linear apesar de não possuir solução analítica evidente (p.e. por transformações integrais). A discretização por diferenças finitas permitiu superar estas dificuldades. Existe uma notável concordância entre os valores experimentais e os previstos no modelo.

Modelagem e simulação de reatores de polimerização em leito fluidizado

O presente trabalho tem como objetivo o desenvolvimento de um modelo dinâmico genérico para reatores de polimerização em leito fluidizado com capacidade de predizer produção, condições de operação e propriedades do produto, como o índice de fluidez e densidade. Este modelo genérico é então aplicado ao caso específico da polimerização de eteno em reatores do processo UNIPOL, utilizando catalisador de óxido de cromo. De modo a obter um modelo validado, os parâmetros cinéticos são estimados utilizando experimentos realizados em uma planta piloto UNIPOL. Também é apresentado um modelo simplificado em relação ao modelo proposto inicialmente, cujo objetivo é remover as interações dos parâmetros das malhas de controle com os parâmetros cinéticos a serem estimados, aproveitando medidas de planta não utilizadas pelo outro modelo. Junto com a parte de estimação de parâmetros é apresentada a metodologia da análise de sensibilidade para o planejamento de experimentos e estimação de parâmetros.

Modelagem e simulação de um reator tubular de alta pressão para produção de PEBD

Foi realizado o estudo do efeito do tipo de nylon nas blendas PP/PP-MA/Ny. Estas blendas foram caracterizadas por calorimetria diferencial de varredura, microscopia eletrônica de varredura, propriedades mecânicas e reológicas. Na análise reológica foram realizados os testes para os nylons N6, N6,9 e N6,10. Para os nylons N6,6 e N6,12 não foi possível fazer as medidas de viscosidade. Os resultados demonstraram que para os teores de 0,34 e 0,46 g% de MA no PP, o nível de dispersão final segue o fator da reatividade dos grupos terminais amino do nylons. Assim, a ordem de aumento de dispersão para as blendas com os nylons foi de N6,9, N6 e N6,10. Este comportamento foi confirmado na análise morfológica. Para o teor de 0,81 g%, a dispersão foi afetada pelo quociente das viscosidades, obedecendo a ordem esperada para os nylons 6,10, 6,9 e 6 e suas blendas.

Modelagem e simulação dos reatores de polimerização em fase líquida do processo spheripol

O processo Spheripol é, atualmente, um dos mais utilizados na produção industrial de polipropileno e seus copolímeros. É composto por dois reatores tipo loop em série e, em alguns casos, um ou dois reatores tipo fase gás também em série. Neste trabalho é apresentado um modelo para os reatores loop do processo Spheripol. Os reatores loop são modelados como CSTR não ideais. O esquema cinético implementado é relativamente amplo, sendo capaz de englobar copolimerizações multisítio. Uma análise qualitativa do modelo é feita utilizando parâmetros cinéticos presentes na literatura e uma análise quantitativa do modelo é feita com base em dados reais de planta industrial. Ainda com base em dados de planta é feita uma estimação de parâmetros e a validação do modelo. Finalmente, o modelo com parâmetros estimados é aplicado num estudo de caso no qual também é criado um modelo de equilíbrio termodinâmico para a mistura reacional baseado na lei de Henry e na equação de estado de Peng-Robinson.

Modelagem e simulação de um reator trickle-bed para hidrogenação de propeno

Neste trabalho foi modelado um reator catalítico heterogêneo, pressurizado, adiabático e provido de reciclo. O leito deste reator está dividido em duas partes, sendo que a primeira opera como um reator trickle-bed, enquanto a segunda opera como um reator de leito fixo, onde os regentes estão em fase líquida. Um Sistema de Hidrogenação de Propeno em operação no Pólo Petroquímico do Sul foi utilizado para fornecer os dados necessários à validação do modelo proposto. Neste sistema, a reação de hidrogenação de propeno, em fase líquida, ocorre sobre um catalisador níquel/óxido de níquel. Uma equação, que representa a taxa da reação no intervalo de temperatura de interesse para o sistema, foi desenvolvida utilizando-se dados experimentais obtidos em um reator batelada. O comportamento termodinâmico foi representado através de dados de equilíbrio do sistema binário hidrogênio-propano. Considerou-se que a fase líquida escoa em plug-flow e a fase gasosa, devido à existência de um selo líquido no interior do reator, foi considerada estagnada.

Modelagem e simulação da câmara de reação de uma caldeira a carvão pulverizado

O presente trabalho é parte integrante de um projeto multidisciplinar cujo objetivo é construir uma ferramenta computacional capaz de simular a operação de uma Usina Termoelétrica tratando, especificamente, da câmara de reação (local da caldeira onde se desenvolve a queima do combustível). Foi desenvolvido um modelo matemático dinâmico, pseudo-homogêneo, unidimensional e, com propriedades físicas variáveis, com o objetivo de descrever a conversão do calor gerado pela combust em energia térmica, que está baseado em balanços de massa e energia. Este modelo foi implementado no simulador gPROMS e seus parâmetros foram estimados utilizando o módulo gEST deste software (método da máxima verossimilhança), com base em dados de operação da Usina e dados cinéticos de literatura. Para a simulação dinâmica da câmara de reação foi utilizado um método de integração implícita de passos múltiplos com ordem e passo variáveis. A câmara de combustão foi dividida em três zonas, sendo realizados, em cada zona, os balanços de massa e energia, adotando-se a hipótese de mistura perfeita. O modelo desenvolvido prevê o perfil de temperatura ao longo da câmara, as correntes de saída e a composição dos gases de combustão em cada zona; incluindo uma previsão aproximada da concentração dos poluentes formados, principalmente monóxido de nitrogênio e dióxido de enxofre. Os resultados apresentados mostraram que o modelo prediz satisfatoriamente a composição do gás de combustão na saída da câmara de reação e que a divisão da câmara em diferentes zonas de combustão é apropriada para avaliar a geração e o aproveitamento energético na caldeira.

Modelagem e simulação do campo de velocidades do vento em tormentas elétricas

Esta tese propõe um modelo para descrever o campo de velocidade do vento originado numa tormenta elétrica. O escoamento descendente que ocorre durante uma tormenta elétrica madura é modelado como um campo de fluxo axial simétrico. A evolução da velocidade do vento com o tempo é considerada por meio de uma função de modulação empírica. As flutuações da velocidade ao redor da velocidade média em evolução são representadas por um processo aleatório do tipo ruído branco em faixa. Todos os parâmetros utilizados no modelo estão relacionados a variáveis meteorológicas, tais como a altura e diâmetro da nuvem cumulonimbus causadora do evento, a diferença de pressão e a distância entre o solo e a base da nuvem. Adicionalmente considera-se que a nuvem é transportada numa direção pré escolhida pelo vento de fundo, que pode ser devido à circulação geral na região ou a um ciclone extra-tropical. O campo de velocidades do vento resulta da superposição de ambos campos. O modelo é validado comparando as velocidades do vento determinadas em eventos simulados, com registros de estações meteorológicas. Adicionalmente, perfis verticais previstos pelo modelo são comparados com vários perfis para tormentas elétricas propostos na literatura. Finalmente, séries de velocidade máximas anuais a altura padrão de 10m foram simuladas para três estações meteorológicas no sul do Brasil, para as quais dados experimentais são conhecidos, em termos do número médio de tormentas elétricas por ano e das funções distribuição de probabilidade que caracterizam o vento de fundo. A comparação com séries de registros máximos existentes é satisfatória, dando credibilidade adicional ao modelo proposto.