American Journal of Mechanical Engineering
Volume 8, 2020 - Issue 1
Website: http://www.sciepub.com/journal/ajme

ISSN(Print): 2328-4102
ISSN(Online): 2328-4110

Article Versions

Export Article

Cite this article

- Normal Style
- MLA Style
- APA Style
- Chicago Style

Review Article

Open Access Peer-reviewed

Md Shahidul Haque, Antara Majumdar, Md. Faisal Kader^{ }, Md. Mahbubur Razzaq

Received March 08, 2020; Revised April 17, 2020; Accepted May 05, 2020

Pulsating Heat Pipe (PHP) has become one of the promising options for passive heat transfer. This technique is especially efficient for cooling of micro and nano sized particles. This paper works on a aluminum made and ammonia charged pulsating heat pipe. The heat pipe consists of aluminum tubing having total 14 turns of 3mm and 4mm as inner and outer diameter. A constant electric power supply of 36 W was provided to run the pulsation mechanism. The experiment was carried out for three filling ratios of 0.4, 0.6 and 0.8 as well as inclination angle of 30°, 45°, 60°, 90° and 180°. The objective of this research is to formulate an empirical correlation between overall heat transfer coefficient and different physical parameter such as heat input, evaporator temperature rise, fill ratios and inclinations. It has been found that the correlation is effective for heat input and filling ratios. Experimental results showed that for these two parameters the correlation converges well.

Closed loop pulsating heat pipe (CLPHP) is a two phase passive heat transfer device that works on continuous forming of liquid slugs and vapor plugs. The device consists of evaporator, condenser and adiabatic section. Self-excitement is created due to temperature difference between evaporator and condenser zone that circulates the annular flow throughout the tubing.

The idea of such kind of passive heat transfer device was first patented by Tesla ^{ 1} in the year 1920. Gaugler ^{ 2} developed the idea to build a heat pipe later in the year 1942. After that, industrially suitable closed loop pulsating heat pipe had been introduced by Akachi ^{ 3} in 1990. Investigations to enhance the performance of pulsating heat pipe had been carried out by changing geometrical shapes, inserting wicks, check valves, working fluid etc ^{ 4, 5, 6, 7, 8, 9, 10}. Researches for the development of mathematical modeling had been carried out to predict the two phase phenomenon of pulsating heat pipe ^{ 11, 12, 13, 14}. In this research work, linear empirical correlation has been tried to develop using physical parameters for the designated setup and working fluid.

A closed loop pulsating heat pipe (CLPHP) (inside volume of 6.264×10^{-5}m^{3}) was made by using aluminum tube of 4.0 mm Outer Diameter (OD), 3.0 mm Inner Diameter (ID) and 8.862 m length to form 14 parallel channels with 13 bends. The CLPHP was made with 3 sections: evaporator, adiabatic and condenser section. The evaporator section is about 210 mm long, the adiabatic section is 120 mm long and the condenser section is 324 mm long. The evaporator section of the heat pipe was placed inside aluminum block (0.35m×0.21m×0.035m) having grooves of same dimension of the heat pipe to avoid any gap between the block and the pipe outside surface to ensure smooth heat flow to the CLPHP. A heating coil was fabricated by winding nichrome wire of 0.25 mm diameter on ceramic bead at a constant interval of 1.50 mm was placed inside the slot of Aluminum block which acted as heat source. Heat was supplied to the working fluid directly from AC power supply with the help of a variac. At different locations on the setup, fourteen K-type thermocouples were attached by insulated tape to monitor the temperature during the experiment. A pressure gauge (maximum range: 300psi) was employed to measure pressure variations at different sections. The experimental setup is illustrated in Figure 1.

The apparatus was placed on a wooden stand to eliminate heat loss by conduction. One gate valve was connected to the heat pipe to remove air as well as to fill the working fluid. The stand can be rotated at different orientations. At the vertical orientation (0°), the evaporator section was below the condenser section which gets reversed for 180° inclination. For the tests under horizontal orientation (90°), all sections were on the same plane. Both evaporator and adiabatic sections were having thermal insulation (0.39m×0.35m×0.075m) using layers of glass wool, while the condenser section was open to the surrounding to get cooled by natural flow of air inside the room. Figure 2 depicts the different orientation of the whole setup.

The tests were conducted between 24°C - 28°C room temperature. At first, the vacuum was created inside the pipe through vacuum pump and the gate valve made closed. The vacuum pressure was monitored by using a vacuum gauge. Liquid ammonia was taken in a syringe which was attached to the gate valve carefully. Under vacuum pressure, ammonia was injected to the pipe taking three different fill ratios as 0.4 or 2.506×10-5m3 (40% of inside volume), 0.6 or 3.758×10-5m3 (60% of inside volume) and 0.8 or 5.011×10-5m3 (80% of inside volume) and the gate valve was closed immediately to resist air insertion.

The thermocouples were first calibrated using a standard thermometer (range: -40°C~1000°C). Temperatures were recorded in boiling water of 100°, melting ice of 0°C and some other temperatures with the standard thermometer and thermocouples. One selector switch was used to monitor the temperature of five and two different points on condenser and adiabatic sections respectively and another selector switch was used to monitor the temperature of seven different points of evaporator section. Knobs were used for selecting desired thermocouples to detect temperature reading on a digital thermometer. The temperature at different sections are monitored and recorded at every 20 minutes until the steady state condition. Below 90° inclination, the system reaches steady state after about 260 minutes while beyond this position it takes more time to reach steady state even though it is not stable. For all the tests, the rise of temperature was monitored and recorded for 260 min. The input electric energy was kept constant at 36 W through variac to the aluminum block. Various parameters were measured for different inclinations at 0°, 30°, 45°, 60°, 90° and 180°, three distinct fill ratios as 40%, 60% and 80%.

Electric heat input is transferred to the aluminum block through nichrome wire heat coil. By conduction heat transfer heat is transferred from aluminum block to the adjacent heat pipes. After that, heat is transferred to the inside working fluid by convection heat transfer This relationship and calculated data table shown by Haque ^{ 15} are provided below,

(1) |

Where, = V^{2}/R_{electrical} = 36 W.

Equation (1) can be rewritten as,

(2) |

Thus gives the amount of total heat loss from the AC power supply to the working fluid.

Here,

(3) |

Where, *M*_{Al}= 7.5 kg, *S*_{Al}= 900 J/kg˚C, = 20 min.

Conduction heat transfer through the PHP wall from the evaporator to the condenser section can be derived from the Fourier’s law of heat conduction,

(4) |

Here, A= nπ(d_{O}²-d_{i}²)/4, n= 14, ΔT= Δx=

Thermal conductivity of Aluminum, K=

Substituting all of these in the equation (4), we get,

(5) |

In the current test set up, n = 14, d_{o} = 0.004 m, d_{i} = 0.003 m, = 0.12 m, =250 W/m °C

Overall Heat Transfer Coefficient,

(6) |

From the analysis of experimental results of closed loop Pulsating Heat Pipe (CLPHP) discussed in this Chapter, the thermal performance of CLPHP depends on various parameters such as heat input, heat transfer coefficient, condenser inclination, evaporator temperature, ambient temperature and fill ratio.

A constant heat input of 36 W was applied in the heating block and temperature rise in various sections was monitored till steady state was reached. This is carried out for three different fill ratios, namely, 0.4, 0.6 and 0.8 and at 0°, 30°, 45°, 60°, 90° and 180° inclinations.

**Fig****ure****3.**Comparison of the rates of increase of temperature with time in the three sections of the CLPHP for fill ratio = 0.4 and inclination = 0° (Vertical position)

**Fig****ure****4.**

**Fig****ure****5.**Comparison of the rates of increase of temperature with time in the three sections of the CLPHP for fill ratio = 0.4 and inclination = 180° (Heat source above the heat sink)

It is seen from the above figures that the difference in temperature between the evaporator and the condenser sections is uniform throughout the experimental range for 0° inclination. But for 90° and 180° inclinations this difference in temperature gradually increases with time. This observation again implies that the heat transfer rate decreases with increase of inclination above 90°.

**Fig****ure****6.**Comparison of the rates of increase of temperature with time in the three sections of the CLPHP for fill ratio = 0.6 and inclination = 0° (Vertical position)

**Fig****ure****7.**Comparison of the rates of increase of temperature with time in the three sections of the CLPHP for fill ratio = 0.6 and inclination = 90° (Horizontal position)

**Fig****ure****8.**Comparison of the rates of increase of temperature with time in the three sections of the CLPHP for fill ratio = 0.6 and inclination = 180° (Heat source above the heat sink)

**Figure****9.**Comparison of the rates of increase of temperature with time in the three sections of the CLPHP for fill ratio = 0.8 and inclination = 0° (Vertical position)

**Figure****10.**Comparison of the rates of increase of temperature with time in the three sections of the CLPHP for fill ratio = 0.8 and inclination = 90° (Horizontal position)

**Figure****11.**Comparison of the rates of increase of temperature with time in the three sections of the CLPHP for fill ratio = 0.8 and inclination = 180° (Heat source above the heat sink).

It is seen from the above Figure 6 to Figure 11, that as with fill ratio =0.4, the difference in temperature between the evaporator and the condenser sections for the other two fill ratio is also uniform throughout the experimental range for 0° inclination. But for 90° and 180° inclinations this difference in temperature gradually increases with time though the rate of increase of temperature difference with 90^{o} is very small. This observation again implies that the heat transfer rate decreases with increase of inclination above 90°.

**Figure****12.**Effect of fill ratio(V/Vmax) on evaporator temperature T_{eva}(°C) in CLPHP at different inclination.

**Figure****13.**Effect of fill ratio(V/Vmax) on condenser temperature T_{ con}(°C) in CLPHP at different inclination.

It is seen from the above Figure 12 to Figure 13 that with fill ratio =0.4 and 0.6 for inclination angle below 90^{o} the evaporator temperature is qualitatively and quantitavely same but for fill ratio 0.8 it is higher.For inclination angle 90° and above evaporator section temperature is much higher for all fill ratios.The condenser sections temperature for all fill ratios is lower for inclination angle 90° and 180°. This observation implies that the heat transfer rate is lower at fill ratio 0.8 compared to fill ratio 0.4 and 0.6. It also decreases with increase of inclination above 90°. It is also evident from Figure 12 that the temperature of evaporator section for fill ratios 0.4 and 0.6 is slightly lower in inclination angle 30°.

For fixed heat input 36W, heat transfer through the PHP has been calculated using equations from chapter 3.

**Figure****14.**Effect of fill ratio(V/Vmax) on Heat transfer rate by working fluid , Q_{php}(Watt) at different inclination.

From Figure 14 above it is observed that for all fill ratios the Heat transfer rate, Q_{php} at inclinations 0°, 30°, 45°and 60° in are very close to each other though it is slightly lower for fill ratio 0.8 . It is also observed that for all fill ratios heat transfer rate, Q_{php} at inclinations 90°and 180° is much lower and it is lowest for fill ratio 0.8.So we can say that at fill ratios 0.4 and 0.6 heat transfer rate, Q_{php }is independent of inclination below 90°.

**Figure****15.**Effect of inclination, θ on heat transfer rate by working fluid, Q_{php}(Watt) of CLPHP for different fill ratio at constant heat input

It is observed from Figure 15 that for all inclination below 90° heat transfer rate by working fluid, Q_{php} for fill ratio 0.4, 0.6 and 0.8 is approximately same though Q_{php} for 0.8 is slightly lower .For inclination 90° and 180° this Q_{php }is lower compared to other inclination. From fig it is also clear that Q_{php }for inclination 180° and fill ratio 0.8 is the lowest.

In order to correlate the data obtained from the present experimental investigation, some dimensionless group such as heat input ratio, temperature difference ratio, ratio of condenser inclination and fill ratio are proposed. To determine the dependency of overall heat transfer ratio on the other dimensionless groups, the following steps are considered.** **

The present experimental investigations were carried out with different heat inputs. From the experimental results, it is found that, with the increase of heat input, the heat transfer rate is increased. For the analysis of heat transfer performance for different heat input, a ratio of heat input to the maximum heat input can be considered as

Thus, the correlation can be assumed as

(7) |

In this case, using conventional techniques, the exponent *n*_{1} was determined by plotting (U/U_{max}) vs. (Q/Q_{max}) from the present experimental data as shown in Figure 16 and the average of the slope was found *n*_{1}_{ }= 0. 7545. So the correlation became:

(8) |

In the heat pipes, the heat is transferred by the evaporation and condensation of working fluid from evaporator and condenser section respectively. So, the heat transfer can be defined as the phase change heat transfer. Another mode of heat transfer can be occurred by the convection which is directly depends to the temperature differential of the evaporator and condenser section. From the experimental date, it is evident that, with the increase of temperature differential of evaporator to condenser section, the overall heat transfer coefficient also increases. From this, a dimensionless group can be assumed as following:

This dimensionless parameter also should include the contribution of heat transfer performance of Looped Parallel Micro Heat Pipes (LPMHP). By including this dimensionless group the correlating equation became,

(9) |

In this case also using conventional techniques, the exponent *n*_{2} of the dimensionless heat transfer ratio was determined by plotting [U/U_{max}] / [(Q/Q_{max})^{0.7545}] vs. (T_{eva}-T_{a})/T_{eva}) of the present experimental data as shown in Figure 17. The average value of the exponent was found to be, *n*_{2 }*=*_{ }-.0575. So the Correlation became:

(10) |

In the heat transfer performance of looped parallel micro heat pipe, the condenser inclination has an effect. For considering the effect of condenser inclination at the transport section, a dimensionless group can be assumed as

Here, θ is the various condenser inclination and θ_{max} is the maximum inclination angle which is at the horizontal position (θ_{max} = 90^{o}). By including this dimensionless group the correlating equation became,

(11) |

In this case also using conventional techniques, exponent n_{3} of the group (1+θ/θ_{max}) was determined by plotting [U/U_{max}] / (Q/Q_{max})^{07545} X [(T_{eva}-T_{a})/T_{eva})]^{.0575}vs. (1+θ/θ_{max}) of the present experimental data as shown in Figure 18. The average value of the exponent was found to be, n_{3 }= 0.0777. So the final Correlation became:

(12) |

The fill ratio V/V_{max} has an effect on the heat transfer performance CLPHP. For considering the effect of fill ratio at the transport section, a dimensionless group can be assumed as V/V_{max}. Here, V is the various fill ratio and V_{max} is the maximum fill ratio in the experiment. By including this dimensionless group, the correlating equation became,

(13) |

In this case also using conventional techniques, exponent n_{4} of the group (V/Vmax) was determined by plotting (U/U_{max})/[(Q/Q_{max})^{07545}X[(T_{eva}-T_{a})/T_{eva})]^{.0575 }X(1+θ/θ_{max})^{.0777} vs (V/V_{max}) of the present experimental data as shown in Figure 19 The average value of the exponent was found to be, n_{4 }= 0.2021. So the final Correlation became:

(14) |

In this final correlation (equation 14), we found that the exponent of [(T_{eva}-T_{a})/T_{eva})^(-.0575)] and (1+θ/θ_{max}) is very close to zero. In other words, it can be conclude that, the effect of [(T_{eva}-T_{a})/T_{eva})^(-.0575)] and (1+θ/θ_{max}) dimensionless groups in the heat transfer characteristic of CLPHP is very insignificant. In this case, if we neglect these two terms, thus the correlation becomes:

(15) |

Then the graphical representation of the simplified correlations for LPMHP is shown in Figure 21. The complete correlation for pulsating heat pipe is given by the following equation:

(16) |

This paper aims at building up an empirical equation considering the physical parameters that directly contributes to heat transfer through the pulsating heat pipe. Heat input, Evaporator temperature rise, inclination and fill ratios are compared to form a linear empirical relationship. From the section of 4.1 and 4.2 it is recommended that the best performance can be achieved at inclination below 90˚ and fill ratios 0.4 and 0.6. It is found that overall heat transfer coefficient varies effectively with the heat input to the PHP and the fill ratios. It is evident that the final correlation converges well for linear relationship.

This work is supported by Bangladesh Army, Department of Mechanical Engineering, Bangladesh Military Academy, Bhatiary, Chittagong, Bangladesh.

A: Cross sectional area of the wall, m^{2}

A_{ci}: Inner surface area of the condenser, m^{2}

A_{ei}: Inner surface area of the evaporator, m^{2}

d_{o}: Outer diameter of the Aluminum pipe, m

d_{i}: Inner diameter of the Aluminum pipe, m

h: Convective heat transfer coefficient ,W/m^{2 }°C

K: Thermal conductivity of metal, W/m.K

*K*_{Al}: Thermal conductivity of Aluminum, W/m.K

Length of the PHP wall between the evaporator and the condenser, m

M_{Al}: Mass of Aluminum block, Kg

n: Number of turns

Heat absorbed by the block, W

Electric heat input, W

Rate of Heat transfer by the working fluid, W

Heat conduction through the pipe wall, W

*S*_{Al}: Specific heat of Aluminum, J/Kg°C

T_{a}: Ambient temperature, °C

T_{con}: Average temperature of conduction section, °C

T_{eva}: Average temperature of evaporator section, °C

T_{Al}: Rise of temperatute of Aluminum block at the time interval Δt, °C

ΔT: Temperature difference between evaporator and condenser, °C

Δt: Time interval, sec

U: Overall heat transfer coefficient, W/m^{2}. °C

V: Fill Ratio.

[1] | N. Tesla, "Valvular Conduit". New York Patent US 1329559 A, 3 February 1920. | ||

In article | |||

[2] | R. S. Gaugler, "HEAT TRANSFER DEVICE". Dayton, Ohio Patent 2,350,348, 6 June 1944. | ||

In article | |||

[3] | H. Akachi, "STRUCTURE OF MICRO-HEAT PIPE". Sagamihara, Japan Patent 5,219,020, 15 August 1991. | ||

In article | |||

[4] | G. M. Grover, "Evaporation-Condenssation Heat Transfer Device". US Patent 3229759, 2 December 1963. | ||

In article | |||

[5] | P. a. R. Dunn, Heat Pipes, Pergamon Press, 1982. | ||

In article | |||

[6] | A. Faghri, "Review and Advances in Heat Pipe Science and Technology," Journal of heat Transfer, vol. 134, no. 12, p. 18 pages, 2012. | ||

In article | View Article | ||

[7] | X. J. L. a. Z. Z. Q. Zhang X. M., "Experimental Study of a Pulsating Heat pipe using FC-72, Ethanol and Water as working fluid," vol. 12, pp. 445-472, 1999. | ||

In article | |||

[8] | B. a. F. Holley, "Analysis of Pulsating Heat Pipe with Capillary Wick and Varying Channel Diameter," International Journal of Heat and Mass Transfer , vol. 48, pp. 2635-2651, 2005. | ||

In article | View Article | ||

[9] | T. Dobson R.T. and Harms, "Lumped Parameter Analysis of Closed and Open Oscillatory Heat Pipes," in 11th International Heat Pipe Conference , Tokyo, Japan , 1999. | ||

In article | |||

[10] | S. K. M. G. Honghai Yang, "Performance characteristics of pulsating heat pipes as integral thermal spreaders," International Journal of Thermal Sciences, vol. 48, pp. 815-824, 2009. | ||

In article | View Article | ||

[11] | Khandekar S. and Groll M., "An Insight into Thermo-Hydraulic Coupling in Pulsating Heat Pipes," Int. J. of Thermal Sciences, vol. 43, no. 1, pp. 13-20, 2004. | ||

In article | View Article | ||

[12] | B. Y. W. T. N. a. O. K. T. Tong, " Closed-loop pulsating heat pipe," Applied Thermal Engineering, vol. 21, p. 1845-1862, 2001. | ||

In article | View Article | ||

[13] | A. F. Y. Z. M.B.Shafi, " Thermal modeling of unlopped and looped pulsating heat pipes," ASME J. Heat Transfer , vol. 123, pp. 1159-1172, 2001. | ||

In article | View Article | ||

[14] | P. T. T. W. P. K. M. M. P Sakulchangsatjatai, "Operation modeling of closed end and closed loop oscillating heat pipe at normal operat-ing condition," Appl. Therm. Eng., vol. 24, pp. 995-1008, 2004. | ||

In article | View Article | ||

[15] | A. M. M. F. K. a. M. M. R. Md Shahidul Haque, "Thermal characteristics of an ammonia-charged closed-loop pulsating heat pipe," Journal of Mechanical Science and Technology, vol. 33, no. 4, pp. 1907-1914, 2019. | ||

In article | View Article | ||

Published with license by Science and Education Publishing, Copyright © 2020 Md Shahidul Haque, Antara Majumdar, Md. Faisal Kader and Md. Mahbubur Razzaq

This work is licensed under a Creative Commons Attribution 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Md Shahidul Haque, Antara Majumdar, Md. Faisal Kader, Md. Mahbubur Razzaq. Empirical Correlation for an Ammonia Charged Aluminum Closed Loop Pulsating Heat Pipe. *American Journal of Mechanical Engineering*. Vol. 8, No. 1, 2020, pp 9-16. http://pubs.sciepub.com/ajme/8/1/2

Haque, Md Shahidul, et al. "Empirical Correlation for an Ammonia Charged Aluminum Closed Loop Pulsating Heat Pipe." *American Journal of Mechanical Engineering* 8.1 (2020): 9-16.

Haque, M. S. , Majumdar, A. , Kader, M. F. , & Razzaq, M. M. (2020). Empirical Correlation for an Ammonia Charged Aluminum Closed Loop Pulsating Heat Pipe. *American Journal of Mechanical Engineering*, *8*(1), 9-16.

Haque, Md Shahidul, Antara Majumdar, Md. Faisal Kader, and Md. Mahbubur Razzaq. "Empirical Correlation for an Ammonia Charged Aluminum Closed Loop Pulsating Heat Pipe." *American Journal of Mechanical Engineering* 8, no. 1 (2020): 9-16.

Share

[1] | N. Tesla, "Valvular Conduit". New York Patent US 1329559 A, 3 February 1920. | ||

In article | |||

[2] | R. S. Gaugler, "HEAT TRANSFER DEVICE". Dayton, Ohio Patent 2,350,348, 6 June 1944. | ||

In article | |||

[3] | H. Akachi, "STRUCTURE OF MICRO-HEAT PIPE". Sagamihara, Japan Patent 5,219,020, 15 August 1991. | ||

In article | |||

[4] | G. M. Grover, "Evaporation-Condenssation Heat Transfer Device". US Patent 3229759, 2 December 1963. | ||

In article | |||

[5] | P. a. R. Dunn, Heat Pipes, Pergamon Press, 1982. | ||

In article | |||

[6] | A. Faghri, "Review and Advances in Heat Pipe Science and Technology," Journal of heat Transfer, vol. 134, no. 12, p. 18 pages, 2012. | ||

In article | View Article | ||

[7] | X. J. L. a. Z. Z. Q. Zhang X. M., "Experimental Study of a Pulsating Heat pipe using FC-72, Ethanol and Water as working fluid," vol. 12, pp. 445-472, 1999. | ||

In article | |||

[8] | B. a. F. Holley, "Analysis of Pulsating Heat Pipe with Capillary Wick and Varying Channel Diameter," International Journal of Heat and Mass Transfer , vol. 48, pp. 2635-2651, 2005. | ||

In article | View Article | ||

[9] | T. Dobson R.T. and Harms, "Lumped Parameter Analysis of Closed and Open Oscillatory Heat Pipes," in 11th International Heat Pipe Conference , Tokyo, Japan , 1999. | ||

In article | |||

[10] | S. K. M. G. Honghai Yang, "Performance characteristics of pulsating heat pipes as integral thermal spreaders," International Journal of Thermal Sciences, vol. 48, pp. 815-824, 2009. | ||

In article | View Article | ||

[11] | Khandekar S. and Groll M., "An Insight into Thermo-Hydraulic Coupling in Pulsating Heat Pipes," Int. J. of Thermal Sciences, vol. 43, no. 1, pp. 13-20, 2004. | ||

In article | View Article | ||

[12] | B. Y. W. T. N. a. O. K. T. Tong, " Closed-loop pulsating heat pipe," Applied Thermal Engineering, vol. 21, p. 1845-1862, 2001. | ||

In article | View Article | ||

[13] | A. F. Y. Z. M.B.Shafi, " Thermal modeling of unlopped and looped pulsating heat pipes," ASME J. Heat Transfer , vol. 123, pp. 1159-1172, 2001. | ||

In article | View Article | ||

[14] | P. T. T. W. P. K. M. M. P Sakulchangsatjatai, "Operation modeling of closed end and closed loop oscillating heat pipe at normal operat-ing condition," Appl. Therm. Eng., vol. 24, pp. 995-1008, 2004. | ||

In article | View Article | ||

[15] | A. M. M. F. K. a. M. M. R. Md Shahidul Haque, "Thermal characteristics of an ammonia-charged closed-loop pulsating heat pipe," Journal of Mechanical Science and Technology, vol. 33, no. 4, pp. 1907-1914, 2019. | ||

In article | View Article | ||