The residential, commercial and office building sector is one of the largest consumers of energy in the world. Yet the current challenge is to reduce the energy consumption and cost caused by these buildings. However, in cities, buildings consume 70% of primary energy . In hot and humid countries, the ambient temperature and relative humidity can reach 41˚C and 84% respectively. In these tropical countries, air conditioning becomes essential to maintain thermal comfort conditions for humans in buildings. These indoor conditions should be maintained around 20˚C - 25˚C with a relative humidity of 50% - 55% according to ASHRAE (2013) and ISO 7730-2005 ; 23˚C - 26˚C in summer and 20˚C - 23.5˚C in winter, according to other standards . Air-conditioning and ventilation systems are among the largest energy consumers in countries or regions with relatively high air temperatures. Air conditioning equipment can consume up to 50% - 80% of the total electricity consumption of residential, public and commercial buildings in both developed and developing countries . In developing countries such as Guinea, the cost of electricity for the consumer from recently completed hydroelectric dams and conventional fossil fuel power plants continues to rise due to high transmission and distribution losses. In order to address this energy crisis and to improve the comfort of the population, it is urgent and important to develop a passive air-conditioning system using a two-phase thermosyphon loop for these buildings in order to reduce their consumption, greenhouse gas emissions and to fight global warming.
So far, the two-phase closed thermosyphon loop is less used in passive air conditioning in the home, but is widely applied in the cooling of electronic equipment. This loop is a natural two-phase flow cycle that can achieve efficient and low-cost heat transfer .
Numerous studies on closed-loop two-phase thermosyphon have been conducted in the literature.
Cao et al.  carried out a study on Variable Conductance Loop Thermosiphon (VCLT) active thermal control for thermal comfort in buildings based on steady state heat transfer modelling. They showed that under normal conditions, opening the control valve from 1 to 0.77 results in a decrease in heat transfer rate from 468.5 W to 71.9 W and an increase in internal flow resistance from 0.002 K/W to 0.305 K/W.
Samba et al.  carried out an experimental study of a two-phase thermosiphon loop using a telecommunication cabinet to cool telecommunication equipment with n-pentane as the working fluid. The results show that the optimum filling ratio is 9.2%.
Han et al.  developed a thermosyphon loop air conditioner combined with a vapour compression refrigerator for cooling a mobile phone base station. They showed that the COP of the thermosyphon loop reaches 9.43 when the outdoor temperature is 20˚C lower than the indoor temperature. Field tests also show that the thermosyphon saves about 34.3% to 36.9% energy compared to the traditional air conditioner.
Zhu et al. , for space heating, proposed a two-phase thermosyphon loop (TPTL) integrated in a building envelope to replace a pump-driven heat injection process system. The results showed a critical impact of the fill rate on the thermal resistance and on the behaviour of the thermosyphon loop during start-up. They also revealed that the fill rate of about 116% is not a fixed value and can be reduced to achieve a higher start-up speed or increased to achieve a lower thermal load resistance.
An experimental study was carried out by Tong et al. , on the self-regulation performance of a two-phase dual evaporator thermosyphon loop based on R744, used in data center air conditioning systems. The experimental analyses show that the mass flow rate of the high-power evaporator is higher than that of the low power evaporator with a very limited self-regulation capacity of the TPTL. In addition, they found that in practical applications with a parallel multi-evaporator TPTL system, a large difference between the heat transfer loads to the evaporators must be avoided.
In the research work of Ding et al. , they studied a heat pipe system compared to a traditional air conditioning system for a computer room. These results showed that the system can maintain the computer room at a temperature of 17˚C to 21˚C for the proper functioning of the computer equipment.
Khodabandeh , made an experimental study of a closed two-phase thermosiphon loop for cooling high heat flux electronic components of a radio base station. The tests were carried out with isobutane (R600a) as working fluid. It proves that in the thermosyphon loop system, the thermal resistances between the heat source/evaporator and the condenser/air dominate over the resistances between the evaporator walls and the condenser walls.
The objective of this study is to numerically model a two-phase thermosyphon loop for passive air conditioning in hot and dry tropical countries such as Guinea in the Mamou region using the finite element method. With high ambient temperatures mainly caused by solar irradiation and the effect of global warming in this country will increase the need and seasonal demand for space cooling in the future.
2. Description of the Physical Model of the System Studied
The two-phase thermosyphon loop (TPTL) is a reliable and energy-efficient heat transfer device based on fluid buoyancy.
In general, thermosyphons work by using a heat transfer fluid with a well-controlled amount injected into the liquid reservoir. The dissipation of heat from the system to be cooled leads to a phase change of the heat transfer fluid under a saturation pressure .
The air-conditioning system studied is a simple, fully self-contained design. It consists of an evaporator, a condenser, a riser and a downcomer (Figure 1). The evaporator is placed in the room to be cooled and the condenser outside the room. However, it should be noted that the evaporator and the condenser are plate heat exchangers. The condenser should be positioned at a higher altitude than the evaporator with some vertical distance . The working fluid absorbs
Figure 1. Schematic of the two-phase thermosyphon loop air conditioning system.
heat from the room and evaporates in the evaporator. The mass of vapour rises in the riser due to the buoyancy forces caused by the differences in density of the fluid in the riser and the downcomer. It condenses in the condenser, releasing its latent heat of condensation. The condensate flows by gravity into the downcomer and back into the evaporator. This cycle is an excellent heat carrier only by natural convection in closed loops and eliminates the need for a conventional pump to circulate the heat transfer fluid.
The material used in the various components of the two-phase thermosyphon loop is copper, with methanol as the working refrigerant.
● Thermo-physical properties of materials
The thermal properties of the materials in the system are assumed to be constant over time (Table 1).
3. Physical Modelling of Heat Transfer in the Air Conditioning System
3.1. Basic Equations
The basic equation governing the evolution of the temperature in the different sections of the air conditioning system over time is written in the following general form :
This Equation (1) will be applied to the various components of the two-phase thermosyphon loop.
3.2. Heat Transfer Balance at the Evaporator and Condenser
The evaporator is the main part of the loop, where the liquid boils. Heat is transferred as heat of vaporisation from the evaporator to the condenser. Figure 2 illustrates the points considered for evaluating the different temperatures.
Table 1. Thermo-physical properties of materials.
Figure 2. The different modes of heat transfer at the evaporator and condenser.
● For external walls
● For internal walls
● At the fluid
The volume flow rate from the evaporator to the riser is calculated by:
3.3. Heat Transfer Balance for Riser and Downcomer
The evaporator is connected to the condenser by pipes called the riser and the downcomer. In the riser (Figure 3(a)), the mass of steam produced at the evaporator rises through this pipe and flows to the condenser where it condenses. The condensate flows by gravity into the downcomer (Figure 3(b)) and returns to the evaporator as a liquid. Figure 3 shows the different points considered for temperature evaluation.
● Equation for the external wall
Figure 3. The different modes of heat transfer for the Riser (a) and for the downcomer (b).
● Equation for the internal wall
● Equation for methanol
For methanol in the vapour state:
For methanol in the liquid state:
3.4. Determination of Heat Transfer Coefficients
● The coefficients of conductive heat transfer through the thickness of these materials are given by the following formula:
● The convective exchange coefficient between the external walls of the condenser and the ambient air is given by the following expression :
● The natural convection heat transfer coefficients between the vertical inner walls between the air in the habitat enclosure and the vertical outer walls of the evaporator are determined using the correlation proposed by Churchill and Chu :
● The heat transfer coefficient between the refrigerant and the inner walls of the evaporator , which assumes nucleated boiling, is calculated by the following relationship:
● The natural convection transfer coefficient during condensation between the fluid and the inner walls of the condenser is calculated by :
● The radiative transfer coefficients between the external walls and the sky and between the external walls and the ground are determined using the following relationships :
4. Coefficient of Performance
➢ Coefficient of thermal performance (COPth)
The efficiency of the machine is expressed by the thermal coefficient of performance given by the relation (22):
: Is the power at the evaporator in (W), which is given by the relation:
: Amount of energy supplied to evaporate the fluid in (W) given by the relation:
5. Results and Discussion
The evaporator was validated from an experimental study by , who performed an experiment on the self-regulation performance of a two-phase thermosiphon loop with double parallel evaporators. They carried out the study with the refrigerant R744 in contrast to our system which operates with methanol as working fluid. The temperature on the surface of the evaporator is measured using Pt1000 sensors. The curve in Figure 4 shows that our results are in good qualitative agreement between our numerical simulation and the chosen experimental model. This result confirms the validity of our simulation model.
Figure 5shows the evolution of the temperature of the fluid, the evaporated
Figure 4. Validation of the evaporator.
Figure 5. Evolution of the fluid temperature, mass (b) and flow rate (a) evaporated at the evaporator.
mass and the flow rate of the vapour at the exit of the evaporator, this system works alone without coupling with a habitat. The simulation is carried out with a surface area of S = 2 m2. The ambient temperature is assumed to be constant over time. The decrease in temperature is justified by the fact that thermal heat is being extracted from the enclosure to be conditioned to the condenser, so we have a heat loss at this level. This drop in temperature shows the feasibility of our thermosiphon system on a large scale for thermal comfort in the home. We notice that the flow rate and the mass increase progressively. The presence of the mass and flow rate of the vapour at the outlet of the evaporator to the condenser shows us a natural circulation of the fluid and the entry of the vapour into the riser (Figure 5(a) and Figure 5(b). The steam produced passes through the riser before entering the condenser. The key role of the evaporator is to ensure superheating at the outlet and sub-cooling at the inlet for better air conditioning.
The condenser was validated from an experimental study by , who conducted experiments on a prototype to explore the operational characteristics of a two-phase thermosyphon loop (TPTL) for heating a building. To measure the temperature and pressure of the TPTL, four four-way connectors coupled to four Pt100 sensors were used. In Figure 6, we observe qualitatively, a comparable curve pattern between the results of our model simulation and the chosen experimental model. This pattern confirms the validity of our simulation model.
Figure 7illustrates the temporal evolution of the fluid temperature, mass and flow rate at the condenser. As the condenser is located outside the building, it exchanges heat with the surrounding environment. The simulation was carried out with a surface S = 0.5 m2 and applying a constant solar flux of 700 W/m2 on the outside walls, which are exposed to the sun. We can see that the temperature increases gradually over time until it reaches a maximum value of about 47˚C. This increase in the condenser temperature is mainly due to the release of heat
Figure 6. Condenser validation.
Figure 7. Evolution of fluid temperature, mass (b) and flow rate (a) condensed at the condenser.
when the methanol vapour condenses on the condenser walls and this is normal because we have a heat gain. We also notice that the increase in flow rate leads to an increase in evaporated mass. After the vapour enters the condenser, it releases latent heat with the exchange of heat from the surrounding environment and the liquid phase of the working heat transfer fluid is obtained at the condenser outlet. The decrease in flow rate (Figure 7(a)) at the condenser justifies an accumulation of the liquid at the downcomer before entering the evaporator to start a new cycle. These developments are comparable to practice.
Figure 8illustrates the influence of the solar flux on the evolution of the fluid temperature at the condenser. We can see that the temperature increases gradually before stabilising because the solar flux reaching the condenser walls has been increased. Indeed, the more the sun’s radiation on the walls of the condenser is important, the more the temperature increases as shown by the curves.
Figure 8. Evolution of the fluid temperature at the condenser with the variation of the solar flux.
The temperature varies according to the solar flux density, reaching a maximum of 55˚C when the flux density is 1000 W/m2 and a minimum of 39˚C when the radiated flux density is 400 W/m2. These evolutions represent the variation of the typical sunshine in the Mamou region during the dry season.
Figure 9illustrates the influence of the evaporator surface area on the temporal evolution of the temperature. We notice that the temperature decreases less with increasing surface area because the air temperature is constant and is fixed at 25˚C. From the curves we can see that the surface area has a very important influence on the temperature, it increases when the surface area is large because the air temperature is constant and fixed at 25˚C. This means that the larger the surface area, the greater the heat exchange, it reaches 12˚C when the surface is S = 0.5 m2 and about 14˚C when the surface is S = 3 m2.
Figure 10shows the influence of the evaporation temperature on the amount of cold produced and the thermal performance coefficient. It can be seen that when the evaporation temperature varies between 270 K and 284 K, the thermal coefficient of performance and the amount of cold produced increase progressively with it, until reaching maximum values of 0.84 and 123.25 W respectively.
Figure 11shows the influence of the condensing temperature on the amount of cold produced and the thermal performance coefficient. It can be seen that by keeping the evaporation temperature constant at 274 K and varying the condensation temperature from 302 K to 320 K, the thermal coefficient of performance of the air conditioning system and the amount of cold produced gradually decrease until reaching minimum values of 0.77 and 114 W respectively.
Figure 12shows the evolution of the fluid flow velocity in the vapour state in the riser and in the liquid state in the downcomer. We see a difference in velocity at the riser and downcomer due to the different density of methanol in the liquid
Figure 9. Evolution of the temperature of the fluid at the evaporator with the variation of the surface area.
Figure 10. Influence of the evaporation temperature on the thermal performance coefficient and the amount of cooling with Tcd = 302 K.
Figure 11. Influence of the condensation temperature on the thermal coefficient of performance and the amount of cooling with Tev = 275 K.
Figure 12. Fluid flow velocity in the riser and downcomer.
and vapour state. The maximum flow velocity value at the downcomer is about 2.2 × 10−3 m/s and that of the vapour at the riser is about 1.73 m/s.
Heat transfer modelling and simulation of a two-phase thermosyphon loop is performed. The simulation results of this study are presented and show that the temperature at the evaporator walls drops from 23˚C to 13˚C, which justifies that the system is capable of cooling. The mass flow rates and masses at the evaporator and condenser increase with temperature. The variation in evaporating and condensing temperature affects the performance of the system. For a constant evaporating and condensing temperature of 2˚C and 29˚C, the COP is 0.77 and 0.84 respectively. It has been found that the solar flux density and the evaporator surface area have a great influence on the air conditioning system. Therefore, to cool a room such as a house, a large evaporator surface is required.
From these results, we can see that the natural phase change two-phase thermosiphon loop has an excellent heat transfer capacity at low cost and without any electricity consumption. This system will therefore be an effective solution in the future for thermal comfort in the home, which is considered to be one of the largest energy consumers in the world. This study will be continued in order to couple with our habitat model to maintain the internal air temperature of the living space by natural convection.
● A two-phase thermosyphon loop is developed to show its feasibility in passive air-conditioning of housing;
● The thermal performance of the loop air conditioning system has been determined;
● The temperature on the walls of the evaporator drops from 25˚C to 13˚C with methanol and can be used to cool the ambient temperature of a room to be air conditioned.
COPth: Coefficient of thermal performance;
Cp: Heat capacity (J/kg·K);
DF: Absorbed solar flux density (W/m2);
d: Diameter (m);
ep: Wall thickness (m);
F: Geometric form factor;
Gn: Incident solar flux density (W/m2);
hxij: Heat transfer coefficient according to the transfer mode x (W/m2·K);
hc: Convection heat transfer coefficient (W/m2·K);
hr: Radiation heat transfer coefficient (W/m2·K);
hcond: Conduction heat transfer coefficient (W/m2·K);
hfg: Latent heat of vaporization (J/kg);
L: Length (m);
m: Mass (kg);
: Mass flow rate (kg/s);
Nu: Nusselt number;
Pu: Unit power (W/m2·K);
P: Vapour saturation pressure (Pa);
Pr: Prandtl number;
Qchaleur: Amount of energy (W);
Qf: Amount of cold (W);
Ra: Rayleigh number;
S: Area (m2);
T: Temperature (K);
t: Time (s);
U: Flow velocity (m/s);
V: Wind velocity (m/s);
w: Mass fraction;
∆x: Step length (m).
: Coefficient of absorptivity;
: Heat source or sink (W/m2);
: Coefficient of mass expansion (K−1);
: Stefan-Boltzmann constant (5.67 × 10−8 Wm2·K−4);
: Emissivity of the wall;
: Density (kg/m3);
: Thermal conductivity (W/m·K);
: Superficial tension (N/m);
: Dynamic viscosity (kg/m·s).
vs: Saturated vapour;
ld: Liquid in downcomer;
id: Inside of the downcomer;
ir: Inside of the riser.
 Huang, Y. and Li, C. (2021) Accurate Heating, Ventilation and Air Conditioning System Load Prediction for Residential Buildings Using Improved Ant Colony Optimization and Wavelet Neural Network. Journal of Building Engineering, 35, Article ID: 101972.
 Opoku, R., Mensah-Darkwa, K. and Samed Muntaka, A. (2018) Techno-Economic Analysis of a Hybrid Solar PV-Grid Powered Air-Conditioner for Daytime Office Use in Hot Humid Climates—A Case Study in Kumasi City, Ghana. Solar Energy, 165, 65-74.
 Charbonneau, J.Y., Douville, S., Schreiber, L. and Lemaire, J. (2004) Confort thermique à l’intérieur d’un établissement. Commission de la santé et de la sécurité du travail du Québec, Direction de la prévention-inspection, Direction des communications.
 Cao, J., Hong, X., Zheng, Z., Asim, M., Wang, Q. and Leung, M.K. (2020) Performance Characteristics of Variable Conductance Loop Thermosyphon for Energy-Efficient Building Thermal Control. Applied Energy, 275, Article ID: 115337.
 Samba, A., Hasna, L.G., Le Masson, S. and Nörterhäuser, D. (2013) Two-Phase Thermosyphon Loop for Cooling Outdoor Telecommunication Equipments. Applied Thermal Engineering, 50, 1351-1360.
 Han, L., Shi, W., Wang, B., Zhang, P. and Li, X. (2013) Development of an Integrated Air Conditioner with Thermosyphon and the Application in Mobile Phone Base Station. International Journal of Refrigeration, 36, 58-69.
 Zhu, L., Yang, Y., Chen, S. and Sun, Y. (2019) Thermal Performances Study on a Façade-Built-In Two-Phase Thermosyphon Loop for Passive Thermo-Activated Building System. Energy Conversion and Management, 199, Article ID: 112059.
 Tong, Z., Liu, X. and Jiang, Y. (2017) Experimental Study of the Self-Regulating Performance of an R744 Two-Phase Thermosyphon Loop. Applied Energy, 186, 1-12.
 Ding, T., He, Z.G., Hao, T. and Li, Z. (2016) Application of Separated Heat Pipe System in Data Center Cooling. Applied Thermal Engineering, 109, 207-216.
 Khodabandeh, R. (2004) Thermal Performance of a Closed Advanced Two-Phase Thermosyphon Loop for Cooling of Radio Base Stations at Different Operating Conditions. Applied Thermal Engineering, 24, 2643-2655.
 Chehade, A., Lemasson, S. and Louhahlia-Gualous, H. (2013) Etude des Performances Thermiques d’une Boucle Diphasique Fonctionnant avec un Nanofluide. 21ème Congrès Français de mécanique, Bordeaux.
 Zhang, P., Wang, B., Shi, W., Han, L. and Li, X. (2015) Modeling and Performance Analysis of a Two-Phase Thermosyphon Loop with Partially/Fully Liquid-Filled Downcomer. International Journal of Refrigeration, 58, 172-185.
 Camara, Y., Chesneau, X. and Kanté, C. (2018) Etude numérique du confort thermique dans un habitat bioclimatique en brique de terre stabilisée pour un climat type de la Guinée. Afrique Science, 14, 238-254.
 Oudrane, A., Aour, B., Zeghmati, B., Chesneau, X. and Hamouda, M. (2017) Analyse numerique de transfert thermique unidimensionnel d’une enceinte habitable desertique. Recueil de Mécanique, 2, 89-102.
 Khalifa, A.J.N. (2001) Natural Convective Heat Transfer Coefficient—A Review: I. Isolated Vertical and Horizontal Surfaces. Energy Conversion and Management, 42, 491-504.
 Naresh, Y. and Balaji, C. (2018) Thermal Performance of an Internally Finned Two Phase Closed Thermosyphon with Refrigerant R134a: A Combined Experimental and Numerical Study. International Journal of Thermal Sciences, 126, 281-293.