NUMERICAL SIMULATION ON THE TWO-PHASE FLOW PATTERN IN THE LOOP HEAT PIPE WITH R-134 A

Y. Park,M. Ha,S. Park
Abstract:This paper discusses the two-phase flow pattern in the loop heat pipe with R-134a. A computational fluid dynamics (CFD) study was carried out using ANSYS FLUENT. VOF model was used to simulate interface between vapor and liquid phase of R134a. A UDF was used to model evaporation and condensation mass transfer between two phases. For the simulation of increase of pressure in the loop heat pipe, the ideal gas law was considered when modelling the density of vapor. The numerically calculated temperatures in this paper and Fadhl’s calculated temperatures and experimentally measured temperatures matched very well [2]. The maximum difference between the calculated and Fadhl’s temperature data is 2.4 %. The bubble figure in the loop heat was observed with time passed in this paper. INTRODUCTION A heat pipe is a two-phase heat transfer device with a highly effective heat transfer rate through evaporating and condensing a fluid that is circulating in a sealed container [1]. The loop heat pipe uses capillary action to remove heat from a source and passively move it to a condenser. In the loop heat pipes, evaporation section and condensation section are existed. The heat enters the evaporation section by the heater and the working fluid absorbs an amount of heat proportional to the latent heat of evaporation, which is sufficient to change the fluid from liquid to vapor [2]. The vapor then flows to the condensation section, where the heat is transferred from vapor to the liquid. In the condensation section, vapor condenses and condensed liquid is then returned to the evaporator due to capillary forces. Loop heat pipe has advantages of being able to provide reliable operation over long distance and the ability to operate against gravity. Many numerical studies are focused on the development of mathematical modelling by calculating thermal resistance of the system or 1D thermal network analysis [3-7]. The computational fluid dynamic (CFD) simulation should be conducted to understand heat and mass transfer mechanism in the heat pipe deeply. A limited number of CFD numerical simulation studies are published because CFD studies regarding multiphase flow with phase change are very computationally expensive compared to single phase problem. NOMENCLATURE  [ 3 kg/m ] Density  [-] Volume fraction  [ kg/m s  ] Viscosity v  [ m/s ] Velocity h [ J/kg ] Enthalpy p [ Pa ] Pressure T [ K ] Temperature E S [ 3 W/m ] Energy source term M S [ 3 kg/m s ] Mass source term k [ 2 2 m /s ] Turbulent kinetic energy  [1/s ] Specific dissipation rate L [ J/kg ] Latent heat Subscripts v Vapor l Liquid sat Saturation There are some published studies regarding multiphase flow simulations in the thermosyphon [1-2, 8-11]. A. Alizadehdakhel et al. [9] modeled a gas/liquid two-phase flow and the simultaneous evaporation and condensation phenomena in a thermosyphon. The volume of fluid (VOF) technique was used to model the interaction between these phases. For modeling of evaporation and condensation in a thermosyphon, additional equations regarding mass and momentum source were calculated proposed by S. C. K. De Schepper [8]. The numerical temperature data were validated with measurement temperature. They concluded that CFD is a useful tool to model and explain the complex flow and heat transfer in a thermosyphon. The density considered in these papers was considered as constant value. However, in case of high pressure problem due to a large amount of heating power, pressure is increased by the supply of heating power continuously and as a result of increase of pressure in the system, the vapor in the system would be compressed. Therefore, the density of vapor would be increased as time passed. There are very limited number of papers considering the increase of vapor density during calculations [12]. In this study, a CFD study was conducted to analyze the effect of the location of heater inserted in the loop heat pipe. For the CFD study, the commercial code ANSYS FLUENT 15.0 was 13th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics
Engineering,Physics
What problem does this paper attempt to address?