Análise da influência de variáveis meteorológicas e da hidrodinâmica sobre os padrões espaciais e temporais da qualidade da água e de florescimentos de cianobactérias e macrófitas em um reservatório urbano

Os reservatórios urbanos estão suscetíveis a uma variedade de interferências antropogênicas que acarretam grande variabilidade espacial e temporal. Contudo, possuem uma dinâmica própria na qual o hidroclima e micro e macro-eventos meteorológicos atuam sobre os processos físicos, químicos e biológicos resultando em respostas particulares de cada corpo de água. No presente estudo a existência de padrões espaciais e temporais na formação de florescimentos de algas, cianobactérias e macrófitas no reservatório Guarapiranga, São Paulo, SP, foi avaliada por meio de experimento de curta escala de tempo durante o evento da entrada de uma frente fria. Foram amostrados 64 pontos em todo o reservatório, e o estudo intensivo de florescimento algal e de cianobactérias em dois ciclos nictemerais, em um ponto selecionado no reservatório. Um modelo tridimensional de hidrodinâmica foi aplicado ao estudo compartimentalizado dos tempos de residência e imagens de satélite foram analisadas para determinação de padrões temporais e espaciais durante períodos de tempo mais amplos. Os resultados revelaram que os períodos mais favoráveis ao surgimento de florescimentos de cianobactérias são geralmente os meses mais quentes, de dezembro e janeiro, ou aqueles em que ocorrem estratificações mais fortes como no fim do inverno, em julho, e após as primeiras chuvas nos meses de setembro e outubro. Existem padrões espaciais recorrentes na formação dos florescimentos, controlados em grande parte pela ação do vento, que no reservatório Guarapiranga é predominantemente nas direções leste e sudeste empurrando os florescimentos na direção da foz dos tributários Embu Mirim e Embu Guaçu e ocasionalmente na direção da foz do rio Parelheiros. As simulações hidrodinâmicas evidenciam as forçantes que determinam os padrões observados e reforçam a importância de se discretizarem os tempos de residência de diferentes compartimentos do reservatório. As séries temporais amplas permitiram a determinação da qualidade da água em cada região e fornecem subsídios para o futuro manejo do reservatório. Como esse comportamento não se restringe ao reservatório Guarapiranga, o tipo de modelagem aqui utilizada pode ser útil para obter informações importantes no processo de planejamento e seleção de medidas para o gerenciamento de reservatórios urbanos tropicais polimíticos, em geral.

Estudo numérico do efeito de blocagem hidrodinâmica em um túnel de cavitação sobre o escoamento ao redor de um submarino.

O projeto de submarinos é um tema que vem se desenvolvendo há muitos anos. No Brasil, vivemos um momento de pleno avanço na área, sendo encabeçado pela Marinha do Brasil com o desenvolvimento do primeiro submarino de propulsão nuclear brasileiro. No decorrer desse projeto, diversos estudos devem ser desenvolvidos, dentre eles, ensaios com modelos em escala, que serão realizados em um túnel de cavitação. Durante estes experimentos, espera-se que surjam diversas fontes de incerteza, que podem comprometer os resultados obtidos, destacando-se o efeito de blocagem hidrodinâmica. Nesse contexto, esta dissertação apresenta uma investigação do fenômeno de blocagem no modelo de submarino DARPA SUBOFF na seção de testes de um túnel de cavitação. Para isso, simulações tridimensionais de média de Reynolds (RANS) foram realizadas para o modelo de submarino posicionado no túnel (sob blocagem) e sob condição de águas profundas (referência), para diferentes ângulos de ataque e número de Reynolds (Re) constante. Coeficientes de força e momento hidrodinâmicos foram obtidos e os métodos de correção de blocagem clássico e WICS foram empregados para obter fatores que pudessem corrigir esses valores. A eficiência de cada método foi avaliada comparando os valores corrigidos com os de referência. Os resultados mostraram que ambos métodos de correção são eficientes na recuperação de coeficientes de força e momento hidrodinâmicos para razões de blocagem inferiores a 10%, mas ressalvas são feitas quanto à utilização de cada um desses métodos, sendo o método clássico indicado para o estudo de geometrias simples e o método WICS para modelos mais complexos. Paralelamente, a topologia do escoamento ao redor do submarino foi estudada através do método dos pontos críticos, para avaliar se a blocagem hidrodinâmica pode causar alterações nas estruturas presentes no escoamento. A análise topológica permitiu avaliar a evolução nas estruturas do escoamento para os ângulos de ataque e Re simulados. Foram capturadas as linhas de separação no corpo e avaliada sua relação com a tensão de cisalhamento ao longo do casco. A blocagem hidrodinâmica não causou alterações drásticas na topologia do escoamento.

Modelação térmica e hidrodinâmica de condensadores para aplicação em sistema de recuperação de energia térmica de veículos automóveis

A recuperação da energia térmica do escape afigura-se como uma vertente com importância relevante para o aumento da eficiência energética dos motores de combustão interna, tendo benefícios em termos económicos e ambientais. O presente trabalho centra-se no estudo e dimensionamento do condensador a integrar num sistema de ciclo de Rankine, tendo-se optado pela seleção de um permutador de placas do tipo “chevron”. O permutador para efeitos de cálculo é dividido em três zonas distintas: i) pré-arrefecedor; ii) condensador e iii) sub-arrefecedor. O presente estudo apresenta uma revisão bibliográfica alargada para o coeficiente de transferência de calor por convecção e perda de carga associada ao escoamento. Para cada zona do permutador são utilizados três fluidos de trabalho distintos: água, etanol e R245fa. Para o fluido refrigerante considerou-se água. Tendo em conta a revisão bibliográfica efetuada, foram implementados modelos para a perda de carga e para a transferência de calor no permutador, atendendo à seleção do fluido de trabalho e ao fluido refrigerante. Os resultados obtidos permitiram verificar que a utilização do fluido de trabalho água, conduz a um condensador mais compacto. Por último é efetuada a otimização do sistema, mediante alteração da temperatura de funcionamento do permutador. Os resultados permitiram verificar a diminuição da temperatura de saída do fluido refrigerante possibilita: i) uma redução do volume do permutador; ii) uma diminuição da perda de carga associada.

Water quality assessment and characterisation of chlorophyll-a variability related to river discharges, within the southeastern Bay of Biscay : evaluation and development of chlorophyll-a algorithms for MODIS and MERIS imagery

222 p. : il.

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.

Distribuição do material particulado em suspensão na plataforma continental e talude superior entre Ubatuba (SP) e o cabo Búzios (RJ)

A pesquisa aqui apresentada amostrou a plataforma continental e talude superior entre Ubatuba (SP) e o cabo Búzios (RJ), em dois projetos de escala local e de meso-escala. A campanha entre Ubatuba e o cabo Frio foi realizada no âmbito do projeto Oceano-Rio: levantamentos oceanográficos integrados ao largo do Estado do Rio de Janeiro, realizado em colaboração com o Ministério de Ciência e Tecnologia (MCT) e com a cooperação da Marinha do Brasil (MB), a bordo do Navio Hidro-oceanográfico Cruzeiro do Sul. Dos dados coletados nesta campanha oceanográfica, foram analisados os dados de CTD, concentração de oxigênio dissolvido, fluorescência, transmissividade, atenuação, retro-espalhamento além de amostras de água para a determinação direta da concentração de material particulado em suspensão (MPS). No experimento de menor escala, entre o cabo Frio e o cabo Búzios, além destes dados, foi realizado um fundeio entre 02 e 30 de julho de 2011, que amostrou a intensidade do eco, as correntes e as ondas. As análises realizadas permitiram observar uma grande abrangência da massa dágua característica da Ressurgência sobre a plataforma continental entre Ubatuba (SP) e cabo Búzios (RJ), inclusive em regiões muito rasas. A ressurgência nem sempre ocorreu associada aos ventos NE, sugerindo que a ocorrência de vórtices tem grande influência sobre o Sistema de Ressurgência de Cabo Frio. O padrão de distribuição de MPS na região é muito influenciado pela ressuspensão dos sedimentos de fundo (aumentada durante as tempestades, quando a coluna dágua passa a apresentar maiores concentrações de MPS), produção fitoplanctônica e aporte das baías costeiras. Outros processos que parecem influenciar este padrão são os vórtices e as ondas internas no talude. A complexidade da interação destas diferentes fontes de MPS e processos associados à sua distribuição ficou evidenciada no processo de conversão dos sinais óticos do transmissômetro e nefelômetro em concentração de MPS. Apesar desta dificuldade associada à grande variabilidade das características do MPS, as análises dos testes de conversão do sinal ótico em concentração de MPS permitem concluir que as melhores conversões podem ser obtidas após a inspeção visual das concentrações de MPS medidos (filtração) e remoção das amostras que não seguem a tendência geral esperada de dispersão sinal ótico versus concentração de MPS.

Circulação e fluxo de material particulado em suspensão no principal canal de acesso à baía de Sepetiba

Entre 06/08/11 e 25/02/12 foram obtidas séries temporais contínuas de intensidade e direção das correntes e intensidade do eco ao longo de toda a coluna dágua, e medições, de temperatura e pressão, próximas ao fundo, na região adjacente ao canal de acesso à baía de Sepetiba (230'16.5"S e 4359'29.4"W, profundidade local de aproximadamente 25 m) através da utilização de um perfilador acústico ADCP (Acoustic Doppler Current Profiler), modelo WorkHorse BroadBand Sentinel (600 kHz, Teledyne-RDI). A partir dos dados adquiridos, observou-se a existência de correntes intensas próximo ao fundo (até 1,02 m/s), principalmente durante os períodos de enchente sob condições de sizígia, que são fortemente influenciadas pela orientação do canal de navegação. A análise das séries temporais das componentes da velocidade mostraram, em conformidade ao relatado por alguns autores, que esta é uma baía cuja circulação é fortemente influenciada pela dinâmica da maré sendo M2, M4, M6 e M8 as principais componentes harmônicas que atuam no sistema. Além disso, observou-se significativa assimetria da maré, sendo os períodos de enchente consideravelmente mais curtos e associados às correntes mais intensas, o que permite concluir que no setor investigado da baía de Sepetiba há dominância de enchente. Outra característica interessante da área de estudo relaciona-se à observação de que ventos intensos de S-SO são responsáveis pelo empilhamento de água no interior da baía, sendo que altura do nível da superfície da água apresenta relação direta com as variações da pressão atmosférica local. Em diversos períodos, quando há atuação de ventos de E-NE é possível encontrar Água Central do Atlântico Sul (ACAS) no interior da baía de Sepetiba, desde a sua entrada principal até as adjacências da ilha Guaíba. Além disso, foram identificadas diversas oscilações de baixa frequência que puderam ser associadas à variação da pressão atmosférica. Oscilações de alta frequência foram associadas à dinâmica da maré, a co-oscilações da maré e ao vento. Em relação à concentração de Material Particulado em Suspensão (MPS), durante os períodos de sizígia foram registradas as maiores concentrações material particulado na coluna dágua. Durante os períodos de sizígia, o fluxo de MPS é mais pronunciado do que durante as quadraturas, sendo que, independente do período da maré, o fluxo cumulativo de material particulado em suspensão é dirigido para o interior da baía de Sepetiba. Considerando o alinhamento dos vetores de velocidade em conformidade a direção preferencial do canal de navegação, tem-se que o fluxo cumulativo de MPS varia entre dirigido para SE na região próxima ao fundo, a dirigido para NE próximo ao topo da coluna dágua, o que sugere a deflexão do movimento das correntes em consequência do atrito. Ao longo de um dia, o fluxo do material particulado é mais pronunciado durante as enchentes, quando há aumento das tensões cisalhantes que atuam sobre o leito da baía redisponibilizando para a coluna dágua o material que estava depositado no fundo.

Previsão do comportamento hidráulico de um rio com base na estimativa de coeficientes que controlam seu escoamento. Estudo de caso: Rio Bengalas, Nova Friburgo-RJ

Nesta tese é realizada a modelagem do comportamento hidráulico dos principais rios que compõem a bacia hidrográfica do Rio Bengalas, localizada no município de Nova Friburgo-RJ, a qual abrange a área mais urbanizada da referida cidade. Para a realização das simulações foi utilizado o Sistema de Modelagem de Águas MOHID, ferramenta MOHID Land. Já para a calibração do modelo foram adotados alguns métodos de otimização, mais precisamente, os algoritmos de Luus- Jaakola (LJ) e Colisão de Partículas (PCA), acoplados ao referido sistema, com o intuito de determinar os principais parâmetros necessários à modelagem de corpos hídricos, bem como suas bacias hidrográficas. Foram utilizados dados topográficos do IBGE disponibilizados pela prefeitura após a elaboração do Plano de Águas Pluviais da região de interesse. Com o modelo devidamente calibrado por meio de dados experimentais, foi realizada a validação do mesmo através da simulação de inundações nesta região. Apesar de técnicas de otimização acopladas à plataforma MOHID terem sido utilizadas pela primeira vez em um rio de montanha, os resultados apresentaram-se importantes e qualitativamente satisfatórios do ponto de vista de auxílio à tomada de decisões, tendo como base a prevenção de danos causados pelas elevações da lâmina dágua que ocorrem frequentemente em Nova Friburgo, como por exemplo, a recente tragédia de janeiro de 2011 ocorrida na Região Serrana do Estado do Rio de Janeiro.

Aplicação do modelo hidrológico SWMM na gestão das águas pluviais urbanas: estudo de caso da bacia hidrográfica do Rio Morto, Rio de Janeiro.

O crescimento da população e dos núcleos urbanos durante o século XX, sobretudo nos países em desenvolvimento, contribuiu para o aumento das áreas impermeáveis das bacias hidrográficas, com impactos importantes nos sistemas de drenagem urbana e na ocorrência de enchentes associadas. As enchentes trazem prejuízos materiais, na saúde e sociais. Recentemente, têm sido propostas práticas conservacionistas e medidas compensatórias, que buscam contribuir para o controle das enchentes urbanas, através do retardo do pico e amortecimento dos hidrogramas. Modelos matemáticos hidrológicos-hidráulicos permitem a simulação da adoção destas medidas de controle, demonstrando e otimizando sua localização. Esta dissertação apresenta os resultados da aplicação do modelo hidrológico Storm Water Management Model (SWMM) à bacia hidrográfica de estudo e representativa do rio Morto localizada em área peri-urbana em Jacarepaguá na cidade do Rio de Janeiro, com área de 9,41 km. O processamento do modelo SWMM foi realizado com o apoio da interface Storm and Sanitary Analysis (SSA), integrada ao sistema AutoCAD Civil 3D. Além da verificação da adequabilidade do modelo à representação dos sistemas hidrológico e hidráulico na bacia, foram desenvolvidos estudos para dois cenários como medidas de controle de enchentes: cenário 1, envolvendo implantação de um reservatório de detenção e, cenário 2, considerando a implantação de reservatórios de águas pluviais nos lotes. Os hidrogramas resultantes foram comparados ao hidrograma resultante da simulação nas condições atuais. Além disso, foram avaliados os custos associados a cada um dos cenários usando o sistema de orçamento da Empresa Rio Águas da PCRJ. Nas simulações foram adotadas a base cartográfica, e os dados climatológicos e hidrológicos previamente observados no contexto do projeto HIDROCIDADES, Rede de Pesquisa BRUM/FINEP, na qual este estudo se insere. Foram representados os processos de geração e propagação do escoamento superficial e de base. Durante o processo de calibração, realizou-se a análise de sensibilidade dos parâmetros, resultando como parâmetros mais sensíveis os relativos às áreas impermeáveis, especialmente o percentual de área impermeável da bacia (Ai). A calibração foi realizada através do ajuste manual de sete parâmetros do escoamento superficial e cinco do escoamento de base para três eventos. Foram obtidos coeficientes de determinação entre 0,52 e 0,64, e a diferença entre os volumes escoados e observados entre 0,60% e 4,96%. Para a validação do modelo foi adotado um evento pluviométrico excepcional observado na cidade em abril de 2010, que à época causou enchentes e grandes transtornos na cidade. Neste caso, o coeficiente de determinação foi igual a 0,78 e a diferença entre volumes foi de 15%. As principais distorções entre hidrogramas observados e simulados foram verificados para as vazões máximas. Em ambos os cenários as enchentes foram controladas. A partir destes estudos, pôde-se concluir que o melhor custo-benefício foi o cenário 2. Para este cenário, foi observado maiores amortecimento e retardo da vazão de pico do hidrograma, igual a 21,51% da vazão simulada para as condições atuais da bacia. Os custos de implantação orçados para os reservatórios de lote ficaram 52% a menos do que o do reservatório de detenção.

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.

Aplicação de múltiplas ferramentas no estudo do transporte de sedimentos na margem interna da restinga da Marambaia (baía de Sepetiba, RJ)

Dados sedimentológicos e hidrodinâmicos foram utilizados para estudar a dinâmica sedimentar da margem interna da restinga da Marambaia (baía de Sepetiba, Rio de Janeiro, SE Brasil). Foram realizadas três campanhas amostrais para coleta de sedimentos ao longo da restinga da Marambaia e na área submersa adjacente. Adicionalmente, em dezembro de 2013, três fundeios simultâneos para aquisição de dados hidrodinâmicos foram conduzidos paralelamente à restinga, ao longo da profundidade média de 6,5 metros. As medições abrangeram dois ciclos de maré (25 horas), incluindo perfilagem com CTDs, perfilagem correntométrica e coleta de amostras de água para determinação das concentrações do material particulado em suspensão. A partir do processamento das amostras de sedimentos em laboratório, foram obtidos os parâmetros estatísticos (média, desvio-padrão e assimetria) para o detalhamento da caracterização sedimentológica da área de estudo, bem como a aplicação de um modelo de tendência direcional de transporte de sedimentos (GisedTrend). A restinga apresenta areias médias e bem selecionadas, enquanto as amostras coletadas na porção submersa, em profundidades maiores que 2 m, passam de areias finas e bem selecionadas a siltes finos e mal selecionados em maiores profundidades. Quatro casos de tendência de transporte foram explorados em um total de 14 possibilidades e os mais significativos foram analisados com o intuito de identificar áreas de acumulação e remoção de sedimentos, podendo estar associados aos processos responsáveis pelo transporte de sedimentos, tais como correntes geradas pelos ventos locais, marés e dispersão de plumas de sedimentos em suspensão. Os dados hidrodinâmicos, durante o período amostrado, permitiram observar o complexo padrão de circulação existente, especialmente na área adjacente à área central da restinga. Dentre as observações realizadas, verificou-se a ocorrência de assimetria da maré, fato já documentado em trabalhos anteriores conduzidos próximo ao principal canal de acesso à baia, com períodos de enchente mais curtos e períodos de vazante mais longos. Diferentemente do que foi reportado nestes trabalhos, os períodos de vazante foram associados às correntes de maior intensidade. As concentrações de material particulado em suspensão situaram-se, em média, entre 10 e 20 mg/L tanto em superfície como próximo ao fundo; no entanto, suas variações temporais não apresentaram qualquer relação com as fases da maré. Possíveis áreas de convergência e divergência da circulação / células de transporte de sedimentos foram identificadas e comparadas aos vetores de tendência de transporte obtidos através do modelo GisedTrend, com resultados satisfatórios.

Uma Metodologia de Estudo de Simulação Tridimensional de Escoamento Turbulento Estratificado no Reservatório de Plantas Hidrelétricas.

Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada.

Experimentos numéricos para deposição de parafinas com modelo Multisólido

Este trabalho objetiva a construção de estruturas robustas e computacionalmente eficientes para a solução do problema de deposição de parafinas do ponto de vista do equilíbrio sólido-líquido. São avaliados diversos modelos termodinâmicos para a fase líquida: equação de estado de Peng-Robinson e os modelos de coeficiente de atividade de Solução Ideal, Wilson, UNIQUAC e UNIFAC. A fase sólida é caracterizada pelo modelo Multisólido. A previsão de formação de fase sólida é inicialmente prevista por um teste de estabilidade termodinâmica. Posteriormente, o sistema de equações não lineares que caracteriza o equilíbrio termodinâmico e as equações de balanço material é resolvido por três abordagens numéricas: método de Newton multivariável, método de Broyden e método Newton-Armijo. Diversos experimentos numéricos foram conduzidos de modo a avaliar os tempos de computação e a robustez frente a diversos cenários de estimativas iniciais dos métodos numéricos para os diferentes modelos e diferentes misturas. Os resultados indicam para a possibilidade de construção de arcabouços computacionais eficientes e robustos, que podem ser empregados acoplados a simuladores de escoamento em dutos, por exemplo.

Inflação em sistemas com fluidos de radiação com viscosidade e dissipação e análise de estabilidade

Este trabalho apresenta um estudo da estabilidade das equações da inflação morna com um fluido de radiação viscoso. A viscosidade do fluido é proveniente do constante decaimento de partículas neste, devido à dissipação do campo escalar da inflação, o ínflaton.Esta viscosidade, que pode ser volumar ou laminar, é tratada em termos de teorias termodinâmicas fora do equilíbrio. Este estudo se limita às equações de fundo da inflação morna, de modo que somente a viscosidade volumar tem um efeito significativo, sendo a viscosidade laminar importante somente no contexto de perturbações cosmológicas. A descrição da viscosidade em termos de uma termodinâmica fora do equilíbrio, porém, não pode ser realizada univocamente, pois a única informação que temos sobre processos irreversíveis é a segunda lei da termodinâmica. Portanto, parte-se em busca de teorias que estejam de acordo com esta lei e que, por argumentos plausíveis, sejam capazes de descrever o comportamento dos fluxos dissipativos próximo ao equilíbrio. O objetivo deste trabalho é estudar a estabilidade da inflação morna viscosa para teorias causais e não causais para o fluido de radiação com viscosidade, de forma que se possa observar o impacto da viscosidade no regime inflacionário e a relevância de se passar a considerar a causalidade. Para o fluido de radiação, as teorias consideradas são a teoria não causal de Eckart e as teorias causais de Israel-Stewart e de Denicol et al (hidrodinâmica dissipativa causal não linear). Obtém-se que as teorias causais, como era de se esperar, além de serem, por definição, consistentes no tocante à finitude da velocidade de propagação dos fluxos dissipativos, tornam o sistema dinâmico estável para valores de viscosidade mais distantes do equilíbrio. Observa-se também, nitidamente, que a teoria de Denicol et al é a mais robusta nesse sentido. Este trabalho, portanto, visa dar continuidade ao estudo dos efeitos não-isentrópicos na inflação, já que, além da dissipação do ínflaton na inflação morna, o impacto da viscosidade tem despertado bastante interesse.

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.