Estudo da difusão translacional em um modelo para água

Apesar da água ser o líquido mais comum na natureza, suas características ainda não estão totalmente explicadas. A relação entre a forma do potencial intermolecular efetivo que representa as interações presentes na água, as várias anomalias existentes na água e a possível existência de dupla criticalidade ainda é uma questão em aberto. Recentemente descobriu-se que a água apresenta, além de anomalias termodinâmicas, anomalia na difusão translacional e rotacional. Mostrou-se que para um modelo SPC/E para água estas anomalias dinâmicas estão conectadas com a temperatura de máxima densidade (TMD) e que as anomalias dinâmicas ocupam uma região maior que a TMD no diagrama de fases p vs. T. Nesta dissertação investigamos a relação entre a anomalia na densidade e a anomalia na difusão translacional da água em um modelo onde as moléculas interagem através de um potencial que contém duas fontes de interação: uma isotrópica atrativa do tipo van der Waals, v, e uma direcional, u, que pode ser nula ou repulsiva. Reproduzimos o diagrama de fases pressão vs. temperatura aonde a anomalia na densidade e a dupla criticalidade se fazem presente. Mostramos que este sistema apresenta anomalia na difusão translacional e que esta se encontra na mesma região do diagrama de fases pressão vs. temperatura em que a TMD está presente. Diferentemente do que ocorre para água SPC/E, a região de anomalia na difusão translacional ocupa um espaço menor, encontrando-se na parte interna da TMD. Sugerimos que a discretização do potencial afeta mais fortemente a anomalia dinâmica do que a termodinâmica e que leva a um encolhimento na região de anomalia na difusão.

Ferramentas para análise dinâmica e estrutural de um reboque de linha leve

Neste trabalho é desenvolvida uma metodologia de projeto para identificar as regiões críticas da estrutura de um reboque de linha leve sendo tracionado em pavimentos do tipo rodovia de baixa qualidade e estrada secundária de muito baixa qualidade. Para tanto, são levantados alguns dados experimentais da estrutura, necessários para a aproximação e simulação dinâmica de um modelo simplificado. A excitação da base é realizada por atuadores que simulam as oscilações verticais de um perfil de estrada, a qual é definida de acordo com os estudos realizados por Dodds e Robson (1973). Isto permite a determinação de um histórico de carregamentos das regiões da estrutura do chassi sob a ação das molas da suspensão. Em seguida, é gerado um modelo estrutural simplificado do reboque em elementos finitos, chamado de global, no qual são determinadas as regiões sob ação das maiores tensões. Tendo identificada a região mais crítica da estrutura, é criado um modelo local desta parte, onde se pode observar a distribuição de tensões com mais detalhe, permitindo a identificação dos pontos de concentração de tensões. Desta forma, com a aplicação do método de análise global-local é possível a obtenção de resultados detalhados quanto aos esforços da estrutura com um menor custo computacional.

Influência de parâmetros moleculares em funções de correlação temporal na dinâmica de solvatação mecânica

No presente trabalho descrevemos nossos resultados relativos à investigação da dinâmica de solvatação mecânica por meio de simulações por dinâmica molecular, respeitando o regime da resposta linear, em sistemas-modelo de argônio líquido com um soluto monoatômico ou diatômico dissolvido. Estudamos sistematicamente a influência dos parâmetros moleculares dos solutos (tamanho, polarizabilidade) e da densidade frente a vários modelos de solvatação. Funções de Correlação Temporal da Energia de Solvatação foram calculadas com relação à correlações de n-corpos (n = 2; 3) distinguindo interações repulsivas e atrativas para ambos os sistemas líquidos. Também obtivemos segundas derivadas temporais dessas funções referindo-se à parcelas translacionais, rotacionais e roto-translacionais na solução do diatômico. Encontramos que funções de correlação temporal coletivas podem ser razoavelmente bem aproximadas por correlações binárias a densidades baixas e, a densidades altas, correlações ternárias tornam-se mais importantes produzindo um descorrelacionamento mais rápido das funções coletivas devido a efeitos de cancelamento parciais. As funções de correlação para interações repulsivas e atrativas exibem comportamentos dinâmicos independentes do modelo de solvatação devido a fatores de escalonamento linear que afetam apenas as amplitudes das dessas funções de correlação temporal. Em geral, os sistemas com grau de liberdade rotacional apresentam tempos de correlação mais curtos para a dinâmica coletiva e tempos de correlação mais longos para as funções binárias e ternárias. Finalmente, esse estudo mostra que os sistemas contendo o diatômico relaxam-se predominantemente por mecanismos translacionais binários em modelos de solvatação envolvendo alterações apenas na polarizabilidade do soluto, e por mecanismos rotacionais atrativos binários em modelos envolvendo alterações no comprimento de ligação.

Um estudo da otimização da geometria de um pára-quedas simplificado

Neste trabalho desenvolve-se um estudo numérico do fluxo de ar em torno da geometria de um pára-quedas tradicional simplificado, para alguns valores de Reynolds. O método baseia-se na solução das equações incompressíveis de Navier- Stokes discretizadas pelo método de diferenças finitas e integradas pelo método de Runge-Kutta. Utiliza-se o método dos contornos virtuais para representar a geometria numa malha cartesiana e o método de otimização não-linear dos poliedros flexíveis para otimização do coeficiente de arraste calculado através do código de dinâmica de fluidos computacional; esteé um método de busca multivariável, onde o pior vértice de um poliedro com n + 1 vérticesé substituído por um novo.

Estudo teórico da solvatação de ânions em solventes orgânicos

O presente trabalho relata a parametrização do modelo de contínuo IEF-PCM para o estudo da solvatação de ânions monovalentes em acetonitrila e N,N-dimetilformamida. Nós propomos duas metodologias de parametrização do modelo de contínuo polarizável para este estudo teórico. A primeira proposta é baseada unicamente no cálculo do termo eletrostático da energia livre e na equiparação deste termo com a energia livre de solvatação experimental. Na segunda proposta de metodologia, a parametrização do modelo de contínuo polarizável foi realizada com o auxílio das simulações de dinâmica molecular e dentro desta segunda proposta nós calculamos todos os termos que contribuem para a energia livre de solvatação, a energia de cavitação. a energia eletrostática e a energia de van der Waals. Com a aplicação da primeira metodologia, é obtido um erro médio absoluto de 2,1 kcal/mol nas energias de solvatação dos ânions em acetonitrila enquanto que o erro médio absoluto calculado para energias de solvatação dos ânions em N,N-dimetilformamida é 2,8 kcal/mol. A aplicação da segunda metodologia resulta num erro médio absoluto de 1,4 kcal/mol nas energias de solvatação dos ânions em acetonitrila e um erro médio absoluto de 1,2 kcal/mol nas energias de solvatação calculadas para os ânions em N,N-dimetilformamida. Estes resultados indicam que o modelo de contínuo polarizável IEF-PCM é adequado para o estudo dos efeitos termodinâmicos nas soluções de ânions nos solventes acetonitrila e N,N-dimetilformamida.

Compressibilidade da matéria nuclear em estrelas de nêutrons

Neste trabalho, são discutidos modelos da hadrodinâmica quântica com aproximação de campo médio aplicados a estrelas de nêutrons. O modelo de Walecka define o ponto de partida para desenvolver o modelo de acoplamento derivativo ajustável. A presente dissertação visa a um estudo detalhado sobre a influência dos parâmetros do modelo ajustável no sistema, analisando seus limites, inclusive quando os parâmetros são iguais a zero ou infinito (modelo exponencial). Esta análise tem o propósito de estabelecer um conjunto de parâmetros que defina um modelo que esteja de acordo com as propriedades fenomenológicas tais como módulo de compressão da matéria nuclear, massa efetiva na saturação da matéria nuclear e também algumas propriedades estáticas globais das estrelas de nêutrons como, por exemplo, massa e raio. Estabelecido o modelo a ser considerado, a autora dessa dissertação introduz, como inovação, a compressibilidade hadrônica como função da densidade. Tradicionalmente, determinam-se propriedades da matéria apenas para a densidade de saturação.

Espectroscopia resolvida no tempo : caracterização de pulsos curtos, dinâmica molecular em líquidos, modelamento de luz incoerente

Nesta tese realizamos três trabalhos que envolvem diferentes aspectos da espectroscopia resolvida no tempo. No primeiro discutimos brevemente a teoria e as principais técnicas para a caracterização de pulsos curtos. Analisamos detalhamente uma destas técnicas e propusemos modi cações que possibilitaram o barateamento dos custos da montagem e, além disso, introduziram novas características que sanaram alguns problemas que a montagem original apresentava e que também possibilitaram uma melhor caracterização da própria técnica. Descrevemos cuidadosamente as condições que devem ser satisfeitas pela geometria dos feixes e pelos componentes da montagem para obter uma caracterização correta dos pulsos curtos. Também apresentamos o procedimento de calibração do sistema. Pulsos com diferentes tempos e funções de fase foram caracterizados e os resultados foram validados por testes de con abilidade da informação recuperada. O trabalho seguinte foi o estudo da dinâmica molecular em líquidos puros e em misturas através da técnica efeito Kerr óptico resolvido no tempo, que é uma técnica do tipo bombeio e prova não ressonante, usando um sistema laser Ti:Sa ra com pulsos de 170 fs, centrados em 800 nm. As moléculas estudadas foram o dissulfeto de carbono (CS2), benzeno (C6H6), alilbenzeno (C9H10) e o poliestireno (PS). A teoria necessária para descrever os resultados medidos foi desenvolvida no regime temporal. O modelo de Debye para a relaxação difusiva da anisotropia orientacional descreve a componente de tempos longa, acima de um picosegundo. Partindo do Hamiltoniano de interação, desenvolvemos a teoria da função resposta linear, chegando a uma expressão para a relaxação da polarizabilidade anisotrópica, necessária para descrever os tempos curtos (subpicosegundos). Além disso, a passagem para o regime espectral utilizando somente dados experimentais, ou seja, sem a necessidade de levar em conta modelos especí cos, também foi discutida. Os resultados mostram que os tempos difusivos tanto nos líquidos puros quanto nas misturas seguem a equação de Debye-Stokes-Einstein que prevê um aumento deste tempo para viscosidades maiores. Os tempos curtos são analisados em termos da componente não difusiva da resposta espectral associada à dinâmica molecular. As alterações do espectro foram quanti cadas e a explicação para as mudanças observadas foi dada em termos das con gurações estruturais de interação que levam a uma alteração do potencial intermolecular dentro do qual as moléculas executam movimentos libracionais. Por último, investigamos a questão do modelamento de pulsos de luz incoerente. Para isto trabalhamos com um laser de corante banda larga sem elemento de seleção espectral intracavidade. A técnica utilizada foi o espalhamento forçado de luz ao qual foi acoplado um modelador composto por uma grade, uma lente e uma máscara que alterava a função de fase espectral relativa entre os feixes formadores da grade transiente na amostra de interesse. Realizamos uma análise detalhada desta montagem e obtivemos expressões para ajustar os dados medidos. Os resultados mostram que a função de correlação pode ser alterada de forma especí ca através da escolha de determinadas funções de fase espectrais.

Estabilidade não linear de um feixe de partículas carregadas sujeito a um campo magnético focalizador

Uma análise da estabilidade não-linear de feixes com seção transversal circular considerando perturbações sem simetria axial é executada. Mostra-se que as oscilações simétricas oficialmente circulares de um feixe podem induzir oscilações não-lineares do tipo anti-simétrica(elípticas), com um conseqüente aumento do tamanho do feixe ao logo de uma direção preferencial. O mecanismo da instabilidade e sua relevância às perdas de partícula no feixe são discutidos.

Modelo unidimensional do atrito em escala atômica : um estudo por dinâmica molecular

o trabalho de Krim e Windom [Phys. Rev B 38, 12184(1988)], que revelou a natureza viscosa do atrito em escala atômica, gerou uma intensa atividade tanto teórica quanto experimental. Contudo, questões fundamentais ainda se matém em aberto com opor exemplo a relação entre o coeficiente de atrito viscoso e a topologia do substrato assim como a dependência com a temperatura do substrato. Neste trabalho apresentamos os resultados, obtidos usando dinâmica molecular, para um modelo unidimensional de um sistema adsorvato/substrato. Pesquisamos diferentes relações de comensuração entre o adsorvato e o substrato assim como também a dependência do coeficiente de atrito fonônico, Nph, com a temperatura. Para todas as configurações estudadas obtivemos que o coeficiente de atrito fonônico depende quadráticamente da amplitude de corrugação do substrato, mas tem uma dependência não trivial com a razão de comensuração adsorvato/substrato. O resultado mais impressionante é que para relações de comensuração entre 0.65 e 0.9 o atrito fonônico é muito fraco. Nosso resultados podem explicar as diferencias encontradas na literatura em relação à magnitude do atrito de origem eletrônico e o de origem fonônico.