Simulação numérica da interação de escoamentos tridimensionais de fluidos compressíveis e incompressíveis e estruturas deformáveis usando o método de elementos finitos

Neste trabalho apresenta-se um algoritmo para a simulação de problemas tridimensionais de interação fluido-estrutura utilizando a técnica de elementos finitos. Um esquema de Taylor-Galerkin de dois passos e elementos tetraédricos lineares são empregados para o fluido, que pode ser compressível ou incompressível. É adotada uma formulação lagrangeana-euleriana arbitrária (ALE), compatível com o movimento da interface fluidoestrutura. Um método ftacionado de correção de velocidade é utilizado para os fluidos incompressíveis. A estrutura é analisada usando elementos triangulares com três nós e seis graus de liberdade por nó (três componentes de deslocamentos e três componentes de rotação). Os efeitos da não-linearidade geométrica são incluídos. O método de Newmark é empregado para integrar no tempo as equações dinâmicas de equilíbrio, usando-se uma descrição lagrangeana atualizada. O sistema de equações alge'bricas é solucionado através do método dos gradientes conjugados e o sistema não-linear, resultante de deslocamentos e rotacões finitas da estrutura, é solucionado com um esquema incremental-iterativo. O código é otimizado para aproveitar as vantagens do processamento vetorial.

Análise de indicadores voltados à tomada de decisão gerencial : um comparativo entre a teoria e a prática

Este trabalho visa abordar um aspecto de grande importância em um sistema de gestão empresarial – os indicadores que devem orientar a tomada de decisão na organização. Os indicadores podem cumprir duas funções administrativas: planejamento e controle. Por suas origens históricas, a ênfase no controle ainda prevalece sobre o planejamento. Porém, o ambiente competitivo dos dias atuais não permite atrasos na tomada de decisão gerencial. Portanto, se os indicadores auxiliarem na função de planejamento – orientando e permitindo inferir sobre possíveis consequências de uma ou outra decisão - estariam possibilitando ações mais focalizadas, com melhores resultados para a empresa. Diminuiriam, assim, eventuais correções no rumo da gestão. Consequentemente, diminuiria também o tempo necessário para o alcance dos resultados planejados. Esta dissertação buscou analisar os pr incipais modelos teóricos de sistemas de indicadores, explicitando quais os indicadores sugeridos e o modo como devem ser tratados. Realizou-se, também, uma pesquisa de campo em 20 empresas do Rio Grande do Sul, a fim de identificar os indicadores que estavam sendo efetivamente utilizados para a tomada de decisão, bem como problemas decorrentes do uso destes. Com isso, pôde-se confrontar os modelos teóricos e a prática atual das empresas pesquisadas, o que permitiu extrair as conclusões do trabalho. Além disso, permitiu verificar a correlação existente entre um uso eficiente de indicadores e o desempenho destas empresas, comparando-se com uma pesquisa realizada com este objetivo em empresas americanas. Como conclusão final, constatou-se a importância central dos indicadores em um sistema de gestão, constituindo uma base sólida na busca por um melhor desempenho empresarial.

Caracterização metalúrgica do aço SAE 1141 microligado ao nióbio por simulação física

Este trabalho analisa, através do emprego de simulador físico, os efeitos da adição do nióbio no aço SAE 1141 e seu comportamento como refinador de grão e endurecedor por precipitação, nas altas temperaturas de forjamento e tempos de reaquecimento aplicados industrialmente. Analisa também qual a influência das deformações e velocidades de resfriamento sobre o tamanho de grão austenítico, microestrutura, dureza, cinética da precipitação e propriedades finais. Na definição dos parâmetros de simulação foram observadas rotinas de produção de peças forjadas comercialmente, a literatura técnica, e possíveis alterações que otimizassem as propriedades do aço em estudo. Para tanto foram definidas as temperaturas de forjamento 1100, 1200 e 1250°C, a taxa de deformação de 1/s, as deformações de 30, 50 e 70%, e as velocidades de resfriamento de 36, 60, 100 e 150°C/min. Os resultados indicaram que a velocidade de resfriamento tem papel preponderante na microestrutura final, que variou de um agregado de ferrita e perlita para bainita (às vezes com alguma martensita) a medida em que a velocidade de resfriamento aumenta. A simulação térmica indicou que tamanho de grão austenítico aumenta a medida em que se eleva a temperatura de reaquecimento e que existe apenas um pequeno efeito ancorador de grão do nióbio nesse quesito, quando não estiver presente alguma deformação Foi observado, através da simulação física que, com a presença de deformação, principalmente para as temperaturas menores de reaquecimento, a precipitação induzida por deformação parece ancorar o crescimento dos grãos recristalizados, resultando num tamanho de grão austenítico menor. Esse trabalho constitui parte integrante do projeto de desenvolvimento da tecnologia dos aços microligados, num esforço conjunto do Laboratório de Metalurgia Física – LAMEF, através do Grupo de Desenvolvimento de Aços Microligados, da Universidade Federal do Rio Grande do Sul, em parceria com a Aços Finos Piratini e DANA – Albarus, com o apoio da CAPES.

Aplicação de um método de implantação da manutenção produtiva total a partir da ótica da teoria das restrições na linha de fabricação de panelas da Tramontina Farroupilha S.A.

Este trabalho tem por objetivo estudar a implantação dos princípios e técnicas da Manutenção Produtiva Total, Teoria das Restrições e Sistema Toyota de Produção em uma indústria metalúrgica visando aumentar a capacidade produtiva da linha de fabricação de panelas de aço inoxidável. O mesmo foi desenvolvido na Tramontina Farroupilha S.A. Após uma fundamentação teórica e a apresentação das características gerais da empresa, é proposto um método de implantação da Manutenção Produtiva Total em dez etapas, aplicado em equipamento gargalo e obtendo um aumento significativo de sua capacidade. As conclusões indicam a possibilidade de implantação deste método à esta atividade industrial bem como a sua utilização em outros ramos industriais, observadas as características específicas de cada organização.

Influência de efeitos não-locais na dispersão de poluentes na camada linite planetária

A equação de difusão-advecção é muito utilizada no campo de estudos da poluição atmosférica na determinação da concentração de poluentes. Uma maneira de solucionar o problema de fechamento desta equação está baseada na hipótese de transporte por gradiente que, em analogia com a difusão molecular, assume que o fluxo turbulento de concentração é proporcional à magnitude do gradiente de concentração média. Neste trabalho, diferentemente do modo tradicional, utiliza-se uma equação genérica para a difusão turbulenta considerando-se que o fluxo mais a sua derivada são proporcionais ao gradiente médio. Desta forma, obtém-se uma equação que leva em conta a assimetria no processo de dispersão de poluentes atmosféricos. Portanto, a proposta do presente trabalho é a obtenção da solução analítica desta nova equação utilizando-se a técnica da Transformada de Laplace, considerando-se a Camada Limite Planetária (CLP) como um sistema multicamadas. Os parâmetros que encerram a turbulência sâo derivados da teoria de difusão estatística de Taylor combinada com a teoria de similaridade convectiva válidos para grandes tempos de difusão. Finalmente, na avaliação da performance deste modelo que considera a assimetria no processo de dispersão de poluentes atmosféricos, utilizam-se os dados experimentais de Copenhagen e Prairie Grass.

Simulação computacional por dinâmica molecular de propriedades de equilíbrio e espectroscópicas de misturas líquidas CS2/C6H6

Propriedades de equilíbrio e espectroscópicas são determinadas para a mistura líquida CS2/C6H6 na temperatura de 298K e nas frações molares de benzeno iguais a 0,25, 0,50 e 0,75, usando simulação computacional por Dinâmica Molecular. As interações intermoleculares são descritas em relaçãao a todos os átomos por um potencial Lennard-Jones (12/6), com as interações eletrostáticas representadas por quadrupolos pontuais, localizados nos centros de massa das mol´eculas. Um segundo potencial também é considerado, onde a distribuição de cargas na molécula é usada para descrever as interações eletrostáticas. Dados termodinâmicos e propriedades estruturais, descritas através de funções de distribuição radial e de correlação angular, são calculados. Os resultados obtidos com os dois potenciais mostram boa concordância com os dados experimentais. A análise das frações molares locais e das energias internas indica que a mistura líquida CS2/C6H6 apresenta comportamento aproximadamente ideal. As funções de distribuição radial e de correlação angular evidenciam uma configuração preferencialmente paralela entre as moléculas de CS2 e C6H6 a distâncias menores em relação a distância do máximo da função de distribução , correspondente a primeira camada de coordenação . Nesta região, a configuração perpendicularé fracamente favorecida em relação as demais orientações. Funções de correlação temporal da polarizabilidade coletiva e da polarizabilidade molecular, correspondentes aos espalhamentos Rayleigh e Raman, respectivamente, são calculadas a partir do modelo de interação dipolo induzido por dipolo de 1a ordem. As funções de correlação são divididas nas contribuições orientacional, induzida por colisão e cruzada, e nas contribuições dos componentes. Propriedades de equilíbrio como anisotropia efetiva e intensidades integradas são determinadas para as misturas e líquidos puros. Os resultados indicam que a correlação orientacional entre as moléculas de CS2 é a maior responsável pelos espectros Rayleigh e espectros Raman deste componente. A maior participação do benzeno ocorre no espalhamento Raman com uma maior correlação na contribuição induzida por colisão. A separação na escala temporal entre a dinâmica reorientacional e induzida por colisão determinada na simulação é menor em relação ao experimento. Tempos de correlação da contribuição induzida por colisão são maiores nas misturas em relação aos líquidos puros no espalhamento Rayleigh, conforme o experimento, e estão relacionados a contribuições significativas das funções de correlação entre os componentes CS2 e benzeno a tempo longo. Funções de correlação temporal do momento dipolar coletivo são determinadas a partir do modelo dipolo induzido por quadrupolo e analisadas em termos de contribuições isotrópicas, anisotrópicas e cruzadas, como também de contribuições dos componentes. A relação entre estes termos é concordante com os parâmetros moleculares de anisotropia e quadrupolo dos componentes. Os resultados da simulação indicam absorbância de excesso para as misturas em uma ordem semelhante a do experimento. O tempo de correlação maior na mistura de fração molar de benzeno igual a 0,25, também encontrado no experimento, está relacionado a elevada contribuição a tempo longo da função de correlação entre dipolos induzidos em moléculas de CS2 por moléculas de C6H6. Os espectros calculados na simulação e os espectros experimentais do espalhamento de luz despolarizado Rayleigh apresentam boa concordância, validando os modelos de potenciais de interação e de indução usados. No entanto, diferenças são encontradas no infravermelho longínquo. Modi cações no modelo simplificado de indução são propostas.

Geração de modelos de co-simulação distribuída para a arquitetura DCB

O aumento na complexidade dos sistemas embarcados, compostos por partes de hardware e software, aliado às pressões do mercado que exige novos produtos em prazos cada vez menores, tem levado projetistas a considerar a possibilidade de construir sistemas a partir da integração de componentes já existentes e previamente validados. Esses componentes podem ter sido desenvolvidos por diferentes equipes ou por terceiros e muitas vezes são projetados utilizando diferentes metodologias, linguagens e/ou níveis de abstração. Essa heterogeneidade torna complexo o processo de integração e validação de componentes, que normalmente é realizado através de simulação. O presente trabalho especifica mecanismos genéricos e extensíveis que oferecem suporte à cooperação entre componentes heterogêneos em um ambiente de simulação distribuída, sem impor padrões proprietários para formatos de dados e para a descrição do comportamento e interface dos componentes. Esses mecanismos são baseados na arquitetura DCB (Distributed Co-Simulation Backbone), voltada para co-simulação distribuída e heterogênea e inspirada nos conceitos de federado (componente de simulação) e federação (conjunto de componentes) que são definidos pelo HLA (High Level Architecture), um padrão de interoperabilidade para simulações distribuídas. Para dar suporte à co-simulação distribuída e heterogênea, esse trabalho descreve mecanismos que são responsáveis pelas tarefas de cooperação e distribuição, chamados de embaixadores, assim como o mecanismo gateway, que é responsável pela interoperabilidade entre linguagens e conversão de tipos de dados. Também é apresentada uma ferramenta de suporte à geração das interfaces de co-simulação, que são constituídas de dois embaixadores configuráveis e um gateway para cada federado, gerado a partir de templates pré-definidos.

Modelagem e simulação da câmara de reação de uma caldeira a carvão pulverizado

O presente trabalho é parte integrante de um projeto multidisciplinar cujo objetivo é construir uma ferramenta computacional capaz de simular a operação de uma Usina Termoelétrica tratando, especificamente, da câmara de reação (local da caldeira onde se desenvolve a queima do combustível). Foi desenvolvido um modelo matemático dinâmico, pseudo-homogêneo, unidimensional e, com propriedades físicas variáveis, com o objetivo de descrever a conversão do calor gerado pela combust em energia térmica, que está baseado em balanços de massa e energia. Este modelo foi implementado no simulador gPROMS e seus parâmetros foram estimados utilizando o módulo gEST deste software (método da máxima verossimilhança), com base em dados de operação da Usina e dados cinéticos de literatura. Para a simulação dinâmica da câmara de reação foi utilizado um método de integração implícita de passos múltiplos com ordem e passo variáveis. A câmara de combustão foi dividida em três zonas, sendo realizados, em cada zona, os balanços de massa e energia, adotando-se a hipótese de mistura perfeita. O modelo desenvolvido prevê o perfil de temperatura ao longo da câmara, as correntes de saída e a composição dos gases de combustão em cada zona; incluindo uma previsão aproximada da concentração dos poluentes formados, principalmente monóxido de nitrogênio e dióxido de enxofre. Os resultados apresentados mostraram que o modelo prediz satisfatoriamente a composição do gás de combustão na saída da câmara de reação e que a divisão da câmara em diferentes zonas de combustão é apropriada para avaliar a geração e o aproveitamento energético na caldeira.

Simulação de reatores de polimerização de estireno : modelagem cinética e termodinâmica

Nesse trabalho, é proposto um modelo fenomenológico para a polimerização de estireno em solução, validado com dados de uma planta industrial. Um modelo matemático para a polimerização do estireno via radicais livres é desenvolvido para predizer o comportamento estacionário e dinâmico do processo contínuo de produção de poliestireno. O modelo inclui a cinética de polimerização, em dois reatores CSTR auto-refrigerados em série, o equilíbrio líquido-vapor na seção de reação e as propriedades do produto final. Os resultados obtidos através do modelo desenvolvido são utilizados para determinar quais são os efeitos da composição de alimentação, temperatura e concentração do iniciador na conversão do monômero, massas molares, polidispersão, formação de dímeros e trímeros, densidade, viscosidade, índice de fluidez, propriedades térmicas e tensão de ruptura. Para o equilíbrio líquido-vapor, um modelo termodinâmico também foi desenvolvido para reatores auto-refrigerados, que usam o calor de evaporação para remover o calor gerado por reações químicas fortemente exotérmicas, e modelado através de uma equação cúbica de estado e uma regra de mistura aplicada ao equilíbrio líquido-vapor a baixas pressões para soluções poliméricas através do emprego de um único parâmetro de interação. O modelo desenvolvido pode ser utilizado para avaliação de estruturas de controle, simulação e otimização do processo de polimerização do estireno.

Particionamento de grafos de aplicações e mapeamento em grafos de arquiteturas heterogêneas

Esta pesquisa visa a modelagem de clusters de computadores, utilizando um modelo analítico simples que é representado por um grafo valorado denominado grafo da arquitetura. Para ilustrar tal metodologia, exemplificou-se a modelagem do cluster Myrinet/SCI do Instituto de Informática da UFRGS, que é do tipo heterogêneo e multiprocessado. A pesquisa visa também o estudo de métodos e tecnologias de software para o particionamento de grafos de aplicações e seu respectivo mapeamento sobre grafos de arquiteturas. Encontrar boas partições de grafos pode contribuir com a redução da comunicação entre processadores em uma máquina paralela. Para tal, utilizou-se o grafo da aplicação HIDRA, um dos trabalhos do GMCPAD, que modela o transporte de substâncias no Lago Guaíba. Um fator importante é o crescente avanço da oferta de recursos de alto desempenho como os clusters de computadores. Os clusters podem ser homogêneos, quando possuem um arquitetura com nós de mesma característica como: velocidade de processamento, quantidade de memória RAM e possuem a mesma rede de interconexão interligando-os. Eles também podem ser heterogêneos, quando alguns dos componentes dos nós diferem em capacidade ou tecnologia. A tendência é de clusters homogêneos se tornarem em clusters heterogêneos, como conseqüência das expansões e atualizações. Efetuar um particionamento que distribua a carga em clusters heterogêneos de acordo com o poder computacional de cada nó não é uma tarefa fácil, pois nenhum processador deve ficar ocioso e, tampouco, outros devem ficar sobrecarregados Vários métodos de particionamento e mapeamento de grafos foram estudados e três ferramentas (Chaco, Jostle e o Scotch) foram testadas com a aplicação e com a arquitetura modeladas. Foram realizados, ainda, vários experimentos modificando parâmetros de entrada das ferramentas e os resultados foram analisados. Foram considerados melhores resultados aqueles que apresentaram o menor número de corte de arestas, uma vez que esse parâmetro pode representar a comunicação entre os processadores de uma máquina paralela, e executaram o particionamento/mapeamento no menor tempo. O software Chaco e o software Jostle foram eficientes no balanceamento de carga por gerarem partições com praticamente o mesmo tamanho, sendo os resultados adequados para arquiteturas homogêneas. O software Scotch foi o único que permitiu o mapeamento do grafo da aplicação sobre o grafo da arquitetura com fidelidade, destacando-se também por executar particionamento com melhor qualidade e pela execução dos experimentos em um tempo significativamente menor que as outras ferramentas pesquisadas.

Anomalia na densidade em um gás de rede com interações competitivas

Apesar de ser um líquido comum na natureza, muitas dúvidas ainda pairam sobre várias características da água. A existência de uma relação comum entre o tipo de potencial intermolecular, criticalidade e as várias formas de anomalia existentes nessa substância ainda é uma questão em aberto, apesar da intensa pesquisa que têm-se feito ao longo dos anos sobre esse assunto. Nesta dissertação, propomos a hipótese de que a anomalia na densidade esta correlacionada à presença de multicriticalidade e que ambos os fenômenos surgem de um potencial de duas escalas. Para dar suporte a esta hipótese, além de trabalhos anteriores, usamos e estudamos um gás de rede com interações que competem (primeiros vizinhos Vi atrativos e segundos vizinhos V2repulsivos). Construímos para este sistema um diagrama de fases J.lV8. T usando dois métodos: aproximação de campo médio e simulações. Encontramos na aproximação de campo médio duas linhas críticas, uma das quais encontra a linha de 1~ ordem separando duas fases líquidas, e um ponto tricrítico, dado que V2/V1 < -0.5. Se V2/V1 > -0.5, a transição líquido-líquido desaparece, dando lugar a apenas duas fases, uma líquido e uma gás, separadas por uma linha de coexistência terminada em um ponto crítico Com a aproximação de campo médio não encontramos anomalia na densidade. Os resultados obtidos com as simulações alteram qualitativamente o diagrama de fases. Tanto as linhas críticas quanto os pontos tricríticos tem suas posições modificadas com relação ao campo médio. Neste caso encontramos um comportamento anômalo na densidade se V2/Vl < -0.5. Concluímos que o potencial de duas escalas competitivas é um ingrediente necessário ao aparecimento de anomalia na densidade e coexistência entre duas fases líquidas. Ainda, mostramos que essa anomalia pode estar associada não apenas a dois pontos críticos, como se espera para a água, mas a uma multicriticalidade em geral, tal como linhas críticas.

Efeitos de ordenamento chiral em supercondutores e sistemas magnéticos desordenados

Este trabalho apresenta uma investigação experimental sobre as propriedades magnéticas e de transporte elétrico de sistemas caracterizados por desordem e frustração. Tais sistemas são amostras granulares do supercondutor de alta temperatura crítica YBa2Cu3O7-δ e amostras do tipo vidro-de-spin e reentrantes das ligas magnéticas diluídas AuMn 8at% e AuFe xat% (x = 8, 12, 15, 18 e 21). No supercondutor granular foram estudados os efeitos de flutuações termodinâmicas na magnetocondutividade nas proximidades da transição supercondutora e a linha de irreversibilidades magnéticas. A transição para o estado de resistência nula é um processo que ocorre em duas etapas. Inicialmente, a transição de pareamento estabiliza a supercondutividade no interior dos grãos. A transição de coerência ocorre em temperaturas inferiores e ativa ligações fracas entre os grãos através de um processo do tipo percolação. O regime que antecede a transição de coerência é caracterizado por flutuações na fase do parâmetro de ordem. A linha de irreversibilidades magnéticas, estudada a partir da magnetoresistência e magnetização DC, revela um comportamento do tipo Almeida-Thouless em baixos campos magnéticos aplicados, seguido de um crossover, em µ0H = 0.1 T, para um comportamento do tipo Gabay-Toulouse. A linha de irreversibilidades é interpretada como sendo uma manifestação experimental de uma transição de vidro chiral. As ligas magnéticas diluídas foram estudadas através do efeito Hall extraordinário. Os resultados mostram claramente que dois termos de sinais contrários contribuem para este efeito nestes sistemas. Em particular, é observada uma anomalia no coeficiente de Hall extraordinário em temperaturas próximas à temperatura de ordenamento dos vidros-de-spin e de “canting” dos reentrantes que não é prevista por nenhuma teoria convencional de efeito Hall. A interpretação dos resultados é feita em termos do modelo de vidro chiral, mostrando a importância de uma contribuição de origem puramente chiral. Esta é a primeira vez que a chiralidade, uma propriedade intrínseca de sistemas desordenados e frustrados, é observada experimentalmente de modo direto.

Decomposição de Calderon e suas aplicações na teoria da regularidade em equações elípticas

Este trabalho tem por objetivo estudar a regularidade de soluções de Equações Diferenciais Parciais Elípticas da forma Lu = f, para f 2 Lp(­), onde p > 1. Para isto, usamos a Decomposição de Calderon-Zygmund e um resultado que é consequência deste, o Teorema da Interpolação de Marcinkiewicz. Além disso, usando quocientes-diferença provamos a regularidade das soluções para o caso p = 2 e L = ¡¢ de uma forma alternativa.

Transporte espermático e resposta inflamatória na égua após a inseminação com diferentes concentrações de espermatozóides.

Normalmente, após a cobertura ou a inseminação artificial de éguas, ocorre uma endometrite aguda transitória em resposta ao sêmen e bactérias no útero. O objetivo deste estudo foi verificar se o transporte espermático e a intensidade da reação inflamatória uterina, 2h, 4h ou 24h após a inseminação com sêmen resfriado, são influenciados pela concentração espermática na dose inseminante. Para tal, foram utilizadas 192 éguas em cio, com folículo dominante ≥35 mm, sem crescimento bacteriano e livres de PMNs aos exames uterinos complementares. As éguas foram distribuídas aleatoriamente em grupos e inseminadas com 20 ml contendo 100x106 (n=30), 500x106 (n=27) ou 1000x106 (n=31) espermatozóides diluídos em solução de 3 ml de plasma seminal e 17 ml de leite desnatado, refrigerado e armazenado por 18 a 22 horas, ou infundidas com 20 ml de plasma seminal (n=33), ou com 20 ml de leite desnatado (n=38). As éguas foram abatidas duas, quatro ou 24h após as inseminações ou infusões. O grupo controle (n=33) não recebeu nenhum tratamento. Os ovidutos foram separados do útero, sendo útero e ovidutos lavados separadamente com PBS. Uma amostra do lavado de cada oviduto foi examinada para contagem de espermatozóides e uma amostra de cada lavado uterino foi utilizada para contagem de leucócitos. Após as lavagens, foi retirada uma amostra de endométrio para exame histopatológico. As éguas inseminadas e infundidas apresentaram reação inflamatória significativamente maior que as éguas do grupo controle, no decorrer das 24 horas. A reação inflamatória foi significativamente maior nas éguas inseminadas que nas infundidas. A reação inflamatória apresentou correlação com a concentração espermática (r=0,389). O número de éguas apresentando espermatozóides nos ovidutos não foi diferente nos grupos inseminados. Concluiu-se que componentes da dose inseminante provocam uma resposta inflamatória, sendo esta tanto mais severa e de resolução mais rápida, quanto maior for a concentração espermática. Por outro lado, até as quatro horas pós-inseminação, o transporte espermático independe da concentração espermática utilizada.