1. Introduction
High-flux electronic devices are dominating the market with the rise of computational power and shrinkage of electronic circuitry. This can be seen almost everywhere, such as with computers, network switches for data centers, laser projectors and high-power LED lighting. To facilitate effective heat transport, the need for high thermal conductivity becomes a must-have as far as effective thermal management is concerned. The mostly used metals like copper or aluminum contain a thermal conductivity in the order of 200–400 W/(mK). To further increase the effective thermal conductivity, the conventional heat pipe that makes use of a two-phase transportation that can raise the equivalent thermal conductivity around 5000 W/(mK) is popularly used in electronic cooling. The maximum heat transfer rate of the heat pipe is mainly limited by the circulation limit of the wick structure. In this regard, typical heat pipes are not suitable for high-flux applications or long-distance transportation.
Unlike conventional heat pipes, the conventional pulsating heat pipe (PHP) is wickless and features long-distance transport capability, yet it can handle multiple heat sources or a heat source with a large area. PHP is bent into a close-loop serpentine shape configuration by a long continuous capillary tube, which consists of several straight pipes and bent tubes, and some working fluid. PHP can be divided into three parts: evaporating section, condensing section and adiabatic section, as shown in
Figure 1. With the dimension of the capillary tube, the working fluid in the tubes forms a train of randomly sized liquid slugs and vapor bubbles because of surface tension. When the evaporating section receives heat to heat up the working fluid and the condensing section dissipations heat to cool down the working fluid, the net pressure difference between the evaporating section and condensing section pushes the liquid slugs and vapor bubbles toward the other section and results in circulation.
With the working fluid is being randomly distributed, the PHP produces unbalanced flow resistance, so the working fluid has different flow patterns under different heating powers. When the power is lower, the net force can only overcome the flow resistance and transport the vapor/liquid slug in a small range where the working fluid oscillates intermittently. When the power is raised to a medium value, the net force can overcome more flow resistance and the working fluid oscillates to and from in a much larger range. When the power is even higher, the net force easily surpasses the flow resistance, thereby generating a comparatively large flow inertia to complete flow cycles in a unitary direction.
Since PHP has no capillary limit, it is suitable for long distance heat transport or high-flux applications. Although PHP has no capillary limit, it has the issue of dry-out in horizontal or anti-gravity conditions due to its symmetrical structure; its performance relies heavily on the gravity-assistance. In fact, PHP with a few turns is not a functional subject to horizontal arrangements [
1]. This is because the working fluid in the evaporation section cannot generate enough net force to overcome the flow resistance after evaporation, and the working fluid stops oscillating accordingly. Therefore, the lower temperature working fluid is unable to flow into the evaporation section, so that the temperature in the evaporation section rises rapidly. This state is called dry-out, and this had been clearly demonstrated in the PHP visualization experiment by Chien et al. [
2].
The working fluid of PHP has a corresponding saturated vapor pressure in the evaporation section and the condensing section. Most of this vapor pressure usually cancels each other out, and the remaining net pressure difference is used to push the working fluid to move to and fro. Asymmetrical forces are generated by randomly distributed working fluid. When the force is sufficient to overcome the flow resistance of the working fluid, the working fluid starts to oscillate in a small range. This status is called startup, and there have been similar temperature curve of startup behaviors in the single loop pulsating heat pipe that is visualized by Khandekar and Groll [
3].
PHPs are passive two-phase devices known for their manufacturing simplicity compared to other passive two-phase devices (i.e., conventional heat pipes). In order to prevent the bent pipes from collapsing on the outside and buckling inside the tube, it has a structural limitation in the radius of curvature. According to the standard JB/T5000.11 and Copper tube handbook [
4], the limitation in the radius of curvature of copper tube is about 2.5–3
Do. Such limitation causes huge heat spreading resistance between the heat spreading layer of the heating block or cooling block and PHP. Therefore, the amount of heat flux on the PHP is decreased appreciably.
To tackle the constraints of horizontal layout or inherited anti-gravity limitations, structural design is a common method, as shown in
Table 1, in which the improvements in PHP can be divided into three parts, namely dual-diameter, external components, and 3D design. Chien et al. [
2] and Yang et al. [
5] conducted visual observation and heat transfer measurements for flat-plate pulsating heat pipes having a non-uniform channel configuration. Tseng et al. [
6,
7] has developed non-uniform cross section area PHPs and double pipe designs to facilitate unbalanced forces in horizontal and anti-gravity applications, respectively. Kwon and Kim [
8] investigated the effect of a dual-diameter tube on the flow and heat transfer characteristics of single-turn pulsating heat pipes which were made of glass capillary tubes with various inner diameters.
Wan et al. [
9] and Feng et al. [
10] also designed the PHP to have a check valve integrated between two heating sections to create a one-way circulatory flow of working fluid. Taslimifar et al. [
11] used the magnetic field in the evaporator and ferrofluid as PHP working fluid. Their results show that the thermal performance of an open loop PHP in both startup and steady state conditions was improved. Vries et al. [
12] designed and tested the single turn PHP having a pair of Tesla-type valves, and successfully demonstrate a 14% decrease of thermal resistance. Pastukhov and Maydanik [
13] studied the PHP having ball-valve or capillary barrier structures and showed that the presence of ball check valves reduced orientation sensitivity and thermal resistance while increased the maximum heating load by 10–14%. Torresin et al. [
14] designed and tested a double condenser aluminum PHP; since the evaporating section is central, the thermal resistance of gravity-anti was same as the gravity-assist condition.
Mameli et al. [
15] tested the startup characteristics in microgravity and proved that the operation is thermally induced using a hybrid loop thermosyphon/pulsating heat pipe. Deng et al. [
16] studied the 3D-PHP, effects of number of turns, filling ratio and length of evaporating section and condensing section on the performance. Tseng et al. [
17] studied the 3D-PHP and used a fin-and-tube heat exchanger as a condenser and found that the thermal module can handle a supplied power up to 1 kW and a maximum heat flux of 15.6 W/cm
2. Hathaway et al. [
18] studied the Uneven-Turn 3D-PHP configuration to create an unbalanced force between the evaporating section and the condensing section, and successfully inverted the heat source.
From the forgoing literature survey, efforts stressed on resolving the problem of anti-gravity startup subject to high flux are still very rare. Hence, the objective of this study is to develop a three-dimensional configuration PHP having non-uniform cross section areas which is based on [
2,
6,
16]. It is a ring-shaped structure containing non-uniform cross section areas, the ring-shaped portions are arranged in order alongside an axis so as to form the rectangle coiled having a compact arrangement structure which helps to prevent a high heat spreading resistance in the heating block and cooling block. Through this configuration, it can offer a much higher operating heat flux and the working fluid automatically flows along a specific direction due to unbalanced flow resistance and unbalanced capillary force between two adiabatic sections caused by non-uniform cross section areas of the adiabatic sections. It will be demonstrated shortly that the proposed non-uniform 3D-PHP is able to operate under anti-gravity situation, which is functional in high-flux applications.
2. Experiment Setup and Data Processing
The schematic of the non-uniform experimental apparatus is shown in
Figure 2. The experimental setup consists of evaporator, condenser, water bath (YOTEC Instrument Corporation, Taiwan, model B20), direct current power supply, magnetic flowmeter, data acquisition module, and a non-uniform 3D-PHP. DC power supply contains a rated output voltage 60 V with accuracy of ±0.12 V; a rated output current 25 A with accuracy being ±0.075 A; and a maximum output power of 1500 W (Good Will Instrument Corporation, Taiwan, model PSU 60-25). A magnetic flowmeter with accuracy 0.35% and measurable range 0–12.8 L/min is used (Yokogawa Electric Corporation, Japan model AXF-005). Temperature acquisition module is used with a maximum measurement range of ±78 mV, sampling rate of 75 samples/s, and a 24-bit analog to digital converter (National Instrument Corporation, United States, model NI-9213). The evaporators were made of copper block (50 mm × 70 mm × 15 mm), and having six electronic heaters with a total equivalent resistance of 3 Ω. The condenser was made of copper block (50 mm × 70 mm × 15 mm) with volumetric flow rate of 1.25 L/min, and an inlet temperature of 50 °C. The objective of this study is to develop a passive solution for PHP upon high-flux applications. The test conditions mainly refer to the specifications from relevant manufacturers of next generation network switches. Note that the startup temperature is comparatively high, but it is not uncommon. In the meantime, the total input electric power is 811 W subject to an area of 67 mm × 67 mm, and the maximum temperature of system must be controlled below than 110 °C in an ambient temperature of 50 °C.
The conventional PHP has a structural limitation in the radius of curvature, and such a limitation causes the minimum distance between adjacent straight pipes to be about 2.5
Do. For uniform area heat source, the heat spreading model can be simplified to the finite isotropic rectangular flux channel mode as shown in
Figure 3. The heat-spreading resistance can be calculated by [
19].
where relative contact area
,
a is half the heating width,
c is half the heat spreading width, relative layer thickness
,
t is the heat-spreading thickness, Biot number
,
h is heat transfer coefficient, and
k is the thermal conductivity of heat-spreading layer.
This study proposes a compact arrangement PHP; the relative width ratio is close to 1, thus the heat spreading resistance is close to 0. Therefore, the compact arrangement PHP is suitable for a uniform area heat source, and the performance can be greater than a conventional PHP.
This assumes
h = 5000 W/(m
2K),
k = 400 W/(mK),
= 0.333,
= 0.1–1.0. The thermal resistance ratio can be calculated using Equation (1), and the results are shown in the
Figure 4. The abscissa is relative ratio
, which refers to the ratio of heat transport channel width to the total spreading width. The ordinate is relative thermal resistance ratio, which represents the ratio of spreading resistance to spreading resistance plus conduction resistance. For conventional PHPs, the relative width ratio
, and the compact arrangement pulsating of this study can attain a value of
. The relative thermal resistance ratio of conventional PHP is 8.1% at
. In this study, the relative ratio of the proposed PHP is
and relative thermal resistance ratio is only 0.6%, meaning a negligible influence of the spreading thermal resistance subject to the present compact arrangement of PHPs.
In this study, a ring-shaped rectangle PHP with non-uniform cross section areas was fabricated and tested, having an evaporating section, a condensing section, and two adiabatic sections, as shown in
Figure 5.
For the present PHP, the outer diameter is 4.0 mm and the inner diameter is 3.4 mm. The heat transport distance between the evaporating section and condensing section is 350 mm, and the distance between two adiabatic sections is 100 mm. The PHP contains 12 turns with a radius of curvature of 20 mm. The tubes of the evaporating section and one side of the adiabatic section were compressed to form an oval configuration with a minor diameter being 2 mm and the hydraulic diameter
Dh ≈ 2.24 mm, as shown in
Figure 6. The total inner volume is about 95,876 mm
3, and methanol is used as the working fluid. Methanol has a lower evaporation temperature and a higher
than water. The filling ratio (FR) is 38.7%.
The volumetric flowrate from thermostat water bath is 1.25 L/min with a supply temperature to maintain the cooling block at 50 °C. The DC power supply is used to regulate the heating power into the heating block ranging from 200 W to 800 W. The experiments were performed with gravity-assisted (bottom heating mode) and anti-gravity (top heating mode) conditions.
A total of 10 OMEGA T-type thermocouples were used, and the locations of the thermocouples and dimensions of the heating block and cooling block are shown in
Figure 7. The lengths of the evaporating section and condensing section of the PHP are both 50 mm, respectively. Thermocouples were denoted as
Te,1–
Te,3 on heating block surface,
Tc,1–
Tc,3 on cooling block surface,
Tae,1–
Tae,2 on the adiabatic section tube surface 30 mm above the evaporating section,
Tac,1–
Tac,2 on the adiabatic section tube surface 30 mm above the condensing section.
Tae,1 and
Tac,1 were attached on the same tube of the circle side.
Tae,2 and
Tac,2 were attached on the same tube of the oval side. The sampling rate of the temperature signal recorder was set to 1.0 s.
The thermal resistance
R of the present PHP was defined as temperature difference between the evaporating section and the condensing section divided by the input power, i.e.,:
where subscript
t of Equation (4) is the
t-th second recorded value of temperature,
refers to the time-average temperature of the last 200 recorded value of temperature,
and
refers to position-average temperature of the heating block and cooling block, respectively.
The combined standard uncertainty in measurement of the thermal resistance can be referred from ISO/IEC GUM [
20] and is calculated by following equation:
where
and
are standard uncertainties and the sensitivity coefficient of each measurement values, respectively.
Substituting Equations (2)–(6) into Equation (7) yields the combined standard uncertainty:
The uncertainty of temperature measurement considered the following items, including the uncertainty inherited from the calibration standard, the rounding error in calculation, the accuracy of the data acquisition device, and the standard deviation of time-average temperature. The calibration standard used a platinum resistance thermometer. The third-party calibration standard was the Fluke Dry-Well Calibrator. The calibration range is 20 °C to 80 °C, step by 20 °C. The calibration contains an expanded uncertainty of 0.04 °C and a coverage factor of 2, and the uncertainty inherited from the calibration standard is:
The temperature data is rounded to one decimal place and the probability distribution was uniform distribution, the uncertainty of rounding error is:
According to the data acquisition device (National Instrument Corporation, United States, model NI-9213) datasheet, the resolution is ±38 µV and accuracy is ±0.77 °C.
The measurement standard deviation of time-average temperature is:
The uncertainty of the temperature is:
The uncertainty of power supply considered the rounding error and accuracy of instrument. The uncertainty of the power supply is:
The uncertainties of each parameter are summarized in
Table 2. For the gravity-assisted condition, under the input power of 200 W, 400 W, 600 W, and 800 W, the combined expanded uncertainties of the thermal resistance at a level of confidence of 95% are ±6.6%, ±3.5%, ±2.7%, and ±2.3%, respectively. For the anti-gravity condition, under the input power of 200 W, 400 W, 600 W, and 800 W, the combined expanded uncertainties of the thermal resistance at a level of confidence of 95% are ±3.8%, ±3.0%, ±2.6%, and ±2.3%, respectively.
3. Results and Discussion
An average temperature variation vs. elapsed time under gravity-assisted condition is shown in
Figure 8. The experimental input power is increased from 200 W to 800 W with a step of 200 W, and the startup temperature of PHP at evaporator is about =80 °C (
Te,avg). Once the PHP is initiated, a pronounced drop of the
Te,avg is encountered (about 10 °C). Subsequently, the steady state prevails at about
Te,avg = 70 °C. On the other hand, the working fluid and the corresponding
Te,avg oscillates with a small amplitude when the input power is less than 400 W. However, appreciable oscillations of the liquid/vapor slug occur to transport heat amid the evaporator and condenser when the input power rises to 600 W. Therefore, the fluctuations of
Te,avg is comparatively large (~6 °C). A further rise of input power to 800 W result in an even pronounced oscillation frequency of the working fluid. However, then the fluctuation of the average temperature
Te,avg is relatively small when compared to an input power of 600 W. This is because the thermal response falls behind the flow oscillation, thereby flow oscillation with higher frequency does not lead to appreciable temperature fluctuation.
An average temperature variation vs. elapsed time under anti-gravity condition is shown in
Figure 9. The experiment input power is increased from 200 W to 800 W which is also in a step of 200 W. Note that the corresponding startup temperature (
Te,avg) is about =90 °C while the steady state is reached at about
Te,avg = 95 °C. When the input power is less than 400 W, the working fluid of the PHP reveals a lower oscillation frequency, therefore the temperature oscillation can be clearly seen in
Figure 9. Once the input power exceeds 600 W, the oscillation frequency of the working fluid also rises to yield a relatively stable
Te,avg.
The temperature variation of the adiabatic section vs. time for input power 200 W to 800 W is shown in
Figure 10. It can be observed that the temperature distributions are different. The left side of the charts show the gravity-assisted condition where the working fluid returns to the evaporator with the aid of gravity. Apparently, no appreciable temperature fluctuation is seen at a power of 200 W, but the oscillation frequency of temperature is increased with the rise of input power until
Q = 600 W, and a further rise to 800 W leads to moderate fluctuation the temperature. The right side of the charts show the anti-gravity operation where the relative relationship of temperature distribution is similar for input power 200 W to 800 W, and the oscillation frequency of temperature increases as the input power increases. This is because the present PHP features an uneven inner diameter for the tubes in the adiabatic section. The circle side of the adiabatic section is
Di = 3.4 mm, which is slightly higher than the suggested upper limit
mm (based on the properties of 80 °C) in the exiting literatures to activate PHP effectively. When the inner diameter of PHP is out of the recommend range, the flow pattern inside the tube can be referred to as the results of Khandekar and Groll [
21]; the circle side of the adiabatic section is unable to sustain the flow-circulation, as shown in the
Figure 11. In this regard, the working fluid returning to the evaporator can only be driven by capillary force and net pressure difference by the oval side, thereby forming a unitary circulation flow.
The corresponding average thermal resistance vs. input power is shown in
Figure 12. When the input power is less than 400 W, the resultant force of the working fluid is insufficient to drive the working fluid stably, thereby yielding unstable oscillation. Hence, the thermal resistance of the anti-gravity condition is higher than that of the gravity-assisted operation. When input power is greater than 600 W, sufficient driving force is able to facilitate a stable oscillation to form a unitary circulation flow. In this regard, the difference in overall thermal resistance between gravity-assisted and anti-gravity operations is rather small. In fact, at an input power of 600 W, the thermal resistance of gravity-assisted and anti-gravity conditions is 0.0883 K/W and 0.0918 K/W, respectively. When the input power is further raised to 800 W, the thermal resistance of gravity-assisted and anti-gravity conditions is 0.0761 K/W and 0.0796 K/W, respectively. In summary, the difference in thermal resistance between gravity-assisted and anti-gravity operations is less than 10% provided that a unitary circulation flow prevails. The results indicate that the proposed PHP is very suitable for high-flux applications where the inertial force of the working fluid may outweigh the influence of gravity.
To further demonstrate the superiority of the proposed PHP, performance is in terms of thermal resistance ratio (
Ranti-gravity/
Rgravity-assisted) and supplied heat flux at the evaporator. The lower thermal resistance ratio indicates that PHP contains better capability to resist the negative influence caused by gravity. The associated performance index is used to compare the proposed PHP with existing literatures [
7,
9,
10,
16,
18,
22,
23] is shown in
Figure 13. Some further details amid the existing studies and the present one are listed in
Table 3. As shown in the figure, the proposed PHP offers the highest heat flux of 22.9 W/cm
2 at a supplied power of 800 W while containing the least thermal resistance ratio (
Ranti-gravity/
Rgravity-assisted) of 1.05 when the input power is 800 W. Apparently, the proposed design is far better than the existing literatures [
7,
9,
10,
16,
18,
22,
23], especially at a high supplied power like 600 W or 800 W, yet it still outperforms most of the existing studies even at a moderate power of 400 W.
For the typical PHP design, a small diameter tube is normally incorporated into the serpentine configuration to make use of the capillary force, yet bends formed in the condenser/evaporator sections normally require carefully mechanical machining to avoid severe deformations. In this regard, bends with large curvature are often used for serpentine PHPs. Unfortunately, this results in a decline in effective contact surface area between evaporator/condenser and PHP circuitry, thereby impairing the benefit of the wickless design of PHP. In essence, despite PHP offering a high-flux transporting characteristic, the thermal module incorporating PHPs normally fail to meet the high-flux demand. The present novel thermal module design offers a 3-D configuration which can cover higher contact surface area of PHP alongside the evaporator. The maximum heat transfer of the pulsating heat pipe is related to the total filling of the working fluid.
This study uses a compact arrangement design, meaning that at the same heat transfer distance more working fluid can be filled, so it contains a higher heat transfer than conventional design. When more working fluid is filled, the higher flow resistance in the pulsating heat pipe will be induced. Therefore, optimal operating input power will increase accordingly and the thermal resistance at lower input power will also increase.