Simulação de sistemas térmicos de potência para geração de energia elétrica

O presente trabalho é dedicado à simulação numérica de sistemas térmicos de potência. O trabalho é iniciado com a modelagem de um ciclo Rankine, dedicado à produção de energia elétrica, para o qual foi elaborado um programa de simulação com a linguagem de programação MATLAB. A partir desse primeiro caso, são apresentados os modelos empregados para representar os diversos componentes que formam o circuito, como o gerador de vapor, a turbina, o condensador e a bomba. Além desses componentes, são introduzidas as equações que representam o escoamento do fluido de trabalho, no caso a água, permitindo assim o cálculo da perda de carga nas diferentes canalizações do circuito, sendo também acoplado o funcionamento da bomba. Essa alternativa pennite uma melhor avaliação do trabalho despendido para operar o sistema. A modelagem do ciclo deixa então de ser exclusivamente tennodinâmica, e passa a incluir aspectos de mecânica de fluidos. Outras variantes desse ciclo simples são também modelados e simulados, incluindo ciclos Rankine regenerativos e com irreversibilidades. As simulações são efetuadas admitindo-se parâmetros de operação, como, potência da turbina, temperatura do vapor d'água na entrada da turbina e pressão do vapor d'água na saída da turbina, com a variante de fixar-se o título do vapor d'água na saída da turbina.

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.

Reação-difusão com difusividade variável para oxidação de silício

Neste trabalho, propomos uma modificação do modelo de reação-difusão (R. M. C. de Almeida et al., Physics Review B, 61, 19 (2000)) incluindo difusividade variável com o objetivo principal de predizer, ou no mínimo descrever melhor, o crescimento de oxido de Si no regime de filmes nos. Estudamos o modelo reação-difusão a coeficiente de difusão, D, fixo e D variável. Estudamos extensivamente o modelo reação-difusão com D fixo caracterizando seu comportamento geral, e resolvendo numericamente o modelo com D variável para um intervalo amplo de relações DSiO2=DSi. Ambos casos apresentam comportamento assintótico parabólico das cinéticas. Obtivemos as equações analíticas que regem o regime assintótico de tais casos. Ambos os modelos apresentam interface não abrupta. Comparações das cinéticas com o modelo linear-parabólico e com dados experimentais foram feitas, e também para as espessuras da interface. Contudo nenhum dos dois modelos de reação-difusãao, com D fixo e com D varáavel, podem explicar a região de filmes nos, que possui taxa de crescimento superior a taxa de crescimento da região assintótica. Porém, ao incluir taxa de reação variável dentro do modelo de reação-difusão com D variável, este aponta para uma solução do regime de filmes nos.

A relação entre orientação para mercado, aprendizagem organizacional e performance

Apesar do grande interesse e da consolidação teórica do tema orientação para mercado e sua influência no desempenho organizacional, uma corrente de autores da área de Marketing tem argumentado que a orientação para mercado é necessária, mas não suficiente para sustentar uma vantagem competitiva de longo prazo. Esta corrente de autores defende a habilidade complementar da organização de se engajar em um processo contínuo de aprendizagem para a sustentação de uma posição distintiva no ambiente competitivo. Visando um aprofundamento deste tema, esta tese teve como objetivo principal a avaliação da relação da orientação para mercado com a performance empresarial, sob a influência de uma postura de aprendizagem organizacional. Para tanto, um modelo teórico, contendo as relações hipotetizadas entre os referidos construtos, foi desenvolvido, testado e ajustado, através da aplicação da técnica de modelagem de equações estruturais. Os resultados do levantamento realizado na indústria Eletro-Eletrônica do Brasil demonstraram uma influência positiva e significativa da orientação para mercado sobre a performance empresarial, além do impacto indireto da orientação para aprendizagem na performance empresarial através da sua forte e positiva influência sobre a orientação para mercado. Resultados complementares, futuras pesquisas e implicações gerenciais foram, ainda, discutidos.

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.

Previsão de falhas no período de garantia a partir da percepção do cliente : um estudo de caso junto a uma empresa de alta tecnologia

A prestação de serviços tem sido um importante diferencial dentro do cenário competitivo atual. Para empresas de manufatura, um desses serviços de destaque é o atendimento à garantia, o qual deve ser desenhado de forma a maximizar a satisfação do cliente, operando dentro de custos otimizados. Para permitir a alocação mais eficiente de recursos aos serviços de atendimento à garantia, é importante entender o comportamento das falhas do produto durante o período de garantia, o qual depende da confiabilidade do produto, do tempo e intensidade de uso e do volume de vendas. Esta dissertação trata da aplicação e adaptação de um modelo apropriado para a predição de falhas durante o período de garantia. Esse modelo poderá auxiliar as organizações a manter controle sobre o custo do atendimento à garantia, auxiliando na melhoria contínua dos serviços prestados.

O desenvolvimento das trejetórias do comportamento delinqüente em adolescentes infratores

O presente trabalho teve como objetivo investigar, através de um delineamento longitudinal retrospectivo e prospectivo, as trajetórias do comportamento delinqüente, em adolescentes infratores brasileiros. Também se investigou variáveis individuais, familiares e sociais, consideradas fatores de risco para a delinqüência. Participaram 200 adolescentes, do sexo masculino, de baixo nível sócio-econômico, que cumpriam medidas sócio-educativas pelo cometimento de atos infracionais. Os adolescentes foram entrevistados individualmente e sua trajetória de vida foi acompanhada durante dois anos. Os resultados mostraram que em 72,5% dos adolescentes investigados ocorreu a desistência do comportamento delinqüente, e que em 26,5% ocorreu a persistência do comportamento delinqüente. Durante a realização da pesquisa dois adolescentes foram mortos. Análises de regressão mostraram que os melhores preditores destas trajetórias foram variáveis contextuais: trabalhar antes do cometimento dos atos infracionais, o número de internações em unidades de proteção durante a infância e a trajetória de inserção no sistema de atendimento à criança e/ou adolescente em situação de risco. Os resultados indicaram também a presença de características pessoais saudáveis nos jovens investigados, tais como um estilo atribucional interno, estratégias de coping adaptativas, altos níveis de auto-estima e baixos índices de depressão. Estes resultados devem ser levados em conta nas estratégias de intervenção com adolescentes em situação de risco social e pessoal para a delinqüência. Também confirmam a importância de intervenções que tenham como foco não somente os jovens, mas também sua comunidade, sua família, sua escola e o planejamento de programas que incrementem oportunidades de emprego.

Determinação dos campos de velocidade e temperatura em ambientes ventilados

No presente trabalho apresenta-se um estudo numérico do escoamento de ar em um ambiente interno tridimensional típico, sujeito a um sistema de ventilação por deslocamento. O estudo numérico consiste na obtenção dos campos de velocidade e temperatura presentes no interior deste ambiente. A determinação dos campos de velocidade e temperatura é necessária na avaliação precisa das condições de conforto térmico, produzidas por um determinado sistema de climatização. Este tipo de estudo não serve somente para determinar a viabilidade da aplicação de um determinado sistema de distribuição de ar, mas através dele também se pode determinar as condições ótimas de uso deste sistema. De forma crescente observa-se o interesse dos projetistas de ar condicionado em aperfeiçoar as suas técnicas de cálculo, de modo que a simulação numérica poderá tornar-se uma técnica usual em pouco tempo. Atualmente duas técnicas estão disponíveis para o estudo do escoamento de ar em ambientes internos: investigação experimental e simulação numérica. Em princípio, medidas diretas dão informações mais realísticas relativas ao escoamento do ar interno, contudo, medições experimentais são muito caras e consomem muito tempo, inviabilizando a sua aplicação caso a caso. Devido às limitações da investigação experimental e à disponibilidade de computadores cada vez mais rápidos, o uso da simulação numérica surge como uma opção considerável nos dias de hoje O método de solução numérica, utilizado no presente trabalho, baseia-se na resolução das equações de conservação de massa, de quantidade de movimento e de energia, em um domínio de cálculo discretizado com o uso do Método dos Volumes Finitos. Os efeitos da turbulência são levados em conta através da utilização do modelo de turbulência K- ε e de um modelo de turbulência simplificado. Resultados bem convergidos para o problema abordado somente são obtidos com a utilização do modelo de turbulência simplificado. As simulações numéricas apresentadas no presente trabalho foram realizadas com o uso do código Fluent. Com base na comparação dos resultados numéricos obtidos no presente trabalho e os resultados experimentais de Yuan et al., 1999a, pode-se concluir que o modelo simplificado é adequado, ao menos provisoriamente, para determinar a distribuição de velocidade e de temperatura no escoamento de ar em ambientes internos.

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.

Simulação numérica de um sistema de armazenamento de calor em estufas de plasticultura

Apresentar um modelo para simular um sistema de armazenamento de calor no solo em estufas para plasticultura é o objetivo do presente trabalho. O sistema consiste num feixe de tubos enterrados no solo. A convecção forçada de ar no seu interior realiza a troca térmica necessária para manter as estufas sob faixas desejadas de temperatura. O objetivo do modelo é investigar os efeitos no calor armazenado e a influência das variáveis, tais como diâmetro, comprimento, espaçamento entre os tubos e a velocidade de ar no canal provocam no sistema. O solo é tratado como um meio difusivo e avalia-se a contribuição do termo de condensação e evaporação da água contida no ar em escoamento nos tubos. A equação da energia é resolvida para o solo e para o ar. Os tubos de seção transversal circular são modelados como tubos de seção transversal quadrada com o objetivo de que as simulações possam ser processadas em coordenadas cartesianas. O programa resolve situações tridimensionais, transientes e emprega o Método dos Volumes Finitos para integrar as equações diferenciais governantes. O modelo original é baseado no modelo de Gauthier et al., 1997, tendo sido os resultados do mesmo foram usados para a validação do presente estudo. Um circuito de água quente é também projetado e apresentado para o aquecimento das estufas. A água circula através de mangueiras sobre o solo e é aquecida por um sistema de queimadores a gás liqüefeito de petróleo ou óleo combustível, transferindo assim calor para o interior da mesma. O projeto de aquecimento foi realizado através de um programa de parceria entre a Ufrgs, Sebrae, Fapergs e a Agropecuária Clarice.

Otimização de geometrias aerodinâmicas utilizando métodos inversos

O objetivo deste trabalho é a obtenção de uma técnica para a modelagem otimizada de corpos submetidos a fluxos de alta velocidade, como aerofólios em escoamentos transônicos e outras geometrias aerodinâmicas. A técnica é desenvolvida através de expansões em séries de Fourier para um conjunto de equações diferenciais com interrelação com as condições de contorno, sendo uma equação para a parte superior e outra para a parte inferior do aerofólio. O método de integração temporal empregado baseia-se no esquema explícito de Runge-Kutta de 5 estágios para as equações da quantidade de movimento e na relação de estado para a pressão. Para a aproximação espacial adota-se um esquema em volumes finitos no arranjo co-localizado em diferenças centrais. Utiliza-se dissipação artificial para amortecer as frequências de alta ordem do erro na solução das equações linearizadas. A obra apresenta a solução de escoamentos bi e tridimensionais de fluidos compressíveis transônicos em torno de perfis aerodinâmicos. Os testes num´ericos são realizados para as geometrias do NACA 0012 e 0009 e asas tridimensionais usando as equações de Euler, para número de Mach igual a 0.8 e ® = 0o. Os resultados encontrados comparam favoravelmente com os dados experimentais e numéricos disponíveis na literatura.

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.

Simulação numérica de processos de solidificação em sistemas binários aplicados à criopreservação de células

A preservação e o armazenamento de células e tecidos têm sido utilizados largamente em pesquisa científica e aplicações clínicas. No entanto, há uma aparente contradição entre o conceito de preservaão e as conclusões baseadas em resultados experimentais que materiais biológicos criopreservados podem ser danificados pelo próprio processo de preservação. A compreensão do processo de solidificação de soluções salinas é fundamental para a proposição de novos protocolos de criopreservação. No presente estudo, o congelamento de uma solução de cloreto de sódio a 1% em massa é simulado. As equações de conservação de massa, momentum, energia, e espécies químicas foram discretizadas e resolvidas numericamente utilizando-se o método dos volumes de controle para um domínio bidimensional que contém a parede da bolsa plástica e a solução salina. A perda de água da célula foi calculada a partir da história de temperatura e concentração durante o processo de solidificação e verificou-se que, dependendo da posição inicial da célula na bolsa, a célula tem probabilidades diferentes de sobreviver durante o processo.

Estudo da instabilidade de feixes de condutores em linhas aéreas de transmissão

A demanda de energia elétrica tem aumentado consideravelmente com o observado crescimento da atividade industrial, exigindo o desenvolvimento de novas tecnologias para aumentar a capacidade de transporte de energia por linha. Neste sentido, as linhas de transmissão de alta voltagem, anteriormente formadas por condutores isolados, podem ter a sua capacidade expressivamente incrementada dispondo-se os condutores em feixes, ou seja, com mais de um condutor por fase. A ação do vento nos condutores em forma de feixes originam problemas vibratórios diferentes daqueles comumente encontrados em condutores isolados. Além de oscilações por galope e desprendimento de vórtices, feixes são suscetíveis a oscilações de baixa freqüência causada pelo efeito de interferência das esteiras entre os condutores. Estas oscilações são de dois tipos: (i) movimento de sub-vão, que é característico do condutor a sotavento, o qual pode-se deslocar independentemente dos outros condutores vizinhos, e (ii) movimento de vão completo, no qual movimentos do feixe ocorrem no vão compreendido entre as torres de suporte. Alguns trabalhos na literatura consideram os movimentos de vão completo em condutores coberto por gelo, mas pouca informação se tem sobre o estudo de oscilações em feixes sem a presença de gelo. Devido a complexidade destes fenômenos, não há ainda critérios claros quanto à estabilidade do feixe para as configurações usuais de linhas deste tipo. Este trabalho tem então como objetivo inicial apresentar as principais características dinâmicas de condutores isolados e em feixes. Além disso, apresentar uma análise dos parâmetros que influenciam os movimentos, dentre os quais são citados: características do terreno, velocidade e turbulência do vento, número e arranjo dos condutores, rugosidade do condutor, espaçamento entre condutores, inclinação do feixe, sistemas de espaçadores e de suspensão, tração do condutor, freqüências naturais e efeito de flecha. Ensaios em túnel de vento são normalmente realizados com o propósito de obter-se dados de tipos específicos de condutores e configurações de feixes de condutores. Posteriormente, estes dados podem ser utilizados na análise teórica dos movimentos dos condutores a fim de prever-se a estabilidade dos mesmos. Neste trabalho foram projetados experimentos do tipo estático, que permitiram a obtenção dos coeficientes aerodinâmicos e suas derivadas em relação aos ângulos de incidência do vento em feixes de cabos. Para isto utilizaram-se duas células de carga com extensômetros, onde em uma delas mediu-se diretamente as forças de arrasto e de sustentação e na outra os momentos de torção. Os experimentos foram conduzidos no túnel de vento Prof. Joaquim Blessmann, da UFRGS. Os modelos eram compostos por condutores rígidos de pequeno comprimento que estavam dispostos isoladamente, em feixe de 02 cabos dispostos lado-a-lado ou em feixes de 04 cabos dispostos na configuração de um quadrado. Os ensaios foram realizados para diferentes espaçamentos entre os cabos (10, 14, 18 e 22 diâmetros), para diversos ângulos de inclinação do feixe (0o a 45o) e ainda para diferentes velocidades e condições de escoamento (suave ou turbulento). Os ensaios foram realizados com cabos lisos e cabos do tipo Rook ACSR 24/7, a fim de analisar-se a influência da rugosidade no comportamento dos modelos. Adicionalmente, foi desenvolvida uma aplicação da análise do problema de instabilidade dinâmica através da utilização de equações linearizadas do movimento de condutores em feixe. A partir destas equações, determinam-se as regiões de instabilidade das oscilações de vão completo, em feixes de dois e quatro condutores. Os coeficientes aerodinâmicos utilizados nestas equações foram aqueles determinados nos ensaios no túnel de vento Prof. Blessmann. Finalmente, o conhecimento da influência dos parâmetros vinculados as características do feixe, vão e vento incidente nos coeficientes aerodinâmicos e na conseqüente estabilidade do feixe, possibilita a determinação de alguns critérios de projeto que garantam maior estabilidade dos feixes de condutores.

Aplicações do método dos elementos discretos em estruturas de concreto

A ciência moderna apresentou significativo avanço a partir do desenvolvimento da análise diferencial. A transformação de equações diferenciais de alta ordem em sistemas de equações algébricas foi possível através do desenvolvimento de métodos numéricos, constituindo este, outro grande avanço. Dentro desses pode-se destacar os métodos de diferenças finitas, dos elementos finitos, dos elementos discretos e mais recentemente, os elementos de contorno. Neste trabalho, faz-se uma contribuição ao desenvolvimento do Método dos Elementos Discretos para aplicações na Mecânica do Contínuo, na Mecânica da Fratura, assim como na determinação do dano em elementos estruturais submetidos a cargas. Neste método, a discretização espacial no modelo se realiza mediante um conjunto de massas ligadas entre se por forças materializadas como um arranjo de barras de treliça com rigidez equivalente ao contínuo que se quer representar, e mediante um esquema de integração explícita, se realiza a integração das equações de movimento no tempo. Verifica-se a validade e a capacidade do método em predizer o efeito de tamanho em elementos de concreto e concreto armado, obtendo-se uma excelente correlação com ensaios encontrados na literatura técnica, além de importantes conclusões a respeito da aplicação de cargas estáticas e dinâmicas, tanto em padrões de fissuração ou ruptura, quanto aos valores limites de resistência dos materiais ou cargas aplicadas, dando-se importância na geração aleatória das propriedades dos materiais mediante o uso do Método de Representação Espectral.