Orientador: Luiz Fernando Milanez / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-19T23:56:56Z (GMT). No. of bitstreams: 1 Santim_ChristianoGarciadaSilva_M.pdf: 4969690 bytes, checksum: ce0fa0a331e41c2ba678acb1f055dad6 (MD5) Previous issue date: 2012 / Resumo: O presente trabalho trata do estudo numérico da fusão de materiais de mudança de fase (PCMs-Phase Change Materials) no interior de uma cavidade retangular e de um cilindro horizontal, usando-se de um aquecimento isotérmico na parede de tais geometrias, utilizando-se de um código CFD. Inicialmente o material sólido encontra-se sub-resfriado. O objetivo é verificar a influência do aquecimento no processo de fusão, bem como o efeito do subresfriamento neste. Algumas hipóteses simplificadoras foram assumidas no modelo matemático. As propriedades termofísicas dos PCMs foram consideradas constantes, exceto para a densidade, onde a aproximação de Boussinesq fora utilizada e, para a fusão do gelo, empregou-se a relação proposta por Gebhart e Mollendorf (1977). O problema foi resolvido usando uma formulação matemática baseada no método da entalpia-porosidade, o qual permite o uso de uma grade espacial fixa. As grades computacionais utilizadas foram selecionadas após verificar-se que os resultados não variavam mais com o refinamento. O fluxo total de calor na superfície interna das paredes também fora obtido para os sistemas, bem como perfis de temperatura, linhas de correntes e perfis das frentes de fusão. Correlações para os tempos totais de fusão dos PCMs em função dos números adimensionais de Stefan e Rayleigh também foram apresentados. Os resultados numéricos foram comparados a resultados experimentais disponíveis na literatura, apresentando boa concordância / Abstract: The present work deals with the numerical study of melting of PCMs inside a rectangular cavity with all isothermally heated walls and inside a horizontal cylinder with isothermally heated wall, utilizing a CFD code. Initially the solid material is sub-cooled. The objective is to verify the influence of the heating in the melting process as well as the sub-cooling effect. Some simplifying hypotheses were assumed in the mathematical model. The thermophysical properties of phase change material were considered constants, except for the density where the Boussinesq aproximation was utilized and, for the ice melting, a relation proposed by Gebhart and Mollendorf (1977) was employed. The problem was solved by using a mathematical formulation based on the enthalpy-porosity method, which allows the use of a fixed spatial grid. The computational grid used was selected after checking that the results did not vary with more refinement. The total heat flux on the inner surfaces was obtained for the systems, as well as several temperature profiles, streamlines and melting front positions. Correlations for the total PCM melting time as a function of the dimensionless numbers of Stefan and Rayleigh were presented. Numerical results were compared with experimental results suitable in the literature showing be in good agreement / Mestrado / Termica e Fluidos / Mestre em Engenharia Mecânica

Orientador : Luiz Fernando Milanez / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-04T00:03:06Z (GMT). No. of bitstreams: 1 DiasJunior_Tito_D.pdf: 2993584 bytes, checksum: de14ac584f2ba9e3b006a79418cf574b (MD5) Previous issue date: 2003 / Resumo: Neste trabalho, o enfoque em resfriamento de componentes eletrônicos é adotado, buscando um entendimento dos fenômenos tridimensionais envolvidos na convecção natural sobre fontes de calor pontuais ou discretas. Para isso utilizou-se simulação numérica por volumes tinitos implementada no FLUENT 6.0. O escoamento, em regime permanente, foi considerado laminar e as propriedades constantes, exceto a densidade no termo de empuxo, sendo a sua variação com a temperatura aproximada levando-se conta a hipótese de Boussinesq. Especial atenção foi dispensada à modelagem do fenômeno, buscando níveis crescentes de complexidade de modo a monitorar os efeitos da condução e da radiação sobre os resultados. Informações quantitativas e qualitativas sobre o escoamento foram obtidas utilizando-se simulações numéricas visando determinar características macroscópicas do escoamento. Foram avaliados detalhadamente os parâmetros que influenciam a estrutura do escoamento, tais como número de Rayleigh, número de Prandtl, razões de aspecto de cavidades e canais, condições de contorno entre outros. Os casos estudados foram elaborados de forma que os resultados e conclusões pudessem preencher lacunas e esclarecer algumas conclusões presentes na literatura / Abstract: In this work, the electronic component cooling is focused, aiming the understanding of the three-dimensional effects present in the natural convection on discrete or pontual heat sources. The finite volume method, implemented in the FLUENT package, was used to simulate the heat transfer processo The steady state flow was considered laminar and the properties were constant, except the density in the buoyancy term, where the Boussinesq approximation was used. Special attention was given to the physical phenomenon, aiming to increasing complexity levels, in order to monitoring the effects of conduction and radiation on the results. Quantitative and qualitative information about the flow were obtained by numerical simulation aiming to determine macroscopical characteristics of the flow. The parameters influencing the flow structure were evaluated, such as Rayleigh number. Prandtl number, cavities and channels aspect ratios, boundary conditions, and others. The studied cases were planned in order to the results and conclusions could fill in the blanks and to elucidate some conclusions present in the literature / Doutorado / Termica e Fluidos / Doutor em Engenharia Mecânica

Orientador: Luiz Felipe Mendes de Moura / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-24T19:39:53Z (GMT). No. of bitstreams: 1 Silva_JoaoBatistaCampos_D.pdf: 10719408 bytes, checksum: c23aba649022079c516a76319b545a1b (MD5) Previous issue date: 1998 / Resumo: Neste trabalho desenvolve-se um método de elementos finitos com formulação por volumes de controle (CVFEM); usando-se elementos quadrilaterais com 9 nós, onde define-se funções de interpolação lineares e quadráticas; para aplicação em simulação de escoamentos transientes incompressíveis de fluidos viscosos, em domínios bidimensionais. O modelo matemático é constituído pelas equações de Navier-Stokes em variáveis primitivas u,v,p. Alguns problemas testes são resolvidos para validação do programa computacional construído com base no modelo numérico. É feita também uma forma de interpolação para os termos convectivos para eliminar coeficientes negativos nas equações discretizadas e evitar resultados não realísticos, quando os números de Reynolds e de Peclet são altos. Esta interpolação ("upwind") é baseada nos fluxos de massa para os volumes de controle num elemento e é conhecida como interpolação ponderada pela massa (MAW- Mass Weigthed Interpolation). Os resultados obtidos nos problemas testes são comparados com resultados de outros modelos numéricos e/ou resultados experimentais. A concordância entre os resultados é satisfatória / Abstract: In this investigation a control volume-finite element method (CVFEM) using a ninenoded finite element with quadratic interpolation functions has been developed for numerical simulation of incompressible, viscous, transient fluid flows. The Navier-Stokes equations in primitive varibles U,V,pare used as the mathematical model. The mass weighted interpolation (MAW) is also implemented for the convective terms of the equations, into the element considered. The obtained results for some test problems agree quite good with results of other numerical methods and/or experimental results / Doutorado / Termica e Fluidos / Doutor em Engenharia Mecânica

Orientador: Amauri Garcia / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-04T03:07:36Z (GMT). No. of bitstreams: 1 Rosa_DanielMonteiro_M.pdf: 7995773 bytes, checksum: 8859b8e785820d46c1b77212d7acaa31 (MD5) Previous issue date: 2004 / Resumo: Os espaçamentos dendríticos são parâmetros microestruturais importantes, resultantes do processo de solidificação. Eles afetam os perfis de microssegregação e governam a formação de segundas fases na região interdendrítica, influenciando, conseqüentemente, as propriedades mecânicas do material. Poucos estudos têm analisado os efeitos da convecção no líquido interdendrítico, bem como a influência da direção de crescimento em relação à gravidade nos espaçamentos dendríticos. Neste trabalho, foi utilizado um dispositivo de solidificação unidirecional vertical descendente, constituído por uma lingoteira de aço inox com diâmetro interno de 56mm, 150mm de comprimento e 10mm de espessura. Na parte superior, foi posicionada uma câmara de refrigeração a água do mesmo material do molde e com espessura útil de 3,0 mm. Após a obtenção dos lingotes e registrados os respectivos perfis térmicos experimentais, foram determinadas as variáveis térmicas de solidificação: coeficiente de transferência de calor metal/molde, velocidades de deslocamento da isoterma liquidus, gradientes térmicos, tempos locais de solidificação e taxas de resfriamento para a solidificação unidirecional descendente de ligas hipoeutéticas do sistema Al-Cu (Al3%Cu, Al5%Cu e Al8%Cu). Estas variáveis térmicas foram confrontadas com as previsões teóricas de um modelo numérico de solidificação e em seguida correlacionadas com os parâmetros microestruturais experimentais (espaçamentos dendríticos primários e secundários). Dessa forma, foram determinadas equações experimentais de crescimento para a solidificação descendente. As previsões teóricas de modelos de crescimento dendrítico representativos da literatura foram conftontadas com os resultados experimentais. Realizou-se também uma análise comparativa destes espaçamentos dendríticos com resultados obtidos para solidificação vertical ascendente de ligas de mesma composição / Abstract: The dendritic spacings are important microstructural parameters, involved in solidification processes. They can affect not only microsegregation profiles but also the formation of secondary phases within interdendritic regions, which influences mechanical properties of cast structures. A small number of studies have been carried out in order to analyze the effects of melt convection within the interdendritic region or to verify the influence of growth direction on dendritic arm spacings. In this work, an experimental set-up of downward solidification with a stainless steel mold, having an internal diameter of 56mm, height 150mm and wall thickness of 10mm, was used. The upper part of the split mold was closed with a water-cooling chamber made of stainless steel, with a wall thickness of 3 mm. Besides, a combined theoretical and experimental approach is developed to quantitatively determine solidification thermal variables such as: transient metallmold heat transfer coefficient, tip growth rate, thermal gradient, tip cooling rate and local solidification time during downward transient solidification of hypoeutectic AI-Cu alloys solidified downwards (Al-3wt%Cu, Al-5wt%Cu and Al-8wt%Cu). These solidification thermal variables are correlated with dendritic spacings (primary and secondary spacings), which have been measured along cross and longitudinal sections of ingots solidified under downward transient heat flow conditions. Predictive theoretical models for dendritic spacings have been compared with experimental observations. A comparison between upward and downward unsteady-state results for dendritic spacings has also been conducted. / Mestrado / Materiais e Processos de Fabricação / Mestre em Engenharia Mecânica

Orientadores: Philippe Remy Bernard Devloo, Sonia Maria Gomes / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-07T13:02:36Z (GMT). No. of bitstreams: 1 Forti_TiagoLuisDuarte_M.pdf: 2510262 bytes, checksum: bb77d518f6e5cbaeafea1a081b587395 (MD5) Previous issue date: 2005 / Resumo: Este trabalho dedica-se ao estudo do método de Galerkin descontínuo aplicado a problemas de convecção-difusão. O método de Galerkin descontínuo (MGD) é um variante do método de elementos finitos tradicional (MEF) em que as funções do espaço de interpolação são descontínuas entre elementos. A motivação para o estudo do MGD vem da mecânica dos fluidos. Muitos problemas de mecânica dos fluidos apresentam solução com fortes gradientes ou descontinuidades. São os problemas de choque e de camada limite. Nessas regiões de forte gradiente ou descontinuidades, o MEF apresenta oscilações na solução numérica. Essas oscilações tornam o método pouco estável, podendo-se obter soluções não-físicas como pressão negativa. O método de Galerkin descontínuo permite evitar ou reduzir essas oscilações. Várias formulações são disponíveis na literatura e algumas delas são tratadas neste trabalho, em especial a formulação de Baumann e Oden. Propõe-se a combinação dos métodos de elementos finitos e Galerkin descontínuo em uma mesma simulação, obtendo-se as vantagens de cada um deles / Abstract: The present work is dedicated to study the discontinuous Galerkin method (DGM) applied to convection-diffusion problems. The DGM is a variant of the so-known finite element method (FEM). In DGM the interpolation space is formed by discontinuous functions between elements, called a broken-space. The motivation of this work comes from fluid mechanics. Many problems in fluid mechanics have discontinuous solutions or high-gradient solutions. They are shock problems or boundary layer problems. On regions of high-gradient or around discontinuities the FEM presents oscillations in the numeric solution. It makes the FEM unstable and it is possible to obtain non-physical solutions like negative pressures. DGM allows to avoid or reduce these oscillations. Some formulations available in the literature are tested in this work, specially that proposed by Baumann and Oden. In this work, we propose a formulation which combines FEM and DGM in the same simulation obtaining the advantages of both methods / Mestrado / Estruturas / Mestre em Engenharia Civil

Orientador : Elias Basile Tambourgi / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-04T01:56:06Z (GMT). No. of bitstreams: 1 Veloso_MariaAuxiliadoraFortini_D.pdf: 6659840 bytes, checksum: 9158a84e3db6fe7278b8e1cc55b8277e (MD5) Previous issue date: 2004 / Resumo: O programa computacional STHIRP-1 (Simulação Termo-Hidráulica de Reatores de Pesquisa), cujos fundamentos são descritos neste trabalho, utiliza os princípios da técnica de subcanais e tem a capacidade de simular, em condições estacionárias e transitórias, os fenômenos térmicos e hidráulicos que ocorrem no núcleo de um reator de pesquisa refrigerado a água sob regime de convecção natural. Os modelos e correlações empíricos necessários para descrição das grandezas do escoamento que não podem ser descritos por relações teóricas foram selecionados de acordo com as características de operação do reator. Apesar de o objetivo primeiro ser o cálculo de reatores de pesquisa, a formulação utilizada para descrever o escoamento do fluido e a condução térmica nos elementos aquecedores é suficientemente geral para estender o uso do programa a aplicações em reatores de potência e a outros sistemas térmicos que tenham as características representadas pelas equações do programa. Para demonstrar a capacidade analítica de STHIRP-1, foram feitas comparações entre resultados calculados e medidos no reator de pesquisa TRIGA IPR-R1 do CDTN/CNEN. Os resultados indicam que o programa reproduz com boa precisão dados experimentais de temperaturas de saída de subcanais. No entanto, resultados experimentais mais consistentes deverão ser usados no futuro para corroborar a validação do programa / Abstract: The STHIRP-1 computer program, which fundamentals are described in this work, uses the principles of the subchannels analysis and has the capacity to simulate, under steady state and transient conditions, the thermal and hydraulic phenomena which occur inside the core of a water-refrigerated research reactor under a natural convection regime. The models and empirical correlations necessary to describe the flow phenomena which can not be described by theoretical relations were selected according to the characteristics of the reactor operation. Although the primary objective is the calculation of research reactors, the formulation used to describe the fluid flow and the thermal conduction in the heater elements is sufficiently generalized to extend the use of the program for applications in power reactors and other thermal systems with the same features represented by the program formulations. To demonstrate the analytical capacity of STHIRP-1, there were made comparisons between the results calculated and measured in the research reactor TRIGA IPR-R1 of CDTN/CNEN. The comparisons indicate that the program reproduces the experimental data with good precision. Nevertheless, in the future there must be used more consistent experimental data to corroborate the validation of the program / Doutorado / Sistemas de Processos Quimicos e Informatica / Doutor em Engenharia Química

Orientador: Carlos Alberto Carrasco Altemani / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-24T09:07:23Z (GMT). No. of bitstreams: 1 Loiola_BrunaRafaella_M.pdf: 2801156 bytes, checksum: 0bbef2296c5a3a938f63b6246554ab2d (MD5) Previous issue date: 2013 / Resumo: A transferência de calor conjugada por convecção forçada-condução de dois aquecedores protuberantes montados na parede inferior (substrato) de um duto retangular foi investigada através de medidas experimentais em laboratório. Potência elétrica era dissipada em resistências embutidas nos aquecedores feitos de alumínio, que eram resfriados em regime permanente por um escoamento forçado de ar no duto. Cada aquecedor transferia calor por convecção diretamente das suas superfícies em contato com o escoamento e por condução para a placa do substrato através da interface comum. Esta condução de calor espalhava na placa do substrato e era transferida por convecção da placa para o escoamento de ar no duto. O resfriamento conjugado dos aquecedores na placa de substrato foi descrito por uma matriz de coeficientes conjugados obtidos por meio de testes experimentais realizados com um aquecedor ativo por vez, em uma faixa do número de Reynolds do escoamento entre 1600 e 6400. Foram utilizadas duas montagens experimentais praticamente idênticas, uma com a placa do substrato de acrílico e a outra, de alumínio. O número de Nusselt adiabático, que é um descritor invariante do processo convectivo de troca de calor, também foi avaliado nos testes com o substrato de acrílico. Os coeficientes conjugados obtidos nestes testes foram então utilizados para prever as temperaturas dos aquecedores em testes adicionais com os dois simultaneamente ativos, sob condições variadas de escoamento e aquecimento. A concordância dessas temperaturas indicou que esses coeficientes são descritores invariantes do resfriamento conjugado dos aquecedores discretos no duto. Simulações numéricas do escoamento e da transferência de calor no duto considerado foram efetuadas com o pacote computacional PHOENICS em um domínio tridimensional semelhante ao duto experimental. Os resultados numéricos obtidos foram comparados com as medidas experimentais / Abstract: Conjugate heat transfer by forced convection-conduction of two protruding heaters mounted on the lower wall (substrate) of a rectangular duct was investigated through experiments in laboratory. Electric power dissipation in the Aluminum heaters¿ embedded resistances was transferred under steady state conditions to forced airflow in the duct. Heat losses from each heater included direct convection from its surfaces in contact with the airflow and conduction to the substrate plate through their common interface. This heat conduction was spread in the substrate plate and eventually transferred by convection from the plate to the airflow in the duct. The heaters¿ conjugate cooling was described by a conjugate matrix with coefficients obtained from experimental tests with a single active heater at a time, in the range of the airflow Reynolds number from 1,600 to 6,400. Two almost identical assemblies were used, one with a Plexiglas plate and the other with Aluminum. The adiabatic Nusselt number, an invariant descriptor of the convective heat transfer, was also evaluated in the tests with the Plexiglas plate. The conjugate coefficients were used to predict the heaters¿ temperatures from additional tests with both active heaters, for arbitrary airflow rates and heaters¿ power dissipation. An agreement of these temperatures indicated that the conjugate coefficients are invariant descriptors of their conjugate cooling. Numerical CFD simulations were performed using the PHOENICS software in a three dimensional domain similar to the experimental assembly. The numerical results were compared to those of the experiments / Mestrado / Termica e Fluidos / Mestra em Engenharia Mecânica

Circuitos de convecção natural ou sistemas de circulação natural são empregados em diversas áreas da engenharia. Reatores nucleares refrigerados a água utilizam circuitos de circulação natural como método passivo de seguranca. Em situações críticas, sem qualquer controle externo, o sistema permanece em segurança por suas próprias características de funcionamento (intrinsecamente seguro). O trabalho proposto consiste em estudar numericamente o circuito de circulação natural de água, localizado no Instituto de Pesquisas Energéticas e Nucleares / Comissão Nacional de Energia Nuclear em São Paulo, por meio do uso de modelos matemáticos, objetivando determinar o padrão do escoamento em condições sem mudança de fase líquido-vapor. A comparação dos resultados de temperatura obtidos por cada um dos modelos de turbulência aos pontos instrumentados no circuito experimental, na condição transitória, revelou desvios significativos nas respostas do modelo de zero equação. Desvios intermediário foram observados nos modelos de transporte da viscosidade turbulenta (EVTE), k - ω, SST e SSG e resultados melhores foram vericados nos modelos k - ε e DES (com significativa superioridade do primeiro modelo). / Natural circulation loops apply to many engineering applications such as: water heating solar energy system (thermo-siphons), thermal management of electrical components (voltage converter), geothermal energy, nuclear reactors, etc. In pressurized water nuclear reactors, known as PWR\'s, the natural circulation loops are employed to ensure passive safety. In critical situations, the heat transfer will occur only by natural convection, without any external control or mechanical devices. This feature is desired and has been considered in modern nuclear reactor projects. This work consists of a numerical study of the natural circulation loop, located at the Instituto de Pesquisas Energeticas e Nucleares / Comissão Nacional de Energia Nuclear in São Paulo, Brazil, in order to establish the ow pattern in single phase conditions. The comparison of numerical results to experiments in transient condition revealed significant deviations for the Zero Equation turbulence model. Intermediate deviations for the Eddy Viscosity Turbulence Equation (EVTE), k - ω, SST e SSG models. And the best results are obtained by the k - ε e DES models (with better results for the k - ε model).

Foram estudados dois sistemas referentes à células a combustível de membrana de troca protônica através de modelos computacionais. Um protótipo de célula foi utilizado para o estudo da dinâmica de fluidos, tendo sido utilizados resultados cedidos para a validação da formulação do escoamento e posterior quantificação da contribuição do mesmo no transporte de espécies nos meios porosos. Este mesmo sistema foi explorado quanto à geometria dos canais de escoamento, sendo possível estabelecer uma hierarquia de perda por sobrepotencial de concentração em função da estequiometria do reagente. Obteve-se que, para as estequiometrias recomendadas em uma célula real, os melhores desempenhos seriam atingidos por geometrias de serpentina simples e serpentina dupla, ou tripla, para o ânodo e cátodo, respectivamente. Uma célula completa também fora estudada, sendo comparada a resultados experimentais obtidos localmente. Foi observado que o modelo utilizado apresenta boa concordância com os aspectos globais do dispositivo real, enquanto que os aspectos locais são largamente mal representados pelo modelo computacional. Este modelo também fora utilizado para quantificar a contribuição do transporte convectivo na distribuição de reagentes, corroborando e estendendo os estudos realizados com o protótipo de célula. Em suma, são sugeridas as possibilidades e desafios de se utilizar modelos computacionais para células a combustível e sistemas similares, neste caso estabelecendo a conexão entre um protótipo e um dispositivo real. / Two systems related to proton-exchange membrane fuel cells were studied through computational models. A cell prototype was used to study fluid dynamics, where the fluid flow formulation was validated against provided experimental results, being later used to quantify its contribution to species transport in the porous media. This same systems was explored regarding its flow channel geometry, where it was possible to establish a hierarchy relative to the concentration overpotential loss as a function of the reactant stoichiometry. It was obtained that, for the recommended stoichiometries for a real cell, the best performance would be attained by a single serpentine and a double, or triple, serpentine geometries for the anode and cathode, respectively. A full cell was also studied, being compared against experimental results obtained locally. It was observed that the model shows good agreement with the global aspects of the real device, while the local aspects are widely misrepresented. This model was also used to quantify the contribution of convection to reactant distribution, corroborating and extending the studies performed with the cell prototype. In all, the possibilities and challenges concerning the use of computational models with fuel cells and similar systems are suggested, particularly by establishing the connection between a prototype and a real device.

Circuitos de convecção natural ou sistemas de circulação natural são empregados em diversas áreas da engenharia. Reatores nucleares refrigerados a água utilizam circuitos de circulação natural como método passivo de seguranca. Em situações críticas, sem qualquer controle externo, o sistema permanece em segurança por suas próprias características de funcionamento (intrinsecamente seguro). O trabalho proposto consiste em estudar numericamente o circuito de circulação natural de água, localizado no Instituto de Pesquisas Energéticas e Nucleares / Comissão Nacional de Energia Nuclear em São Paulo, por meio do uso de modelos matemáticos, objetivando determinar o padrão do escoamento em condições sem mudança de fase líquido-vapor. A comparação dos resultados de temperatura obtidos por cada um dos modelos de turbulência aos pontos instrumentados no circuito experimental, na condição transitória, revelou desvios significativos nas respostas do modelo de zero equação. Desvios intermediário foram observados nos modelos de transporte da viscosidade turbulenta (EVTE), k - ω, SST e SSG e resultados melhores foram vericados nos modelos k - ε e DES (com significativa superioridade do primeiro modelo). / Natural circulation loops apply to many engineering applications such as: water heating solar energy system (thermo-siphons), thermal management of electrical components (voltage converter), geothermal energy, nuclear reactors, etc. In pressurized water nuclear reactors, known as PWR\'s, the natural circulation loops are employed to ensure passive safety. In critical situations, the heat transfer will occur only by natural convection, without any external control or mechanical devices. This feature is desired and has been considered in modern nuclear reactor projects. This work consists of a numerical study of the natural circulation loop, located at the Instituto de Pesquisas Energeticas e Nucleares / Comissão Nacional de Energia Nuclear in São Paulo, Brazil, in order to establish the ow pattern in single phase conditions. The comparison of numerical results to experiments in transient condition revealed significant deviations for the Zero Equation turbulence model. Intermediate deviations for the Eddy Viscosity Turbulence Equation (EVTE), k - ω, SST e SSG models. And the best results are obtained by the k - ε e DES models (with better results for the k - ε model).