Vibrações naturais em um sistema de interação viga-água sem o efeito de onda de superfície livre

Neste trabalho é apresentado um estudo sobre o comportamento dinâmico de sistemas flexíveis com interação viga-água. Especificamente, é escolhida uma barragem freqüentemente encontrada em problemas de engenharia de irrigação e represas. O modelo matemático utilizado para descrever tal fenômeno tem como variáveis principais a pressão hidrodinâmica e os deslocamentos horizontais da barragem. Uma vez formuladas as equações governantes de tal sistema, que resultam em uma equação de onda para a pressão e uma equação de Euler-Bernoulli para os deslocamentos, é utilizada a técnica de separação de variáveis para proceder à sua solução. Por simplicidade, assumiu-se que no modelo não havia o efeito de onda da superfície livre. Primeiro, foi resolvida a equação para a pressão, sendo calculadas as suas componentes temporais e espaciais. Depois, foi resolvida a equação de Euler-Bernoulli junto com a condição de interface. Foram simulações para o material específico, e apresentados os resultados conseguidos. No final do trabalho são enunciadas as conclusões pertinentes.

Vibrações naturais em um sistema de interação viga-água incluindo o efeito de onda de superfície

Este trabalho tem como objetivo o estudo modal de um sistema estruturafluído, modelado pela equação de Euler-Bernoulli, para uma viga elástica, sujeita a pressão da água, e por ondas de superfície livre. O sistema acoplado possui condições para o domínio sólido (viga fixa-livre), para o domínio fluído (impermeabilidade e rigidez inferior), com ondas de superfície, e de interfacec fluido-estrutura (condições de continunidade na deflexão, ângulo de rotação, força de corte interno e momento de curvatura). Para deteerminar as vibrações livres e deslocamento no seco e no molhado, utiliza-se o método espectral para eliminar a dependência oscilatória temporal e concentrar-se na determinação dos modos através do estudo de problemas de contornos espaciais. Os modos podem ser calculados com o uso da base clássica de Euler ou da base dinâmica gerada pela resposta impulso. Foram feitas simulações para um material específico, e apresentados os resultados obtidos.

Solução de um problema não linear do tipo reação-difusão na modelagem dispersão de insetos

Neste trabalho tratamos da solução de um problema não linear do tipo tração-difusão, na modelagem de dispersão de insetos. Começamos estabelecendo uma lei de conservação e a partir desta, deduzimos algumas equações importantes para o desenvolvimento do nosso estudo, tais como a equação de convecção, de difusão e simultaneamente convecção e difusão. Se considerarmos uma escala de tempo que possibilite a adição ou retirada de indivíduos no meio, conforme seja considerada reprodução, migração ou morte, podemos acrescentar ao processo difusivo um termo de reação, obtendo então, a equação do tipo reação-difusão. Se o temp de reação for deendendee da densidade populacional e do tipo logístico, obtém-se a equação de Fischer. Dessa equação abordamos alguns aspectos, tais como, determinação dos estados estacionários, análise da estabillidade dos mesmos, representação gráfica no plano de fase e por último investigamos a existência de solução do tipo onda viajante. Abordamos, também, alguns exemplos apresentados na literatura, envolvendo equação da difusão com coeficiente constante e com coeficiente dependente da densidade populacional. Além disso, apresentamos o resultados obtidos com a modelagem em tempo discreto, a partir de um trabalho experimental com besouros marcados para o experimento e depois liberados Banks et al (1985) , em que os autores admitiram uma variação temporal e a partir dos dados obtidos fizeram uma estimativa para o coeficiente de difusão D (t), bem como para o coeficiente de decaimento α(t) do termo de reação linear em u. Construimos curvas que se ajustam a essas alternativas e apresentamos esses coeficientes em versão continua D (t) e α(t), dependentes da variável tempo t. Através de uma abordagem numérica, os modelos foram comparadas da variável tempo t. Através de uma abordagem numérica, os modelos foram comparados para diversos casos, usando diferentes combinação de D constante e D variando no tempo, a constante e a variando no tempo. Além disso, analisamos tambén, o efeito da substituição do coeficiente de difusão D constante por D(t) na equação de Fisher.

Solução da equação de transporte multidimensional em geometria cartesiana e meio infinito usando derivada fracionária

Neste trabalho, foi construída uma forma integral para a solução das equações de transporte em uma, duas e três dimensões, considerando o núcleo de espalhamento de Klein-Nishina, espalhamento isotrópico e o núcleo de espalhamento de Rutherford, respectivamente, seguindo a mesma idéia proposta em trabalhos recentes, nos quais foi construída uma solução para a equação de transporte de nêutrons em geometria cartesiana, usando derivada fracionária. A metodologia consiste em igualar a derivada fracionária do fluxo angular à equação integral, determinar a ordem da derivada fracionária comparando o núcleo da equação integral com o da definição de Riemann-Liouville. Essa formulação foi aplicada ao cálculo de dose absorvida. São apresentadas soluções geradas a partir do emprego do método da derivada fracionária e comparadas a resultados disponíveis na literatura.

Estudo, proposta e previsão de desempenho de um filtro de energia para experimentos com alta resolução de comprimento de onda em fontes pulsantes de nêutrons

Em fontes de nêutrons por spallation para pesquisa de materiais, o comprimento de onda dos nêutrons é geralmente determinado pelos tempos de percurso (TOF) dos nêutrons desde a fonte até o detector. A precisão atingível é limitada pelo fato de o tempo de emissão do pulso característico do sistema alvo/moderador ser diferente de zero (a situação ideal fictícia seria a emissão se ocorresse na forma de um impulso). “Moderadores acoplados” (elementos usados para produzir feixes de alta intensidade com nêutrons de baixa energia) apresentam um decaimento de intensidade em função do tempo muito longo, ao longo de todo o espectro usado nos experimentos. Por este motivo, “moderadores desacoplados”, os quais produzem feixes com intensidade mais reduzida, são freqüentemente usados para instrumentos que requerem alta resolução. Neste trabalho, propusemos e analisamos uma nova técnica de filtragem dinâmica de feixes de nêutrons polarizados de baixa energia para experimentos que utilizam TOF na determinação do comprimento de onda. O dispositivo consiste de um sistema ótico polarizador/analisador e um inversor de spin seletivo em energia, o qual funciona por ressonância espacial do spin. Variando a condição de ressonância em sincronia com a estrutura temporal do pulso de nêutrons (através do controle de campos magnéticos), o filtro pode ajustar a resolução de energia (ou de comprimento de onda) de pulsos de banda larga em tais experimentos, separando os nêutrons com a correta relação “TOF/comprimento de onda” dos demais Um método para o cálculo de desempenho do sistema foi apresentado em um Trabalho Individual (TI) (PARIZZI et al., 2002 - i), do qual se fará aqui uma breve revisão além de amplo uso na otimização dos parâmetros a serem ajustados para o projeto do filtro. Os resultados finais mostram que ganhos consideráveis em resolução podem ser obtidos com a aplicação desta técnica em experimentos de reflectometria por tempo de percurso, sem que para tal seja necessário comprometer a intensidade do feixe usado pelo mesmo instrumento quando operado em um modo de baixa resolução, dando ao usuário do instrumento a opção de escolher a relação ótima entre intensidade e resolução para seu experimento. Como parte da conclusão desta dissertação, é apresentada uma proposta de parâmetros para a construção deste tipo de filtro e previsão de desempenho da configuração proposta, baseada no software de modelamento desenvolvido.

Solução da equação de transporte de fótons para uma placa plana heterogênea, modelo de multigrupo com núcleo de espalhamento de Klein-Nishina

Neste trabalho o método LTSN é utilizado para resolver a equação de transporte de fótons para uma placa plana heterogênea, modelo de multigrupo, com núcleo de espalhamento de Klein-Nishina, obtendo-se o fluxo de fótons em valores discretos de energia. O fluxo de fótons, juntamente com os parâmetros da placa foram usados para o cálculo da taxa de dose absorvida e do fator de buildup. O método LTSN consiste na aplicação da transformada de Laplace num conjunto de equações de ordenadas discretas, fornece uma solução analítica do sistema de equações lineares algébricas e a construção dos fluxos angulares pela técnica de expansão de Heaviside. Essa formulação foi aplicada ao cálculo de dose absorvida e ao fator de Buildup, considerando cinco valores de energia. Resultados numéricos são apresentados.

Estudo e solução da equação de transporte de nêutrons bidimensionais pelo método LTSn para elevadas ordens de quadraturas angulares : LTSn2D - Diag e LTSn2D - DiagExp

O principal objetivo dessa tese consiste em determinar uma solução numéricada equação bidimensional do transporte de nêutrons para elevadas ordens de quadratura angular. Diagonalizando a matriz de transporte LTSN bidimensional , construímos dois algoritmos que se diferenciam pela forma de representar os termos de fuga transversal, que surgem nas equações LTSN integradas transversalmente. Esses termos no método LTSN2D − Diag são expressos como combinação linear dos autovetores multiplicados por exponenciais dos respectivos autovalores. No método LTSN2D − DiagExp os termos de fuga transversal são representados por uma função exponencial com constante de decaimento heuristicamente identificada com parâmetros materiais característicos do meio. A análise epectral desenvolvida permite realizar a diagonalização. Um estudo sobre o condicionamento é feito e também associamos um número de condicionamento ao termo de fuga transversal. Definimos os erros no fluxo aproximado e na fórmula da quadratura, e estabelecemos uma relação entre eles. A convergência ocorre com condições de fronteira e quadratura angular adequadas. Apresentamos os resultados numéricos gerados pelos novos métodos LTSN2D − Diag e LTSN2D − DiagExp para elevadas ordens de quadratura angular para um problema ilustrativo e comparamos com resultados disponíveis na literatura.

A equação de Poisson-Boltzmann em regiões com fronteira irregular

Propomos uma idealização da situação em que uma macromolécula é ionizada em um solvente. Neste modelo a área da superfície da molécula é suposta ser grande com respeito a seu diâmetro. A molécula é considerada como um dielétrico com uma distribuição de cargas em sua superfície. Utilizando as condições de transmissão, a distribuição de Boltzmann no solvente e resultados recentes sobre espaços de Sobolev no contexto de espaços métricos, bem como de integração sobre superfícies irregulares, o problema é formulado em forma variacional. Resultados clássicos do cálculo de variações permitem a resolução analítica do problema.

Solução da equação de condução de calor bidimensional, em meios multicompostos, pelos métodos nodal, com parâmetros concentrados, e a técnica da transformada de Laplace

Neste trabalho, desenvolvemos uma metodologia semi-analítica para solução de problemas de condução de calor bidimensional, não-estacionária em meios multicompostos. Esta metodologia combina os métodos nodal, com parâmetros concentrados, e a técnica da transformada de Laplace. Inicialmente, aplicamos o método nodal. Nele, a equação diferencial parcial que descreve o problema é integrada, transversalmente, em relação a uma das variáveis espaciais. Em seguida, é utilizado o método de parâmetros concentrados, onde a distribuição de temperatura nos contornos superior e inferior é substituída pelo seu valor médio. Os problemas diferenciais unidimensionais resultantes são então resolvidos com o uso da técnica da transformada de Laplace, cuja inversão é avaliada numericamente. O método proposto é usado na solução do problema de condução de calor, em paredes de edificações. A implementação computacional é feita, utilizando-se a linguagem FORTRAN e os resultados numéricos obtidos são comparados com os disponíveis na literatura.

O problema de Dirichlet para a equação das superfícies mínimas em domínios não necessariamente convexos

Estudamos o problema de Dirichlet para a equação das superfícies mínimas em domínios limitados do plano. Provamos a existência e unicidade de gráficos mínimos sobre domínios limitados e não necessariamente convexos, com valores no bordo satisfazendo uma condição que denominamos condição da declividade limitada generalizada a qual, usando cilindros no lugar de planos, generaliza a condição clássica da declividade limitada. Com este resultado, dado um domínio limitado e suave qualquer do plano, conseguimos obter cotas explícitas para a norma C2 de dados no bordo deste domínio que garantem a existência de solução ao correspondente problema de Dirichlet.

Formulação semi-analítica para a equação transformada resultante da aplicação da GITT em problemas difusivos-advectivos

Neste trabalho se propõe um avanço para a Técnica Transformada Integral Generalizada, GITT. O problema transformado, usualmente resolvido por subrotinas numéricas, é aqui abordado analiticamente fazendo-se uso da Transformada de Laplace. Para exemplificar o uso associado destas duas transformadas integrais, resolvem-se dois problemas. Um de concentração de poluentes na atmosfera e outro de convecção forçada com escoamento laminar, entre placas planas paralelas, com desenvolvimento simultâneo dos perfis térmico e hidrodinâmico. O primeiro é difusivo, transiente e com coeficientes variáveis. Sua solução é obtida de forma totalmente analítica. Além de mostrar o uso da técnica, este exemplo apesar de ter coeficientes variáveis, é resolvido com o auxílio de um problema de autovalores associado com coeficientes constantes. No segundo, obtém-se a solução da Equação da Energia analiticamente. Já a Equação da Conservação do Momentum é linearizada e resolvida de forma iterativa. A solução de cada iteração é obtida analiticamente.

Solução analítica da equação de difusão-advenção pelo método GILTT aplicada à dispersão de poluentes atmosféricos

O objetivo deste trabalho é obter uma nova solução analítica para a equação de advecção-difusão. Para tanto, considera-se um problema bidimensional difusivo-advectivo estacionário com coeficiente de difusão turbulenta vertical variável que modela a dispersão de poluentes na atmosfera. São utilizados três coeficientes difusivos válidos na camada limite convectiva e que dependem da altura, da distância da fonte e do perfil de velocidade. A abordagem utilizada para a resolução do problema é a técnica da Transformada Integral Generalizada, na qual a equação transformada do problema difusivo-advectivo é resolvida pela técnica da Transformada de Laplace com inversão analítica. Nenhuma aproximação é feita durante a derivação da solução, sendo assim, esta é exata exceto pelo erro de truncamento. O modelo ´e avaliado em condições moderadamente instáveis usando o experimento de Copenhagen. Apresentam-se os resultados numéricos e estatísticos, comparando os resultados obtidos com dados experimentais e com os resultados da literatura. O modelo proposto mostrou-se satisfatório em relação aos dados dos experimentos difusivos considerados.

Modelos cinéticos da equação linearizada de Boltzmann e um problema de transferência de calor em microescala

Neste trabalho, um problema de transferência de calor da dinâmica de gases rarefeitos, causado pela diferença de temperaturas nas superfícies de um canal, é abordado. O problema é formulado através dos modelos cinéticos BGK, S e Gross-Jackson da equação linearizada de Boltzmann e resolvido, de forma unificada, pelo método analítico de ordenadas discretas (método ADO). Resultados numéricos para as perturbações de densidade e temperatura e também para o fluxo de calor são apresentados e comparados, mostrando que não se pode dizer que algum dos três modelos seja uma melhor aproximação da solução aos resultados da equação linearizada de Boltzmann.

Modelagem estratigráfica de clinoformas deposicionais : construção e aplicação de um modelo computacional baseado em mecânica estatística

Sob a premissa de que a geometria do perfil deposicional das clinoformas pode conter informações sobre a composição textural dos sedimentos que as compõem, bem como sobre a energia da coluna d’água, foi desenvolvido um modelo computacional para a análise teórica dos fatores controladores do tipo de curvatura dos foresets das clinoformas, e, por extensão dos taludes submarinos. Um modelo análogo de suspensão dinâmica foi implementado experimentalmente com um programa em código Matlab, denominado MAXWELL, e é classificado como um modelo estratigráfico, bidimensional, analítico, determinístico, dinâmico, parcialmente dedutivo e parcialmente baseado em regras. Contém um modelo de escoamento de fluido em linhas de corrente, e trata indiretamente a tensão de cisalhamento no domínio de um sistema fechado bifásico líquido-vapor análogo, a partir a proposta de analogias entre os processos de evaporação e de suspensão de sedimentos. É uma abordagem baseada na competência de transporte do fluxo d’água, pois considera temperatura e velocidade de fluxo combinado onda-corrente como variáveis controladoras análogas da intensidade de transferência de energia. Os processos deposicionais marinhos são reduzidos a um fenômeno da superfície deposicional, que é considerada análoga a uma interface líquidovapor dentro de um sistema fechado. A equação de distribuição estatística de velocidades moleculares de Maxwell é usada para calcular a proporção de moléculas na fase líquida, considerada análoga à proporção de sedimentos depositados, bem como a proporção na fase vapor, tomada como análoga à proporção de sedimentos em suspensão. O estado energético do sistema é parametrizado por três potenciais: energia interna, como função do tamanho de grão (areia, silte ou argila); energia térmica do meio, como função da energia hidrodinâmica, e energia gravitacional, como função do declive topográfico As simulações indicam que os principais fatores controladores do perfil deposicional das clinoformas, e por extensão, dos taludes submarinos em geral, são a energia hidrodinâmica da coluna d’água e a granulometria (ou coesão) do sedimento, que podem ser consideradas dois parâmetros comutáveis, isto é, grãos grossos ou coesos podem produzir sobre a geometria das clinoformas o mesmo efeito da baixa energia hidrodinâmica, e vice-versa. Com base no fator de decaimento da energia hidrodinâmica com o aumento da batimetria, foram definidas três ordens de grandeza da intensidade da energia da coluna d’água: baixa energia (10–1), alta x energia (10-2) e muito alta energia (10-3). Com base nesse critério, foram caracterizados quatro tipos de perfis deposicionais de clinoformas: dois tipos sigmoidais e dois tipos exponenciais. Os tipos sigmoidais podem ser de alta energia ou de muito alta energia, e distinguem-se pela granulometria do sedimento e pela distribuição de declividades em relação à dimensão na direção horizontal. Os perfis de tipo exponencial podem ser tanto de baixa energia, quanto de alta energia. A subida relativa do nível do mar afeta o tipo geométrico do perfil deposicional apenas de forma indireta, a partir da retenção da fração grossa do influxo sedimentar na plataforma ou topset. Os principais fatores controladores do ângulo de declividade do perfil das clinoformas são a granulometria do sedimento, a energia hidrodinâmica, a profundidade d’água da bacia e o desvio lateral da corrente de fundo. A dedução da litofácies a partir da geometria das clinoformas parece promissora apenas para os tipos teóricos sigmoidais, que são distintos na forma e no conteúdo sedimentar.

Solução LTSn da equação adjunta de transporte de nêutrons com fonte arbitrária para elevada ordem de quadratura numa placa homogênea

O objetivo deste trabalho consiste em estender o método LTSN à solução do problema adjunto de transporte de nêutrons. A solução adjunta é interpretada fisicamente como uma função importância que designa a capacidade de contribuição de cada cela do espaço de fase para um funcional resposta. A derivação desta interpretação, através do princípio variacional, está sucintamente apresentada. Surgida da necessidade de generalização da fonte adjunta, também propõe-se uma nova formulação LTSN capaz de resolver problemas de transporte, tanto direto quanto adjunto, com fonte arbitrária, para elevada ordem de quadratura em geometria de placa. Esta nova formulção inspira-se na propriedade de invariância de projeção dos meios isotrópicos mas também é válida para os meios anisotrópicos. Todos os resultados apresentados pelas simulações numéricas de problemas adjuntos são calculados pela nova formulação LTSN e são comparados ou com a definição de função importância ou pelas relações de reciprocidade ou pelo código ANISN.