Simulação e controle de um sistema de suspensão simplificado

As aplicações da mecânica vibratória vêm crescendo significativamente na análise de sistemas de suspensões e estruturas de veículos, dentre outras. Desta forma, o presente trabalho desenvolve técnicas para a simulação e o controle de uma suspensão de automóvel utilizando modelos dinâmicos com um, dois e três graus de liberdade. Na obtenção das equações do movimento para o sistema massa-mola-amortecedor, o modelo matemático utilizado tem como base a equação de Lagrange e a segunda lei de Newton, com condições iniciais apropriadas. A solução numérica destas equações é obtida através do método de Runge-Kutta de 4ª ordem, utilizando o software MATLAB. Para controlar as vibrações do sistema utilizou-se três métodos diferentes de controle: clássico, LQR e alocação de pólos. O sistema assim obtido satisfaz as condições de estabilidade e de desempenho e é factível para aplicações práticas, pois os resultados obtidos comparam adequadamente com dados analíticos, numéricos ou experimentais encontrados na literatura, indicando que técnicas de controle como o clássico podem ser simples e eficientes.

Critérios para locação amostral baseados em simulação estocástica

Em planejamento de lavra, freqüentemente as campanhas de sondagem são complementadas por falta de segurança nas estimativas. Na maioria dos casos, essas sondagens adicionais são locadas sem o uso de qualquer procedimento matemático que indique os locais mais carentes de informação. Essa dissertação propõe uma metodologia fundamentada em simulação estocástica para tratar esse problema. Com a utilização desse método, é possível que seja considerada a variabilidade entre as informações, além da configuração espacial apresentada. Essa metodologia permite gerar várias estimativas equiprováveis do atributo em questão, reproduzindo o histograma e o variograma das informações originais. A variabilidade apresentada pelos conjuntos de n valores equiprováveis, em cada bloco discretizado da jazida, é medida por meio do cálculo de índices de incerteza. Esses índices são calculados após o processo de simulação estocástica por meio da compilação das n estimativas realizadas para cada bloco. Os índices de incerteza considerados nessa dissertação são: (i) variância condicional; (ii) coeficiente de variação condicional e (iii) intervalo entre quartis. A partir da classificação desses índices de incerteza, são identificados os setores que apresentam as maiores variabilidades em suas estimativas. Com essa classificação é possível a confecção de mapas de incerteza que auxiliam na visualização dos setores/blocos mais variáveis da jazida Nesses setores, são adicionadas informações oriundas da escolha aleatória de um dos n valores equiprováveis disponíveis para cada bloco simulado. Por meio dos índices de incerteza e de acordo com a metodologia, são propostas informações adicionais na área considerada, que contribuem na diminuição das incertezas globais e locais. Essa metodologia usa técnicas extremamente robustas, com embasamento teórico muito consistente, atacando diretamente os setores/blocos que apresentam as maiores flutuabilidades nas estimativas. Um estudo de caso em uma jazida de carvão verifica a aplicabilidade do método. Esse estudo, considera a variabilidade nas estimativas de espessura de uma camada carvão.

Interligação entre o planejamento de longo, médio e curto prazo com o uso do pacote computacional MSPROJECT

A divisão do processo de planejamento e controle da produção (PCP) em diferentes níveis hierárquicos é um dos mecanismos básicos para lidar com a incerteza. Por outro lado, a falta de integração entre estes níveis é apontada como uma dos obstáculos mais importantes à eficácia deste processo, uma vez que dificulta que as decisões tomadas nos diferentes níveis gerenciais sejam consistentes entre si. Este problema é particularmente importante em empreendimentos de construção, os quais, em função de sua complexidade e da incerteza envolvida, em geral requerem freqüentes modificações nos planos, podendo estas ser consideravelmente trabalhosas. O objetivo geral do presente trabalho consiste em propor diretrizes para a interligação do processo de planejamento e controle nos níveis de longo, médio e curto prazo, através da utilização de pacotes computacionais tradicionais, assim como contribuir para a consolidação e refinamento do modelo de PCP proposto por Formoso et al. (1999). O trabalho realizado foi dividido em três grandes etapas. A primeira delas compreendeu a realização de uma pesquisa bibliográfica. A segunda etapa consistiu no desenvolvimento de dois estudos de caso, realizados em empresas de construção do Estado da Bahia, uma delas atuante nos mercados de obras públicas e a outra em incorporações imobiliárias para a classe média-alta. Por fim, na terceira etapa foram analisados os resultados e propostas as diretrizes para a interligação dos diferentes níveis do PCP Entre as conclusões, pode-se destacar que a eficácia dos planos produzidos em pacotes computacionais depende fortemente da existência de um plano de ataque da obra, o qual deve ser discutido com a gerência da produção e explicitado graficamente, de forma a aumentar a transparência de processos. Quando o plano de longo prazo não é elaborado adequadamente, levando em conta as reais necessidades da gerência de produção, o esforço despendido na realização de alterações tende a ser muito grande, dificultando a utilização deste plano para controle em tempo real. O estudo também sugere um conjunto de indicadores que permitem avaliar o cumprimento de prazo e a aderência entre planos de diferentes níveis gerenciais. Estes indicadores são obtidos com relativa facilidade, a partir da utilização do pacote computacional MSProject. Finalmente, o estudo oferece algumas contribuições relativas a critérios para a definição de pacotes de trabalho, à medida que os planos necessitam serem detalhados nos níveis de médio e curto prazo.

Modelagem e simulação de um reator trickle-bed para hidrogenação de propeno

Neste trabalho foi modelado um reator catalítico heterogêneo, pressurizado, adiabático e provido de reciclo. O leito deste reator está dividido em duas partes, sendo que a primeira opera como um reator trickle-bed, enquanto a segunda opera como um reator de leito fixo, onde os regentes estão em fase líquida. Um Sistema de Hidrogenação de Propeno em operação no Pólo Petroquímico do Sul foi utilizado para fornecer os dados necessários à validação do modelo proposto. Neste sistema, a reação de hidrogenação de propeno, em fase líquida, ocorre sobre um catalisador níquel/óxido de níquel. Uma equação, que representa a taxa da reação no intervalo de temperatura de interesse para o sistema, foi desenvolvida utilizando-se dados experimentais obtidos em um reator batelada. O comportamento termodinâmico foi representado através de dados de equilíbrio do sistema binário hidrogênio-propano. Considerou-se que a fase líquida escoa em plug-flow e a fase gasosa, devido à existência de um selo líquido no interior do reator, foi considerada estagnada.

Otimização de pré-formas e matrizes em problemas bidimensionais de forjamento

Este trabalho apresenta uma sistemática para realizar a otimização numérica de pré-formas e de matrizes em problemas de forjamento axissimétricos e em estado plano de deformações. Para este fim, desenvolveu-se um código computacional composto basicamente de três módulos: módulo de pré-processamento, módulo de análise e módulo de otimização. Cada um destes foi elaborado acrescentando rotinas em programas comerciais ou acadêmicos disponíveis no GMAp e no CEMACOM. Um programa gerenciador foi desenvolvido para controlar os módulos citados no processo de otimização. A abordagem proposta apresenta uma nova função objetivo a minimizar, a qual está baseada em uma operação booleana XOR (exclusive or) sobre os dois polígonos planos que representam a geometria desejada para o componente e a obtida na simulação, respectivamente. Esta abordagem visa eliminar possíveis problemas geométricos associados com as funções objetivo comumente utilizadas em pesquisas correlatas. O trabalho emprega análise de sensibilidade numérica, via método das diferenças finitas. As dificuldades associadas a esta técnica são estudadas e dois pontos são identificados como limitadores da abordagem para problemas de conformação mecânica (grandes deformações elastoplásticas com contato friccional): baixa eficiência e contaminação dos gradientes na presença de remalhamentos. Um novo procedimento de diferenças finitas é desenvolvido, o qual elimina as dificuldades citadas, possibilitando a sua aplicação em problemas quaisquer, com características competitivas com as da abordagem analítica Malhas não estruturadas são tratadas mediante suavizações Laplacianas, mantendo as suas topologias. No caso de otimização de pré-formas, o contorno do componente a otimizar é parametrizado por B-Splines cujos pontos de controle são adotados como variáveis de projeto. Por outro lado, no caso de otimização de matrizes, a parametrização é realizada em termos de segmentos de reta e arcos de circunferências. As variáveis de projeto adotadas são, então, as coordenadas das extremidades das retas, os raios e centros dos arcos, etc. A sistemática é fechada pela aplicação dos algoritmos de programação matemática de Krister Svanberg (Método das Assíntotas Móveis Globalmente Convergente) e de Klaus Schittkowski (Programação Quadrática Sequencial – NLPQLP). Resultados numéricos são apresentados mostrando a evolução das implementações adotadas e o ganho de eficiência obtido.

Análise numérica de uma escavação de grande porte em Porto Alegre/RS: caso de obra

Neste trabalho, apresenta-se um estudo numérico, através do Método dos Elementos Finitos, do efeito de instalação de uma parede diafragma com tirantes, seguida de uma grande escavação em solo residual de Granito Independência. As simulações numéricas são comparadas a deslocamentos medidos em uma obra localizada no Bairro Moinhos de Vento, na cidade de Porto Alegre, Rio Grande do Sul - Brasil. Para realizar esta análise utilizou-se o programa computacional Plaxis. Este programa se destina especificamente a análises da deformabilidade e estabilidade de projetos de engenharia geotécnica, cuja simulação requer o uso de métodos numéricos que consideram linearidade e não-linearidade constitutiva, bem como dependência do tempo. Foi utilizado nas simulações o Modelo Elástico-Perfeitamente Plástico com Superfície de Ruptura de Mohr- Coulomb. Os parâmetros do solo foram determinados a partir de ensaios de SPT e CPT, e balizados através de dados obtidos na literatura. Os resultados das análises numéricas são comparados com os resultados medidos in situ, durante a execução da escavação. Foram feitas duas análises distintas, uma chamada Classe A, onde os parâmetros estimados através dos ensaios SPT e CPT foram adotados como padrão e utilizados na simulação da obra Na segunda análise, denominada análise Classe C, correspondente a uma avaliação dos valores de ângulo de atrito (φ), coesão (c) e Módulo de Elasticidade (E), procurando-se avaliar a sensibilidade das previsões. Os resultados da análise numérica Classe A apresentam uma boa aproximação dos resultados medidos em 3 dos 5 perfis analisados, sendo que as simulações reproduziram de forma qualitativa a tendência das curvas experimentais. Parâmetros de resistência e deformabilidade têm um grande influência no comportamento do modelo numérico. A variação de parâmetros de projeto, como a força de ancoragem nos tirantes e a rigidez da parede diafragma também foram testadas. Os resultados demonstram que a influência destes parâmetros é significativa na medida dos deslocamentos.

Ferramenta computacional para a medição de campos de velocidade utilizando processamento digital de imagens

A proposta deste trabalho, consiste na elaboração de uma ferramenta computacional para a medição de campos de velocidades em escoamentos com baixas velocidades (< 0,5 m/s) utilizando o processamento digital de imagens. Ao longo dos anos, inúmeras técnicas foram desenvolvidas com este objetivo. Para cada tipo de aplicação, uma técnica se aplica com maior ou menor eficiência do que outras. Para o caso de estudos em fluídos transparentes, onde o escoamento pode ser visualizado, técnicas que utilizam processamento digital de imagens vêm ganhando um grande impulso tecnológico nos últimos anos. Este impulso, é devido a fatores como: câmaras vídeo filmadoras de última geração, dispositivos de aquisição de imagens e componentes de processamento e armazenamento de dados cada vez mais poderosos. Neste contexto, está a velocimetria por processamento de imagens de partículas cuja sigla é PIV (particle image velocimetry). Existem várias formas de se implementar um sistema do tipo PIV. As variantes dependem, basicamente, do equipamento utilizado. Para sua implementação é necessário, inicialmente, um sistema de iluminação que incide em partículas traçadoras adicionadas ao fluido em estudo. Após, as partículas em movimento são filmadas ou fotografadas e suas imagens adquiridas por um computador através de dispositivos de captura de imagens. As imagens das partículas são então processadas, para a obtenção dos vetores velocidade. Existem diferentes formas de processamento para a obtenção das velocidades. Para o trabalho em questão, devido às características dos equipamentos disponíveis, optou-se por uma metodologia de determinação da trajetória de partículas individuais, que, apesar de limitada em termos de módulo de velocidade, pode ser aplicada a muitos escoamentos reais sob condições controladas Para validar a ferramenta computacional desenvolvida, imagens ideais de partículas foram simuladas como se estivessem em escoamento, através do deslocamento conhecido de vários pixels. Seguindo o objetivo de validação, foi utilizada ainda uma imagem real de partículas, obtida com o auxílio de um plano de iluminação de luz coerente (LASER) e câmaras de vídeo tipo CCD. O programa desenvolvido foi aplicado em situações de escoamento real e os resultados obtidos foram satisfatórios dentro da escala de velocidades inicialmente presumida.

Simulação da irrigação por inundação e da drenagem nos solos de várzea do Rio Grande do Sul

Os solos de várzea são responsáveis por uma parcela significativa da economia do Rio Grande do Sul. As principais atividades de exploração desses solos são o cultivo do arroz irrigado por inundação e a pecuária. Tendo em vista o método de irrigação adotado, o cultivo do arroz é acompanhado de grandes consumos de água e energia. As más condições de drenagem, associadas à uma série de outros fatores, praticamente limitam a atividade agrícola à essa cultura, e fazem com que a pecuária seja quase uma imposição. Tendo esse quadro como referência foram definidos os objetivos desse trabalho, através dos quais buscaram-se algumas respostas que possam auxiliar na definição de um modelo de exploração mais adequado ao grande potencial que esses solos apresentam. Inicialmente o modelo agroidrológico SWAP foi adaptado às particularidades da irrigação por inundação e dos sistemas de drenagem que têm seus drenos implantados com subsoladores. Criou-se então uma nova versão, denominada SWAP_INU. Na seqüência o modelo foi validado a partir da comparação de seus resultados com observações feitas em áreas da Estação Experimental do Arroz do IRGA. Nessas áreas os perfis de solo têm como característica a presença de horizonte B superficial e impermeável. Uma série de simulações permitiram que fossem confrontadas as lâminas de irrigação e profundidades do lençol freático registradas nos experimentos e calculadas pelo modelo. No caso específico da irrigação, as lâminas observadas e calculadas também foram comparadas com as recomendadas pelo IRGA, quando pode-se constatar que o limite inferior do intervalo recomendado é compatível com as observações e cálculos feitos. Além disso, verificou-se que nesse tipo de perfil as características do horizonte A têm pouca importância sobre os resultados finais. Em outra série de simulações procurou-se relacionar a drenagem com a possibilidade de mecanização. Foram criados cenários com diferentes condições de drenagem superficial e subsuperficial, além da presença de drenos livres. Os resultados mostraram uma dependência muito grande entre os sistemas convencionais, e um desempenho bastante satisfatório da drenagem não convencional. Por fim, os resultados de um modelo simplificado, IRRIGA, foram confrontados com os do modelo SWAP_INU. A proximidade entre os valores calculados indicou um ótimo desempenho do modelo simplificado na estimativa da necessidade total de água para a irrigação por inundação.

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.

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.

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.

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.

Dimensionamento de equipes baseado em modelos de previsão, simulação e alocação : caso de uma empresa do setor elétrico

O bom dimensionamento de equipes contribui para o aumento do nível dos serviços prestados pelas empresas, com o menor custo possível. Uma alternativa para abordar a questão foi dimensionar as equipes de eletricistas, de uma empresa do setor elétrico, (utilizando técnicas de previsão de demanda, de simulação e de alocação) para atender de forma otimizada, a demanda variável das atividades prestadas - fornecimento de energia. Um equilíbrio entre a demanda por serviços e a capacidade de execução da empresa evitaria longas filas de espera dos clientes e servidores (eletricistas) ociosos. Cinco etapas forma cumpridas: fase exploratória, coleta de dados, previsão de demanda e simulação do processo e alocação do recurso. Na primeira houve um entendimento de como chegava o pedido do serviço na empresa até a finalização da ordem de serviço. Na coleta de dados foram levantados aproximadamente 80 tipos diferentes de atividades desenvolvidas pelos eletricistas e classificadas de acordo com a prioridade de urgência, prazos de atendimento dos serviços e afinidade de execução das tarefas. Nesta etapa ainda foram coletados os volumes de serviços gerados e tempos médios de deslocamento e execução das atividades. Na terceira etapa foi utilizado um software de previsão de demanda chamado Forecast Pro, possibilitando a escolha automática do modelo de previsão mais apropriado para a série histórica em estudo. Na quarta etapa, foi utilizado um software de simulação de processos chamado Arena. Desenvolveu-se um modelo do processo real com os respectivos dados de entrada dos serviços, tempos de deslocamento e execução e número de equipes. Na última etapa, utilizando a ferramenta Solver do Excel otimizou-se o número de equipes. Um dos resultados da ação foi obter vários cenários com a variação do número de equipes e seus respectivos tempos médios de atendimento, sem causar nenhum dano para a empresa, podendo assim ser analisado qual o melhor cenário para ser implementado na companhia, minimizando o problema.

Arquitetura para um ambiente de grade computacional direcionado às aplicações distribuídas, móveis e conscientes do contexto da computação pervasiva

Neste início de década, observa-se a transformação das áreas de Computação em Grade (Grid Computing) e Computação Móvel (Mobile Computing) de uma conotação de interesse emergente para outra caracterizada por uma demanda real e qualificada de produtos, serviços e pesquisas. Esta tese tem como pressuposto a identificação de que os problemas hoje abordados isoladamente nas pesquisas relativas às computações em grade, consciente do contexto e móvel, estão presentes quando da disponibilização de uma infra-estrutura de software para o cenário da Computação Pervasiva. Neste sentido, como aspecto central da sua contribuição, propõe uma solução integrada para suporte à Computação Pervasiva, implementada na forma de um middleware que visa criar e gerenciar um ambiente pervasivo, bem como promover a execução, sob este ambiente, das aplicações que expressam a semântica siga-me. Estas aplicações são, por natureza, distribuídas, móveis e adaptativas ao contexto em que seu processamento ocorre, estando disponíveis a partir de qualquer lugar, todo o tempo. O middleware proposto, denominado EXEHDA (Execution Environment for Highly Distributed Applications), é adaptativo ao contexto e baseado em serviços, sendo chamado de ISAMpe o ambiente por este disponibilizado. O EXEHDA faz parte dos esforços de pesquisa do Projeto ISAM (Infra-Estrutura de Suporte às Aplicações Móveis Distribuídas), em andamento na UFRGS. Para atender a elevada flutuação na disponibilidade dos recursos, inerente à Computação Pervasiva, o EXEHDA é estruturado em um núcleo mínimo e em serviços carregados sob demanda. Os principais serviços fornecidos estão organizados em subsistemas que gerenciam: (a) a execução distribuída; (b) a comunicação; (c) o reconhecimento do contexto; (d) a adaptação; (e) o acesso pervasivo aos recursos e serviços; (f) a descoberta e (g) o gerenciamento de recursos No EXEHDA, as condições de contexto são pró-ativamente monitoradas e o suporte à execução deve permitir que tanto a aplicação como ele próprio utilizem essas informações na gerência da adaptação de seus aspectos funcionais e não-funcionais. O mecanismo de adaptação proposto para o EXEHDA emprega uma estratégia colaborativa entre aplicação e ambiente de execução, através da qual é facultado ao programador individualizar políticas de adaptação para reger o comportamento de cada um dos componentes que constituem o software da aplicação. Aplicações tanto do domínio da Computação em Grade, quanto da Computação Pervasiva podem ser programadas e executadas sob gerenciamento do middleware proposto.

Simulação numérica da combustão turbulenta de gás natural em câmara cilíndrica

O presente trabalho apresenta uma modelagem detalhada de processos de combustao turbulentos para um jato concentrico de combustıvel e ar. A modelagem é fundamentada nas equacões de conservacão de massa, de quantidade de movimento, de energia e de espécies quımicas. A turbulencia é resolvida pela utilizacão do modelo k- padrão. Dois modelos de reacões quımicas são apresentados. O modelo SCRS – Simple Chemically-Reacting Systems, que assume taxas instantâneas de reacões quımicas. Também é abordado o modelo E-A – Eddy Breakup - Arrhenius, que assume taxas finitas de reacões quımicas. A radiacão térmica, fenômeno de grande importância devido as altas temperaturas alcancadas em processos de combustão,é modelada através do Método das Zonas. O modelo da soma ponderada de gases cinzas – WSGGM,é usado para determinar o espectro de emissão e absorcão dos gases no processo. Para a solucão destas equacões diferenciais, juntamente com os modelos de turbulência, de reaçõoes químicas e radiação térmica, faz-se o uso do Método dos Volumes Finitos. Para validar a modelagem apresentada resolve-se o processo de combustão em uma câmara cilíndrica. A câmara de combustão usada áa mesma abordada no First Workshop on Aerodynamics of Steady State Combustion Chambers and Furnaces, organizado pela ERCORTAC - European Research Community On Flow Turbulence And Combustion, em outubro de 1994, que apresenta dados experimentais de temperatura e concentração das espécies químicas para várias posições de interesse no interior da câmara. Utiliza-se o gás natural como combustível e o ar atmosférico como oxidante. O processo de combustão sem pré-mistura é resolvido para a condição de excesso de combustível de 5 % para ambos os modelos, onde o gás natural é injetado por um duto circular central, e o ar atmosférico por um orifício anular externo a esse duto, no mesmo plano Uma reação química não estagiada é assumida para o modelo SCRS. Para o modelo E-A duas situações são resolvidas: combustão não estagiada, com uma etapa global de reação química; e reação quımica estagiada, com duas etapas globais. Os resultados obtidos com o modelo SCRS para a distribuição de temperaturas, em termos de tendências gerais, são razoáveis. Já as concentrações de espécies químicas não apresentam dados satisfatórios para este modelo. Para o modelo E-A os resultados apresentam boa concordância com os dados experimentais, principalmente para a situação em que o processo de combustão é assumido em duas etapas globais. ´E analisado em detalhe o papel desempenhado pela transferencia de calor por radiacao, com meio participante. Para melhor verificar as trocas de calor, assume-se uma camara de combustao cilındrica com paredes d’agua. A injecao do combustıvel e do oxidante e feita atraves de um queimador central, semelhante ao usado para validar a modelagem, porem com dois orifıcios concentricos para injecao de combustıvel. Nesta situação o efeito do turbilhonamento (swril), assumido como 20 % da velocidade axial de entrada, sobre a injecao de ar e computado atraves da condicao contorno da equacao de conservacao da quantidade de movimento angular. Nesta fase apenas o modelo E-A, com duas etapas globais de reacoes quımicas, e considerado, ja que o mesmo apresenta os melhores resultados. O processo de combustao e simulado com e sem a presenca da radiacao termica. Verifica-se que a presenca da radiacao termica homogeneiza a temperatura dos gases no interior da camara. Com isso verifica-se tambem alterações nas taxas de reacoes quımicas, modificando a magnitude das fracoes das especies quımicas Quando a radiacao termica e considerada efeitos de extinção local da chama sao verificados nas regioes de temperaturas mais altas, diminuindo o consumo de oxigenio e aumentando a producao de monoxido de carbono, caracterizando assim uma combustao incompleta. Em algumas situacoes tem-se uma variacao de temperatura de ate 500 K, a montante da chama. A radiacao termica tambem aumenta a taxa de transferencia de calor dos gases quentes para as paredes da camara, e desta para o seu exterior. Com os resultados obtidos a partir desta modelagem e possıvel determinar o perfil da zona de combustao, a distribuicao de concentracoes de especies quımicas, o campo de velocidades e as taxas de transferencia de calor para as paredes da camara de combustao, total, por conveccao superficial e por radiacao. Estes resultados sao de extrema importancia para prever a performance de camaras de combustao, assim como auxiliar na sua otimizacao.