Propagation properties of beam generated by Gaussian mirror resonator

Based on scalar diffraction theory, propagation properties of beam generated by Gaussian mirror resonator were investigated. Explicit expressions for the field distribution of the beam that propagate through a paraxial optical ABCD system were derived and numerical examples were given to illustrate our analytical results. (c) 2005 Elsevier B.V. All rights reserved.

Possible antigravity regions in F (R) theory?

We construct an F(R) gravity theory corresponding to the Weyl invariant two scalar field theory. We investigate whether such F (R) gravity can have the antigravity regions where the Weyl curvature invariant does not diverge at the Big Bang and Big Crunch singularities. It is revealed that the divergence cannot be evaded completely but can be much milder than that in the original Weyl invariant two scalar field theory. (C) 2014 The Authors. Published by Elsevier B.V.

电子束蒸发氧化锆薄膜的粗糙度和光散射特性

利用电子束蒸发工艺，以Ag层为衬底，沉积了中心波长为632．8nm的氧化锆（ZrO2）薄膜，膜层厚度在80-480nm范围内变化．研究了不同厚度样品的粗糙度变化规律和表面散射特性．结果发现，随着膜层厚度的逐渐增加．其表面均方根（RMS）粗糙度和总积分散射（TIS）均呈现出先减小后增大的趋势．利用非相关表面粗糙度的散射模型对样品的TIS特性进行了理论计算，所得结果与测量结果相一致．

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

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.

Estudo da quantificação de incertezas para o problema de contaminação de meios porosos heterogêneos

As técnicas de injeção de traçadores têm sido amplamente utilizadas na investigação de escoamentos em meios porosos, principalmente em problemas envolvendo a simulação numérica de escoamentos miscíveis em reservatórios de petróleo e o transporte de contaminantes em aquíferos. Reservatórios subterrâneos são em geral heterogêneos e podem apresentar variações significativas das suas propriedades em várias escalas de comprimento. Estas variações espaciais são incorporadas às equações que governam o escoamento no interior do meio poroso por meio de campos aleatórios. Estes campos podem prover uma descrição das heterogeneidades da formação subterrânea nos casos onde o conhecimento geológico não fornece o detalhamento necessário para a predição determinística do escoamento através do meio poroso. Nesta tese é empregado um modelo lognormal para o campo de permeabilidades a fim de reproduzir-se a distribuição de permeabilidades do meio real, e a geração numérica destes campos aleatórios é feita pelo método da Soma Sucessiva de Campos Gaussianos Independentes (SSCGI). O objetivo principal deste trabalho é o estudo da quantificação de incertezas para o problema inverso do transporte de um traçador em um meio poroso heterogêneo empregando uma abordagem Bayesiana para a atualização dos campos de permeabilidades, baseada na medição dos valores da concentração espacial do traçador em tempos específicos. Um método do tipo Markov Chain Monte Carlo a dois estágios é utilizado na amostragem da distribuição de probabilidade a posteriori e a cadeia de Markov é construída a partir da reconstrução aleatória dos campos de permeabilidades. Na resolução do problema de pressão-velocidade que governa o escoamento empregase um método do tipo Elementos Finitos Mistos adequado para o cálculo acurado dos fluxos em campos de permeabilidades heterogêneos e uma abordagem Lagrangiana, o método Forward Integral Tracking (FIT), é utilizada na simulação numérica do problema do transporte do traçador. Resultados numéricos são obtidos e apresentados para um conjunto de realizações amostrais dos campos de permeabilidades.

Population structure, long-term connectivity, and effective size of mutton snapper (Lutjanus analis) in the Caribbean Sea and Florida Keys

Genetic structure and average long-term connectivity and effective size of mutton snapper (Lutjanus analis) sampled from offshore localities in the U.S. Caribbean and the Florida Keys were assessed by using nuclear-encoded microsatellites and a fragment of mitochondrial DNA. No significant differences in allele, genotype (microsatellites), or haplotype (mtDNA) distributions were detected; tests of selective neutrality (mtDNA) were nonsignificant after Bonferroni correction. Heuristic estimates of average long-term rate of migration (proportion of migrant individuals/generation) between geographically adjacent localities varied from 0.0033 to 0.0054, indicating that local subpopulations could respond independently of environmental perturbations. Estimates of average longterm effective population sizes varied from 341 to 1066 and differed significantly among several of the localities. These results indicate that over time larval drift and interregional adult movement may not be sufficient to maintain population sustainability across the region and that there may be different demographic stocks at some of the localities studied. The estimate of long-term effective population size at the locality offshore of St. Croix was below the minimum threshold size considered necessary to maintain the equilibrium between the loss of adaptive genetic variance from genetic drift and its replacement by mutation. Genetic variability in mutton snapper likely is maintained at the intraregional level by aggregate spawning and random mating of local populations. This feature is perhaps ironic in that aggregate spawning also renders mutton snapper especially vulnerable to overexploitation.

A 4500-year time series of Pacific cod (Gadus macrocephalus) size and abundance: archaeology, oceanic regime shifts, and sustainable fisheries

A 4500-year archaeological record of Pacific cod (Gadus macrocephalus) bones from Sanak Island, Alaska, was used to assess the sustainability of the modern fishery and the effects of this fishery on the size of fish caught. Allometric reconstructions of Pacific cod length for eight prehistoric time periods indicated that the current size of the nearshore, commercially fished Pacific cod stocks is statistically unchanged from that of fish caught during 4500 years of subsistence harvesting. This finding indicates that the current Pacific cod fishery that uses selective harvesting technolog ies is a sustainable commercial fishery. Variation in relative Pacific cod abundances provides further insights into the response of this species to punctuated changes in ocean climate (regime shifts) and indicates that Pacific cod stocks can recover from major environmental perturbations. Such palaeofisheries data can extend the short time-series of fisheries data (<50 yr) that form the basis for fisheries management in the Gulf of Alaska and place current trends within the context of centennial- or millennial-scale patterns.

Hartle's model within the general theory of perturbatiye matchings: the change in mass

Hartle's model provides the most widely used analytic framework to describe isolated compact bodies rotating slowly in equilibrium up to second order in perturbations in the context of General Relativity. Apart from some explicit assumptions, there are some implicit, like the "continuity" of the functions in the perturbed metric across the surface of the body. In this work we sketch the basics for the analysis of the second order problem using the modern theory of perturbed matchings. In particular, the result we present is that when the energy density of the fluid in the static configuration does not vanish at the boundary, one of the functions of the second order perturbation in the setting of the original work by Hartle is not continuous. This discrepancy affects the calculation of the change in mass of the rotating star with respect to the static configuration needed to keep the central energy density unchanged.

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.