Energy Conversion and Management: Sciencedirect
Energy Conversion and Management: Sciencedirect
Energy Conversion and Management: Sciencedirect
A R T I C LE I N FO A B S T R A C T
Keywords: A key challenge in the deployment of practical latent heat energy storage systems employing phase change
Latent heat energy storage (LHES) materials is the inherent low thermal conductivity of these materials. The present research is motivated by the
Phase change material (PCM) need to intensify the buoyancy-driven convection flow in the phase change material to enhance the thermal
Perforated fin performance of the system. In this paper, for the first time, the effect of applying perforated fins on the thermal
Natural convection
performance enhancement of a vertical shell and tube latent heat energy storage heat exchanger is experi-
mentally investigated and the results are compared with those of the unfinned and solid finned heat exchangers
as the base cases. Lauric acid as the phase change material is placed in the shell side and the water is passed
through the inner tube. The shells of the heat exchangers were made of transparent Plexiglas tubes to enable the
visual comparison of the melting processes. The fins and tubes were made of copper. The melting process of
phase change material is studied under different inlet water flow rates (0.5 and 1 l/min) and temperatures (55
and 65 °C). The experimental results showed that the time-averaged Nusselt number of the perforated finned
heat exchanger is about 30% higher than that of the solid finned heat exchanger due to the minor hindering
effect of the perforated fins on the development of the convection flows. Moreover, the total melting time of the
perforated finned heat exchanger is about 7% lower than that of the solid finned heat exchanger.
1. Introduction batteries [12,13], refrigeration [14] and waste heat recovery [15].
However, a significant downside of PCMs is their low thermal con-
The continuous increase in energy demand and global warming due ductivity, which decreases the heat transfer rate and restricts the de-
to the greenhouse gas emissions have motivated intensive research for ployment of PCMs in large-scale LHTES units [16]. In order to enhance
efficient use of energy and development of energy storage systems [1]. the thermal performance of the LHTES, many techniques have been
Thermal energy storage (TES) which stores heat in a material and re- established to improve the thermal conductivity of PCMs such as using
leases it when it is needed is one of the efficient techniques to reduce fins [17,18], porous structures with high thermal conductivity [19,20],
the gap between energy supply and demand. There are many practical heat pipes [21] and nanoparticles [22].
applications for TES including solar thermal systems, waste heat re- Among the above-mentioned enhancement techniques, the utiliza-
covery, district heating and cooling, thermal power plants and build- tion of high conductivity fins is one of the most efficient, cost-effective
ings [2]. There are three kinds of TES: sensible heat storage, latent heat and dependable methods for heat transfer augmentation in LHTES
storage, and reversible thermochemical reaction. Among different types systems. However, fins may have an undesirable influence on the de-
of TES, the use of phase change materials (PCMs) for latent heat velopment of natural convection flows in the liquefied PCM. The re-
thermal energy storage (LHTES) systems is receiving growing attention striction of natural convection flow suppresses the positive effect of fins
due to the absorb and release of thermal energy at an almost constant on heat transfer enhancement. Therefore, the geometric parameters of
temperature during the solid–liquid phase change processes [3]. the fins have been the subject of numerous studies to develop the
PCMs have a broad variety of applications such as solar heating LHTES systems with higher efficiency, lower weight, and smaller size
systems [4,5], temperature regulation of photovoltaic panels [6,7], [23]. Shell and tube heat exchangers (HXs) are attractive candidates for
buildings [8,9], electronics cooling [10,11], temperature control of the integration with practical LHTES units due to their simplicity in design,
⁎
Corresponding author.
E-mail address: [email protected] (B. Kamkari).
https://doi.org/10.1016/j.enconman.2020.112679
Received 23 November 2019; Received in revised form 1 March 2020; Accepted 2 March 2020
0196-8904/ © 2020 Elsevier Ltd. All rights reserved.
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
good effectiveness and ease of manufacturing. Hence, the focus of this tube LHTES unit with annular fins. It was found that by adding the
investigation is to study the solid–liquid phase change rate in the finned number of fins, the conduction heat transfer through the PCM is im-
shell and tube HXs. There are several studies in the literature on the proved but the formation of natural convection flows is limited due to
effect of fin geometry including longitudinal fins [24–33], pin fins [34], the decrease in fin spaces. Therefore, an optimum fin number was in-
helical fins [35], and annular fins [36–39] on the thermal behavior of troduced beyond which any increase reduced the heat transfer rate.
the PCM in the horizontally or vertically oriented shell and tube HXs. Parsazadeh and Duan [38] numerically investigated the melting process
Darzi et al. [26] numerically investigated the influence of adding of nano-PCM in a vertical shell and tube HX with angled annular fins. It
the number of longitudinal fins on the melting and solidification times was noticed that the fins with upward angles are more favorable to
of n-eicosane in a horizontal shell and tube HX. It was concluded that improve the heat transfer rate due to the natural convection in-
both the melting and solidification times reduce when the fins are im- tensification. Also, the results indicated that the thermal conductivity
plemented. However, the fins were found to be more effective during enhancement with nanoparticles was not able to overcome the sup-
the solidification process. This was attributed to the hampering effect of pression of natural convection caused by higher viscosity of the nano-
fins on the evolution of buoyancy-driven flows during the melting PCM which leaded to a longer melting time. Shahsavar et al. [39] nu-
process. The same thermal behavior was also reported by Kamkari and merically evaluated the effects of fin arrangement and diameter on
Shokouhmand [27] for the melting of PCM in a finned rectangular melting and solidification processes in a vertical finned shell and tube
LHTES enclosure. Yuan et. al [28] studied the impact of the installation LHTES system. It was found that the optimized non-uniform fin array
angle of the longitudinal fins on the melting of PCM in a horizontal reduces the melting time by 23.9% compared to the uniform fin dis-
double pipe LHTES unit. It was concluded that the angle of fins should tribution. Also, increasing the fin diameter beyond a certain value
be selected so that it minimizes the obstruction effect of fins on the showed an adverse effect due to suppressing the natural convection
development of natural convection flows. Sodhi et al. [29] numerically effect.
analyzed the melting and solidification of PCM in horizontal conical Based on the above literature review, it is found that applying the
shell-based LHTES systems. The optimum cone angle was obtained and annular fins on the inner tube of a vertical shell and tube LHTES sys-
also the effects of longitudinal fins with different geometries on the tems effectively improves the heat transfer rate. However, the fins re-
charging and discharging process were investigated. Mahdi and Nsofor strict the development of convection flows inside the liquid PCM.
[30] numerically investigated the solidification of PCM in a horizontal Therefore, the conflict between the enhancement in heat diffusion
triplex-tube LHTES system by applying longitudinal fins and nano- through the PCM and the simultaneous decrease in the buoyancy effect
particles and validated with experimental data. It was proposed that the due to adding fins should be considered by the designers through the
utilization of fins alone results in a better enhancement than applying selection of the optimum fin geometry to maximize the fin effect.
either nanoparticles alone or both of the nanoparticles and fins. The Some previous studies have already shown that the use of perfo-
solidification time was reduced by up to 55% using fins alone, 8% using rated fins instead of solid fins may enhance the single-phase convective
nanoparticles alone and 30% using the combination of fins and nano- heat transfer coefficient [40–42]. Awasarmol and Pise [43] experi-
particles. In another study, Mahdi et al. [33] numerically investigated mentally investigated the natural convection heat transfer enhance-
the simultaneous charging and discharging processes in a triplex-tube ment from perforated rectangular fin arrays using air as the heat
heat exchanger. A comparison between the fin system and that con- transfer fluid. It was concluded that perforated fins can increase the
taining nanoparticles showed that the insertion of fins with the re- heat transfer coefficient by up to 32% while decreasing the fin mass by
commended structure is more favorable for achieving higher thermal 30%. The same concept can be used to improve the thermal perfor-
performance. mance of LHTES heat exchangers, which has not been investigated
Besides the use of longitudinal fins, some researches have been before. The motivation of the present study is the enhancement of the
conducted to enhance the thermal performance of shell and tube LHTES convective heat transfer coefficient between the liquid PCM and the
systems using different fin geometries including annular, pin and he- heat transfer surface by using the perforated annular fins instead of the
lical shapes [35]. Tay et al. [34] numerically compared the effective- solid fins. To the best of authors’ knowledge, this is the first time that
ness of annular and pin fins in LHTES systems. The analysis showed that the effect of applying perforated fins on the melting process of PCM has
using the annular fins is a more effective heat transfer enhancement been investigated. For this purpose, an experimental setup has been
method than applying the pin fins. Ogoh and Groulx [36] investigated constructed to measure the instantaneous heat transfer and energy
the melting of paraffin wax in a vertical shell and tube HX with dif- storage rates in unfinned, solid finned and perforated finned HXs during
ferent number of annular fins. Yang et al. [37] conducted a numerical the melting process of the PCM under the different working conditions.
investigation on the melting process of paraffin in a vertical shell and The results of the perforated finned HX are compared with those of the
2
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
unfinned and solid finned HXs as the based cases. The findings provide Sufficient time was given to each liquid layer to become solid and then
useful information on the performance improvement of the LHTES the next liquid layer was added to the annulus. The same procedure was
systems. The total weight of the perforated fins is 16% less than that of repeated to fill the annulus with PCM.
the solid fins, which leads to a lighter and more cost-effective HX while Prior to the onset of each experiment, in order to provide a constant
enhancing the thermal performance of the system initial thermal condition, cold water at 22 °C flowed from the cooling
circulator system into the tube of the HX until all thermocouples in the
2. Experimental setup and procedure PCM reached the HTF inlet temperature. The HTF inlet temperature and
water flow rate were set at desired values. The tests were initiated by
2.1. Experimental setup flowing the hot water from the constant temperature water bath to the
HX. The temperatures of the inlet/outlet HTFs and PCM were recorded
An experimental setup was built to measure the thermal perfor- by the data logger at a time interval of 2 sec. Melting photographs were
mance of different shell and tube LHTES HXs. Fig. 1 shows a schematic taken by removing the front half of the insulation at predefined time
diagram and a photograph of the experimental setup used in this study. intervals of 15 min. The tests proceeded until no solid PCM was seen in
Three different shell and tube HXs: unfinned, solid finned and perfo- the annulus.
rated finned HXs were fabricated to be used as the LHTES systems The aforementioned experimental procedure was used to investigate
(Fig. 2). The shells of the HXs were made of the transparent Plexiglas the heat transfer rate and energy storage in the unfinned, solid finned
tubes to provide the visual observation of the melting process. The and perforated finned HXs. The HXs were tested under HTF inlet tem-
inner and outer diameters of the shell were 60 and 80 mm, respectively. peratures of 55 and 65 °C and HTF flow rates of 0.5 and 1 l/min. At HTF
A copper tube was mounted concentrically inside the shell. The inner inlet temperature of 55 °C, the flow rates of 0.5 and 1 l/min are asso-
and outer diameters of the tube were 8 and 10 mm, respectively. Fifteen ciated with Reynolds number of 2600 and 5200, respectively. At HTF
copper fins of 50 mm diameter and 1 mm thickness were welded to the inlet temperature of 65 °C, the flow rates of 0.5 and 1 l/min are asso-
inner tube and spaced equally to fabricate the solid finned HX. The ciated with the Reynolds numbers of 3000 and 6000, respectively. The
perforated finned HX consists of fifteen perforated copper fins with the non-dimensional numbers of Rayleigh and Stefan numbers corre-
same size as the solid fins. Every perforated fin has 24 circular holes sponding to the charging temperatures are calculated using the fol-
with a diameter of 4 mm. The height of all HXs is 500 mm. Water as the lowing formulas [47]:
heat transfer fluid (HTF) flows through the inner tube and the PCM fills
gβ (Tw − Tm) L3
the annulus between the shell and tube. Ra =
να (1)
Fig. 3 shows the distribution of thermocouples used to measure the
temperature of the PCM. Six T–type thermocouples were placed inside Cp, s (Tm − To) + Cp, l (Tw − Tm)
the PCM at two different heights of the HX. Three thermocouples were Ste∗ =
hSL (2)
placed at each height at three different radial distances from the ver-
tical axis of the HX. In order to prevent the displacement of thermo-
couples during the melting and solidification processes, a plastic needle
with a diameter of 1.3 mm was used as a support in each row of the
thermocouples. Two thermocouples were mounted at the inlet and
outlet of the tube to measure the water temperature. All thermocouples
were calibrated with an accuracy of ± 0.1 °C and connected to the
channels of the multiplexer. The multiplexer was the interface between
the data acquisition device and the thermocouples. The computer re-
corded the temperatures of the thermocouples at predefined time in-
tervals. Flexible silicone tubes and valves were used to connect the HX
to the cooling circulator and constant temperature hot water bath
systems.
In order to minimize the heat loss from the HX to the ambient air, a
closed-cell Ethylene Propylene Diene Monomer (EPDM) insulation tube
with a thickness of 60 mm was wrapped around the HX. The tem-
perature of the insulation surface was measured by an infrared ther-
mometer to evaluate the heat loss from the HX. The maximum tem-
perature difference between the outer surface of the insulation and
ambient air at the end of the melting process was less than 1 °C. This
small temperature difference indicates that the heat loss from the HX is
negligible and can be ignored in the calculations.
The PCM used in this experimental study is lauric acid (dodecanoic
acid), which is a fatty acid, with 99% purity (Panreac Chemical). It has
favorable characteristics such as, non-toxicity, good chemical stability
as well as medium-temperature phase change, which qualify it as a
promising PCM for medium-temperature thermal storage applications
[44,45]. The thermophysical properties of lauric acid are listed in
Table 1.
The lauric acid in the solid state, which is a white powder, was Fig. 1. (a) Schematic diagram, (b) Photograph of the experimental setup. (1)
heated in an oven to become liquid. In order to avoid the formation of Heat exchanger (2) Data acquisition system (3) Personal computer (4) PID
air pockets in the PCM, the liquid PCM was poured layer-by-layer into temperature controller (5) Constant temperature hot water bath (6) Cooling
the annulus of the HXs. The height of each layer was about 10 mm. circulator (7) Circulation Pump (8) Channel multiplexer.
3
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
Fig. 2. Photograph and the schematic views of the unfinned, solid finned and perforated finned HXs: (a) Photograph of the HXs, (b) Schematic views and dimensions
of the HXs (all dimensions are in mm).
4
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
3. Data reduction
The instantaneous heat transfer rate from the HTF to the PCM is
equal to the decrease in the sensible heat of the HTF flowing through
the tube. By measuring the temporal temperatures of the inlet and
outlet HTFs, the heat transfer rate (Q̇ (t ) ) is obtained by the following
equation:
where ṁ is the mass flow rate of the HTF, Cp, w is the specific heat
Fig. 3. Distribution of the thermocouples. capacity of the HTF, Tin and Tout are the HTF inlet and outlet tempera-
tures, respectively.
Table 1 To further evaluate the effect of applying fins on the heat transfer
Thermophysical properties of lauric acid and copper [45,46]. rate from the HTF to the PCM, the time-averaged heat transfer rate over
the melting time (Q̇ (t ) ) is calculated for different HXs by using the
Property Lauric acid Copper
following equations:
Specific heat capacity solid/liquid (kJ/kg K) 2.18/2.39 0.381 ttotal
1
Thermal conductivity solid/liquid (W/m K) 0.16/0.14 387.6
Q¯̇ = ∫ Q̇ (t ) dt
Density solid/liquid (kg/m3) 940/885 8978 ttotal 0 (4)
Melting temperature range (°C) 43.5–48.2
Latent heat of fusion (kJ/kg) 187.21
Kinematic viscosity (m2/s) -6 where ttotal is the total melting time of the PCM.
6.2 × 10
Prandtl 93.8 In order to evaluate the heat transfer rate per unit area of the heat
transfer surface (tube and fins surfaces), the time-averaged heat flux
(Q̄″) is defined as follows:
In the above equations, considering the high thermal conductivity of
the copper material and small temperature difference between the HTF Q¯̇
Q¯ ″ =
inlet and outlet temperatures, the heat transfer surface temperature (Tw ) A (5)
is estimated as the average temperature of the HTF inlet and outlet
where A is the total heat transfer area of the HX. For the finned HXs, A
temperatures. The Rayleigh numbers corresponding to charging
is equal to the sum of the areas of the tube and fins.
Table 2
Compactness factors and potential of energy storage ratios for the finned HXs.
HX type Compactness factor Potential of energy storage ratio at Tin = 55 °C Potential of energy storage ratio at Tin = 65 °C
Unfinned HX 1 1 1
Solid finned HX 97.9% 99.2% 99.1%
Perforated finned HX 98.1% 99.3% 99.2%
5
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
⎡ π (Df − Dt )
2 πDp2 ⎤ compared to the ultimate cumulative energy of the HX measured by Eq.
A = πDt L + ⎢ −n× × 2N (5).
⎣ 4 4 ⎥ ⎦ (6)
Qstored = mpcm Cp, s (Tm − Ti ) + mpcm hsl + mpcm Cp, l (Tf − Tm)+
where Df , Dt and Dp are diameter of fin, outer diameter of tube and the
mshell Cp, shell (Tf − Ti ) + mcopper Cp, copper (THTF − Ti ) (9)
perforation diameter. n and N are the number of perforations on each
fin and fin number, respectively. where m is the mass of the PCM in the LHTES HX, Cp, s and Cp, l are the
The time-averaged Nusselt number is calculated using the following specific heat capacities of the solid and liquid PCM, respectively. Ti and
equation [24]: Tf are the initial and final temperatures of the PCM, respectively. mshell
ttotal ttotal is the mass of the shell, mcopper is the mass of the copper tube and fins,
1 1 Q̇ (t ) Lc
¯ =
Nu
ttotal
∫ Nu (t ) dt = ttotal ∫ A (Tw − Tm) KPCM
dt Cp, shell and Cp, copper are the specific heat capacities of the Plexiglas and
0 0 (7) copper, respectively. In Eq. (5), the first term on the right-hand side
where the characteristic length is Lc = Ds − Dt and Ds is the inner expresses the sensible heat gained by the solid PCM to reach its melting
diameter of the shell. point. The second term represents the latent heat storage of the PCM
during the phase change process. The sensible heating of the liquid PCM
is expressed by the third term. The fourth and the fifth terms estimate
3.2. Energy storage
the sensible heat gained by the shell of the HX, copper fins, and tube,
respectively. The calculations showed that the difference between the
Based on the energy conservation law, the energy stored in the PCM
maximum theoretical energy storage (Eq. (5)) and the measured value
is equal to the energy input by the HTF minus the heat loss from the HX
(Eq. (4)) is less than 6% proving the confidence in energy storage cal-
to the ambient air. Considering that the heat loss from the HX is neg-
culations.
ligible, the cumulative energy storage at any time (Qstored (t ) ) is obtained
by summing the input power over the desired time interval as follows:
3.3. Fin effectiveness
t
Qstored (t ) = ∫ Q̇ (t ) dt
The overall fin effectiveness is another significant factor to quantify
0 (8)
the HX thermal performance. Two different definitions for the fin ef-
To confirm the validity of energy storage calculations, the maximum fectiveness are presented: classical fin effectiveness and perforated fin
theoretical energy storage capacity of the HX determined by Eq. (4) is effectiveness. Classical fin effectiveness is calculated as the ratio of
Fig. 4. Sequential photographs of the melting process in unfinned, solid finned and perforated finned HXs.
6
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
Q¯̇finned HX
εfin,classical =
Q¯̇unfinned HX (10)
where Q¯̇finned HX and Q¯̇unfinned HX represent the time-averaged heat
transfer rates from the finned and unfinned HXs, respectively.
In order to quantify the superiority of perforated fins over the solid
fins, the perforated fin effectiveness is defined as the improvement in
time-averaged heat transfer rate in the presence of perforated fins re-
lative to the solid fins as follows.
Q¯̇perforated fin HX
εfin,perforated =
Q¯solid
̇ fin HX (11)
where Q¯̇perforated fin HX and Q¯solid
̇ fin HX are the time-averaged heat transfer
rates from the perforated and solid finned HXs, respectively.
7
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
enhance the natural convection flows while in the solid finned HX, the
solid fins hinder the upward convection flow.
8
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
accumulates more sensible heat at the upper part of the HX. Fig. 7(b)
illustrates the heat transfer rates in the unfinned and solid finned HXs
when the HTF inlet temperature is 65 °C. The average values of the heat
transfer rates in the unfinned and solid finned HXs, during the con-
vection dominated melting, are 0.01 and 0.04 W, respectively which
are lower than those of Fig. 7(a). As it is expected, the less temperature
difference between the HTF and PCM leads to a lower heat transfer rate.
Fig. 8 compares the heat transfer rates between the solid and per-
forated finned HXs at HTF inlet temperatures of 65 and 55 °C. Fig. 8(a)
shows that during the initial stage of the melting process, the heat
transfer rates of the solid finned HX are greater than those of the per-
forated finned HX. The superiority of the solid fins over the perforated
fins during the early stage of the melting process is attributed to the
higher rate of heat conduction from the fin surface to the PCM due to
the larger heat transfer area of the solid fins compared to the perforated
ones. However, during the convection dominated melting period, the
heat transfer rate of the perforated finned HX surpasses the solid finned
HX due to the higher intensity of natural convection heat transfer
caused by the flow of melted PCM through the perforations. The pre-
dominance of the perforated finned HX continues in the third stage of
Fig. 7. Comparison of heat transfer rates and temperature histories between the
unfinned and solid finned HXs at different HTF inlet temperatures: (a) 65 °C, (b)
55 °C.
Fig. 7 compares the heat transfer rates between the unfinned and
solid finned HXs at two different HTF inlet temperatures of 65 and
55 °C. Fig. 7(a) clearly shows that the heat transfer rates in the solid
finned HX are higher than those of the unfinned HX due to the larger
heat transfer area. This difference is more pronounced during the early
stage of the melting process when the conduction heat transfer is
dominant. As an instance, the values of the heat transfer rates in the
unfinned and solid finned HXs after 5 min from the onset of the melting
process are 0.03 and 0.18 W, respectively. During the convection
dominated melting period, the average values of the heat transfer rates
in the unfinned and solid finned HXs are 0.02 and 0.07 W, respectively.
This thermal behavior indicates that, although the fins greatly improve
the heat transfer rate over the whole melting process, they are more
effective during the early stage of melting when the conduction is the
main heat transfer mechanism. Moreover, rapid increases in tempera-
ture histories are observed sooner in the solid finned HX compared to
the unfinned HX implying the faster melting rate in the former one
compared to the latter one. Also, the final temperatures of thermo-
couples located in section B of the unfinned HX show higher values than Fig. 8. Comparison of heat transfer rates and temperature histories between the
those of the solid finned HX revealing that the unfinned HX solid and perforated finned HXs at different HTF inlet temperatures: (a) 65 °C.
(b) 55 °C.
9
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
the melting. Temperature histories also show that initially, tempera- transferred from the HTF to the PCM. The increase in HTF inlet tem-
tures of the PCM in the solid finned HX increase at faster rates com- perature accelerates the initial conduction melting of PCM followed by
pared to those of the perforated finned HX. Thereafter, by increasing sooner onset of natural convection flow in the annulus leading to a
the liquid fraction, the temperature values of the perforated finned HX shorter melting time.
surpass the solid finned HX indicating that the heat is transferred at a Increasing the HTF flow rate enhances the convective heat transfer
faster rate to the solid–liquid interface due to the higher intensity of coefficient between the HTF and the internal surface of the tube, which
natural convection flows in the liquid PCM of the former compared to contributes to increasing the overall heat transfer rate to the PCM. In
the latter. For melting with HTF inlet temperature of 55 °C similar addition, increasing the HTF flow rate decreases the temperature dif-
trends to those presented in Fig. 8(a) are observed in Fig. 8(b). ference between the inlet and outlet HTF temperatures, which provides
Fig. 9 compares the time-averaged heat transfer rates between the a more uniform temperature distribution on the tube surface with a
finned and unfinned HXs at different HTF inlet temperatures and flow higher magnitude leading to a faster melting rate of the PCM. However,
rates. It shows that the time-averaged heat transfer rate greatly im- the decrease in melting time with the increase in the HTF flow rate is
proves by adding the solid fins to the inner tube of the HX. Using the not at the same degree of significance as the increase of the HTF inlet
perforated fins also further increases the heat transfer rate compared to temperature.
the solid finned HX. This reveals that a vertical shell and tube LHTES Table 3 quantitatively presents the melting time reductions in the
system with perforated fins can perform better than a solid finned unfinned and solid finned HXs under the different working conditions.
LHTES system while having less weight. Moreover, it is seen that the For comparison, the unfinned HX with HTF inlet temperature of 55 °C
time-averaged heat transfer rate grows by increasing the HTF inlet and a flow rate of 0.5 l/min is considered as the base case. The highest
temperature or flow rate because of the enhancement in natural con- melting time reduction is 90% obtained by the solid finned HX with a
vection. However, the effect of raising the HTF inlet temperature is HTF inlet temperature of 65 °C and a flow rate of 1 l/min. It is inter-
more remarkable than increasing the flow rate. esting to note that the enhancement in melting time reduction with the
Fig. 10 compares the time-averaged heat fluxes between the finned increase in HTF inlet temperature or flow rate is more remarkable for
and unfinned HXs at different HTF inlet temperatures and flow rates. the unfinned HX when the values are compared to those of the solid
Unlike Fig. 9, showing the increase in time-averaged heat transfer rate finned HX. For instance, at a flow rate of 1 l/min, the value of melting
by adding the solid fins to the HX, Fig. 10 shows that the time-averaged time reduction for the unfinned HX varies from 24% to 60% (36% in-
heat flux of the solid finned HX is less than that of the unfinned HX. The crease) when the HTF inlet temperature rises from 55 to 65 °C, while for
contradiction between the increase in heat transfer rate and decrease in the solid finned HX at the same working condition, the melting time
the heat flux by adding the solid fins to the HX is related to the re- reduction changes slightly from 82% to 90% (8% increase). This
striction of natural convection flows in the annulus imposed by the solid thermal behavior implies that the increment in heat transfer rate with
fins. However, the time-averaged heat flux of the perforated finned HX the increase in HTF inlet temperature or flow rate is more remarkable in
is about 35% higher than that of the solid finned HX due to the minor the unfinned HX compared to the solid finned HX.
hindering effect of the perforated fins on the development of the con- Fig. 13 compares the melting time of PCM between the solid and
vection flows. This reveals that per unit area of the heat transfer surface perforated finned HXs. The figure shows that for all the HTF inlet
of the perforated finned HX contributes more effectively in heat transfer temperatures and flow rates, melting times of the PCM in the perforated
enhancement compared to the solid finned HX. finned HX are less than those of the solid finned HX. The decrease in
Fig. 11 shows the time-averaged Nusselt numbers of the unfinned, melting time is more pronounced at higher HTF inlet temperature while
solid finned and perforated finned HXs under different HTF inlet tem- the effect of increasing the HTF flow rate on decreasing the melting
peratures and flow rates. The hampering effect of the solid fins on the time is less significant.
evolution of natural convection flows reduces the convective heat At a HTF inlet temperature of 65 °C, the decrease in melting time by
transfer coefficient between the heat transfer surface and liquid PCM using the perforated fins instead of the solid fins are 9.5% and 8.7% for
leading to a decrease in the Nusselt number. However, it is worth the HTF flow rates of 0.5 and 1 l/min, respectively. However, when the
mentioning that the increase in the heat transfer area outweighs the
reduction in convective heat transfer coefficient leading to a higher
heat transfer rate of the solid finned HX compared to the unfinned HX
as shown in Fig. 9. When the perforated fins are applied instead of the
solid fins, the Nusselt number increases because of the minor ham-
pering effect of the perforated fins on the development of the convec-
tion flows compared to the solid fins. The small reduction in the heat
transfer area due to the perforations is offset by the enhancement of the
convective heat transfer coefficient causing the higher heat transfer rate
of the perforated finned HX compared to the solid finned HX.
Also, Fig. 11 shows that increasing the HTF inlet temperature or
inlet flow rate increases the Nusselt number. At an inlet flow rate of 1 l/
min, rising the HTF inlet temperature from 55 to 65 °C provides 15.3%
and 10.5% enhancement in Nusselt number in the solid and perforated
finned HXs, respectively.
10
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
Fig. 10. Time-averaged heat fluxes of the unfinned, solid finned and perforated
finned HXs under different HTF inlet temperatures and flow rates.
Fig. 12. Comparison of total melting times of the PCM in the unfinned and solid
finned HXs under different HTF inlet temperatures and flow rates.
Table 3
Melting time reductions in the unfinned and solid finned HXs.
HTF inlet temperature HTF flow rate Melting time
(°C) (l/min) reduction (%)
Unfinned HX 55 0.5 0
55 1 24
65 0.5 51
65 1 60
Solid finned HX 55 0.5 78
55 1 82
65 0.5 88
65 1 90
11
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
Fig. 15. Comparison of cumulative energy storage with the time between the
solid and perforated finned HXs under different flow rates and a HTF inlet
temperature of 65 °C.
Fig. 13. Comparison of total melting times of the PCM in the solid and perfo-
rated finned HXs under different HTF inlet temperatures and flow rates.
Table 4
Melting time reductions in the solid and perforated finned HXs.
HX type HTF inlet HTF flow Melting time
temperature (°C) rate (l/min) reduction (%)
Fig. 16. Comparison of the classical fin effectiveness of the solid and perforated
finned HXs under different HTF inlet temperatures and flow rates.
Fig. 16 presents the classical fin effectiveness for both the solid and
perforated finned HXs at different HTF inlet temperatures and flow
rates. It is found that, for all cases, the classical fin effectiveness of the
perforated finned HX are greater than those of the solid finned HX in-
dicating that perforated fins are more effective for improving the
Fig. 14. Comparison of cumulative energy storage with time between the un- melting heat transfer rate compared to the solid fins. This is attributed
finned and solid finned HXs under different flow rates and a HTF inlet tem- to the existence of the perforations, which provide paths for the flow of
perature of 65 °C. the liquid PCM through the fins and accelerate the convective heat
transfer rate between the fins and the liquid PCM. The figure also
12
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
• At a fixed working condition, the time-averaged heat transfer rate of [9] Xu T, Chiu JN, Palm B, Sawalha S. Experimental investigation on cylindrically
macro-encapsulated latent heat storage for space heating applications. Energy
the perforated finned HX is about 10% higher than that of the solid Convers Manag 2019;182:166–77. https://doi.org/10.1016/j.enconman.2018.12.
finned HX. 056.
• The time-averaged Nusselt number of the solid finned HX is about [10] Yang X-H, Tan S-C, He Z-Z, Liu J. Finned heat pipe assisted low melting point metal
13
R. Karami and B. Kamkari Energy Conversion and Management 210 (2020) 112679
PCM heat sink against extremely high power thermal shock. Energy Convers Manag 10.1016/j.enconman.2019.03.022.
2018;160:467–76. https://doi.org/10.1016/J.ENCONMAN.2018.01.056. [30] Mahdi JM, Nsofor EC. Solidification enhancement of PCM in a triplex-tube thermal
[11] Farzanehnia A, Khatibi M, Sardarabadi M, Passandideh-Fard M. Experimental in- energy storage system with nanoparticles and fins. Appl Energy 2018;211:975–86.
vestigation of multiwall carbon nanotube/paraffin based heat sink for electronic https://doi.org/10.1016/j.apenergy.2017.11.082.
device thermal management. Energy Convers Manag 2019;179:314–25. https://doi. [31] Khan Z, Khan Z, Tabeshf K. Parametric investigations to enhance thermal perfor-
org/10.1016/J.ENCONMAN.2018.10.037. mance of paraffin through a novel geometrical configuration of shell and tube latent
[12] Qin P, Liao M, Zhang D, Liu Y, Sun J, Wang Q. Experimental and numerical study on thermal storage system. Energy Convers Manag 2016;127:355–65. https://doi.org/
a novel hybrid battery thermal management system integrated forced-air convec- 10.1016/j.enconman.2016.09.030.
tion and phase change material. Energy Convers Manag 2019;195:1371–81. [32] Liu C, Groulx D. Experimental study of the phase change heat transfer inside a
https://doi.org/10.1016/j.enconman.2019.05.084. horizontal cylindrical latent heat energy storage system. Int J Therm Sci
[13] Huang YH, Cheng WL, Zhao R. Thermal management of Li-ion battery pack with the 2014;82:100–10. https://doi.org/10.1016/j.ijthermalsci.2014.03.014.
application of flexible form-stable composite phase change materials. Energy [33] Mahdi JM, Lohrasbi S, Ganji DD, Nsofor EC. Simultaneous energy storage and re-
Convers Manag 2019;182:9–20. https://doi.org/10.1016/j.enconman.2018.12.064. covery in the triplex-tube heat exchanger with PCM, copper fins and Al2O3 nano-
[14] Bista S, Hosseini SE, Owens E, Phillips G. Performance improvement and energy particles. Energy Convers Manag 2019;180:949–61. https://doi.org/10.1016/j.
consumption reduction in refrigeration systems using phase change material (PCM). enconman.2018.11.038.
Appl Therm Eng 2018;142:723–35. https://doi.org/10.1016/j.applthermaleng. [34] Tay NHS, Bruno F, Belusko M. Comparison of pinned and finned tubes in a phase
2018.07.068. change thermal energy storage system using CFD. Appl Energy 2013;104:79–86.
[15] Soni V, Kumar A, Jain VK. Performance evaluation of nano-enhanced phase change https://doi.org/10.1016/j.apenergy.2012.10.040.
materials during discharge stage in waste heat recovery. Renew Energy [35] Rozenfeld A, Kozak Y, Rozenfeld T, Ziskind G. Experimental demonstration, mod-
2018;127:587–601. https://doi.org/10.1016/j.renene.2018.05.009. eling and analysis of a novel latent- heat thermal energy storage unit with a helical
[16] Qureshi ZA, Ali HM, Khushnood S. Recent advances on thermal conductivity en- fin. Int J Heat Mass Transf 2017;110:692–709. https://doi.org/10.1016/j.
hancement of phase change materials for energy storage system: A review. Int J ijheatmasstransfer.2017.03.020.
Heat Mass Transf 2018;127:838–56. https://doi.org/10.1016/J. [36] Ogoh W, Groulx D. Effects of the number and distribution of fins on the storage
IJHEATMASSTRANSFER.2018.08.049. characteristics of a cylindrical latent heat energy storage system : a numerical study.
[17] Gasia J, Maldonado JM, Galati F, De Simone M, Cabeza LF. Experimental evaluation Heat Mass Transf 2012;48:1825–35. https://doi.org/10.1007/s00231-012-1029-3.
of the use of fins and metal wool as heat transfer enhancement techniques in a [37] Yang X, Lu Z, Bai Q, Zhang Q, Jin L, Yan J. Thermal performance of a shell-and-tube
latent heat thermal energy storage system. Energy Convers Manag 2019;184:530–8. latent heat thermal energy storage unit: Role of annular fins. Appl Energy
https://doi.org/10.1016/j.enconman.2019.01.085. 2017;202:558–70. https://doi.org/10.1016/J.APENERGY.2017.05.007.
[18] Wang P, Yao H, Lan Z, Peng Z, Huang Y, Ding Y. Numerical investigation of PCM [38] Parsazadeh M, Duan X. Numerical study on the effects of fins and nanoparticles in a
melting process in sleeve tube with internal fins. Energy Convers Manag shell and tube phase change thermal energy storage unit. Appl Energy
2016;110:428–35. https://doi.org/10.1016/j.enconman.2015.12.042. 2018;216:142–56. https://doi.org/10.1016/j.apenergy.2018.02.052.
[19] Tauseef-ur-Rehman HM, Ali MM, Janjua U, Sajjad W-M. Yan, A critical review on [39] Shahsavar A, Goodarzi A, Mohammed HI, Shirneshan A, Talebizadehsardari P.
heat transfer augmentation of phase change materials embedded with porous ma- Thermal performance evaluation of non-uniform fin array in a finned double-pipe
terials/foams. Int J Heat Mass Transf 2019;135:649–73. https://doi.org/10.1016/J. latent heat storage system. Energy 2020. https://doi.org/10.1016/j.energy.2019.
IJHEATMASSTRANSFER.2019.02.001. 116800.
[20] Sardari PT, Grant D, Giddings D, Walker GS, Gillott M. Composite metal foam/PCM [40] Abu-Hijleh BAK. Natural convection heat transfer from a cylinder with high con-
energy store design for dwelling space air heating. Energy Convers Manag ductivity permeable fins. J Heat Transfer 2003;125:282–8. https://doi.org/10.
2019;201. https://doi.org/10.1016/j.enconman.2019.112151. 1115/1.1532013.
[21] Tiari S, Mahdavi M, Qiu S. Experimental study of a latent heat thermal energy [41] Suryawanshi SD, Sane NK. Natural Convection Heat Transfer From Horizontal
storage system assisted by a heat pipe network. Energy Convers Manag Rectangular Inverted Notched Fin Arrays. J Heat Transfer 2009;131:082501https://
2017;153:362–73. https://doi.org/10.1016/j.enconman.2017.10.019. doi.org/10.1115/1.3109993.
[22] Leong KY, Abdul Rahman MR, Gurunathan BA. Nano-enhanced phase change ma- [42] AlEssa AH, Al-Widyan MI. Enhancement of natural convection heat transfer from a
terials: A review of thermo-physical properties, applications and challenges. J fin by triangular perforation of bases parallel and toward its tip. Appl Math Mech
Energy Storage 2019;21:18–31. https://doi.org/10.1016/j.est.2018.11.008. (English Ed) 2008;29:1033–44. https://doi.org/10.1007/s10483-008-0807-x.
[23] Abdulateef AM, Mat S, Abdulateef J, Sopian K, Al-Abidi AA. Geometric and design [43] Awasarmol UV, Pise AT. An experimental investigation of natural convection heat
parameters of fins employed for enhancing thermal energy storage systems: A re- transfer enhancement from perforated rectangular fins array at different inclina-
view. Renew Sustain Energy Rev 2018;82:1620–35. https://doi.org/10.1016/j.rser. tions. Exp Therm Fluid Sci 2015;68:145–54. https://doi.org/10.1016/j.
2017.07.009. expthermflusci.2015.04.008.
[24] Kamkari B, Shokouhmand H, Bruno F. Experimental investigation of the effect of [44] Desgrosseilliers L, Whitman CA, Groulx D, White MA. Dodecanoic acid as a pro-
inclination angle on convection-driven melting of phase change material in a rec- mising phase-change material for thermal energy storage. Appl Therm Eng
tangular enclosure. Int J Heat Mass Transf 2014;72:186–200. https://doi.org/10. 2013;53:37–41. https://doi.org/10.1016/j.applthermaleng.2012.12.031.
1016/j.ijheatmasstransfer.2014.01.014. [45] Shokouhmand H, Kamkari B. Experimental investigation on melting heat transfer
[25] Kamkari B, Amlashi HJ. Numerical simulation and experimental verification of characteristics of lauric acid in a rectangular thermal storage unit. Exp Therm Fluid
constrained melting of phase change material in inclined rectangular enclosures. Int Sci 2013;50:201–12. https://doi.org/10.1016/j.expthermflusci.2013.06.010.
Commun Heat Mass Transf 2017;88:211–9. https://doi.org/10.1016/j. [46] Tian LL, Liu X, Chen S, Shen ZG. Effect of fin material on PCM melting in a rec-
icheatmasstransfer.2017.07.023. tangular enclosure. Appl Therm Eng 2020;167. https://doi.org/10.1016/j.
[26] Rabienataj Darzi AA, Jourabian M, Farhadi M. Melting and solidification of PCM applthermaleng.2019.114764.
enhanced by radial conductive fins and nanoparticles in cylindrical annulus. Energy [47] Robak CW, Bergman TL, Faghri A. Enhancement of latent heat energy storage using
Convers Manag 118 2016:253–63. https://doi.org/10.1016/j.enconman.2016.04. embedded heat pipes. Int J Heat Mass Transf 2011;54:3476–84. https://doi.org/10.
016. 1016/j.ijheatmasstransfer.2011.03.038.
[27] Kamkari B, Shokouhmand H. Experimental investigation of phase change material [48] Kline SJ, McClintock FA. Describing uncertainties in single sample experiments.
melting in rectangular enclosures with horizontal partial fins. Int J Heat Mass Mech Eng 1953;75:3–8.
Transf 2014;78:839–51. https://doi.org/10.1016/j.ijheatmasstransfer.2014.07. [49] Kamkari B, Groulx D. Experimental investigation of melting behaviour of phase
056. change material in finned rectangular enclosures under different inclination angles.
[28] Yuan Y, Cao X, Xiang B, Du Y. Effect of installation angle of fins on melting char- Exp Therm Fluid Sci 2018;97:94–108. https://doi.org/10.1016/J.
acteristics of annular unit for latent heat thermal energy storage. Sol Energy EXPTHERMFLUSCI.2018.04.007.
2016;136:365–78. https://doi.org/10.1016/j.solener.2016.07.014. [50] Rathod MK, Banerjee J. Thermal performance enhancement of shell and tube Latent
[29] Sodhi GS, Jaiswal AK, Vigneshwaran K, Muthukumar P. Investigation of charging Heat Storage Unit using longitudinal fins. Appl Therm Eng 2015;75:1084–92.
and discharging characteristics of a horizontal conical shell and tube latent thermal https://doi.org/10.1016/j.applthermaleng.2014.10.074.
energy storage device. Energy Convers Manag 2019;188:381–97. https://doi.org/
14