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.

Neste trabalho e apresentado um avanço na tecnica GILTT(Generalized Integral and Laplace Transform Technique) solucionando analiticamente um sistema de EDO's(Equações Diferenciais Ordinarias) de segunda ordem resultante da transformação pela GITT(Generalized Integral Transform Technique). Este tipo de problema usualmente aparece quando esta tecnica é aplicada na solução de problemas bidimensionais estacionários. A principal idéia consiste na redução de ordem do problema transformado em outro sistema de EDO's lineares de primeira ordem e a solução analítica deste problema, pela técnica da transformada de Laplace. Como exemplo de aplicação é resolvida a equação da energia linear bidimensional e estacionária. São apresentadas simulações numéricas e comparações com resultados disponíveis na literatura.

Neste trabalho apresentamos uma solução analítica para a equação de cinética unidimensional de difusão de nêutrons, para o modelo de dois grupos de energia, nas geometrias cartesiana e cilíndrica, pelo método espectral. A ideia básica desta metodologia consiste na expansão da solução para os fluxos de nêutrons rápido e térmico, assim como, a solução para as concentrações de precursores de nêutrons atrasados, em séries de autofunções adequadamente escolhidas para a geometria considerada, ou seja, funções seno para a cartesiana e Bessel para a cilíndrica. Dessa maneira, substituem-se estas expansões na equação de cinética, integrando a equação resultante multiplicada pelas autofunções apropriadas à geometria estudada e, ao usar a propriedade da ortogonalidade, obtemos uma equação diferencial matricial linear de primeira ordem com solução conhecida. Assim, ao lançar mão das autofunções adequadas, mostramos a generalidade desta metodologia para solução deste tipo de problema nas geometrias consideradas. Por fim, apresentamos simulações numéricas e comparações com resultados da literatura para os fluxos de nêutrons e concentrações de nêutrons atrasados. / In this work, we report an analytical solution for the kinetic onedimensional neutron di usion equation for the two-group energy model, in cartesian and cylindrical geometry, by the spectral method. The basic idea of the proposed methodology, relies on the expansion of the fast and thermal neutron uxes, as well, the delay neutron precursors concentrations, in a series of eigenfunctions, properly selected for the geometry considered, we mean sine functions for cartesian geometry and Bessel functions for cylindrical geometry. Replacing these expansions in the kinetic equation, taking moments and using the ortogonality property, we come out with a linear rst order matrix di erential equation with a well known solution. Actually, we have shown the generality of the proposed methodology to work out this kind of problem for the geometries considered using properly the eigenfunctions. Finally, we present numerical simulations and comparisons against literature results for the neutron uxes and delayed neutron precursors concentrations.

Neste trabalho apresentamos uma solução analítica para equações difusivas unidimensionais em geometria cilíndrica da Teoria geral de Perturbação em um cilindro homogêneo pela transformada de Hankel. Apresentamos soluções analíticas para o problema de fonte fixa. Resolvemos também um caso monoenergético em um cilindro heterogêneo utilizando uma fomulação recursiva e também usando a Transformada de Hankel. Foi obtida também uma solução analítica dependente do tempo utilizando a Transformada Finita de Hankel. Conhecendo o fluxo de nêutrons, exceto por constantes de integração aplicamos condições de contorno e de interface, após avaliar estas constantes de integração, obtemos a formulação final para o fluxo nestes casos. Os resultados obtidos neste trabalho foram comparados com a literatura, bem como algumas aproximações especificadas e devidamente explicadas ao longo deste. / In this work we present an analytical solution for the one-dimensional diffusion equation in cylindrical geometry from Perturbation Theory inside a homogeneous cylinder by the Hankel Transform. We present analytical solutions for fixed source problem. We solved also a monoenergetic case inside a heterogeneous cylinder using a recursive formulation together with the Hankel Transform. We obtained an analytical solution for the time-dependent diffusion equation using the Finite Hankel Transform. Once the flux of neutrons except for a constant of integration, we apply boundary and interface conditions. After evaluating this constants we obtain the final flux. The obtained results of this work are compared with the literature, as well as some specific approximations are explain throughout this work.

Existe um grande interesse no entendimento de escoamentos de fluidos não-Newtonianos através de contrações abruptas. Este interesse se origina na importância destes escoamentos em processos poliméricos e no freqüente uso do escoamento com contração abrupta como um problema teste para validação de métodos numéricos. As duas propriedades reológicas mais freqüentemente apontadas e estudadas nos fluidos não-Newtonianos são: a viscosidade dependente da taxa de cisalhamento e a natureza viscoelástica destes fluidos. O caráter viscoelástico do fluido é descrito pelo número de Deborah. Na sua definição usual, o número de Deborah envolve as propriedades elásticas do fluido, isto é, o coeficiente primário de tensão normal ou tempo de relaxação terminal originado dos experimentos com tensões de relaxação. Para caracterizar a viscosidade pseudoplástica dos fluidos não-Newtonianos, utiliza-se a equação viscosidade de Carreau, a qual emprega uma viscosidade para baixas taxas de cisalhamento, um tempo característico igual ao recíproco da taxa de cisalhamento na qual se dá o início da pseudoplasticidade, e a inclinação da região power-law da viscosidade para altas taxas de cisalhamento. Ao número de Deborah, pode-se definir o número de Carreau de modo que ele possa ser utilizado para medir o efeito pseudoplástico nos escoamentos não-Newtonianos. A metodologia de elementos finitos utilizada, conhecida como método Galerkin Mínimos Quadrados (GLS), supera as dificuldades enfrentadas pelo método de Galerkin clássico em escoamentos sujeitos a altos números de Reynolds, a saber: a compatibilização dos subespaços de velocidade e pressão – satisfazendo deste modo a condição denominada de Babuška-Brezzi – e as oscilações espúrias devido a natureza assimétrica da aceleração convectiva da equação de momentum. O método GLS, adicionando termos malha-dependentes à formulação clássica de Galerkin, aumenta a formulação de Galerkin sem, contudo, prejudicar sua consistência. Esta dissertação objetiva a modelagem mecânica e simulação numérica via o método de Galerkin Mínimos-Quadrados, de escoamentos planares e axissimétricos de fluidos pseudoplásticos de Carreau através de contrações abruptas. / There is a great interest in understanding the flow of non-Newtonian fluids through sudden contractions. This interest arises from the importance of these flows in polymeric processes and in the frequent use of the flow with sudden contraction as a problem test for validation of numeric methods. The two most frequently mentioned and studied rheological properties of non-Newtonian fluids are the shear rate dependent viscosity and the viscoelastic nature of these fluids. The Deborah number describes the viscoelastic characteristic of the fluid. The Deborah number, in its basic definition, involves the elastic properties of the fluid, i.e., the primary coefficient of normal tension or time of terminal relaxation, originating experiments with relaxed tensions. The pseudoplastic viscosity of non-Newtonian fluids, usually the Carreau viscosity equation, is characterized by the employment of a low rate shear viscosity, a characteristic time - equal to the reciprocal of the shear rates in the beginning of the pseudoplasticity - and the slope of the power-law region of the viscosity for high shear rates. The Carreau number, with the Deborah number, can be used to measure the pseudoplastic effect in non-Newtonian flows. The finite elements method employed in this work, known as the Galerkin least-squares method (GLS), improve on the difficulties faced by the classic Galerkin method in flows subjected to high Reynolds numbers, i.e., the compatibility of sub-spaces of velocity and pressure (the Babuška-Brezzi condition) and the spurious oscillations due to the asymmetric nature of the convective acceleration of the momentum equation. The GLS method, adding mesh-dependents terms to the classic formulation of Galerkin, increases its stability without harm its consistency. This work has as main objective the mechanical modeling and numerical simulation, via the GLS method, of planar and axisymmetric flow of Carreau pseudoplastic fluids through an abrupt 4:1 contraction. In order to take into account the inertia and pseudoplastic effects, the Reynolds and Carreau numbers and power-law coefficient have been investigated within a range of 0<Re<100 and 0<Cu<100 and 0,2<n<1,0, respectively. All numerical results shown to be physically relevant and in accordance with the literature.

Neste trabalho apresentamos uma solução analítica para equações difusivas unidimensionais em geometria cilíndrica da Teoria geral de Perturbação em um cilindro homogêneo pela transformada de Hankel. Apresentamos soluções analíticas para o problema de fonte fixa. Resolvemos também um caso monoenergético em um cilindro heterogêneo utilizando uma fomulação recursiva e também usando a Transformada de Hankel. Foi obtida também uma solução analítica dependente do tempo utilizando a Transformada Finita de Hankel. Conhecendo o fluxo de nêutrons, exceto por constantes de integração aplicamos condições de contorno e de interface, após avaliar estas constantes de integração, obtemos a formulação final para o fluxo nestes casos. Os resultados obtidos neste trabalho foram comparados com a literatura, bem como algumas aproximações especificadas e devidamente explicadas ao longo deste. / In this work we present an analytical solution for the one-dimensional diffusion equation in cylindrical geometry from Perturbation Theory inside a homogeneous cylinder by the Hankel Transform. We present analytical solutions for fixed source problem. We solved also a monoenergetic case inside a heterogeneous cylinder using a recursive formulation together with the Hankel Transform. We obtained an analytical solution for the time-dependent diffusion equation using the Finite Hankel Transform. Once the flux of neutrons except for a constant of integration, we apply boundary and interface conditions. After evaluating this constants we obtain the final flux. The obtained results of this work are compared with the literature, as well as some specific approximations are explain throughout this work.

Neste trabalho e apresentado um avanço na tecnica GILTT(Generalized Integral and Laplace Transform Technique) solucionando analiticamente um sistema de EDO's(Equações Diferenciais Ordinarias) de segunda ordem resultante da transformação pela GITT(Generalized Integral Transform Technique). Este tipo de problema usualmente aparece quando esta tecnica é aplicada na solução de problemas bidimensionais estacionários. A principal idéia consiste na redução de ordem do problema transformado em outro sistema de EDO's lineares de primeira ordem e a solução analítica deste problema, pela técnica da transformada de Laplace. Como exemplo de aplicação é resolvida a equação da energia linear bidimensional e estacionária. São apresentadas simulações numéricas e comparações com resultados disponíveis na literatura.

Existe um grande interesse no entendimento de escoamentos de fluidos não-Newtonianos através de contrações abruptas. Este interesse se origina na importância destes escoamentos em processos poliméricos e no freqüente uso do escoamento com contração abrupta como um problema teste para validação de métodos numéricos. As duas propriedades reológicas mais freqüentemente apontadas e estudadas nos fluidos não-Newtonianos são: a viscosidade dependente da taxa de cisalhamento e a natureza viscoelástica destes fluidos. O caráter viscoelástico do fluido é descrito pelo número de Deborah. Na sua definição usual, o número de Deborah envolve as propriedades elásticas do fluido, isto é, o coeficiente primário de tensão normal ou tempo de relaxação terminal originado dos experimentos com tensões de relaxação. Para caracterizar a viscosidade pseudoplástica dos fluidos não-Newtonianos, utiliza-se a equação viscosidade de Carreau, a qual emprega uma viscosidade para baixas taxas de cisalhamento, um tempo característico igual ao recíproco da taxa de cisalhamento na qual se dá o início da pseudoplasticidade, e a inclinação da região power-law da viscosidade para altas taxas de cisalhamento. Ao número de Deborah, pode-se definir o número de Carreau de modo que ele possa ser utilizado para medir o efeito pseudoplástico nos escoamentos não-Newtonianos. A metodologia de elementos finitos utilizada, conhecida como método Galerkin Mínimos Quadrados (GLS), supera as dificuldades enfrentadas pelo método de Galerkin clássico em escoamentos sujeitos a altos números de Reynolds, a saber: a compatibilização dos subespaços de velocidade e pressão – satisfazendo deste modo a condição denominada de Babuška-Brezzi – e as oscilações espúrias devido a natureza assimétrica da aceleração convectiva da equação de momentum. O método GLS, adicionando termos malha-dependentes à formulação clássica de Galerkin, aumenta a formulação de Galerkin sem, contudo, prejudicar sua consistência. Esta dissertação objetiva a modelagem mecânica e simulação numérica via o método de Galerkin Mínimos-Quadrados, de escoamentos planares e axissimétricos de fluidos pseudoplásticos de Carreau através de contrações abruptas. / There is a great interest in understanding the flow of non-Newtonian fluids through sudden contractions. This interest arises from the importance of these flows in polymeric processes and in the frequent use of the flow with sudden contraction as a problem test for validation of numeric methods. The two most frequently mentioned and studied rheological properties of non-Newtonian fluids are the shear rate dependent viscosity and the viscoelastic nature of these fluids. The Deborah number describes the viscoelastic characteristic of the fluid. The Deborah number, in its basic definition, involves the elastic properties of the fluid, i.e., the primary coefficient of normal tension or time of terminal relaxation, originating experiments with relaxed tensions. The pseudoplastic viscosity of non-Newtonian fluids, usually the Carreau viscosity equation, is characterized by the employment of a low rate shear viscosity, a characteristic time - equal to the reciprocal of the shear rates in the beginning of the pseudoplasticity - and the slope of the power-law region of the viscosity for high shear rates. The Carreau number, with the Deborah number, can be used to measure the pseudoplastic effect in non-Newtonian flows. The finite elements method employed in this work, known as the Galerkin least-squares method (GLS), improve on the difficulties faced by the classic Galerkin method in flows subjected to high Reynolds numbers, i.e., the compatibility of sub-spaces of velocity and pressure (the Babuška-Brezzi condition) and the spurious oscillations due to the asymmetric nature of the convective acceleration of the momentum equation. The GLS method, adding mesh-dependents terms to the classic formulation of Galerkin, increases its stability without harm its consistency. This work has as main objective the mechanical modeling and numerical simulation, via the GLS method, of planar and axisymmetric flow of Carreau pseudoplastic fluids through an abrupt 4:1 contraction. In order to take into account the inertia and pseudoplastic effects, the Reynolds and Carreau numbers and power-law coefficient have been investigated within a range of 0<Re<100 and 0<Cu<100 and 0,2<n<1,0, respectively. All numerical results shown to be physically relevant and in accordance with the literature.