A presente tese envolve um estudo experimental da ebulição convectiva no interior de um dissipador de calor baseado em multi-microcanais. Resultados experimentais para perda de pressão e coeficiente de transferência de calor foram levantados para os hidrocarbonetos R600a (isobutano), R290 (propano) e R1270 (propileno), fluidos com reduzido GWP (Global Warming Potential) e ODP (Ozone Depletion Potential) nulo. O desempenho termo-hidráulico destes fluidos foi avaliado em um dissipador de calor de cobre, contendo cinquenta canais paralelos com seção transversal retangular de 123x494 µm2 , 15 mm de comprimento e área de base de 15x15 mm2. Os experimentos foram realizados para fluxos de calor de até 400 kW/m2, velocidade mássica variando entre 165 e 823 kg/m2s, graus de sub-resfriamento do líquido na entrada da seção de testes de 5, 10 e 15°C e temperaturas de saturação de 21 e 25°C. Os dados experimentais foram amplamente analisados e discutidos, focando o efeito do fluido refrigerante. Oscilações dos sinais de temperatura e pressão foram analisadas parametricamente visando caracterizar efeitos de instabilidades térmicas. Adicionalmente, realizou-se análise comparativa de desempenho dos refrigerantes baseada na 2ª Lei da Termodinâmica. Os dados para hidrocarbonetos foram comparados com resultados de trabalhos prévios para o refrigerante R134a levantados na mesma seção de testes e utilizando a mesma bancada experimental. A partir destes dados, conclui-se que os hidrocarbonetos proporcionam coeficientes de transferência de calor superiores ao R134a. Em geral, o coeficiente de transferência de calor apresenta a seguinte ordem decrescente: R290, R1270, R600a e R134a. No entanto, o R290 necessitou superaquecimentos da parede superiores ao R1270 para iniciar o processo de ebulição. O refrigerante R1270 proporcionou perdas de pressão totais inferiores aos demais fluidos segundo a seguinte ordem decrescente: R600a, R134a, R290 e R1270. O refrigerante R1270 apresentou frequências de oscilação inferiores na temperatura da câmara de saída. Baseado na análise de desempenho da 2ª Lei da Termodinâmica, conclui-se que, as irreversibilidades devido ao processo de transferência de calor foram predominantes quando comparadas àquelas devido à perda de pressão. Através desta análise também constatou-se o melhor desempenho para o refrigerante R290. / The present thesis concerns an experimental study on flow boiling inside a microchannel array. Experimental results for two-phase pressure drop and heat transfer coefficient were acquired for the hydrocarbons R600a (isobutane), R290 (propane) and R1270 (propylene). These fluids present low Global Warming Potential (GWP) and null Ozone Depletion Potential (ODP). The cooling performance of these hydrocarbons were evaluated for a copper heat sink containing fifty parallel microchannels. The microchannels are rectangular with cross section of 123x494 µm2, 15 mm length and a footprint area of 15x15 mm2. The experimental evaluation was performed in a test facility located at the Laboratory of Thermal and Fluid Engineering of School of Engineering of São Carlos, University of Sao Paulo. The experiments were performed for heat fluxes up to 400 kW/m2, mass velocities from 165 to 823 kg/m2s, degrees of liquid subcooling at the test section inlet of 5, 10 and 15°C and saturation temperatures of 21 and 25°C. The experimental data were carefully analyzed and discussed focusing on the effects of the fluid on the heat sink thermal hydraulic performance. Fluctuations in the temperature and pressure were analyzed parametrically in order to evaluate thermal instability effects. Additionally, an exergy analysis was performed to evaluate the refrigerant efficiency during convective evaporation. Subsequently, the parametric effects and performance of hydrocarbons were compared with previous results for refrigerant R134a obtained in the same test facility and under the same experimental conditions. The refrigerant R290 provided heat transfer coefficients higher than R600a and R1270. However, R290 needed a degree of wall superheating for the onset of nucleate boiling higher than R1270. Based on the exergy analysis it was concluded that, the irreversibility associated to the heat transfer process are predominant compared with the irreversibility due to the pressure drop. According to the Second Law analyses it was also concluded R290 as the fluid providing the best performance.

A presente tese envolve um estudo experimental da ebulição convectiva no interior de um dissipador de calor baseado em multi-microcanais. Resultados experimentais para perda de pressão e coeficiente de transferência de calor foram levantados para os hidrocarbonetos R600a (isobutano), R290 (propano) e R1270 (propileno), fluidos com reduzido GWP (Global Warming Potential) e ODP (Ozone Depletion Potential) nulo. O desempenho termo-hidráulico destes fluidos foi avaliado em um dissipador de calor de cobre, contendo cinquenta canais paralelos com seção transversal retangular de 123x494 µm2 , 15 mm de comprimento e área de base de 15x15 mm2. Os experimentos foram realizados para fluxos de calor de até 400 kW/m2, velocidade mássica variando entre 165 e 823 kg/m2s, graus de sub-resfriamento do líquido na entrada da seção de testes de 5, 10 e 15°C e temperaturas de saturação de 21 e 25°C. Os dados experimentais foram amplamente analisados e discutidos, focando o efeito do fluido refrigerante. Oscilações dos sinais de temperatura e pressão foram analisadas parametricamente visando caracterizar efeitos de instabilidades térmicas. Adicionalmente, realizou-se análise comparativa de desempenho dos refrigerantes baseada na 2ª Lei da Termodinâmica. Os dados para hidrocarbonetos foram comparados com resultados de trabalhos prévios para o refrigerante R134a levantados na mesma seção de testes e utilizando a mesma bancada experimental. A partir destes dados, conclui-se que os hidrocarbonetos proporcionam coeficientes de transferência de calor superiores ao R134a. Em geral, o coeficiente de transferência de calor apresenta a seguinte ordem decrescente: R290, R1270, R600a e R134a. No entanto, o R290 necessitou superaquecimentos da parede superiores ao R1270 para iniciar o processo de ebulição. O refrigerante R1270 proporcionou perdas de pressão totais inferiores aos demais fluidos segundo a seguinte ordem decrescente: R600a, R134a, R290 e R1270. O refrigerante R1270 apresentou frequências de oscilação inferiores na temperatura da câmara de saída. Baseado na análise de desempenho da 2ª Lei da Termodinâmica, conclui-se que, as irreversibilidades devido ao processo de transferência de calor foram predominantes quando comparadas àquelas devido à perda de pressão. Através desta análise também constatou-se o melhor desempenho para o refrigerante R290. / The present thesis concerns an experimental study on flow boiling inside a microchannel array. Experimental results for two-phase pressure drop and heat transfer coefficient were acquired for the hydrocarbons R600a (isobutane), R290 (propane) and R1270 (propylene). These fluids present low Global Warming Potential (GWP) and null Ozone Depletion Potential (ODP). The cooling performance of these hydrocarbons were evaluated for a copper heat sink containing fifty parallel microchannels. The microchannels are rectangular with cross section of 123x494 µm2, 15 mm length and a footprint area of 15x15 mm2. The experimental evaluation was performed in a test facility located at the Laboratory of Thermal and Fluid Engineering of School of Engineering of São Carlos, University of Sao Paulo. The experiments were performed for heat fluxes up to 400 kW/m2, mass velocities from 165 to 823 kg/m2s, degrees of liquid subcooling at the test section inlet of 5, 10 and 15°C and saturation temperatures of 21 and 25°C. The experimental data were carefully analyzed and discussed focusing on the effects of the fluid on the heat sink thermal hydraulic performance. Fluctuations in the temperature and pressure were analyzed parametrically in order to evaluate thermal instability effects. Additionally, an exergy analysis was performed to evaluate the refrigerant efficiency during convective evaporation. Subsequently, the parametric effects and performance of hydrocarbons were compared with previous results for refrigerant R134a obtained in the same test facility and under the same experimental conditions. The refrigerant R290 provided heat transfer coefficients higher than R600a and R1270. However, R290 needed a degree of wall superheating for the onset of nucleate boiling higher than R1270. Based on the exergy analysis it was concluded that, the irreversibility associated to the heat transfer process are predominant compared with the irreversibility due to the pressure drop. According to the Second Law analyses it was also concluded R290 as the fluid providing the best performance.

A presente dissertação trata de um estudo teórico-experimental sobre escoamento monofásico e bifásico em um dissipador de calor baseado em microcanais. Este tipo de dissipador de calor tem sido usado para a intensificação da troca de calor em sistemas compactos e de alto desempenho. A intensificação da troca de calor promovida pelo escoamento em microcanais é acompanhada de um incremento na perda de pressão, portanto o estudo destes dois parâmetros é essencial para o entendimento dos fenômenos relacionados e fundamental para o desenvolvimento de ferramentas de projeto para dissipadores de calor baseados em microcanais. Inicialmente, um levantamento bibliográfico extenso sobre a ebulição convectiva em microcanais de reduzido diâmetro foi realizado. Este estudo da literatura trata de critérios de transição entre micro- e macro-escala, padrões de escoamento, métodos de previsão do coeficiente de transferência de calor e perda de pressão. Atenção específica foi dada a estudos de dissipadores de calor baseados em microcanais. Com base nesta análise da literatura, uma bancada experimental foi confeccionada para que dados experimentais de transferência de calor e perda de pressão pudessem ser levantados a partir de um dissipador de calor de microcanais. O dissipador de calor fabricado para este estudo é constituído de 50 microcanais retangulares dispostos paralelamente com 15 mm de comprimento, 100 µm de largura, 500 µm de profundidade e espaçados entre si de 200 µm. Experimentos foram executados para o R134a, velocidades mássicas de 400 a 1500 kg/m²s, título de vapor máximo de 0,35 e fluxos de calor de até 310 kW/m². Como conclusão deste trabalho observa-se perda de pressão elevada em relação aos valores fornecidos pelos métodos de previsão da literatura e um coeficiente de transferência de calor próximo ao estimado pelo modelo de três zonas proposto por Thome et al. (2004). / This study presents a theoretical and experimental investigation on single and two-phase flows in a microchannel based heat sink. Multi-microchannel heat sinks are able of dissipating extremely high heat fluxes under confined conditions. Such characteristics have attracted the attention of academia and industry and actually several studies are being carried out in order to evaluate and optimize such devices. Initially, an extensive investigation of the literature concerning convective boiling in micro-scale channels was performed. This literature review covers transitional criteria between micro- and macro-scale flow boiling, two phase flow patterns, heat transfer coefficient and pressure drop during convective boiling. Special attention was given to studies concerning microchannels based heat sinks. Based on this investigation, an experimental facility was built for performing heat transfer and pressure drop measurements during single-phase flow and flow boiling in microchannel based heat sinks. For this study, a microchannel based heat sink was also manufactured. The heat sink contains 50 rectangular parallel microchannels, 15 mm long, 100 µm wide by 500 µm deep and separated by 200 µm walls. Experiments were performed for R134a, mass velocity of 400-1500 kg/m²s, maximum vapor quality of 0,35 and heat fluxes up to 310 kW/m². The database obtained in the present study was compared against pressure drop and heat transfer coefficient prediction methods from the literature. It was found that no one method is accurate in predicting heat sink pressure drop while heat transfer coefficient results were accurately predicted by the 3-zone model proposed by Thome et al. (2004).

A presente dissertação trata de um estudo teórico-experimental sobre escoamento monofásico e bifásico em um dissipador de calor baseado em microcanais. Este tipo de dissipador de calor tem sido usado para a intensificação da troca de calor em sistemas compactos e de alto desempenho. A intensificação da troca de calor promovida pelo escoamento em microcanais é acompanhada de um incremento na perda de pressão, portanto o estudo destes dois parâmetros é essencial para o entendimento dos fenômenos relacionados e fundamental para o desenvolvimento de ferramentas de projeto para dissipadores de calor baseados em microcanais. Inicialmente, um levantamento bibliográfico extenso sobre a ebulição convectiva em microcanais de reduzido diâmetro foi realizado. Este estudo da literatura trata de critérios de transição entre micro- e macro-escala, padrões de escoamento, métodos de previsão do coeficiente de transferência de calor e perda de pressão. Atenção específica foi dada a estudos de dissipadores de calor baseados em microcanais. Com base nesta análise da literatura, uma bancada experimental foi confeccionada para que dados experimentais de transferência de calor e perda de pressão pudessem ser levantados a partir de um dissipador de calor de microcanais. O dissipador de calor fabricado para este estudo é constituído de 50 microcanais retangulares dispostos paralelamente com 15 mm de comprimento, 100 µm de largura, 500 µm de profundidade e espaçados entre si de 200 µm. Experimentos foram executados para o R134a, velocidades mássicas de 400 a 1500 kg/m²s, título de vapor máximo de 0,35 e fluxos de calor de até 310 kW/m². Como conclusão deste trabalho observa-se perda de pressão elevada em relação aos valores fornecidos pelos métodos de previsão da literatura e um coeficiente de transferência de calor próximo ao estimado pelo modelo de três zonas proposto por Thome et al. (2004). / This study presents a theoretical and experimental investigation on single and two-phase flows in a microchannel based heat sink. Multi-microchannel heat sinks are able of dissipating extremely high heat fluxes under confined conditions. Such characteristics have attracted the attention of academia and industry and actually several studies are being carried out in order to evaluate and optimize such devices. Initially, an extensive investigation of the literature concerning convective boiling in micro-scale channels was performed. This literature review covers transitional criteria between micro- and macro-scale flow boiling, two phase flow patterns, heat transfer coefficient and pressure drop during convective boiling. Special attention was given to studies concerning microchannels based heat sinks. Based on this investigation, an experimental facility was built for performing heat transfer and pressure drop measurements during single-phase flow and flow boiling in microchannel based heat sinks. For this study, a microchannel based heat sink was also manufactured. The heat sink contains 50 rectangular parallel microchannels, 15 mm long, 100 µm wide by 500 µm deep and separated by 200 µm walls. Experiments were performed for R134a, mass velocity of 400-1500 kg/m²s, maximum vapor quality of 0,35 and heat fluxes up to 310 kW/m². The database obtained in the present study was compared against pressure drop and heat transfer coefficient prediction methods from the literature. It was found that no one method is accurate in predicting heat sink pressure drop while heat transfer coefficient results were accurately predicted by the 3-zone model proposed by Thome et al. (2004).

A pesquisa realizada envolveu a avaliação experimental dos efeitos do fluido e da orientação do escoamento no desempenho de um dissipador de calor baseado na ebulição convectiva em microcanais. Estes dissipadores de calor são usados como uma nova aplicação para a refrigeração dos novos dispositivos eletrônicos que geram altas taxas de calor. Efetuou-se inicialmente uma extensa pesquisa bibliográfica sobre o escoamento monofásico e a ebulição convectiva em microcanais e em multi-microcanais através da qual levantou-se os principais métodos de previsão do coeficiente de transferência de calor e da perda de pressão. Então, utilizando o aparato experimental desenvolvido durante o mestrado de Do Nascimento (2012) avaliou-se a transferência de calor e perda de pressão de um dissipador de calor baseado em multi-microcanais paralelos. O dissipador de calor avaliado possui 50 microcanais retangulares dispostos paralelamente com 15 mm de comprimento, 100 µm de largura, 500 µm de altura e espaçados de 200 µm. Ensaios experimentais foram executados para o R245fa, fluido de baixa pressão utilizado em ciclos frigoríficos de baixa pressão, e R407C, fluido de alta pressão usado para conforto térmico, temperatura de saturação de 25 e 31°C, velocidades mássicas de 400 a 1500 kg/m²s, graus de subresfriamento do líquido de 5, 10 e 15°C, título de vapor máximo de até 0,38, fluxos de calor de até 350 kW/m², e para 3 orientações diferentes do dissipador de calor, horizontal, vertical com os canais alinhados horizontalmente e vertical com escoamento ascendente. Os resultados obtidos foram parametricamente analisados e comparados com métodos da literatura. Coeficientes de transferência de calor médios de até 35 kW/m² °C foram obtidos. Resultados adquiridos para o R245fa e R407C foram inferiores aos levantados por Do Nascimento (2012) para o R134a utilizando o mesmo dissipador. O fluido R407C apresentou frequências e amplitudes de oscilações inferiores aos fluidos R134a e R245fa. Nenhum método para o coeficiente de transferência de calor e perda de pressão proporcionou previsões satisfatórias dos dados experimentais. O modelo Homogêneo com viscosidade da mistura bifásica dada por Cicchitti et al. (1960) apresentou as melhores previsões da perda de pressão, já para o coeficiente de transferência de calor, os métodos de Bertsch et al. (2009) e Liu e Winterton (1991) apresentaram as melhores previsões. O dissipador com sua base posicionada horizontalmente fornece coeficientes de transferência de calor superiores enquanto sua base na vertical e escoamento ascendente verificam-se perdas de pressão inferiores. Imagens do escoamento bifásico foram obtidas com uma câmera de alta velocidade e analisadas. / This study presents an experimental investigation on the effect of the fluid and the footprint orientation on the performance of a heat spreader based on flow boiling inside micro-scale channels. This heat spreader is used in an electronics cooling application with high-power density. Initially an extensive investigation of the literature concerning single-phase and two-phase flow inside a single microchannels and multi-microchannels was performed. In this literature review the leading predictive methods for heat transfer coefficient and pressure drop are described. The experimental study was carried out in the apparatus developed by Do Nascimento (2012). The heat sink evaluated in the present study is comprised of fifty parallel rectangular microchannels with cross-sectional dimensions of 100 µm width and of 500 µm depth, and total length of 15 mm. The fins between consecutive microchannels are 200 µm thick. Experimental tests were performed for R245fa, low-pressure fluid used in low pressure refrigeration cycles, and R407C, high-pressure fluid used for heat comfort, saturation temperature of 25 and 31°C, mass velocities from 400 to 1500 kg/m² s, degrees of subcooling of the liquid of 5, 10 and 15°C, outlet vapor quality up to 0.38, heat fluxes up to 350 kW/m², and for the following footprint heat sink orientations: horizontal, vertical with the microchannels aligned horizontally and vertical with upward flow. The results were parametrically analyzed and compared again the predictive methods from literature. Average heat transfer coefficients up to 35 kW/m² °C were obtained. The results for R134a from Do Nascimento (2012) for the same heat sink presented heat transfer coefficients higher than R245fa and R407C. The fluid R407C presented oscillation of the temperature due to thermal instability effects with lower frequency and amplitude lower than R134a, and R245fa. None predictive method provided satisfactory heat transfer coefficient and pressure drop predictions of the experimental data. The Homogeneous model with the viscosity given by Cicchitti et al. (1960) provided the best pressure drop prediction while the heat transfer coefficient was best predicted by Bertsch et al. (2009) and Liu and Winterton (1991). The horizontal orientation of the footprint provided the highest heat transfer coefficients while the vertical footprint orientation with upward flow the lowest pressure drops. Images of the two-phase flow were obtained with a high-speed camera and analyzed.

A pesquisa realizada envolveu a avaliação experimental dos efeitos do fluido e da orientação do escoamento no desempenho de um dissipador de calor baseado na ebulição convectiva em microcanais. Estes dissipadores de calor são usados como uma nova aplicação para a refrigeração dos novos dispositivos eletrônicos que geram altas taxas de calor. Efetuou-se inicialmente uma extensa pesquisa bibliográfica sobre o escoamento monofásico e a ebulição convectiva em microcanais e em multi-microcanais através da qual levantou-se os principais métodos de previsão do coeficiente de transferência de calor e da perda de pressão. Então, utilizando o aparato experimental desenvolvido durante o mestrado de Do Nascimento (2012) avaliou-se a transferência de calor e perda de pressão de um dissipador de calor baseado em multi-microcanais paralelos. O dissipador de calor avaliado possui 50 microcanais retangulares dispostos paralelamente com 15 mm de comprimento, 100 µm de largura, 500 µm de altura e espaçados de 200 µm. Ensaios experimentais foram executados para o R245fa, fluido de baixa pressão utilizado em ciclos frigoríficos de baixa pressão, e R407C, fluido de alta pressão usado para conforto térmico, temperatura de saturação de 25 e 31°C, velocidades mássicas de 400 a 1500 kg/m²s, graus de subresfriamento do líquido de 5, 10 e 15°C, título de vapor máximo de até 0,38, fluxos de calor de até 350 kW/m², e para 3 orientações diferentes do dissipador de calor, horizontal, vertical com os canais alinhados horizontalmente e vertical com escoamento ascendente. Os resultados obtidos foram parametricamente analisados e comparados com métodos da literatura. Coeficientes de transferência de calor médios de até 35 kW/m² °C foram obtidos. Resultados adquiridos para o R245fa e R407C foram inferiores aos levantados por Do Nascimento (2012) para o R134a utilizando o mesmo dissipador. O fluido R407C apresentou frequências e amplitudes de oscilações inferiores aos fluidos R134a e R245fa. Nenhum método para o coeficiente de transferência de calor e perda de pressão proporcionou previsões satisfatórias dos dados experimentais. O modelo Homogêneo com viscosidade da mistura bifásica dada por Cicchitti et al. (1960) apresentou as melhores previsões da perda de pressão, já para o coeficiente de transferência de calor, os métodos de Bertsch et al. (2009) e Liu e Winterton (1991) apresentaram as melhores previsões. O dissipador com sua base posicionada horizontalmente fornece coeficientes de transferência de calor superiores enquanto sua base na vertical e escoamento ascendente verificam-se perdas de pressão inferiores. Imagens do escoamento bifásico foram obtidas com uma câmera de alta velocidade e analisadas. / This study presents an experimental investigation on the effect of the fluid and the footprint orientation on the performance of a heat spreader based on flow boiling inside micro-scale channels. This heat spreader is used in an electronics cooling application with high-power density. Initially an extensive investigation of the literature concerning single-phase and two-phase flow inside a single microchannels and multi-microchannels was performed. In this literature review the leading predictive methods for heat transfer coefficient and pressure drop are described. The experimental study was carried out in the apparatus developed by Do Nascimento (2012). The heat sink evaluated in the present study is comprised of fifty parallel rectangular microchannels with cross-sectional dimensions of 100 µm width and of 500 µm depth, and total length of 15 mm. The fins between consecutive microchannels are 200 µm thick. Experimental tests were performed for R245fa, low-pressure fluid used in low pressure refrigeration cycles, and R407C, high-pressure fluid used for heat comfort, saturation temperature of 25 and 31°C, mass velocities from 400 to 1500 kg/m² s, degrees of subcooling of the liquid of 5, 10 and 15°C, outlet vapor quality up to 0.38, heat fluxes up to 350 kW/m², and for the following footprint heat sink orientations: horizontal, vertical with the microchannels aligned horizontally and vertical with upward flow. The results were parametrically analyzed and compared again the predictive methods from literature. Average heat transfer coefficients up to 35 kW/m² °C were obtained. The results for R134a from Do Nascimento (2012) for the same heat sink presented heat transfer coefficients higher than R245fa and R407C. The fluid R407C presented oscillation of the temperature due to thermal instability effects with lower frequency and amplitude lower than R134a, and R245fa. None predictive method provided satisfactory heat transfer coefficient and pressure drop predictions of the experimental data. The Homogeneous model with the viscosity given by Cicchitti et al. (1960) provided the best pressure drop prediction while the heat transfer coefficient was best predicted by Bertsch et al. (2009) and Liu and Winterton (1991). The horizontal orientation of the footprint provided the highest heat transfer coefficients while the vertical footprint orientation with upward flow the lowest pressure drops. Images of the two-phase flow were obtained with a high-speed camera and analyzed.

[pt] Diante dos constantes avanços da tecnologia os dispositivos eletrônicos vêm passando por um processo de miniaturização, ao mesmo tempo em que sustentam um aumento de potência. Essa tendência se mostra um desafio para seu gerenciamento térmico, uma vez que os sistemas de resfriamento típicos para eletrônicos utilizam ar como fluido de trabalho, e o seu baixo coeficiente de transferência de calor limita sua capacidade de atender às necessidades térmicas da indústria atual. Nesse sentido, o resfriamento bifásico tem sido considerado uma solução promissora para fornecer resfriamento adequado para dispositivos eletrônicos. Circuitos de termossifão bifásico combinam a tecnologia de resfriamento bifásico com sua inerente natureza passiva, já que o sistema não requer uma bomba para fornecer circulação para seu fluido de trabalho, graças às forças da gravidade e de empuxo. Um dissipador de calor de microcanais, localizado bem em cima do dispositivo eletrônico, dissipa o calor gerado. Isto o torna uma solução de baixo custo e energia. Além disso, ter um circuito de termossifão operando com um refrigerante de baixo GWP, como o R-1234yf, resulta em baixo impacto para o meio ambiente, uma vez que é um refrigerante ecologicamente correto e o sistema tem baixo ou nenhum consumo de energia. Este trabalho fornece um modelo numérico detalhado para a simulação de um circuito de termossifão bifásico, operando em condições de regime permanente. O circuito compreende um evaporador (chip e dissipador de calor de micro-aletas), um riser, um condensador refrigerado a água de tubo duplo e um downcomer. Equações fundamentais e constitutivas foram estabelecidas para cada componente. Um método numérico de diferenças finitas, 1-D para o escoamento do fluido por todos os componentes do sistema, e 2-D para a condução de calor no chip e evaporador foi empregado. O modelo foi validado com dados experimentais para o refrigerante R134a, mostrando uma discrepância em relação ao fluxo de massa em torno de 6 por cento, para quando o sistema operava sob regime dominado pela gravidade. A pressão de entrada do evaporador prevista apresentou um erro relativo máximo de 4,8 por cento quando comparada aos resultados experimentais. Além disso, a maior discrepância da temperatura do chip foi inferior a 1 grau C. Simulações foram realizadas para apresentar uma comparação de desempenho entre o R134a e seu substituto ecologicamente correto, R1234yf. Os resultados mostraram que quando o sistema operava com R134a, ele trabalhava com uma pressão de entrada no evaporador mais alta, assim como, com um fluxo de massa mais alto. Por causa disso, o R134a foi capaz de manter a temperatura do chip mais baixa do que o R1234yf. No entanto, essa diferença na temperatura do chip foi levemente inferior a 1 grau C, mostrando o R1234yf como comparável em desempenho ao R134a. Além disso, o fator de segurança da operação do sistema foi avaliado para ambos os refrigerantes, e para um fluxo de calor máximo do chip de 33,1 W/cm2, R1234yf mostrou um fator de segurança acima de 3. Isso significa que o circuito de termossifão pode operar com segurança abaixo do ponto crítico de fluxo de calor. Dada a investigação sobre a comparação de desempenho dos refrigerantes R134a e R1234yf, os resultados apontaram o R1234yf como um excelente substituto ecologicamente correto para o R134a, para operação em um circuito de termossifão bifásico. / [en] Given the constant advances in technology, electronic devices have been going through a process of miniaturization while sustaining an increase in power. This trend proves to be a challenge for thermal management since commonly electronic cooling systems are air-based, so that the low heat transfer coefficient of air limits its capacity to keep up with the thermal needs of today s industry. In this respect, two-phase cooling has been regarded as a promising solution to provide adequate cooling for electronic devices. Two-phase thermosyphon loops combine the technology of two-phase cooling with its inherent passive nature, as the system does not require a pump to provide circulation for its working fluid, thanks to gravity and buoyancy forces. A micro-channel heat sink located right on top of the electronic device dissipates the heat generated. This makes for an energy and cost-efficient solution. Moreover, having a thermosyphon loop operating with a low GWP refrigerant such as R-1234yf results in low impact for the environment since it is an environmentally friendly refrigerant, and the system has low to none energy consumption. This work provides a detailed numerical model for the simulation of a two-phase thermosyphon loop operating under steady-state conditions. The loop comprises an evaporator (chip and micro-fin heat sink), a riser, a tube-in-tube water-cooled condenser and a downcomer. Fundamental and constitutive equations were established for each component. A finite-difference method, 1-D for the flow throughout the thermoysphon s components and 2-D for the heat conduction in the evaporator and chip, was employed. The model was validated against experimental data for refrigerant R134a, showing a mass flux discrepancy of around 6 percent for when the system operated under gravity dominant regime. The predicted evaporator inlet pressure showed a maximum relative error of 4.8 percent when compared to the experimental results. Also, the chip temperature s largest discrepancy was lower than 1 C degree. Simulations were performed to present a performance comparison between R134a and its environmentally friendly substitute, R1234yf. Results showed that when the system operated with R134a, it yielded a higher evaporator inlet pressure as well as a higher mass flux. Because of that, R134a was able to keep the chip temperature lower than R1234yf. Yet, that difference in chip temperature was slightly lower than 1 C degree, showing R1234yf as comparable in performance to R134a. In addition, the safety factor of the system s operation was evaluated for both refrigerants, and for a maximum chip heat flux of 33.1 W/cm2, R1234yf showed a safety factor above 3. This means the thermosyphon loop can operate safely under the critical heat flux. Given the investigation on the performance comparison of refrigerants R134a and R1234yf, results pointed to R1234yf being a great environmentally friendly substitute for R134a for the two-phase thermosyphon loop.