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.

Redução de dimensionalidade em equilíbrio de fases fluidas

Atualmente, existem modelos matemáticos capazes de preverem acuradamente as relações entre propriedades de estado; e esta tarefa é extremamente importante no contexto da Engenharia Química, uma vez que estes modelos podem ser empregados para avaliar a performance de processos químicos. Ademais, eles são de fundamental importância para a simulação de reservatórios de petróleo e processos de separação. Estes modelos são conhecidos como equações de estado, e podem ser usados em problemas de equilíbrios de fases, principalmente em equilíbrios líquido-vapor. Recentemente, um teorema matemático foi formulado (Teorema de Redução), fornecendo as condições para a redução de dimensionalidade de problemas de equilíbrios de fases para misturas multicomponentes descritas por equações de estado cúbicas e regras de mistura e combinação clássicas. Este teorema mostra como para uma classe bem definidade de modelos termodinâmicos (equações de estado cúbicas e regras de mistura clássicas), pode-se reduzir a dimensão de vários problemas de equilíbrios de fases. Este método é muito vantajoso para misturas com muitos componentes, promovendo uma redução significativa no tempo de computação e produzindo resultados acurados. Neste trabalho, apresentamos alguns experimentos numéricos com misturas-testes usando a técnica de redução para obter pressões de ponto de orvalho sob especificação de temperaturas.

Modelagem, simulação e otimização operacional de unidades de síntese de metanol

A indústria de processos químicos tem sofrido consideráveis transformações devido ao acirramento da competitividade. Importantes progressos tecnológicos têm sido atingidos através de técnicas de modelagem, simulação e otimização visando o aumento da lucratividade e melhoria contínua nos processos industriais. Neste contexto, as plantas de metanol, um dos mais importantes produtos petroquímicos, podem ser destacadas. Atualmente, a principal matéria-prima para obtenção de metanol é o gás natural. A produção do metanol é caracterizada por três etapas: geração de gás de síntese, conversão do gás de síntese em metanol (unidade de síntese ou loop de síntese) e purificação do produto na especificação requerida. Os custos fixos e variáveis da unidade de síntese são fortemente dependentes das variáveis operacionais, como temperatura, pressão, razão de reciclo e composição da carga. Desta forma, foi desenvolvido um conjunto de modelos e algoritmos computacionais para representar matematicamente unidades de síntese de metanol. O modelo apresenta operações unitárias associadas aos seguintes equipamentos: divisores de correntes, misturadores de correntes, compressores, trocadores de calor, vasos de flash e reatores. Inicialmente, foi proposto um simulador estacionário, que serviu como base para um pseudo-estacionário, o qual contempla a desativação do catalisador por sinterização térmica. Os simuladores foram criados segundo uma arquitetura seqüencial modular e empregou-se o método de substituição sucessiva para a convergência dos reciclos. O estudo envolveu dois fluxogramas típicos, um constituído por reatores adiabáticos em série, enquanto o outro constituído por um reator tipo quench. Uma análise do efeito das principais variáveis operacionais foi realizada para o reator e para o loop de síntese. Estudou-se também o efeito da desativação do catalisador ao longo do tempo. Uma ferramenta de otimização operacional foi empregada para alcançar a máxima produção manipulando as injeções de carga fria na entrada dos leitos catalíticos. Formulou-se também um problema de maximização do lucro em função da otimização de paradas da unidade para troca do catalisador. Os resultados obtidos apontam que a ferramenta desenvolvida é promissora para a compreensão e otimização da unidade objeto deste estudo

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.

Impacto do uso de biocombustíveis na qualidade do ar da cidade do Rio de Janeiro um estudo de simulação

A frota de veículos movidos por biocombustíveis tem aumentado nos últimos dez anos e com isso as emissões para a atmosfera nas cidades têm sofrido alterações. Devido a esta problemática, estudos com modelos de previsão de emissões são o foco deste trabalho. Alguns modelos, homologados por agências de regulamentação de alguns países ou comunidades, servem de base para as análises de risco com auxílio de simulação.O objetivo geral deste trabalho é avaliar o impacto do uso de biocombustíveis na qualidade do ar na cidade do Rio de Janeiro e estudar cenários de qualidade do ar, em função do aumento da quantidade de biodiesel e álcool adicionada ao diesel e gasolina. Nesta dissertação, monóxido de carbono, óxidos de nitrogênio, compostos orgânicos voláteis e ozônio são os principais poluentes estudados para a cidade do Rio de Janeiro. A avaliação de cenários foi realizada empregando um modelo de qualidade do ar com base no modelo de trajetórias OZIPR e no modelo químico SAPRC. Os resultados demonstram que o aumento do uso de biodiesel diminui a concentração de ozônio na atmosfera em relação ao caso base estudado, em 10,23% utilizando a mistura BE Diesel, em 5,28% utilizando a mistura B20 e apenas 0,33% utilizando a mistura B10. Isso de fato acontece pois o biodiesel possuiu mais oxigenados na sua estrutura. O estudo revelou o aumento nas emissões e concentrações de NOx na troposfera da cidade do Rio de Janeiro, para as mistura de BE- Diesel, B20 e o B10 em relação ao caso base, foram de 5,66%, 2,83% e 0% respectivamente. Em especial na cidade do Rio de janeiro existe um consumodo O3 pelo NO presente na atmosfera urbana. A redução na concentração de CO para as misturas de BE- Diesel, B20 e o B10 em relação ao caso base, foram de 13,11%, 6,56% e 4,12% respectivamente. O aumento da disponibilidade de oxigênio presente no biodiesel produz uma melhor queima do combustível, contribuindo para uma menor concentração de CO na troposfera

Uma Simulação computacional do passeio aleatório simples

Em 1828 foi observado um fenômeno no microscópio em que se visualizava minúsculos grãos de pólen mergulhados em um líquido em repouso que mexiam-se de forma aleatória, desenhando um movimento desordenado. A questão era compreender este movimento. Após cerca de 80 anos, Einstein (1905) desenvolveu uma formulação matemática para explicar este fenômeno, tratado por movimento Browniano, teoria cada vez mais desenvolvida em muitas das áreas do conhecimento, inclusive recentemente em modelagem computacional. Objetiva-se pontuar os pressupostos básicos inerentes ao passeio aleatório simples considerando experimentos com e sem problema de valor de contorno para melhor compreensão ao no uso de algoritmos aplicados a problemas computacionais. Foram explicitadas as ferramentas necessárias para aplicação de modelos de simulação do passeio aleatório simples nas três primeiras dimensões do espaço. O interesse foi direcionado tanto para o passeio aleatório simples como para possíveis aplicações para o problema da ruína do jogador e a disseminação de vírus em rede de computadores. Foram desenvolvidos algoritmos do passeio aleatório simples unidimensional sem e com o problema do valor de contorno na plataforma R. Similarmente, implementados para os espaços bidimensionais e tridimensionais,possibilitando futuras aplicações para o problema da disseminação de vírus em rede de computadores e como motivação ao estudo da Equação do Calor, embora necessita um maior embasamento em conceitos da Física e Probabilidade para dar continuidade a tal aplicação.

Perfil dos HPA prioritários na exaustão de veículo a diesel, no combustível diesel utilizado durante os ensaios de emissão veicular e no óleo lubrificante do motor

As concentrações na exaustão e os fatores de emissão dos hidrocarbonetos policíclicos aromáticos (HPA) prioritários de um veículo a diesel e as suas respectivas concentrações no diesel usado durante os ensaios de emissão veicular foram determinados com a finalidade de estimar a contribuição dos HPA provenientes do combustível nas emissões. Os produtos da combustão foram coletados diretamente nas emissões brutas do escapamento, utilizando um sistema de amostragem a volume constante sem diluição dos gases da exaustão. Os HPA associados ao MP foram amostrados de forma estratificada, utilizando um impactador em cascata MOUDI e filtros de fibra de vidro como substratos, e os HPA em fase gasosa foram amostrados usando cartuchos de amberlite XAD-2. A concentração dos HPA no óleo lubrificante do motor também foi monitorada ao longo do tempo até a sua troca após 12.000 km de uso. Após a extração e tratamento das amostras, a identificação e quantificação dos HPA foram realizadas, utilizando cromatografia de fase gasosa acoplada à espectrometria de massas (CG-EM) com injetor de grande volume de vaporização com a temperatura programável (PTV-LVI). Cinco variáveis do PTV-LVI foram otimizadas, utilizando planejamento de experimentos, o que permitiu obter limites de detecção menores do que 2,0 g L-1. Somente 7 dos 16 HPA prioritários foram identificados na exaustão: NAP, ACY, ACE, FLU, FEN, FLT e PYR. Os ensaios de emissão veicular foram realizados com o veículo em modo estacionário, sem aplicação de carga e com baixa velocidade de rotação do motor (1500 rpm), utilizando um diesel com menor teor de enxofre (10 mg kg-1) e com 5% v/v de biodiesel. Esses fatores possivelmente contribuíram para reduzir as emissões dos outros 9 HPA a valores abaixo dos limites de detecção do método desenvolvido. Aproximadamente 80% da massa dos HPA totais associados ao MP estavam presentes em partículas com tamanho entre 1,0 m e 56 nm, e aproximadamente 4,5% estavam presentes em partículas menores do que 56 nm. Partículas menores que 2,5 m são facilmente inaladas e depositadas no trato respiratório e na região alveolar, justificando a preocupação com relação às emissões de HPA associados a partículas provenientes da exaustão veicular de motores a diesel. Somente 5 dos 7 HPA identificados na exaustão foram detectados no diesel: NAP, ACY, FLU, FEN e PYR. A razão entre os fatores de emissão (g L-1diesel) dos HPA na exaustão e suas respectivas concentrações do diesel (g L-1) variaram de 0,01 0,02 a 0,05 0,029, dependendo do HPA. Esses valores indicam que pelo menos 95 a 99% dos HPA identificados no diesel foram destruídos e/ou transformados em outros compostos durante a combustão, e/ou foram retidos no reservatório do óleo lubrificante. Por outro lado, os HPA que tiveram maiores concentrações no diesel também apresentaram maiores fatores de emissão, o que sugere que os HPA provenientes do diesel possuem uma contribuição significativa para as emissões dos HPA totais. O perfil dos HPA prioritários no óleo lubrificante mostrou-se semelhante ao perfil dos HPA no diesel e nas emissões totais, onde o NAP, FEN e PYR foram os HPA majoritários

Avaliação do possível impacto das técnicas de MEOR (Microbial Enhanced Oil Recovery) no fator de recuperação das reservas de petróleo e gás do Brasil

Os métodos tradicionais de estimular a produção de petróleo, envolvendo a injeção de água, vapor, gás ou outros produtos, estabeleceram a base conceitual para novos métodos de extração de óleo, utilizando micro-organismos e processos biológicos. As tecnologias que empregam os processos de bioestimulação e bioaumentação já são amplamente utilizadas em inúmeras aplicações industriais, farmacêuticas e agroindustriais, e mais recentemente, na indústria do petróleo. Dada a enorme dimensão econômica da indústria do petróleo, qualquer tecnologia que possa aumentar a produção ou o fator de recuperação de um campo petrolífero gera a expectativa de grandes benefícios técnicos, econômicos e estratégicos. Buscando avaliar o possível impacto de MEOR (microbial enhanced oil recovery) no fator de recuperação das reservas de óleo e gás no Brasil, e quais técnicas poderiam ser mais indicadas, foi feito um amplo estudo dessas técnicas e de diversos aspectos da geologia no Brasil. Também foram realizados estudos preliminares de uma técnica de MEOR (bioacidificação) com possível aplicabilidade em reservatórios brasileiros. Os resultados demonstram que as técnicas de MEOR podem ser eficazes na produção, solubilização, emulsificação ou transformação de diversos compostos, e que podem promover outros efeitos físicos no óleo ou na matriz da rocha reservatório. Também foram identificadas bacias petrolíferas brasileiras e recursos não convencionais com maior potencial para utilização de determinadas técnicas de MEOR. Finalmente, foram identificadas algumas técnicas de MEOR que merecem maiores estudos, entre as técnicas mais consolidadas (como a produção de biossurfatantes e biopolímeros, e o controle da biocorrosão), e as que ainda não foram completamente viabilizadas (como a gaseificação de carvão, óleo e matéria orgânica; a dissociação microbiana de hidratos de gás; a bioconversão de CO2 em metano; e a bioacidificação). Apesar de seu potencial ainda não ser amplamente reconhecido, as técnicas de MEOR representam o limiar de uma nova era na estimulação da produção de recursos petrolíferos existentes, e até mesmo para os planos de desenvolvimento de novas áreas petrolíferas e recursos energéticos. Este trabalho fornece o embasamento técnico para sugerir novas iniciativas, reconhecer o potencial estratégico de MEOR, e para ajudar a realizar seu pleno potencial e seus benefícios.

Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.

Determinação experimental e modelagem termodinâmica do equilíbrio de fases em sistemas com CO2, n-hexano, n-hexadecano e tetralina

A descoberta de petróleo na camada de Pré-Sal possibilita a geração de ganhos em relação à dependência energética do país, mas também grandes desafios econômicos e tecnológicos. Os custos de extração são maiores devido a vários fatores como a exigência de equipamentos de exploração que suportem elevadas pressões, altas temperaturas e grandes concentrações de gases ácidos, tais quais, dióxido de carbono (CO2) e sulfeto de hidrogênio (H2S). Uma das principais preocupações com o CO2 é evitar liberá-lo para a atmosfera durante a produção. Com a modelagem termodinâmica de dados de equilíbrio de sistemas envolvendo CO2 supercrítico e hidrocarbonetos é possível projetar equipamentos utilizados em processos de separação. A principal motivação do trabalho é o levantamento de dados de equilíbrio de fases de sistemas compostos de CO2 e hidrocarbonetos, possibilitando assim prever o comportamento dessas misturas. Os objetivos específicos são a avaliação do procedimento experimental, a estimação e predição dos parâmetros de interação binários para assim prever o comportamento de fases dos sistemas ternários envolvendo CO2 e hidrocarbonetos. Duas metodologias foram utilizadas para obtenção dos dados de equilíbrio: método estático sintético (visual) e método dinâmico analítico (recirculação das fases). Os sistemas avaliados foram: CO2 + n-hexano, CO2 + tetralina, CO2 + n-hexadecano, CO2 + n-hexano + tetralina e CO2 + tetralina + n-hexadecano à alta pressão; tetralina + n-hexadecano à baixa pressão. Para o tratamento dos dados foi utilizada equação de estado cúbica de Peng-Robinson e a regra de mistura clássica

Avaliação técnico-econômica da implementação de um sistema de cultivo de microalgas na usina termelétrica Barbosa Lima Sobrinho com vista à biofixação de CO2

Este estudo teve como finalidade levantar dados para uma avaliação das alternativas tecnológicas (cultivos de microalgas e reflorestamento) para a biofixação de CO2 da atmosfera próxima à usina termelétrica; tendo sido utilizada como referência a Usina Barbosa Lima Sobrinho. Já existe um projeto de avaliação do efeito do reflorestamento na fixação do CO2 nesta usina e, neste trabalho, foi avaliada a alternativa do cultivo de microalgas. Uma pesquisa inicial foi feita na literatura para verificar qual a espécie de microalga seria a mais adequada para ser utilizada no estudo, tendo sido a espécie Chlorella sp. a selecionada. Posteriormente os sistemas de cultivo de microalgas mais comumente empregados no mercado foram levantados e foi selecionado o cultivo em tanques abertos como referência para a modelagem do processo. Utilizando os dados da termelétrica e da literatura foi possível estimar a quantidade de CO2 que será retirada da atmosfera caso um sistema de cultivo seja efetivamente instalado na usina termelétrica. Uma análise econômica foi realizada para determinar a viabilidade do projeto. Os resultados indicam que a utilização deste tipo de tecnologia é promissora

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.

Simulação numérica do processo de alteamento de áreas de deposição de resíduos pelo método a montante.

A previsão do comportamento de resíduos constituiu-se em um desafio geotécnico, uma vez que estes materiais apresentam uma resposta distinta dos materiais usualmente encontrados em depósitos naturais. A análise dos recalques da fundação, decorrentes da sobrecarga imposta pelo alteamento, é complexa, tendo em vista que o adensamento de resíduos pressupõe grandes deformações, invalidando o uso de teorias clássicas de adensamento. Atualmente, no Brasil, a técnica de disposição de resíduos de bauxita prevê uma operação inicial de lançamento no interior de lagos artificiais, em forma de polpa. Após o esgotamento do lago e ressecamento do resíduo, inicia-se o lançamento pelo método a montante. Neste método, a polpa é lançada sobre o resíduo pré-existente, que se encontra em processo de adensamento. O presente trabalho tem como objetivo reproduzir numericamente o comportamento de áreas de resíduos durante a etapa de alteamento a montante. A pesquisa tem como enfoque 2 áreas de resíduos de bauxita. Uma delas encontra-se em fase de reabilitação e dispõe de instrumentação de campo (recalques e deslocamentos horizontais). A outra se encontra em fase de operação do alteamento e dispõe de dados experimentais. Desta forma, a metodologia consistiu na reprodução numérica do processo de alteamento da área instrumentada e comparação dos resultados com a instrumentação de campo, com objetivo de avaliar o modelo numérico e os parâmetros do resíduo. Posteriormente, realizou-se a previsão do comportamento do resíduo de fundação da área em fase de alteamento. Os parâmetros geotécnicos foram definidos a partir de um extenso programa de ensaios de campo e laboratório, executado no local em estudo, fazendo-se uso de um tratamento estatístico dos dados experimentais. Os resultados numéricos mostraram a potencialidade do programa na previsão do comportamento de áreas de resíduos durante o alteamento a montante, com previsões de recalques e deslocamentos horizontais coerentes com a instrumentação de campo.

Emissão de CO2 do solo e sua correlação com a rizosfera de diferentes paisagens de áreas mineradas do município de Santo Antônio de Pádua-RJ

Objetivando avaliar o comportamento das emissões de CO2 do solo em áreas mineradas do município de Santo Antônio de Pádua-RJ e sua correlação com a rizosfera, este presente estudo utilizou uma câmara fechada com sensor de infravermelho em três diferentes paisagens, a saber: A-1 (área reflorestada há 10 anos), A-2 (área desmatada) e A-3 (área em processo de recuperação). Em cada área foi instalada três cilindros para efeitos de repetição. O monitoramento foi realizado durante os meses representativos de cada estação do ano de 2013, sendo a análise realizada durante dois dias consecutivos. Concomitantemente as coletas de CO2 foram realizadas coletas de dados de temperatura e umidade do solo, sendo também avaliadas informações pedológicas através das análises de granulometria, porosidade, pH, carbono orgânico e matéria orgânica. Informações meteorológicas e microclimatológicas também foram extraídas através de uma estação meteorológica automática e através de sensores portáteis. Os resultados permitem concluir que existe uma variação sazonal dos fluxos de CO2, havendo uma tendência de máximos de emissão durante o verão e de mínimas durante o inverno, sendo o outono e a primavera marcados por valores medianos. A correlação das áreas entre os dois dias monitoramento indicam que as emissões foram semelhantes ao da análise em dia anterior, apresentando uma correlação significativa a 5% para A-1 e A-2 e de 1% para A-3. A-1 e A-2 apresentaram emissões de CO2 mais homogêneas que A-3, havendo, entretanto, um maior fluxo de CO2 durante o verão para todas as áreas. Os dados de MOS, COS e pH demonstraram não haver uma correlação direta com as emissões de CO2. Os dados de porosidade e densidade, porém, apontam para uma possível correlação com as menores emissões de CO2 em A-3 devido a menor porosidade e maior densidade de seus solos. A temperatura do solo foi a variável que mais se correlacionou com as emissões de CO2, havendo um índice igual a r =0,68 para A-1 e de r =0,74 para A-2, sendo que em A-3 esta correlação não foi significativa. A temperatura do ar demonstrou uma correlação somente na área descampada de A-2. No que se refere à correlação da umidade do solo não houve correlações diretas significativas, sendo que somente houve uma correlação negativa (r=-0,50) significativa a 5% em A-3 com a umidade do ar. O diagnóstico ambiental das áreas de monitoramento revela que estas possuem baixos indicadores de qualidade, sendo afetados também pela escassez hídrica da região durante oito meses do ano. A-1 apresenta os melhores indicadores biológicos, químicos e microclimáticos, seguidos por A-2 e A-3 que apresentam diversas deficiências e problemas em termos de estrutura e atividade biológica dos solos. Neste estudo permite-se concluir que áreas próximas, porém com características distintas podem produzir diferentes padrões de emissão de CO2, dificultando, portanto, estimativas globais de emissão de CO2. Os elementos mais associados às emissões de CO2 parecem estar relacionados à temperatura do solo e do ar, umidade do ar e estrutura do solo, havendo, entretanto, outros fatores que podem estar indiretamente relacionados e que exercem diferentes influências de acordo com o ambiente analisado.

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.