Journal of Electronics Cooling and Thermal Control
Vol.4 No.1(2014), Article ID:43507,11 pages DOI:10.4236/jectc.2014.41003

Operating Characteristics of Multiple Evaporators and Multiple Condensers Loop Heat Pipe with Polytetrafluoroethylene Wicks

Sho Okutani1, Hosei Nagano1, Shun Okazaki2, Hiroyuki Ogawa2, Hiroki Nagai3

1Department of Aerospace Engineering, Nagoya University, Nagoya, Japan

2Japan Aerospace Exploration Agency, Sagamihara, Japan

3Department of Aerospace Engineering, Tohoku University, Sendai, Japan

Email: nagano@nuae.nagoya-u.ac.jp

Copyright © 2014 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 16 December 2013; revised 16 January 2014; accepted 23 January 2014

Abstract

This paper presents fabrication and testing of a multiple-evaporator and multiple-condenser loop heat pipe (MLHP) with polytetrafluoroethylene (PTFE) porous media as wicks. The MLHP has two evaporators and two condensers in a loop heat pipe in order to adapt to various changes of thermal condition in spacecraft. The PTFE porous media was used as the primary wicks to reduce heat leak from evaporators to compensation chambers. The tests were conducted under an atmospheric condition. In the tests that heat loads are applied to both evaporators, the MLHP was stably operated as with a LHP with a single evaporator and a single condenser. The relation between the sink temperature and the thermal resistance was experimentally evaluated. In the test with the heat load to one evaporator, the heat transfer from the heated evaporator to the unheated evaporator was confirmed. In the heat load switching test, in which the heat load is switched from one evaporator to another evaporator repeatedly, the MLHP could be stably operated. The loop operation with the large temperature difference between the heat sinks was also tested. From this result, the stable operation of the MLHP in the various conditions was demonstrated. It was also found that a flow regulator which prevents the uncondensed vapor from the condensers is required at the inlet of the common liquid line when one condenser has higher temperature and cannot condense the vapor in it.

Keywords:Loop Heat Pipe; Multiple Evaporators; Multiple Condensers; Thermal Control; Two-Phase Heat Transfer

1. Introduction

As sophistication and diversification of recent space mission, the thermal requirements for spacecraft are getting stricter. Under such background, a heat transport device with no electric power, named loop heat pipe (LHP) is under development. LHPs are two phase heat transfer devices which utilize the capillary forces developed in porous media and latent heat at phase change of fluid [1] [2] . LHPs have been used for thermal control of spacecraft such as GLAS, SWIFT, and MAXI [3] -[5] . However, future space mission requires much more strict thermal control. For example, planetary probes and lunar rovers are exposed to the large changes of external thermal environment. The amount of heats generated from the instruments changes with the situation. Additionally, usable electric power for thermal control is limited. Therefore more functional thermal control device is required. This research focuses on multiple-evaporator and multiple-condenser loop heat pipes (MLHPs) which have multiple evaporators and condensers in a LHP. The concept of the MLHP has been proposed and several researches on the MLHP have been conducted. Bienert et al. [6] proposed the concept of the dual evaporators LHP, and the feasibility of the dual LHP was demonstrated with titanium wicks and water as the working fluid. Ku et al. [7] [8] have developed a miniature MLHP with titanium wicks and ammonia as the working fluid, and comprehensive test programs were executed in order to apply a miniature MLHP to ST-8 program. However, the researches on the MLHP are still limited, and complex behaviors in the MLHP have not been well cleared. As a result, the technology of the MLHP has not been demonstrated in orbit. The final goal of this study is to figure out operations and behaviors of the MLHP under various conditions, and to establish the detailed model of the MLHP in order to develop a high functional MLHP. In the past, authors have developed a single LHP with a PTFE wick and demonstrated a good performance of the LHP [9] . In this paper, a small MLHP with PTFE wicks was newly designed and fabricated, and basic performance of the MLHP was experimentally evaluated in an atmospheric condition. Firstly, the basic characteristics and advantages of the MLHP for spacecraft thermal control were explained. Secondly, the specification of the MLHP which was designed and fabricated, and test conditions were explained. Next, the test results of the MLHP including the heat loads applied to both evaporators, to only one evaporator, and switched from one evaporator to another evaporator were presented. Finally, the test result when the sink temperature on the condensers was different under a constant heat load was presented.

2. Features and Principle of MLHP

Figure 1 shows a schematic view of a MLHP. The MLHP consists of two evaporators, two condensers, vapor lines, liquid lines and a working fluid. The evaporators and the condensers are connected in parallel with vapor lines and liquid lines. Each evaporator has a reservoir, so called compensation chamber (CC) and pours wicks.

The principle of the MLHP is as follows; In the case that heat loads applied to both evaporators (Figure 1(a)), the vapor generated in each evaporator flows to the vapor lines and confluents at the common vapor line, then flows into each condenser. The working fluid is condensed in the condensers, converges at the common liquid line, and flows to each evaporator. The mass flow rate through each evaporator is determined by a value of heat load applied to each evaporator. Temperature of each CC is determined by the ambient and the sink temperatures, the flow rate of the returning liquid from the liquid lines, and heat leak from the evaporator. If the total pressure drop of the loop exceeds the maximum capillary forces generated in a wick, the vapor penetrates the wick from the evaporator to the CC, then the operating temperature increases, and finally, exceeds the temperature limit.

In the case that the uneven heat loads are applied to each evaporator, for instance, the heat load to evaporator 1 (Evap. 1) is larger than that to evaporator 2 (Evap. 2), the CC at Evap. 1 side (CC1) is cooled and filled with cold returning liquid since the mass flow rate through Evap. 1 is larger than that through Evap. 2. In this case, the fluid in the CC at Evap. 2 side (CC2) is two-phase state, and the CC2 controls the operating temperature. When the heat load is applied to only one evaporator (e.g. Evap. 1) (Figure 1(b)), the flow direction of the working fluid in Evap. 2 is reverse (reverse flow) of the direction in the cases that heat loads applied to both

(a)(b)

Figure 1. Schematic of a MLHP. (a) Heat loads to both evaporators; (b) Heat load to one evaporator.        

evaporators (forward flow). The vapor generated in Evap. 1 is flows to the condensers through the vapor line. On the other hand, a part of the vapor flows into Evap. 2 through the vapor line, and the condensed in Evap. 2. The liquid flows into Evap. 1 through the liquid line. If the liquid from Evap. 2 has higher temperature than the returning liquid from the condensers, CC2 temperature is higher than that of CC1, and CC2 controls the operating temperature of the loop. If the total pressure drop of the loop exceeds the maximum capillary forces generated in the wick, however, the vapor penetrates the wick from Evap. 1 to CC1, then CC1 temperature becomes higher, and CC1 starts controlling the operating temperature of the loop.

Features and advantages of MLHPs are summarized as follows:

1) By having two evaporators, MLHPs can be used for cooling multiple heat sources or a heat source with a large thermal footprint.

2) By having two condensers, MLHPs can dissipate heats for two areas. For instance, in the case that one radiator is high-temperature and cannot condense the vapor in it, the other radiator can be used for dissipating heat.

3) Because of the heat transfers between the evaporators through the common vapor line, MLHPs can be used for keeping multiple instruments in a similar temperature. As a result, especially in a cold case, the electric power to maintain the instrument’s temperature can be saved.

3. Experiment System and Test Conditions

3.1. Experiment System

The MLHP was designed by using mathematical model for a single evaporator LHP [10] -[12] , on condition that both evaporators are heated and heat is transferred in forward direction, that is, the heat is transferred from each evaporator to the common vapor line, the condensers, and the common liquid line. To reduce the heat leak from the evaporators to the CCs and to lower the operating temperature, PTFE porous media is used as primary wicks which have lower thermal conductivity than traditional metal wicks [9] . Figure 2 illustrates 1) overview of the experimental MLHP which consisted of two parallel evaporators, two parallel condensers, a common vapor transport line and a common liquid return line, and 2) evaporator-CC section with temperature measurement points. The configurations of the MLHP are presented in Table 1. Each evaporator was made of stainless steel tubing with an outer diameter (O.D.) of 12 mm and a length of 70 mm. Each CC was made of a stainless steel tube of 36 mm in O.D. × 45 mm in length. The vapor line and liquid line, each 770 mm long, were made of stainless steel tubes with O.D.s of 3.18 mm, respectively. Each condenser was made of a stainless steel tube with an O.D. of 3.18 mm and length of 500 mm. Pure acetone (99.5%) is selected as a working fluid. The fluid inventory is determined as 60 cc. The size of the compensation chamber and the fluid inventory were decided from following procedure; Even in the cold case, that is, when no heat load is applied to the evaporator and the condenser sink temperature is at its minimum, the LHP is charged such that some liquid is left in the compensation chamber when the rest of the loop is completely flooded. The fluid inventory must satisfy the following relation under the cold case:

(1)

Table 1. Configurations of the MLHP.                                                                     

(a)(b)

Figure 2. Experimental MLHP. (a) Overview; (b) Evaporator-CC section.

where M is the mass of working fluid, and are liquid density and vapor density, respectively, of the working fluid at the minimum temperature, Vloop is the loop total volume excluding the compensation chamber, VCC is the compensation chamber volume, and β is the fraction of compensation chamber volume occupied by the liquid.

Another requirement is that even in the hot case, that is, when a maximum heat load is applied to the evaporator and the condenser sink temperature is at its maximum, some vapor space is available in the compensation chamber when the condenser is fully utilized. The fluid inventory must satisfy the following relation under the hot case:

(2)

where, are liquid density and vapor density, respectively, of the working fluid at the maximum temperature, VLL is the volume of the liquid line, Vpw is the void volume of the primary wick, Vgr is the volume of the grooves, VVL is the volume of the vapor line, and Vcond is the volume of the condenser lines, and α is the void fraction of compensation chamber volume.

In this study, the values of α and β are determined to be about 0.1. Once values of α and β are determined, the compensation chamber volume and the fluid inventory can be calculated from the above equations. Finally, the fluid charge must be checked against the upper limit allowed by the loop. An upper limit exists because the loop must be able to contain all the liquid volume at the maximum non-operating temperature in order to prevent bursting due to the hydrostatic pressure. Thus the following constraint applies:

(3)

where is the liquid density at the maximum nonoperating temperature. From Equations (1) and (2), the volume of the CC, VCC1 (= VCC2 = VCC) is determined as follows:

(4)

A 200-gram aluminum heater blocks were attached to the evaporators to simulate the instrument mass. Four cartridge heaters capable of delivering 1 W to 200 W were inserted into each thermal mass to provide heat loads to the evaporators. 60 T-type thermocouples are used for measurement of temperature for the MLHP. The circles and numbers in Figure 2 mark the temperature measurement points and thermocouple numbers. The condenser tubes are connected to aluminum plates, which were attached to copper cold plates. Each condenser cold plate was cooled convectively by a separate chiller. The two chillers allowed temperatures of the two condenser plates to vary independently. The applied power was obtained from the multiplication of the measured voltage and measured current. A data acquisition system consisting of a data logger, a personal computer, and a screen monitor was used to collect and store temperature and power data every second. Forty T type thermocouples were used to monitor the temperatures. The accuracies of measurement for absolute temperature was ±1˚C. Labview software was used for the command and control of the test conditions. The MLHP was set horizontally. The MLHP is insulated from environment during the tests. All tests were conducted under atmospheric conditions.

3.2. Test Conditions

Table 2 lists test conditions. Following four kinds of tests were conducted:

1) Test with heat loads to both evaporators to confirm the behavior of the MLHP. This test was conducted in the conditions that sink temperatures was varied from 10˚C to 50˚C.

2) Test with heat load to one evaporator to confirm the heat transfer between the evaporators.

3) Test to examine the response to the switching the heat load from one evaporator to the other.

4) Test with sink temperature changes under the constant heat load to investigate the operation with the large temperature difference between the sinks.

4. Test Results

4.1. Test Result of the Case 1

In the test Case 1, the heat loads applied to both evaporators were increased in 5 W increments from 0 W/0 W to 50 W/50 W in the condition that the sinks temperature was kept constant. The temperature profile during test case 1 in the sinks temperature 20˚C/20˚C is shown in Figure 3. Numbers in parentheses in Figure 3 correspond to the thermocouple numbers in Figure 2. The MLHP started 10 minutes after the heat loads of 5 W/5 W were applied, and was stably operated until 40 W/40 W was loaded. In the heat loads of 45 W/45 W, the CCs temperature began to rise, and in the heat loads of 50 W/50 W, the large temperature increasing were observed at the evaporators. This is because that the total pressure drop of the loop was over the maximum capillary forces of the wick, and the vapor penetrates the wick to the CC from the evaporator. There was the same tendency in the results of the sink temperatures from 10˚C/10˚C to 40˚C/40˚C. In the case that the sinks temperature was kept at 50˚C/50˚C, however, because the vapor was not condensed completely in the condenser, and the temperature of the returning liquid was high, the CCs and evaporators temperatures became high.

Table 2. Test conditions.

Figure 3. Temperature profile during test case 1.                       

From these results, the thermal resistance between the evaporator and the condenser Revap_cond is evaluated for each sink temperatures by Equation (5).

(5)

where, is the average temperature of the evaporator casing (TC1-4) and is the average temperature of the condenser (TC20-31).

Figure 4 shows the thermal resistances for each sink temperature. Only Evap. 1 side was shown about the argument on thermal resistance, since the same tendency was also obtained at Evap. 2 side.

It is found that as the sink temperature becomes higher, the thermal resistance becomes smaller. In the case that the sink temperature was from 10˚C to 40˚C, the thermal resistance becomes smaller as the heat loads become higher. This indicates that the operating temperature does not change comparatively with the heat loads. When the sink temperature was 40˚C, and the heat loads was 45 W/45 W, the thermal resistance takes the value of 0.58 K/W. In the case that the sink was 50˚C, and the heat loads were more than 20 W/20 W, the thermal resistance increased by rising of the operating temperature.

is divided into three parts and expressed by Equation (6).

(6)

where, Revap is the thermal resistance of the evaporator. It is evaluated from the evaporator average temperature and the vapor line (TC9) inlet temperature by Equation (7).

(7)

Rvl is the thermal resistance of the vapor line. It is evaluated from the vapor line inlet temperature (TC9) and the vapor line outlet temperature (TC16) by the following equation;

(8)

Rcond is the thermal resistance of the condenser. It is evaluated from the vapor line outlet temperature (TC16) and the condenser average temperature (TC20-25) by the following equation;

(9)

Figure 5 shows Revap, Rvl and Rcond as a function of the heat load in the case that sink temperatures were 10˚C/10˚C and 40˚C/40˚C respectively.

Figure 5(a) reveals that Rcond is the largest in the case the sinks temperatures were 10˚C/10˚C, especially large for low heat loads. This is because the most part of the condensers is used for sub-cooling the condensed liquid.

Figure 4. Thermal resistance for sink temperatures.                       

(a)(b)

Figure 5. Thermal resistances. (a) Sink 10˚C/10˚C; (b) Sink 40˚C/40˚C.     

Even when the heat loads were 50 W/50 W, the condensation finished at the first half of the condenser. For this reason, Rcond becomes larger. In the case the sinks were 40˚C/40˚C, on the other hand, Rcond becomes smaller. It is because the temperature difference between the evaporator and the condenser becomes smaller. In order to reduce the thermal resistance at lower sink temperature, it is required to shorten the condenser length or reduce thermal conductance between the condenser tube and condenser plates. For heat dissipation in space application, the thermal resistance is not necessarily important parameter since heat dissipation is conducted passively, and no electric power such as a fun is required. However, in order to keep instrument’s temperature in a narrow range by reducing subcooling, it is required to introduce thermo-electric converter or heat exchange blocks between the vapor and the liquid lines. Therefore thermal design to reduce thermal resistance as a basic characteristic is required.

Next, the thermal resistance of the evaporator Revap and the thermal resistance of the vapor line Rvl are considered. Revap and Rvl at sink temperature of 30˚C/30˚C are shown in Figure 6.

The thermal resistance of the evaporator is valued between 0.1 and 0.3 K/W. When the heat loads were less than 25 W/25 W, Revap increased with the increasing of the heat loads. However, in the heat loads more than 30 W/30 W, Revap decreased with the increasing of the heat loads and Rvl began to increase. From Figure 3, it is also found that the vapor line temperature rose when more than 30 W/30 W was applied. Thus it is considered that the situation of the heat transfer in the evaporator changes in heat loads of between 25 W/25 W and 30 W/30 W.

Figure 7 shows the temperatures of the evaporator casing (TC1-4), the vapor line inlet (TC9) and the vapor line confluence (TC13).

When the heat loads were from 10 W/10 W to 25 W/25 W, the temperature of the vapor line inlet agrees with the temperature of the vapor line confluence. When the heat loads were more than 30 W/30 W, however, the temperature difference between the vapor line inlet and the confluence exists. This means that the saturated vapor was generated in the evaporator and super-heated in the vapor grooves. For this reason, more than 30 W/30 W,

Figure 6. Thermal resistance of evaporator.                           

Figure 7. Temperatures for heat loads.                                

Revap decreases, and Rvl increases.

The heat used for super-heat is evaluated by Equation (10).

(10)

where, Tvl_co is the temperature at the confluence of the vapor line. From this equation, the value was estimated as 0.8 W in 40 W/40 W heat loads.

Based on the above-described results, the following thing is considered. For heat loads from 10 W to 25 W, the thermal resistance in the evaporator increases with the expansion of vaporization zone at the wick outside. For the heat loads of more than 30 W, the vapor super heating in the wick grooves becomes active by increasing of the thermal resistance from heating area to the evaporating zone. From this reason, Revap decreases and Rvl increases.

4.2. Test Result of the Case 2

In the test case 2, the heat load was applied to Evap. 2 only. The temperature profile during test case 2 is shown in Figure 8. The heat load of 5 W was applied to Evap. 2, and then vapor flowed to Evap. 2 side vapor line after 8 minutes. This means the startup of the MLHP. As soon as the heat load of 10 W was applied to Evap. 2, an increase of the temperature of Evap. 1, which had no heat load, was observed. This is because the vapor generated in Evap. 2 flowed into Evap. 1 through the vapor line. It is demonstrated that the heat load to Evap. 2 is usable to warm up Evap. 1.

In the heat load of 0 W/25 W, CC2 had slightly higher temperature than CC1, and controlled the operating temperature. In the heat load of 0 W/50 W, CC2 temperature began to rise. In the heat load of 0 W/55 W, CC2 temperature exceeded CC1 temperature. Accordingly, it is considered that the vapor reaches to CC2 through the wick at this moment. This means that Evap. 2 reached the capillary limit. In the heat load of 0 W/W/60 W, Evap. 2 temperature reached to 100˚C, and the test was terminated. The amount of heat transferred to Evap. 1 from Evap. 2 was calculated by Equation (11).

(11)

Figure 8. Temperature profile during test case 2.                       

where, Tvap is the vapor temperature in the evaporator—in this calculation, the vapor line inlet temperature (TC10) is used. Assuming that the heat transfer coefficient in Evap. 1 is 100 W/m2K, Qevap is estimated as 0.8 W in the heat load from 0 W/20 W to 0 W/50 W. The reason for this small Qevap is that the heat dissipation to the ambient from Evap. 1 was restricted by the thermal insulation. It is possible to increase the amount of heat transferred between the evaporators by absorbing more heat from Evap. 1.

4.3. Test Result of the Case 3

In the test case 3, the heat loads to each evaporators were switched under the condition that total heat load to the MLHP was kept at 50 W. The temperature profile during test case 3 is shown in Figure 9. At first, the heat loads of 25 W/25 W were applied and the MLHP started up. Secondly, the heat loads was changed into 50 W/0 W from 25 W/25 W. The temperature of Evap. 2 began to fall. However, it was stable at around 55˚C because the vapor flowed into Evap. 2 from Evap. 1, and the vapor flow direction change was occurred with the change of the heat loads. In this time, steady temperatures of Evap. 1 and Evap. 2 were 79.4˚C and 54.5˚C respectively. Thirdly, the heat loads was changed into 0 W/50 W from 50 W/0 W. Then the flow direction of the working fluid in the evaporators was passively changed, namely, became forward flow in Evap. 2 and reverse flow in Evap. 1.

When the heat load is applied to one evaporator, the CC of the evaporator with no heat load controls the operating temperature of the loop, and it is passively switched with changing of the heat loads. When the heat loads change to 25 W/25 W again, it is confirmed that the flow direction switches to forward.

In this way, when the heat loads was largely changed, the state of inner flow changed passively and the MLHP could be operated stably.

4.4. Test Result of the Case 4

In the test case 4, the heat loads to both evaporators were kept at 30 W/30 W, condenser 1 sink temperature was raised from −20˚C to 80˚C at hourly intervals, while condenser 2 sink temperature was kept at −20˚C. The temperature profile during test case 4 is shown in Figure 10. The steady temperature of each part of the MLHP when the sink temperatures are −20˚C/−20˚C, 40˚C/−20˚C and 80˚C/−20˚C are shown in Figure 11. When the sink temperatures were −20˚C/−20˚C, the MLHP was operated stably (Figure 11(a)). The operating temperature stayed about the same when the sink temperatures were less than +20˚C/−20˚C as shown in Figure 10. When sink1 temperature was raised to +40˚C, however, the temperatures of the both evaporators and the both condensers began to rise (Figure 11(b)). It is considered that the working fluid through condenser 1 was not condensed completely and condensed in the liquid line by meeting the sub-cooled liquid from condenser 2. When sink1 temperature was raised to +80˚C, it seems that the vapor reached the liquid line outlet, that is, the vaper was condensed from the liquid line to the CC. However the MLHP could keep operating by releasing heats to the ambient from the liquid line and the CCs (Figure 11(c)). In the case the sink temperature is high, to prevent the increasing of the operating temperature, a flow regulator which prevents the vapor to reach the CCs will be required.

Figure 9. Temperature profile during test case 3.                       

Figure 10. Temperature profile during test case 4.                      

(a)(b)(c)

Figure 11. Steady temperatures. (a) Sink: −20˚C/−20˚C; (b) Sink: +40˚C/−20˚C; (c) Sink: +80˚C/−20˚C.                

5. Conclusions

The MLHP with PTFE wicks was newly fabricated and tested. In the tests that heat loads are applied to both evaporators, the MLHP was stably operated up to 40 W/40 W. The relation between the thermal resistance and the sink temperature was also cleared. In the test with heat load to one evaporator, the MLHP was stably operated up to 0 W/50 W, and the heat transfer between the evaporators was confirmed. In the heat load switching test, the MLHP could stably be operated. The operation with the large temperature difference between the heat sinks was tested. From this result, it was shown that retuned liquid or vapor from the condenser with higher temperature increased the operating temperature of the loop. It was also found that a flow regulator which prevents the uncondensed vapor from the condensers is required at the inlet of the common liquid line.

Hereafter, solving these problems, improvements of the steady state model, thermal vacuum tests will be conducted.

References

  1. Ku, J. (1999) Operating Characteristics of Loop Heat Pipes. Proceedings of the 29th International Conference on Environmental Systems, Paper No.1999-01-2007.
  2. Maydanik, Y.F. (2005) Review Loop Heat Pipes. Applied Thermal Engineering, 25, 635-657. http://dx.doi.org/10.1016/j.applthermaleng.2004.07.010
  3. Baker, C. and Grob, E. (2001) System Accommodation of Propylene Loop Heat Pipes for the Geoscience Laser Altimeter System (GLAS) Instrument. Proceedings of the 31st International Conference on Environmental Systems, SAE Paper No. 2001-01-2263.
  4. Choi, M. (2004) Thermal Vacuum/Balance Test Results of Swift BAT with Loop Heat Pipe Thermal System. Proceedings of the 2nd IECEC, AIAA Paper No. 2004-5683.
  5. Nagai, H. and Ueno, S. (2005) Performance Evaluation of Doublecondenser Loop Heat Pipe onboard Monitor of All-Sky X-Ray Image (MAXI) in Thermal Vacuum Testing. Proceedings of the 35th International Conference on Environmental Systems, SAE Paper No. 2005-01-2939.
  6. Bienert, W.B., Wolf, D.A., Nikitkin, M.N., Maydanik, Y.F., Fershtater, Y., Vershinin, S. and Gottschlich, J.M. (1997) The Proof-Of-Feasibility of Multiple Evaporator Loop Heat Pipes. Proceedings of the Sixth European Symposium on Space Environmental Control Systems, 393-398.
  7. Ku, J., Ottenstein, L., Douglas, D., Pauken, M. and Birur, G. (2005) Miniature Loop Heat Pipe with Multiple Evaporators for Thermal Control of Small Spacecraft. Government Microcircuit Applications and Critical Technology Conference, Paper No.183.
  8. Nagano, H. and Ku, J. (2007) Capillary Limit of a Multiple-Evaporator and Multiple-Condenser Miniature Loop Heat Pipe. Journal of Thermophysics and Heat Transfer, 21, 694-701. http://dx.doi.org/10.2514/1.26151
  9. Nagano, H., Fukuyoshi, F., Ogawa, H. and Nagai, H. (2011) Development of an Experimental Small Loop Heat Pipe with Polytetrafluoroethylene Wicks. Journal of Thermophysics and Heat Transfer, 25, 547-552. http://dx.doi.org/10.2514/1.T3614
  10. Nagano, H., Fukuyoshi, F., Nagai, H. and Ogawa, H. (2008) Study on Thermal Characteristics of a Small Loop Heat Pipe. Journal of Thermal Science and Technology, 3, 355-367. http://dx.doi.org/10.1299/jtst.3.355
  11. Nishikawara, M., Nagano, H. and Kaya, T. (2013) Transient Thermo-Fluid Modeling of Loop Heat Pipes and Experimental Validation. Journal of Thermophysics and Heat Transfer, 27, 641-647. http://dx.doi.org/10.2514/1.T3888
  12. Kaya, T. and Hoang, T.T. (1999) Mathematical Modeling of Loop Heat Pipes and Experimental Validation. Journal of Thermophysics and Heat Transfer, 13, 314-320. http://dx.doi.org/10.2514/2.6461