Analytical Investigation of the Heat-Transfer Limits of a Novel Solar Loop-Heat Pipe Employing a Mini-Channel Evaporator

This paper presents an analytical investigation of heat-transfer limits of a novel solar loop-heat pipe developed for space heating and domestic hot water use. In the loop-heat pipe, the condensate liquid returns to the evaporator via small specially designed holes, using a mini-channel evaporator. The study considered the commonly known heat-transfer limits of loop-heat pipes, namely, the viscous, sonic, entrainment, boiling and heat-transfer limits due to the two-phase pressure drop in the loop. The analysis considered the main factors that affect the limits in the mini-channel evaporator: the operating temperature, mini-channel aspect ratio, evaporator length, evaporator inclination angle, evaporator-to-condenser height difference and the dimension of the holes. It was found that the entrainment is the main governing limit of the system operation. With the specified loop design and operational conditions, the solar loop-heat pipe can achieve a heat-transport capacity of 725 W. The analytical model presented in this study can be used to optimise the heat-transfer capacity of the novel solar loop-heat pipe.


Introduction
Loop-heat pipes (LHPs) are two-phase heat-transfer (evaporation/condensation) devices that are able to transfer large amounts of heat over long distances due to a capillary or gravitational structure [1,2]. The major advantages of LHPs compared to heat pipes (HPs) are an ability to operate against gravity and a higher maximum heat-transport capacity [1]. The maximum heat-transport capacity for an LHP system is governed by different limits, including the viscous, sonic, entrainment, boiling and capillary limits [3]. The minimum value of these limits represents the maximum heat-transport capacity of the LHP system, which is the maximum amount of heat the LHP system could transfer. This maximum heat-transfer capacity represents the most significant performance characteristic of an LHP [3]. Due to the above advantages, LHPs are ideally suitable for use in solar collection systems for hot water and space heating use, which allow heat to be collected by an evaporator outside a building, and transfer it to water flowing across the heat exchanger. Therefore, in recent years, some numerical and experimental studies have been performed to study the heat-transfer limit and thermal performance of conventional LHPs applied to solar collectors [4][5][6][7]. Existing solar LHPs are not technically mature; there are still opportunities to enhance their maximum heat-transport capacity and thermal performance [3]. An efficient way to enhance the thermal performance of solar LHPs could be the use of mini-channel heat pipes (MCHPs) in the evaporator.

LHP System Description and Operation
The loop-heat pipe studied is described in Figure 1. The system is filled with the refrigerant R134a instead of water, which is chemically incompatible with the aluminium mini-channel. In fact, the water interaction with the aluminium can lead to a corrosion reaction without natural inhibition. This corrosion can generate non-condensable gas and then causes the breakage of the internal vacuum, which causes the non-functioning of the heat pipe [15]. The system is composed of the following main elements (Figure 1): the evaporator mini-channel, the condenser (shell and tube heat exchanger), the vapour transportation line (vtl), the liquid transportation line (ltl), the vapour header (vh), the liquid header (lh) and the liquid collector at the bottom of the evaporator. The fluid evaporated in the mini-channel is collected by the vapour header and transferred towards the heat exchanger via the vapour transportation line. After condensation at the heat exchanger, the liquid returns to the water header and enters the evaporator through small holes (Figures 1 and 2). The natural fluid circulation in the loop-heat pipe system is governed by the pressure head between the condenser and the bottom of the evaporator. The solar collector is composed of 20 mini-channel heat pipes. Each mini-channel is composed of 10 ports (Figure 2). The condensed liquid enters into each port via four small holes with 0.75 mm diameter. The holes were placed on one side of each mini-channel port ( Figure 2). They do not exist at the opposite side of the wall to avoid blocking the vapour flow because of a higher gravitational effect when the evaporator is inclined. Holes for condensate liquid return into the mini-channel   Table 1 summarises the first set of the LHP system components. These values will be used to experimentally study the thermal performance of the system.    Table 1 summarises the first set of the LHP system components. These values will be used to experimentally study the thermal performance of the system.  Table 1 summarises the first set of the LHP system components. These values will be used to experimentally study the thermal performance of the system.

Mathematical Equations for Heat-Transfer Limits
Five main heat-transfer limits that affect the system operation are considered, namely, the viscous limit Q VL , sonic limit Q SL , entrainment limit Q EL , boiling limit Q BL and pressure-drop limit Q PL . These limits depend on the working fluid, the heat pipe dimensions, and the heat pipe operating temperature t v . The minimum value of these limits determines the maximum heat-transfer capacity of the LHP system. These limits are defined as following:

Viscous Limit Q VL
The viscous limit occurs at low operating temperatures, when the viscous forces are larger than the pressure gradients [1]. This results in an insufficient pressure available to drive the vapour, and leads to the evaporator dry-out. The viscous limit at the evaporator can be expressed [16,17]: By considering the liquid thickness in the micro-channel port, the vapour area A v can be expressed as the following: where 'a' is the mini-channel port width, 'b' is the mini-channel port height and δ lf is the liquid film thickness. This liquid film is formed by the condensate liquid from the four holes in the mini-channel ( Figure 2). According to Imura et al. [18], the average film thickness with a current flow may be approximated to that without one because the average film thickness becomes a little larger. The liquid film thickness δ lf can be approximated as follows: and δ lf = 0.369 3µ 2 l ρ 2 l g 1/2 Re 1/2 for Re > 400, where Re is the Reynolds number, µ l the liquid dynamic viscosity, g the gravitational acceleration. The Reynolds number is given by: U l is the superficial velocity of the liquid film. Since the liquid flow through the holes of the microchannel depends of the driving force of the LHP, which is the pressure head, it is suitable to link it with the superficial flow. To consider this effect, the superficial velocity is linked to the pressure head by the following equation (it is simply expressed as the liquid flow throughout a hole caused by the pressure head): where A h is the hole's section, H hp−he is the pressure head (height difference between the top of the evaporator and the heat exchanger), d h is the hole's diameter and D lh is the diameter of liquid header. The mass velocity U l is assumed the same for the four holes, and then the film thickness is assumed the same along the adjacent wall. A l is the liquid film section expressed as following: where 'a' is the mini-channel port width. Cd [-], the discharge coefficient of the flow from the liquid head to the hole, is expressed as following [19]: The liquid thickness is calculated iteratively by assuming that the liquid mass flow from the liquid header is equal to the mass flow of the liquid film. Equation (1) can be applied to the other components of the loop, for example, the vapour transportation line Q v, vtl , the vapour header Q v, vh and the heat exchanger Q v, he . Thus, the minimum among these items would be the viscous limit of the system. Q VL = min Q VL, vtl , Q VL, vh , Q VL, he , Q VL, e (9)

Sonic Limit Q SL
The sonic limit is due to the fact that at low vapour densities the corresponding vapour flow rate in the heat pipe may result in very high vapour speed, and a choked vapour flow can occur in the vapour passage [1,16]. The sonic limit at the evaporator can be expressed as: Equation (10) can be applied to the other components of the loop, for example, the vapour transportation line Q SL, vtl , the vapour header Q SL, vh and the heat exchanger Q SL, he . Thus, the minumum among these items would be the sonic limit of the system. Q SL = min Q SL, vtl , Q SL, vh, Q SL, he , Q SL, e

Entrainment Limit Q EL
The entrainment limit refers to the case of high shear forces that are developed when the vapour passes over the liquid in the counter-flow direction, and liquid droplets may be entrained with Energies 2018, 11, 148 6 of 18 the vapour towards the condenser. The entrainment limit exists in two parts of the system, in the evaporator and in the condenser. The entrainment limit for a wickless heat pipe can be expressed as follows [17]: where f 1 ( Figure A1), f 2 and f 3 ( Figure A2) are variables given in [16] and presented in Appendix A.
For the heat exchanger, the entrainment limit is given by: where D he,i is the heat-exchanger's internal diameter, δ w is the liquid thickness (Appendix C (A11)) and D w is the condensed liquid hydraulic diamater. The smaller of Q EL,hp and Q EL,he represents the entrainment limit of the system, as shown:

Boiling Limit Q BL
The boiling limit occurs when the applied evaporator heat flux is sufficient to cause nucleate boiling. This creates vapour bubbles that partially block the liquid return and can lead to evaporator dry-out. The boiling limit can occur at the evaporator and the condenser. The boiling limit for the evaporator can be expressed as follows [20]: where Ku = 0.14 1 − Tv The boiling limit for the heat exchanger can be expressed as follows [13,14]: where λ w is the thermal conductivity of the water and r b is the critical radius of bubble generation and is assumed to be 2.54 × 10 −7 m [14]. The smaller of Q BL,he and Q BL,he represents the boiling limit of the heat pipe system, as shown:

Pressure-Drop Limit Q PL
The heat-transfer limit due to the pressure-drop limit Q PL is reached when the condensate return level reaches the end of the condenser, and then any further increase of the heat-transfer rate causes the condensate to block part of the condenser, increasing the overall thermal resistance of the heat pipe [21]. The balance of the pressure drop is crucial for the loop operation; the total pressure head ∆P g must be superior to the total pressure drop ∆P t in the loop: where ∆P t = ∆P vtl,v + ∆P ltl,l + ∆P vh,v + ∆P lh,l + ∆P e + ∆P cond . ∆P vtl,v is the vapour pressure drop in the vapour transportation line and ∆P ltl,l is the liquid pressure drop in the liquid transportation lines. ∆P vh,v is the pressure drop in the vapour header, ∆P lh,l is the pressure drop in the liquid header, ∆P e is the two phase pressure drop in the evaporator and ∆P cond (Pa) is the two-phase pressure drop in the condenser. By solving Equation (19), the heat-transfer limit due to the pressure drop Q PL can be obtained. The items in Equation (19) can be addressed as follows.
The gravitational pressure ∆P g : ∆P g is the gravitational pressure caused by the height difference between the top of the heat exchanger and the bottom of the mini-channel evaporator, and can be expressed as [13,14]: where H he is the height of the heat exchanger and H he,hp is the height difference between the outlet of the heat exchanger and the top of the evaporator. L hp is the evaporator length and β is the evaporator inclination.
The pressure drop in the evaporator ∆P e : The two-phase pressure drop in the mini-channel evaporator is expressed as follows [22]: where ∆P tp,F is the frictional pressure drop, ∆P tp,G (Pa) is the gravitational and ∆P tp,A (Pa) accelerational pressure drops as given in Appendix B.
The pressure drops in the vapour transportation line ∆P vtl,v and the vapour header ∆P vh,v : The vapour pressure drop in the vapour transportation line and the vapour header could be expressed as [23,24]: where Q PL is the pressure-drop limit, estimated at~3.6. f vtl,v , and C vtl are given in Appendix C (A10). The vapour pressure drops in the vapour header ∆P vh,v can be expressed by Equation (22) similarly.
Two-phase pressure drop in the condenser ∆P cond : The Muller-Steinhagen and Heck correlation [25] is used to estimate the two-phase pressure drop in the condenser. The correlation is given as: where Λ and dP dz LO are given in Appendix D.
Total vapour pressure drop ∆P v : Total vapour pressure drop could be written as follows: Liquid pressure drop in the liquid header ∆P lh : The liquid pressure drop in liquid header can be expressed by the following equation: Liquid pressure drop in the liquid transportation line ∆P ltl,l can be expressed by the same form as Equation (25).

Algorithm for the Computer Model
The modelling approach used to calculate the heat-transfer limits of the LHP has been validated experimentally for a mini-channel heat pipe [26], and the mini-channel evaporator model used has been validated by much experimental data in the literature [22]. The algorithm used to compute the five heat-transfer limits is described as follows. Figure 3 presents the chart of the computer model.

i.
Given the geometry of the solar LHP, technical parameters are presented in Table 1. ii.
Given a certain operating temperature, the thermodynamic properties of the refrigerant are estimated. iii.
Calculating the viscous limits at an appropriate region in the operation, and taking the minimum value as the viscous limit by Equations (1) and (9). iv.
Calculating the sonic limit at an appropriate region in the operation, and taking the minimum value as the sonic limit by Equations (10) and (11). v.
Calculating the entrainment limits at an appropriate region in the operation, and taking the minimum value as the entrainment limit by Equations (12)- (14). vi.
Calculating the boiling limits at an appropriate region in the operation, and taking the minimum value as the boiling limit by Equations (15), (17) and (18). vii.
Running a numerical iteration to calculate the pressure-drop limit.
(a) Given the initial value of Q PL . (b) Gravity pressure is estimated ∆P g by Equation (20). (c) Total pressure drop in the loop is estimated ∆P t by Equation (19).
.5% (error allowance), heat balance is achieved and real value of Q PL can be obtained. viii.
Taking the minimum value, the five limits as the governing limit of the system operation. ix.
Program stops. Taking the minimum value, the five limits as the governing limit of the system operation. ix.
Program stops.

The Impact of the Operating Temperature of the Mini-Channel Evaporator
By using the computer model, the influence of the operating temperature on the different heat-transfer limits was investigated. Table 2 presents the variation of the five heat-transfer limits (viscous, sonic, entrainment, boiling and pressure). For each given operating temperature, the governing heat-transfer limit (in bold in the table) is the minimum value of the five heat-transfer limits. At a given operating temperature, this minimum value represents the maximum heat-transport capacity of the system.

The Impact of the Operating Temperature of the Mini-Channel Evaporator
By using the computer model, the influence of the operating temperature on the different heat-transfer limits was investigated. Table 2 presents the variation of the five heat-transfer limits Energies 2018, 11, 148 9 of 18 (viscous, sonic, entrainment, boiling and pressure). For each given operating temperature, the governing heat-transfer limit (in bold in the table) is the minimum value of the five heat-transfer limits. At a given operating temperature, this minimum value represents the maximum heat-transport capacity of the system.  Table 2 shows that the viscous and sonic limits (Q VL and Q SL ) increase with the operating temperature t v ; however, the entrainment limit, the boiling limit and the pressure-drop limit decrease when the temperature is increased. For each operating temperature, the governing heat-transfer limit is the entrainment limit. The viscous limit Q VL occurs in the vapour transportation line. The viscous limit reflects the difficulty of the vapour flow because of higher viscous forces, and is related to the Reynolds number. When the temperature increases, the Reynolds number, which is related to the working fluid thermodynamic properties, rises. Therefore, the vapour pressure difference becomes larger to overcome the viscous forces, and thus facilitates the vapour flow, and consequently the viscous limit increases. The sonic limit Q SL also appears in the evaporator. The sonic limit increases with the vapour flow rate. Increasing the operating temperature would cause a smaller vapour density and a higher vapour flow rate. Therefore, the sonic limit of the system would be higher. The sonic limit increases and tends to a limit where any increase in temperature does not affect the vapor speed and tends to be constant. Otherwise, the increased vapour flow, because of the temperature rise, increases the possibility to entrain droplets from the evaporator towards the vapour transportation line, and thus the entrainment limit Q EL that occurs in the mini-channel evaporator becomes lower. As with the entrainment limit, the boiling limit decreases when there is an increase in the operating temperature, because of the increase of the evaporation rate that can lead to the dry-out. This increased operating temperature increases the total pressure drop in the loop ∆P t because of higher flow velocity, and as a result, the heat-transfer limit Q PL due to the pressure drop decreases. Table 3 presents the effect of the aspect ratio on heat-transfer limits when the operating temperature t v is held at 50 • C and the height from the top of the absorber to the bottom of the condenser H hp,he is 0.6 m. The aspect ratio (b/a) is the ratio of the mini-channel port height 'b' by its width 'a'. Table 3 shows that all the heat-transfer limits increase with the aspect ratio. The increase of the aspect ratio by increasing 'b' on one hand creates higher vapour transporting space A v and thus higher evaporation rate, and on the other hand, leads to higher Bond number Bo (gravitational forces become more important compared to tension forces); as a result, the entrainment of the liquid droplets decreases and consequently the entrainment limit Q EL increases. Low aspect ratios lead to decreased entrainment limits because the higher vapour flow, the small vapour space and the lower importance of gravitational forces (lower Bond number Bo) favour the entrainment of liquid droplets towards the condenser. As the aspect ratio increases, the vapour transportation space and the importance of gravitational forces (higher Bond number Bo) increase, enabling more liquid to return along the heat pipe, which causes higher boiling limitation Q BL . As the vapour space increases, the shear forces decrease and the entrainment decreases slightly, tending to a constant value. Increasing the heat pipe ratio leads to a higher Reynolds number and thus lower viscous forces, and therefore, the LHP viscous limit Q VL increases. The increase of the aspect ratio causes lower vapour speed, and as a result, the sonic limit increases. The sonic limit Q SL first appears at the evaporator; as the vapour spaces increase, the sonic limit appears in the vapour transportation line and becomes constant. The fact that the aspect ratio increases allows more space for the vapour flow and a lower pressure drop, meaning that there is a higher heat-transfer limit Q PL due to the pressure drop.  Table 4 presents the effect of condenser-to-absorber height difference H hp-he ( Figure 1) on heat-transfer limits when the operating temperature t v is held at 50 • C and the aspect ratio at 0.58 (0.001/0.0017). It shows that the entrainment is the governing heat-transfer limit and it increases with an increase of the condenser-to-absorber height difference. The viscosity and sonic limits decrease with the height difference because of the reduction of the vapour space due the increase of the liquid film thickness δ lf . Figure 4 presents the evolution of the liquid thickness in function of the height H hp,he .

The Impact of the Condenser-to-Absorber Height Difference
The pressure limit remained unchanged because the impact height variation on the pressure drop was not considered.  Table 4 presents the effect of condenser-to-absorber height difference Hhp-he ( Figure 1) on heat-transfer limits when the operating temperature tv is held at 50 °C and the aspect ratio at 0.58 (0.001/0.0017). It shows that the entrainment is the governing heat-transfer limit and it increases with an increase of the condenser-to-absorber height difference. The viscosity and sonic limits decrease with the height difference because of the reduction of the vapour space due the increase of the liquid film thickness δ lf . Figure 4 presents the evolution of the liquid thickness in function of the height Hhp,he. The pressure limit remained unchanged because the impact height variation on the pressure drop was not considered.  An approximate linear relation is found and the increase of the height difference leads to an increase of the maximum heat-transfer capacity ( Figure 5). This effect can be explained by the influence

Impact of the Individual Hole's Dimension
For the loop-pipe system pipe with a 0.58 aspect ratio, an evaporator-to-condenser height difference Hhp,he of 0.6 m and the operating temperature of 50 °C, the impact of the individual hole's diameter on the heat-transfer limits was investigated. As shown in Table 5, the hole diameter ( Figure  2) has an influence on viscous limit, sonic limit and the entrainment limit. Increasing the hole's diameter increases the liquid film thickness in the micro-channel and increased the possibility to entrain the liquid towards the condenser, and so the entrainment limit becomes lower. This increase of liquid film thickness decreases the vapour space and then causes the decrease of the sonic and viscous limits.

Impact of the Evaporator Inclination Angle
While maintaining an aspect ratio (b/a) of 0.58, an operational temperature tv of 50 °C, an evaporator length Lhp of 1.9 m, an evaporator-to-condenser height difference Hhp,he at 0.6 m and keeping the other LHP parameters the same, a change of the evaporator inclination angle β would influence the heat-transfer limits of the LHP. Table 6 shows that the entrainment limit is the governing limit, and increases with the evaporator inclination angle, while the other three limits all remain the same expect for the boiling limit, which increases and attains a maximum at 60° and then decreases. The entrainment limit increases between 0° to 60°, attains a maximum at 60° and afterwards decreases between 60-90° ( Figure 6). The behaviour of this curve reflects the evolution of different parameters f1 and f3 of Equation (12) with the Bond number (Appendix A). Equation (12), used to estimate the entrainment limit of the wickless heat pipe (thermosyphon), is empirical, and many experimental studies show that the maximal heat-transport capacity occurs when the wickless heat pipe is at 60-80° to the horizontal, not vertical [17].

Impact of the Individual Hole's Dimension
For the loop-pipe system pipe with a 0.58 aspect ratio, an evaporator-to-condenser height difference H hp,he of 0.6 m and the operating temperature of 50 • C, the impact of the individual hole's diameter on the heat-transfer limits was investigated. As shown in Table 5, the hole diameter ( Figure 2) has an influence on viscous limit, sonic limit and the entrainment limit. Increasing the hole's diameter increases the liquid film thickness in the micro-channel and increased the possibility to entrain the liquid towards the condenser, and so the entrainment limit becomes lower. This increase of liquid film thickness decreases the vapour space and then causes the decrease of the sonic and viscous limits.

Impact of the Evaporator Inclination Angle
While maintaining an aspect ratio (b/a) of 0.58, an operational temperature t v of 50 • C, an evaporator length L hp of 1.9 m, an evaporator-to-condenser height difference H hp,he at 0.6 m and keeping the other LHP parameters the same, a change of the evaporator inclination angle β would influence the heat-transfer limits of the LHP. Table 6 shows that the entrainment limit is the governing limit, and increases with the evaporator inclination angle, while the other three limits all remain the same expect for the boiling limit, which increases and attains a maximum at 60 • and then decreases. The entrainment limit increases between 0 • to 60 • , attains a maximum at 60 • and afterwards decreases between 60-90 • (Figure 6). The behaviour of this curve reflects the evolution of different parameters f 1 and f 3 of Equation (12) with the Bond number (Appendix A). Equation (12), used to estimate the entrainment limit of the wickless heat pipe (thermosyphon), is empirical, and many experimental studies show that the maximal heat-transport capacity occurs when the wickless heat pipe is at 60-80 • to the horizontal, not vertical [17].  Figure 6. Impact of the evaporator inclination angle on the entrainment limit.
A higher inclination angle led to higher hydrostatic force, which improved the system capability in transporting the condensed liquid film flow, and thus led to increased heat flux while the other limits remained constant.

Impact of the Evaporator Length
While keeping the operational temperature, evaporator diameter, evaporator inclination angle, evaporator-to-condenser height difference and other LHP parameters constant, a change in the evaporator length can influence the heat-transfer limits of the LHP. Table 7 shows that by increasing the evaporator length, the viscous limit decreases while the boiling limit varies in the opposite trend. The viscous limit trend can be explained by the fact that the increase of the heat-pipe length increases the vapour distance, which could reduce the vapour volume flow rate and lead to a lower viscous limit. The sonic limit remains constant while the pressure limit increases slightly with evaporator length because of the pressure-drop decrease. The entrainment limit was found to be the main governing heat-transfer limit. The entrainment limit is constant because there is no interrelation (Equation (12)) between the entrainment limit and the heat-pipe length. Otherwise, the boiling limit became larger as the increased evaporator length led to an increased vapour space in Figure 6. Impact of the evaporator inclination angle on the entrainment limit.
A higher inclination angle led to higher hydrostatic force, which improved the system capability in transporting the condensed liquid film flow, and thus led to increased heat flux while the other limits remained constant.

Impact of the Evaporator Length
While keeping the operational temperature, evaporator diameter, evaporator inclination angle, evaporator-to-condenser height difference and other LHP parameters constant, a change in the evaporator length can influence the heat-transfer limits of the LHP. Table 7 shows that by increasing the evaporator length, the viscous limit decreases while the boiling limit varies in the opposite trend. The viscous limit trend can be explained by the fact that the increase of the heat-pipe length increases the vapour distance, which could reduce the vapour volume flow rate and lead to a lower viscous limit. The sonic limit remains constant while the pressure limit increases slightly with evaporator length because of the pressure-drop decrease. The entrainment limit was found to be the main governing heat-transfer limit. The entrainment limit is constant because there is no interrelation (Equation (12)) between the entrainment limit and the heat-pipe length. Otherwise, the boiling limit became larger as the increased evaporator length led to an increased vapour space in the LHP, evaporation rate of the fluid and operating temperature.

Conclusions
This paper presented a numerical investigation of the heat-transfer limits for a novel solar LHP system using a mini-channel heat-pipe evaporator. Six main limits of the LHP, namely, the entrainment, viscous, boiling, sonic and pressure limits, were considered. The influence of the main mini-channel evaporator parameters (aspect ratio, length, inclination angle), the LHP operating temperature and the condenser-to-absorber height difference on the limits were assessed. The study showed that, for all varied parameters, the entrainment limit was the main governing heat-transfer limit of the novel solar LHP because of higher vapour flow and shear forces in the mini-channel evaporator. With the considered design of the LHP system and operational conditions, the system can achieve a heat-transfer limit of 725 W. Operation of the system with higher operating temperature enables higher heat-transfer capacity of the system. Similarly, the aspect ratio should be relatively larger and the condenser-to-absorber height difference should be low to avoid a higher liquid film thickness that could decrease the maximum heat-transport capacity of the system. The inclination of the mini-channel evaporator should enhance the maximum heat-transfer capacity of the system. The increase of the mini-channel length would lead an increased maximum heat-transfer capacity of the system. Finally, the computer model presented in this study can be used to optimise the heat-transfer capacity of the novel solar LHP. where K p = P v g (ρ l −ρ v )σ 0.5 .
The factor f3 is a function of the inclination angle of the wickless heat pipe and the Bond number: Figure A2. Variation of factor f3 with the inclination angle β and the Bond number [27,28].
The Bond number for the mini-channel port can be expressed as follows: (A3) Figure A2. Variation of factor f 3 with the inclination angle β and the Bond number [27,28].
The Bond number for the mini-channel port can be expressed as follows: