Issue 
Mechanics & Industry
Volume 21, Number 3, 2020



Article Number  309  
Number of page(s)  10  
DOI  https://doi.org/10.1051/meca/2020021  
Published online  27 April 2020 
Regular Article
An investigation on effect of EDL on heat transfer of micro heat pipe with square and triangular cross section
^{1}
Department of Mechanical Engineering, Science and Research Branch, Islamic Azad University, Tehran, Iran
^{2}
School of Mechanical Engineering, College of Eng. University of Tehran, Tehran, Iran
^{*} email: ho.shokouhmand@gmail.com
Received:
8
December
2017
Accepted:
9
February
2020
Electronic industries have always been trying to improve the efficiency of electronic devices with small dimensions through thermal management of this equipment, thus increasing the use of small thermal sinks. In this study micro heat pipes with triangular and square cross sections have been manufactured and tested. One of the main objectives is to obtain an understanding of micro heat pipes and their role in energy transmission with electrical double layer (EDL). Micro heat pipes are highly efficient heat transfer devices, which use the continuous evaporation/condensation of a suitable working fluid for twophase heat transport in a closed system. Since the latent heat of vaporization is very large, heat pipes transport heat at small temperature difference, with high rates. Because of variety of advantage features these devices have found a number of applications both in space and terrestrial technologies. The theory of operation micro heat pipes with EDL is described and the micro heat pipe has been studied. The temperature distribution have achieved through five thermocouples installed on the body. Water and different solution mixture of water and ethanol have used to investigate effect of the electric double layer heat transfer. It was noticed that the electric double layer of ionized fluid has caused reduction of heat transfer.
Key words: heat transfer / micro heat pipe / electrical double layer / condensation / evaporation
© AFM, EDP Sciences 2020
1 Introduction
Components, mechanical and electrical devices with micronsize are becoming more common in commercial applications and scientific research. Rapid progress in microelectric has created a basic revolution in computers, electronics devices, communications systems and many other tools. Lightweight, compact and highperformance of microsystems will have important applications in transportation, the environment and space research, biochemical and other industries of chemical process. The main advantages of these systems and microfluids are including high rates of heat transfer, mass transfer and chemical reactions and reduced size of the samples [1].
With increasing heat dissipation of microelectronic devices and reduced size of them, thermal management becomes a more important element than the electronic product design. So, long life and reliable performance of parts may be achieved by effective control of device operating temperature with ranges adjustment by engineering design. Heat sinks are devices that facilitate the release of heat from hot surfaces that usually, transfer the heat generated by the components to the cooler environment. The primary objective of the heat sink is keeping the device temperature in the lower the maximum temperature allowed by the manufacturer.
Mainly, micro devices have a different behavior than macro devices. That it can be attributed to several factors, among this is electrokinetic effects. Electrokinetic transport applications in micro and nanodimension are important. Electric double layer (EDL) in the containing ion flows has significant effects on Heat transfer and fluid flow inside the micro heat pipes. So that the usual equations of heat transfer and flow in the channels about micro heat pipe loses its credibility [2].
Induced surface charge at the interface between solid and liquid phase causes the formation of an electrical double layer and the nonuniform distribution of fluid ions on the surface. In transfer of an electrolyte by pressuredriving way inside a micropipe, electrical double layer formation and movement of loads in it along with applied fluid creates electric current which is called streaming.
In the past decades, a number of cooling system designs including simple flat thermal loser and more complex, such as twophase heat sink cooling microjet have been proposed and widely have been studied experimentally and theoretically. Knight et al. [3] worked in the field of allowed optimization microchannel heat sink by reducing the aspect ratio, the ratio of the width of the blade channel, Nusselt number and type of flow and volumetric rate of flow. Groll et al. [4] have investigated number of articles about the micro heat pipes for thermal control of electronic equipment. They have compared and made microclassification more than 100 articles about sinks heat pipe based on their performance, their design, and their thermal physical parameters. Applied fluids and wall material is discussed and combination possibilities are offered. Le Berre et al. [5] accomplished experimental study on the performance of an array of micro heat pipe for different charge of liquid in different experimental conditions. Micro heat pipe dimensions are 20 mm × 20 mm and contain 27 parallel channels with a triangular cross section channel, width and depth of 350–500, respectively.
2 Material and method
2.1 Governing equation
Micro heat pipes function is much like conventional heat pipes. Vaporizing of working fluid in the evaporator result in heat absorption from the heat source and Steam move toward condenser and in there condense and release its latent heat. Micro heat pipes, contrary to conventional heat pipes, use grooves or sharp edges to return the fluid from the evaporator to the condenser, Instead of using porous media (Fig. 1).
Solution area has two phases of Vapor and liquid. In addition, electrical area which is existed has influence on solution area. To track the interface between the liquid and vapor phase Lagrangian and Eulerian methods is available. Lagrangian method scope is limited to flows in which changing in surface form are not extreme. Lagrangian method is less accurate than Eulerian methods. But it can use to simulate flows with complex interface (here vapour–liquid). One of the most famous Eulerian method is VOF or volume of fluid method that defined by volume fraction of fluid in the cell network computing solutions and it's expressed by α. If the cell not containing fluid, the amount of α is zero. And if it is completely fluid filled, the amount of α is one. And if part of it is filled, α is between zero and one. In this way, the contact area between the liquid and vapor phases will be achieved as shown in Figure 2 [2,6,7].
Fig. 1
Schematic of the micro heat pipe function and utilizing groove to return the fluid from the evaporator to the condenser [16]. 
2.2 NavierStokes equations
The conservation equations in the singlefield representation must account for differences in material properties of individual fluids as well as surface tension acting at the interface. For the entire flow field, the continuity and momentum equations are expressed as:(1) (2) where v is velocity, τ is the rate of strain tensor, p is pressure. The last term in equation (2) represents the source of momentum due to surface tension effects.
The stress tensor τ is expressed as:(3)where μ is the dynamic viscosity of the “mixture” and ∇v ^{ T } is the transpose of the(4) (5)where T is vector transpose of matrix and µ is dynamic viscosity of the mixture and at the end the momentum equation which is a suitable equation for micro pipes would be [8].
2.3 VOF equation
VOF equation based on transmission of the volume fraction α is written as follows:(6)and physical properties of the mixture would be as follows:(7) (8)
However, according to both liquid and vapor phase, Eulerian equation VOF method is based on the volume fraction of α and the velocity of liquid and vapor phases are as follows [6,9,10]:(9) (10)
The indexes l and g respectively represent the phase velocity of the liquid and vapor areas. According to the Expression, a weight expression for total velocity based on velocity of the liquid and vapor phases can be provided as follows:(11)
The general VOF equation for the twophase vaporliquid area is as follows:(12)
2.4 Energy equation
In the above equation, the term S_{h} is any source of heat such as radiation and k_{eff} is effective thermal conductivity which it were calculated for the cell contains two phases as equations (4) and (5), and T is the temperature and is calculated for the liquid and vapor phases using the following formula:(14)
2.5 Conductive heat transfer (wall)
The wall has conductive heat transfer and its governing equation is [6]:(15)where q″′ is the heat generated by the heat source and k_{s} is wall heat transfer coefficient (solid) and C_{p} is specific heat.
2.6 The effect of body force
If φ be the external electric potential on micro heat pipe, the body force which is caused by that is written as follows [6,11]:(16) (17) (18)
In the above equations, k is inverse of Debye length and F is Faraday constant and R is the gases constant and T is the absolute temperature and c_{i,0} the mass concentration of ion i and z_{i} is iValence of particle and ψ_{1} is potential of the wall flux. Electric potential Ψ_{2} and volumetric charge density ρ_{e} of the electric double layer, caused by liquid and solid contacts with the condition of symmetry, are expressed by the PoissonBoltzmann equation as follows [12,13].(19) (20)where ε is the dielectric constant and n_{0} is bulk ionic concentration and z is ions valence and e is protons flux and k_{b} is Boltzmann constant and T is absolute temperature. Body force, which is caused by electric double layer between the liquid and solid phases, is obtained from the following equation(21)
and E is the electroviscous force. Body force, which is caused by electric double layer between the liquid and vapor phases, is as follows [14]:(22)where U_{sum} is concerned to intermolecular potential and z_{i} is the charge number and e is initial charge and V_{e} is the electrostatic potential and value of ψ_{3} achieved by following formula [20](23)
In the above equation, µ^{S} is the chemical potential of ions at the interface and q is the value of ion charge. In the present research, apart from the gravity, there are electric forces that affect the momentum. That all these forces are:(24)
3 Setup used in research
Making of micro heat pipe was significantly complex because of its small dimension. In this research, the micro copper pipe with triangular and square cross section were used which were prepared from commercially cylindrical pipe which was available. Due to the need to determine the exact temperature in the desired areas, cross section and the pipe diameter before deformation was brought to the desired dimensions by using heat and tension (Figs. 3 and 4; Tab. 1).
After the design and manufacturing of micro pipes, quality of working fluid charge is very important. Presence of noncondensable gases leads to disorder of fluid flow inside the heat pipe and disrupt the performance of heat pipes. At worst condition completely stopped the flow of fluid and the performance of heat pipe and heat pipe acts as a conductor.
To avoid this situation, in manufacturing process of required samples, firstly the end of micro heat pipe with a length longer than samples in above table was completely blocked and suction was fulfilled from other side. Then working fluid is injected as amount as 20% of volume of required micro pipe and leakage test was conducted after severing additional section and blocking the beginning of micro pipe.
Test equipment for thermal performance of MHP consists of a vacuum chamber and a heat transfer system to cool the condenser. Evaporator consists of a heater with a 0.36 diameter wire as heat source. This heater is settled around the body of copper pipe for providing a constant thermal load and is connected and linked to the evaporator wall via an external link. The condenser is cooled by the circulation of water flow. Vacuum chamber which is made of acrylic material to minimize heat loss to the environment was used that is shown in Figure 5.
K type Thermocouples (0.08 mm) are located by soldering in two section on the evaporator wall, one point in the adiabatic section, and two points in the condenser. The locations of the thermocouples are shown in Schematic of test set up as X that these dimensions are for the MHP with a 50 mm total length (Figs. 5 and 6). An adapter has been used to set the sensors (Fig. 7).
MHP is made of oxygenfree copper and its wall thickness is between 0.27 and 0.28 mm pure water and also ordinary water and different solution mixture of water and ethanol as a working fluid is used, which has a relatively large surface tension and Specific latent heat between 30 and 160 °C operating temperature. Filling ratio of fluid is 20% of the internal volume of the MHP in which steam is not possible, may occur in the condenser. Heat transfer by changing working fluid phase could not be completed in blocked fluid region. The temperature of end of the condenser is 5–20 °C less than the area adjacent to the condenser, and heat transfer rate of MHP due to blocking of the liquid has dropped. Thus, inactive amount of working fluid due to the closure of the design phase MHP should be considered for the calculation of filling ratio of fluid. The details of blocking mechanism have been not reported in any paper and still remain an underdeveloped subject.
Fig. 3
Cross section of micro heat pipe (triangular and square curves type). 
Fig. 4
Special tool designed to form a pipe to the triangular and square crosssection. 
Characteristics of examples of micro heat pipes in experimental test.
Fig. 5
Location and installation of thermocouple. 
Fig. 7
Micro heat pipe adapter. 
4 Working fluid
In this study, water and different solution mixture of water and ethanol were used as working fluid. Considering that the dimensions, very small and the weight of the nonfluids, were considered to be extremely important to obtain a solution with the desired concentration, first a certain volume of water was injected into a laboratory balloon by a pipette. Then the balloon was set to the laboratory scale and set to zero on zero and using the pipette, ethanol was added to the desired concentration. The balloon containing water and ethanol was produced in a magnetic stirring for five minutes (Fig. 8).
Fig. 8
Steps for preparing the water solution and ethanol. 
5 Heat transfer characteristics in micro heat pipe
Heat pipes can transfer a lot of heat via a little temperature difference between the evaporator and condenser. In general, measuring the temperature difference between the evaporator and condenser is one of the test methods to check whether there is a noncompressible gas in heat pipe or a heat pipe works well after the end of the production process. In the micro heat pipe, the production process requires highprecision technologies. Due to the existence of noncondensable gases or pollutants, even trace amounts may be harmful to the performance of heat pipes.
Figure 9 shows temperature distribution in the axial length of 50 mm in micro heat pipe with triangular and square cross section. Tested MHPs have curved triangular and square cross section with a filling ratio of 20% of internal total volume. In the condenser, heat was dispersed via conduction heat transfer. Temperature after a steady state is averaged for at least one minute to reduce minor fluctuations error. As shown in Figure 9, MHP wall temperature by increasing the heat load increases. This means that thermal equilibrium, in other word the isotherm properties of MHP from the evaporator to condenser is fulfilled well. The temperature difference between the evaporator and condenser in different thermal loads is shown in the next curves. However, the temperature difference between the evaporator and condenser, when the working fluid was changed to water and also mixture of water and ethanol was more than that for pure water. This is because the electric double layer leads to disruption of the flow inside the pipe. This due to the small dimensions of the micro heat pipe and also ions of the working fluid influenced performance of micro pipe clearly and result in greater temperature difference between the two section of evaporator and the condenser.
Fig. 9
Temperature distributions in the axial length of 50 mm in micro pipe with. 
6 Software modeling
In addition to experimental tests, micro heat pipe was simulated by the software. It should be noted that the effects of electrical double layer was not considered in this modeling. In simulating the thermodynamic properties of working fluid, the initial conditions and geometry of the micro heat pipeline were considered in accordance with the specimen. The results obtained from the software and the experiments were compared (Figs. 10–12).
Fig. 10
Meshing that used for boundary conditions. 
Fig. 11
Thermal contour for input current. 
Fig. 12
3D model of MHP. 
7 Error analysis
Due to the accuracy of measurement equipment and its effect on the error rate of computational quantities, including Reynolds number, Nusselt number, and heat transfer coefficient, the analysis of these errors is necessary. In this study, Coleman and Steele method is used to calculate the uncertainty of variables. According to this method, if the quantity of R is the function of the variables x_{1}, x_{2} then the uncertainties are derived from the following equations:(25) (26)
The accuracy of the measurement devices in the experimental system is given in the following table presents the uncertainty of the computational quantities (Tabs. 2 and 3).
Accuracy of measurement devices.
Uncertainty of computational variables.
8 Discussion and conclusion
The effect of the cross section shape on the function of micro heat pipes with triangle and square crosssection with thermal loads of 1–4 watts has been studied theoretically and experimentally. Triangular and square micro heat pipes with thin structures were manufactured and tested. Water and different solution mixture of water and ethanol were used as working fluid. MHP which is charged with pure water as working fluid represented good isotherm properties over pipe length, and the temperature difference between the evaporator and condenser was below 16 ^{°}C. Effect of gradient angle on thermal performance was low and thermal characteristics from low heating mode to high heating mode were stable. In use of ordinary water, wall temperature has been higher. This reflects the effect of an electric double layer in the use of ordinary water. MHP which is made in the present study showed a high capacity for heat dissipation and also represented that in used Small scale pipes, the impact of Electrical double layer negatively affected the performance of heat pipes by 20%. Thus it will be able to be used widely in integrated electronic devices as a cooling unit. Triangular Micro heat pipes with the same geometry and working fluid were tested as a reference to investigate the setup and comparing the results with previous reference and studies [15] and the accuracy of regulation of micro heat pipes were confirmed (Fig. 13).
To test the effect of EDL on heat transfer of micro heat pipe, with caring to the influence of alcohol concentration in a solution of water and alcohol, we used various solution alcohol and water in experimental tests and numerical study.
The difference between the experimental and theoretical results in the input power of 1 watt and the solution of water and 0.018 (gr) ethanol was about 0.7 degrees in triangular MHP But this difference in square cross section was far less; of course, the highest dispute in the adiabatic zone was 4% in triangular cross section (Fig. 14).
As expected, the effect of EDL increased by adding more drops of alcohol to water. By increasing the concentration of ethanol to 3 times in the solution (0.054 gr), the difference between theoretical and experimental results was significantly different in both triangular and square MHPs, so that we observed an average of 5% of the difference (Fig. 15). At the maximum concentration of the tested solution, 10 drops of ethanol in solution (18 gr), the theoretical results close to the results of low concentrations of ethanol were obtained. While the experimental results showed a significant difference between the results at different concentrations, this was due to the EDL effects (Fig. 16).
It was observed the difference between experimental results and extracted data from the theoretical solution have been increased by increasing input power to 4 watt, so in this state the minimum of difference was 2% and maximum in the worst state was 17% (Figs. 17–20). But the difference between experimental and theoretical results in MHP with square crosssectional area due to decreasing effects of EDL because of the geometry of corners was lower than that with triangular cross section.
In both triangular and square MHP, high amount of difference in the EDL influence of 1 or 2 watt input power was not seen. While in the 4 watt input power in addition of differences between experimental and theoretical results because of the existence of EDL factor, temperature difference between evaporator and condenser was around 27%; in other words, the existence of EDL has a negative effect on efficiency of both square and triangular cross section MHP. It was noticed that these effects were observed in triangular section due to its geometry more than the square section. Corners in square samples manufactured and tested in this study were normal, but in samples with curved corners such as triangular sample, it is expected that the micro heat pipe efficiency and also the effect of electrical double layer increase. This is due to one additional corner and the increase of the capillary pressure. It is suggested that this issue will be discussed in future researches.
Fig. 13
Validation of sample thermal resistance tested in this case with previous record and researches. 
Fig. 14
Compare experimental and theoretical graph of wall temperature, water and 0.018 g ethanol as working fluid, 1 W. 
Fig. 15
Compare experimental and theoretical graph of wall temperature, water and 0.054 g ethanol as working fluid, 1 W. 
Fig. 16
Compare experimental and theoretical graph of wall temperature, water and 0.18 g ethanol as working, 1 W. 
Fig. 17
Compare experimental and theoretical graph of wall temperature, water and 0.018 g ethanol as working, 4 W. 
Fig. 18
Compare experimental and theoretical graph of wall temperature, water and 0.054 g ethanol as working, 4 W. 
Fig. 19
Compare experimental and theoretical graph of wall temperature, water and 0.18 g ethanol as working, 1 W. 
Fig. 20
Compare experimental and theoretical graph of wall temperature, water and 0.054 g ethanol as working, 2 W. 
Nomenclature
v_{l} : Velocity of the liquid
v_{r} : The difference between the velocity of the liquid phase and the vapor phase velocity v_{l} –v_{g}
f_{1} : The force caused by an external electric potential
f_{2} : Body force generated by an electrical double layer between the liquid and solid
f_{3} : Body force generated by an electrical double layer between the liquid and vapor
k : The curvature (viewed from the vapor phase)
k : The inverse of Debye length
n^{⇀} : The unit normal vector of the interface (outward from the gas to liquid phases)
S_{h} : Any heating source including radiation and
k_{eff} : Effective thermal conductivity
q ^{‴} : At generated by the thermal source
k_{s} : Thermal conductivity coefficient of the wall
c_{i,0} : Bulk concentration of the ith ion
n_{0} : Bulk ionic concentration,
Greek Symbol
ρ_{e} : Volumetric charging of the double electrical layer
σ : Coefficient of contact surface tension
φ : External electrical potential
ψ_{1} : Potential of the wall recharge
ρ_{e} : Volumetric charging of the double electrical layer
μ^{S } : Chemical potential of ions at the contact level
References
 M. Rahmat, P. Hubert, Twophase simulations of micro heat pipes, Comput. Fluids 39 , 451–460 (2010) [Google Scholar]
 G. Mohiuddin Mala, D. Li, J.D. Dale, Heat transfer and fluid flow in microchannels, Int. J. Heat Mass Transfer 40 , 3079–3088 (1997) [CrossRef] [Google Scholar]
 R.W. Knight, D.J. Hall, J.S. Goodling, R.C. Jaeger, Heat sink optimization with application to microchannels, IEEE Trans. Compon. Hybrids Manufact. Technol. 15 , 832–842 (1992) [CrossRef] [Google Scholar]
 M. Groll, M. Schneider, V. Sartre, M. Chaker Zaghdoudi, M. Lallemand, Thermal control of electronic equipment by heat pipes, Rev. Génér. Therm. 37 , 323–352 (1998) [CrossRef] [Google Scholar]
 M.L. Berre, G. Pandraud, P. Morfouli, M. Lallemand, The performance of micro heat pipes measured by integrated sensors, J. Micromech. Microeng. 16 , 1047–1050 (2006) [Google Scholar]
 M. Fallah Abbasi, H. Shokouhmand, Experimental investigation on effect of EDL on heat transfer of micro heat pipe, Microgr. Sci. Technol. 31 , 317–326 (2019) [CrossRef] [Google Scholar]
 Y. Lei, Z. Chen, J. Shi, Analysis of condensation heat transfer performance in curved triangle microchannels based on the volume of fluid method, Microgr. Sci. Technol. 29 , 433–443 (2017) [CrossRef] [Google Scholar]
 K. Fukagata, N. Kasagi, P. Uaarayaporn, T. Himeno, Numerical simulation of gasliquid twophase flow and convective heat transfer in a micro tube, Int. J. Heat Fluid Flow 28 , 72–82 (2007) [Google Scholar]
 G. Černe, S. Petelin, I. Tiselj, Coupling of the interface tracking and the twofluid models for the simulation of incompressible twophase flow, J. Comput. Phys. 171 , 776–804 (2001) [Google Scholar]
 E. Berberović, N.P. van Hinsberg, S. Jakirlić, I.V. Roisman, C. Tropea, Drop impact onto a liquid layer of finite thickness: dynamics of the cavity evolution, Phys. Rev. E 79 , 036306 (2009) [Google Scholar]
 P. Zhang, H.H. Qiu, Investigation of the patterned surface modification on 3D vortex flow generation in a micropipe, J. Micromech. Microeng. 18 , 115030 (2008) [Google Scholar]
 L. Ren, D. Li, W. Qu, Electroviscous effects on liquid flow in microchannels, J. Colloid Interface Sci. 233 , 12–22 (2001) [Google Scholar]
 L. Gong, J.K. Wu, L. Wang, K. Chao, Periodical streaming potential and electroviscous effects in microchannel flow, Appl. Math. Mech. 29 , 715–724 (2008) [CrossRef] [Google Scholar]
 H. Nakajima, Mass Transfer, Open access peerreviewed Edited Volume (November 4th 2011) (2011) [Google Scholar]
 H. Shokouhmand, A.A.R. Bahrami, Study of the effects of electrokinetic field on heat transfer through rectangular microchannels, J. Fac. Eng. 39, 773–784 (2006) [Google Scholar]
 B.Q.J.A.W. Qian, MEMSBased Microheat Pipes (June 15th 2016) [Google Scholar]
 P. Krata, The impact of sloshing liquids on ship stability for various dimensions of partly filled tanks, TransNav, Int. J. Mar. Navigat. Saf. Sea Transp. 7 , 481–489 (2013) [CrossRef] [Google Scholar]
Cite this article as: M. Fallah Abbasi, H. Shokouhmand, M. Khayat, An investigation on effect of EDL on heat transfer of micro heat pipe with square and triangular cross section, Mechanics & Industry 21, 309 (2020)
All Tables
All Figures
Fig. 1
Schematic of the micro heat pipe function and utilizing groove to return the fluid from the evaporator to the condenser [16]. 

In the text 
Fig. 2
Contact area of two phases in real mode (above) and by Eulerian VOF method (bottom) [17]. 

In the text 
Fig. 3
Cross section of micro heat pipe (triangular and square curves type). 

In the text 
Fig. 4
Special tool designed to form a pipe to the triangular and square crosssection. 

In the text 
Fig. 5
Location and installation of thermocouple. 

In the text 
Fig. 6
Schematic of test set up [6]. 

In the text 
Fig. 7
Micro heat pipe adapter. 

In the text 
Fig. 8
Steps for preparing the water solution and ethanol. 

In the text 
Fig. 9
Temperature distributions in the axial length of 50 mm in micro pipe with. 

In the text 
Fig. 10
Meshing that used for boundary conditions. 

In the text 
Fig. 11
Thermal contour for input current. 

In the text 
Fig. 12
3D model of MHP. 

In the text 
Fig. 13
Validation of sample thermal resistance tested in this case with previous record and researches. 

In the text 
Fig. 14
Compare experimental and theoretical graph of wall temperature, water and 0.018 g ethanol as working fluid, 1 W. 

In the text 
Fig. 15
Compare experimental and theoretical graph of wall temperature, water and 0.054 g ethanol as working fluid, 1 W. 

In the text 
Fig. 16
Compare experimental and theoretical graph of wall temperature, water and 0.18 g ethanol as working, 1 W. 

In the text 
Fig. 17
Compare experimental and theoretical graph of wall temperature, water and 0.018 g ethanol as working, 4 W. 

In the text 
Fig. 18
Compare experimental and theoretical graph of wall temperature, water and 0.054 g ethanol as working, 4 W. 

In the text 
Fig. 19
Compare experimental and theoretical graph of wall temperature, water and 0.18 g ethanol as working, 1 W. 

In the text 
Fig. 20
Compare experimental and theoretical graph of wall temperature, water and 0.054 g ethanol as working, 2 W. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.