Neste trabalho, estendemos, de forma analítica, a formulação LTSN à problemas de transporte unidimensionais sem simetria azimutal. Para este problema, também apresentamos a solução com dependência contínua na variável angular, a partir da qual é estabelecido um método iterativo de solução da equação de transporte unidimensional. Também discutimos como a formulação LTSN é aplicada na resolução de problemas de transporte unidimensionais dependentes do tempo, tanto de forma aproximada pela inversão numérica do fluxo transformado na variável tempo, bem como analiticamente, pela aplicação do método LTSNnas equações nodais. Simulações numéricas e comparações com resultados disponíveis na literatura são apresentadas. / In this work, we extend, in an analytical way, the LTSNformulation to one-dimensional transport problems without azimuthal symmetry. For this problem, we also present the solution with continuous dependence on the angular variable, from which an iterative method of the one-dimensional transport equation solution is established. We also discuss how the LTSN formulation is applied to the resolution of the one-dimensional time dependent transport problem, both in an approximate way, by the numerical inversion of the transformed time flux, and analytically, by the application of the LTSN method into the nodal equations. Numerical simulations and comparisons with results available in the literature are presented.

Durante a realização de projetos de sistemas de refrigeração assistidos por sistemas de armazenagem de energia, na forma latente, é necessário que se saiba de antemão as taxas de transferência de energia dos mesmos. Na literatura soluções para este problema são escassos e não abordam todas as geometrias utilizadas para o armazenamento de energia latente, inexistindo correlaçõesmatemáticas que permitam aos projetistas prever o comportamento da fusão e da solidificação do material de mudança de fase (MMF) e as taxas de transferência de energia térmica. Na solução destes problemas são usualmente adotadas hipóteses simplificativas para amenizar o grau de complexidade da descrição matemática do problema Uma destas hipóteses constitui-se em desconsiderar o processo advectivo na transferência de calor no MMF, em sua fase líquida. Outra hipótese adotada, quando o MMF é a água pura, consiste em utilizar um comportamento linear para a massa específica ou um comportamento médio, consistindo em um afastamento da realidade. Este trabalho descreve esforços no sentido de aproximar estas soluções cada vez mais da realidade fisica do problema. São apresentadas simulações para a mudança de fase da água pura levando em conta a inversão da massa específica, em geometria polar, na forma transiente, utilizando o método de Volumes Finitos para solução das equações governantes. As soluções apresentadas, otimizadas quanto às malhas espacial e temporal, envolvem condições de contorno de primeira e terceira espécies Como resultado das soluções são apresentados o Nusselt local, ao longo da parede do cilindro além do Nusselt médio. Correlações matemáticas para o número de Nusselt médio durante a fusão e para o volume de MMF fundido, são propostos e apresentados. A evolução geométrica da fronteira de mudança de fase também é avaliada para as duas condições de contorno e diversos números de Rayleigh. Comparações entre o volume fundido considerando apenas o processo difusivo e o volume fundido utilizando a convecção são realizadas. Para a verificação da independ~encia de malha aplica-se o índice de convergência de malha (GCI) . Algumas comparações qualitativas com soluções experimentais disponíveis na literatura são apresentadas.

O principal objetivo dessa tese consiste em determinar uma solução numéricada equação bidimensional do transporte de nêutrons para elevadas ordens de quadratura angular. Diagonalizando a matriz de transporte LTSN bidimensional , construímos dois algoritmos que se diferenciam pela forma de representar os termos de fuga transversal, que surgem nas equações LTSN integradas transversalmente. Esses termos no método LTSN2D − Diag são expressos como combinação linear dos autovetores multiplicados por exponenciais dos respectivos autovalores. No método LTSN2D − DiagExp os termos de fuga transversal são representados por uma função exponencial com constante de decaimento heuristicamente identificada com parâmetros materiais característicos do meio. A análise epectral desenvolvida permite realizar a diagonalização. Um estudo sobre o condicionamento é feito e também associamos um número de condicionamento ao termo de fuga transversal. Definimos os erros no fluxo aproximado e na fórmula da quadratura, e estabelecemos uma relação entre eles. A convergência ocorre com condições de fronteira e quadratura angular adequadas. Apresentamos os resultados numéricos gerados pelos novos métodos LTSN2D − Diag e LTSN2D − DiagExp para elevadas ordens de quadratura angular para um problema ilustrativo e comparamos com resultados disponíveis na literatura.

Substâncias carcinogênicas são diariamente lançadas no meio ambiente, tanto na atmosfera quanto em corpos hídricos. Estas são capazes de ligar-se a proteínas constituintes de tecidos vivos produzindo um carcinoma, cuja probabilidade de formação depende da afinidade do poluente com os grupos funcionais presentes nos substratos protéicos. Contudo, os mecanismos pelos quais ocorre a formação do carcinoma não estão totalmente esclarecidos. Alguns modelos baseados em propriedades moleculares foram formulados na tentativa de prever quais serão os mecanismos, compostos intermediários e produtos finais de reação. Porém, esses modelos apresentam sérias limitações por não levarem em conta a dinâmica do processo reativo. Para que se possa estimar os mecanismos, é preciso detectar a formação ou ruptura de ligações ao longo do tempo, o que torna necessário utilizar modelos transientes. O presente trabalho apresenta um modelo que resolve a equação de Schrödinger dependente do tempo para verificar qual o mecanismo de reação entre uma substância carcinogênica e um aminoácido. A simulação do cenário transiente proposto requer um baixo tempo de processamento e possibilita uma fácil interpretação dos resultados obtidos.

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.