Propagação de ondas em teorias alternativas da gravitação

Uma forma de generalizar a teoria de Einstein da gravitação é incorporar na lagrangiana termos que dependem de escalares formados com os tensores de Ricci e Riemann, tais como (Ricci)2, ou (Riemann)2. Estas teorias tem sido estudadas intensamente nos últimos anos, já que elas podem ser usadas para descrever a expansão acelerada do universo no modelo cosmológico standard. Entre os desfios de modificar a teoria de Einstein, se encontra o de limitar a ambiguidade na escolha da dependência da lagrangiana com os escalares antes mencionados. A proposta desta dissertação é a de colocar limites sobre as possíveis lagrangianas impondo que as ondas (isto é, perturbações lineares) se propaguem no vácuo sem que apareça, shocks.

O pré-aquecimento do universo: da ressonância paramétrica à turbulência

A inflação consegue dar conta de uma série de problemas do Modelo padrão da Cosmologia, preservando ainda o sucesso do modelo do Big Bang. Na sua versão mais simples, a inflação é controlada por um campo escalar, o ínflaton, que faz com que o universo se expanda exponencialmente. Após, o ínflaton decai e ocorre o reaquecimento do universo. Contudo, alguns autores apontam a existência de uma fase intermediária, chamada de pré-aquecimento. O decaimento do ínflaton possui uma rica dinâmica não-linear. No primeiro estágio, a ressonância paramétrica promove o crescimento exponencial de alguns modos do ínflaton. Isto altera a dinâmica do modo homogêneo do ínflaton, promovendo uma reestruturação das cartas de ressonâncias da equação de movimento dos modos perturbativos. Desta forma, ocorre a transferência de energia para estes modos, até que o universo termaliza. Esta transferência de energia é típica de um sistema turbulento. Por se tratar de uma evolução não-linear, torna-se conveniente a implementação computacional de métodos numéricos. Neste contexto, os métodos espectrais têm se mostrado uma excelente ferramenta para estudar este tipo de sistema. Esta dissertação apresenta os resultados do esquema numérico desenvolvido para o modelo com potencial quártico, que será a base para os demais estudos a serem desenvolvidos. Como mostrado, este esquema é extremamente preciso e eficiente.

Implementação paralela de um código de elementos finitos em 2D para as Equações de Navier-Stokes para fluidos incompressíveis com transporte de escalares.

O estudo do fluxo de água e do transporte escalar em reservatórios hidrelétricos é importante para a determinação da qualidade da água durante as fases iniciais do enchimento e durante a vida útil do reservatório. Neste contexto, um código de elementos finitos paralelo 2D foi implementado para resolver as equações de Navier-Stokes para fluido incompressível acopladas a transporte escalar, utilizando o modelo de programação de troca de mensagens, a fim de realizar simulações em um ambiente de cluster de computadores. A discretização espacial é baseada no elemento MINI, que satisfaz as condições de Babuska-Brezzi (BB), que permite uma formulação mista estável. Todas as estruturas de dados distribuídos necessárias nas diferentes fases do código, como pré-processamento, solução e pós-processamento, foram implementadas usando a biblioteca PETSc. Os sistemas lineares resultantes foram resolvidos usando o método da projeção discreto com fatoração LU por blocos. Para aumentar o desempenho paralelo na solução dos sistemas lineares, foi empregado o método de condensação estática para resolver a velocidade intermediária nos vértices e no centróide do elemento MINI separadamente. Os resultados de desempenho do método de condensação estática com a abordagem da solução do sistema completo foram comparados. Os testes mostraram que o método de condensação estática apresenta melhor desempenho para grandes problemas, às custas de maior uso de memória. O desempenho de outras partes do código também são apresentados.

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.

Estudo de níveis de ozônio troposférico usando métodos de estatística univariada e multivariada para duas localidades da cidade do Rio de Janeiro

Com cada vez mais intenso desenvolvimento urbano e industrial, atualmente um desafio fundamental é eliminar ou reduzir o impacto causado pelas emissões de poluentes para a atmosfera. No ano de 2012, o Rio de Janeiro sediou a Rio +20, a Conferência das Nações Unidas sobre Desenvolvimento Sustentável, onde representantes de todo o mundo participaram. Na época, entre outros assuntos foram discutidos a economia verde e o desenvolvimento sustentável. O O3 troposférico apresenta-se como uma variável extremamente importante devido ao seu forte impacto ambiental, e conhecer o comportamento dos parâmetros que afetam a qualidade do ar de uma região, é útil para prever cenários. A química das ciências atmosféricas e meteorologia são altamente não lineares e, assim, as previsões de parâmetros de qualidade do ar são difíceis de serem determinadas. A qualidade do ar depende de emissões, de meteorologia e topografia. Os dados observados foram o dióxido de nitrogênio (NO2), monóxido de nitrogênio (NO), óxidos de nitrogênio (NOx), monóxido de carbono (CO), ozônio (O3), velocidade escalar vento (VEV), radiação solar global (RSG), temperatura (TEM), umidade relativa (UR) e foram coletados através da estação móvel de monitoramento da Secretaria do Meio Ambiente (SMAC) do Rio de Janeiro em dois locais na área metropolitana, na Pontifícia Universidade Católica (PUC-Rio) e na Universidade do Estado do Rio de Janeiro (UERJ) no ano de 2011 e 2012. Este estudo teve três objetivos: (1) analisar o comportamento das variáveis, utilizando o método de análise de componentes principais (PCA) de análise exploratória, (2) propor previsões de níveis de O3 a partir de poluentes primários e de fatores meteorológicos, comparando a eficácia dos métodos não lineares, como as redes neurais artificiais (ANN) e regressão por máquina de vetor de suporte (SVM-R), a partir de poluentes primários e de fatores meteorológicos e, finalmente, (3) realizar método de classificação de dados usando a classificação por máquina de vetor suporte (SVM-C). A técnica PCA mostrou que, para conjunto de dados da PUC as variáveis NO, NOx e VEV obtiveram um impacto maior sobre a concentração de O3 e o conjunto de dados da UERJ teve a TEM e a RSG como as variáveis mais importantes. Os resultados das técnicas de regressão não linear ANN e SVM obtidos foram muito próximos e aceitáveis para o conjunto de dados da UERJ apresentando coeficiente de determinação (R2) para a validação, 0,9122 e 0,9152 e Raiz Quadrada do Erro Médio Quadrático (RMECV) 7,66 e 7,85, respectivamente. Quanto aos conjuntos de dados PUC e PUC+UERJ, ambas as técnicas, obtiveram resultados menos satisfatórios. Para estes conjuntos de dados, a SVM mostrou resultados ligeiramente superiores, e PCA, SVM e ANN demonstraram sua robustez apresentando-se como ferramentas úteis para a compreensão, classificação e previsão de cenários da qualidade do ar

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.

Dinâmica da condensação de Bose-Einstein em gases fracamente interagentes

A presente dissertação estuda com detalhes a evolução temporal fora do equilíbrio de um condensado de Bose-Einstein homogêneo diluído imerso em um reservatório térmico. Nós modelamos o sistema através de um campo de Bose escalar complexo. É apropriado descrever o comportamento microscópico desse sistema por meio da teoria quântica de campos através do formalismo de Schwinger-Keldysh. Usando esse formalismo, de tempo real a dinâmica do condensado é solucionada por um grupo de equações integro-diferencial auto consistente, essas são solucionadas numericamente. Estudamos também o cenário quench, e como a densidade do gás e as interações entre as flutuações tem o efeito de provocar as instabilidades nesse sistema. Aplicamos esse desenvolvimento para estudar o comportamento de duas espécies homogêneas de um gás de Bose diluído imerso em um reservatório térmico.

Modelo de Ginzburg-Landau a partir da teoria de campos a temperatura finita

Neste trabalho, utilizamos o formalismo de teorias quânticas de campos a temperatura finita, tal como desenvolvidas por Matsubara, aplicado a uma hamiltoniana de N campos escalares com autointeração quártica a N grande. Obtém-se uma expressão, na primeira aproximação quântica, para o coeficiente do termo quadrático da hamiltoniana ("massa quadrada"), renormalizado, como função da temperatura. A partir dela, estudamos o processo de quebra espontânea de simetria. Por outro lado, a mesma hamiltoniana é conhecida como modelo de Ginzburg-Landau na literatura de matéria condensada, e que permite o estudo de transições de fase em materiais ferromagnéticos. A temperatura é introduzida através do termo quadrático na hamiltoniana, de forma linear: é proporcional à diferença entre a variável de temperatura e a temperatura crítica. Tal modelo, porém, possui validade apenas na regi~ao de temperaturas próximas à criticalidade. Como resultado de nossos cálculos na teoria de campos a temperatura finita, observamos que, numa faixa de valores em torno da temperatura crítica, a massa quadrática pode ser aproximada por uma relação linear em relação à variável de temperatura. Isso evidencia a compatibilidade da abordagem de Ginzburg-Landau, na vizinhança da criticalidade, com respeito ao formalismo de campos a temperatura finita. Discutimos também os efeitos causados pela presença de um potencial químico no sistema.

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.

Reducing energy demand: what are the practical limits?

Concern over the global energy system, whether driven by climate change, national security, or fears of shortage, is being discussed widely and in every arena but with a bias toward energy supply options. While demand reduction is often mentioned in passing, it is rarely a priority for implementation, whether through policy or through the search for innovation. This paper aims to draw attention to the opportunity for major reduction in energy demand, by presenting an analysis of how much of current global energy demand could be avoided. Previous work led to a "map" of global energy use that traces the flow of energy from primary sources (fuels or renewable sources), through fuel refinery, electricity generation, and end-use conversion devices, to passive systems and the delivery of final energy services (transport, illumination, and sustenance). The key passive systems are presented here and analyzed through simple engineering models with scalar equations using data based on current global practice. Physically credible options for change to key design parameters are identified and used to predict the energy savings possible for each system. The result demonstrates that 73% of global energy use could be saved by practically achievable design changes to passive systems. This reduction could be increased by further efficiency improvements in conversion devices. A list of the solutions required to achieve these savings is provided.

Experiments in constituting a taste panel for canned and frozen prawn products

With a view to constituting a taste panel in the laboratory for detecting the flavor changes in canned and frozen prawn, three methods of panel selection (scalar scoring method, range and deviation method and triangular method) were tried. Out of the three, triangular method was found to be suitable for panel formation in canned and frozen prawn. Using this method a panel of six members was formed for detecting flavor changes in the two products.

Investigation of the losses in the stator duct spacers of large A.C. motors using the A -ip 3-D eddy current formulation

This paper describes how the A -if) formulation may be applied to determine the losses in the stator duct spacers of large a.c. motors. The model is described in terms of its geometry and boundary conditions. The novel aspects of the application of the formulation to this problem are explained. These include the modelling of fixed currents sources (the stator windings), the location of the necessary cut surfaces and the determination of their magnetic scalar potential differences, and the implementation of periodic boundary conditions for vector variables. Results are presented showing how the duct spacer losses vary with load, and with the relative permeability of the spacer material. The effects of modelling iron nonlinearity, of both the spacer and the steel laminations, are also illustrated. © 1996 IEEE.

Secondary flows and loss caused by blade row interaction in a turbine stage

A study of the three-dimensional stator-rotor interaction in a turbine stage is presented. Experimental data reveal vortices downstream of the rotor which are stationary in the absolute frame - indicating that they are caused by the stator exit flowfield. Evidence of the rotor hub passage vortices is seen, but additional vortical structures away from the endwalls, which would not be present if the rotor were tested in isolation, are also identified. An unsteady computation of the rotor row is performed using the measured stator exit flowfield as the inlet boundary condition. The strength and location of the vortices at rotor exit are predicted. A formation mechanism is proposed whereby stator wake fluid with steep spanwise gradients of absolute total pressure is responsible for all but one of the rotor exit vortices. This mechanism is then verified computationally using a passive-scalar tracking technique. The predicted loss generation through the rotor row is then presented and a comparison made with a steady calculation where the inlet flow has been mixed out to pitchwise uniformity. The loss produced in the steady simulation, even allowing for the mixing loss at inlet, is 10% less than that produced in the unsteady simulation. This difference highlights the importance of the time-accurate calculation as a tool of the turbomachine designer.

Three-dimensional direct numerical simulations of autoignition in turbulent non-premixed flows with simple and complex chemistry

3D Direct Numerical Simulations (DNS) of autoignition in turbulent non-premixed flows between fuel and hotter air have been carried out using both 1-step and complex chemistry consisting of a 22 species n-heptane mechanism to investigate spontaneous ignition timing and location. The simple chemistry results showed that the previous findings from 2D DNS that ignition occurred at the most reactive mixture fraction (ξMR) and at small values of the conditional scalar dissipation rate (N|ξMR) are valid also for 3D turbulent mixing fields. Performing the same simulation many times with different realizations of the initial velocity field resulted in a very narrow statistical distribution of ignition delay time, consistent with a previous conjecture that the first appearance of ignition is correlated with the low-N content of the conditional probability density function of N. The simulations with complex chemistry for conditions outside the Negative Temperature Coefficient (NTC) regime show behaviour similar to the single-step chemistry simulations. However, in the NTC regime, the most reactive mixture fraction is very rich and ignition seems to occur at high values of scalar dissipation. Copyright © 2006 by ASME.

An evaluation of volumetric interest points

This paper presents the first performance evaluation of interest points on scalar volumetric data. Such data encodes 3D shape, a fundamental property of objects. The use of another such property, texture (i.e. 2D surface colouration), or appearance, for object detection, recognition and registration has been well studied; 3D shape less so. However, the increasing prevalence of depth sensors and the diminishing returns to be had from appearance alone have seen a surge in shape-based methods. In this work we investigate the performance of several detectors of interest points in volumetric data, in terms of repeatability, number and nature of interest points. Such methods form the first step in many shape-based applications. Our detailed comparison, with both quantitative and qualitative measures on synthetic and real 3D data, both point-based and volumetric, aids readers in selecting a method suitable for their application. © 2011 IEEE.