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.

Identificação de danos estruturais utilizando dados no domínio do tempo provenientes de ensaios de vibração

O presente trabalho aborda o problema de identificação de danos em uma estrutura a partir de sua resposta impulsiva. No modelo adotado, a integridade estrutural é continuamente descrita por um parâmetro de coesão. Sendo assim, o Modelo de Elementos Finitos (MEF) é utilizado para discretizar tanto o campo de deslocamentos, quanto o campo de coesão. O problema de identificação de danos é, então, definido como um problema de otimização, cujo objetivo é minimizar, em relação a um vetor de parâmetros nodais de coesão, um funcional definido a partir da diferença entre a resposta impulsiva experimental e a correspondente resposta prevista por um MEF da estrutura. A identificação de danos estruturais baseadas no domínio do tempo apresenta como vantagens a aplicabilidade em sistemas lineares e/ou com elevados níveis de amortecimento, além de apresentar uma elevada sensibilidade à presença de pequenos danos. Estudos numéricos foram realizados considerando-se um modelo de viga de Euler-Bernoulli simplesmente apoiada. Para a determinação do posicionamento ótimo do sensor de deslocamento e do número de pontos da resposta impulsiva, a serem utilizados no processo de identificação de danos, foi considerado o Projeto Ótimo de Experimentos. A posição do sensor e o número de pontos foram determinados segundo o critério D-ótimo. Outros critérios complementares foram também analisados. Uma análise da sensibilidade foi realizada com o intuito de identificar as regiões da estrutura onde a resposta é mais sensível à presença de um dano em um estágio inicial. Para a resolução do problema inverso de identificação de danos foram considerados os métodos de otimização Evolução Diferencial e Levenberg-Marquardt. Simulações numéricas, considerando-se dados corrompidos com ruído aditivo, foram realizadas com o intuito de avaliar a potencialidade da metodologia de identificação de danos, assim como a influência da posição do sensor e do número de dados considerados no processo de identificação. Com os resultados obtidos, percebe-se que o Projeto Ótimo de Experimentos é de fundamental importância para a identificação de danos.

Simulador MATLAB para Sistemas Ópticos WDM Amplificados.

Neste trabalho é apresentado um simulador MATLAB para sistemas ópticos WDM amplificados baseado na solução das equações não lineares de Schrödinger acopladas, pelo método de Fourier de passo alternado. Este simulador permite o estudo da propagação de pulsos em fibras ópticas, considerando dispersão cromática, efeitos não lineares como automodulação de fase e modulação de fase cruzada e atenuação, prevendo também o emprego de amplificadores ópticos a fibra dopada com Érbio (EDFAs). Através de simulações numéricas, foi explorada a técnica de otimização do posicionamento de um EDFA ao longo de um enlace óptico, sem repetidores, que objetiva a redução dos custos de implantação de sistemas ópticos, seja pela diminuição da potência do transmissor ou pela relaxação da exigência de sensibilidade do receptor. Além disto, pode favorecer um aumento na capacidade do sistema, através do aumento do alcance ou da taxa de transmissão. A concordância dos resultados obtidos com os disponíveis na literatura confirmam a validade da técnica, bem como a versatilidade e robustez do simulador desenvolvido.

Simulações de problemas inversos com aplicações em engenharia nuclear usando técnicas de transporte de partículas neutras monoenergéticas na formulação unidimensional de ordenadas discretas

Neste trabalho, três técnicas para resolver numericamente problemas inversos de transporte de partículas neutras a uma velocidade para aplicações em engenharia nuclear são desenvolvidas. É fato conhecido que problemas diretos estacionários e monoenergéticos de transporte são caracterizados por estimar o fluxo de partículas como uma função-distribuição das variáveis independentes de espaço e de direção de movimento, quando os parâmetros materiais (seções de choque macroscópicas), a geometria, e o fluxo incidente nos contornos do domínio (condições de contorno), bem como a distribuição de fonte interior são conhecidos. Por outro lado, problemas inversos, neste trabalho, buscam estimativas para o fluxo incidente no contorno, ou a fonte interior, ou frações vazio em barras homogêneas. O modelo matemático usado tanto para os problemas diretos como para os problemas inversos é a equação de transporte independente do tempo, a uma velocidade, em geometria unidimensional e com o espalhamento linearmente anisotrópico na formulação de ordenadas discretas (SN). Nos problemas inversos de valor de contorno, dado o fluxo emergente em um extremo da barra, medido por um detector de nêutrons, por exemplo, buscamos uma estimativa precisa para o fluxo incidente no extremo oposto. Por outro lado, nos problemas inversos SN de fonte interior, buscamos uma estimativa precisa para a fonte armazenada no interior do domínio para fins de blindagem, sendo dado o fluxo emergente no contorno da barra. Além disso, nos problemas inversos SN de fração de vazio, dado o fluxo emergente em uma fronteira da barra devido ao fluxo incidente prescrito no extremo oposto, procuramos por uma estimativa precisa da fração de vazio no interior da barra, no contexto de ensaios não-destrutivos para aplicações na indústria. O código computacional desenvolvido neste trabalho apresenta o método espectronodal de malha grossa spectral Greens function (SGF) para os problemas diretos SN em geometria unidimensional para gerar soluções numéricas precisas para os três problemas inversos SN descritos acima. Para os problemas inversos SN de valor de contorno e de fonte interior, usamos a propriedade da proporcionalidade da fuga de partículas; ademais, para os problemas inversos SN de fração de vazio, oferecemos a técnica a qual nos referimos como o método físico da bissecção. Apresentamos resultados numéricos para ilustrar a precisão das três técnicas, conforme descrito nesta tese.

Análise matemática de soluções descontínuas de leis de conservação hiperbólicas e resoluções numéricas para a captura de ondas de choque em escoamentos multifásicos em meios porosos

O processo de recuperação secundária de petróleo é comumente realizado com a injeção de água ou gás no reservatório a fim de manter a pressão necessária para sua extração. Para que o investimento seja viável, os gastos com a extração precisam ser menores do que o retorno financeiro obtido com a produção de petróleo. Objetivando-se estudar possíveis cenários para o processo de exploração, costuma-se utilizar simulações dos processos de extração. As equações que modelam esse processo de recuperação são de caráter hiperbólico e não lineares, as quais podem ser interpretadas como Leis de Conservação, cujas soluções são complexas por suas naturezas descontínuas. Essas descontinuidades ou saltos são conhecidas como ondas de choque. Neste trabalho foi abordada uma análise matemática para os fenômenos oriundos de leis de conservação, para em seguida utilizá-la no entendimento do referido problema. Foram estudadas soluções fracas que, fisicamente, podem ser interpretadas como ondas de choque ou rarefação, então, para que fossem distinguidas as fisicamente admissíveis, foi considerado o princípio de entropia, nas suas diversas formas. As simulações foram realizadas nos âmbitos dos escoamentos bifásicos e trifásicos, em que os fluidos são imiscíveis e os efeitos gravitacionais e difusivos, devido à pressão capilar, foram desprezados. Inicialmente, foi feito um estudo comparativo de resoluções numéricas na captura de ondas de choque em escoamento bifásico água-óleo. Neste estudo destacam-se o método Composto LWLF-k, o esquema NonStandard e a introdução da nova função de renormalização para o esquema NonStandard, onde obteve resultados satisfatórios, principalmente em regiões onde a viscosidade do óleo é muito maior do que a viscosidade da água. No escoamento bidimensional, um novo método é proposto, partindo de uma generalização do esquema NonStandard unidimensional. Por fim, é feita uma adaptação dos métodos LWLF-4 e NonStandard para a simulação em escoamentos trifásicos em domínios unidimensional. O esquema NonStandard foi considerado mais eficiente nos problemas abordados, uma vez que sua versão bidimensional mostrou-se satisfatória na captura de ondas de choque em escoamentos bifásicos em meios porosos.

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.

Toxicidade aguda sobre Daphnia similis decorrente de contaminação aquática recente e antiga por gasolina, etanol e suas misturas: simulações em escala de bancada.

A contaminação de ambientes aquáticos decorrente de acidentes com gasolina, álcool combustível e misturas binárias representa um risco crescente, tendo em vista as projeções do setor para os próximos 50 anos. O objetivo do presente estudo foi avaliar a toxicidade aguda da Gasolina C, Gasolina P e álcool combustível isoladamente e em misturas binárias, assim como de suas respectivas Frações Solúveis em Água (FSA) e Frações Dispersas em Água (FDA) sobre Daphnia similis. O estudo ainda incluiu a avaliação da toxicidade aguda remanescente na matriz água de uma contaminação antiga (intemperismo) com a Gasolina C. Paralelamente, foram conduzidos ensaios de toxicidade aguda com amostras ambientais (água subterrânea, superficial e elutriato a partir de sedimentos) de uma área alagada com histórico de contaminação antiga. O cultivo e os ensaios com D. similis foram de acordo com a NBR 12.713 (2009). Tanto a gasolina C quanto a P foram extremamente tóxicas para os organismos, apresentando valores médios de CE50% em 48 h de 0,00113% e 0,058% respectivamente. As diferenças entre os resultados obtidos com a Gasolina C e aqueles obtidos com suas frações FSA e FDA foram significativas (p < 0,05), sendo que não houve diferença significativa entre a toxicidade aguda da FSA e da FDA (p < 0,05). Os resultados obtidos com os ensaios com Gasolina P e FDA não apresentaram diferenças significativas entre si (p < 0,05), mas, foram significativamente diferentes daqueles obtidos com FSA (p < 0,05). Os resultados dos ensaios de toxicidade aguda com misturas binárias sugeriram efeito menos que aditivo (antagonismo). Os resultados da simulação de uma contaminação antiga demonstraram redução acentuada da toxicidade para D. similis ao longo de apenas 28 dias. Entretanto, com relação aos ensaios com as amostras ambientais da área com histórico de contaminação, apesar da ausência ou baixa toxicidade nas amostras de água superficial (sugerindo intemperismo), toxicidade alta foi observada em amostras de água subterrânea e no elutriato de sedimentos, sugerindo condições de adsorção aos sedimentos com alto teor de argila e/ou aprisionamento dos compostos em zona saturada.

Análise tensão deformação de uma barragem de concreto em solo residual preponderantemente anisotrópico.

Estudos de fundações de obras de arte, como barragens de concreto-gravidade e barragens de terra, devem contemplar todos os esforços atuantes no maciço de fundação, principalmente as tensões e as deformações esperadas durante todo o processo construtivo e no período pós-construtivo. Quando essas estruturas são apoiadas sobre rochas de boa sanidade, a escolha do barramento geralmente favorece a implantação de barragens de concreto. Entretanto, quando os maciços de fundação são formados por solos, a opção técnica geralmente mais bem aceita é quanto à utilização de barragens de terra. Em ambos os casos, as análises de estabilidade e de deformação são desenvolvidas por métodos analíticos bem consolidados na prática. Nas condições mais adversas de fundação, seja em rochas ou em solos com marcante anisotropia e estruturas reliquiares herdadas da rocha-mãe, a utilização de modelos constitutivos anisotrópicos em análises por elementos finitos propicia simulações mais realistas dessas feições estruturais, contribuindo para o seu melhor conhecimento. O presente trabalho teve por objetivo demonstrar a utilização de um modelo constitutivo anisotrópico no estudo da fundação da Barragem San Juan, localizada na República Dominicana, que foi concebida como uma estrutura tipo concreto-gravidade apoiada sobre solos residuais jovens altamente anisotrópicos. Nessa obra, apesar de sua pequena altura, a presença marcante de descontinuidades ensejou um estudo mais detalhado do comportamento tensão-deformação da fundação, levando em conta o levantamento detalhado da atitude das descontinuidades presentes no maciço e os resultados de ensaios de resistência em planos paralelos e normais às descontinuidades. Para a estimativa de deformações, os respectivos módulos de Young (Es) foram estimados com base em correlações com a resistência à penetração SPT desenvolvidas nesta dissertação, a partir de um estudo estatístico baseado em várias publicações disponíveis na literatura. As análises numéricas por elementos finitos foram desenvolvidas através do programa Plaxis 2D, utilizando-se como modelo constitutivo aquele denominado The Jointed Rock Model, que é particularmente recomendado para análises de estabilidade e deformação de materiais anisotrópicos. Os resultados das análises numéricas foram comparados com as análises de equilíbrio limite elaboradas para o projeto executivo da referida barragem, pelo programa Slope, utilizando o método rigoroso de Morgenstern e Price, que se mostrou conservador. Os resultados das análises numéricas mostraram sua inequívoca versatilidade para a escolha de opções de reforço da fundação, através de dentes que objetivavam o aumento das condições de estabilidade da barragem.

Avaliações econômicas do uso da vacina contra o Papilomavírus Humano (HPV) em meninas adolescentes: uma revisão sistemática

O câncer de colo do útero persiste como um importante problema de saúde em todo o mundo, em particular nos países em desenvolvimento. Duas vacinas contra o papilomavirus humano (HPV) encontram-se atualmente disponíveis e aprovadas para uso em meninas adolescentes, antes do início da vida sexual: uma bivalente, contra os sorotipos 16 e 18 e outra quadrivalente, contra os sorotipos 6, 11, 16 e 18. Estes imunobiológicos têm por objetivo induzir uma imunidade contra o papilomavírus e, desta forma, atuar na prevenção primária do câncer do colo de útero. As avaliações econômicas podem ser um elemento que auxiliem nos processos de tomada de decisão sobre a incorporação da vacina em programas de imunização nacionais. Estas avaliações foram o objeto central deste trabalho, que teve como objetivo sintetizar as evidências procedentes de uma revisão sistemática da literatura de estudos de avaliação econômica da utilização da vacina contra o HPV em meninas adolescentes e pré-adolescentes. Foi realizada uma busca na literatura nas bases MEDLINE (via Pubmed), LILACS (via Bireme) e National Health Service Economic Evaluation Database (NHS EED) ate junho de 2010. Dois avaliadores, de forma independente, selecionaram estudos de avaliação econômica completa, que tivessem como foco a imunização para HPV em mulheres com as vacinas comercialmente disponíveis direcionada à população adolescente. Após a busca, 188 títulos foram identificados; destes, 39 estudos preencheram os critérios de elegibilidade e foram incluídos na revisão. Por tratar-se de uma revisão de avaliações econômicas, não foi realizada uma medida de síntese dos valores de relação incremental entre custos e efetividade. Os 39 artigos incluídos envolveram 51 avaliações econômicas em 26 países. Predominaram estudos de custo-utilidade (51%). Do ponto de vista da perspectiva da análise, predominou o dos sistemas de saúde (76,4%). A maioria dos trabalhos (94,9%) elegeu meninas, com idade entre 9 e 12 anos, como sua população alvo e desenvolveu simulações considerando imunidade para toda a vida (84,6%). Os modelos utilizados nos estudos foram do tipo Markov em 25 análises, de transmissão dinâmica em 11 e híbridos em 3. As análises de sensibilidade revelaram um conjunto de elementos de incerteza, uma parte significativa dos quais relacionados a aspectos vacinais: custos da vacina, duração da imunidade, necessidade de doses de reforço, eficácia vacinal e cobertura do programa. Estes elementos configuram uma área de especial atenção para futuros modelos que venham a ser desenvolvidos no Brasil para análises econômicas da vacinação contra o HPV.

Estudo do impacto do complexo petroquímico COMPERJ na atmosfera

Os modelos de dispersão de ar têm sido utilizados amplamente para investigação de padrões de dispersão, comportamento de emissões, estimativas de potenciais riscos a saúde humana, elaboração de propostas de gerenciamento ambiental e, também, na previsão de impacto da qualidade do ar. Alguns modelos, homologados por agência de regulação de alguns países ou comunidades, servem de base para as análises de risco com auxílio de simulação. Somente após este tipo de análise uma unidade fabril terá direito à sua instalação e operação nestes países. Nesta dissertação, serão abordados os principais poluentes presentes em uma indústria petroquímica básica, uma revisão sobre os principais tipos de modelos existentes no mercado e um estudo de caso será realizado empregando os modelos AERMOD e OZIPR/SAPRC. A indústria petroquímica básica a servir de modelo será o Complexo Petroquímico do Rio de Janeiro, que deverá ser instalado no município de Itaboraí e com operação prevista para o ano de 2012. De acordo com as simulações realizadas neste trabalho, o poluente NOx apresentou os resultados mais críticos violando em algumas áreas os padrões primários e secundários de emissão. Diante deste fato, o ozônio se tornou um poluente secundário importante a ser analisado. E para sua simulação, premissas tiveram que ser tomadas, devido a ausência de dados, criando cenários que apresentaram resultados díspares: Ora violando os limites, ora se mantendo abaixo deles. Apesar disso, estes cenários apontaram a mesma solução para minimizar esta formação de ozônio: controlar as emissões de compostos orgânicos voláteis.

Um método de matriz resposta com esquema iterativo de inversão parcial por região para problemas unidimensionais de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas

Um método de matriz resposta (RM) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas de transporte de nêutrons monoenergéticos e com fonte fixa, em geometria unidimensional na formulação de ordenadas discretas (SN). O método RM com esquema iterativo de inversão parcial por região (RBI) 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 SN, afora os erros de arredondamento da aritmética finita computacional. Desenvolvemos um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em qualquer ponto 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.

Modelagem do comportamento dinâmico de passarelas tubulares em aço e mistas (aço-concreto)

A experiência dos engenheiros estruturais e os conhecimentos adquiridos pelo uso de materiais e novas tecnologias, têm ocasionado estruturas de aço e mistas (aço-concreto) de passarelas cada vez mais ousadas. Este fato tem gerado estruturas de passarelas esbeltas, e consequentemente, alterando os seus estados de limite de serviço e último associados ao seu projeto. Uma consequência direta desta tendência de projeto é o aumento considerável das vibrações das estruturas. Portanto, a presente investigação foi realizada com base em um modelo de carregamento mais realista, desenvolvido para incorporar os efeitos dinâmicos induzidos pela caminhada de pessoas. O modelo de carregamento considera a subida e a descida da massa efetiva do corpo em cada passo. A posição da carga dinâmica também foi alterada de acordo com a posição do pedestre sobre a estrutura e a função do tempo gerada, possui uma variação espacial e temporal. O efeito do calcanhar do pedestre também foi incorporado na análise. O modelo estrutural investigado baseia-se em uma passarela tubular (aço-concreto), medindo 82,5m. A estrutura é composta por três vãos (32,5 m, 20,0 m e 17,5 m, respectivamente) e dois balanços (7,5 m e 5,0 m, respectivamente). O sistema estrutural é constituído por perfis de aço tubular e uma laje de concreto, e é atualmente utilizada para travessia de pedestres. Esta investigação é realizada com base em resultados experimentais, relacionando a resposta dinâmica da passarela com as obtidas via modelos de elementos finitos. O modelo computacional proposto adota as técnicas de refinamento de malha, usualmente presente em simulações pelo método de elementos finitos. O modelo de elementos finitos foi desenvolvido e validado com resultados experimentais. Este modelo de passarela tubular permitiu uma avaliação dinâmica completa, investigando especialmente ao conforto humano e seus limites de utilização associados à vibração. A resposta dinâmica do sistema, em termos de acelerações de pico, foi obtida e comparada com os valores limites propostos por diversos autores e padrões de projeto. As acelerações de pico encontradas na presente análise indicou que a passarela tubular investigada apresentou problemas relacionados com o conforto humano. Por isso, foi detectado que este tipo de estrutura pode atingir níveis de vibrações excessivas que podem comprometer o conforto do usuário na passarela e especialmente a sua segurança.

Simulação biofisicoquímica em reservatório tropical polimítico

O reservatório do Lobo, localizado no estado de São Paulo, é um sistema dinâmico no qual se desenvolve um ciclo diurno de estratificação e mistura, de modo similar ao que tem sido observado em outros lagos tropicais. Utilizou-se simulação 3D computacional com os softwares ELCOM (Estuary and Lake Computer Model) acoplado ao CAEDYM (Computacional Aquatic Ecosystem Dynamics Model), ambos desenvolvidos pelo CWR (Center for Water Research) da Universidade da Austrália. Foram realizadas cinco simulações: Piloto Primavera baseada em dados reais da estação no ano primavera no reservatório para o ano de 2007; Primavera-P em que as concentrações de fósforo total, fosfato inorgânico e fosfato total dissolvido foram aumentadas em 100% no reservatório (coluna de água e sedimento) e nos rios tributários; Primavera-V na qual a intensidade dos ventos foi aumentada em 50%; Primavera-T onde a temperatura da água (reservatório e tributários) e do ar foram aumentadas em 10C e, Primavera-X, onde a temperatura da água (reservatório e tributários) e do ar sofreu aumento em 10C, as concentrações de fósforo total, fosfato inorgânico e fosfato total dissolvido foram aumentadas em 100% e a velocidade do vento aumentada em 50%. A concentração de clorofila a foi representada pelos grupos cianobactérias e clorofíceas. O espaço de tempo das simulações representou 90 dias. As clorofíceas apresentaram maior desenvolvimento populacional do que as cianobactérias em todas as simulações. No reservatório, a mistura vertical é ocasionada diariamente pelo vento ou por processos convectivos causados pela perda de calor no corpo de água. A oxigenação do reservatório é maior com a ocorrência de ventos e de grupos fotossintéticos. As concentrações totais de fósforo e nitrogênio apresentaram aumento em todas as simulações.

Variação temporal de uma comunidade fitoplanctônica do reservatório de APM-Manso através de modelagem ecológica tridimensional.

A ecologia de reservatórios, que são ecossistemas complexos, dinâmicos e artificiais, vem assumindo destaque no Brasil. O objetivo deste trabalho foi avaliar a viabilidade da aplicação, no reservatório de APM-Manso, de um modelo ecológico tridimensional em estudos sobre a dinâmica fitoplanctônica, simulando a variação temporal do fitoplâncton para cenários distintos de carga de nutrientes. O modelo CAEDYM foi acoplado ao ELCOM e simulação foi realizada em duas etapas: uma hidrodinâmica e outra ecológica. Escolheu-se para as simulações o período de cinco meses, a partir de 1 de setembro de 2005. Foram construídos dois cenários de simulação, o primeiro contendo os valores reais de carga de nutrientes dos principais rios contribuintes medidos em campo, e o segundo com redução na carga nutricional destes rios, simulando um possível processo de substituição de áreas florestadas por áreas de pastagem na bacia do rio Manso. A comunidade fitoplanctônica simulada apresentou rápidas respostas à disponibilidade nutricional do ambiente, e os resultados obtidos corroboraram com diversas teorias sobre as estratégias adaptativas e sobre as dinâmicas algais. Dentre as classes simuladas, Bacillariophyceae e Cryptophyceae se mostraram mais sensíveis às reduções de carga, enquanto Chrolophyceae e Cyanophyceae, apesar de terem suas biomassas reduzidas, sofreram menos com o impacto, sugerindo estarem mais adaptadas à limitação de nutrientes. Os picos chuvosos influenciaram positivamente as taxas de crescimento das Bacillariophyceae apenas no Cenário 1, uma vez que a limitação por nutrientes foi mais decisiva para esta classe no Cenário 2. Observou-se em ambas as simulações uma tendência de substituição na dominância de Cyanophyceae por Chlorophyceae.

Avaliação de correlações para propriedades físico-químicas de frações de petróleo

O petróleo é uma mistura complexa consistindo em um número muito grande de hidrocarbonetos. A descrição completa de todos os hidrocarbonetos existentes nessas misturas é inviável experimentalmente ou consome tempo excessivo em simulações computacionais. Por esta razão, uma abordagem molecular completa para cálculos de propriedades dessas misturas é substituído por uma abordagem pseudo-componente ou por correlações entre propriedades macroscópicas. Algumas dessas propriedades são utilizadas de acordo com a regulamentação de venda de combustíveis, e.g., para gasolina. Dependendo do esquema de refino e do óleo cru utilizado para produção desse combustível, uma larga variedade de valores é encontrada para as propriedades de correntes de processo que compõe o combustível final. A fim de planejar com precisão adequada a mistura dessas correntes, modelos devem estar disponíveis para o cálculo preciso das propriedades necessárias. Neste trabalho, oito séries de combustíveis brasileiros e duas séries de combustíveis estrangeiros foram analisadas: frações de gasolina, querosene, gasóleo e diesel. As propriedades analisadas para as frações são: número de octano, teor de aromáticos, teor de enxofre, índice de refração, densidade, ponto de fulgor, ponto de fluidez, ponto de congelamento, ponto de névoa, ponto de anilina, pressão de vapor Reid e número de cetano. Diversas correlações foram avaliadas e os melhores desempenhos foram destacados, permitindo uma estimação precisa das propriedades do combustível avaliado. Um processo de re-estimação de parâmetros foi aplicado e novos modelos foram ajustados em comparação com os dados experimentais. Esta estratégia permitiu uma estimativa mais exata das propriedades analisadas, sendo verificada por um aumento considerável no desempenho estatístico dos modelos. Além disso, foi apresentado o melhor modelo para cada propriedade e cada série