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 : 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: Marcelo Moreira Ganzarolli / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-22T20:18:14Z (GMT). No. of bitstreams: 1 Avelar_AnaCristina_M.pdf: 5757345 bytes, checksum: ab66436b252ec7b8b48d4756dd9e7f42 (MD5) Previous issue date: 1997 / Resumo: Os constantes avanços tecnológicos em eletrônica e informática tem tornado os sistemas eletrônicos cada vez mais compactos, aumentando-se a quantidade de calor a ser removida dos componentes e placas de circuito impresso. Por este motivo, atualmente são exigidos sistemas de dissipação de calor altamente eficientes. Este estudo analisa teórica e experimentalmente a transferência de calor em canais verticais formados por placas de circuito impresso dispostas paralelamente, resfriadas por convecção natural e propõe uma modelagem, baseada em relações existentes na literatura, que busca prever a distribuição de temperaturas em pontos significativos nos canais e placas em função da potência dissipada e da distância entre as placas. Este estudo visa também analisar os efeitos do aquecimento não-uniforme das placas que formam o canal. Na simulação teórica, devido à pequena espessura da placa o gradiente de temperatura ao longo da espessura da mesma foi desprezado e resolveu-se numericamente a equação de transferência de calor em coordenadas cartesianas bidimensionais e em regime permanente. Equacionou-se balanços de energia para os componentes e para o ar no canal e o problema foi resolvido numericamente através de um programa computacional. As placas de circuito impresso utilizadas nos testes experimentais, concebidas exclusivamente para fins de estudos térmicos, possuem uma base de epóxi com 25 resistores discretamente distribuídos sobre sua superfície de 200x164mm. Realizou-se testes com aquecimento uniforme e não-uniforme, variando-se a potência por componentes nas placas. Variou-se também as potências por placas e as distâncias entre as mesmas. Testou-se as potências de 2, 4,6 e 8 W e as distâncias entre placas de 12,24, e 48mm. Verificou-se boa concordância entre os resultados numéricos e experimentais, principalmente para a menor distância entre placas, 12 mm, onde a diferença entre os resultados teóricos e experimentais foi muito pequena / Abstract: The constant technological advances in electronics and computations have made the electronic system increasingly more compact, thus increasing the amount of heat to be removed ITomthe components and printed circuit boards. For this reason, highly efficient system of heat removal are presently required. This study analyses both theoretically and experimentallyheat transfer in an array of vertical parallel printed circuit boards, cooled off by natural convection and proposes a modeling, based on correlation found in literature which tries to predict temperature distribution in significant places in the channels and boards as regards dissipated power and distance among boards. This study also aims at analyzing nonuniform heating effets of the channel made boards. In the theoritical simulation, due to the board' s small thickness, the temperature gradient across the board has been neglected, and the equation of heat transfer in two-dimensional Cartesian was numericallysolved, and on a steady state. Energy balances for the componentes were formulated, and the problem was numerically solved by a computer programo The printed circuit boards used in the experimental tests, manufactured speciallyfor heat trasnfer studies have an epoxy basis with 25 resistors discretely distributed on 200 x 164 mm. Uniform and non-uniform heating tests were performed, thus a variation of power and distance among them had a variation. The 2, 4, 6 and 8 W power and distance among 12, 24 and 48 mm boards were tested. Good agreements between numerical and experimentalresults were observed, mainlyfor the smaller distance among boards, 12 mm, the differencesbetween the theoretical and experimentalresults were small / Mestrado / Termica e Fluidos / Mestre em Engenharia Mecânica

Orientador: Luiz Fernando Milanez / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-23T05:12:48Z (GMT). No. of bitstreams: 1 Malateaux_EzequielCosta_M.pdf: 98402739 bytes, checksum: 8bb545a108688271a9e6d884f1c8ac49 (MD5) Previous issue date: 2013 / Resumo: A transferência de calor em cilindros horizontais isolados é de grande importância na indústria tanto no caso em que se pretende minimizar a dissipação térmica (quando dutos isolados são usados no transporte de fluidos), como quando se pretende maximizá-la (no caso de cabos de eletricidade revestidos, por exemplo). Em ambas as situações, convecção natural e radiação térmica estão presentes, podendo ser consideradas como os principais mecanismos envolvidos de dissipação térmica. Normalmente o isolamento térmico é circular e concêntrico ao cilindro. No entanto, o coeficiente convectivo e a temperatura da superfície não são constantes ao redor do cilindro, levando a crer que a adoção de isolantes concêntricos pode não ser o melhor método de se distribuir o isolante. Assim, se um isolamento excêntrico for adotado, a temperatura superficial deve apresentar variação ainda maior, alterando o escoamento ao redor do cilindro e, portanto, as perdas térmicas. No presente trabalho é apresentada uma análise numérica da dissipação térmica em cilindros horizontais isotérmicos isolados excentricamente expostos a convecção natural (Pr= 0,715) e a radiação. A condução na camada de isolante foi resolvida analiticamente (através do uso do sistema de coordenadas bicilíndricas) e integrada numericamente. A partir da temperatura obtida para a superfície externa do isolante, o software PHOENICS foi utilizado para resolver o problema de convecção natural, fornecendo um novo perfil para o coeficiente convectivo. Para cada caso, algumas iterações foram necessárias para se alcançar a convergência. Os parâmetros analisados são a emissividade da superfície externa, a razão entre os raios externo e interno do isolante, a razão entre as condutividades térmicas do isolante e do ar, o número de Rayleigh (definido de forma a ser independente dos outros parâmetros) e a excentricidade do isolante. A dissipação e a temperatura superficial local são mostradas graficamente, assim como os campos de temperatura e velocidade do escoamento. Uma equação é sugerida para a determinação da influência da excentricidade sobre o calor total dissipado em função da razão entre as resistências térmicas de convecção natural-radiação e a resistência térmica de condução, para determinada razão de raios e excentricidade / Abstract: Heat transfer in horizontal insulated cylinders is of importance in industry as when one wants to reduce the heat dissipation (when insulated ducts are used to transport fluids), as when one wants to enhance it (in the case of covered electric cables, for example). In both cases, free convection and thermal radiation are present and can be considered as the main mechanisms of heat dissipation. Generally, the insulating layer is round and is arranged concentrically in the cylinder. However, the heat transfer coefficient and the surface temperature are not constant around the cylinder, raising the point that the use of concentric insulation may not be the best way to arrange the insulating layer. Thus, if an eccentric layer is used, the surface temperature should have a larger variation, changing the flow around the cylinder and, therefore, the heat dissipation. A numerical analysis of the heat dissipation in horizontal isothermal eccentrically insulated cylinders exposed to free convection (Pr= 0,715) and radiation is presented. The conduction through the insulating layer was solved analytically (using the bicylindrical coordinate system) and integrated numerically. From the solution for the surface temperature, the PHOENICS software was used to solve the problem of free convection, giving a new profile for the heat transfer coefficient. For each case, some iterations were necessary to reach convergence. The parameters analyzed were the emissivity of the external surface, the ratio between the outer and the inner radius of the insulating layer, the ratio between the insulation and the air thermal conductivities, the Rayleigh number (defined in such a way to be independent of the other parameters) and the eccentricity of the insulation. The heat dissipation and the surface temperature are shown graphically, as well as the temperature and velocity flow fields. An equation is suggested to define the importance of the eccentricity over the total heat as a function of the ratio between the free convection-radiation thermal resistance and the conduction thermal resistance, for a given ratio of radius and eccentricity / Mestrado / Termica e Fluidos / Mestre em Engenharia Mecânica

Orientador: Kamal Abdel Radi Ismail / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-24T07:44:11Z (GMT). No. of bitstreams: 1 Martinelli_ValdisleyJose_M.pdf: 3830645 bytes, checksum: fa60238d6afea2e97d6f6e2886a91fc4 (MD5) Previous issue date: 1998 / Resumo: O objetivo deste trabalho é estudar o desenvolvimento de transferência de calor laminar por convecção natural em canais de placas planas verticais sujeito a aquecimento assimétrico. São analisadas condições de contorno com fluxo de calor uniforme nas paredes das placas e temperatura da parede uniforme. Soluções para escoamento em regime permanente são obtidos para o ar para diferentes fluxos de calor e diferentes temperaturas nas paredes (acima da temperatura do fluido na entrada do canal); bem como para diferentes geometrias do canal (variação da distância entre as placas (b)). As equações foram escritas em diferenças finitas utilizando variáveis primitivas e resolvidas implicitamente, usando discretização central em todos os termos das equações governantes, exceto para o caso de fluxo de calor que foi utilizado diferença unilateral de três níveis nos termos da derivada da temperatma. Os resultados indicam que existe uma distância ótima entre as placas onde ocorre a máxima troca de calor entre as placas e o fluido. Os resultados foram comparados com a literatura, tanto trabalhos numéricos como experimentais / Abstract: This work studies the developm.entof laminar free convection in vertical flat parallel plates with asymmetric heating. Thermal boundary conditions of uniform wall heat fluxes and uniform wall temperature are considered. Solutions of the flow on steady state are obtained for air at different wall heat fluxes and wall temperature differences (above the temperature of the fluid at the channel) and also different geometries of the channel (distance between the plates (b)). The equations were written in finite difference using primitive variables and they were solved by implicitesqueme all the terms of the governing equations, and by means of a three point derivative formula the terms wall temperature at heat fluxes. The results show there is one optimal distance between the plates to maximize the heat transfer. The results were compared with the literature, both numerical and experimentalworks / Mestrado / Termica e Fluidos / Mestre em Engenharia Mecânica

Orientador : Amauri Garcia / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-04T03:13:50Z (GMT). No. of bitstreams: 1 Spinelli_JoseEduardo_D.pdf: 12716637 bytes, checksum: 79131dc23b1f9de6e4d77299eeff6973 (MD5) Previous issue date: 2005 / Resumo: 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 nas variáveis térmicas e estruturais durante a solidificação. Neste trabalho, foi utilizado um dispositivo de solidificação unidirecional vertical descendente, constituído por uma lingoteira de aço inoxidável com diâmetro interno de 60 mm, 157 mm de comprimento e 9 mm 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 e taxas de resfriamento para a solidificação unidirecional descendente de ligas hipoeutéticas do sistema Sn-Pb (Sn5%Pb, Sn15%Pb e Sn20%Pb). Estas variáveis térmicas são confrontadas com as previsões teóricas de um modelo numérico de solidificação e, em seguida, correlacionadas com os parâmetros estruturais: espaçamentos dendríticos primários, secundários e transição colunar/equiaxial (TCE). Dessa forma, são determinadas equações experimentais de crescimento para a solidificação descendente e verificado um critério de previsão da TCE. Realiza-se também uma análise comparativa dos presentes resultados com aqueles obtidos para solidificação vertical ascendente de ligas de mesma composição. Os resultados comparativos mostram que o efeito convectivo estimula a ocorrência da TCE e ainda é responsável por uma sensível redução dos espaçamentos dendríticos primários. As previsões teóricas de modelos de crescimento dendrítico representativos da literatura são confrontadas com os resultados experimentais / Abstract: Only a few studies have reported influences of interdendritic convection and growth direction on dendrite arm spacings. A downward directional solidification apparatus was used having a stainless steel split mold with an internal diameter of 60 mm, height 157 mm and a 9 mm wall thickness. The upper part of the split mold was closed with a water-cooling chamber made of stainless steel, with a wall thiclmess of 3 mm. A combined theoretical and experimental approach is developed to quantitatively determine the solidification thermal parameters: transient heat transfer coefficients, tip growth rates, thermal gradients and cooling rates during downward unsteady state solidification ofhypoeutectic Sn-Pb alloys. The resulting thermosolutal convection can start in the melt both within the interdendritic region and ahead of the dendrite array. The experimental results have shown that melt convection may be causing pileup of fractioned dendritic arms, which must stimulate the CET occurrence. The results have supported a criterion recently proposed based on a critical tip cooling rate. For upward unidirectional condition this critical value was found to be about 0.014 K/s for hypoeutectic Sn-Pb alloys. In the present study, in conditions of downward solidification, melt convection seems to favor the structural transition, which is anticipated and occurs for a critical cooling rate of about 0.03 K/s, for any of three hypoeutectic alloys experimentally examined. Primary dendritic arm spacings have been affected by the direction of growth, decreasing in conditions of downward vertical solidification when compared with those grown vertically upwards. The main predictive theoretical models for dendritic spacings are compared with the experimental observations / Doutorado / Materiais e Processos de Fabricação / Doutor em Engenharia Mecânica

CAPES / Vários estudos numéricos investigaram cavidades fechadas sob o efeito da convecção natural preenchidas com fluidos newtonianos generalizados (FNG) nos últimos anos pelas aplicações diretas em trocadores de calor compactos, no resfriamento de sistemas eletrônicos e na engenharia de polímeros. Neste trabalho é realizada a investigação numérica do processo de convecção natural de fluidos de lei de Potência e de Bingham em cavidades fechadas, aquecidas lateralmente e preenchidas com meios heterogêneos e bloco centrado. O meio heterogêneo é constituído de blocos sólidos, quadrados, desconectados e condutores de calor. Como parâmetros são utilizados a faixa de Rayleigh de 104 à 107, índice de potência n de 0, 6 à 1, 6, número de Bingham de 0, 5 até Bimax , sendo investigado da influência do número de Prandtl para cada modelo de fluido. Nas cavidades com meio heterogêneo são utilizadas as quantidades de blocos de 9, 16, 36 e 64, mantendo-se a razão entre a condutividade térmica do sólido e do fluido κ = 1. Para as cavidades com bloco centrado, são utilizados os tamanhos adimensionais de 0, 1 à 0, 9 com κ = 0, 1; 1 e 10. A modelagem matemática é realizada pelas equações de balanço de massa, de quantidade de movimento e de energia. As simulações são conduzidas no programa comercial ANSYS FLUENT R . Inicialmente são resolvidos problemas com fluidos newtonianos em cavidade limpa, seguida de cavidade preenchida com meio heterogêneo e posteriormente bloco centrado para validação da metodologia de solução. Na segunda etapa é realizada o estudo com os modelos de fluidos de lei de Potência e de Bingham seguindo a mesma sequência. Os resultados são apresentados na forma de linhas de corrente, isotermas e pelo número de Nusselt médio na parede quente. De maneira geral, a transferência de calor na cavidade é regida pelo número de Rayleigh, tamanho e condutividade térmica dos blocos, pelo índice de potência para o modelo de lei de Potência e do número de Bingham para o modelo de Bingham. O número de Prandtl tem grande influência nos dois modelos de fluidos. O meio heterogêneo reduz a transferência de calor na cavidade quando interfere na camada limite térmica para ambos os fluidos, sendo feita uma previsão analítica para o fluido de lei de Potência. Para bloco centrado, a interferência na camada limite com fluido de lei de Potência também foi prevista analiticamente. A transferência de calor aumentou com bloco de baixa condutividade térmica e pouca interferência e com bloco de alta condutividade térmica e grande interferência, para ambos os fluidos. / Many studies have been carried out in square enclosures with generalized Newtonian fluids with natural convection in past few years for directly applications in compact heat exchangers, cooling of electronics systems and polymeric engineering. The natural convection in square enclosures with differently heated sidewalls, filled with power-law and Bingham fluids in addition with heterogeneous medium and centered block are analyzed in this study. The heterogeneous medium are solid, square, disconnected and conducting blocks. The parameters used are the Rayleigh number in the range 104 - 107 , power index n range of 0, 6 - 1, 6, Bingham number range of 0, 5 - Bimax , being the influence of Prandtl number investigated for each fluid model. The number of blocks for heterogeneous medium are 9, 16, 36 and 64, keeping constant solid to fluid conductive ratio, κ = 1. For enclosures with centered block are used the nondimensional block size from 0, 1 to 0, 9, with solid to fluid conductive ratio in range κ = 0, 1; 1 and 10. Mathematical modeling is done by mass, momentum and energy balance equations. The solution of equations have been numerically solved in ANSYS FLUENT R software. Firstly, numerical solutions for validation with Newtonian fluids in clean enclosures are conducted, followed by enclosures with heterogeneous medium and centered block. Subsequently, numerical solutions of power-law and Bingham fluids with same enclosures configurations are conducted. The results are reported in the form of streamlines, isotherms and average Nusselt number at hot wall. In general, the heat transfer process in enclosure is governed by Rayleigh number, size and thermal conductivity of the blocks, power index n for power-law fluid and Bingham number for Bingham fluid. Both fluid models are very sensitive with Prandtl number changes. Heterogeneous medium decrease heat transfer in enclosure when affects thermal boundary layer for both fluid models. One analytical prediction was made for power-law fluid. An increase in heat transfer occurs with low thermal conductivity block and few interference and with high thermal conductivity block and great interference, for both fluids.

CAPES / Vários estudos numéricos investigaram cavidades fechadas sob o efeito da convecção natural preenchidas com fluidos newtonianos generalizados (FNG) nos últimos anos pelas aplicações diretas em trocadores de calor compactos, no resfriamento de sistemas eletrônicos e na engenharia de polímeros. Neste trabalho é realizada a investigação numérica do processo de convecção natural de fluidos de lei de Potência e de Bingham em cavidades fechadas, aquecidas lateralmente e preenchidas com meios heterogêneos e bloco centrado. O meio heterogêneo é constituído de blocos sólidos, quadrados, desconectados e condutores de calor. Como parâmetros são utilizados a faixa de Rayleigh de 104 à 107, índice de potência n de 0, 6 à 1, 6, número de Bingham de 0, 5 até Bimax , sendo investigado da influência do número de Prandtl para cada modelo de fluido. Nas cavidades com meio heterogêneo são utilizadas as quantidades de blocos de 9, 16, 36 e 64, mantendo-se a razão entre a condutividade térmica do sólido e do fluido κ = 1. Para as cavidades com bloco centrado, são utilizados os tamanhos adimensionais de 0, 1 à 0, 9 com κ = 0, 1; 1 e 10. A modelagem matemática é realizada pelas equações de balanço de massa, de quantidade de movimento e de energia. As simulações são conduzidas no programa comercial ANSYS FLUENT R . Inicialmente são resolvidos problemas com fluidos newtonianos em cavidade limpa, seguida de cavidade preenchida com meio heterogêneo e posteriormente bloco centrado para validação da metodologia de solução. Na segunda etapa é realizada o estudo com os modelos de fluidos de lei de Potência e de Bingham seguindo a mesma sequência. Os resultados são apresentados na forma de linhas de corrente, isotermas e pelo número de Nusselt médio na parede quente. De maneira geral, a transferência de calor na cavidade é regida pelo número de Rayleigh, tamanho e condutividade térmica dos blocos, pelo índice de potência para o modelo de lei de Potência e do número de Bingham para o modelo de Bingham. O número de Prandtl tem grande influência nos dois modelos de fluidos. O meio heterogêneo reduz a transferência de calor na cavidade quando interfere na camada limite térmica para ambos os fluidos, sendo feita uma previsão analítica para o fluido de lei de Potência. Para bloco centrado, a interferência na camada limite com fluido de lei de Potência também foi prevista analiticamente. A transferência de calor aumentou com bloco de baixa condutividade térmica e pouca interferência e com bloco de alta condutividade térmica e grande interferência, para ambos os fluidos. / Many studies have been carried out in square enclosures with generalized Newtonian fluids with natural convection in past few years for directly applications in compact heat exchangers, cooling of electronics systems and polymeric engineering. The natural convection in square enclosures with differently heated sidewalls, filled with power-law and Bingham fluids in addition with heterogeneous medium and centered block are analyzed in this study. The heterogeneous medium are solid, square, disconnected and conducting blocks. The parameters used are the Rayleigh number in the range 104 - 107 , power index n range of 0, 6 - 1, 6, Bingham number range of 0, 5 - Bimax , being the influence of Prandtl number investigated for each fluid model. The number of blocks for heterogeneous medium are 9, 16, 36 and 64, keeping constant solid to fluid conductive ratio, κ = 1. For enclosures with centered block are used the nondimensional block size from 0, 1 to 0, 9, with solid to fluid conductive ratio in range κ = 0, 1; 1 and 10. Mathematical modeling is done by mass, momentum and energy balance equations. The solution of equations have been numerically solved in ANSYS FLUENT R software. Firstly, numerical solutions for validation with Newtonian fluids in clean enclosures are conducted, followed by enclosures with heterogeneous medium and centered block. Subsequently, numerical solutions of power-law and Bingham fluids with same enclosures configurations are conducted. The results are reported in the form of streamlines, isotherms and average Nusselt number at hot wall. In general, the heat transfer process in enclosure is governed by Rayleigh number, size and thermal conductivity of the blocks, power index n for power-law fluid and Bingham number for Bingham fluid. Both fluid models are very sensitive with Prandtl number changes. Heterogeneous medium decrease heat transfer in enclosure when affects thermal boundary layer for both fluid models. One analytical prediction was made for power-law fluid. An increase in heat transfer occurs with low thermal conductivity block and few interference and with high thermal conductivity block and great interference, for both fluids.