© 2021 IIETA. This article is published by IIETA and is licensed under the CC BY 4.0 license (http://creativecommons.org/licenses/by/4.0/).
OPEN ACCESS
Pulsating heat pipe is one of the prominent technology for thermal management of electronic devices. It consists of three sections namely evaporator, adiabatic and condenser section. PHP is a two phase passive device having efficient and quick ability of transferring heat from evaporator section to condenser section. At first an 8 turn pulsating heat pipe of closed loop ends (CLPHP) with copper tube capillary dimensions is investigated experimentally for different fill ratios and for different inclinations by varying range of heat inputs. Different working fluids viz Water, Acetone, Ethanol and Methanol are considered for the experimentation. One of the recent analytical technology for modelling of CLPHPs is Artificial Neural Network (ANN) approach. The analytical models are having limited scope of applicability and they are simple in nature. The present paper describes Validation of experimental data by training prediction model ANN with available data. Three input nodes such as input heat, fill ratio and angle of inclination and one output node corresponding to PHP that is thermal resistance are considered. The feed forward neural network (FFNN) architecture is adopted for predictions. By using the physical phenomena of the system modelling are clearly known for obtaining feasible results which is main function of ANN. The predicted data validates experimental data in a satisfactory range and the results are found to be in good agreement with in the range of ± 10 percent.
angle of inclination, ANN model, CLPHP, fill ratio (FR), heat input, thermal resistance (Rth), working fluid
Pulsating heat pipes are mainly developed for control of electronic devices heat dissipation rates. It’s a novel wickless heat pipe with phase change phenomena and invented by Akachi and Polášek [1].^{. }PHP consists of Evaporator where heat input is supplied to the device, adiabatic section and condenser sections at which the vapour condenses and becomes liquid i.e. rejection of heat take place. By pulsating motion heat transfer takes place from one end to other end. Surface tension and gravitational forces plays an important role in the working mechanism of heat transfer. In order to reduce the force due to surface tension over the force due to gravitation, the diameter of the tube is selected as capillary dimensions. Due to capillary dimensions the fluid distributes itself in the tube. This causes liquid slugs and vapour bubbles in the tube. As the physical phenomena occurs in PHP are complex, the application and scope of analytical models are limited in nature. The operating mechanism of PHP is mainly depending on physical, geometrical and operational parameters which influences its performance. Simulation studies on pulsating heat pipes were conducted by Jokar et al. [2]. And natural algorithms for simulation was applied. This simulator clearly explains complex behaviour of the PHP by using nonlinear structure. Numerical model to investigate the pressure losses in the tube with respect to number of turns was developed by Sakulchangsatjatai et al. [3], and Mameli et al. [4]. In this, the correlations of heat transfer are developed for predicting the PHP performance. Using ANN and regression analysis, Prediction models for pulsating heat pipe were conducted by Patel and Mehta [5]. All the literature data related to experiments and 9 variables are considered as input to develop prediction models are collected and developed Feedforward backpropagation model to verify the data. On chaotic flow numerical investigations are carried out by Pouryoussefi and Zhang [6]. For analysis evaporative heating and condenser cooling constant wall temperatures are applied. Fully connected feed forward multilayer ANN is used to train with test data and then experimental data is validated by Khandekar et al. [7]. ANN consists of two input nodes and single output node. The results concluded that if sufficient data is available ANN approach seems to be a promising approach. Various nonlinear models based on fluctuations in the wall temperatures to understand the complex flow distribution in the PHP are implemented by Qu et al. [8], and Song and Xu [9]. A neural network of two layers to predict the PHP behaviour is studied by Shokouhmand et al. [10]. The results obtained are appropriate for predicting effective parameters trend on performance of PHP and they are in concurrent with available data. To predict the thermal performance of a PHP, Han et al. [11] carried out experiments. The thermal performance of PHP is influenced by operational parameters (heat flux applied, PHP orientation and fill ratio), inner diameter, number of turns as geometrical parameters and pure fluids, binary mixtures and nano fluids are as working fluids. The correlation for startup heat flux and startup mechanism of PHP is theoretically described by Qu and Ma [12]. The results explains that the influence of working fluids makes the PHP’s internal heat and mass transfer phenomena more complex.
From the literature it is observed that closed loop pulsating heat pipe prediction models are developed with limited fluids at a specified fill ratio and orientation only. The main objective of the present work is to carry out a comprehensive work on water, ethanol, methanol and acetone fluids with 50%, 60% and 75% fill ratios at 0degree, 45degree and 90degree positions. The experimental values are compared with prediction values. The feed forward neural network is used over other methods because in this method only in one direction that is in the forward direction the information always moves. The connections between the nodes do not form a cycle that is it does not go backwards. The transfer of nodes is from input to output nodes through the hidden nodes.
To understand the pulsating motion of liquid slugs and vapour bubbles of fluid in tube and also effect of parameters influence on CLPHP many prediction models are developed. But the number of parameters prediction is very much limited. To predict the thermal performance of PHP the correlations of heat transfer mechanism including all the parameters are so far not available. Keeping this in view to develop prediction model for thermal performance of PHP the present research work is carried out by using ANN approach. ANN approach is one of the modelling techniques used to predict the data where the understanding of mechanism in a device is difficult.
The design and functioning of human brain and nervous system is the main motivation for the design of ANN device. It’s a datadriven machine learning tool. It’s neither an algorithm nor actual hardware. All the processing elements are interconnected and makes different layers. These layers are useful to obtain a relationship between inputs and outputs, learning and memorizing the data. Usually, the structure of ANN consists of three layers shown in Figure 1. For receiving data input layer, for processing hidden layer and to send computed information an output layer. The artificial neurons are represented by nodes of each layer. Through weights the interconnection between the nodes and its preceding layer is possible. To match with the required output through the learning process by adjusting of associated weights the desired output is obtained. To model complex functions all, the individual neurons are connected with each other. The output vector is obtained when input vector is accepted by network.
Figure 1. Architecture of artificial neural network [7]
By ANN always it is not possible to get prediction of new data. This is due to overtraining the data. To avoid this data is set into different parts. For training process almost 70% of data is assigned. In the remaining 30%, 15% of data us used for testing and the remaining 15% is for validation. For developing accurate model more data is required. Therefore, training data is used for this purpose. To train the network the learning algorithm is used. The algorithm measures actual and desired output vector difference. The error resulted uses an iterative method for calculating weights and to improve network such a way to reduce error. Several times this process is repeated till the error is in the required level so that it becomes valid model for prediction. To reduce the error to check generalization and accuracy the validation data is used. The final behaviour of network is examined by validating data.
The experimental work is carried out on an 8 turn closed loop pulsating heat pipe made up of copper. Copper tube with 2 mm and 3.1 mm capillary dimensions and of length 264 mm was bent into 8 Uturns into major section of CLPHP such as evaporator, adiabatic and condenser section. The evaporator box of size 330X90X55 mm are placed at the bottom and the condenser of size 320X95X55 mm are placed at the top to insert the parallel tubes. The CLPHP performance is mainly depending upon the parameters like working fluid, volume fill ratio, inclination angle, number of turns, inner diameter of the tube, heat input supplied etc. Fill ratio is defined as the ratio of liquid volume to the total volume of PHP. It varies from 0% to 100% i.e. CLPHP tubes without any fluid (bare tubes) or empty device means no fluid in the tube is referred as 0% filling. In this case the heat transfer in CLPHP is only by pure conduction mode. A fully filled tube refers to 100% filling ratio. In this case there is no pulsating effect in the tube. Due to buoyancy the liquid circulates in the tube and buoyancy helps for heat transfer. Therefore, between 10% and 90% fill ratios only CLPHP can function in a pulsating mode. The exact value of fill ratio is mainly depending upon working fluids, construction of PHP and other parameters that effect PHP. Different mass flow rates are maintained for different heat inputs. The mass flow rate depends upon the properties such as its specific heat, density etc. of working fluid considered.
Heat is supplied at evaporator section of CLPHP. The liquid converts into vapour at evaporator and gets condensed in the condenser section. Amount of heat supplied is mainly depending upon the type of fluid used and its boiling point. The orientation of CLPHP also another important parameter which effects the performance. Different inclinations varying from 0 degree to 180 degrees is possible. In all these cases the physics that is underlined for understanding of heat transfer mechanism is not so clear and is differs with respect to its operation.
The schematic diagram of an 8 turn CLPHP experimental setup is as shown in Figure 2. The experiment is conducted for working fluids namely Water, Ethanol, Methanol and Acetone. By incorporating thermocouples at evaporator and condenser sections, the temperatures were recorded and heat input is supplied at bottom end of evaporator section by using variac. The ratio of difference in temperature with its heat input gives the performance parameter thermal resistance, the required output.
Figure 2. Schematic of experimental setup
The experimental results are noted for all the working fluids considered at different fill ratios (50%, 60% and 75%) and also at different orientations (0degree, 45degrees and 90degree). These results are validated by using ANN.
To develop the Artificial Neural Network (ANN) model a neural network toolbox in MATLAB was used This model was developed by considering Heat input (20100 W) depending upon the working fluid, Fill ratio (50%, 60% and 75%) and inclination angle (0degree, 45degrees and 90degree) as input parameters and thermal resistance as output parameter. The following Tables 14 gives different process parameters considered for different working fluids namely Water, Acetone, Ethanol and Methanol (Considered for the present study).
Table 1. Process parameters and their levels for water CLPHP
Process parameters 
Levels 

1 
2 
3 

Heat Input 
20 
60 
100 
Fill Ratio (%) 
50 
60 
75 
Inclination angle 
0 
45 
90 
Table 2. Process parameters and their levels for acetone CLPHP
Process parameters 
Levels 

1 
2 
3 

Heat Input 
20 
40 
60 
Fill Ratio (%) 
50 
60 
75 
Inclination angle 
0 
45 
90 
Table 3. Process parameters and their levels for ethanol CLPHP
Process parameters 
Levels 

1 
2 
3 

Heat Input 
20 
40 
60 
Fill Ratio (%) 
50 
60 
75 
Inclination angle 
0 
45 
90 
Table 4. Process parameters and their levels for methanol CLPHP
Process parameters 
Levels 

1 
2 
3 

Heat Input 
20 
40 
60 
Fill Ratio (%) 
50 
60 
75 
Inclination angle 
0 
45 
90 
Initially, into the neural network the input data were fed and then to achieve the output it was preceded for simulation in Figure 3. By using 3^{3} factorial designs simulations are carried out taking into account of maximum and minimum limits of the fluid considered. A total of 27 experiments investigations were considered to predict the data using ANN architecture namely Feed Forward Neural Network (FFNN).
Figure 3. Feed forward back propagation ANN
The Tables 5a and 5b describes the experimental plan for CLPHP modeling. Different operating parameters that affect the CLPHP performance are considered in different combinations. The combination includes heat input, fill ratio and angle of inclination for different fluids Water, Acetone, Ethanol and Methanol.
Table 5a. Experimental plan for modeling water and acetone CLPHP
Water 
Acetone 

Q 
FR 
Angle 
Q 
FR 
Angle 
20 
50 
0 
60 
50 
0 
60 
50 
0 
60 
50 
45 
100 
50 
0 
60 
50 
90 
20 
50 
45 
60 
60 
0 
60 
50 
45 
60 
75 
0 
100 
50 
45 
60 
60 
45 
20 
50 
90 
60 
60 
90 
60 
50 
90 
60 
75 
45 
100 
50 
90 
60 
75 
90 
20 
60 
0 
40 
50 
0 
60 
60 
0 
40 
50 
45 
100 
60 
0 
40 
60 
0 
20 
60 
45 
40 
50 
90 
60 
60 
45 
40 
60 
45 
100 
60 
45 
40 
75 
0 
20 
60 
90 
40 
60 
90 
60 
60 
90 
40 
75 
45 
100 
60 
90 
40 
75 
90 
20 
75 
0 
20 
50 
0 
60 
75 
0 
20 
60 
0 
100 
75 
0 
20 
50 
45 
20 
75 
45 
20 
75 
0 
60 
75 
45 
20 
60 
45 
100 
75 
45 
20 
50 
90 
20 
75 
90 
20 
60 
90 
60 
75 
90 
20 
75 
45 
100 
75 
90 
20 
75 
90 
Table 5b. Experimental plan for modeling ethanol and methanol CLPHP
Ethanol 
Methanol 

Q 
FR 
Angle 
Q 
FR 
Angle 
60 
50 
0 
60 
75 
45 
60 
50 
45 
60 
50 
0 
60 
50 
90 
60 
50 
45 
60 
60 
0 
60 
50 
90 
60 
60 
45 
60 
60 
0 
60 
75 
0 
60 
75 
0 
60 
75 
45 
60 
60 
45 
60 
75 
90 
60 
60 
90 
60 
60 
90 
60 
75 
90 
40 
50 
0 
40 
50 
0 
40 
50 
45 
40 
50 
45 
40 
50 
90 
40 
60 
0 
40 
75 
0 
40 
50 
90 
40 
60 
0 
40 
60 
45 
40 
75 
90 
40 
60 
90 
40 
75 
45 
40 
75 
0 
40 
60 
45 
40 
75 
45 
40 
60 
90 
40 
75 
90 
20 
50 
0 
20 
50 
0 
20 
50 
45 
20 
60 
0 
20 
60 
0 
20 
75 
0 
20 
60 
45 
20 
50 
45 
20 
75 
0 
20 
60 
45 
20 
60 
90 
20 
50 
90 
20 
50 
90 
20 
60 
90 
20 
75 
45 
20 
75 
45 
20 
75 
90 
20 
75 
90 
By using TRAINLM and LEARNGDM the training and Adaption learning functions are considered as shown in Figure 4.
Figure 4. Creating network using ANN
Figure 5. Neural network training
For every neural network model the following model is adopted.
The output thermal resistance is determined by using experimental values for all considerations adopted for experimentation. These resistance values are validated by using ANN. The error can be calculated as percentage error by using the following formula.
Error (%) =$\frac{(\text { Rth }(\text { prediction })\text { Rth }(\text { Experimental }))}{\text { Rth(Experimental) }}$
5.1 For working fluid water
The following Figures 68 show the comparison of water at different fill ratio and at different inclinations by considering both ANN and Experimental result.
5.1.1 At 50% fill ratio
The percentage error observed for all the readings regarding CLPHP operating with working fluid water when operating at 50%FR at different inclination angles was below ±5% range.
Figure 6. Comparison of experimental Rth with predicted ANN at 50% FR and 0,45,90degree inclinations water CLPHP
5.1.2 At 60% fill ratio
The percentage error observed for all the readings regarding PHP operating with working fluid water when operating at 60%FR at different inclination angles was below ±10% range.
Figure 7. Comparison of experimental Rth with predicted ANN at 60% FR and 0,45,90degree inclinations for water CLPHP
5.1.3 At 75% fill ratio
The percentage error observed for all the readings regarding PHP operating with working fluid water when operating at 75%FR at different inclination angles was below ±13% range. The following graph Figure 8 shows these variation.
Among all the working conditions of CLPHP with water, when it is operating at 50% fill ratio is having an error of less value i.e ±5% compared with 60% and 75% fill ratios.Thermal resistance value are high for higher filling ratio compared to lower filling ratio. At low heat transfer limit, higher thermal resistance is observed and performance of CLPHP is very poor. For high filling ratios, sufficient bubbles are not developed to provide pumping action in the fluid. So at this FR drastic deterioration of performance occurs.The formation of bubbles will results in better heat transfer rates and hence performnce improvement.
Figure 8. Comparison of experimental Rth with predicted ANN at 75% FR and 0,45,90degree inclinations for water CLPHP
5.2 For working fluid acetone
Acetone has low value of latent heat.At a given temperature and presure quck evaporation of liquid takes place has it has low latent heat value i.e 520.56 kJ/kg.The velocities of liquid slugs wil be increased and results in improvement of heat ransfer.The physical and chemical properties of fluid are: boilin point is 56.2℃, density 748.5 kg/m^{3} in liqud form, specific heat 2.28 kJ/kgK in liquid form, 1.385 kJ/kgK for vapor acetone and its density is 19.09X10^{3} N/m. The following Figure 9 shows comparison of Acetone fluid at different fill ratios and at different inclination angles by considering experimental results and ANN predictions.
When the CLPHP with acetone as working fluid was operating at 50% fill ratio and at 0degree inclination the error obtained was ±6% and for 45degree and 90degree the percentage of errors are ±1% and ±1% respectively.
When the CLPHP with working fluid acetone was operating at 60% fill ratio and at 0degree inclination the error obtained was ±2% and for 45degree and 90 degree the percentage of errors are ±3% and ±5% respectively as shown in Figure 10.
Figure 9. Comparison of experimental Rth with predicted ANN at 50% FR and 0,45,90degree inclinations for acetone CLPHP
Figure 10. Comparison of experimental Rth with predicted ANN at 60% FR and 0,45,90degree inclinations for acetone CLPHP
Figure 11. Comparison of experimental Rth with predicted ANN value at 75% FR and 0,45,90degree inclinations for acetone CLPHP
When CLPHP with working fluid acetone was operating at 75% fill ratio and at 0degree inclination the error obtained was ±9%. where as for 45degree and 90 degree the percentage of errors are ±1% and ±6% respectively shown in Figure 11. And observed that 5% error is observed for acetone flud at 60% fill ratio and 90degree inclination position.
5.3 For working fluid ethanol
The following set of graphs (Figures 12, 13 and 14) shows the comparative analysis of experimentally obtained thermal resistance values with ANN predictions values obtained by training in nntool for PHP working with ethanol as working fluid. The heat input is considered in the range 2080 watts and fill ratios considered are 50%, 60% and 75% at inclinations of 0degree, 45degree and 90degree positions.
Figure 12. Comparison of experimental Rth with predicted ANN at 50% FR and 0,45,90degree inclinations for ethanol CLPHP
When CLPHP is operating with Ethanol working fluid at 50% filling ratio, (i) at 0degree position the error in thermal resistance obtained is ±10% (ii) at 45degree position inclination the error is ±11% (iii) and at an inclination angle 90degree it was observed as ±3%.
Figure 13. Comparison of experimental Rth with predicted ANN at 60% FR and 0,45,90degree inclinations for ethanol CLPHP
When CLPHP is operating with Ethanol working fluid at 60% filling ratio, (i) at 0degree position the error in thermal resistance obtained is ±3% (ii) at 45degree position inclination the error is ±8% (iii) and at an inclination angle 90degree it was observed as ±9%.
Figure 14. Comparison of experimental Rth with predicted ANN at 75% FR and 0,45,90degree inclinations for ethanol CLPHP
When CLPHP is operating with Ethanol working fluid at 75% filling ratio, (i) at 0degree position the error in thermal resistance obtained is ±2% (ii) at 45degree position inclination the error is ±1% (iii) and at an inclination angle 90degree it was observed as ±3%.
5.4 For working fluid methanol
The experiments are carried out for CLPHP by considering methanol as working fluid. The temperatures of evaporator and condenser at different locations are recorded by using data logger and hence thermal resistance was found by varying heat input. Fill ratios of 50%, 60% and 75% are considered. All the experiments are carried at 0 –degree, 45degree and 90degree inclinations of PHP. The results obtained are validated with ANN tool. The following group of curves (Figures 15, 16 and 17) represents the variations of Rth in both experimental values and also predicted ANN values.
Figure 15. Comparison of experimental Rth with predicted ANN at 50% FR and 0,45,90degree inclinations for methanol CLPHP
When methanol CLPHP is considered at 50% fill ratio the error percentage obtained experimentally in comparison with predicted ANN thermal resistance was different from one position to other position. At 0degree inclination position the error is in the range of ±8% whereas for 45 and 90degree positions it is in the range of ±4% only. This is due to the gravity that acts in 0degree position.
Figure 16. Comparison of experimental Rth with predicted ANN Rth at 60% FR and 0,45,90degree inclinations for methanol CLPHP
For CLPHP operating at 60% fill ratio the error percentage of thermal resistance in comparison with experimental and predicted ANN values is at 0degree and 90–degree inclination positions ±3% and at 45degree position observed as ±9%.
Figure 17. Comparison of experimental Rth with predicted ANN at 75% FR and 0,45,90degree inclinations for methanol CLPHP
For CLPHP operating at 75% fill ratio the error percentage of thermal resistance obtained in comparison with experimental and predicted ANN values is at 0degree±2%, at 45degree ±9%. And at 90 degrees it is 0.2 percent.
[1] Akachi, H., Polášek, F. (1996). Pulsating Heat Pipes. In: Heat Pipe Symposium, Melbourne, Australia, 208217.
[2] Jokar, A., Godarzi, A.A., Saber, M., Shafii, M.B. (2016). Simulation and optimization of a pulsating heat pipe using artificial neural network and genetic algorithm. Heat and Mass Transfer, 52(11): 24372445. https://doi.org/10.1007/s0023101617598
[3] Sakulchangsatjatai, P., Terdtoon, P., Wongratanaphisan, T., Kamonpet, P., Murakami, M. (2004). Operation modeling of closedend and closedloop oscillating heat pipes at normal operating condition. Applied Thermal Engineering, 24(7): 9951008. https://doi.org/10.1016/j.applthermaleng.2003.11.006
[4] Mameli, M., Marengo, M., Zinna, S. (2012). Numerical model of a multiturn closed loop pulsating heat pipe: Effects of the local pressure losses due to meanderings. International Journal of Heat and Mass Transfer, 55(4): 10361047. https://doi.org/10.1016/j.ijheatmasstransfer.2011.10.006
[5] Patel, V.M., Mehta, H.B. (2016). Artificial neural network modeling of a closed loop pulsating heat pipe. World Academy of Science, Engineering, and TechnologyInternational Journal of Mechanical, Aerospace, Industrial, Mechatronics, and Manufacturing Engineering, 10: 17541757.
[6] Pouryoussefi, S.M., Zhang, Y. (2016). Numerical investigation of chaotic flow in a 2D closedloop pulsating heat pipe. Applied Thermal Engineering, 98: 617627. https://doi.org/10.1016/j.applthermaleng.2015.12.097
[7] Khandekar, S., Cui, X., Groll, M. (2002). Thermal performance modeling of pulsating heat pipes by artificial neural network. In Proceedings of 12th International Heat Pipe Conference, pp. 215219.
[8] Qu, J., Wu, H., Cheng, P., Wang, X. (2009). Nonlinear analyses of temperature oscillations in a closedloop pulsating heat pipe. International Journal of Heat and Mass Transfer, 52(1516): 34813489. https://doi.org/10.1016/j.ijheatmasstransfer.2009.03.012
[9] Song, Y.X., Xu, J.L. (2009). Chaotic behaviour of pulsating heat pipes. International Journal of Heat and Mass Transfer, 52(1314): 29322941. https://doi.org/10.1016/j.ijheatmasstransfer.2009.02.030
[10] Shokouhmand, H., Gharib, N., Bahrami, H. (2006). Modeling of closed loop pulsating heat pipes by neural networks. Engineering Systems Design and Analysis, 42487: 131136. https://doi.org/10.1115/ESDA200695417
[11] Han, X., Wang, X., Zheng, H., Xu, X., Chen, G. (2016). Review of the development of pulsating heat pipe for heat dissipation. Renewable and Sustainable Energy Reviews, 59: 692709. https://doi.org/10.1016/j.rser.2015.12.350
[12] Qu, W., Ma, H.B. (2007). Theoretical analysis of startup of a pulsating heat pipe. International Journal of Heat and Mass Transfer, 50(1112): 23092316. https://doi.org/10.1016/j.ijheatmasstransfer.2006.10.043