No description available.

[pt] Neste trabalho foi feita um análise teórica, da influência da relação densidade-temperatura, nas proximidades da densidade máxima, sobre a transferência de calor por convecção natural. Para este estudo foi considerado o sistema formado por uma esfera isotérmica imersa em água: Usando as simplificações da camada limite e a transformação de similaridade, as equação de conservação de massa, momentum e energia, foram reduzidas a sistemas de duas equações diferencias ordinárias, não lineares, de condições de contorno. Estas equações diferenciais simultâneas as quais descrevem os campos de velocidade e temperatura da superfície da esfera, como da temperatura da água. A solução destas equações dão dois tipos de regimes de fluxo, o primeiro o usual fluxo unidirecional e o segundo o bidirecional no qual existe fluxo inverso. Em ambos regimens a direção do fluxo depende tanto da temperatura da superfície da esfera, como da temperatura da água. O objetivo principal deste estudo foi obter a variação do coeficiente de transmissão de calor, o qual depende tanto da temperatura de superfície da esfera, como da temperatura da água. Tendo em vista que a densidade máxima da água ocorre na temperatura de 3,98 graus Celsius , e a temperatura da esfera de 0 graus Celsius a 35 graus Celsius. Os resultados numéricos foram obtidos com o emprego dos computadores digitais IBM-1130 e /370 do Rio Datacentro da Puc. / [en] In this work a theoretical analysis was made on the influence of the temperature density relationship on natural convetion heat transfer in the region of maximum desity. An Isothermal sphere immersed in water was considered in this study. Using boundary layer simplifications and similarity transformationhs, thecontinity, momentum, and energy equations which are non linear and depend on the boudary conditions. These symultaneous differential equations, which describe the velocity and temperature of the sphere as well as the water. The sotution of these equations gives twotypes of flow regimes; the first, the common unidirectional one., and the second, a bidirectional one, in which there is flow reversal. In both of theases regimes of the sphere and water. The principle objective oh this study was to abtain the variation in the heat transfer coefficient wich is dependent on the velocity field, wich in turn is dependent upo the temperatures of the sphere and water. Since the maximum desity of water occurs at 3,98 Celsius degrees, the temperature of the water in this study was varied between 0 Celsius degrees and 20 Celsius degrees, while the temperature of the sphere was varied between 0 Celsius degrees and 35 Celsius degrees. Numerical results were obtain with the use of the IBM-1130 and 370 computors at Rio Datacentro of PUC.

No description available.

<p>Solar thermosiphons integrated into the thermal envelop of buildings has been studied for their potential to take advantage of solar energy in heating buildings. The annual performance of solar thermosiphons cannot currently be predicted with the correlations from previous research. Also, no work has been done on the supply mode of a solar thermosiphon even though it has the potential to provide heating and fresh ventilation air. An investigation was done with the goal of developing a numerical model that could predict the performance of the supply mode of a solar thermosiphon. The numerical model included infrared thermal radiation and conduction through the glass, phenomenon which had not been used in previous numerical models. To validate the numerical model a novel steady state experiment was developed. This experiment included radiation as the heat source and the ability to vary geometric lengths. The performance parameters of mass flow rate and thermal efficiency were comparable between the numerical predictions and experimental results. However, due to uncertainties in the current experimental setup, full validation of the numerical model was not possible. These uncertainties would have to be addressed before the numerical model that was developed can be fully validated and used for generating correlations. After consideration of practical implementation constrains, it was shown that it was easier to implement the indoor air curtain mode of a solar thermosiphon than the supply mode. The indoor air curtain mode provides the same amount of energy from solar radiation to heat a building as the supply mode of a solar thermosiphon.</p> / Master of Applied Science (MASc)

no / The coupling between natural convection and conduction within rectangular enclosure was investigated numerically. Three separate heat sources were flush mounted on a vertical wall and an isoflux condition was applied at the back of heat sources. The governing equations were solved using control volume formulation. A modified Rayleigh number and a substrate/fluid thermal conductivity ratio were used in the range 10^4 −10^7 and 10−10^3 respectively. The investigation was extended to examine high thermal conductivity ratio values. The results illustrated that, when Rayleigh number increased the dimensionless heat flux and local Nusselt number increased and the boundary layers along hot, cold and horizontal walls were reduced significantly. An opposite behaviour for the thermal spreading in the substrate and the dimensionless temperature, were decreased for higher Rayleigh number. Moreover, the thermal spreading in the substrate increased for higher substrate conductivity, which affected the temperature level. However the effect of the substrate is negligible when the thermal conductivity ratio higher than 1,500. / The full text of book chapters are not available for self deposit under the Publisher's copyright restrictions.

yes / A numerical study of natural convection heat transfer in water filled cavity has been examined in 3-D for single phase liquid cooling system by using an array of parallel plate fins mounted to one wall of a cavity. The heat generated by a heat source represents a computer CPU with dimensions of 37.5∗37.5mm mounted on substrate. A cold plate is used as a heat sink installed on the opposite vertical end of the enclosure. The air flow inside the computer case is created by an exhaust fan. A turbulent air flow is assumed and k-ε model is applied. The fins are installed on the substrate to enhance the heat transfer. The applied power energy range used is between 15 - 40W. In order to determine the thermal behaviour of the cooling system, the effect of the heat input and the number of the parallel plate fins are investigated. The results illustrate that as the fin number increases the maximum heat source temperature decreases. However, when the fin number increases to critical value the temperature start to increase due to the fins are too closely spaced and that cause the obstruction of water flow. The introduction of parallel plate fins reduces the maximum heat source temperature by 10% compared to the case without fins. The cooling system maintains the maximum chip temperature at 64.68°C when the heat input was at 40W that is much lower than the recommended computer chips limit temperature of no more than 85°C and hence the performance of the CPU is enhanced.

Computational fluid dynamics was used to simulate solar chimney power plants and investigate modeling techniques and expected energy output from the system. The solar chimney consists of three primary parts: a collector made of a transparent material such as glass, a tower made of concrete located at the center of the collector, and a turbine that is typically placed at the bottom of the tower. The collector absorbs solar radiation and heats the air below, whereby air flows inward towards the tower. As air exits at the top of the tower, more air is drawn below the collector repeating the process. The turbine converts pressure within the flow into power. The study investigated three validation cases to numerically model the system properly. Modeling the turbine as a pressure drop allows for the turbine power output to be calculated while not physically modeling the turbine. The numerical model was used to investigate air properties, such as velocity, temperature, and pressure. The results supported the claim that increasing the energy into the system increased both the velocities and temperatures. Also, increasing the turbine pressure drop decreases the velocities and increases the temperatures within the system. In addition to the numerical model, analytical models representing the vertical velocity without the turbine and the maximum power output from a specific chimney were used to investigate the effects on the flow when varying the geometry. Increasing the height of the tower increased the vertical velocity and power output, and increasing the diameter increased the power output. Dimensionless variables were used in a regression analysis to develop a predictive equation for power output. The predictive equation was tested with new simulations and was shown to be in very good agreement. / Master of Science

Experimental observations indicate the presence of attached, gravity induced, horizontal buoyant currents above large area fires. Their driving mechanism is indirect and resembles the one observed above heated horizontal plates. Classic plume modelling is satisfactory for providing information for the flow far from the source. In dealing with large areas and directing attention to the flow close to the source, the classic plume theory should fail because the radial pressure gradient that is responsible for the driving of the flow is squeezed in the long and thin classic plume assumption. For this we propose a new plume structure for the description of the buoyant flow above a circular region of large radius L as “The flow field must be divided into three regions. A region where the flow is predominantly horizontal and attached to the surface, a transition region from horizontal to vertical where separation of the attached current takes place, and a region where vertical flow is established and classic plume theory can be applied”. A model for the description of the gross properties of the horizontal currents is developed under the term “horizontal plume”. The modified Richardson number for the horizontal plume a, being analogous to the radius of the large area, is studied asymptotically in the limit a → ∞ and second order uniformly valid semi-analytical solutions are obtained. The hot plate experiment was set up in order to test the model and facilitate its improvement. A chapter is dedicated to the data analysis coming from thermocouple readings and visualisation of the flow using particle image velocimetry.In the remainder of this thesis two classic problems of laminar natural convection are revisited. That of the first order laminar boundary layer above an isothermal circular plate of radius a and the first order laminar boundary layer above the semi- infinite plate inclined to horizontal. In both cases allowances to variable property effects were made through the introduction of a nondimensional parameter λT, with its value set to zero implying the assumption of the Boussinesq approximation. For the circular plate, fourth order series solutions were obtained valid at the edge of the plate where the effects of λT and Prandtl number Pr are studied. Furthermore a finite difference scheme for the numerical solution of the nonsimilar partial integro- differential equation was developed using the Keller Box method and compared with results obtained from the commercial finite element software COMSOL Multiphysics 3.5a. For the semi-infinite plate, fourth order series approximations valid at the edge of the plate were obtained, while an extensive analysis for the effect of λT, Pr and inclination parameter σ was performed on the flow. Positions of the separation points when the inclination is negative (σ < 0) as a function of Pr and λT were recovered.

O método de Lattice-Boltzmann (MLB) vem ganhando destaque nas últimas décadas pela sua capacidade de solução de escoamentos complexos como escoamentos multifásicos e multicomponentes, meios porosos e magnetohidrodinâmicos. Também existem extensões do método para a solução de problemas de advecção-difusão, que permitem acoplamento dessa metodologia para a solução numérica do conjunto de equação Navier-Stokes-Fourier. No trabalho presente, os objetivos iniciais podem ser resumidos no estudo dos conceitos básicos necessários para entender a derivação do método a partir da teoria cinética e estudo do método de Lattice Boltzmann, com operador de colisão mais simples conhecido como operador BGK, para problemas bidimensionais térmicos e fluidodinâmicos. A implementação numérica do método foi realizada em linguagem C, Matlab e CUDA C, com foco na solução de cinco problemas incompressíveis e laminares em regime permanente, conhecidos na literatura: equação viscosa de Burgers, escoamento de Poiseuille com transferência de calor, convecção natural em uma cavidade quadrada, regimes de convecção natural, forçada e mista em uma cavidade com uma das fronteiras móvel, e por fim, convecção forçada em uma cavidade ventilada, com uma entrada e uma saída. Além disso, um pequeno estudo no tempo computacional utilizando três implementações distintas foram testadas: implementação em série, uso de interpolações entre malhas grosseiras como condição inicial para malhas mais refinadas, e por fim, a adição da implementação do código em paralelo. Os ganhos de tempo entre a primeira e segunda estratégia ficaram entre 1.5 e 6, ao passo que o código paralelizado mostrou-se entre 20 e 25 vezes mais rápido que a segunda estratégia testada, comprovando o benefício de utilizar o processamento em paralelo em unidades gráficas. Os resultados obtidos para os problemas foram comparados com outros trabalhos da literatura, mostrando boa concordância para os primeiros quatro problemas estudados. Para a cavidade ventilada, diferenças relativas de até 15.7% no coeficiente de troca de calor e de até 28.38% para o coeficiente de queda de pressão adimensional foram observadas. Análises a respeito dos termos de erros do método apresentado, e a utilização de outras metodologias com método de Lattice Boltzmann, como por exemplo, o emprego de outros operadores de colisão, para trazerem mais estabilidade e precisão, podem elucidar melhor as divergências observados entre o trabalho presente e outros trabalhos da literatura. / Lattice-Boltzmann Method (LBM) has gained attention over last decades due to its performance in solving complex flows such as multiphase and multicomponent flow, porous media and magnetohydrodynamics. There are also further techniques which makes LBM able to solve advection-difussion problems, which allows coupling this methodology to solve Navier-Stokes-Fourier equations. In this work, initial goals are, in a few words, studying main concepts required to understand numerical method origins from kinectic thery, and studying the method implemetantion to solve 2D fluid dynamic and thermal problems. Numerical implemetation was written in C, CUDA C and Matlab languages, keeping in focus five different cases of laminar incompressible flows in steady state: viscous Burgers\' equation, thermal Poiseuille, natural convection in square cavities, natural, forced and mixed convection in a lid driven cavity, and, finally, forced convection in a ventilated cavity, with one inlet and outlet ports. Besides that, a quick study regarding impact in simulation time of three distinct approaches was conducted. First approach consists in implementing a serial code, the second in using interpolation for coarser grids results as initial condition for finer meshes, and the last one is adding parallelized code implementation. Time gains between first and second approach range from 1.5 to 6, while parallelized code was able to converge from 20 to 25 times faster than second approach, confirming great benefits in using graphics processing units. Results obtained from numerical solutions of problems were compared with other works from literature, and a good agreement among them was observed, specially in the first four problems studied. In ventilated cavity problem, relative differences up to 15.7% in heat transfer coefficient and up to 28.38% in dimensionless pressure drop coefficient were observed. A further error terms analysis of method used in this work, and use of distinct approaches, such as different collision operators, in order to increase numerical solution stability and accurarcy, may shed a light on results divergences observed between this work and others from literature.

O coletor solar de tubos evacuados une uma alta absortividade de radiação solar a um ótimo grau de isolamento térmico. Estas características, aliadas a um custo relativamente baixo, fazem deste tipo de coletor o mais utilizado no mundo. Por isso, diversos tipos de abordagens estão sendo adotadas para descrever seu comportamento térmico. Nesta linha, este trabalho foi desenvolvido através de abordagem experimental e teórica, sendo a última subdividida em numérica por CFD e analítica. A abordagem experimental contou com a construção de uma bancada para medida de temperatura e radiação, em um coletor de 24 tubos evacuados acoplados em um reservatório de 178 L. A abordagem por CFD utilizou um modelo tridimensional transiente. Com o modelo numérico validado, utilizando resultados experimentais, simularam-se diferentes condições de operação, em termos de ângulo de inclinação, fluxo de calor sobre o coletor, tamanho do reservatório e temperatura de entrada da água. Em relação à abordagem analítica, esta é dividida em: modelo de irradiância, modelo do tubo e modelo do reservatório. O modelo de irradiância determina a distribuição da radiação solar ao longo da circunferência do tubo. Parte desta radiação é absorvida pelo coletor e transferida para água. Esta quantidade é determinada com o modelo do tubo, que é baseado no método de resistências térmicas. O modelo do reservatório descreve o comportamento térmico da água em seu interior, tanto em aquecimento quanto em resfriamento, analisando a interação energética com o coletor e com o meio externo. O desenvolvimento do modelo do tubo passa pela avaliação da vazão mássica entre o tubo e o reservatório, além disso, o coeficiente de transferência de calor por convecção no interior do tubo deve ser determinado. Tais variáveis são determinadas a partir de uma correlação para o número de Reynolds, a qual foi obtida com resultados da avaliação por CFD e é função dos números de Nusselt, Grashof e Prandtl. A correlação proposta apresentou bom ajuste com os resultados numéricos. Com a bancada de ensaio experimental foram feitas medidas de temperatura da água no reservatório ao longo de alguns dias. Para as mesmas condições do experimento, a temperatura média da água no reservatório foi estimada com resultados da integração dos modelos de irradiância, do tubo e do reservatório. A diferença entre os resultados experimental e teórico foi de 4,8% para a energia acumulada. / The evacuated tube solar collector combines high solar radiation absorptivity to a great thermal insulation degree. These characteristics, combined with a relatively low cost, make this type of collector the most used in the world. Therefore, various types of approaches are being adopted to describe its thermal behavior. In this way, this work was developed through experimental and theoretical approaches, the latter being subdivided into numeric, by CFD, and analytical analysis. For the experimental approach a test bench was built. The tests was carried on a solar collector with 24 evacuated tubes coupled to a 178 L tank, measuring temperature and solar radiation. The CFD approach used a transient three-dimensional model. After the numerical model validation using experimental data, simulations was carried over different operating conditions in terms of angle, heat flux on the collector, tank size and water inlet temperature. The analytical approach is divided into: irradiance model, tube model and tank model. The irradiance model determines the irradiance distribution of solar radiation along the circumference of the tube. Part of this radiation is absorbed by the collector and transferred to water, this amount is determined with the tube model, using the thermal resistance method. The tank model describes the thermal behavior of inside water, both in heating and in cooling, analyzing energy interaction with the collector and the external environment. The development of the tube model involves the assessment of the mass flow rate between the tube and the tank, furthermore the convection heat transfer coefficient inside the tube must be determined. These variables are determined from a correlation for the Reynolds number, which was obtained with evaluation results by CFD. Proposed Reynolds number is a function of the Nusselt, Prandtl and Grashof numbers. The correlation presented a good agreement with the numerical results. Using the experimental test bench the water temperature was measured into the tank over a few days. For the same experiment conditions, the average temperature of the water into the tank was estimated by results of integration of irradiance, tube and tank models. The difference between the experimental and theoretical results was 4.8% for the accumulated energy.