Aproximações de Galerkin clássico e Galerkin/mínimos quadrados de escoamentos viscosos incompressíveis : aplicação à medição de fluxos

Desde a antigüidade a medição do escoamento dos fluidos tem sido uma marca de nossa civilização, ajudando a predizer a fertilidade das terras e o consumo d’água em fontes e aquedutos. Nos nossos dias, a área de medição de fluxo está bem estabelecida e ainda desperta grande interesse nas linhas de pesquisa da mecânica dos fluidos experimental e computacional. Em particular, o estudo da medição de fluxo com elementos intrusivos, tais como placas de orifício, é de grande interesse dado o preço baixo do medidor, e sua boa precisão quando comparada à sua simplicidade de instalação e manutenção. Esta dissertação tem como objetivo o estudo da aplicação de elementos finitos aos escoamentos de fluidos viscosos - via aproximação clássica de Galerkin e Galerkin/mínimos-quadrados (GLS) – com particular ênfase na aproximação das equações de Navier-Stokes incompressível no escoamento newtoniano através de um canal obstruído por uma placa de orifício. Inicialmente, são apresentadas as dificuldades do método de Galerkin clássico na aproximação de escoamentos incompressíveis; ou seja, através da simulação de escoamentos viscosos bem conhecidos - como o escoamento no interior de uma cavidade e através de uma expansão súbita - fica evidenciada a restrição imposta pela condição de Babuška-Brezzi quando da escolha dos subespaços aproximantes do campo de velocidade e pressão Como alternativa às patologias do método de Galerkin clássico, esta dissertação emprega a metodologia de Galerkin/mínimos-quadrados na simulação acima mencionada da placa de orifício, a qual permite o uso de elementos de igual-ordem para aproximar velocidade e pressão e capturar de maneira estável escoamentos sujeitos a altos números de Reynolds. Os testes computacionais realizados se apresentaram fisicamente realistas quando comparados com a literatura e dados experimentais, sendo todos desenvolvidos no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) do Departamento de Engenharia Mecânica da Universidade Federal do Rio Grande do Sul.

Planejamento de perturbações para a identificação de modelos dinâmicos

A identificação de modelos é determinante no sucesso das modernas técnicas de controle avançado de processos. Um modelo para o sistema pode ser obtido através de modelagem rigorosa, baseada em equações governantes do sistema ou através da modelagem empírica e estimação de parâmetros. Embora mais rápida e fácil, a modelagem empírica necessita de alguns testes de identificação nos quais as variáveis manipuladas são variadas de modo que resultem em variações nas variáveis controladas. Os testes de identificação podem apresentar custos muito elevados tendo em vista que o sistema pode sair de seu ponto normal de operação, gerando produtos com folga de especificação. Este fato ocorre porque usualmente as perturbações aplicadas nas variáveis manipuladas nas indústrias de processos são independentes umas das outras, aumentando a duração do teste de identificação. Desta forma, neste trabalho foi desenvolvida uma nova metodologia de projeto de perturbações simultâneas para a identificação de modelos dinâmicos baseada na direcionalidade do sistema, com o objetivo de fornecer dados mais ricos para se capturar corretamente o comportamento multivariável do sistema e manter o processo no ponto de operação normal. As perturbações são projetadas conforme as características de um modelo simplificado do processo, ou pré-modelo. Este modelo inicial é obtido essencialmente de dados históricos de planta, selecionados através de uma sistemática análise de correlação desenvolvida neste trabalho A metodologia proposta é composta de duas partes: a primeira parte diz respeito à análise dos dados históricos de planta para obtenção de informações prelimirares as quais são utilizadas no planejamento de perturbações, tais como amplitude do ruído de medida, correlação entre as variáveis de processo, constante de tempo do sistema e matriz de ganhos. E a segunda parte consiste no cálculo da amplitude das perturbações baseado nos resultados da primeira etapa do planejamento. Para sistemas mal-condicionados verificou-se que as perturbações planejadas pela metodologia removem menos a planta de seu ponto de operação gerando resultados mais consistentes em relação às perturbações tradicionais. Já para sistemas bem-condicionados, os resultados são semelhantes. A metodologia foi aplicada em uma unidade piloto experimental e numa unidade de destilação da PETROBRAS, cujos resultados apontam pouca remoção dos sistemas do ponto de operação e modelos consistentes. A validação dos modelos também foi contemplada na dissertação, uma vez que foi proposto um novo critério de validação que considera a derivada dos dados de planta e a do modelo e não apenas os dados de planta e os dados da simulação das saídas do modelo.

Análise de transformações geométricas para o georreferenciamento de imagens do satélite CBERS-1

Neste trabalho serão descritos métodos e técnicas na análise de transformações geométricas para o georreferenciamento de imagens do satélite CBERS I, utilizando o sensor CCD na região de Porto Alegre, com a utilização das transformações afim e equações projetivas e com o uso de pontos de apoio coletados com receptor GPS. Os resultados experimentais obtidos com as transformações afim e equações projetivas são animadores, recomendando-se então continuar os estudos para as imagens do CBERS-I que pode ser um atalho importante para atualizar a cartografia regional brasileira, pois neste caso do sensor CCD se vislumbra a possibilidade de se gerarem cartas imagens nas escalas de 1:100.000 e 1:50.000.

Simulador em regime transiente de alguns equipamentos de um ciclo de Rankine regenerativo

Este trabalho apresenta um modelamento computacional dos diversos componentes de um ciclo de Rankine regenerativo em operação transiente, utilizando as leis fundamentais de balanço de massa e energia, além de algumas equações constitutivas. O cálculo das propriedades termodinâmicas foi realizado através de rotinas desenvolvidas no Grupo de Estudos Térmicos e Energéticos (GESTE) da UFRGS. Os componentes com comportamento inercial são o condensador e o tanque de água de alimentação. Os demais respondem instantaneamente. O método de solução das equações algébricas é seqüencial e o avanço no tempo é totalmente implícito. Primeiramente o sistema é colocado em regime permanente, a partir do qual é introduzida uma mudança de carga na turbina em um determinado intervalo de tempo. Algumas variáveis prescritas foram obtidas a partir de dados reais de uma usina termelétrica. A partir desse modelamento foi possível investigar o comportamento dinâmico do condensador e do tanque de alimentação e de que maneira estes efeitos eram percebidos pelos demais componentes. Os resultados obtidos demonstram que ambos possuem um comportamento dinâmico, sendo que o condensador apresenta uma constante de tempo menor que o tanque de alimentação. O condensador, por rejeitar o calor necessário, acaba por perceber rapidamente as oscilações introduzidas no ciclo Além disto é apresentado um método para determinação de coeficientes globais de troca térmica em condensadores e regeneradores. Os resultados obtidos neste trabalho foram utilizados na modelagem desses componentes durante a construção de um simulador em regime permanente da UTE Presidente Médici - Fase B. Também foi realizada uma pesquisa sobre os diversos programas simuladores de plantas de potência, assim como uma discussão das vantagens introduzidas por simuladores dinâmicos e em regime permanente nas decisões de engenharia.

Equação de burgers propriedades e comportamento assintótico

No presente trabalho, obtemos e analisamos diversas propriedades das soluções u(·, t) da equação de difusão linear (equação do calor em meios unidimensionais homogêneos) ut = μuxx x 2 R, t > 0 correspondentes a estados iniciais u(x, 0) = u0(x), com u0 2 Lp(R), para algum 1 p < 1; bem como da equação de Burgers ut + cuux = μuxx x 2 R, t > 0 onde c, μ são constantes dadas, sendo c 6= 0 e μ > 0 e ainda assumindo u(x, 0) = u0(x) com u0 2 Lp(R) para 1 p < 1, e limitado. Estudamos também a equação mais geral da forma ut + f(u)x = μuxx x 2 R, t > 0 discutindo várias propriedades importantes das soluções, associadas a estados iniciais u0 2 Lp(R) \ L1(R) para algum 1 p < 1. Em particular, examinamos o comportamento de ku(·, t)kLr(R), p r 1, para t >> 1, e diversas propriedades relacionadas.

A prática do marketing de relacionamento e a retenção de clientes: um estudo aplicado em um ambiente de serviços

Este trabalho aborda a estratégia da retenção de clientes por meio da prática do marketing de relacionamento no ambiente de serviços. É proposto um Modelo Teórico, com suas respectivas hipóteses, o qual é comparado a dois Modelos Rivais, identificados a partir da percepção de executivos de uma empresa com atuação proeminente no setor de serviços no Brasil, a Ticket, uma empresa do Grupo Accor. Nesse contexto, os construtos Satisfação de Clientes, Valor, Reputação do Provedor de Serviços e Confiança são considerados antecedentes da Retenção de Clientes. A análise central dos dados foi implementada por técnicas de estatística multivariada e pela Modelagem de Equações Estruturais. Os resultados indicam que o Modelo Teórico Reespecificado apresenta o melhor ajuste em relação aos demais modelos avaliados, embora precise de um maior amadurecimento teórico e empírico, até pelo ineditismo do modelo integrado. O trabalho propicia evidências relevantes, comprovando que o valor percebido pelos clientes é um construto antecedente da confiança depositada no provedor de serviços, relação esta mediada pela reputação do provedor de serviços; que a reputação do provedor de serviços é importante no contexto relacional em estudo; e que a retenção de clientes é positivamente influenciada pela confiança depositada no provedor de serviços.

Efeito do nível e fonte de fibra sobre a concentração e a utilização da energia metabolizável de dietas para frangos de corte em crescimento

A utilização de ingredientes fibrosos pode resultar efeitos deletérios sobre a energia metabolizável (EM) e líquida (EL) da dieta para frangos que dependem do nível, fonte e estrutura macromolecular da fibra. No presente trabalho foram desenvolvidos dois experimentos com frangos em crescimento, com os objetivos de avaliar os efeitos do nível e fonte de fibra sobre a disponibilidade e partição da EM, testar diversas medidas da fibra em equações de predição da concentração dietética de EM e EL e avaliar e aperfeiçoar metodologias de estudo da partição metabólica da energia. No experimento I, utilizando abate comparativo e regressão múltipla com frangos ISA em gaiolas individuais, foram observados efeitos significativos do nível e fonte de fibra sobre a disponibilidade e eficiência de retenção da EM (K0). As melhores medidas isoladas do efeito da fibra foram a parede celular insolúvel em água e a fibra total da dieta, e as melhores equações de predição obtidas para a EM e a EL contiveram combinações de uma medida da quantidade de fibra, como as fibra bruta, detergente ácido ou uma das citadas anteriormente, com as medidas de viscosidade e/ou de retenção de água. No experimento II, utilizando calorimetria indireta e modelagem com frangos Ross alojados em grupo, foi demonstrado que o efeito do nível de fibra sobre o K0 em dietas com baixa viscosidade pode estar relacionado com o aumento no efeito térmico do alimento (ETA) de curto prazo. A perda de calor total (61% dos 360 kcal EM/kg0,75 consumidos) e partida em metabolismo basal (52%), atividade física (19%) e ETA total (29%) e de longo prazo (16%), não foram alterados pelo nível de fibra. O peso do trato gastrintestinal das aves aumentou com o aumento no nível de fibra nos dois experimentos Foram obtidos coeficientes médios para a EL de mantença de 97,5 e 115,2 kcal/kg0,75/dia, respectivamente para os experimentos I e II, e o K0 variou de 0,638 (alta fibra de trigo e far. soja) a 0,708 (controle milho e far. soja) no experimento I, e foi de 0,706 (dieta normal) e 0,714 (alta fibra) no experimento II.

Desenvolvimento de modelos simplificados para o estudo da extrusão de polímeros

Neste trabalho, desenvolveram-se modelos matemáticos simplificados para o cálculo de duas variáveis importantes no estudo da extrusão de polímeros: vazão mássica (M&) e pressão na saída da extrusora (Psaída), em função das propriedades dos materiais e das condições de operação do equipamento. Podem-se utilizar esses modelos como ferramentas simples para a definição de critérios de ajustes que se devem fazer em uma extrusora de parafuso único para obter-se o desempenho desejado quando se alimenta o equipamento com um novo material. Para desenvolverem-se os modelos simplificados, utilizaram-se dados experimentais da extrusão de poliestireno (PS) e de polipropileno (PP), bem como resultados preditos por um programa computacional de simulação de extrusão disponível comercialmente. Mediram-se os dados experimentais de vazão mássica e de pressão na saída da extrusora em um equipamento de parafuso único de 45 mm de diâmetro. Realizaram-se esses testes, variando-se a velocidade de rotação do parafuso de 70 a 100 rpm para ambos os polímeros. No primeiro conjunto de simulações, utilizou-se o simulador Flow 2000 (Compuplast Inc.) para ajustarem-se os valores preditos de M& e de Psaída aos dados obtidos experimentalmente através da estimação dos fatores de atrito barril-polímero tanto para o PP quanto para o PS. Posteriormente, realizou-se um planejamento de experimentos, do tipo fatorial fracionado , para obter-se um segundo conjunto de simulações, considerando-se as propriedades dos materiais (reológicas e térmicas) e as condições de operação da extrusora (velocidade de rotação do parafuso e perfil de temperatura nas zonas de aquecimento da extrusora) como fatores de investigação. Com as novas simulações no Flow 2000, ajustaram-se os parâmetros dos modelos simplificados aos valores de vazão mássica e de pressão na saída da extrusora preditos no simulador. Elaboraram-se os modelos simplificados levando-se em conta as interações entre os fatores cujos efeitos consideraram-se significativos nas análises de variância (ANOVA). Obteve-se um modelo linear com 37 termos para o cálculo da vazão mássica e um modelo linear com 41 termos para o cálculo da pressão na saída da extrusora. Posteriormente, aplicou-se uma técnica de regressão multivariável para selecionar apenas os termos importantes dessas 1402IV2− XVI equações, conduzindo a um modelo linear com 10 termos para o cálculo da vazão mássica e a um modelo com 6 termos para o cálculo da pressão na saída da extrusora. Conseguiu-se boa concordância entre os dados experimentais e os valores preditos quando se aplicaram os modelos simplificados.

População estelar em galáxias ativas versus não ativas : o papel das novas gerações de estrelas

Esta tese apresenta o estudo da população estelar nuclear e extranuclear de uma amostra de galáxias ativas próximas (trinta e sete galáxias Seyfert 2 e vinte e quatro rádio-galáxias); e de uma amostra de controle de galáxias não ativas (onze elípticas: sete lenticulares e dezoito espirais). Foram para isto utilizados espectros óticos de fenda longa com boa razão sinal-ruído obtidos em telescôpios com aberturas de 1;5m a 4m. As regiões amostradas nas extrações correspomlem nas galáxias a 100-2000 parsecs (pc), com um valor mediano de 800 pc. A fim de verificar se existe relação entre a presença de formação estelar recente e a atividade nuclear; foram determinadas as idades das populações estelares da região nuclear e circumnuclear destas galáxias; utilizando o método de síntese espectral. Também foi verificada a infiuência do núcleo ativo sobre a população circunmclear; através do estudo da diluição das larguras equivalentes nucleares em relação àquelas de fora do núcleo. Um dos aspectos mais relevantes do presente estudo é a inclusão de uma amostra de controle de galáxias não ativas em um nÚmero comparável ao de galáxias ativas. O estudo conjunto de uma amostra de controle serviu para quantificar as diferenças obtidas dos espectros - em particular no contínuo e na população estelar das galáxias devido a existência da atividade nuclear. Os principais resultados são os que seguem: Uma grande fração de galáxias Seyfert 2 apresenta formação estelar recente (com idades iguais ou inferiores a 100 milhões de anos), tanto no núcleo quanto fora dele. Por outro lado, as rádio-galáxias em geral são dominadas por populações velhas e de idade intermediária (10 e 1 bilhão de anos); as populações jovens são significativas em apenas 10% destes objetos Em vários aspectos, galáxias Seyfert 2 e rádio-galáxias apresentam características diferentes. As primeiras apresentam resultados da síntese bem mais diversificados do que as últimas; tanto em termos de contribuição das populações de diferentes idades, quanto em relação ao comportamento dos gradientes de população estelar e dos avermelhamentos internos das galáxias. As galáxias Seyfert 2 apresentam uma diversidade de populações bastallte grande; e em geral essas populações são muito diferentes daquelas encontradas em galáxias não ativas lenticulares e espirais. Apesar das rádio-galáxias apresentarem uma pequena diversidade de populações estelares e, aparentemente, populações semelhantes àquelas encontradas em galáxias não ativas elípticas e lenticulares, elas têm populações levemente mais jovens do que as encontradas nas galáxias não ativas, sob a forma de uma contribuição maior da população de 1 bilhão de anos. Nenhuma rádio-galáxia Fanaroff-Riley tipo I apresenta contribuição significativa de populações mais jovem; do que 1 bilhão de anos: enquanto que em uma ou duas rádio-galáxias FRI isso ocorre. Esse resultàdo sugere que existe uma pequena diferença entre a população estelar das rádio-galáxias FRI e FGII estudadas. Essa diferença precisa ser melhor estudada, através de uma amostra de rádio-galáxias maior, para se verificar se não é originada apenas pelo pequeno número de objetos analisados em cada grupo Para as galáxias Seyfert 2, os resultados encontrados são consistentes com um cenário evolutivo 1 onde uma interação provocaria a queda do gás na direção do núcleo provocando um ou mais episódios de formação estelar e iniciando também a atividade nuclear. Enquanto a formação estelar é dominante, as assinaturas da interação são ainda visíveis e o espectro apresenta as características de uma população jovem e linhas de emissão intemediárias entre aquelas de galáxias Seyfert 2 e de galáxias "Starlmrst". A partir do momento em que o episóidio de formação estelar se enfraquece, o espectro passa a ser dominado por características de uma população mais velha e com linhas de emissão típicas de uma galáxia Seyfert. Dentro do cenário acima, a diversidade de populações encontrada pode então ser explicada como sendo devida a diferentes estágios evolutivos da interação; além disso: um outro fator importante parece ser a quantidade de gás disponível; se esta quantidade for muito pequena, pode não ocorrer o disparo de formação estelar. No caso das rádio-galáxias, as interações que teriam originado a atividade nuclear parecem ter ocorrido há mais tempo (1 bilhão de anos atrás para a maioria dos objetos estudados), sugerindo um maior intervalo de tempo entre a interação e o disparo da atividade rádio do que entre a interação e a atividade Seyfert.

Desenvolvimento de produtos em subsidiárias de empresas multinacionais no Brasil

Esta tese é um estudo sobre os papéis das subsidiárias localizadas no Brasil na estratégia de desenvolvimento de produto de empresas multinacionais. O principal objetivo da tese é identificar papéis estratégicos assim como desenvolver e validar um modelo para analisá-los. Esse modelo está alicerçado em um fundamento teórico composto de duas abordagens complementares – a organizacional, que trata das relações hierárquicas e de dependência entre as unidades de uma empresa, e a econômica, que se refere aos mecanismos que visam atingir uma maior eficiência por meio da redução de custos. O modelo teórico é desdobrado em um modelo conceitual de análise de papéis estratégicos que considera os diferentes modos de coordenação das atividades de desenvolvimento de produto pela multinacional, a posição da subsidiária focal em custos de desenvolvimento de produto e os diferentes tipos de interação com parceiros locais. O “modo de coordenação” é de fundamental importância neste estudo e diz respeito às forças que influenciam a distribuição das atividades inovadoras pelas diferentes unidades da empresa multinacional. O método usado para o levantamento de dados foi o de pesquisa survey pela Internet com uma amostra de 146 unidades de desenvolvimento de produtos. A análise de dados foi feita combinando duas técnicas complementares: primeiro, foi usada uma técnica de modelagem de equações estruturais (PLS) para validar o modelo e identificar a estrutura de relações entre os principais construtos; em seguida, foi usada a técnica de análise de conglomerados, a qual fornece mais riqueza de detalhes porque possibilita identificar papéis estratégicos com baixa representação na amostra. Os resultados da análise dão suporte ao modelo, mas contestam-no parcialmente, o que pode ser explicado pela situação particular de um país em desenvolvimento. Ademais, foram identificados cinco grupos significativamente distintos de papéis estratégicos, os “adaptadores locais”, “os inovadores nascentes”, os “inovadores locais”, os “inovadores para mercados emergentes” e os “inovadores globais”. A caracterização desses papéis estratégicos permite tirar conclusões sobre o grau de integração das subsidiárias nas redes globais de inovação das multinacionais, assim como no sistema de inovação do país hospedeiro. A tipologia e o modelo podem servir como base para o desenho de estratégias e políticas de C&T que visem aumentar a integração global e local das atividades de desenvolvimento de produto realizadas em subsidiárias de empresas multinacionais.

Confiança e risco percebido em relacionamentos comprador-vendedor organizacional

Os relacionamentos entre compradores e vendedores organizacionais têm como um dos principais fatores de sucesso o estabelecimento de confiança entre as partes. Também presente na maioria dos relacionamentos, o risco percebido é aplicado explicita ou implicitamente em muitas definições de confiança apresentadas na literatura. A relação entre os construtos da confiança e do risco percebido não é clara, no entanto. Em alguns momentos o risco percebido é concebido como um antecedente do estabelecimento da confiança e em outros a confiança impacta sobre o risco percebido. Uma das hipóteses levantadas neste estudo é que os dois construtos se relacionariam de forma recíproca. Desta forma, o objetivo desta pesquisa é desenvolver e testar um modelo de relações causais recíprocas entre o risco percebido e a confiança em um contexto de relacionamentos organizacionais. Para isso, a relação proposta foi analisada, através da modelagem de equações estruturais, a partir da visão dos clientes de empresas que oferecem o serviço de transporte rodoviário de cargas fracionadas. Além de apresentar a relação entre os dois construtos principais do estudo, o modelo ainda avaliou a relação de ambos com o valor percebido e as intenções de recompra. Os resultados permitiram a aceitação do modelo de relações recíprocas entre risco percebido e confiança, sendo que o impacto do risco percebido sobre a confiança foi negativo (quanto mais risco percebido, menos confiança), porém o impacto da confiança sobre o risco percebido foi positivo (quanto mais se confia, mais se assume o risco).

Um guia de política monetária no Brasil durante o regime de metas de inflação

O objetivo deste trabalho é abordar os principais mecanismos utilizados pelo Banco Central do Brasil (BCB) no processo de decisão acerca da condução da política monetária no período durante o regime de metas de inflação. O BCB analisa e reavalia seus modelos econométricos, que buscam captar inter-relações dos principais agregados da economia brasileira, em conjunto com informações externas aos modelos e julgamentos subjetivos de integrantes do conselho. Portanto, inicialmente será discutida a conjuntura histórica que levou o país a adotar o regime de metas de inflação, assim como sua aplicabilidade. Na segunda parte, são feitas as estimações das equações IS e Phillips que mostram que a taxa de juros afeta o hiato do produto com defasagem de um trimestre e que o produto se correlaciona positivamente com a inflação com defasagem de mais um trimestre. Num terceiro momento, são feitas estimações de uma função de reação do BCB, sugerindo que a instituição leva em consideração o desvio da expectativa de inflação do mercado com relação a meta, o comportamento do hiato do produto e a trajetória dos preços livres em seu processo de decisão.

Desilgualdaddes de Harnack

Neste trabalho provamos a Desigualdade de Harnack para solu»c~oes posi- tivas de equa»c~oes diferenciais parciais el¶³pticas de segunda ordem usando as itera»c~oes de JÄurgen Moser. A partir disso, mostramos que estas solu»c~oes s~ao fun»c~oes de HÄolder e estudamos o seu comportamento no in¯nito. Aplicamos este resultado para provar que uma superf¶³cie m¶³nima, gr¶a¯co de uma fun»c~ao de¯nida fora de um cubo em Rn e com derivadas limitadas, aproxima-se de um plano em Rn+1.

Simulação numérica da transferência simultânea de energia e umidade através do solo em um sistema trocador-armazenador de calor

O presente trabalho pesquisa o fenômeno da transferência simultânea de calor e umidade em solos insaturados que envolvem dutos enterrados. Estes dutos e o solo envolvente compõem o chamado sistema trocador-armazenador de calor no solo, que é muitas vezes utilizado como fonte complementar de aquecimento em estufas solares agrícolas. O funcionamento do sistema trocador-armazenador de calor é baseado na energia armazenada no solo durante os horários de máxima insolação, sendo que durante a noite parte desta energia é recuperada. Durante o dia o ar quente do meio interno da estufa solar é bombeado para dentro do feixe de tubos enterrados, que devolvem o ar mais frio na outra extremidade. Por outro lado, durante a noite o ar mais frio do meio interno da estufa é bombeado para dentro dos dutos, os quais efetuam troca térmica com o solo envolvente, que neste horário possui a temperatura mais elevada do sistema. O objetivo geral deste trabalho é resolver o problema transiente periódico e tridimensional da transferência simultânea de calor e umidade em solos não-saturados, que compõem o sistema trocador-armazenador de calor, utilizando a simulação numérica. Como objetivos secundários têm-se: o melhoramento do sistema de troca térmica, a quantificação da parcela de calor transportado pela difusão da umidade no solo e a análise dos campos de temperatura e de conteúdo de umidade, sendo que a análise dos campos de umidade permite verificar se existe formação de frentes de secagem significativas na vizinhança dos dutos durante os sucessivos períodos de aquecimento e resfriamento a que o sistema é submetido. Na resolução do problema em questão é empregado o modelo clássico de Philip e De Vries para a transferência simultânea de calor e massa em meios porosos insaturados Neste modelo, as equações de conservação de energia e massa obtidas trazem explicitamente as influências combinadas dos gradientes de temperatura e de conteúdo de umidade nos processos de transporte de calor e umidade. O sistema de equações diferenciais governantes do problema em questão é resolvido numericamente utilizando o Método dos Volumes Finitos e na discretização destas equações é usada uma integração temporal totalmente implícita. Todas as propriedades difusivas e termofísicas empregadas são consideradas variáveis com a temperatura e o conteúdo de umidade. Os dutos de seção transversal circular do sistema trocadorarmazenador de calor no solo, são modelados como dutos de seção transversal quadrada de área equivalente para que coordenadas cartesianas possam ser utilizadas nos modelos analisados. Neste trabalho são simulados quatro modelos computacionais associados ao sistema trocador-armazenador de calor no solo. Estes modelos são compostos por: um duto isolado, um duto com convecção, dois dutos isolados e dois dutos com convecção. A variação da temperatura do ar na entrada do(s) escoamento(s), assim como a temperatura do meio ambiente, para os modelos com convecção, é dada por uma senóide com uma amplitude de 14 ºC. No modelo de um duto isolado, são realizadas simulações utilizando várias combinações dos parâmetros do modelo em questão e os resultados, assim obtidos, são comparados com aqueles encontrados na literatura Visando melhorar o sistema de troca térmica dos modelos computacionais investigados, são selecionados valores e intervalos de valores recomendados para os parâmetros do modelo de um duto isolado. Para este modelo, com um diâmetro de 0,1 m, são escolhidos valores (ou intervalos de valores) recomendados: de 4 m/s para a velocidade do escoamento interno dentro do duto, de 0,25 para o conteúdo de umidade do solo, de 5 até 20 metros para o comprimento do duto e de 0,20 até 0,30 m para a distância entre centros do dutos. As simulações dos quatro modelos computacionais realizadas utilizando as várias combinações dos valores recomendados para os parâmetros destes modelos, mostrou que não há diferença significativa entre os valores de calor volumétrico armazenado no solo empregando a resolução acoplada das equações de energia e de massa e a resolução da equação da temperatura. Mesmo para os modelos de um e de dois dutos com convecção a diferença percentual encontrada foi insignificante. Finalmente, são apresentados e analisados os campos de temperatura e de conteúdo de umidade para os quatro modelos computacionais avaliados. Os perfis de temperatura e de conteúdo de umidade em diferentes horários mostraram que, durante o dia, o solo absorve calor dos escoamentos internos de ar e, uma vez que, junto à superfície dos dutos tem-se regiões de maior temperatura, há, conseqüentemente, uma migração da umidade nestas regiões. Durante a noite, ocorre o contrário, o solo fornece calor aos escoamentos dentro dos dutos, e, desta forma, as regiões próximas aos dutos apresentam níveis de conteúdo de umidade superiores ao inicial. Ainda, os perfis de conteúdo de umidade para todas as situações analisadas mostraram que, não há formação de frentes de secagem significativas nas proximidades dos dutos que compõem os quatro modelos computacionais avaliados.

Controle de temperatura em reatores de polimerização em fase gasosa

Nos reatores de polimerização em fase gasosa, o controle de temperatura em uma estreita faixa de temperatura é importantíssimo para assegurar que a temperatura na zona de reação seja mantida acima do ponto de orvalho dos reagentes, assim como abaixo do ponto de fusão do polímero. Outro ponto importante é que se a temperatura estiver em malha aberta, estes reatores estão propensos à instabilidade e ao surgimento de ciclos limites, entretanto podem ser facilmente estabilizados com um controlador de temperatura apropriado. Neste trabalho é realizado um estudo detalhado do comportamento dinâmico de um reator industrial de polimerização em leito fluidizado, com a localização de pontos de bifurcação e estabilização do sistema com controlador PID projetado via otimização no domínio da freqüência. O modelo completo do sistema foi implementado no software MATLAB® para ajustar seus parâmetros com os dados da planta industrial, realizar simulações dinâmicas e analisar as estruturas de controle. Comparações entre as simulações e os dados da planta industrial mostram que as predições do modelo proposto são satisfatórias para reproduzir o comportamento da planta industrial e as propriedades do polímero. A adição de um sistema de recirculação de gás e de um sistema de troca térmica ao modelo do reator torna possível o estudo dos comportamentos dinâmicos e a localização dos pontos de bifurcação, os quais não poderiam ser reproduzidos somente com o modelo do reator. O modelo desenvolvido foi implementado no AUTO®, um software para problemas de continuação e bifurcação em equações diferenciais ordinárias.