Transição à turbulência na camada de mistura estavelmente estratificada utilizando simulação numérica direta e simulação de grandes escalas

A transição à turbulência em uma camada de mistura estavelmente estratificada é de grande interesse para uma variedade de problemas geofísicos e de engenharia. Esta transição é controlada pela competição entre o cisalhamento do escoamento de base e as forças de empuxo, devido à estratificação em densidade do ambiente. Os efeitos do empuxo atuam no escoamento reduzindo a taxa de crescimento das perturbações e retardando a transição à turbulência, enquanto o cisalhamento fornece energia cinética ao escoamento. O presente trabalho investiga a natureza da transição à turbulência em uma camada de mistura temporal estavelmente estratificada através de Simulação Numérica Direta (DNS) e Simulação de Grandes Escalas (LES). O propósito da investigação é analisar o efeito da estratificação estável no desenvolvimento da instabilidade de Kelvin-Helmholtz (K-H) e na formação dos vórtices longitudinais, que se formam após a saturação dos turbilhões primários de K-H. Além deste propósito, é examinado, utilizando de DNS e LES, o desenvolvimento das instabilidades secundárias de K-H na camada baroclínica. Os testes numéricos tridimensionais são realizados com diferentes tipos de condições iniciais para a flutuação de velocidade transversal, enquanto uma condição forçada é usada para as outras duas componentes de flutuação de velocidade. Em particular, o efeito do comprimento transversal do domínio de cálculo é testado empregando diferentes comprimentos, enquanto são usadas as mesmas dimensões para a direção longitudinal e vertical. As simulações bidimensionais mostram que o aumento da estratificação inibe o processo de emparelhamento, reduz a troca de energia entre os turbilhões de K-H e o escoamento, atenua a instabilidade de K-H e diminui o fluxo vertical de massa. A instabilidade secundária do tipo K-H é identtificada na camada baroclínica para Re ¸ 500 quando há o processo de emparelhamento dos vórtices simulados. Na simulação a Re = 500, a instabilidade secundária de K-H aparece tanto para Ri = 0.07 (fraca estratificação) como para Ri = 0.167 (forte estratificação). Os resultados tridimensionais demonstram que os vórtices longitudinais são claramente formados na camada a Ri = 0. Por outro lado, nos casos estratificados os vórtices são enfraquecidos, devido ao gradiente longitudinal de densidade, que diminui a vorticidade nos turbilhões de K-H enquanto aumenta na região entre eles.

Análise numérica de escoamentos turbulentos tridimensionais empregando o método de elementos finitos e simulação de grandes escalas

O presente trabalho apresenta o estudo e implementação de um algoritmo numérico para análise de escoamentos turbulentos, tridimensionais, transientes, incompressíveis e isotérmicos, através da Simulação de Grande Escalas, empregando o Método de Elementos Finitos. A modelagem matemática do problema baseia-se nas equações de conservação de massa e quantidade de movimento de um fluido quase-incompressível. Adota-se um esquema de Taylor-Galerkin, com integração reduzida e fórmulas analíticas das funções de interpolação, para o elemento hexaédrico de oito nós, com funções lineares para as componentes de velocidade e constante no elemento para a pressão. Para abordar o problema da turbulência, emprega-se a Simulação de Grandes Escalas, com modelo para escalas inferiores à resolução da malha. Foram implementados o modelo clássico de Smagorinsky e o modelo dinâmico de viscosidade turbulenta, inicialmente proposto por Germano et al, 1991. Uma nova metodologia, denominada filtragem por elementos finitos independentes, é proposta e empregada, para o processo de segunda filtragem do modelo dinâmico. O esquema, que utiliza elementos finitos independentes envolvendo cada nó da malha original, apresentou bons resultados com um baixo custo computacional adicional. São apresentados resultados para problemas clássicos, que demonstram a validade do sistema desenvolvido. A aplicabilidade do esquema utilizado, para análise de escoamentos caracterizados por elevados números de Reynolds, é discutida no capítulo final. São apresentadas sugestões para aprimorar o esquema, visando superar as dificuldades encontradas com respeito ao tempo total de processamento, para análise de escoamentos tridimensionais, turbulentos e transientes .

Simulação de grandes escalas de escoamentos incompressíveis com transferência de calor e massa por um método de elementos finitos de subdomínio

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Análise de escoamentos incompressíveis utilizando simulação de grandes escalas e adaptação de malhas

No presente estudo, são apresentadas soluções numéricas de problemas de Engenharia, na área de Dinâmica dos Fluidos Computacional, envolvendo fluidos viscosos, em escoamentos incompressíveis, isotérmicos e não isotérmicos, em regime laminar e turbulento, podendo envolver transporte de massa. Os principais objetivos deste trabalho são a formulação e a aplicação de uma estratégia de adaptação automática de malhas e a inclusão de modelos de viscosidade turbulenta, integrados com um algoritmo utilizado para simular escoamentos de fluidos viscosos bi e tridimensionais, no contexto de malhas não estruturadas. O estudo é dirigido no sentido de aumentar o conhecimento a respeito das estruturas de escoamentos turbulentos e de estudar os efeitos físicos no transporte de quantidades escalares propiciando, através de técnicas de adaptação automática de malhas, a obtenção de soluções numéricas precisas a um custo computacional otimizado. As equações de conservação de massa, de balanço de quantidade de movimento e de quantidade escalar filtradas são utilizadas para simular as grandes escalas de escoamentos turbulentos e, para representar as escalas submalha, são utilizados dois modelos de viscosidade turbulenta: o modelo de Smagorinsky clássico e o modelo dinâmico. Para obter soluções numéricas com precisão, é desenvolvida e implementada uma estratégia de adaptação automática de malhas, a qual é realizada simultaneamente e interativamente com a obtenção da solução. O estudo do comportamento da solução numérica é fundamentado em indicadores de erro, com o propósito de mapear as regiões onde certos fenômenos físicos do escoamento ocorrem com maior intensidade e de aplicar nestas regiões um esquema de adaptação de malhas. A adaptação é constituída por processos de refinamento/desrefinamento e por um processo de suavização laplaciana. Os procedimentos para a implementação dos modelos de viscosidade turbulenta e a estratégia de adaptação automática de malhas são incorporados ao código computacional de elementos finitos tridimensionais, o qual utiliza elementos tetraédricos lineares. Aplicações de escoamentos de fluidos viscosos, incompressíveis, isotérmicos e não isotérmicos em regime laminar e turbulento são simuladas e os resultados são apresentados e comparados com os obtidos numérica ou experimentalmente por outros autores.

Simulação numérica de escoamentos de fluidos pelo método de elementos finitos baseado em volumes de controle em malhas não estruturadas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação numérica de escoamentos de fluidos com efeitos de rotação

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Simulação numérica do escoamento de jatos impingentes axissimétricos confinados

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Simulação numérica de escoamentos de fluidos pelo método de elementos finitos de mínimos quadrados

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Simulação numérica de escoamentos de fluídos pelo método de elementos finitos baseado em volumes de controle com a técnica de passo fracionado

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Efeito das escalas de simulação sobre os padrões deposicionais de fluxos turbidíticos

Este trabalho apresenta simulações físicas de correntes de densidade não-conservativas desenvolvendo-se num canal bidimensional de grande porte (5,00x0,40x1,00m) e declividade variável. A metodologia empregada foi desenvolvida em estudos anteriores (Manica, 2002 e Ávila, 2003) para ensaios em canal de pequeno porte e adaptada para o canal de maior porte, objetivando gerar subsídios para o estudo dos efeitos da escala de simulação sobre os padrões deposicionais de correntes de turbidez. As correntes foram formadas por uma mistura de água e carvão com massa específica em torno de 1023 kg/m3. A granulometria dos carvões utilizados ficou entre 0,062-0,297mm para o Carvão 207 e entre 0,062-0,177mm para o Carvão 205. Os experimentos foram divididos em três etapas. Na primeira etapa, foram simuladas correntes simples com o objetivo de realizar a transposição da metodologia de simulação física de correntes de turbidez do canal de pequeno porte para o canal de grande porte. Os dados obtidos foram comparados com os resultados pré-existentes do canal de pequeno porte provenientes dos estudos de Ávila (2003). Na segunda etapa foram simuladas correntes consecutivas, as quais foram compostas de um total de quatro fluxos por ensaio. Novamente, os dados obtidos foram comparados com os resultados pré-existentes dos estudos de Ávila (2003). Finalmente, na terceira etapa, foram simuladas correntes simples sobre canal com declividades de –0,5° e 2,0°. As correntes 00000geradas durante a segunda e a terceira etapas também tiveram seus parâmetros geométricos (alturas da cabeça e do corpo) e hidrodinâmicos (velocidades e acelerações da cabeça) comparados com dados pré-existentes de simulações com correntes conservativas (Fabian, 2002) Os resultados obtidos demonstraram que os depósitos gerados durante a primeira etapa foram os que mais se assemelharam aos dados pré-existentes de Ávila (2003) no que se refere às distribuições longitudinais e verticais. Também se observou que as correntes nãoconservativas não apresentaram variações significativas das alturas da cabeça e do corpo para as declividades ensaiadas. Tais observações acerca dos depósitos e da hidrodinâmica sugerem que, para que se garanta a eliminação de efeitos de escala, devem ser mantidas constantes: a massa específica do sedimento empregado na mistura, as condições de admissão da mesma e também o valor do número de Froude Densimétrico. Também há indícios de que variações nos parâmetros geométricos das correntes podem estar associadas com a massa específica do sedimento empregado para a simulação e com as condições de admissão da mistura.

Inflação estocástica não-isentrópica

Em modelos inflacionários não-isentrópicos, a contribuição para o espectro de potência é essencialmente proveniente das flutuações térmicas. Esta é a situação oposta a da inflação fria, onde as flutuações de origem quântica fornecem toda contribuição para o espectro. Pouca ou nenhuma importância tem sido dada ao regime intermediário, onde as flutuações quânticas e térmicas são comparáveis. Neste trabalho, tendo como bases a inflação não-isentrópica e a inflação estocástica de Starobinsky, propomos um quadro geral onde é possível tratar de maneira conjunta, explícita e transparente tanto a contribuição de origem quântica quanto a de origem térmica para o espectro de potência do inflaton.O espectro de potência geral obtido reproduz, nos limites apropriados, todos os resultados caracteríssticos tanto da inflação fria, quanto da inflação não-isentrópica. Com o objetivo de checar a consistência e a viabilidade do modelo, foram usados os típicos potenciais polinomiais característicos da inflação caótica. Apesar destes potenciais já estarem praticamente descartados pelas observações no contexto da inflação fria, surpreendentemente pudemos constatar que efeitos dissipativos e de temperatura são capazes de restaurar a compatibilidade dos mesmos com os parâmetros cosmológicos inferidos através dos dados do nono ano do WMAP. Através da inserção de tais efeitos na dinâmica de grandes escalas do inflaton, estendemos ainda alguns resultados relacionados ao cenário conhecido como inflação eterna.

Análise da distribuição de um simulador multinível

A validação de projetos de sistemas eletrônicos pode ser feita de diversas maneiras, como tem sido mostrado pelas pesquisas em síntese automática e verificação formal. Porém, a simulação ainda é o método mais utilizado. O projeto de um sistema digital típico pode ser desenvolvido em diversos níveis de abstração, como os níveis algorítmico, lógico ou analógico. Assim, a simulação também deve ser executada em todos esses níveis. A simulação apresenta, contudo, o inconveniente de não conseguir conciliar uma alta acuracidade de resultados com um baixo tempo de simulação. Quanto mais detalhada é a descrição do circuito, maior é o tempo necessário para simulá-lo. O inverso também é verdadeiro, ou seja, quanto menor for a acuracidade exigida, menor será o tempo dispendido. A simulação multinível tenta conciliar eficiêencia e acuracidade na simulação de circuitos digitais, propondo que partes do circuito sejam descritas em diferentes níveis de abstração. Com isso, somente as partes mais críticas do sistema são descritas em detalhes e a velocidade da simulação aumenta. Contudo, essa abordagem não é suficiente para oferecer um grande aumento na velocidade de simulação de grandes circuitos. Assim, surge como alternativa a aplicação de técnicas de processamento distribuído à simulação multinível. Os aspectos que envolvem a combinação dessas duas técnicas são abordados nesse trabalho. Como plataforma para os estudos realizados, optou-se por usar duas ferramentas desenvolvidas nessa Universidade: os simuladores do Sistema AMPLO e o Sistema Operacional HetNOS. São estudadas técnicas de sincronização em sistemas distribuídos, fundamentais para o desenvolvimento dos simuladores e, finalmente, são propostas alternativas para a distribuição dos simuladores. É realizada, ainda, uma análise comparativa entre as versões propostas.

Anisotropias da radiação cósmica de fundo e vínculos em modelos com decaimento do vácuo

Many astronomical observations in the last few years are strongly suggesting that the current Universe is spatially flat and dominated by an exotic form of energy. This unknown energy density accelerates the universe expansion and corresponds to around 70% of its total density being usually called Dark Energy or Quintessence. One of the candidates to dark energy is the so-called cosmological constant (Λ) which is usually interpreted as the vacuum energy density. However, in order to remove the discrepancy between the expected and observed values for the vacuum energy density some current models assume that the vacuum energy is continuously decaying due to its possible coupling with the others matter fields existing in the Cosmos. In this dissertation, starting from concepts and basis of General Relativity Theory, we study the Cosmic Microwave Background Radiation with emphasis on the anisotropies or temperature fluctuations which are one of the oldest relic of the observed Universe. The anisotropies are deduced by integrating the Boltzmann equation in order to explain qualitatively the generation and c1assification of the fluctuations. In the following we construct explicitly the angular power spectrum of anisotropies for cosmologies with cosmological constant (ΛCDM) and a decaying vacuum energy density (Λ(t)CDM). Finally, with basis on the quadrupole moment measured by the WMAP experiment, we estimate the decaying rates of the vacuum energy density in matter and in radiation for a smoothly and non-smoothly decaying vacuum

El potencial de servicios multimodales en el corredor de transporte Asunción - Montevideo

Incluye Bibliografía

Energia escura acoplada

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)