Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu

CMB as a probe of new physics and old times

Cosmic birefringence (CB)---a rotation of photon-polarization plane in vacuum---is a generic signature of new scalar fields that could provide dark energy. Previously, WMAP observations excluded a uniform CB-rotation angle larger than a degree.

In this thesis, we develop a minimum-variance--estimator formalism for reconstructing direction-dependent rotation from full-sky CMB maps, and forecast more than an order-of-magnitude improvement in sensitivity with incoming Planck data and future satellite missions. Next, we perform the first analysis of WMAP-7 data to look for rotation-angle anisotropies and report null detection of the rotation-angle power-spectrum multipoles below L=512, constraining quadrupole amplitude of a scale-invariant power to less than one degree. We further explore the use of a cross-correlation between CMB temperature and the rotation for detecting the CB signal, for different quintessence models. We find that it may improve sensitivity in case of marginal detection, and provide an empirical handle for distinguishing details of new physics indicated by CB.

We then consider other parity-violating physics beyond standard models---in particular, a chiral inflationary-gravitational-wave background. We show that WMAP has no constraining power, while a cosmic-variance--limited experiment would be capable of detecting only a large parity violation. In case of a strong detection of EB/TB correlations, CB can be readily distinguished from chiral gravity waves.

We next adopt our CB analysis to investigate patchy screening of the CMB, driven by inhomogeneities during the Epoch of Reionization (EoR). We constrain a toy model of reionization with WMAP-7 data, and show that data from Planck should start approaching interesting portions of the EoR parameter space and can be used to exclude reionization tomographies with large ionized bubbles.

In light of the upcoming data from low-frequency radio observations of the redshifted 21-cm line from the EoR, we examine probability-distribution functions (PDFs) and difference PDFs of the simulated 21-cm brightness temperature, and discuss the information that can be recovered using these statistics. We find that PDFs are insensitive to details of small-scale physics, but highly sensitive to the properties of the ionizing sources and the size of ionized bubbles.

Finally, we discuss prospects for related future investigations.

Towards understanding the mixing characteristics of turbulent buoyant flows

This work proposes a new simulation methodology in which variable density turbulent flows can be studied in the context of a mixing layer with or without the presence of gravity. Specifically, this methodology is developed to probe the nature of non-buoyantly-driven (i.e. isotropically-driven) or buoyantly-driven mixing deep inside a mixing layer. Numerical forcing methods are incorporated into both the velocity and scalar fields, which extends the length of time over which mixing physics can be studied. The simulation framework is designed to allow for independent variation of four non-dimensional parameters, including the Reynolds, Richardson, Atwood, and Schmidt numbers. Additionally, the governing equations are integrated in such a way to allow for the relative magnitude of buoyant energy production and non-buoyant energy production to be varied.

The computational requirements needed to implement the proposed configuration are presented. They are justified in terms of grid resolution, order of accuracy, and transport scheme. Canonical features of turbulent buoyant flows are reproduced as validation of the proposed methodology. These features include the recovery of isotropic Kolmogorov scales under buoyant and non-buoyant conditions, the recovery of anisotropic one-dimensional energy spectra under buoyant conditions, and the preservation of known statistical distributions in the scalar field, as found in other DNS studies.

This simulation methodology is used to perform a parametric study of turbulent buoyant flows to discern the effects of varying the Reynolds, Richardson, and Atwood numbers on the resulting state of mixing. The effects of the Reynolds and Atwood numbers are isolated by looking at two energy dissipation rate conditions under non-buoyant (variable density) and constant density conditions. The effects of Richardson number are isolated by varying the ratio of buoyant energy production to total energy production from zero (non-buoyant) to one (entirely buoyant) under constant Atwood number, Schmidt number, and energy dissipation rate conditions. It is found that the major differences between non-buoyant and buoyant turbulent flows are contained in the transfer spectrum and longitudinal structure functions, while all other metrics are largely similar (e.g. energy spectra, alignment characteristics of the strain-rate tensor). Also, despite the differences noted between fully buoyant and non-buoyant turbulent fields, the scalar field, in all cases, is unchanged by these. The mixing dynamics in the scalar field are found to be insensitive to the source of turbulent kinetic energy production (non-buoyant vs. buoyant).

Relativistic velocity - potential hydrodynamics and stellar stability

The equations of relativistic, perfect-fluid hydrodynamics are cast in Eulerian form using six scalar "velocity-potential" fields, each of which has an equation of evolution. These equations determine the motion of the fluid through the equation

U_ʋ=µ^-1 (ø,_ʋ + αβ,_ʋ + ƟS,_ʋ).

Einstein's equations and the velocity-potential hydrodynamical equations follow from a variational principle whose action is

I = (R + 16π p) (-g)^1/2 d⁴x,

where R is the scalar curvature of spacetime and p is the pressure of the fluid. These equations are also cast into Hamiltonian form, with Hamiltonian density –T₀⁰ (-g^oo)^-1/2.

The second variation of the action is used as the Lagrangian governing the evolution of small perturbations of differentially rotating stellar models. In Newtonian gravity this leads to linear dynamical stability criteria already known. In general relativity it leads to a new sufficient condition for the stability of such models against arbitrary perturbations.

By introducing three scalar fields defined by

ρ ᵴ = ∇λ + ∇x(xi + ∇xɣi)

(where ᵴ is the vector displacement of the perturbed fluid element, ρ is the mass-density, and i, is an arbitrary vector), the Newtonian stability criteria are greatly simplified for the purpose of practical applications. The relativistic stability criterion is not yet in a form that permits practical calculations, but ways to place it in such a form are discussed.

Produção de oscillons durante restaurações e quebras espontâneas de simetrias

Na natureza há vários fenômenos envolvendo transições de fase com quebra ou restauração de simetrias. Tipicamente, mudanças de fase, são associadas com uma quebra ou restauração de simetria, que acontecem quando um determinado parâmetro de controle é variado, como por exemplo temperatura, densidade, campos externos, ou de forma dinâmica. Essas mudanças que os sistemas sofrem podem levar a formação de defeitos topológicos, tais como paredes de domínios, vórtices ou monopolos magnéticos. Nesse trabalho estudamos particularmente mudanças de fase associadas com quebras ou restaurações dinâmicas de simetria que estão associadas com formação ou destruição de defeitos do tipo paredes de domínio em modelos de campos escalares com simetria discreta. Nesses processos dinâmicos com formação ou destruição de domínios, estudamos a possibilidade de formação de estruturas do tipo oscillons, que são soluções não homogêneas e instáveis de campo, mas que podem concentrar nelas uma quantidade apreciável de energia e terem uma vida (duração) suficientemente grande para serem de importância física. Estudamos a possibilidade de formação dessas soluções em modelos de dois campos escalares interagentes em que o sistema é preparado em diferentes situações, com a dinâmica resultante nesses sistemas estudada numa rede discreta.

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.

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.

Désintégration du faux vide médiée par des kinks en 1+1 dimensions

Dans ce mémoire, on étudie la désintégration d’un faux vide, c’est-à-dire un vide qui est un minimum relatif d’un potentiel scalaire par effet tunnel. Des défauts topologiques en 1+1 dimension, appelés kinks, apparaissent lorsque le potentiel possède un minimum qui brise spontanément une symétrie discrète. En 3+1 dimensions, ces kinks deviennent des murs de domaine. Ils apparaissent par exemple dans les matériaux magnétiques en matière condensée. Un modèle à deux champs scalaires couplés sera étudié ainsi que les solutions aux équations du mouvement qui en découlent. Ce faisant, on analysera comment l’existence et l’énergie des solutions statiques dépend des paramètres du modèle. Un balayage numérique de l’espace des paramètres révèle que les solutions stables se trouvent entre les zones de dissociation, des régions dans l’espace des paramètres où les solutions stables n’existent plus. Le comportement des solutions instables dans les zones de dissociation peut être très différent selon la zone de dissociation dans laquelle une solution se trouve. Le potentiel consiste, dans un premier temps, en un polynôme d’ordre six, auquel on y rajoute, dans un deuxième temps, un polynôme quartique multiplié par un terme de couplage, et est choisi tel que les extrémités du kink soient à des faux vides distincts. Le taux de désintégration a été estimé par une approximation semi-classique pour montrer l’impact des défauts topologiques sur la stabilité du faux vide. Le projet consiste à déterminer les conditions qui permettent aux kinks de catalyser la désintégration du faux vide. Il appert qu’on a trouvé une expression pour déterminer la densité critique de kinks et qu’on comprend ce qui se passe avec la plupart des termes.

Exercises B

Exercises and solutions in LaTex

Exercises B pdf

Exercises and solution in PDF

VACUUM POLARIZATION EFFECTS ON QUASINORMAL MODES IN ELECTRICALLY CHARGED BLACK HOLE SPACE-TIMES

We investigate the influence of vacuum polarization of quantum massive fields on the scalar sector of quasinormal modes in spherically symmetric black holes. We consider the evolution of a massless scalar field on the space-time corresponding to a charged semiclassical black hole, consisting of the quantum-corrected geometry of a Reissner-Nordstrom black hole dressed by a quantum massive scalar field in the large mass limit. Using a sixth order WKB approach we find a shift in the quasinormal mode frequencies due to vacuum polarization.

Condições de energia de hawking e ellis e a equação de raychaudhuri

In the Einstein s theory of General Relativity the field equations relate the geometry of space-time with the content of matter and energy, sources of the gravitational field. This content is described by a second order tensor, known as energy-momentum tensor. On the other hand, the energy-momentum tensors that have physical meaning are not specified by this theory. In the 700s, Hawking and Ellis set a couple of conditions, considered feasible from a physical point of view, in order to limit the arbitrariness of these tensors. These conditions, which became known as Hawking-Ellis energy conditions, play important roles in the gravitation scenario. They are widely used as powerful tools for analysis; from the demonstration of important theorems concerning to the behavior of gravitational fields and geometries associated, the gravity quantum behavior, to the analysis of cosmological models. In this dissertation we present a rigorous deduction of the several energy conditions currently in vogue in the scientific literature, such as: the Null Energy Condition (NEC), Weak Energy Condition (WEC), the Strong Energy Condition (SEC), the Dominant Energy Condition (DEC) and Null Dominant Energy Condition (NDEC). Bearing in mind the most trivial applications in Cosmology and Gravitation, the deductions were initially made for an energy-momentum tensor of a generalized perfect fluid and then extended to scalar fields with minimal and non-minimal coupling to the gravitational field. We also present a study about the possible violations of some of these energy conditions. Aiming the study of the single nature of some exact solutions of Einstein s General Relativity, in 1955 the Indian physicist Raychaudhuri derived an equation that is today considered fundamental to the study of the gravitational attraction of matter, which became known as the Raychaudhuri equation. This famous equation is fundamental for to understanding of gravitational attraction in Astrophysics and Cosmology and for the comprehension of the singularity theorems, such as, the Hawking and Penrose theorem about the singularity of the gravitational collapse. In this dissertation we derive the Raychaudhuri equation, the Frobenius theorem and the Focusing theorem for congruences time-like and null congruences of a pseudo-riemannian manifold. We discuss the geometric and physical meaning of this equation, its connections with the energy conditions, and some of its several aplications.

Two field BPS solutions for generalized Lorentz breaking models

In this work we present nonlinear models in two-dimensional space-time of two interacting scalar fields in the Lorentz and CPT violating scenarios. We discuss the soliton solutions for these models as well as the question of stability for them. This is done by generalizing a model recently published by Barreto and collaborators and also by getting new solutions for the model introduced by them.

Disclosing the generic behavior of topological solutions: An orbit-based approach

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

Nonpolynomial potentials with deformable topological structures

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