Desprendimento de vórtices e controle em esteira de cilindros por simulação numérica direta

O desprendimento de vórtices é responsável por vibrações prejudiciais a estruturas e pelo aumento de seus coeficientes de arraste. Portanto o correto entendimento e o controle do desprendimento de vórtices é de grande interesse de modo a diminuir custos de manutenção de estruturas estáticas e custos com combustíveis de estruturas móveis, por exemplo. O caráter tridimensional do desprendimento de vórtices a baixos números de Reynolds, entre 100 e 300, é aqui estudado através de Simulações Numéricas Diretas (DNS) do escoamento ao redor de cilindros circulares em duas e três dimensões. O controle do desprendimento de vórtices foi estudado nesta dissertação, através do uso de placas de separação do escoamento em simulações numéricas bidimensionais para números de Reynolds 100, 160 e 300, mostrando que para estes números de Reynolds é possível reduzir as freqüências de desprendimento, e em alguns casos até suprimi-las. A simulação tridimensional para Re = 300 mostrou um escoamento aparentemente bidimensional que afeta sensivelmente a freqüência de desprendimento quando em comparação o caso 2D. Para a realização das simulação numéricas foi utilizado um esquema compacto de diferenças finitas em associação com um método de representação virtual de obstáculos imersos em escoamentos, que permite a representação de geometrias complexas sobre uma malha cartesiana.

Resolução de sistemas de equações lineares através de métodos de decomposição de domínio

A paralelização de métodos de resolução de sistemas de equações lineares e não lineares é uma atividade que tem concentrado várias pesquisas nos últimos anos. Isto porque, os sistemas de equações estão presentes em diversos problemas da computação cientí ca, especialmente naqueles que empregam equações diferenciais parciais (EDPs) que modelam fenômenos físicos, e que precisam ser discretizadas para serem tratadas computacionalmente. O processo de discretização resulta em sistemas de equações que necessitam ser resolvidos a cada passo de tempo. Em geral, esses sistemas têm como características a esparsidade e um grande número de incógnitas. Devido ao porte desses sistemas é necessária uma grande quantidade de memória e velocidade de processamento, sendo adequado o uso de computação de alto desempenho na obtenção da solução dos mesmos. Dentro desse contexto, é feito neste trabalho um estudo sobre o uso de métodos de decomposição de domínio na resolução de sistemas de equações em paralelo. Esses métodos baseiam-se no particionamento do domínio computacional em subdomínios, de modo que a solução global do problema é obtida pela combinação apropriada das soluções de cada subdomínio. Uma vez que diferentes subdomínios podem ser tratados independentemente, tais métodos são atrativos para ambientes paralelos. Mais especi camente, foram implementados e analisados neste trabalho, três diferentes métodos de decomposição de domínio. Dois desses com sobreposição entre os subdomínios, e um sem sobreposição. Dentre os métodos com sobreposição foram estudados os métodos aditivo de Schwarz e multiplicativo de Schwarz. Já dentre os métodos sem sobreposição optou-se pelo método do complemento de Schur. Todas as implementações foram desenvolvidas para serem executadas em clusters de PCs multiprocessados e estão incorporadas ao modelo HIDRA, que é um modelo computacional paralelo multifísica desenvolvido no Grupo de Matemática da Computação e Processamento de Alto Desempenho (GMCPAD) para a simulação do escoamento e do transporte de substâncias em corpos de águas.

Análise dos principais fatores intervenientes no comportamento hidrodinâmico do banhado do Taim (RS)

Os ecossistemas conhecidos como banhados são corpos d’água rasos de grande importância ecológica, embora apenas recentemente isso começou a ser reconhecido. Além de servirem de habitat para uma diversidade de espécies de seres vivos e constituírem áreas de transição entre ambientes aquáticos e terrestres, os banhados ainda são responsáveis pela atenuação de cheias, melhoria da qualidade da água e recarga/descarga de águas subterrâneas. Ao lado da oscilação natural do nível da água, conhecido como hidroperíodo, a circulação da água desempenha um papel de fundamental importância na estrutura e no funcionamento dos banhados. Tais ecossistemas são caracterizados por pequenas lâminas d’água e pela grande presença de macrófitas aquáticas, que são sensíveis a alterações no nível e no padrão de circulação da água. Isso tem ocorrido no banhado do Taim (RS), objeto deste estudo, em decorrência da ação antrópica, como a captação de água para orizicultura e a operação das comportas de saída para a lagoa Mirim. Neste trabalho, aplicou-se um modelo hidrodinâmico bidimensional horizontal (modelo IPH-A), associado a um algoritmo de secagem/inundação, ao sistema formado pela lagoa Mangueira e pelo banhado do Taim, procurando analisar os principais fatores intervenientes no comportamento hidrodinâmico do banhado. Foram avaliados os efeitos decorrentes da ação de cisalhamento do vento na superfície da água, da resistência ao escoamento e da proteção contra a ação do vento proporcionados pelas macrófitas aquáticas e a interação com a lagoa Mangueira. A resistência da vegetação ao escoamento foi representada na modelagem pela adoção de um coeficiente de Manning superior ao das áreas consideradas sem vegetação e, para representar o efeito protetor contra a ação do vento, considerou-se uma redução da tensão de cisalhamento na superfície da água nas áreas com vegetação. Verificou-se, através das simulações numéricas, que os fatores analisados têm papel fundamental sobre a hidrodinâmica do Taim, estando fortemente associados entre si. Em particular, os resultados mostraram que a distribuição espacial da vegetação norteia o padrão geral da circulação da água no banhado e sua ocorrência, principalmente na interface com a lagoa, diminui a influência desta sobre o Taim, significativamente. Outro fato relevante é que, conforme o nível da água, diferentes regiões do banhado ficam secas, resultando padrões hidrodinâmicos distintos, que também diferem quanto à influência que sofrem do vento e da lagoa Mangueira.

Formulação semi-analítica para a equação transformada resultante da aplicação da GITT em problemas difusivos-advectivos

Neste trabalho se propõe um avanço para a Técnica Transformada Integral Generalizada, GITT. O problema transformado, usualmente resolvido por subrotinas numéricas, é aqui abordado analiticamente fazendo-se uso da Transformada de Laplace. Para exemplificar o uso associado destas duas transformadas integrais, resolvem-se dois problemas. Um de concentração de poluentes na atmosfera e outro de convecção forçada com escoamento laminar, entre placas planas paralelas, com desenvolvimento simultâneo dos perfis térmico e hidrodinâmico. O primeiro é difusivo, transiente e com coeficientes variáveis. Sua solução é obtida de forma totalmente analítica. Além de mostrar o uso da técnica, este exemplo apesar de ter coeficientes variáveis, é resolvido com o auxílio de um problema de autovalores associado com coeficientes constantes. No segundo, obtém-se a solução da Equação da Energia analiticamente. Já a Equação da Conservação do Momentum é linearizada e resolvida de forma iterativa. A solução de cada iteração é obtida analiticamente.

Estudo teórico-experimental do processo de microfiltração de partículas de sílica em suspensão aquosa

Cada vez mais muitos processos industriais convencionais de separação, tais como destilação, centrifugação e extração com solventes, entre outros, estão sendo substituídos pelos processos de separação por membranas em virtude desses apresentarem vantagens em relação aos processos convencionais, além de serem considerados uma tecnologia limpa. Entre estas vantagens pode-se citar a economia de energia, pois não envolvem mudança de fase; a possibilidade de processamento de substâncias termolábeis, pois podem operar à temperatura ambiente; a maior seletividade; a facilidade de scale – up, pois são modulares e não extensivos em mão de obra, entre outros. Uma característica limitante destes processos é a existência de uma camada de elevada concentração de partículas ou moléculas que se forma próxima à superfície da membrana e que ocasiona uma redução no fluxo permeado, podendo modificar as características de retenção da membrana. Este fenômeno é conhecido como polarização por concentração. O entendimento dos mecanismos de polarização é considerado de grande importância neste campo e o estudo e validação de modelos matemáticos preditivos torna-se necessário para se obter sucesso na aplicação dos processos com membranas no meio industrial O objetivo geral deste trabalho consiste na análise teórico-experimental do processo de microfiltração tangencial com suspensão aquosa de partículas de sílica. O objetivo específico é a utilização de dados experimentais para a estimação de um parâmetro do modelo matemático que descreve o perfil de concentração na superfície da membrana ( membrana). A configuração tangencial implica que a alimentação, contendo as partículas a serem removidas ou concentradas, flui paralelamente à superfície da membrana reduzindo, assim, o fenômeno de polarização por concentração. Os experimentos foram efetuados em uma unidade de bancada de microfiltração tangencial, também conhecido como fluxo cruzado, em um módulo para membrana plana com canal de escoamento retangular sob condições de escoamento laminar. Os seguintes parâmetros foram variados durante os experimentos: pressão transmembrana, velocidade tangencial e concentração de partículas de sílica Os dados experimentais foram utilizados em uma modelagem matemática, baseada nas equações fundamentais da mecânica dos fluidos e da transferência de massa, que considera as propriedades físicas da suspensão, viscosidade dinâmica e difusividade mássica, variáveis com a concentração de sílica para se estimar a concentração máxima de sílica junto à superfície da membrana. Esta modelagem engloba os parâmetros operacionais, bem como as características da membrana nesta única variável. Com a finalidade de se verificar como a concentração máxima sobre a membrana é afetada pelas condições operacionais e, conseqüentemente, como esta afeta o fluxo permeado, a concentração máxima sobre a superfície da membrana é apresentada em função das condições operacionais variadas durante os experimentos de microfiltração. Finalmente, após a estimação desta variável, pôde-se concluir que ela é influenciada pelas condições operacionais típicas do processo de microfiltração tangencial, tais como velocidade tangencial de escoamento, pressão através da membrana e concentração da alimentação.

Estudo teórico-experimental dos efeitos de tamanho e de velocidade de deformação em estruturas de aço

o objetivo principal da pesquisa é analisar a influência conjunta do tamanho e da velocidade de deformação (ou velocidade de alongamento) nas propriedades mecânicas do aço, procurando determinar a possível existência de acoplamento produzido com a variação conjunta dos dois parâmetros. A pesquisa inclui a revisão bibliográfica, a análise teórica do assunto e um estudo experimental com 130 corpos de prova de diferentes tamanhos, extraídos duma única chapa grossa de aço, ensaiados com diferentes velocidades de deformação, obtendo resultados de distintas variáveis (resistência à tração, deformação de ruptura, energia consumida por unidade de volume, deformação correspondente à resistência à tração e redução da área na zona de escoamento). As principais conclusões da análise teórica são: a) Levando em consideração que o conceito de comprimento característico do material é de fundamental importância para interpretar o efeito de tamanho, se determina a forma geral da função f que define o fenômeno, para materiais com um e dois comprimentos característicos. b) Por analogía, para analisar o efeito produzido pela velocidade de deformação se pode supor a existência de uma, o mais de uma, velocidade de deformação característica do material, resultando uma função f que tem a mesma forma geral que no efeito de tamanho. c) Os atuais modelos que interpretam os efeitos de tamanho e de velocidade de deformação podem ser expressados na forma das expressões gerais encontradas. d) Também se determina a forma geral da função.f que define o efeito conjunto de tamanho e de velocidade de deformação quando não existe acoplamento das duas variáveis e no caso geral de um material com um comprimento e uma velocidade de deformação característica Os principais resultados do estudo experimental são: a) A variação das propriedades quando varia o tamanho, a velocidade de deformação, a posição na espessura e a região da qual é extraído o corpo de prova. b) Existe acoplamento dos dois efeitos claramente significativo para a resistência à tração e a deformação correspondente à resistência à tração, apenas significativo para a deformação de ruptura e a energia consumida por unidade de volume e não significativo para a redução da área na zona de escoamento. c) Do estudo experimental da resistência à tração resulta que o acoplamento dos dois efeitos é tal que quando aumenta a velocidade de deformação e diminui o tamanho o efeito conjunto se acentua e quando diminui a velocidade de deformação e aumenta o tamanho o efeito conjunto se atenua. d) Com base nos critérios existentes que analisam os dois efeitos por separado é deíinido um novo critério para a resistência à tração, incorporando termos que permitem interpretar o acoplamento dos dois efeitos, obtendo um melhor resultado que reduz a soma de quadrados dos resíduos em 20 % . As conclusões da tese proporcionam maior informação teórica e experimental dos efeitos de tamanho e de velocidade de deformação nas propriedades mecânicas dos materiais que podem contribuir para adotar decisões na consideração destes efeitos nas estruturas.

Análise da estabilidade do corte ferroviário do KM 112 da linha Cacequi-Rio Grande, RS e proposta de solução

Esta Dissertação de Mestrado tem por finalidade apresentar a Análise da Estabilidade de Um Corte Ferroviário e Proposta de Solução para Estabilização de Taludes de um segmento da linha férrea Cacequi – Rio Grande (RS), Sub-trecho entre São Gabriel e Bagé, localizado entre os km 111 + 600 e km 112 + 800. Este trecho de variante foi implantado por volta de 1975, em substituição ao traçado da linha velha que comprometia o tráfego das composições ferroviárias. Tal segmento está assente sobre a formação Rio Bonito, constituída de paraortoconglomerados, arenitos quartzosos, arenitos feldspáticos, siltitos, folhelhos carbonosos, diamictitos e carvão. O corte em questão é tipo caixão, de altura variável, em torno de 12,00 m na sua parte mais alta. Quando da sua implantação foi concebido com valetões de pé de corte aberto e recobertos com leiva. Desde sua construção, os taludes apresentaram forte erosão superficial e instabilidade geral, ocorrendo, no ano de 1982, uma ruptura circular no lado esquerdo, a qual provocou um deslocamento de cerca de 4,0m da linha férrea da sua posição original em direção ao talude do lado direito (sentido da quilometragem). A massa de solo envolvida apresenta deslocamentos mais evidentes junto ao pé do talude, provocando uma manutenção constante corretiva do nivelamento da via férrea. Também se observa que os valetões estão assoreados, o que ocasiona obstruções e provoca a retenção d´água superficial e o encharcamento do solo, contribuindo na instabilidade do talude Além da instabilidade global, existem instabilidades localizadas (escorregamentos superficiais) que são devidos à erosão superficial dos taludes em razão do inadequado direcionamento das águas de escoamento superficial e da fraca proteção vegetal nos taludes. Nas sondagens foi descoberta uma fina camada de argilito cinza-escuro numa cota muito próxima da cota em que foi definida a plataforma. Esta camada provocou a mobilização da resistência ao cisalhamento residual, ocasionando a instabilização do talude. A proposta de estabilização prevê, a um baixo custo, melhorar em aproximadamente 70% as condições locais, mediante a execução de uma trincheira drenante no pé do talude, em toda a extensão do corte, a instalação de drenos sub-horizontais na região mais afetada e a correção dos sistemas de drenagem existentes, acompanhados de um gerenciamento da execução dos serviços. O monitoramento posterior definirá a necessidade ou não da execução de um muro de contenção em solo reforçado na região mais crítica.

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.

Previsão de vazão em tempo real no rio Uruguai com base na previsão meteorológica

A operação de reservatórios para geração de energia, ou controle de cheias é definida em função dos volumes afluentes que são resultantes das chuvas que ocorrem sobre a bacia. Devido à aleatoriedade e às próprias incertezas envolvidas na ocorrência das precipitações e vazões; a produção de energia, a segurança das barragens e o controle das cheias à montante e jusante ficam comprometidas. Para que as incertezas sejam reduzidas é necessário o aprimoramento das previsões de vazões de afluência em tempo real. A previsão em tempo real pode se realizada com base na vazão de postos de montante e jusante, na precipitação observada e, ou, na precipitação prevista. A previsão de precipitação é necessária para aumentar a antecipação da previsão e melhoria de resultados para tempos futuros além do tempo de concentração da bacia. Esta pesquisa tem como objetivo a avaliação do ganho da previsão de vazão com uso integrado de previsão de precipitação através de uso de um modelo meteorológico regional (meso-escala) com um modelo hidrológico distribuído. Os resultados do modelo meteorológico regional foram fornecidos pelo Laboratório de Planejamento Energético (LabPlan) da Universidade Federal de Santa Catarina (UFSC), onde está sendo utilizado, de forma operacio nal, o Modelo Numérico Regional ARPS (Advanced Regional Prediction System). O modelo hidrológico de transformação chuva-vazão utilizado é um modelo distribuído com discretização em módulos para grandes bacias - MGB (Modelo de Grandes Bacias). O estudo de caso foi realizado na bacia do rio Uruguai até a Usina Hidrelétrica de Machadinho, cuja área de drenagem é de, aproximadamente, 32.000 km2. Diversos cenários de previsão foram simulados. Para o período de 2001 e 2002 foi feita a análise das previsões de eventos isolados, segundo a disponibilidade de dados de previsão meteorológica. Para o período de 2003, durante 6 meses, foi feita a análise das previsões contínuas. Para este período, através de algumas estatísticas, avaliou-se o ganho hidrológico obtido, em termos de vazão prevista com utilização do modelo hidrológico chuva -vazão considerando chuva futura zero e considerando a previsão da chuva com modelo meteorológico regional. Para o período de 2001 a 2003 avaliou-se, também, a importância da rede de pluviógrafos para previsão em tempo real. Formas de atualização simples das variáveis de estado foram testadas e mostraram significativa melhora das previsões. Os resultados da previsão por eventos mostraram ganhos significativos na previsão de vazão quando a previsão de chuva foi incorporada. Já no período de previsão contínua o mesmo não foi observado, porém este período foi bastante seco com poucos eventos de cheia prejudicando a análise do uso das previsões de chuva no modelo hidrológico para previsão. A análise da importância da rede de pluviógrafos destacou a região sul da bacia como a região mais importante em termos de geração de escoamento rápido ao reservatório de Machadinho. Além disso, uma análise simplificada mostrou que uma rede de pluviógrafos distribuídos na bacia, segundo as recomendações da Organização Meteorológica Mundial (OMM), poderia reduzir em aproximadamente 25% o erro padrão nas previsões de vazão com 12 horas de antecedências em Machadinho.

Blendas de polipropileno com elastômeros termoplásticos tipo SBS e SEBS

Neste trabalho, foram preparadas blendas de polipropileno (PP) com elastômeros termoplásticos (TPE) tipo SBS (poliestireno-bloco-polibutadieno-blocopoliestireno) e tipo SEBS (poliestireno-bloco-poli(etileno-co-butileno)-blocopoliestireno) a fim de se avaliar a influência do tipo e concentração de elastômero nas propriedades do polímero final. Foram utilizados como matrizes de polipropileno um homopolímero (PP-H) e um copolímero aleatório de propileno-co-etileno (PP-R). A fim de avaliar o efeito da presença de um agente nucleante na matriz PP-R, algumas blendas de PP-R/TPE foram preparadas utilizando-se dibenzilideno sorbitol (DBS). Os resultados mostraram que o TPE teve efeito nucleante na cristalização do PP, aumentando a temperatura de cristalização (Tc) e estreitando o pico de fusão, sendo que o SBS teve um maior efeito do que o SEBS em ambas matrizes, resultando em maiores valores de Tc. Embora o efeito de nucleação possa mudar a cristalinidade da matriz na blenda, as propriedades mecânicas foram mais suscetíveis às alterações na morfologia ou dispersão do TPE nas matrizes de PP. Estudos de microscopia eletrônica de varredura (MEV) mostraram que as morfologias das blendas com TPE apresentaram-se diferenciadas, dependendo do tipo e principalmente da quantidade adicionada de TPE em cada matriz. Na matriz PP-H observou-se um maior número de domínios de TPE do que nas matrizes PP-R e PP-RN O tamanho médio dos domínios elastoméricos nas blendas PP/SEBS foram menores e mais bem dispersos do que nas blendas PP/SBS, resultando em um material com melhor resistência ao impacto à baixa temperatura. Ensaios de tensão-deformação mostraram que a adição dos TPES diminui a tensão no escoamento e o módulo de Young para todas matrizes, aumentando o alongamento na ruptura de forma mais significativa no caso das blendas PP-R/TPE. Ambos TPES tiveram igual influência nas propriedades de tração, especialmente com a matriz PPR, mostrando que ambos elastômeros podem ser usados como modificadores de impacto, em diferentes formulações que atendam aplicações específicas de mercado.

Alternativas para a superação de obstáculos relativos à criação de valor, aumento da competitividade e acordos comerciais Internacionais na negociação comercial do Brasil com a China: o caso da soja

China and Brazil have got commercial potentials and strategic partnerships are desired from both sides, but in order to guarantee long term, stable relationship these countries must overcome barriers and obtain mutual gains. The objective of this study was to determine how one can overcome obstacles related to value creation, enhance of competitiveness and international deals in terms of Brazil-China commercial negotiation specifically for soy beans. The CEO of one of the biggest and most important international trading-companies has been contacted and mentioned problems faced by his company in terms of brazillian soy beans sale to chinese clients. Problems such as: the break up of legally based contracts, demanding public pressure, insufficient infra-structure in Brazil for production dispatching, cultural differences, lack of trust, power asymmetry were related and has been analyzed on this study. The perspective used on this study refers to elements, processes and forms of negotiation between these two continental giants.

Caracterização de escoamentos turbulentos transientes usando a transformada de ondaletas

Este trabalho apresenta a aplicação da transformada de ondaletas, como ferramenta de análise para o estudo de escoamentos turbulentos transientes e não homogêneos que podem ocorrer nas situações reais de engenharia, onde fenômenos transientes e não contínuos estão normalmente presentes. O estudo experimental da turbulência em túnel de vento habitualmente pressupõe que o escoamento seja estacionário; assim, fenômenos transientes são estudados como uma sucessão de situações estacionárias intermediárias. Isto é necessário porque a ferramenta clássica para o estudo experimental da turbulência, a análise de Fourier, só se aplica a fenômenos estacionários, pois seus resultados se referem a comportamentos de conjunto e as singularidades do sinal não aparecem na análise. Para estes escoamentos onde a transformada de Fourier não se aplica ou não apresenta resultados satisfatórios, a transformada de ondaletas, entre outras possibilidades, é a ferramenta matemática que vem sendo mais freqüentemente utilizada a partir da última década. São apresentados os fundamentos matemáticos, bem como uma breve história da transformada de ondaletas, da transformada de Fourier e da estatística aplicada á turbulência. Para estudar a transformada de ondaletas e buscar a melhor forma de aplicá-la ao estudo da turbulência, foram analisados três escoamentos distintos: na esteira de um cilindro, em turbulência homogênea e em um banco de tubos de configuração quadrangular. Os resultados experimentais foram obtidos por anemometria de fio quente. As medições foram feitas em um túnel de vento, partindo do repouso até atingir o regime permanente. Os dados experimentais obtidos foram analisados utilizando técnicas de ondaletas, usando várias ondaletas e várias transformadas de ondaletas diferentes. Os resultados são comparados entre si e com a transformada de Fourier. Também foram realizados ensaios em regime permanente para várias velocidades do escoamento entre o repouso e a velocidade nominal do ventilador, com o auxílio de um inversor de freqüência. Os resultados para a variação do número de Strouhal com a velocidade, obtidos da análise destes dados pela transformada de Fourier mostraram boa concordância com os obtidos da análise de ondaletas das aquisições transientes. Os resultados mostram que a transformada de ondaletas é uma ferramenta valiosa para a análise dos dados experimentais obtidos nos problemas investigados, qualificando e complementando a análise de Fourier onde esta se aplica e substituindo-a quando o fenômeno for não estacionário.

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.

Modelagem estratigráfica de clinoformas deposicionais : construção e aplicação de um modelo computacional baseado em mecânica estatística

Sob a premissa de que a geometria do perfil deposicional das clinoformas pode conter informações sobre a composição textural dos sedimentos que as compõem, bem como sobre a energia da coluna d’água, foi desenvolvido um modelo computacional para a análise teórica dos fatores controladores do tipo de curvatura dos foresets das clinoformas, e, por extensão dos taludes submarinos. Um modelo análogo de suspensão dinâmica foi implementado experimentalmente com um programa em código Matlab, denominado MAXWELL, e é classificado como um modelo estratigráfico, bidimensional, analítico, determinístico, dinâmico, parcialmente dedutivo e parcialmente baseado em regras. Contém um modelo de escoamento de fluido em linhas de corrente, e trata indiretamente a tensão de cisalhamento no domínio de um sistema fechado bifásico líquido-vapor análogo, a partir a proposta de analogias entre os processos de evaporação e de suspensão de sedimentos. É uma abordagem baseada na competência de transporte do fluxo d’água, pois considera temperatura e velocidade de fluxo combinado onda-corrente como variáveis controladoras análogas da intensidade de transferência de energia. Os processos deposicionais marinhos são reduzidos a um fenômeno da superfície deposicional, que é considerada análoga a uma interface líquidovapor dentro de um sistema fechado. A equação de distribuição estatística de velocidades moleculares de Maxwell é usada para calcular a proporção de moléculas na fase líquida, considerada análoga à proporção de sedimentos depositados, bem como a proporção na fase vapor, tomada como análoga à proporção de sedimentos em suspensão. O estado energético do sistema é parametrizado por três potenciais: energia interna, como função do tamanho de grão (areia, silte ou argila); energia térmica do meio, como função da energia hidrodinâmica, e energia gravitacional, como função do declive topográfico As simulações indicam que os principais fatores controladores do perfil deposicional das clinoformas, e por extensão, dos taludes submarinos em geral, são a energia hidrodinâmica da coluna d’água e a granulometria (ou coesão) do sedimento, que podem ser consideradas dois parâmetros comutáveis, isto é, grãos grossos ou coesos podem produzir sobre a geometria das clinoformas o mesmo efeito da baixa energia hidrodinâmica, e vice-versa. Com base no fator de decaimento da energia hidrodinâmica com o aumento da batimetria, foram definidas três ordens de grandeza da intensidade da energia da coluna d’água: baixa energia (10–1), alta x energia (10-2) e muito alta energia (10-3). Com base nesse critério, foram caracterizados quatro tipos de perfis deposicionais de clinoformas: dois tipos sigmoidais e dois tipos exponenciais. Os tipos sigmoidais podem ser de alta energia ou de muito alta energia, e distinguem-se pela granulometria do sedimento e pela distribuição de declividades em relação à dimensão na direção horizontal. Os perfis de tipo exponencial podem ser tanto de baixa energia, quanto de alta energia. A subida relativa do nível do mar afeta o tipo geométrico do perfil deposicional apenas de forma indireta, a partir da retenção da fração grossa do influxo sedimentar na plataforma ou topset. Os principais fatores controladores do ângulo de declividade do perfil das clinoformas são a granulometria do sedimento, a energia hidrodinâmica, a profundidade d’água da bacia e o desvio lateral da corrente de fundo. A dedução da litofácies a partir da geometria das clinoformas parece promissora apenas para os tipos teóricos sigmoidais, que são distintos na forma e no conteúdo sedimentar.

Solução da equação advectivo-difusiva utilizando regras para manipulação de exponenciais de operadores diferenciais e simetrias de Lie : aplicações em engenharia ambiental

Vários métodos analíticos, numéricos e híbridos podem ser utilizados na solução de problemas de difusão e difusão-advecção. O objetivo deste trabalho é apresentar dois métodos analíticos para obtenção de soluções em forma fechada da equação advectivo-difusiva em coordenadas cartesianas que descreve problemas de dispersão de poluentes na água e na atmosfera. Um deles é baseado em regras de manipulação de exponenciais de operadores diferenciais, e o outro consiste na aplicação de simetrias de Lie admitidas por uma equação diferencial parcial linear. Desenvolvem-se regras para manipulação de exponenciais de operadores diferenciais de segunda ordem com coeficientes constantes e para operadores advectivo-difusivos. Nos casos em que essas regras não podem ser aplicadas utiliza-se uma formulação para a obtenção de simetrias de Lie, admitidas por uma equação diferencial, via mapeamento. Define-se um operador diferencial com a propriedade de transformar soluções analíticas de uma dada equação diferencial em novas soluções analíticas da mesma equação. Nas aplicações referentes à dispersão de poluentes na água, resolve-se a equação advectivo-difusiva bidimensional com coeficientes variáveis, realizando uma mudança de variáveis de modo a reescrevê-la em termos do potencial velocidade e da função corrente correspondentes ao respectivo escoamento potencial, estendendo a solução para domínios de contornos arbitrários Na aplicação referente ao problema de dispersão de poluentes na atmosfera, realiza-se uma mudança de variáveis de modo a obter uma equação diferencial parcial com coeficientes constantes na qual se possam aplicar as regras de manipulação de exponenciais de operadores diferenciais. Os resultados numéricos obtidos são comparados com dados disponíveis na literatura. Diversas vantagens da aplicação das formulações apresentadas podem ser citadas, a saber, o aumento da velocidade de processamento, permitindo a obtenção de solução em tempo real; a redução da quantidade de memória requerida na realização de operações necessárias para a obtenção da solução analítica; a possibilidade de dispensar a discretização do domínio em algumas situações.