Estudo de fluxos numéricos para modelos magneto-hidrodinâmicos utilizando o código Flash

A teoria magneto-hidrodinâmicos permite a estruturação de modelos computacionais, designados modelos MHDs, que são uma extensão da dinâmica dos fluidos para lidar com fluidos eletricamente carregados, tais como os plasmas, em que se precisa considerar os efeitos de forças eletromagnéticas. Tais modelos são especialmente úteis quando o movimento exato de uma partícula não é de interesse, sendo que as equações descrevem as evoluções de quantidades macroscópicas. Várias formas de modelos MHD têm sido amplamente utilizadas na Física Espacial para descrever muitos tipos diferentes de fenômenos de plasma, tais como reconexão magnética e interações de ventos estelares com diferentes objetos celestiais. Neste trabalho, o objetivo é analisar o comportamento de diversos fluxos numéricos em uma discretização de volumes finitos de um modelo numérico de MHD usando um esquema de malha entrelaçada sem separação direcional considerando alguns casos testes. Para as simulações, utiliza-se o código Flash, desenvolvido pela Universidade de Chicago, por ser um código de amplo interesse nas simulações astrofísicas e de fenômenos no espaço próximo à Terra. A metodologia consiste na inclusão de um fluxo numérico, permitindo melhoria com respeito ao esquema HLL.

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.

Biomassa arbórea em fragmentos da Mata Atlântica Fluminense.

2008

Decomposição de resíduos de milho e soja em função do tempo e do manejo do solo.

A evolução da decomposição dos resíduos culturais é um aspecto muito complexo, ainda pouco conhecido e estudado, porém de grande importância na compreensão de processos como a formação da matéria orgânica e a manutenção da fertilidade e sustentabilidade dos solos. A sua velocidade de decomposição determina o tempo de permanência da cobertura morta na superfície do solo, sendo variável conforme a constituição dos diferentes tecidos vegetais. O objetivo deste trabalho foi a obtenção de equações matemáticas, que expliquem a evolução da dinâmica de decomposição de resíduos de milho e soja, em semeadura direta e convencional, utilizando- se dados de um ensaio de longa duração (12 anos) conduzido no campo experimental da Embrapa Soja, em Londrina, PR. Foram obtidas equações exponenciais, que mostraram semelhanças na velocidade de decomposição dos resíduos de milho e soja e níveis de decomposição mais elevados em semeadura convencional,tanto para o milho quanto para a soja.

Decomposição de resíduos de aveia e trigo em função do tempo e do manejo do solo.

O uso excessivo de máquinas e implementos agrícolas causou sérios problemas de compactação, erosão, redução dos níveis de matéria orgânica e fertilidade dos solos. Neste cenário, estudos sobre a evolução da decomposição de resíduos de culturas são importantes para a compreensão de processos como o da formação da matéria orgânica, crucial para a manutenção da fertilidade e sustentabilidade dos solos. Foram conduzidos estudos da dinâmica da decomposição de resíduos de aveia e trigo, em semeadura direta e convencional, utilizando-se dados de 13 anos de um ensaio conduzido no campo experimental da Embrapa Soja, em Londrina, PR. Foram obtidas equações exponenciais, que mostraram a melhor representação do processo de decomposição ao longo do tempo. A velocidade de decomposição foi elevada nas primeiras semanas, para ambas as culturas, sendo observados níveis de decomposição mais elevados em semeadura convencional quando em comparação com o plantio direto, tanto para a aveia quanto para o trigo.

Um modelo determinístico para avaliação do controle biológico de praga de citros.

O objetivo desse trabalho é apresentar um modelo matemático compartimental e determinístico para analisar a dinâmica populacional do inseto-praga Phyllocnistis citrella em interação com seus inimigos naturais Galeopsomyia fausta (parasito nativo) e Ageniaspis citricola (parasito exótico). Com o estudo analítico do sistema de equações diferenciais que descreve o modelo, avalia-se a eficiência do controle biológico da P.citrella através da possibilidade de coexistência entre as espécies e da competição entre os dois inimigos naturais.

Avaliação da composição estrutural e biomassa viva acima do solo, em florestas sob efeito da fragmentação na Amazônia Ocidental.

A finalidade deste estudo foi avaliar o impactos do processo de fragmentação florestal sobre a estrutura horizontal do fragmento e biomassa viva acima do solo, nos ambientes de borda e interior. Para esta avaliação foram considerados os seguintes parâmetros: uso atual da terra; composição florestal; mensuração da biomassa vegetal viva acima do solo, por meio de equações alométricas ajustadas para a região e estimativa do estoque de carbono.

Cálculo diferencial e integral e tecnologias digitais: perspectivas de exploração no sofware GeoGebra

A ênfase algébrica dada ao longo do tempo nos cursos de Cálculo Diferencial e Integral não oportunizou que tratamentos gráficos e numéricos fossem privilegiados, visto a ausência de softwares que possibilitassem uma abordagem diferenciada aos conceitos inerentes a esta disciplina (Richit, 2010, Guimarães, 2001). Contudo, iniciativas no mundo inteiro têm dedicado esforços e desenvolvido softwares que possibilitam explorações qualitativamente diferentes para conceitos de Cálculo a partir de representações gráficas, numéricas ou algébricas envolvendo visualização, a simulação, o aprofundamento do pensamento matemático, conjecturas e validações, etc. Deste modo, a incorporação das tecnologias digitais na aula de Cálculo remove um pouco o fardo da manipulação algébrica, possibilitando a transição entre a ação física (interação do estudante com a tecnologia) e a representação matemática de um conceito. Assim, a proposta de oficina aqui apresentada objetiva explorar conceitos de Cálculo (Funções, Limites, Derivadas e Integrais) em uma perspectiva de investigação com o software GeoGebra.

Uma abordagem da lógica fuzzy no ensino médio

O uso cada vez mais intenso de conceitos subjetivos - provenientes da lógica fuzzy - aplicados a problemas reais nos motivou a desenvolver procedimentos para a introdução destas idéias no âmbito do ensino médio. O projeto propõe inicialmente o estudo de conjuntos fuzzy que podem ser entendidos com exemplos - de variação populacional, de controle de pragas e de epidemias. Posteriormente, usar as “operações fuzzy” Sup e Inf em produtos de matrizes para realizar diagnósticos e avaliações subjetivas. As situações abordadas já estão na literatura (Barros e Bassanezi, 2006), entretanto não como fonte para o Ensino Médio. Um dos objetivos principais deste trabalho é contrapor a crença de exatidão da matemática clássica com os resultados provenientes de lógica subjetiva, utilizando conceitos apropriados para os estudantes destas séries: teoria dos conjuntos, relações e funções, matrizes, equações de diferenças e outros.

A influência do terrorismo nos comportamentos de viagem internacionais

A partir de uma amostra de 600 turistas internacionais que circulam em Portugal, Espanha e Itália, este estudo identifica as principais os conceitos chave relacionados com o terrorismo, a percepção de risco, envolvimento e motivação para a segurança dos turistas internacionais. Diferentes níveis de preocupação relativamente à segurança pode influenciar as decisões dos turistas. No seu processo de decisão, os turistas avaliam vários factores, nomeadamente, o nível de risco ou de segurança que consideram nos destinos (Sonmez, 1998). Os turistas adoptam uma atitude protectora alterando os seus comportamentos durante os processos de decisão, substituindo os destinos que consideram inseguros por outros associados a uma maior segurança (Gu & Martin, 1992; Mansfeld, 1996). O terrorismo exacerbado pelos media tem efeitos graves nas receitas dos destinos turísticos (Taylor, 2006). Através da publicidade negativa, um destino turístico que experiencia um incidente terrorista pode ver a sua reputação danificada e a actividade turística severamente comprometida (Sonmez, 1998). Inclusivamente, a imageme negativa de um destino pode ser generalizada e pode também afectar outros países ou regiões por períodos de tempo indeterminados (Taylor, 2006). Um modelo de equações estruturais revela que os turistas são motivados para adquirir informação sobre o terrorismo nos media, nomeadamente mostram atenção e interesse sobre essas notícias e esse facto influencia directamente o seu risco percebido. A percepção de risco influencia directamente o envolvimento dos turistas no planeamento da viagem, especificamente a procura de informação antes da viagem e quando estão no destino. A percepção de risco e o envolvimento dos turistas influencia a percepção da importância da segurança.A discussão foca as implicações deste modelo para a teoria e para as instituições e organizações turísticas. São igualmente apresentadas recomendações para os gestores e promotores dos destinos e para os gestores das organizações turísticas. Direcções futuras de investigação são igualmente apresentadas.

Matérias-primas da plataforma do Mondego para cerâmica

O presente trabalho divulga os resultados dos estudos levados a efeito nas matérias-primas argilosas que se integram tipologicamente na argila comum, ocorrentes na designada Plataforma do Mondego, Centro de Portugal, na região entre Miranda do Corvo, a SW, e Tábua, a NE. Para tal realizou-se a cartografia superficial dos sedimentos continentais do Cretácico e do Terciário, aí preservados, estabeleceram-se as colunas sedimentares, a nível regional e a nível local, nas áreas de maior potencial reconhecido e, caracterizaram-se as matérias-primas argilosas amostradas, em termos de composição, textura e aptidão cerâmica, objectivando a definição das zonas das jazidas com maior interesse industrial. A informação obtida e compilada, relativa ao caulino e ao feldspato é também apresentada, com vista a uma percepção da potencialidade destas matériasprimas cerâmicas ocorrentes na área. O estudo de caracterização dos materiais argilosos investigados teve como base 53 amostras obtidas em seis regiões da Plataforma do Mondego aqui definidas por motivo de organização do trabalho, segundo os critérios geográfico e geotectónico. As formações de Côja e de Campelo, do Terciário, constituem as unidades onde ocorrem litótipos produtivos, em termos de matéria-prima para a Cerâmica de Barro Vermelho ou Cerâmica de Construção. Relativamente às características texturais e composicionais, em síntese, apresentam-se os factos relevantes seguintes: A matéria-prima argilosa existente na área estudada materializa, granulometricamente, na maioria das amostras, silte, caindo os níveis amostrados com maior percentagem de argila, no domínio do silte argiloso. O défice em fracção argila implica limitações quanto à possibilidade de diversificação de produtos cerâmicos fabricados com estas matérias-primas tal-qual. Os diferentes métodos analíticos utilizados na caracterização mineralógica dos materiais amostrados confirmaram uma composição em termos dos minerais argilosos, consistindo de ilite/mica (em geral, o mais abundante), caulinite e esmectite, interestratificados e clorite. Os minerais não argilosos são quartzo (predominante), feldspato (sobretudo potássico) e hematite, com uma representação baixa. A mineralogia da fracção inferior a 2μm das amostras, não difere das amostras totais, salvo no teor mais elevado dos minerais argilosos e acentuada redução dos minerais não argilosos. Os resultados da análise química por fluorescência de raios X das amostras integrais correlacionam-se com as características mineralógicas observadas através das técnicas analíticas utilizadas. No respeitante às propriedades e comportamento cerâmico verifica-se: Os parâmetros relacionados com a plasticidade indicam que parte das pastas elaboradas com estes materiais argilosos têm uma trabalhabilidade aceitável, mas existem problemas de conformação e acentuada retracção num número significativo de amostras, devido à elevada plasticidade da maioria das amostras. A extrusão é satisfatória a óptima. A RMF e a retracção em seco assumem valores, respectivamente, moderados a baixos e moderados, embora seja necessário ter em conta o procedimento de extrusão dos provetes, sem dispositivo de vácuo. Todas as amostras foram sujeitas a cozedura a 900ºC, e um conjunto seleccionado foi cozido a 1000ºC e a 1100ºC. As fases mineralógicas ocorrentes após cozeduras a 900ºC e 1100ºC foram identificadas num conjunto de amostras, tendo-se evidenciado a coerência dessas fases, com a mineralogia das amostras em seco. Após cozedura a 900ºC, os valores de RMF das amostras satisfazem geralmente os valores mínimos, exigidos para o fabrico de tijolo, abobadilha e, com alguma frequência, de telha, como já se verificava com os valores daquela propriedade em seco. Os valores de retracção seco-cozido são em geral, modestos. A capacidade de absorção de água é maioritariamente elevada. A formação de vidro, sobretudo, condiciona o comportamento destas propriedades por cozedura dos provetes a 1000ºC e a 1100ºC. A coloração predominante em cru das matérias-primas argilosas amostradas é amarelo acastanhado a castanho avermelhado. Após cozedura a 900ºC, há um acentuado escurecimento e incremento no grau de vermelho. As cozeduras a 1000ºC e 1100ºC promovem escurecimento gradual, com ligeira influência na cor. O comportamento dos provetes após as cozeduras cerâmicas revelou-se homogéneo a cada uma das respectivas temperaturas, não se registando também defeitos significativos, nem eflorescências. A análise das amostras em termos composicionais e tecnológicos permitiu destacar as principais características e aspectos distintivos das matériasprimas, nas diferentes regiões definidas objectivando as diferentes potencialidades cerâmicas. Nesta abordagem comparativa foram consideradas só as amostras dos campos silte e silte arenoso, por serem aquelas com maior interesse para a Cerâmica de Construção. Nas colunas sedimentares das regiões de Tábua e Santa Quitéria constata-se a ocorrência de dois ritmos de sedimentação, que embora assumam características específicas em cada região, têm aspectos composicionais e tecnológicos em comum, traduzindo melhor aptidão cerâmica os ritmos inferiores. A conjugação da cartografia realizada com os estudos laboratoriais permite concluir que as regiões de Tábua e de Santa Quitéria serão as que têm maior potencial por explorar, em matéria-prima para Cerâmica de Construção, apesar da primeira já ser intensamente explorada. Na região de Tábua, as amostras têm como fases mineralógicas principais ilite e quartzo na mesma proporção média (35%) e caulinite (média=19%) que regista enriquecimento significativo na fracção argila (média=38%). As argilas desta região registam a cor em cru mais vermelha e pH mais ácido observados. O ritmo de sedimentação inferior, com esmectite e interestratificados e ligeiramente menos quartzoso, apresenta melhores propriedades cerâmicas. A composição mineralógica média das amostras da região de Santa Quitéria é próxima daquela da região de Tábua, mas menos caulinítica, em especial na sequência inferior, na qual o teor médio de caulinite na fracção argila (7%) é o mais baixo observado. Na região de Côja – Arganil, a actividade extractiva é significativa na bacia de Côja. Aqui, a exploração de novas áreas potenciais é condicionada pela cobertura conglomerática e por estruturas tectónicas e não tectónicas relacionadas com comportamento plástico. Na restante área desta região, os recursos argilosos são penalizados por material areno-conglomerático. As amostras desta região distinguem-se das restantes a nível textural pela maior fracção areia e mineralogicamente pela presença de clorite, teor reduzido de caulinite e elevado de feldspato. As potencialidades em barro vermelho na região de Sanguinheda não serão significativas, pois a Formação de Côja é predominantemente arcósica e a Formação de Campelo pouco espessa e conglomerática. As argilas com melhor aptidão cerâmica foram amostradas na região de Miranda do Corvo – Lousã. Na composição, estas amostras são as que contêm maior fracção de argila, maior teor de ilite e caulinite e menor de argilas expansivas. Registam os melhores valores nas características tecnológicas, nomeadamente a RMF e absorção de água. A consistência dos grupos amostrais e ritmos definidos em termos composicionais e tecnológicos é corroborada pelas técnicas de análise estatística multivariada aplicadas, que os identificam. As condições de amostragem na região de Tábua, permitiram a elaboração de uma coluna tipológica, na qual, a partir de uma caracterização expedita de amostras é possível identificar a sua afinidade com os ritmos argilosos definidos e, consequentemente, a sua situação na coluna sedimentar regional e aptidão cerâmica. A cartografia dos recursos argilosos potenciais, elaborada à escala 1:25.000, constitui um dos objectivos principais deste trabalho e nela constam unidades litológicas, nas quais são diferenciadas unidades argilosas tendo também em consideração a tipologia e guias mineralógicos resultantes da caracterização das amostras. As características composicionais destas unidades denotam diferente aptidão cerâmica e, portanto, permitem salientar as zonas com maior interesse económico.

Exact and numerical solutions for diffusio-reaction equations with convection term

Esta dissertação estuda essencialmente dois problemas: (A) uma classe de equações unidimensionais de reacção-difusão-convecção em meios não uniformes (dependentes do espaço), e (B) um problema elíptico não-linear e paramétrico ligado a fenómenos de capilaridade. A Análise de Perturbação Singular e a dinâmica de Hamilton-Jacobi são utilizadas na obtenção de expressões assimptóticas para a solução (com comportamento de frente) e para a sua velocidade de propagação. Os seguintes três métodos de decomposição, Adomian Decomposition Method (ADM), Decomposition Method based on Infinite Products (DIP), e New Iterative Method (NIM), são apresentados e brevemente comparados. Adicionalmente, condições suficientes para a convergência da solução em série, obtida pelo ADM, e uma aplicação a um problema da Telecomunicações por Fibras Ópticas, envolvendo EDOs não-lineares designadas equações de Raman, são discutidas. Um ponto de vista mais abrangente que unifica os métodos de decomposição referidos é também apresentado. Para subclasses desta EDP são obtidas soluções numa forma explícita, para diferentes tipos de dados e usando uma variante do método de simetrias de Bluman-Cole. Usando Teoria de Pontos Críticos (o teorema usualmente designado mountain pass) e técnicas de truncatura, prova-se a existência de duas soluções não triviais (uma positiva e uma negativa) para o problema elíptico não-linear e paramétrico (B). A existência de uma terceira solução não trivial é demonstrada usando Grupos Críticos e Teoria de Morse.

Caracterização de novas membranas de titanossilicatos microporosos por ensaios de permeabilidade

O interesse crescente das membranas inorgânicas deve-se à potencial aplicação em novas áreas de investigação e da indústria, e em alternativa a operações mais convencionais. Em particular, as membranas de titanossilicatos oferecem vantagens importantes sobre as de zeólitos, pois podem ser sintetizadas sem agentes estruturantes orgânicos, para evitar a calcinação subsequente usualmente responsável por defeitos irreversíveis, exibem novas possibilidades de substituição isomórfica da matriz, permitindo um ajuste mais fino das propriedades catalíticas e de adsorção, e são capazes de separar misturas com base em diferenças de afinidade e tamanho molecular (efeito de peneiro). Os objectivos principais deste trabalho foram: i) a caracterização dinâmica de membranas do tipo zeolítico sintetizadas no Laboratório Associado CICECO, realizando-se experiências de permeação com gases puros e misturas; ii) o desenvolvimento e validação de novos modelos para a transferência de massa multicomponente através de membranas porosas pela abordagem de Maxwell-Stefan, tendo em conta os mecanismos específicos encontrados, particularmente a contribuição por difusão superficial; e iii) a modelação dos pontos experimentais medidos, bem como dados compilados da literatura. De forma a realizar os ensaios de permeação, desenhou-se, montou-se e testou-se uma instalação experimental. Para gases puros, os objectivos principais foram a medição de permeâncias a temperatura constante, por variação da pressão transmembranar r ( ΔP ), e de permeâncias a temperatura programada, conduzidas a ΔP constante. Seguidamente, calcularam-se as selectividades ideais. Em relação a misturas, a determinação de selectividades reais requer as fracções molares no permeado e no retido. Na globalidade, estudaram-se três suportes diferentes (aço inoxidável e α − alumina) e dezanove membranas de AM-3, ETS-10, ZSM-5 e zeólito 4A, utilizando-se H2, He, N2, CO2, e O2. A primeira avaliação exploratória da qualidade das membranas foi feita permeando azoto à temperatura ambiente. Assim, permeâncias superiores a 10−6 mol/m2s.Pa evidenciavam defeitos grosseiros, levando-nos a efectuar cristalizações adicionais sobre as primeiras camadas. Este procedimento foi implementado com oito membranas. Um trabalho experimental mais detalhado foi conduzido com cinco membranas. Membranas com curvas permeância-temperatura ( Π −T ) decrescentes indicam tipicamente transporte viscoso e de Knudsen, i.e. meso e macrodefeitos. Por exemplo, a membrana nº 3 de AM-3 exibiu este comportamento com H2, He, N2 e CO2 puros. A contribuição de Knudsen foi confirmada pela relação linear encontrada entre as permeâncias e o inverso da raiz quadrada da massa molar. O mecanismo viscoso foi também identificado, pois as permeâncias eram inversamente proporcionais à viscosidade do gás ou, atendendo a equações do tipo de Chapman-Enskog, directamente proporcionais a 2 0.5 k d M (onde k d é o diâmetro cinético e M a massa molar). Um comportamento de permeação distinto observou-se com a membrana nº 5 de AM-3. As permeâncias registadas a temperatura programada eram aproximadamente constantes para o N2, CO2 e O2, enquanto com o H2 cresciam significativamente. Conjuntamente elas evidenciam a ocorrência de macro, meso e microdefeitos intercristalinos. O transporte gasoso activado através dos microporos compensa o impacto diminuidor dos meso e macroporos. Ao contrário do N2, CO2 e O2, o pequeno diâmetro do hidrogénio torna-lhe possível permear através dos microporos intracristalinos, o que lhe adiciona um mecanismo de transferência responsável por esse crescimento. No que respeita à difusão superficial, o sistema CO2/ZSM-5 pode ser tomado como um exemplo paradigmático. Uma vez que este zeólito adsorve o CO2, as permeâncias diminuem com o crescimento de ΔP , em virtude de as concentrações no sólido aumentarem de forma não linear e tenderem para a saturação. Os resultados contrastantes obtidos com azoto realçam ainda mais o mecanismo superficial, pois o N2 não é adsorvido e as permeâncias medidas são constantes. Globalmente, as selectividades ideais calculadas ( α* ) variam de cerca de 1 a 4.2. Este parâmetro foi também utilizado para discriminar as melhores membranas, uma vez que baixos valores de α* denotam o escoamento viscoso não-selectivo típico de macrodefeitos. Por exemplo, o H2/CO2 na membrana nº 3 de AM-3 apresentou α* = 3.6 − 4.2 para 40–120ºC, enquanto que na membrana nº 5 de AM-3 originou α* = 2.6 − 3.1. Estes resultados corroboraram as observações anteriores, segundo as quais a membrana nº 5 era melhor do que a nº 3. Alguns ensaios foram realizados com membranas saturadas com água para aumentar a selectividade: as medições mostraram claramente uma melhoria inicial seguida de uma redução consistente de α* com o aumento da temperatura, devido à remoção das moléculas de água responsáveis pela obstrução de alguns poros. Em relação às selectividades reais de misturas contendo hidrogénio, devem ser realizadas mais experiências e a quantificação do hidrogénio deve ser melhorada. No que concerne à modelação, novos factores termodinâmicos de Maxwell- Stefan foram derivados para as isotérmicas mono e multicomponente de Nitta, Langmuir-Freundlich e Toth, tendo sido testadas com dados de equilíbrio e de permeação da literatura. (É importante realçar que só estão publicadas equações para Langmuir e Dual-Site Langmuir de componentes puros e misturas). O procedimento de validação adoptado foi exigente: i) as isotérmicas multicomponente foram previstas a partir das de gás puro; ii) os parâmetros de difusão dos componentes puros foram ajustados a dados de permeação de cada gás; iii) depois, as difusividades cruzadas de Maxwell- Stefan foram estimadas pela relação de Vignes; finalmente, v) as novas equações foram testadas usando-se estes parâmetros, tendo sido capazes de estimar com sucesso fluxos binários. Paralelamente ao enfoque principal do trabalho, derivou-se um novo modelo para permuta iónica em materiais microporosos baseado nas equações de Maxwell-Stefan. Este foi validado com dados experimentais de remoção de Hg2+ e Cd2+ de soluções aquosas usando ETS-4. A sua capacidade preditiva foi também avaliada, sendo possível concluir que se comporta muito bem. Com efeito, conseguiram-se boas previsões com parâmetros optimizados a partir de conjuntos de dados independentes. Este comportamento pode ser atribuído aos princípios físicos sólidos da teoria de Maxwell-Stefan.

Phase equilibria modeling for biofuels production

Em consequência de uma série de problemas ambientais, económicos e políticos relacionados com o uso de combustíveis convencionais, vários países estão agora a focar as suas atenções em combustíveis alternativos. O biodiesel está na linha da frente das alternativas ao petróleo no sector dos transportes, sendo considerado uma opção a curto prazo visto que o seu preço é competitivo e não são necessárias mudanças nos motores para implementar o seu uso. De entres os possíveis processos de produzir biodiesel, a reacção de transesterificação com catálise básica é o método preferido. Depois da reacção são sempre necessários processos de purificação de modo ao biodiesel produzido cumprir os standards definidos para os combustíveis alternativos, reduzindo problemas de motor e consequentemente aumentando a sua aceitação por parte dos consumidores. De entre as especificações encontram-se o conteúdo em água, em álcool e em glicerol. Ser-se capaz de descrever correctamente o equilíbrio de fases de sistemas que são de interesse para os processos de purificação de biodieseis numa gama alargada de condições termodinâmicas é uma condição necessária para uma correcta simulação do processo industrial, de modo a se atingir uma elevada produtividade a baixos custos de operação. O uso de moléculas oxigenadas como combustíveis representa uma alteração significativa em termos da termodinâmica de soluções. Para combustíveis baseados em petróleo as equações de estado cúbicas e os modelos clássicos de coeficientes de actividade mostraram ser apropriados, no entanto para combustíveis novos como o biodiesel, sendo mais complexos do ponto de vista das interacções intermoleculares com formação de dipolos e pontes de hidrogénio, são necessários modelos termodinâmicos mais complexos para descrever essas interacções. Neste trabalho a CPA EoS (Cubic-Plus-Association Equation of State) será desenvolvida de modo a permitir uma descrição adequada dos equilíbrios líquido-vapor e líquido-líquido para uma serie de sistemas binários e multicomponentes contendo água, ácidos gordos, ésteres de ácidos gordos, glicerol e álcoois.

Hydraulics of stepped chutes: experimental-numerical-theoretical study

O estudo de escoamentos turbulentos em descarregadores em degraus tem sido um desafio para os investigadores. A macro-rugosidade do leito, a ondulação da superfície livre, a intermitência da localização da secção inicial de entrada de ar e o escoamento bi-fásico a jusante da secção inicial de entrada de ar fazem com que a caracterização do escoamento deslizante sobre turbilhões em descarregadores em degraus não seja simples. Actualmente, é possível combinar técnicas de medição fiáveis com simulações numéricas e análise teórica. Nesta dissertação, o estudo experimental baseia-se em resultados experimentais obtidos em duas instalações experimentais: a instalação A, do Laboratório Nacional de Engenharia Civil (LNEC), dotada de um descarregador em degraus com declive de 1V:0.75H e a instalação B, do Instituto Superior Técnico (IST), dotada de um descarregador em degraus com declive de 1V:2H. Uma sonda de ar, um tubo de Pitot modificado e vários hidrómetros permitiram o estudo do escoamento deslizante sobre turbilhões mono- e bifásico ao longo do descarregador em degraus da instalação A. Na bacia de dissipação de energia, a observação visual do escoamento foi auxiliada pelas leituras efectuadas em tomadas de pressão localizadas na soleira da bacia e pelas leituras da altura do escoamento efectuadas em réguas graduadas localizadas nas paredes da bacia. Na instalação B foram utilizados um tubo de Pitot e vários hidrómetros para estudar a região não arejada do escoamento deslizante sobre turbilhões no descarregador em degraus. O estudo numérico é baseado em simulações numéricas da região não arejada do escoamento deslizante sobre turbilhões sobre descarregadores em degraus com o código comercial de CFD FLOW-3D® de modo a reproduzir as condições ensaiadas experimentalmente. As simulações beneficiaram da técnica de blocos múltiplos (multi-block) num sistema de coordenadas cartesianas, da determinação da superfície livre pelo método TruVOF e da utilização de dois modelos de turbulência: os modelos k- e RNG k-. Por último, o estudo teórico consistiu em desenvolver um modelo simplificado 1D para determinar as características hidráulicas principais do trecho não arejado do escoamento deslizante sobre turbilhões em descarregadores em degraus. O modelo foi desenvolvido a partir das equações de Navier-Stokes, conjuntamente com resultados experimentais e numéricos. Os resultados apresentados nesta dissertação contribuem para o conhecimento do escoamento deslizante sobre turbilhões em descarregadores em degraus, nomeadamente na região não arejada, na secção inicial de entrada de ar e na região arejada. A hidráulica dos dissipadores de energia a jusante de descarregadores em degraus, em particular das bacias tipo III do USBR, é também objecto de estudo. Em relação à região não arejada do escoamento em descarregadores com declive acentuado, e com base em resultados experimentais e numéricos, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, da espessura da camada limite, da concentração média de ar, do coeficiente de energia cinética, da dissipação de energia, do factor de resistência e do coeficiente n da fórmula de Manning. São ainda propostas expressões adimensionais para a energia cinética turbulenta e sua dissipação. Para declives moderados, são propostas expressões para estimar o desenvolvimento da altura equivalente de água, do coeficiente de energia cinética e da energia específica residual. São ainda apresentados valores do expoente 1/N da expressão adimensional da distribuição de velocidades, quer para descarregadores com declive acentuado quer com declive moderado. Em conformidade com outros estudos centrados em escoamentos de parede e com derivações teóricas, para a região não arejada do escoamento em descarregadores em degraus de acentuado declive, observa-se que o factor de resistência depende da macro-rugosidade criada pelos degraus e da geometria da secção transversal e que o coeficiente n da fórmula de Manning aumenta com a rugosidade. A descrição estatística da turbulência do escoamento é igualmente explorada, contribuindo para o conhecimento da estrutura do escoamento. Observou-se que para números de Reynolds rugoso não superiores a 6.8x104 a energia cinética turbulência e a sua dissipação cumprem leis de semelhança. Estas expressões adimensionais estão de acordo com os resultados obtidos por outros autores para escoamentos completamente desenvolvidos em canais abertos e no escoamento em rios com leito de gravilha. Em acréscimo, a taxa de dissipação de energia, quer para descarregadores de declive acentuado quer de moderado declive, é baixa. Por último, observa-se que os valores da média temporal da concentração de ar entre 0 e 1 medidos na região não arejada do escoamento dizem respeito não só ao ar capturado entre ondas de água, na zona de ondulação da superfície livre, mas também ao ar emulsionado no escoamento, i.e., sob a forma de bolhas de ar, quando perto da secção média inicial de entrada de ar, devido à diferença entre localizações instantânea e média temporal. Foram revistas metodologias e fórmulas para estimar a localização da secção inicial de entrada de ar e apresentadas expressões para estimar a concentração média de ar e a altura equivalente de água nessa secção. Relativamente à região de escoamento arejado em descarregadores em degraus com declive acentuado, os resultados experimentais apresentados nesta dissertação permitiram estimar a influência da definição da superfície livre nos parâmetros hidráulicos da região do escoamento arejado e estimar a máxima elevação do escoamento nesta região do escoamento. Com base nos resultados experimentais obtidos na bacia de dissipação de energia do tipo III do USBR localizada a jusante do descarregador em degraus da instalação A, observou-se que os perfis da altura piezométrica e da altura do escoamento tendem a seguir o perfil recomendado pelo USBR para bacias tipo III. A excepção ocorre à entrada da bacia, onde as alturas piezométricas apresentadas nesta dissertação excedem largamente as apresentadas pelo USBR. É ainda observado que, tal como entre as bacias tipo I e tipo III do USBR, o ressalto hidráulico estabiliza muito mais rapidamente numa bacia tipo III a jusante de um descarregador em degraus do que uma bacia tipo I a jusante do mesmo descarregador em degraus. Finalmente, observa-se que os blocos de amortecimento a colocar no descarregador não têm influência visível nos resultados da altura piezométrica nem da altura do escoamento ao longo da bacia. Relativamente às simulações numéricas do escoamento não arejado, a proximidade entre resultados experimentais e numéricos permite validar o modelo teórico e a integração numérica usados no FLOW-3D®. As simulações desenvolvidas também mostraram que o modelo de turbulência k- permite representar as características do escoamento não arejado em descarregadores em degraus, uma vez que não foram observadas diferenças significativas entre as simulações com este modelo e com o modelo RNG k-. Finalmente, observou-se que o modelo de entrada de ar usado no FLOW-3D® é válido para estimar a localização da secção inicial de entrada de ar. Por último, a proximidade entre os resultados obtidos da aplicação do modelo teórico desenvolvido no âmbito desta dissertação e os resultados experimentais indica que as hipóteses e simplificações consideradas no desenvolvimento do modelo são adequadas.