Modelagem numérica da evolução da linha de costa das praias localizadas a oeste da cidade de Fortaleza, Ceará : trecho compreendido entre o rio Ceará e a praia do Cumbuco

Há aproximadamente meio século, as praias situadas a sotamar do Porto do Mucuripe, em Fortaleza, vem sofrendo intensos processos erosivos, creditados em grande parte à construção e ampliação deste porto. O fato é que o acentuado crescimento urbano da capital cearense ocasionou a fixação de dunas e a quebra do fluxo longitudinal de sedimentos em seu litoral, resultando no recuo da linha de costa e na necessidade de intervenção antrópica por meio de obras rígidas que viessem a garantir a preservação da infra-estrutura existente nos trechos mais afetados. Como conseqüência da fixação das praias, o suprimento de material sedimentar passou a ficar retido, enquanto que o potencial de transporte das ondas se preservou. A quebra deste equilíbrio dinâmico acarretou a transferência dos processos erosivos para as praias adjacentes, o que tornou-se um problema cada vez maior, pois as soluções adotadas nestas praias eram idênticas às anteriores. As conseqüências deste processo para uma cidade como Fortaleza, onde o turismo é uma das principais fontes de renda, são graves, dado que como resultado final, encontramos longos trechos de praias com a balneabilidade comprometida e perda de qualidade visual. O litoral situado a oeste da capital é limitado à direita pela foz do Rio Ceará e à esquerda por um promontório rochoso, onde situa-se a Ponta do Pecém. Este trecho compreende aproximadamente 30 km de praias arenosas, com granulometria média e fina, e com ondas incidindo sobre a costa de forma obliqua, o que as torna o principal mecanismo de transporte de sedimentos. A ocupação urbana concentra-se principalmente nas praias mais próximas a Fortaleza, onde observa-se ainda, o afloramento de rochas de praia e grande perda de material sedimentar, fornecendo indícios da transferência dos processos erosivos da orla marítima da capital para estas praias. Com a conclusão das obras do Porto do Pecém e de um pólo industrial que visa desfrutar da localização estratégica deste porto, é natural que ocorra uma intensificação nos processos de ocupação urbana das praias próximas à área. Tal constatação motivou um trabalho de modelagem da dinâmica desta zona com o objetivo de nortear um plano de uso e ocupação das áreas localizadas próximas à praia, de forma que se possa prever o comportamento da linha de costa e evitar que sejam repetidos certos equívocos como a construção em zonas de forte dinâmica e a fixação das fontes primárias de fornecimento de sedimentos, que são as dunas frontais. Dada a disponibilidade de dados, bons processadores e aos custos significativamente reduzidos da modelagem numérica, adotou-se o pacote GENESIS – RCPWAVE, que além de ser de domínio público, é a base do sistema de modelagem de linha de costa adotado pelo CERC (Coastal Engineering Research Center), U.S.A., para aplicações em costa aberta, em regiões sujeitas às intervenções humanas. A calibração do modelo se fez considerando as linhas de praia medidas em 1974 pela DHN e em 2001 com o uso de GPS. Os dados de onda utilizados foram obtidos por um ondógrafo direcional do tipo Waverider, instalado a uma profundidade de 18 metros nas proximidades da Ponta do Pecém. Os dados relativos ao modelo conceitual dos processos predominantes na região, como: contribuições externas, variação granulométrica e variações sazonais de perfis foram obtidos de levantamentos bibliográficos de trabalhos anteriores. Por último, informações relativas às estruturas existentes e seu comportamento, ao afloramento de formações rochosas e o último levantamento da linha de praia, foram obtidas através de trabalhos de campo. De uma forma geral, o comportamento previsto pelo modelo mostrou-se semelhante ao observado nos diferentes levantamentos. Considerando-se as limitações dos processos envolvidos no levantamento de dados, onde tanto a carta da DHN quanto o mapeamento por satélite estão sujeitos a imprecisões e ainda, que a série de dados confiáveis de ondas para a região possuía apenas dois anos, é importante notar que, em linhas gerais, a formulação matemática do modelo representou satisfatoriamente os processos envolvidos. Os resultados fornecidos possibilitam a extrapolação da evolução da linha de costa e indicam pontos de provável recuo ou avanço da praia, norteando a sua ocupação. A ferramenta gerada proporciona ainda a avaliação do impacto de intervenções por meio de estruturas rígidas ou engordamento de praia ao longo do tempo e gera uma estimativa dos valores de deriva litorânea para os diferentes trechos de praia, possibilitando avaliar os efeitos das intervenções nas praias adjacentes.

Utilização de imagens de satélite e modelagem numérica para o estudo da dispersão de poluentes nas usinas termoelétricas de Charqueadas e São Jerônimo

O objetivo geral deste trabalho é desenvolver estudos relacionados à dispersão de poluentes, considerando a queima de carvão para geração de energia nas Usinas Termoelétricas de Charqueadas e São Jerônimo. O período de estudo foi do dia 17 a 23 de junho de 2003. Neste período houve a passagem de um sistema frontal. Sendo possível avaliar a dispersão dos poluentes em condições pré-frontal, frontal e pós-frontal. Para simular o comportamento dos poluentes neste período, foi utilizada uma subrotina de dispersão acoplada ao modelo RAMS (Regional Atmospheric Modeling System). Os resultados mostraram que nos dias classificados como pré-frontal e pós-frontal as concentrações do material particulado, dióxido de enxofre e do óxido de nitrogênio, atingiram seus valores máximos, pelo fato da umidade relativa do ar estar bastante baixa em torno de 60%, pressão atmosférica da ordem de 1021 hPa e a intensidade dos ventos fraca. No dia classificado como frontal, as concentrações estavam praticamente nulas, devido à passagem do sistema frontal que causou a queda na pressão atmosférica, aumento da umidade relativa e também pelo fato da ocorrência de precipitação atmosférica neste dia. As comparações dos resultados simulados, com os dados observados na estação de qualidade do ar mostraram-se satisfatórios. Com exceção do dia 22 de junho, que apresentou uma diferença da ordem de 90%.

Desenvolvimento de uma metodologia para planejamento urbano, utilizando técnicas de sensoriamento remoto, modelagem numérica do terreno (MNT) e sistemas de informações geográficas (SIG)

o monitoramento da expansão das áreas urbanas e a análise da sua interação com o meio físico têm sido um grande desafio para os técnicos de planejamento urbano. No Brasil, em especial, dada a velocidade com que o fenômeno se processa e graças a um crescimento desordenado das cidades nas últimas décadas, esses estudos, que envolvem um elevado número de informações, tem exigido decisões e diagnósticos urbanos cada vez mais rápidos. Esta dissertação propõe uma metodologia para o planejamento racional do uso do solo urbano através do emprego integrado de tecnologias recentes como Sistema de Informações Geográficas (SIG), Modelagem Numérica do Terreno (MNT) e Sensoriamento Remoto através de imagens orbitais. Para isso, são implementados no SIG desenvolvido pelo INPE dados provenientes de cartas topográficas, de mapas temáticos do meio físico e de imagens orbitais LANSAT/TM da região estudada. A partir desses dados iniciais são geradas, também num SIG, outras informações com objetivo de estudar a evolução da área urbana, identificar áreas com suscetibilidade preliminar à erosão laminar, áreas com restrição ao uso urbano e áreas de eventos perigosos e riscos. o trabalho apresenta inicialmente uma revisão bibliográfica sobre a aplicação de Sensoriamento Remoto, Modelagem Numérica do Terreno (MNT) e Sistema de Informações Geográficas (SIG) em estudos urbanos. Segue-se a conceituação e aspectos teóricos dessas três ferramentas básicas utilizadas. A metodologia propriamente dita traz os planos de informações originais e as suas respectivas fontes de informações, os processos de classificação de imagens digitais empregados e os modelos de cruzamentos desenvolvidos para um SIG. A área teste escolhida é a sub-bacia do Arroio Feijó, localizada na região metropolitana de Porto Alegre, na porção centro-leste do Estado do Rio Grande do Sul. A região é caracterizada por uma elevada densidade populacional, pela presença de áreas inundáveis e pela ocorrência de processos eroslVOS. Os resultados mostram que a metodologia proposta é adequada e eficiente para agilizar as atividades de planejamento urbano, subsidiando a elaboração de Planos Diretores de Desenvolvimento Integrado e orientando o crescimento das cidades para regiões mais favoráveis. Além disso, contribui para a prevenção de parcela dos riscos e problemas geotécnicos relacionados ao meio físico nas áreas urbanas.

Estudo para a padronização das dimensões de unidades de alvenaria estrutural no Brasil através do uso da coordenação modular

Este estudo apresenta uma proposta para a padronização das dimensões de unidades de alvenaria estrutural no Brasil. Para isso, primeiramente são discutidos os benefícios da aplicação da teoria da coordenação modular na alvenaria estrutural, através da relação da coordenação modular com a arquitetura e da síntese dos seus princípios fundamentais. Tendo sempre o foco na padronização do componente modular e no seu uso em projetos coordenados modularmente. Em seguida, é apresentada a situação da alvenaria estrutural no Brasil e na Alemanha, observando principalmente o uso da coordenação modular, as características das unidades, do projeto e da execução nesse processo construtivo. Destacamse os problemas enfrentados no Brasil, onde adaptações do projeto arquitetônico de edifícios em concreto armado para projetos em alvenaria estrutural aumentam os custos e problemas no processo construtivo. A principal razão para isto é que as unidades de alvenaria não são coordenadas modularmente e necessitam rigorosa padronização. Assim, finalmente, a partir dos aspectos condicionantes para as dimensões das unidades de alvenaria estrutural é possível alcançar o objetivo central deste trabalho que é a sugestão à indústria brasileira de uma série de medidas preferidas para a padronização das dimensões de unidades de alvenaria estrutural no Brasil, apresentando os possíveis arranjos de paredes e elaborando um projeto piloto coordenado modularmente, sem que haja a restrição da liberdade criativa dos arquitetos.

Modelagem numérica aplicada aos fenômenos de escoamento e mistura em modelos físicos de panelas de aciaria

Esse trabalho tem por objetivo geral o desenvolvimento de uma metodologia de modelagem numérica que represente o escoamento e o fenômeno de mistura em um modelo físico de panela de aciaria de base elíptica. Os objetivos específicos do trabalho são: o estudo dos coeficientes das forças de interação líquido-bolha, dos modelos de turbulência e da mudança do formato da base da panela de circular para elíptica. O escoamento e o fenômeno de mistura foram calculados através do método de Volume Finitos baseado em Elementos por meio do software CFX5.7 da Ansys. Dados da literatura e ensaios em modelo físico, realizados em laboratório, auxiliaram na validação dos modelos numéricos. O estudo dos coeficientes das forças de não-arrasto mostrou que os resultados da distribuição de ar ao longo da altura do banho mudam com a variação dos coeficientes. No final, coeficientes para 3 configurações de panelas em diferentes vazões de ar foram encontrados. Com relação ao estudo dos modelos de turbulência, observou-se que para a solução do escoamento e do fenômeno de mistura em uma panela de base circular, o k-ε é o modelo de turbulência mais indicado. Por outro lado, para uma panela de base elíptica, o modelo RSM mostrou-se o mais adequado. Finalmente, com relação ao estudo da mudança do formato da base da panela, observou-se que os tempos de mistura de uma panela de base elíptica são maiores que uma de base circular e aumentam à medida que a vazão de ar diminui.

Análise interativa de situações de reforço e recuperação de peças de concreto armado, empregando o método dos elementos finitos

O objetivo deste trabalho é desenvolver um programa computacional, baseado no método dos elementos finitos, para simular situações de reforço e recuperação de peças de concreto armado. A linguagem Matlab é a ferramenta utilizada para a elaboração do programa. É uma linguagem de alta performance para computação técnica. Integra computação, visualização e programação em um fácil ambiente para uso, onde problemas e soluções são expressos em familiar notação matemática. A característica principal deste programa é de permitir alterações na modelagem numérica durante a análise do problema, podendo-se retirar ou acrescentar elementos da estrutura, aumentar ou diminuir rigidezes, alterar materiais, etc, viabilizando-se, assim, uma avaliação das diversas etapas ou hipóteses dos procedimentos de recuperação ou reforço de estruturas. O programa permite a mudança de vinculações do elemento estrutural, durante a análise. Além disto, é permitido não só incrementos de forças como incrementos de deslocamentos impostos à estrutura. O programa computacional utiliza duas etapas distintas, baseadas em um modelo elasto-viscoplástico. Na primeira etapa, faz-se a determinação da resposta da estrutura ao longo do tempo. Considera-se, nesta etapa, que o material possui comportamento viscoelástico. Na segunda, busca-se determinar a resposta da estrutura para um carregamento instantâneo, considerando-se o material com comportamento elastoplástico Para melhor representar o comportamento do concreto, quando submetido a carregamento de longa duração, utilizou-se um modelo de camadas superpostas. A armadura é introduzida no modelo como uma linha de material mais rígido dentro do elemento de concreto. Considera-se aderência perfeita entre o concreto e o aço. Assim, os deslocamentos ao longo das barras de aço são determinados a partir dos deslocamentos nodais dos elementos de concreto. Para verificar a precisão do programa computacional desenvolvido, comparam-se os resultados numéricos com resultados experimentais, cujos ensaios foram realizados no Laboratório de Ensaios e Modelos Estruturais do Programa de Pós-Graduação em Engenharia Civil da Escola de Engenharia da Universidade Federal do Rio Grande do Sul, e no Laboratório de Estruturas da Escola de Engenharia de São Carlos- USP. Nessas comparações, obtiveram-se ótimas aproximações entre os resultados numéricos e experimentais.

Análise numérica de uma torre de resfriamento de grande porte

Este trabalho tem como objetivo principal determinar a troca de calor em uma torre de resfriamento de grande porte, do tipo seca. A força motriz deste sistema é o empuxo resultante da diferença de temperatura entre o ar interior e exterior da torre. A partir de uma certa velocidade o vento externo também desempenha um papel importante. A convecção natural ocorre em regime turbulento, devido às grandes dimensões. A modelagem numérica é feita pelo Método dos Volumes Finitos com malhas criadas sobre Elementos Finitos, através dos softwares comerciais Fluent e CFX. Os resultados são obtidos para diferentes potências dissipadas no trocador de calor. A partir dos resultados são realizados estudos sobre a viabilidade de mudanças na geometria da torre, a fim de se obter melhorias no escoamento, bem como uma troca de calor mais eficiente. Faz-se ainda uma regressão a fim de obter correlações entre a potência dissipada no trocador de calor e a vazão mássica de ar necessária para a troca. Estas correlaçõoes são de utilidade, para aplicações posteriores, na implantação de um simulador para o processo de troca de calor nesta torre. O trocador de calor também é enfocado a fim de se determinar o coeficiente global de troca térmica, a fim de disponibilizá-lo para aplicação em outros trabalhos

Aplicação de um modelo de dispersão atmosférico na região de Candiota-RS e sua relação com as condições meteorológicas

Este trabalho tem por objetivo, estudar a dispersão de poluentes na região da Usina Termelétrica Presidente Médici (UTPM), localizada no município de Candiota-RS, utilizando a modelagem numérica (mais especificamente o modelo ISCST - Industrial Surce Complex Term Short Term) e relacionar o seu comportamento com as condições meteorológicas. Esse modelo de dispersão utilizado é recomendado pela EPA (Environmental Protection Agency) para tratamento de dispersão de poluentes emitidos por fontes industriais como refinarias, termelétricas, siderúrgicas, etc. O período de estudo foi de 19/01/03 à 28/01/03, quando se realizou a primeira campanha experimental, pertencente ao projeto "Estudo Ambiental Aplicando a Modelagem Numérica na Região de Candiota-RS". Os resultados simulados pelo modelo foram comparados com os dados das concentrações de Material Particulado (MP) e Dióxido de Enxofre (SO 2 ), medidos através dos pontos receptores em duas estações de qualidade do ar próximas à Usina Termoelétrica. Uma correlação entre as concentrações de poluentes modelados e os dados coletados foi realizada, assim como, uma comparação com os padrões de qualidade do ar estipulados CONAMA (Conselho Nacional do Meio Ambiente). Os resultados denotam que o modelo ISCST subestima as concentrações de MP e SO 2 observadas nas estações de qualidade do ar, o que provavelmente ocorre, devido à presença de outras fontes de emissão (como indústrias de cimento e calcário) nas proximidades das estações. A situação com a melhor correlação nos dez dias, entre modelo e amostrador, apresentou-se com 90%. Em apenas um dos dez dias, foi excedido o padrão secundário (100 m g/m 3 ) com o poluente SO 2 , caracterizado a qualidade do ar da região como boa para esse período.

Estudo comparativo sobre a permeabilidade das alvenarias em blocos cerâmicos e alvenarias em blocos de concreto

A umidade constitui-se em uma das principais manifestações patológicas que afligem os usuários de edificações. Procura-se, neste trabalho, abordar as diversas características dos materiais, bem como, os fenômenos e propriedades físicas que agem no sentido de propiciar a penetração de umidade nas alvenarias. Neste trabalho, analisa-se o desempenho de alvenarias compostas por unidades de blocos cerâmicos e blocos de concreto, quanto à penetração de umidade. São apresentados métodos de ensaios concebidos para avaliação da estanqueidade em paredes e, dentre estes, três métodos de ensaio usualmente utilizados pelo meio técnico são empregados na análise da capacidade estanque, no sentido de se proceder a um estudo comparativo do desempenho das alvenarias de blocos cerâmicos e de concreto, com e sem revestimento.

Situações atmosféricas que propiciam as marés meteorológicas

O objetivo central deste trabalho consiste, a partir de um caso de ciclogênese ocorrido na costa leste do América do Sul, caracterizar as condições sinóticas favoráveis ao desenvolvimento de ciclones marítimos que possam gerar ressacas. Uma ciclogênese desenvolveu-se a leste da Argentina dia 03 de maio de 2001. Neste mesmo dia sobre o continente, existia um centro de alta pressão, que em conjunto com este ciclone marítimo formava um pista de vento de quadrante sul, paralelo ao continente, deixando o mar revolto, aumentando nível no extremo sul do Brasil, como registrado pela tábua de maré da cidade de Rio Grande, no estado do Rio Grande de Sul. No dia 04 de maio, à leste do estado de São Paulo, outro ciclone começa a se formar, intensificando a pista de vento pré existente. A metodologia utilizada afim de atingir os objetivos propostos, foi a da modelagem numérica. Utilizou-se modelo Regional Atmospheric Modeling System (RAMS). Foram utilizados como entrada dados de outros dois modelos, modelo Global e Modelo ETA. Para ambos os experimentos usando o RAMS, observou-se que o ciclone do dia 04 teve sua gênese em torno de 25ºS e 45ºW, se deslocando para sudeste. Na análise de mesoescala, foi observado que durante o período estudado os ventos sempre tiveram uma componente de sul e que os mais intensos ventos ocorreram em 28ºS e 48ºW, pois nesta região haviam maiores valores de gradiente de pressão, fluxo de calor sensível e fluxo de calor latente, já que é também nesta região onde se verifica a maior temperatura da superfície do mar.

Solução de um problema não linear do tipo reação-difusão na modelagem dispersão de insetos

Neste trabalho tratamos da solução de um problema não linear do tipo tração-difusão, na modelagem de dispersão de insetos. Começamos estabelecendo uma lei de conservação e a partir desta, deduzimos algumas equações importantes para o desenvolvimento do nosso estudo, tais como a equação de convecção, de difusão e simultaneamente convecção e difusão. Se considerarmos uma escala de tempo que possibilite a adição ou retirada de indivíduos no meio, conforme seja considerada reprodução, migração ou morte, podemos acrescentar ao processo difusivo um termo de reação, obtendo então, a equação do tipo reação-difusão. Se o temp de reação for deendendee da densidade populacional e do tipo logístico, obtém-se a equação de Fischer. Dessa equação abordamos alguns aspectos, tais como, determinação dos estados estacionários, análise da estabillidade dos mesmos, representação gráfica no plano de fase e por último investigamos a existência de solução do tipo onda viajante. Abordamos, também, alguns exemplos apresentados na literatura, envolvendo equação da difusão com coeficiente constante e com coeficiente dependente da densidade populacional. Além disso, apresentamos o resultados obtidos com a modelagem em tempo discreto, a partir de um trabalho experimental com besouros marcados para o experimento e depois liberados Banks et al (1985) , em que os autores admitiram uma variação temporal e a partir dos dados obtidos fizeram uma estimativa para o coeficiente de difusão D (t), bem como para o coeficiente de decaimento α(t) do termo de reação linear em u. Construimos curvas que se ajustam a essas alternativas e apresentamos esses coeficientes em versão continua D (t) e α(t), dependentes da variável tempo t. Através de uma abordagem numérica, os modelos foram comparadas da variável tempo t. Através de uma abordagem numérica, os modelos foram comparados para diversos casos, usando diferentes combinação de D constante e D variando no tempo, a constante e a variando no tempo. Além disso, analisamos tambén, o efeito da substituição do coeficiente de difusão D constante por D(t) na equação de Fisher.

Análise numérica de problemas de fratura elástica e elastoplástica

Este trabalho apresenta o método dos elementos finitos em conjunto com métodos numéricos especificos para a solução de problemas de fratura. Esta é uma poderosa ferramenta para a análise de fraturas e soluções confiáveis são obtidas para problemas complexos de Engenharia tanto no campo linear como no não-linear. O elemento finito. implementado é do tipo isoparamétrico quadrâtico da família Serendipity. Com dois graus de liberdade por nó, permite discretizar em estado plano de tensão ou deformação estruturas com geometrias bastante variadas. Para a análise linear são implementadas quatro técnicas consagradas para a avaliação do fator de intensidade de tensão no modo I de fratura: extrapolação de doslocamentos (usando malha somente com elementos convencionais e malha mesclada com elementos especiais), taxa de liberação de energia de defermação, extensão virtual da trinca e o método da integral J, descartando-se neste caso a hipótese de descarregamento. A linguagem de programação adotada é o FORTRAN 77. A entrada de dados é feita por intermédio de arquivos previamente preparados. Os resultados obtidos são confrontados com resultados experimentais e computacionais fornecidos por outros programas. Analisam-se placas, estruturas de uso na indústria e simulam-se ensaios como o corpo de prova de flexão em três pontos e o corpo de prova de tensão. compacto.

Modelo numérico para análise à flexão de elementos estruturais com protensão aderente e não aderente

Na protensão não aderente, a armadura permanece livre para se movimentar ao longo de seu perfil em todas as seções com exceção das de ancoragem. Não há aderência entre concreto e armadura, e a hipótese da compatibilidade de deformações entre o aço e concreto não é aplicável, tornando inviável o desenvolvimento de uma solução analítica. Visando colaborar para a maior compreensão do comportamento à flexão das estruturas com protensão não aderente e para o desenvolvimento de critérios nacionais de projeto, um modelo numérico foi implementado utilizando o elemento finito do tipo híbrido para pórticos planos. Nesta formulação, a equação para as solicitações ao longo do elemento é a função de interpolação, e as forças as variáveis interpoladas. Como esta função resulta das condições de equilíbrio, sem hipóteses arbitrárias, o método é considerado exato para forças e curvaturas. Elementos longos são possíveis, de maneira que um único elemento finito pode ser utilizado para um vão de viga ou pilar, reduzindo o esforço computacional. O caráter exato da formulação contribui para a boa modelagem dos cabos não aderentes, já que a tensão nestas armaduras depende das curvaturas de todas as seções do elemento. O modelo numérico proposto prevê a não linearidade geométrica, carregamentos cíclicos e a construção composta. Relações constitutivas já consolidadas na literatura são empregadas para os materiais. A cadeia de Maxwell é utilizada para representar o comportamento reológico do concreto e do aço de protensão, respeitando as características de envelhecimento de cada material. Inúmeros exemplos são apresentados, permitindo a comparação entre resultados numéricos e experimentais Uma análise paramétrica foi realizada, caracterizando o desempenho do modelo numérico frente a variações nos parâmetros de entrada. Discutem-se, ainda, dois critérios de ruptura para a utilização do modelo numérico no estudo do comportamento à flexão das estruturas com protensão não aderente. Os resultados indicam o bom desempenho do modelo numérico e a sua adequação para a realização de pesquisas sobre o assunto.

Desenvolvimento de um plano agregado de produção para um sistema agroindustrial

Esta dissertação traz como tema principal o planejamento da produção, com enfoque no planejamento agregado da produção. Tal planejamento de longo prazo é o elo entre a alta administração da empresa e a manufatura. Através dele, as metas estratégicas da empresa são comunicadas por meio de uma linguagem comum ao corpo tático/operacional, visando coordenar os esforços e buscar comprometimento para o alcance das metas. O desenvolvimento de um plano agregado de produção é concretizado através do uso de algoritmos de pesquisa operacional, visando otimizar o sistema, considerando-se os fatores restritivos. O planejamento agregado da produção definirá o grau de utilização dos recursos de produção e servirá de base para elaboração dos planejamentos de produção mais detalhados, que definirão como melhor utilizar os recursos disponibilizados num horizonte de planejamento mais curto e sob menores incertezas de mercado. Esta dissertação traz como aplicação prática a modelagem de um plano agregado de produção para uma agroindústria de grande porte e diversas plantas produtoras, no segmento de produtos suínos.

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.