Uma implementação geral de elementos de contorno e sua aplicação a problemas de flexão de placas

o presente trabalho aborda a aplicação do método dos elementos de contorno (MEC) para solução de problemas de flexão linear e geometricamente não-linear de placas semiespessas. Os modelos de placa empregados consideraram a influência do cisalhamento através de teorias de primeira ordem, especificamente as de Mindlin e Reissner. Uma formulação integral unificada dos modelos de placa utilizados é desenvolvida para o operador de Navier do problema, onde foram mantidos alguns termos de ordem superior no tensor deformação de Green. A formulação integral do problema de membrana acoplado ao de flexão é igualment desenvolvida, levando a um sistema de equações integrais não-lineares que descreve completamente problemas de placas que envolvem grandes deslocamentos. Estas equações podem ser particularizadas para problemas de flexão linear e estabilidade elástica. Tendo em vista a necessidade de se considerar derivadas dos deslocamentos translacionais, as equações integrais correspondentes ao gradiente dos deslocamentos também foram deduzidas, caracterizando uma formulação hipersingular. o método empregado para solução numérica do sistema de equações integrais foi o método direto dos elementos de contorno. Um tratamento das integrais fortemente singulares presentes nas equações foi realizado, baseado em expansões assint6ticas dos núcleos. Deste procedimento resulta uma abordagem regularizada que emprega apenas quadraturas padrão de Gauss-Legendre.

Métodos de simulação para ciclos de Rankine

O presente trabalho é dedicado ao estudo de métodos de simulação para ciclos de Rankine. O trabalho é iniciado com a modelagem de um ciclo de Rankine simples e segue evoluindo para configurações mais complexas tal como o ciclo de Rankine com reaquecimento e regeneração. São adotadas as considerações mais convencionais da prática de projeto de centrais termelétricas cujos sistema térmicos baseiam-se no ciclo de Rankine, incluindo-se queda de pressão em tubulações do circuito além de outras perdas. Em seguida, são estabelecidas as expressões matemáticas que possibilitam a determinação das propriedades termodinâmicas da água em seus mais diversos estados ao longo do ciclo. Por último, são desenvolvidos métodos de simulação, chamados neste trabalho de Substituição Sucessiva e Bloco Único, que caracterizam-se pela resolução simultânea do conjunto de equações algébricas dos ciclos elaborados. As simulações são efetuadas através de programas escritos na linguagem Fortran. Os métodos de simulação são aplicados para a obtenção dos resultados considerados mais importantes na análise de sistemas térmicos de potência, tais como rendimento térmico do ciclo, título na saída da turbina, vazões mássicas pelo sistema, potência nas bombas e calor trocado no gerador de vapor e no condensador Na maioria das simulações, estes resultados apresentam-se como funções da: (1) potência elétrica requerida, eficiência isentrópica e pressões na turbina; (2) eficiência térmica, pressão e temperatura no gerador de vapor; (3) pressão e grau de subresfriamento do líquido saturado no condensador e (4) eficiência isentrópica das bombas. São obtidos os mesmos resultados para os métodos de simulação utilizados. O método da Substituição Sucessiva apresentou menor tempo computacional, principalmente para configurações de ciclo mais complexas. Uma aplicação alternativa do método de Bloco Único demonstrou ser inconveniente para ciclos de configurações mais complexas devido ao elevado tempo computacional, quando todas as equações de cálculo das propriedades termodinâmicas são incluídas no sistema de equações a ser resolvido. Melhores rendimentos térmicos e título na saída da turbina foram obtidos para configurações de ciclo de Rankine com reaquecimento e regeneração.

Filtragem de imagens com preservação das bordas usando a Transformada Wavelet

A filtragem de imagens visando a redução do ruído é uma tarefa muito importante em processamento de imagens, e encontra diversas aplicações. Para que a filtração seja eficiente, ela deve atenuar apenas o ruído na imagem, sem afetar estruturas importantes, como as bordas. Há na literatura uma grande variedade de técnicas propostas para filçtragem de imagens com preservação de bordas, com as mais variadas abordagens, deentrte as quais podem ser citadas a convolução com máscaras, modelos probabilísticos, redes neurais, minimização de funcionais e equações diferenciais parciais. A transformada wavelet é uma ferramenta matemática que permite a decomposição de sinais e imagens em múltiplas resoluções. Essa decomposição é chamada de representação em wavelets, e pode ser calculada atrravés de um algorítmo piramidal baseado em convoluções com filtros passa-bandas e passa-baixas. Com essa transformada, as bordas podem ser calculadas em múltiplas resoluções. Além disso, como filtros passa-baixas são utilizados na decomposição, a atenuação do ruído é um processo intrínseco à transformada. Várias técnicas baseadas na transformada wavelet têm sido propostas nos últimos anos, com resultados promissores. Essas técnicas exploram várias características da transformada wavelet, tais como a magnitude de coeficientes e sua evolução ao longo das escalas. Neste trabalho, essas características da transformada wavelet são exploradas para a obtenção de novas técnicas de filtragem com preservação das bordas.

Análise numérica de escoamentos turbulentos tridimensionais empregando o método de elementos finitos e simulação de grandes escalas

O presente trabalho apresenta o estudo e implementação de um algoritmo numérico para análise de escoamentos turbulentos, tridimensionais, transientes, incompressíveis e isotérmicos, através da Simulação de Grande Escalas, empregando o Método de Elementos Finitos. A modelagem matemática do problema baseia-se nas equações de conservação de massa e quantidade de movimento de um fluido quase-incompressível. Adota-se um esquema de Taylor-Galerkin, com integração reduzida e fórmulas analíticas das funções de interpolação, para o elemento hexaédrico de oito nós, com funções lineares para as componentes de velocidade e constante no elemento para a pressão. Para abordar o problema da turbulência, emprega-se a Simulação de Grandes Escalas, com modelo para escalas inferiores à resolução da malha. Foram implementados o modelo clássico de Smagorinsky e o modelo dinâmico de viscosidade turbulenta, inicialmente proposto por Germano et al, 1991. Uma nova metodologia, denominada filtragem por elementos finitos independentes, é proposta e empregada, para o processo de segunda filtragem do modelo dinâmico. O esquema, que utiliza elementos finitos independentes envolvendo cada nó da malha original, apresentou bons resultados com um baixo custo computacional adicional. São apresentados resultados para problemas clássicos, que demonstram a validade do sistema desenvolvido. A aplicabilidade do esquema utilizado, para análise de escoamentos caracterizados por elevados números de Reynolds, é discutida no capítulo final. São apresentadas sugestões para aprimorar o esquema, visando superar as dificuldades encontradas com respeito ao tempo total de processamento, para análise de escoamentos tridimensionais, turbulentos e transientes .

Modelagem física e matemática dos processos de extração de óleo de soja em fluxos contracorrente cruzados

Este trabalho visa desenvolver um modelo físico e matemático geral para os processos de extração sólido-líquido em fluxos contracorrente cruzados (CCC) que são utilizados na indústria de alimentos. Levam-se em consideração os processos principais (o transporte de massa entre as fases, difusão e convecção) envolvidos por todo o campo de extração, com uma abordagem bidimensional evolutiva, incluindo as zonas de carregamento, drenagem e as bandejas acumuladoras. O modelo matemático é formado por equações diferenciais parciais que determinam a alteração das concentrações nas fases poro e “bulk” em todo o campo de extração e equações diferenciais ordinárias (que refletem as evoluções das concentrações médias nas bandejas). As condições de contorno estabelecem as ligações entre os fluxos CCC da micela e matéria-prima e consideram, também, a influência das zonas de drenagem e carregamento. O algoritmo de resolução utiliza o método de linhas que transforma as equações diferenciais parciais em equações diferenciais ordinárias, que são resolvidas pelo método de Runge-Kutta. Na etapa de validação do modelo foram estabelecidos os parâmetros da malha e o passo de integração, a verificação do código com a lei de conservação da espécie e um único estado estacionário. Também foram realizadas a comparação com os dados experimentais coletados no extrator real e com o método de estágios ideais, a análise da influência de propriedades da matéria-prima nas características principais do modelo, e estabelecidos os dados iniciais do regime básico (regime de operação) Foram realizadas pesquisas numéricas para determinar: os regimes estacionário e transiente, a variação da constante de equilíbrio entre as fases, a variação do número de seções, a alteração da vazão de matéria-prima nas características de um extrator industrial e, também foram realizadas as simulações comparativas para diferentes tipos de matéria-prima (flocos laminados e flocos expandidos) usados amplamente na indústria. Além dessas pesquisas, o modelo também permite simular diferentes tipos de solventes. O estudo da capacidade de produção do extrator revelou que é necessário ter cuidado com o aumento da vazão da matéria-prima, pois um pequeno aumento desta pode causar grandes perdas de óleo tornando alto o custo da produção. Mesmo que ainda seja necessário abastecer o modelo com mais dados experimentais, principalmente da matéria-prima, os resultados obtidos estão em concordância com os fenômenos físico-químicos envolvidos no processo, com a lei de conservação de espécies químicas e com os resultados experimentais.

Modos em vigas com secção transversal de variação linear

O objetivo principal deste trabalho é a obtenção dos modos e as freqüências naturais de vigas de variação linear e em forma de cunha, com condições de contorno clássicas e não-clássicas, descritas pelo modelo estrutural de Euler-Bernoulli. A forma dos modos foi determinado com o uso das funções cilíndricas. No caso forçado se considera uma força harmônica e se resolve o problema pelo método espectral, utuilizando o software simbólico Maple V5. Realiza-se uma análise comparativa dos resultados obtidos com os resultados existentes na literatura para vigas uniformes.

Simulação numérica de esteiras em transição utilizando o método dos contornos virtuais

Simulações Numéricas são executadas em um código numérico de alta precisão resolvendo as equações de Navier-Stokes e da continuidade para regimes de escoamento incompressíveis num contexto da turbulência bidimensional. Este código utiliza um esquema compacto de diferenças finitas de sexta ordem na aproximação das derivadas espaciais. As derivadas temporais são calculadas usando o esquema de Runge-Kuta de terceeira ordem com baixo armazenamento. Tal código numérico fornece uma representação melhorada para uma grande faixa de escalas de comprimento e de tempo. As técnicas dos contornos imersos acopladas ao método dos contornos virtuais permitem modelar escoamentos não-estacionários sobre geometrrias complexas, usando simplesmente uma malha Cartesiana uniforme. Por meio de procedimentos de aproximação/interpolação, as técnicas dos contornos imersos (aproximação Gaussiana, interpolação bilinear e redistribuição Gaussiana), permitem a representação do corpo sólido no interior do campo de escoamento, com a superfície não coincidindo com a malha computacional. O método dos contornos virtuais, proposto originalmente por Peskin, consiste, basicamente, na imposição na superfície e/ou no interior do corpo, de um termo de força temporal acrescentando às equações do momento. A aplicação deste campo de força local leva o fluido ao repouso na superfície do corpo, permitindo obter as condições de contorno de não-deslizamento e de não penetração de fluido na parede. A análise das oscilações induzidas no escoamento-contorno pelo processo de desprendimento de vórtices na esteira do cilindro circular e de geometria retangulares na incidência, para números de Reybolds variando de 40 a 400, confirma a eficiência computacional e a aplicabilidade das técncias implementadas.

Um assistente de feedback para o serviço de filtragem do software direto

Este trabalho descreve a especificação e implementação do protótipo Assistente de Feedback que ajuda os usuários a ajustarem os parâmetros do serviço de filtragem de mensagens vindas do correio eletrônico de sistemas como o Direto. O Assistente de Feedback é instalado no computador do usuário do Direto para monitorar suas preferências representadas pelas ações aplicadas nas mensagens do correio eletrônico. O trabalho apresenta, ainda, uma revisão bibliográfica sobre os conceitos gerais de probabilidades, redes Bayesianas e classificadores. Procura-se descrever as características gerais dos classificadores, em especial o Naive Bayes, sua lógica e seu desempenho comparado a outros classificadores. São abordados, também, conceitos relacionados ao modelo de perfil de usuário e o ambiente Direto. O Naive Bayes torna-se atraente para ser utilizado no Assistente de Feedback por apresentar bom desempenho sobre os demais classificadores e por ser eficiente na predição, quando os atributos são independentes entre si. O Assistente de Feedback utiliza um classificador Naive Bayes para predizer as preferências por intermédio das ações do usuário. Utiliza, também, pesos que representarão a satisfação do usuário para os termos extraídos do corpo da mensagem. Esses pesos são associados às ações do usuário para estimar os termos mais interessantes e menos interessantes, pelo valor de suas médias finais. Quando o usuário desejar alterar os filtros de mensagens do Direto, ele solicita ao Assistente de Feedback sugestões para possíveis exclusões dos termos menos interessantes e as possíveis inclusões dos termos mais interessantes. O protótipo é testado utilizando dois métodos de avaliação para medir o grau de precisão e o desempenho do Assistente de Feedback. Os resultados obtidos na avaliação de precisão apresentam valores satisfatórios, considerando o uso de cinco classes pelo classificador do Assistente de Feedback. Os resultados dos testes de desempenho permitem observar que, se forem utilizadas máquinas com configurações mais atualizadas, os usuários conseguirão receber sugestões com tempo de respostas mais toleráveis.

Serviços ao cliente e marketing de relacionamento no setor hoteleiro de Porto Alegre

Nos últimos vinte anos, diversos fatores vêm moldando o panorama econômico e mercadológico, entre os quais destacam-se a globalização, tecnologia e competição cada vez mais acirrada. Neste cenário, a oferta de serviços ao cliente e as relações comerciais com foco no longo prazo, baseados em benefícios mútuos, figuram como importantes fontes de vantagem competitiva sustentável. Seguindo as tendências de pesquisas contemporâneas, grandes atenções por parte da academia têm-se voltado para o desenvolvimento de modelos teóricos que visam explicar o comportamento das dimensões determinantes para o estabelecimento e manutenção de trocas relacionais. Ao mesmo tempo, os serviços ao cliente vêm sendo estudados como uma atividade de relevância estratégica para potenciais de diferenciação. Este estudo objetiva examinar as correlações existentes entre a oferta de serviços ao cliente e os atributos de Marketing de Relacionamento nas transações entre redes hoteleiras de Porto Alegre e seus clientes organizacionais. A mensuração do Marketing de Relacionamento é feita através de constructos como dependência do cliente, nível de comparação com clientes alternativos, investimentos no relacionamento, troca de informações, confiança e comprometimento. O nível de serviços ao cliente é obtido através de indicadores gerados em entrevistas em profundidade com informantes-chave dos hotéis e alguns de seus clientes organizacionais. Os resultados obtidos através da utilização de modelagem de equações estruturais sugerem bons índices de ajustamento para a estrutura teórica proposta sobre o comportamento dos atributos de Marketing de Relacionamento e serviços ao cliente, apontando para a confirmação inicial das hipóteses de pesquisa.

Uma metodologia unificada no domínio tempo para sistemas concentrados, discretos e distribuídos

A resposta impulso é utilizada como ferramenta padrão no estudo direto de sistemas concentrados, discretos e distribuídos de ordem arbitrária. Esta abordagem leva ao desenvolvimento de uma plataforma unificada para a obtenção de respostas dinâmicas. Em particular, as respostas forçadas dos sistemas são decompostas na soma de uma resposta permanente e de uma resposta livre induzida pelos valores iniciais da resposta permanente. A teoria desenvolve-se de maneira geral e direta para sistemas de n-ésima ordem, introduzindo-se a base dinâmica gerada pela resposta impulso na forma padrão e normalizada, sem utilizar-se a formulação de estado, através da qual reduz-se um sistema de ordem superior para um sistema de primeira ordem. Considerou-se sistemas de primeira ordem a fim de acompanhar-se os muitos resultados apresentados na literatura através da formulação de espaço de estado. Os métodos para o cálculo da resposta impulso foram classificados em espectrais, não espectrais e numéricos. A ênfase é dada aos métodos não espectrais, pois a resposta impulso admite uma fórmula fechada que requer o uso de três equações características do tipo algébrica, diferencial e em diferenças Realizou-se simulações numéricas onde foram apresentados modelos vibratórios clássicos e não clássicos. Os sistemas considerados foram sistemas do tipo concentrado, discreto e distribuído. Os resultados da decomposição da resposta dinâmica de sistemas concentrados diante de cargas harmônicas e não harmônicas foram apresentados em detalhe. A decomposição para o caso discreto foi desenvolvida utilizando-se os esquemas de integração numérica de Adams-Basforth, Strömer e Numerov. Para sistemas distribuídos, foi considerado o modelo de Euler-Bernoulli com força axial, sujeito a entradas oscilatórias com amplitude triangular, pulso e harmônica. As soluções permanentes foram calculadas com o uso da função de Green espacial. A resposta impulso foi aproximada com o uso do método espectral.

Estudo da variabilidade genética e ecológica de comunidades de Drosophila em regiões de Mata Atlântica de ilhas do continente de Santa Catarina

Foi realizado um estudo da dinâmica de assembléias de Drosofilídeos em oito amostras insulares e continentais de Santa Catarina através de dados coletados em várias visitas ao longo de dois anos. Dentre os resultados obtidos está a estimativa do grau de diversidade destas assembléias. Nossas coletas mostraram que a predição de qual espécie será dominante, num determinado período amostrado, é razoavelmente possível. A análise dos índices de diversidade nos indica que o Morro da Lagoa é o ponto de menor diversidade específica, seguido de Ratones Grande. Contudo, os dois pontos têm um grande número de espécies diferentes, sendo a sua diversidade baixa em função da alta dominância do subgrupo willistoni neles encontrada. Analisando o componente S (número de espécies) da diversidade, nas ilhas pontos de coleta, percebe-se que a ilha maior (Ilha de Santa Catarina) tem realmente um maior número de espécies coletadas – 46 no ponto A além de 10 espécies diferentes coletadas no ponto D (56 espécies no total) – do que as ilhas menores (42, 44, 40 e 50), o que corrobora a teoria da biogeografia de ilhas. No continente, a curva espécie/área se comportou da mesma forma que nas ilhas se considerarmos a diversidade como um todo. Realmente estes pontos foram uns dos que apresentaram maior diversidade, principalmente o ponto F, com um H’ de 2,22, que se manteve com Mata Atlântica Primária até o final do período de coletas, sendo portanto o ponto mais preservado de todos utilizados e, teoricamente, o que apresentava maior diversidade de nichos ecológicos para serem ocupados Outros Drosofilídeos como Zaprionus, Zygotricha, Gitona, Cladochaeta, Diathoneura, Micodrosophila, Leucophenga e Amiota foram coletados. Embora nosso interesse preliminar fosse apenas o gênero Drosophila, a inclusão destes outros gêneros em nosso estudo visou uma maior compreensão das possíveis associações que podem ocorrer entre eles e espécies de Drosophila. Os dois primeiros gêneros foram mais freqüentes nas nossas coletas. Zaprionus indianus foi considerada uma espécie invasora, pois surgiu com freqüências baixíssimas que aumentaram gradualmente nas coletas subsequentes, superando em freqüência as espécies nativas. Isto confirma o caráter generalista e polifágico deste tipo de espécie. Neste trabalho, é relatado o primeiro registro do Gênero Zaprionus (Diptera, Drosophilidae) para o Estado de Santa Catarina, na região litorânea central que inclui as Ilhas de Santa Catarina, Arvoredo, Ratones Grande, Ratones Pequeno e Campeche. Drosophila roerhae, D. unipunctata, D. schineri, D. bifilum, D. fuscolineata, D. meridionalis, D. neosaltans, D. bocainoides e D. platitarsus foram pela primeira vez registradas para a região Sul do Brasil, aumentando, portanto o limite meridional de suas distribuições. Como um ponto de partida para estudar o polimorfismo para inversões cromossômicas em D. neocardini, foi construído um fotomapa de referência dos cromossomos politênicos de glândulas salivares de larvas de terceiro estágio. Pelo menos 258 indivíduos (aproximadamente três núcleos por glândula) de sete diferentes localidades (Sertão do Peri, Ilha do Arvoredo, Serra do Tabuleiro, Ilha de Ratones Grande, Ilha de Ratones Pequeno, Morro da Lagoa da Conceição e Ilha do Campeche, todos no Estado de Santa Catarina) foram analisados e fotomicrografias foram obtidas, até se chegar a um consenso sobre a identidade dos elementos cromossômicos. Uma nova inversão no braço cromossômico IIIL foi registrada e denominada de IIILA. A variabilidade cromossômica encontrada nas espécies de Drosophila do grupo cardini em todas as localidades também foi pesquisada, e foi comparada visando contribuir para uma melhor compreensão da evolução destas comunidades. Analisando o polimorfismo cromossômico de D. polymorpha encontramos nove inversões diferentes pela primeira vez descritas. Uma das inversões novas foi encontrada no cromossomo X, duas outras foram encontradas para o braço IIL; quatro foram catalogadas para o braço cromossômico IIIR e duas inversões novas foram achadas no braço cromossômico IIIL. Com relação ao polimorfismo, em D. neocardini foi encontrada apenas uma nova inversão no braço IIIL e para D. cardinoides uma nova inversão no braço IIIL. O estudo discute as implicações ecológicas e evolutivas deste tipo de polimorfismo, para um maior entendimento da evolução deste grupo de espécies.

Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Relação entre orientação para serviços ao cliente e performance das organizações

O estudo do conceito de orientação para o mercado e suas implicações têm despertado interesse de acadêmicos e práticos nos últimos anos. Esta pesquisa trata da utilização do conceito aplicado aos serviços ao cliente, traduzido na orientação para serviço ao cliente. Os modelos de orientação para serviços ao cliente, propostos por alguns teóricos, sugerem que, com a utilização desses serviços como um componente estratégico de diferenciação das empresas, seriam atingidos níveis superiores de performance organizacional. Os serviços ao cliente são o principal componente desses modelos, uma vez que apresentam algumas características únicas, como a sua aplicabilidade nos constructos de orientação para o mercado e a possibilidade de ser mensurado e comparado aos concorrentes. Esta pesquisa propôs-se a verificar a aplicação de um modelo de orientação para serviços ao cliente e relacioná-lo com a performance das organizações. Para tanto, foram pesquisadas as empresas da indústria eletroeletrônica de todo país, analisando-se questões como os serviços prestados, o atendimento de pedidos, os produtos oferecidos, entre outros. Esses indicadores foram comparados com variáveis de performance como: participação de mercado, receita de vendas, lucratividade, etc. Através da técnica de modelagem de equações estruturais, verificou-se que existe uma relação causal entre a orientação para os serviços ao cliente e a performance das empresas, demonstrando que a adoção de tal orientação traduzse em um melhor desempenho das empresas.

Paralelizações de métodos numéricos em clusters empregando as bibliotecas MPICH, DECK e Pthread

Este trabalho tem como objetivo desenvolver e empregar técnicas e estruturas de dados agrupadas visando paralelizar os métodos do subespaço de Krylov, fazendo-se uso de diversas ferramentas e abordagens. A partir dos resultados é feita uma análise comparativa de desemvpenho destas ferramentas e abordagens. As paralelizações aqui desenvolvidas foram projetadas para serem executadas em um arquitetura formada por um agregado de máquinas indepentes e multiprocessadas (Cluster), ou seja , são considerados o paralelismo e intra-nodos. Para auxiliar a programação paralela em clusters foram, e estão sendo, desenvolvidas diferentes ferramentas (bibliotecas) que visam a exploração dos dois níveis de paralelismo existentes neste tipo de arquitetura. Neste trabalho emprega-se diferentes bibliotecas de troca de mensagens e de criação de threads para a exploração do paralelismo inter-nodos e intra-nodos. As bibliotecas adotadas são o DECK e o MPICH e a Pthread. Um dos itens a serem analisados nestes trabalho é acomparação do desempenho obtido com essas bibliotecas.O outro item é a análise da influência no desemepnho quando quando tulizadas múltiplas threads no paralelismo em clusters multiprocessados. Os métodos paralelizados nesse trabalho são o Gradiente Conjugação (GC) e o Resíduo Mínmo Generalizado (GMRES), quepodem ser adotados, respectivamente, para solução de sistemas de equações lineares sintéticos positivos e definidos e não simétricas. Tais sistemas surgem da discretização, por exemplo, dos modelos da hidrodinâmica e do transporte de massa que estão sendo desenvolvidos no GMCPAD. A utilização desses métodos é justificada pelo fato de serem métodos iterativos, o que os torna adequados à solução de sistemas de equações esparsas e de grande porte. Na solução desses sistemas através desses métodos iterativos paralelizados faz-se necessário o particionamento do domínio do problema, o qual deve ser feito visando um bom balanceamento de carga e minimização das fronteiras entre os sub-domínios. A estrutura de dados desenvolvida para os métodos paralelizados nesse trabalho permite que eles sejam adotados para solução de sistemas de equações gerados a partir de qualquer tipo de particionamento, pois o formato de armazenamento de dados adotado supre qualquer tipo de dependência de dados. Além disso, nesse trabalho são adotadas duas estratégias de ordenação para as comunicações, estratégias essas que podem ser importantes quando se considera a portabilidade das paralelizações para máquinas interligadas por redes de interconexão com buffer de tamanho insuficiente para evitar a ocorrência de dealock. Os resultados obtidos nessa dissertação contribuem nos trabalhos do GMCPAD, pois as paralelizações são adotadas em aplicações que estão sendo desenvolvidas no grupo.

Escadas de spin integráveis

Neste trabalho definimos três modelos de escadas de spin integráveis novos que correspondem a variações de um modelo de escada de spin baseado na simetria SU(4). Os modelos são exatamente solúveis através do método do ansatz de Bethe e as equações do ansatz de Bethe, os autovalores de energia e o gap de spin são derivados e propriedades físicas interessantes são discutidas. Inicialmente apresentamos um modelo de escada de spin integrável que possui um parâmetro livre além do acomplamento ao longo dos degraus. Determinamos a dependência do parâmetro anisotrópico na transição de fase entre uma região com gap e outra sem gap. Nós também mostramos que o modelo é um caso especial de uma Hamiltoniana mais geral que possui três parâmetros livres. A susceptibilidade magnética em função da temperatura é obtida numericamente e sua dependência no parâmetro anisotrópico é determinada explicitamente. Uma comparação entre o gap de spin obtido através da curva de susceptibilidade magnética e aquele obtido das equações do ansatz de Bethe é feita e uma boa concordância encontrada. A conexão com alguns compostos é apresentada e mostramos que os nossos resultados ajustam bem a curva da susceptibilidade magnética dos compostos KCuCI3, CU2(C5H12N2hC14e (C5H12NhCuBr4. A seguir nós propomos dois tipos diferentes de modelos integráveis com impurezas. Mostramos em ambos os casos que uma transição de fase entre uma região com gap e outra sem gap ocorre para um valor crítico do acoplamento ao longo dos degraus. Além disso, a dependência das impurezas na transição de fase é determinada explicitamente. Em um dos modelos o gap diminui com o aumento da intensidade da impureza A. E, fixando a intensidade de impureza A, é observada uma redução do gap com o aumento da concentração de impurezas. Este resultado está qualitativamente de acordo com resultados experimentais.