Melting Performance Enhancement of Triplex Tube Latent
Published: Last Edited:
Disclaimer: This essay has been submitted by a student. This is not an example of the work written by our professional essay writers. You can view samples of our professional work here.
Any opinions, findings, conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of UK Essays.
Thermal Storage Using FinsNanoPCM Technique
Ammar M. Abdulateef^{1}'*, Sohif Mat^{1}, Jasim Abdulateef^{2}
1 Solar Energy Research Institute, University Kebangsaan Malaysia, Bangi, Selangor, Malaysia department of Mechanical Engineering, University of Diyala, 32001 Diyala, Iraq
Latent heat thermal energy storage (LHTES) systems using phase change material (PCM) could have lower heat transfer rates during charging/discharging processes due to its low inherence of the thermal conductivity. In this study, heat transfer enhancement using internal longitudinal fins employing PCM first and nanoPCM secondly in a large triplex tube heat exchanger (TTHX) was investigated by Fluent 15 software numerically. The results showed the thermal conductivity of pure PCM (0.2 W/m.K) could be enhanced to 25% by dispersing 10% alumina (AEO3) as a nanoparticle. However, the melting time is reduced to 12% as compared with the PCM only therefore, a longitudinal finsnanoPCM technique achieved a complete PCM melting shortly (218 minutes). Consequently, the simulation results have been validated and illustrated a good agreement with the PCM and nanoPCM experimentally.
Keywords: phase change material, triplex tube heat exchanger, melting time, longitudinal fins, nanoparticle
The major emphasis associated with most of the solar devices application is the continuous power generation during cloud transients and nondaylight hours. Thermal energy storage (TES) systems especially the latent heat thermal energy storage (LHTES) systems offer possibility to store higher amounts of thermal energy in comparison with sensible heat thermal energy storage (SHTES) systems. However, most the phase change materials (PCM) that used as storage media in the LHTES systems suffers from the low thermal conductivity (0.2 W/m.K), it often leads to uncompleted melting/solidification process and significant temperatures difference within the PCM, which in some cases can cause a material failure and system overheating.
Many researchers studied the different kinds of heat exchangers used in the LHTES systems with (PCM). Among these, concentric cylinder, shell and tube, and triplex tube heat exchanger (TTHX) [1, 2]. Most of these have been proved a high efficient for minimum
ISSN: 236789921 volume. Agyenim et al. [3] have been presented a significant comparison for three experimental configurations, a concentric tube system with no fins and augmented with circular and longitudinal fins. The system with longitudinal fins gave the most performance with increasing thermal response during charging and reduced subcooling in the melt during discharging. Further, the melting performance enhancement of a small scale TTHX used in LHTES system has received a significant interest by [4, 5] where numerical and experimental investigations have been made using longitudinal fins technique only to improve the melting time of simple PCM. It can be seen, longitudinal fins are most common extended surfaces have been considered in TES systems. In addition, when a triplex tube heat exchanger (TTHX) is used, the heat transfer area is also extended to the PCM and thermal performance is enhanced respect to cylinder or shell and tube heat exchanger. On the other hand, the unloading latent thermal storage, the solidliquid interface moves away from the heat transfer surface and the heat flux decreases because of increasing the thermal resistance of the growing layer of the molten/solidified medium. This effect could be reduced by a technique of dispersing high thermal conductivity nanoparticles. The PCM melting dispersed with various volumetric concentrations of alumina (AhO3) that is heated from one side of a square enclosure is investigated numerically [6]. Wang et al. [7] improved thermal properties of paraffin wax by the addition of (TiO2) as a nanoparticle successfully without any surfactant.
The biggest challenge that is faced to investigate for both of PCM and nanoPCM was a large triplex tube heat exchanger (TTHX). Therefore, the contribution in the heat transfer rate between the PCM and the HTF are augmented using internal longitudinal fins first and dispersing a high conductivity material such as alumina (Al2O3) secondly to be formed with longitudinal fins as finsnanoPCM technique to produce the biggest demand thermal energy stored that is required for application in air conditioning systems.
The physical configurations of the TTHX model for two cases (1) pure PCM and (2) nanoPCM are elucidated in Fig. 1. It consists of inner tube, middle tube, and outer tube that have 38.1 mm, 190.5 mm, and 250 mm in radius and 3 mm thickness, respectively with eight internal longitudinal fins each one has 121 mm long and 2 mm thickness. The inner tube and middle tube are made from copper and outer tube from steel. The water is used as HTF to transfer the heat by conviction to the walls and by conduction to the PCM or nanoPCM. The heat transfer during the PCM melting process is based on the both sides heating method where the heat is supplied from both inner and outer tubes during the charging process. The minimum temperature has been required to operate the PCMLHTES system was approximately 90 Â°C.
The PCM melting numerical model is solved using Ansys Fluent 15 software based on the enthalpyporosity technique and the finite volume method [8]. The model is drawn and meshed in a two dimensions( r, 9) as well as boundary layers and zone types are defined using
ISSN: 236789922 International Journal of Theoretical and Applied Mechanics
Gambit 2.4.6 software. The grids size number of the numerical model for internal longitudinal fins was calculated to 56200 as illustrated in Fig. 2.
Fig. 1. Physical configurations of the TTHXinternal longitudinal fins.
Fig. 2. Distribution of the grids size number in the middle tube of TTHXinternal longitudinal fin.
For the numerical analysis of the thermal process, the following assumptions are made: (1) the melting is Newtonian and incompressible; (2) the flow in the melting process is laminar, unsteady with negligible viscous dissipations; (3) the thermophysical properties of the HTF and PCM are independent on the temperature; (4) the heat transfer is both of conduction and of convection controlled. The effect of natural convection during the charging process is considered by invoking the Boussinesq approximation that is valid for the density variations of buoyancy force, otherwise the effect is ignored. The density variation is defined as follow: p=Pi/(J3(TTl) + 1) (1)
Volume 2, 2017 Ammar M. Abdulateef et al. International Journal of Theoretical and Applied Mechanics ttp://www.iaras.org/iaras/journals/ijtam
where pi is the PCM density at the melting temperature at Tt and ft is the thermal expansion coefficient. The temperature distribution and viscous incompressible flow are solved by using the NavierStokes and thermal energy equations, respectively. The continuity, momentum, and thermal energy equations as follows [9].
The continuity equation: dt(p) + di(pui) = 0(2)
The momentum equation:
dt(pUi)+ dj(pui uj) = pdjj Uidip + pgt +
Si(3)
The energy equation:
dt(ph) + dt(pAH) + di(pui h) = di(kdiT) (4) where p is the density of the PCM, ut is the fluid velocity, p is the dynamic viscosity, p is the pressure, g is the gravity acceleration, k is the thermal conductivity and h is a sensible enthalpy. The sensible enthalpy equation:
T
h = href + f^CpAT(5)
The total enthalpy H equation:
H = h +AH(6)
where href is the reference enthalpy at the reference temperature Tref, Cp is the specific heat, AH is the latent heat content of thePCM that changes between zero (solid) and L (liquid), y is the liquid fraction, which is generated during the phase change between the solid and liquid state when the temperature is Tt > T > T_{s}, which can be written as:
y = AH/L y = 0 y=l
(7)
if T <TS if T>T,
Y =
TTe
if T* <T <T,
TiT_{s}
From equation (3) the source term St is:
(8)
Si = C(ly)'
where C(1 y)
y^{3}+s
2 ui
Y^{3}+£
(9)
is the "porosity function"
U
defined by Brent et al. [10]. C is a constant describes how sharply the velocity is reduced to zero when the material solidifies. This constant varies between 10^{4} and 10^{7} (10^{5} is considered), and £ is a small (0.001) to prevent division by zero.
2.3. Boundary and initial conditions
At the initial time, the PCM was in a solid state and the temperature reached to 27 oC. A constant
temperature of the tube wall represented the HTF temperature [11, 12] that was at approximately 90 Â°C.The boundary conditions as follows:
Both sides heating method: at r = rt^ T = Thtf(10)
at r = r_{m} ^ T = Thtf(11)
Initial temperature of the model: at t = 0 ^ T = Tini(12)
In case of nanoP CM, we have considere d the same conservation equations, boundary, and initial conditions mentioned above.
2.4. Thermophysical properties
Table 1 describes the thermophysical properties of materials are used [4], the thermophysical properties of the nanoPCM are calculated [13]: The density equation:
Pnpcm^{0}Pnp + (1^{0)}Ppcm(13)
The sp ec ific heat cap acity e quati on:
_ <P(pcv)
C,
p,npcm
Pnpcm
The late nt h e at equatio n :
_ (l  $)(.pL)pcm
(14)
j=(15)
npcm^{(15)}
Pnpcm
The dynamics viscosity of nanoPCM is
given by [14]:
Pnpcm= 0.983e^{(12959}Â®p_{pcm}(16)
The effective thermal conductivity of the
nanoPCM, which includes the effects of particle
size (dnp), particle volume fraction (0), and
temperature dependence as well as propertie s of
the base PCM. The particle subject to Brownian
motion is also given by [14]:
Knp) 0
Knp + 2Kpcm 2 jj^pcm
^{npCm}K_{np}+2K_{pcm}+2(K_{pcm}K_{np}) 0 P^{cm}
+
5 x 1 0 4 y_{k} g0p_{pC}mcp,_{pcm}J_{Pnpdnp} f(T' 0)
(17)
where B is the Boltzmann constant (1.381 x 1023 J/K) and yk = 8.4407(1000)10â„¢4. f(T, 0) = (2.8217 x 1020+ 3.917 x 103) + (3.0669 x 1020 3.91123 x
^{T}ref
10^{3})(18)
where Tref is the reference temperature = 273 K.
We have evaluated in the equation (17), the effects of nanoparticle diameter (dnp = 20 nm), nanoparticle volume fraction (0 = 10%), and the reference temperature (Tref = 237 K).
ISSN: 23678992
3
Volume 2, 2017
Ammar M. Abdulateef et al.
Table 1. Thermophysical properties of PCM, copper, and alumina (AI2O3).
Properties 
PCM (RT82) 
Copper 
A^O3 
Density, solid, ps (Kg/m3) 
950 
8978 
3600 
Density, lquid, pi (Kg/m3) 
770 
 
 
Specific heat, Cpi , Cps (J/kgK) 
2000 
381 
765 
Latent heat of fusion, L (J/kg) 
176000 
 
 
Dynamic viscosity, p (kg/m.s) 
0.03499 
 
 
Melting temperature, Tm (K) 
350.15 358.15 
2345 

Thermal conductivity, K fW/m.K) 
0.2 
387.6 
36 
Thermal expansion coefficient, (1/K) 
0.001 
 
 
A schematic diagram of the LHTES system apparatus is illustrated in Fig. 3. The middle tube of TTHX is filled with 100 kg PCM first. The present numerical model for PCM and nanoPCM has been validated experimentally with PCM as illustrate in Fig. 4. A comparison resulted was not exceeded in percentage errors of 3% and showed a good agreement with an experimental test for two cases. Moreover, the average temperature of the PCM was 27 Â°C when melting process started and the HTF charging temperature by both sides heating method [4] was 90 Â°C with an experimental mass flow rate 37.5 L/min.
Fig. 3. Schematic diagram of experimental apparatus of LHTES system, which includes; 1. Evacuated tube solar collectors (ETSC), 2. Flow meter, 3. Triplex tube heat exchanger (TTHX), 4. Thermocouple Jtype, 5. Sensor (water), 6. Internal longitudinal fin, 7. Pressure vessel tank, 8. Pump, 9. Data acquisition, 10. Computer, 11. Water storage tank, 12. Electrical heater, 13. Pipes, 14. Valve two ways, 15. Valve three ways.
ISSN: 236789924
International Journal of Theoretical and Applied Mechanics
Time (min)
Fig. 4. Validation of an experimental and numerical model
The isothermal contours of the PCM in TTHX with internal fins at different times (10, 60, 120, and 247 min) are elucidated in Fig. 5. firstly, heat transfer occurred between the hot wall of the tube and solid surface of the PCM by conduction, which dominated the melting process at the early stage and caused a very thin layer of the liquid that is surrounded the longitudinal fin surface and hot wall of the tube while the rest of the PCM remained solid without any phase change because of the effects of natural convection were limited. After 10 minutes, small convection cells are formed between the fins wall and subsequently expanded to the middle tube. Over time, cells convection emerged and facilitated the formation of the large convection cells at 60 minutes that are expanded to the bottom part of tube at 120 minutes because heat transfers by fins. The full PCM melting was accomplished at 247 minutes.
The thermophysical properties of the nanoPCM with various volumetric concentrations of the alumina (AfO3) are calculated using equations (1317). It is found that, the specific heat and latent heat of the nanoPCM are lower than the pure PCM whereas the thermal conductivity and dynamic viscosity of the nanoPCM are higher than the pure PCM, see Table 2. This variation in
Volume 2, 2017
Liquid fraction
International Journal of Theoretical and Applied Mechanics
Ammar M. Abdulateef et al.http://www.iaras.org/iaras/journals/ijtam
the thermal conductivity and dynamic viscosity agree well with the results that reported in [6]. Moreover, augmenting the alumina nanoparticle (AhO3) volume concentrations caused to reduce the PCM melting time, see Fig. 6. Consequently, the PCM with 10% alumina (AhO3) is considered
sssssssssssssssasssss
10 min60 min
120 min247 min
Fig. 5. Isothermal contours of the PCM in TTHX longitudinal fins.
Table 2. Variation of the thermal conductivity and dynamic viscosity of nanoPCM.
Volumetric concentration <H%) 
Thermal conductivity k (W/m.K) 
Dynamic viscosity g (kg/m.s) 
Simple PCM 
0.2 
0.03499 
NanoPCM (1% M2O3 ) 
0.206 
0.0121161 
NanoPCM (4% M2O3) 
0.225 
0.0485 
NanoPCM (7% M2O3) 
0.245 
0.084812 
NanoPCM (10% M2O3 ) 
0.265 
0.121161 
1.2
0100200300
Time (min)
Fig. 6. Effect of the nanoparticle concentrations.
4.2.1. Nanoparticleinternal fins technique
The isothermal contours of the finsnanoPCM technique in TTHX at different times (10, 60, 120, and 218 min) are shown in Fig. 7. A significant reduction in time was observed by dispersing 10% nanoparticle to the PCM when the absorbed energy was stored to the required load under the effects of both sides heating method, which is augmented the conduction heat transfer rate. Therefore, the full melting of the PCM is completed at 218 minutes. Consequently, the nanoparticle plays a significant role in the melting rate enhancement where the thermal conductivity of simple PCM (0.2 W/m.K) could be enhanced to 25% significantly that is caused to increase the conduction heat transfer.
10 min60 min
120 min218 min
Fig. 7. Isothermal contours of the finsnanoPCM technique.
Figure 8 illustrates liquid fraction vs. melting time for the PCM and nanoPCM in TTHX internal longitudinal fins. As shown, the PCM melting time is reduced using nanoPCM to 12% as compared to the PCM only. The PCM melting retardation was reduced because of augmenting the thermal conductivity of PCM effectively.
ISSN: 23678992
5
Volume 2, 2017
Liquid fraction
Ammar M. Abdulateef et al.
International Journal of Theoretical and Applied Mechanics
Consequently, the model of finsnanoPCM is considered the most efficient technique to achieve the PCM melting shortly (218 min).
Fig. 8. Liquid fraction vs. melting time for the PCM and nanoPCM in TTHXintemal longitudinal fins.
Heat transfer enhancement for a large triplex tube heat exchanger (TTHX) has been represented the biggest challenge in LHTES system. The results showed the thermal conductivity of simple PCM (0.2 W/m.K) could be enhanced to 25% by dispersing 10% alumina and the melting time is reduced to 12% as compared with the PCM only. Consequently, the model of finsnanoPCM has been considered the most efficient technique based on both sides heating method to achieve the PCM melting shortly (218 min). However, the numerical results have validated and showed a good agreement with the PCM and nanoPCM experimentally.
BBoltzmann constant (J/K)
Cmushy zone constant (kg/m^{3}s)
Cpspecific heat (J/kg.K)
gi gravity acceleration in the idirection (m/s^{2})
Henthalpy (J/kg)
HTFheat transfer fluid
Llatent heat fusion (J/kg)
kthermal conductivity (W/m.K)
ppressure (Pa)
Tmmelting temperature (^{o}C or K)
uvelocity component (m/s)
Simomentum source term in the idirection (Pa/m)
pfluid density (kg/m^{3})
yliquid fraction
Pthermal expansion coefficient (1 /K)
Zcorrection factor
The authors gratefully appreciate a financial support that provided by Solar Energy Research Institute(SERI),University Kebangsaan Malaysia (UKM), Malaysia.
 H. Niyas, P. Muthukumar, Performance analysis of latent heat storage systems, International Journal of Scientific & Engineering Research 4 (2013) 22295518.
 Y.L. Jian, Numerical and experimental investigation for heat
transfer in triplex concentric tube with phase change material for thermal energy storage, Solar Energy 32
 85977.
 F. Agyenim, P. Eames, M. Smyth, A comparison of heat transfer enhancement in a medium temperature thermal energy storage heat exchanger using fins, Solar Energy 83
 15091520.
 S. Mat, A.A. AlAbidi, K. Sopian, M.Y. Sulaiman, A.T. Mohammad, Enhance heat transfer for PCM melting in triplex tube with internalexternal fins, Energy Conversion and Management 74 (2013) 223236.
 A.A. AlAbidi, S. Mat, K. Sopian, M.Y. Sulaiman, A.T. Mohammad, Heat transfer enhancement for PCM thermal energy storage in triplex tube heat exchanger, Heat Transfer Engineering, vol. 37, pp. 705712, 2016.
 A.V. Arasu, A.S. Mujumdar, Numerical study on melting of paraffin wax with Al2O3 in a square enclosure, International Communications in Heat and Mass Transfer 39 (2012) 816.
 J. Wang, H. Xie, Z. Guo, L. Guan, Y. Li, Improved thermal properties of paraffin wax by the addition of TiO2 nanoparticles, Applied Thermal Engineering (2014) 17.
 S.V. Patankar, Numerical heat transfer and fluid flow, McGraw Hill, New York, 1980.
 A.A.R. Darzi, M. Farhadi, K. Sedighi, Numerical study of melting inside concentric and eccentric horizontal annulus, Appl Math Model 36 (2012) 40804086.
 A.D. Brent, V.R.Voller, K.J. Reid, Enthalpyporosity technique for melting convectiondiffusion phase change: application to the melting of a pure metal, Numer Heat Transfer 13 (1988) 297318.
 C. Guo, W. Zhang, Numerical simulation and parametric study on new type of high temperature latent heat thermal energy storage system, Energy Convers Management 49 (2008) 27919.
 M.J. Hosseini, A.A. Ranjbar, K. Sedighi, M. Rahimi, A combined experimental and computational study on the melting behavior of a medium temperature phase change storage material inside shell and tube heat exchanger, International Communications in Heat and Mass Transfer 39 (2012) 14161424.
[ 1 3 ] A.P. Sasmito, J.C. Kurnia, A.S. Mujumdar, Numerical evaluation of laminar heat transfer enhancement in nanofluid flow in coiled square tubes, Nanoscale Research Letters 6 (2011) 376.
[14] R.S. Vajjha, D.K. Das, PK. Namburu, Numerical study of fluid dynamic and heat transfer performance of Al_{2}O_{3} and CuO nanofluids in the flat tubes of a radiator, International Journal of Heat Fluid Flow 31 (2010 ) 613621.
ISSN: 23678992
6
Volume 2, 2017
Cite This Essay
To export a reference to this article please select a referencing stye below: