Este trabalho visa o uso da função de transferência, a qual relaciona distribuição de temperatura e fluxo de calor, na comparação no domínio freqüência, entre o modelo de difusão usual (parabólico) e um modelo ondulatório (hiperbólico) que inclue o efeito de propagação do calor não instantâneo, sendo avaliados os casos de meio semi-infinito e finito. Para o caso de meio semi-infinito, são determinadas as expressões para as características de amplitude e de fase, considerando tanto a abordagem parabólica quanto hiperbólica. É observada a relação entre estas duas abordagens, mostrando que a abordagem parabólica é uma caso particular da abordagem hiperbólica, podendo ser obtida através de um processo de limite envolvendo o tempo de relaxação r . Para o caso de meio finito, são determinadas as expressões para as caracteríısticas de amplitude de ambas as faces da placa unidimensional, considerando tanto a abordagem parabólica quanto hiperbólica. Estas expressões são transformadas para a forma adimensional quando então são deduzidas as expressões correspondentes das características de amplitude. Mais uma vez, todos os resultados para o caso parabólico podem ser determinados a partir dos resultados do caso hiperbólico, através de um processo de limite envolvendo o tempo de relaxação r São apresentados resultados numéricos referentes às características de amplitude, onde é apontada a existência de uma freqüência limite, acima da qual a diferença entre os dois modelos, do tipo parabólico ou hiperbólico, aumenta rapidamente. Também é apresentada uma forma alternativa de cálculo da distribuição de temperatura transiente que faz uso da função de transferência do sistema.

No trabalho proposto é apresentado um método analítico para a obtenção da solução da equação advectivo-difusiva que rege a propagação de poluentes em corpos hídricos. Este método utiliza transformações conformes para efetuar o mapeamento de um corpo hídrico com geometria arbitrária em um semi-plano, fornecendo uma solução em forma fechada expressa em termos da função corrente e do potencial velocidade.

Este trabalho tem por objetivo simular e analisar uma usina termelétrica a carvão em várias condições de funcionamento. A usina simulada neste trabalho é a AVV 1 localizada em Copenhague, Dinamarca. A AVV 1 é uma usina de geração de potência e aquecimento distrital que pode funcionar em várias condições de carga. A simulação da usina supracitada foi tema de um concurso de simuladores proposto no congresso ECOS 2003 realizado em Copenhague, Dinamarca. Para a realização deste trabalho foi construído um programa na linguagem FORTRAN 90. Cada componente da usina é modelado através de equações de balanço de massa e energia, e o sistema completo tem sua solução obtida pelo método de substituição sucessiva. Para viabilizar essa solução é necessário também implementar uma rotina de cálculo de propriedades do fluido de trabalho. No caso estudado, o fluido de trabalho da usina é a água e a formulação utilizada para o cálculo de suas propriedades nos diversos estados é a IAPWS IF-97. A usina é simulada em dois modos de operação: modo de condensação, onde ocorre apenas geração de eletricidade, e em modo de contrapressão, onde há geração de eletricidade e aquecimento distrital, conforme nomenclatura sugerida pela organização do concurso No modo de condensação, são feitas quatro séries de simulações variando a carga de 100% a 40%. Cada série contém um conjunto de hipóteses quanto à variação das eficiências isoentrópicas e pressões das turbinas em função da vazão mássica. No modo de contrapressão, a usina é simulada funcionando com 100% da carga. O programa desenvolvido calcula as propriedades para qualquer ponto de trabalho ao longo da planta, assim como a eficiência da mesma, a potência gerada, e todas as vazões mássicas pertinentes. Além disso, é feita também uma análise exergética da planta. A simulação demonstrou que a planta possui uma eficiência global de 42,02% com uma geração de 250,2 MW em 100% de carga no modo de condensação. Nessas mesmas condições, do ponto de vista exergético, a eficiência encontrada é de 37,21%. No modo de contrapressão, a usina apresenta uma eficiência exergética de 40,19% com um aproveitamento energético de 90,55%. Por fim, é possível também verificar a comportamento da eficiência da planta e a variação de água de resfriamento do condensador com a carga. Os resultados gerados são próximos àqueles encontrados pelos diversos pesquisadores que abordaram o problema.

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.

Esta dissertação tem por objetivo propor melhorias no desempenho térmico de condensadores evaporativos a partir do estudo de dois modelos desenvolvidos para representar a transferência de calor e a efetividade deste tipo de equipamento, chamados de Ajuste Global e Psicrométrico. São elaborados dois aplicativos computacionais para a simulação dos modelos. Dados experimentais obtidos em uma bancada de testes são usados para o ajuste dos modelos. São apresentados gráficos comparativos do erro relativo de cada modelagem com relação ao experimental. É verificado que os dois modelos apresentam boa exatidão e podem ser utilizados para representar o equipamento. A técnica de análise de incertezas é aplicada no intuito de verificar a influência das diferentes grandezas no cálculo da taxa de transferência de calor do condensador. A temperatura de condensação do fluido refrigerante, juntamente com as vazões de amônia e do ar são os principais parâmetros que apresentam uma maior incerteza, necessitando mais cuidado em sua aquisição. Uma tabela de correção de capacidades é elaborada a partir do modelo Psicrométrico para prever o comportamento do condensador operando em faixas de condições de temperatura de condensação e temperatura de bulbo úmido do ar, diferentes das condições nominais. São propostas modificações na atual estrutura física do condensador evaporativo a fim de melhorar seu desempenho térmico, baseadas no estudo dos modelos e também dos diferentes ensaios realizados na bancada de testes. / The present work aims to propose improvements in the thermal performance of evaporative condensers from the study of two mathematical models developed to represent its heat transfer and effectiveness, called here Global adjustment and Psychometric model. It was elaborated two programs using EES to simulate both models. Experimental data was collected for the adjustment of both approaches. It was presented several graphs comparing the error of each model to the experimental. It was verified that both models presented good accuracy and they may be used to represent the equipment. Uncertainty analysis was realized to verify which variable affects mostly the heat transfer rejected by the evaporative condenser. The ammonia temperature of condensation, the air and ammonia mass flows are the variables presenting dominant influence on the uncertainty being necessary good care in the measurement of those values. A correction table was elaborated using Psychometric model to predict heat transfer in the EC operating in different ammonia temperature of condensation and air wet bulb temperature, different from the nominal ones. Concluding the study, some modifications are proposed in the current structure of the evaporative condenser aiming to improve its performance, based on both models and also the tests realized.

Nessa dissertação apresenta-se o desenvolvimento de um modelo espectral para o decaimento da turbulência na Camada Limite Convectiva (CLC). A equação dinâmica para o espectro de energia é obtida a partir das equações de Navier-Stokes [Rinze, 1975J, na qual o termo de produção de energia por efeito mecânico foi desprezado. Durante o período de transição dia-noite foi considerado que o fluxo de calor da superfície é interrompido instântaneamente, o que permite desconsiderar o termo de produção ou perda de energia por efeito térmico. Para parametrizar o termo de transferência de energia cinética, para um fluxo turbulento e isotrópico, foi utilizado o modelo de Heisenberg, que descreve o decaimento da turbulência como um processo em que os turbilhões de um certo tamanho transferem energia por efeito de uma viscosidade cinemática turbulenta, aos turbilhões menores através de interações em cascata. Para resolver a equação diferencial não linear, que descreve o decaimento da turbulência na Camada Limite Convectiva (CLC), fixa-se o parâmetro "z" e considera-se que o espectro de energia pode ser dividido como um produto de duas funções. Neste trabalho, foi utilizado o modelo de Kristensen [1989] para obter o espectro inicial tridimensional.Foi calculada a componente vertical do coeficiente difusão, considerando o espectro vertical, que foi obtido usando-se uma função peso, a qual informa a contribuição de cada componente na formação do espectro tridimensional. Finalmente, os resultados obtidos no trabalho foram confrontados com os dados de LES - Large Eddy Simulation [Nieuwstadt e Brost, 1986] existentes na literatura. / This work presents the development of a spectral model for the decay of turbulence in the Convective Boundary Layer (CBL). The dynamic equation for the energy spectrum is obtained from the Navier-Stokes equations [Hinze, 1975J, in which the term of energy production by mechanical effect was disregarded. During the period of transition from day to night it was considered that the flow of heat from the surface is instantaneously interrupted, which aIlowed the dismissal of the term of production or loss of energy instantaneously. For the purpose of parameters for the kinetic energy term of transference in a turbulent isotropic flow, the Heisenberg model was used, which describes the decay of turbulence as a process in which certain sized eddy transfer energy to smaIler eddy by means of turbulent kinetic viscosity through cascading interactions. For the purpose of solving the non-linear differential equation describing the decay of turbulence in the Convective Boundary Layer (CBL), we established the "z" parameter and considered that the energy spectrum can be divided as a product of two functions. This paper uses the Kristensen model [1989] for obtaining the initial tridimensional spectrum. The vertical component of the diffusion coefficient was calculated considering the vertical spectrum, which in turn was obtained by means of a weight function, which informs the contribution of each component in the whole of the tridimensional spectrum.FinaIly, the results obtained were confronted with the LES - Large Eddy Simulation [Nieuwstadt e Brost, 1986] existing in literature, for validation.

Neste trabalho apresentamos uma solução analítica para problemas de transferência radiativa em uma placa para névoa e nuvem, respectivamente, com grau de anisotropia 82 e 299, pela nova versão do método LT SN que usa a decomposição espectral para transformar a matriz LTSN em uma matriz diagonal. Por esse procedimento, resolvemos de uma forma direta, a transformada de Laplace de um sistema simbólico de equações lineares, bem como a inversão do fluxo angular transformado. Isto representa que uma solução analítica é determinada com uma significante redução do tempo computacional. Para validar a nova formulação LTSN, são apresentadas simulações numéricas e comparações com os resultados disponíveis na literatura. / In this work we report an analytical solution for radiative transfer problems in a slab for Haze and Cloud, respectively, with anisotropy of degree 82 and 299, by the new version of the LT SN spectral decomposition method based in a transformation of the LT SN matrix in a diagonal matrix. By this procedure is a straightforward task to solve the Laplace transform symbolic linear system as well to invert the angular fluxo It turns out that an analytical solution is determined, with a significant reduction of the computational time. To validate the new formulation, we also present numerical simulations and comparisons with results available in the literature.

No presente trabalho foi desenvolvido um modelo para o decaimento da energia cinética turbulenta em uma camada limite planetária convectiva. O ponto de partida foi a equação de balanço para a função espectro de energia cinética turbulenta. Para parametrizar o termo de tranferência inercial foi considerado que a energia é transferida de grandes para pequenos turbilhões até a escala de Kolmogorov, onde a energia é dissipada na forma de calor. O termo de fonte ou perda de energia cinética por efeito térmico foi parametrizado considerando que se pode escrevê-lo como o produto de duas funções, onde o decaimento no tempo pode ser expresso por uma função cosseno. O termo de produção de turbulência por efeito mecânico foi parametrizado a partir de análise dimensional e considerando-se a teoria de similaridade de Monin-Obukov. Foi observado que a função espectro de energia decai de forma mais lenta quando se considera a fonte de turbulência mecânica do que quando predomina a turbulência gerada por convecção. A energia cinética turbulenta obtida pelo modelo proposto foi comparada com resultados de LES (Large Eddy Simulation). / In this present work a model for decline of the turbulent kinetic energy in the Convective Boundary Layer was developed. The starting point was the budget equation of the turbulent kinetic energy spectrum. To parametrize the term of inertial transfer, it was considered that the energy is transferred of great to small eddies until the scale of Kolmogorov, where the energy is dissipated in the heat form. In the buoyancy effect term was considered that it can written as the commodity of two functions, where the decline in the time can be express as a cosine function. The turbulence production term generated by shear of the wind (mechanics turbulence) was parametrized from the dimensional analysis and considering it the theory of the similarity of Monin-Obukhov. It was observed that the energy spectrum function to decays slower when is considered the mechanic source of what when predominates the turbulence generated by convection. The turbulent kinetic energy gotten by considered model was compared with results of LES (Large Eddy Simulation).

O objetivo deste trabalho é simular numericamente o escoamento de um fluido viscoplástico através de um cilindro posicionado sobre duas placas planas paralelas, utilizando um algoritmo numérico baseado no método de elementos finitos. A modelagem mecânica do problema é baseada nas equações de conservação de massa e conservação da quantidade de movimento, acopladas ao modelo constitutivo recentemente proposto por Souza Mendes e Dutra [Souza Mendes e Dutra, 2004], por vezes mencionado apenas como fluido SMD ou simplesmente viscoplástico. A modelagem mecânica foi aproximada pelo método de Galerkin Mínimos-Quadrados (GLS), construído para superar a falha numérica presente no método de Galerkin clássico. De forma prévia avalia-se o escoamento de um fluido Newtoniano em uma cavidade forçada, sob condições de inércia desprezível (Re=1) e em condições inerciais onde o número de Reynolds é excursionado entre valores iguais a 100, 400 e 1000. Em seguida, avaliase sob esta mesma geometria, duas outras situações reológicas distintas, pseudoplasticidade e dilatância, nestes casos considerando-se o modelo constitutivo definido para um fluido Newtoniano generalizado (GNL), tendo também como principal objetivo à avaliação do código numérico utilizado no presente trabalho, frente os resultados obtidos na literatura. Por fim aborda-se a análise do escoamento do fluido viscoplástico sob a geometria proposta. Avalia-se nesta condição, variações dos valores do índice de potência n, o valor do número de salto J e principalmente as condições de inércia do escoamento, através do número de Reynolds característico para a função SMD. Em todas as situações analisadas, apresentam-se os resultados numéricos, através das isobáricas de pressão, isolinhas de velocidade e isoregiões de tensão (caracterizando regiões de maior e menor rigidez do escoamento), sempre com a principal finalidade de fornecer importantes informações sobre o fenômeno físico envolvido no escoamento do fluido viscoplástico proposto, além da caracterização da topologia presente em cada condição, podendo na maioria dos casos avaliados serem corroborados com a literatura disponível. / The objective of this text is to simulate numerically the draining of a viscoplastic fluid through a cylinder located between two plain parallel plates, using a numerical algorithm based in the method of finite elements. The modeling mechanics of the problem is based on the equations of mass conservation and movement conservation, connected to the constitutive model recently considered by Souza Mendes and Dutra [Souza Mendes and Dutra, 2004], for times mentioned only as SMD fluid or simply viscoplastic. The modeling mechanics was approached by the method of Galerkin Least-Square (GLS), constructed to surpass the present numerical imperfection of the classic Galerkin method. Previously the text perform an analysis of the Newtonian fluid draining through a lid-driven cavity, under conditions of worthless inertia (Re=1) and in inertial conditions where the Reynolds number is varied between values of 100, 400 and 1000. After that, over the same geometry, two other distinct rheology situations are evaluated, shear-thinning and shear-thickening, in these cases considering itself the defined constitutive generalized Newtonian fluid (GNL) model, also having as main objective to evaluate the numerical code used in this present work, comparing with the results gotten in literature. Finally it is approached analysis of the draining of the viscoplastic fluid under proposed mentioned geometry. It is evaluated in this condition, variations of the values of the power-law index n, value of the jump number J and the conditions of inertia of the draining, through the characteristic of the Reynolds number specified for the SMD function. In all the analyzed situations the numerical results are showing through the isobaric lines, velocity profiles and the stress regions (characterizing regions of greater and minor rigidity of the draining), always with the main purpose to supply important information on the involved physical phenomena in the draining of the considered viscoplastic fluid, beyond the characterization of the present topology in each condition, being able in the majority of the evaluated cases to be corroborated with available literature.

A utilização de canais compostos está bastante presente nos mais diversos campos da engenharia. Na engenharia mecânica, por exemplo, os problemas estão ligados à utilização destas estruturas no interior de trocadores de calor, reatores nucleares e até mesmo em equipamentos eletrônicos. Equipamentos nos quais o maior desafio é a maximização das taxas de transferência de calor. A acuracidade das relações que determinam a transferência de calor em trocadores de calor ou reatores nucleares, passa por uma análise dinâmica destes equipamentos, visto que suas características geométricas podem ser responsáveis pela formação de estruturas coerentes dentro do escoamento principal, que podem vir a serem importantes fontes de perturbações. A proposta deste trabalho é investigar experimentalmente as características do escoamento turbulento em alguns tipos de canais compostos fechados. Utilizando como técnica experimental a anemometria de fio quente. Nesse sentido, também foi investigado a influência dos parâmetros geométricos da seção, tanto na distribuição dessas quantidades quanto nas características dinâmicas do escoamento. No primeiro tipo o fluido é forçado a escoar entre duas placas paralelas de profundidade "p" e separadas por uma distância "d". Essas placas foram fixadas a uma das paredes laterais de um canal aerodinâmico. Na segunda seção as placas foram substituídas por feixes de seções retangulares, também fixados a mesma parede. Esses feixes também eram separados por uma distância "d". Na terceira o canal aerodinâmico foi dividido em dois. Os feixes da seção precedente foram movimentados para o centro, dando origem a dois subcanais principais conectados por uma fenda de largura "d" e profundidade "p". A quarta seção é uma simplificação das seções compostas por feixe de barras. Consiste em duas semi-circunferencias de raio externo igual a 66 mm. As semi-circunferencias foram postas uma em frente a outra, separadas também por uma distancia "d" As mudanças nos parâmetros geométricos das seções, profundidade, "p", largura da fenda, "d" e comprimento da seção "L", deram origem a 20 vinte diferentes seções. Os resultados mostraram, em todas as seções, o desenvolvimento de uma camada de mistura na direção principal do escoamento. Quanto as características dinâmicas do escoamento, estas parecem ser reguladas pelos perfis médios de velocidade. / Compound channels are present in many engineering areas, for example, in nuclear reactors, in cooling systems of electronic devices and in supply water channels. On the mechanical engineering field, the problems are linked to heat exchangers, nuclear reactor fuel elements and even electronic devices, where, the challenge is the rise in heat exchange rates. However, a good prediction on the heat transfer rates in compound channels, pass through the accuracy of the relations used: the heat exchange coefficient in heat exchangers or in nuclear reactors and correlations for the friction factors in a supplying channel. The dynamic analysis of these equipments, considering fluid structure interaction, is also important, since their geometric characteristics can be responsible for the formation of coherent structures on the flow that can become important sources of disturbances. The purpose of this work is the experimental investigation of the characteristics of the turbulent flow in compound channels and the influence of the geometric parameters of the section in its dynamic characterization. By using hot wire anemometry techniques, mean quantities and Reynolds stresses distribution were investigated in four kinds of compound channels. In the first one the fluid flows between two parallel plates, both attached on the lateral wall of aerodynamic channel. These plates have as important dimensional parameters their depth, "p", and width, "d". In second configuration the plates were replaced by rectangular bars. These ones also attached on the same lateral wall, apart from each other by a distance "d". In the third test section the aerodynamic channel was spitted in two. The bars attached on the lateral wall were moved to the middle channel, giving rise to subcahnnels connected by a narrow gap, with depth "p" and width "d". The last one was a rod bundle simplification. By using two semicircles, placed in front of one each other, it was possible to create a compound channel. Again, the topology gave rise two main subchannels connected by a narrow gap. Many changes were performed in the geometric parameters during the experiments. By modifying depth, "p", width "d" and the length, "L", enable us to get twenty different test sections. The results showed a mixing layer developing in the streamwise direction. As regards mean, normal and Reynolds stresses quantities, the experiments depicted self-preserving values when made dimensionless by local scales. According to results showed here, dynamical features are mainly ruled by mean axial velocity profile characteristics.