Evaluation of Melting Mechanism and Natural Convection Effect in a Triplex Tube Heat Storage System with a Novel Fin Arrangement

: In this research, a numerical analysis is accomplished aiming to investigate the effects of adding a new design ﬁns arrangement to a vertical triplex tube latent heat storage system during the melting mechanism and evaluate the natural convection effect using Ansys Fluent software. In the triplex tube, phase change material (PCM) is included in the middle tube, while the heat transfer ﬂuid (HTF) ﬂows through the interior and exterior pipes. The proposed ﬁns are triangular ﬁns attached to the pipe inside the PCM domain in two different ways: (1) the base of the triangular ﬁns is connected to the pipe, (2) the tip of the triangular ﬁns is attached to the pipe and the base part is directed to the PCM domain. The height of the ﬁns is calculated to have a volume equal to that of the uniform rectangular ﬁns. Three different cases are considered as the ﬁnal evaluation toward the best case as follows: (1) the uniform ﬁn case (case 3), (2) the reverse triangular ﬁn case with a constant base (case 12), (3) the reverse triangular ﬁn case with a constant height (case 13). The numerical results show that the total melting times for cases 3 and 12 increase by 4.0 and 10.1%, respectively, compared with that for case 13. Since the PCM at the bottom of the heat storage unit melts slower due to the natural convection effect, a ﬂat ﬁn is added to the bottom of the heat storage unit for the best case compared with the uniform ﬁn cases. Furthermore, the heat storage rates for cases 3 and 12 are reduced by 4.5 and 8.5%, respectively, compared with that for case 13, which is selected as the best case due to having the lowest melting time (1978s) and the highest heat storage rate (81.5 W). The general outcome of this research reveals that utilizing the tringle ﬁns enhances the thermal performance and the phase change rate.


Introduction
The overconsumption of fossil fuels has resulted in increasing the level of greenhouse gas emissions and the depletion of the ozone layer [1]. Therefore, renewable sources of energy have been considered as an appropriate alternative to fossil fuels [2,3]. Solar energy, as one of the most prospective sources of energy, benefits from many advantages over fossil fuels, including its safety and friendliness to the environment, which is freely available in abundance [4]. However, there is a mismatch between energy supply and demand exists in solar energy systems that can be well managed using thermal energy storage (TES) systems [5]. Moreover, energy storage can be used in passive thermal control applications, which is highly important in a net-zero emission environment in the near future [6][7][8]. TES systems have been widely used in various applications, such as buildings, solar heating storage, refrigeration, solar power plants, etc., [9][10][11][12].
Energy storage methods are grouped into three categories, including sensible thermal energy, latent heat, and thermochemical energy [13]. Materials with long-term stability and high specific heat capacity can improve the performance of sensible heat storage systems [14]. In latent heat energy storage systems, energy is released or absorbed when the phase change process occurs. Applying the phase change materials (PCMs) in the latent heat energy storage systems is known as one of the most efficient methods to store thermal energy. The amount of energy stored by PCMs during the phase change process depends on the mass and latent heat of PCMs [15]. However, the conductivity of the PCMs applied as the storage media in the TES systems is considerably low (k ≤ 0.2 W/m K). Therefore, material failure and system overheating may occur due to a high-temperature difference within the PCM and incomplete melting and solidification processes. Improving the thermal performance of the TES systems is of great importance due to their widespread utilization in solar systems, drying technology, waste heat recovery, electronic cooling, and air conditioning [16]. The thermal performance of PCMs can be improved efficiently via various techniques proposed by researchers, including the implementation of fins with high thermal performance [17][18][19][20], nano-particles [21][22][23][24], multiple-PCMs [25,26], nanoencapsulating PCMs [27,28], high conducting particles [29,30], use of porous media [31][32][33][34], imposing magnetic fields [27,35], and improving the convection heat transfer [28,36]. Further, improving the performances of the energy storage and heat exchanger have been investigated in different works [37][38][39] Many researchers investigated the feasibility of using various types of fins to improve the thermal performance of TES systems both numerically and experimentally [40][41][42][43]. Applying fins in the TES systems to enhance the thermal performance of the system benefits from some advantages, i.e., low construction costs, high efficiency, and easy fabrication [44]. An experimental energy storage system using a horizontal triplex tube heat exchanger (TTHX) in the presence of longitudinal/triangular fins was developed, tested, and analyzed by Abdulateef et al. [45] to improve the melting process of PCMs. They reported that applying internal, internal-external and external triangular fins resulted in notable improvement of 11, 12, and 15%, respectively, compared with those of the cases with longitudinal fins. Guo et al. [46] assessed the melting process of a PCM in a novel LHTES system with an angled fin design. Senthil et al. [47] studied the effects of adding fins with variable lengths on the melting behavior of a PCM in a vertical solar system. Mahdi et al. [48] enhanced the melting rate of PCM in a triplex tube heat exchanger (TTHX) by applying a novel fin configuration, compared with that of using nanoparticles. They reported that the melting process was accelerated via using long fins in the regions where the conduction heat transfer was dominant. Furthermore, better thermal performance was achieved via applying fewer and relatively shorter fins at the upper half of the system. Abdulateef et al. [49] studied the feasibility of using a fin-nanoparticle combination to accelerate the melting and solidification of a PCM in a TTHX. Fin parameters were optimized to improve the heat transfer rate inside the PCM. They concluded that the lowest solidification and melting times were achieved for the system with an external triangular fins-nanoparticles combination, considering the fins' number of 8, fins' length of 141 mm, and fins' aspect ratio of 18% [50].
Sun et al. [51] investigated the impacts of using circular fins with a staggered arrangement to intensify the thermal response of a PCM inside a vertical triple-tube heat exchanger with a counter-current configuration. They reported that the PCM melting rate increases significantly due to the staggered arrangement of fins. Hassan et al. [5] studied the advancements in thermal properties of PCMs and thoroughly discussed the terms of enhancement in melting and solidification rates. Sivalakshmi et al. [52] studied the influences of helical fins addition on the thermal performance of a double pipe heat exchanger. To evaluate the performance of the system, various parameters, including heat transfer coefficient, average heat transfer rate, and heat exchanger effectiveness in the plain inner pipe, were assessed and compared with those of the case with helical fins installed over the inner pipe. It was concluded that the installation of fins over the inner pipe of the double-pipe heat exchanger increases the heat transfer coefficient significantly. Abdulateef et al. [53] proposed a fins-nanoparticle combination as an augmentation technique to improve the melting process of a PCM in a TTHX-based TES system. They applied paraffin (RT82) and Alumina (Al2O3) nanoparticles to study the thermal performance of TTHX numerically and experimentally. They deduced that applying the fins-nanoparticle combination results in an improvement in the melting process of the PCM compared with that of the case without nanoparticles. Safari et al. [54] investigated the feasibility of using bifurcated and straight fins inside a shell-and-tube heat exchanger to enhance the thermal performance of a PCM. They also analyzed the effects of fin's length and thickness on the thermal performance of the heat exchanger while the mass of fins was assumed to be constant for all cases. It was concluded that the thermal performance of the case with a bifurcated fin was superior to that of the case with a straight fin. Furthermore, fin length played an outstanding role to determine the optimum fin arrangement. Masoumpour-Samakoush et al. [55] presented a novel combination of triangular and rectangular fins to study the melting process of a PCM in a cavity numerically. The impacts of various parameters, including the schematic of internal fins, hot wall temperature, height of triangular fins, and various kinds of PMCs on the melting phenomenon, were evaluated. The numerical results showed that increasing the number of fins results in a shorter melting time and consequently a higher energy charge rate. Ghalambaz et al. [56] utilized variable-length fins to improve the heat transfer during the melting. The length of the fins could be reduced in an ascending or descending arrangement. They reported that the fin's aspect ratio is one of the most important parameters that can accelerate the melting process. Najim et al. [57] investigated the impacts of circular fins on accelerating the thermal response of a PCM in a vertical triple-tube heat exchanger. The impacts of improving the fin geometrical features, as well as the arrangement of the fins in various spatial zones of heat exchangers were determined. They also conducted a parametric analysis to find the best fin's structure, considering the transient behavior of the melting process and rates of heat storage. Punniakodi and Senthil [58] comprehensively reviewed the container geometry and orientation of a PCM in solar thermal systems. It was concluded that variations in fins' arrangement, as well as the direction of the heat transfer flow results in a faster melting process.
Since the tringle fins in the triple-tube PCM were not investigated in this field previously, what are the impacts of the different orientations and sizes of the tringle fins in the thermal energy storage system? This paper performed a numerical analysis of the performance of a vertical triple-tube latent heat storage system to investigate the impacts of adding triangular fins for heat transfer enhancement, compared with that using rectangular fins. The motivation of this work is to receive a better understanding of how different shapes of the fins added affect the melting performance of the system. The triangular fins are added to the inner and the middle tubes in both direct and upturned connections. The main contribution of this work is presenting a new architecture of thermal energy storage systems at the design stage by combining the effects of the shape of the fins and orientations of the fins. Firstly, two scenarios are considered for the added fins: (1) the base of the triangular fins is connected to the pipe (2) the tip of the triangular fins is attached to the pipes while the base is directed to the PCM domain. Then, various dimensions for the base of the triangular fins are studied for both the above-mentioned orientations of the fins. The height of the fins is calculated to have a surface area equal to that of the uniform rectangular fins. After determining the best dimension via comparing the liquid fraction and temperature counters, time of melting process, and heat storage values of various cases, the effect of adding a flat plate to the bottom of the heat storage unit is evaluated to improve the thermal performance of the proposed cases toward a higher heat transfer rate from the HTF to the PCM and thus lower melting time.

System Description
In this study, a vertical finned triple-tube PCM heat storage is examined considering five triangular fins connected to the interior and exterior sides of the pipe inside the middle tube. Thus, the total number of fins is considered to be ten. The triangular fins are connected in two ways, i.e., (1) the base of the triangular fin is connected to the pipe, (2) the base of the fin is located inside the PCM domain, and the pointed part of the fins is connected to the pipe. Note that the PCM is filled in the middle tube, while the heat transfer fluid (which is water) flows through the interior and exterior pipes. Furthermore, the configuration of the HTF in the inner and outer tubes is supposed to be counter-current according to the previous research in the literature [59]. In other words, the HTF flows through the outer tube in the gravity direction, while the HTF in the inner pipe is in the opposite gravity direction. The Reynolds number of the HTF flow is 2000 in both inner and outer tubes (velocity of 0.0269 m/s), which results in the flow rate of 33.8 and 169 cm 3 /s, respectively.
The schematics of the proposed systems, including the no-fin case, uniform-fin case, triangular finned (TF) case, reverse triangular finned (RTF) case, and the RTF case with the added fin are shown in Figure 1. It should be noted that the flow is considered axisymmetric because of the characteristics of the problem, as well as the lack of circumferential flow variations in the studied system. The triangular fins' faces are perpendicular to the axis of the PCM tube. The boundary conditions, as well as the dimensions of the system, are shown in Figure 1a for the no-fin case. Note that the dimension for the PCM domain is selected based on various studies on triplex tube latent heat storage systems in the literature [22,57]. and temperature counters, time of melting process, and heat storage values of various cases, the effect of adding a flat plate to the bottom of the heat storage unit is evaluated to improve the thermal performance of the proposed cases toward a higher heat transfer rate from the HTF to the PCM and thus lower melting time.

System Description
In this study, a vertical finned triple-tube PCM heat storage is examined considering five triangular fins connected to the interior and exterior sides of the pipe inside the middle tube. Thus, the total number of fins is considered to be ten. The triangular fins are connected in two ways, i.e., (1) the base of the triangular fin is connected to the pipe, (2) the base of the fin is located inside the PCM domain, and the pointed part of the fins is connected to the pipe. Note that the PCM is filled in the middle tube, while the heat transfer fluid (which is water) flows through the interior and exterior pipes. Furthermore, the configuration of the HTF in the inner and outer tubes is supposed to be counter-current according to the previous research in the literature [59]. In other words, the HTF flows through the outer tube in the gravity direction, while the HTF in the inner pipe is in the opposite gravity direction. The Reynolds number of the HTF flow is 2000 in both inner and outer tubes (velocity of 0.0269 m/s), which results in the flow rate of 33.8 and 169 cm 3 /s, respectively.
The schematics of the proposed systems, including the no-fin case, uniform-fin case, triangular finned (TF) case, reverse triangular finned (RTF) case, and the RTF case with the added fin are shown in Figure 1. It should be noted that the flow is considered axisymmetric because of the characteristics of the problem, as well as the lack of circumferential flow variations in the studied system. The triangular fins' faces are perpendicular to the axis of the PCM tube. The boundary conditions, as well as the dimensions of the system, are shown in Figure 1a for the no-fin case. Note that the dimension for the PCM domain is selected based on various studies on triplex tube latent heat storage systems in the literature [22,57].   Different values of 2.5, 3, 3.5, and 4 mm are considered for the base of the triangular fin in comparison with the uniform-fin case, where the dimension of the fins is 2 × 5 mm, as well as the no-fin case. It should be noted that the minimum value of the base of the triangular fins is considered 2.5 mm, since, for the lower values (≤2), the height of the triangular fins is obtained higher than 10 mm, which is not possible according to the diameter of the middle tube. The height of the triangular fins is then calculated, considering the equal area for the fins in comparison with the uniform-fin case. Therefore, the PCM mass is considered similar in all the studied cases. The dimensions of the fins with the description and name of each proposed case are presented in Table 1. In the previous paper by the authors [57], it was shown that, due to the natural convection effect in a vertical heat storage system during the phase change process, the melting process occurs in the upper region of the heat storage system faster compared to the bottom section of the heat storage unit. Thus, it was recommended to add a rectangular fin to the bottom of the heat exchanger to increase the heat transfer rate. Thus, in this paper, for the best case, a fin is added at the bottom, and its thermal performance is compared with that of the uniform-fin case with added fin, as shown in Table 1. The dimensions of the added fin are 1 × 20 mm, which is considered similar to all the added fin cases. For the best case with the added fin, two configurations are evaluated, i.e., (1) the base of the triangular fins is kept constant, and the height is then calculated; (2) the height of the triangular fins is kept constant, and the base is then calculated.
Paraffin RT-35 is used in the present study. The melting temperature of this PCM is suitable for HVAC applications in buildings. The properties of RT-35 are presented in Table 2.

Mathematical Modeling
To simulate the phase change phenomenon in the PCM domain, the enthalpy-porosity approach is selected to be used in this study. Density variations of the PCM result in the natural convection effect in the domain, which was determined using the Boussinesq approximation. The fluid flow was assumed to be Newtonian, transient, and incompressible. Furthermore, the flow regime was laminar. The effects of viscous dissipation are neglected in the study, while the boundary condition at the walls is assumed to be a no-slip boundary condition. The boundary condition at the outer walls of the heat exchanger is adiabatic, while high-quality insulation is considered for the walls of the PCM domain at the top and bottom of the heat exchanger. The following assumptions are considered in this study as follows [48]:

3.
Consider variable density while the other physical properties of PCM are constant.

4.
Apply Boussinesq approximation for density variation in the liquid PCM.
Based on the assumption made, the governing equations of continuity, momentum, and energy can be developed as follows [60]: It should be noted that the term ( → S ) in Equation (2) is defined to consider the effect of phase transition interpreted as the velocity damping of the Darcy law [61] as follows: It should be mentioned that the value of Am called the mushy zone was assumed to be 105, based on the previous studies available in the literature [38,39]. In order to analyze the phase change phenomenon accurately, the PCM liquid fraction (λ) is defined as follows [40]: In order to estimate the value of (S L ) in Equation (3), Equation (6) is defined as follows: Regarding the melting process, Equation (7) can be applied to estimate the rate of energy stored: where E ini and E end stand for the total energy amounts of the PCM at the initial and final times of the phase-transition course, respectively, which are the summation of sensible heat (MC p dT), and latent heat (ML f ) of the PCM. t m is the melting time.

Numerical Analysis
To evaluate the thermofluidic characteristics of PCM during the melting phenomenon, the SIMPLE methodology is used combined with a Green-Gauss cell-based meshing technique in the ANSYS-FLUENT environment. The QUICK differencing technique is applied to solve the momentum and energy equations, while the pressure corrected equations are used numerically using the PRESTO technique. After achieving a universal pre-selection, the under-relaxation values for liquid fraction, velocity components, energy equation, and pressure correction equation are assumed to be 0.5, 0.3, 1.0, and 0.3, respectively. To stop the iterative solutions of the governing equations at each time step, the convergence criteria for the continuity, momentum, and energy equations are assumed to be 10 −4 , 10 −4 , and 10 −6 , respectively.

Spatial and Temporal Discretization
To apply the grid and time step size independency examinations performed for the uniform-fin case, the corresponding tests are accomplished applying the time step size of 0.2 with different cell numbers of 28,500, 43,000, and 81,620. For various cell numbers, variations of liquid fraction versus time are illustrated in Figure 2. As can be seen, values of liquid fraction for the grid sizes of 43,000 and 81,620 are almost equal (the maximum difference at the identical time is less than 0.2%). Therefore, the grid size value is set to 43,000 for further analysis. To analyze the impacts of the time step size on the melting time, the melting time amounts for different time step sizes are calculated and shown in Table 3. According to Table 3, the melting times for various time step size values are almost identical. The difference between the time step sizes of 0.1 s and 0.2 s is around 0.13%. Thus, the value of the time step size is considered 0.2 s for the next studies.

Model Validation
To verify the accuracy of the numerical scheme applied in the present study, the numerical results of the present research are validated against those of Mat et al. [62]. The same boundary conditions and performance parameters of the reference system are applied in the present model. Mat et al. [62] conducted a set of numerical and experimental investigations to explore the melting performance of a PCM (RT58) in a finned doubletube TES system. The objective of their research was to evaluate the impacts of adding staggered fins to the exterior and interior walls of the PCM shell. The data for the average PCM temperature and liquid fraction from the present study are compared to the numer-

Model Validation
To verify the accuracy of the numerical scheme applied in the present study, the numerical results of the present research are validated against those of Mat et al. [62]. The same boundary conditions and performance parameters of the reference system are applied in the present model. Mat et al. [62] conducted a set of numerical and experimental investigations to explore the melting performance of a PCM (RT58) in a finned double-tube TES system. The objective of their research was to evaluate the impacts of adding staggered fins to the exterior and interior walls of the PCM shell. The data for the average PCM temperature and liquid fraction from the present study are compared to the numerical and experimental findings reported by Mat et al. [62] as shown in Figure 3. The statistical validation gave a maximum percentage error of 1.4%. Thus, the present model can be adopted for exploring the thermal characteristics of the PCM-based triple-tube system in the presence of triangular fins. It should be noted that the experimental uncertainty in the study of Mat et al. [62] is ±0.15 • C (around 0.2%) for the temperature measurement.

Results and Discussion
In this section, the effects of different orientations and dimensions of the triangular fins on the charging process are discussed. Then, the performance of the determined optimized system is compared with those of the uniform-fin and the no-fin cases. The final analysis investigates the effects of adding a flat fin to the bottom of the heat exchanger on the melting rate, temperature distribution, and heat storage rate. It should be mentioned that the major objective of this work is to detect the capability of the triangular fins to provide better storage operation in the thermal unit. Besides, as stated in the introduction, numerous works on double and triple-tube PCM-based heat storage units in the presence of circular or longitudinal fins have been extensively examined to enhance the efficiency of heat storage units. This study offers an additional step ahead, compared with other general works in the literature to determine an efficient finned triple-tube heat storage unit in the presence of triangular fins.

Effects of Using Triangular Fins on the Charging Phenomenon
Using fins in the heat storage unit enhances the performance of the unit due to increasing the heat transfer surface area. Increasing the average thermal conductivity of the system due to utilizing a metal fin with high conductivity also plays a major role in improving the performance. Furthermore, the fins are also capable of transferring heat to the depth of the PCM from both sides, including the inner and outer sides of the PCM domain. The fins have a dominant role in free convection behavior inside the domain as they form a barrier against the molten PCM movement. Five pairs of triangular fins are attached to the shell part of the system on the inner and the outer sides. Various cases studied in this research are shown in Table 1. The no-fin case (case 1) is considered as the reference case. The uniform-fin case with the dimensions of 2 mm × 5 mm (base × height) is considered as case 2, while case 3 stands for the uniform-fin case with an added fin to the bottom of the system. Case 4 includes five opposite triangular fins with dimensions of 2.5 mm and 8 mm (base and height) in a horizontal direction (TF-base 2.5 case). In cases 5, 6, and 7, the dimensions of the fins are changed by increasing the base value, while the fin height is

Results and Discussion
In this section, the effects of different orientations and dimensions of the triangular fins on the charging process are discussed. Then, the performance of the determined optimized system is compared with those of the uniform-fin and the no-fin cases. The final analysis investigates the effects of adding a flat fin to the bottom of the heat exchanger on the melting rate, temperature distribution, and heat storage rate. It should be mentioned that the major objective of this work is to detect the capability of the triangular fins to provide better storage operation in the thermal unit. Besides, as stated in the introduction, numerous works on double and triple-tube PCM-based heat storage units in the presence of circular or longitudinal fins have been extensively examined to enhance the efficiency of heat storage units. This study offers an additional step ahead, compared with other general works in the literature to determine an efficient finned triple-tube heat storage unit in the presence of triangular fins.

Effects of Using Triangular Fins on the Charging Phenomenon
Using fins in the heat storage unit enhances the performance of the unit due to increasing the heat transfer surface area. Increasing the average thermal conductivity of the system due to utilizing a metal fin with high conductivity also plays a major role in improving the performance. Furthermore, the fins are also capable of transferring heat to the depth of the PCM from both sides, including the inner and outer sides of the PCM domain. The fins have a dominant role in free convection behavior inside the domain as they form a barrier against the molten PCM movement. Five pairs of triangular fins are attached to the shell part of the system on the inner and the outer sides. Various cases studied in this research are shown in Table 1. The no-fin case (case 1) is considered as the reference case. The uniform-fin case with the dimensions of 2 mm × 5 mm (base × height) is considered as case 2, while case 3 stands for the uniform-fin case with an added fin to the bottom of the system. Case 4 includes five opposite triangular fins with dimensions of 2.5 mm and 8 mm (base and height) in a horizontal direction (TF-base 2.5 case). In cases 5, 6, and 7, the dimensions of the fins are changed by increasing the base value, while the fin height is reduced to keep the cross-section area of the fin constant (10 mm 2 for each fin) in all the cases. Figure 4 shows the liquid fraction development for cases 4 to 7 at various times, including 300, 600, 1200, 1800, and 2400 s. For all cases, the PCM was melted at the attached region of the wall, firstly. During the first 300 s, the adjacent PCM to the inner and outer HTF channel walls is melted. The molten area expands gradually during the melting. All the cases almost show the same behavior, and the solid part is confined between the two adjacent fins. Within 1200 s, the solid parts between the fins are separated, the solid portions remain on the opposite fins, and the area between the two opposite fins is melted. Over time, more PCM is liquefied, forming a thicker liquid layer, which is surrounded by the solid part. The molten PCM is concentrated at the top of the heat exchanger because of the natural convection effect. At the bottom of the system, the PCM is still solid. In t = 1800 s, the percentage of the PCM melted for case 4 (the case with the highest fin height (8 mm)) and case 7 (the case with the lowest fin height (5 mm)) are 87 and 83%, respectively. It can be explained that in case 4, which has the largest surface area, more heat is transferred to the PCM compared to the case with the shortest fins' height. Because of the natural convection effect, circulation of the liquid PCM is formed in the clockwise direction. On the one hand, the triangular fins act as a barrier against this movement, which limits the free convection effect. On the other hand, the slope orientation of the triangular fins reduces this barrier and helps the molten PCM to slip on the side of the fin. The lesser the slope (case with the highest fin height), the better the barrier to keep the PCM inside the zone between the two pairs of adjacent fins. The higher slope of the fins side (the case with the shortest fin height) allows better circulation of the PCM through the zones, which helps the liquified PCM to move faster to the upper regions, and the solid PCM gathers at the bottom of the system. Figure 5 illustrates the temperature distribution within 2400 s for cases 4 to 7. Because of the short length of the system, the HTF channels almost remain at a fixed temperature (shown in red color) during the melting process. For all the cases, the temperature increases in the adjacent region to the HTF channel walls and around the triangular fins. The temperature of the fins does not reach the thermal equilibrium, as the fins receive heat from the base side of the triangular fin and releases heat to the PCM from the other side. Therefore, the fins temperature records higher and lower values compared to the PCM and HTF, respectively. The closest area to the center of the PCM domain has the lowest temperature as this area is far from the wall and the natural convection effect is not strong enough in this area. The temperature increases gradually within 2400 s and the average temperatures reach 45.3, 45.9, 47.4, and 46.5 • C, for cases 4, 5, 6, and 7, respectively. Due to remaining the solid PCM at the bottom of the system, the bottom temperature of the heat exchanger is the lowest compared with that of the rest of the system. It should be noted that all the cases have almost the same behavior and almost the same average temperature. Figure 6a shows the development of the liquid fraction for cases 4 to 7 during the entire melting process. The charging process occurs faster during the first 300 s due to the conduction heat transfer with the solid PCM, and then the melting process is reduced slightly due to the increasing heat transfer rate by convection compared with that of conduction. Between t = 2100 and 2200 s, the curves become a plateau and more than 96% of the PCM of all the cases is melted. Case 4 has the fastest phase change process due to having the largest surface area of the fin, while case 7 has the slowest phase change phenomenon. The average temperature noticeably increases during the first 200 s because of the conduction heat transfer based on Figure 6b. Additionally, the average temperature of case 4 TF base 2.5 is the highest due to having the largest surface area. The temperature of the PCM reaches the thermal equilibrium within 3000 s. limits the free convection effect. On the other hand, the slope orientation of the triangular fins reduces this barrier and helps the molten PCM to slip on the side of the fin. The lesser the slope (case with the highest fin height), the better the barrier to keep the PCM inside the zone between the two pairs of adjacent fins. The higher slope of the fins side (the case with the shortest fin height) allows better circulation of the PCM through the zones, which helps the liquified PCM to move faster to the upper regions, and the solid PCM gathers at the bottom of the system.     Figure 5 illustrates the temperature distribution within 2400 s for cases 4 to 7. Because of the short length of the system, the HTF channels almost remain at a fixed temperature (shown in red color) during the melting process. For all the cases, the temperature increases in the adjacent region to the HTF channel walls and around the triangular fins.  Figure 6a shows the development of the liquid fraction for cases 4 to 7 during the entire melting process. The charging process occurs faster during the first 300 s due to the conduction heat transfer with the solid PCM, and then the melting process is reduced slightly due to the increasing heat transfer rate by convection compared with that of conduction. Between t = 2100 and 2200 s, the curves become a plateau and more than 96% of the PCM of all the cases is melted. Case 4 has the fastest phase change process due to having the largest surface area of the fin, while case 7 has the slowest phase change phenomenon. The average temperature noticeably increases during the first 200 s because of the conduction heat transfer based on Figure 6b. Additionally, the average temperature of case 4 TF base 2.5 is the highest due to having the largest surface area. The temperature of the PCM reaches the thermal equilibrium within 3000 s.

Effects of Using Reverse Triangular Fins on the Charging Process
In this part of the study, the impacts of triangular fin orientations on the charging process are investigated. In other words, the directions of the triangular fins are reversed and the tip part of the triangular fins is attached to the inner and the outer walls of the HTF channels, as shown in Figure 7. As mentioned before, the longer fins provide a higher surface area compared with the shorter fin. A larger surface area accelerates the melting process due to transferring more heat to the PCM. Although the point connection limits the amount of the heat transfer rate from the channel of the fin and then to the PCM; however, the heat is still transferred enough to the fin (the fins in Figures 7 and 8 are

Effects of Using Reverse Triangular Fins on the Charging Process
In this part of the study, the impacts of triangular fin orientations on the charging process are investigated. In other words, the directions of the triangular fins are reversed and the tip part of the triangular fins is attached to the inner and the outer walls of the HTF channels, as shown in Figure 7. As mentioned before, the longer fins provide a higher surface area compared with the shorter fin. A larger surface area accelerates the melting process due to transferring more heat to the PCM. Although the point connection limits the amount of the heat transfer rate from the channel of the fin and then to the PCM; however, the heat is still transferred enough to the fin (the fins in Figures 7 and 8 are presented in white color to be observed clearly). The fins form a strong barrier that limits the circulation of the molten PCM. This behavior helps the molten PCM to circulate inside each zone rather than the whole system. Within 600 s, the solid PCM is still connected through the zones. The solid parts finally separate at t = 1200 s and they are arranged on the pairs of opposite fins. The solid parts shrink during melting and, through time, most of the PCM is melted except for the region at the bottom part of the system. Case 8 shows a higher performance, compared with the other cases, as the entire zone (98.8%) is melted except for the region at the bottom of the unit. It should be noted that small portions of the solid part remain in different places of the domain for the other cases. each zone rather than the whole system. Within 600 s, the solid PCM is still connected through the zones. The solid parts finally separate at t = 1200 s and they are arranged on the pairs of opposite fins. The solid parts shrink during melting and, through time, most of the PCM is melted except for the region at the bottom part of the system. Case 8 shows a higher performance, compared with the other cases, as the entire zone (98.8%) is melted except for the region at the bottom of the unit. It should be noted that small portions of the solid part remain in different places of the domain for the other cases.   Figure 8 illustrates the temperature contours of cases 8, 9, 10, and 11. As the heat is transferred to the PCM, the PCM domain is warmed and the temperature increases gradually. The longer fins (higher surface area) transfer more heat to the PCM resulting in faster melting and higher temperature. The temperature rises at the top of the unit initially  Figure 9 shows the streamlines of cases 4 and 8 at various times. Regarding case 4, the liquid PCM rarely circulates within the zones, especially at the early stage of the melting process (300-600 s). Due to the slope direction of the fins toward the PCM domain, more movement is achieved among the zones. Because of the slope direction of the  Figure 8 illustrates the temperature contours of cases 8, 9, 10, and 11. As the heat is transferred to the PCM, the PCM domain is warmed and the temperature increases gradually. The longer fins (higher surface area) transfer more heat to the PCM resulting in faster melting and higher temperature. The temperature rises at the top of the unit initially due to the convection heat transfer. Furthermore, an increase in the temperature at the top of each zone is also observed. Regarding the shortest fin, the average temperature reaches 41 • C within 2400 s, while the corresponding value of the longest fin is 46 • C. Figure 9 shows the streamlines of cases 4 and 8 at various times. Regarding case 4, the liquid PCM rarely circulates within the zones, especially at the early stage of the melting process (300-600 s). Due to the slope direction of the fins toward the PCM domain, more movement is achieved among the zones. Because of the slope direction of the reversed fins toward the HTF channel (case 8), the liquid PCM is transferred hardly through the zones between the two pairs of the opposite fins. In other words, in the reverse triangular fins, a strong barrier is formed against transferring the molten PCM through the zone, causing all the PCM to just circulate in their own zones. The influence of the reverse triangular fin on the melting process is investigated in Figure 10a for cases 8 to 11. Different dimensions of the fins are considered in this part. The longest fin provides a faster charging process due to having the largest surface area, resulting in transferring more heat to the PCM domain. The effect of the fins height on the melting procedure appears at a constant rate until 2400 s (the PCM liquid fraction is 99% for case 8 and 96% for case 11). Then, the melting rate reduces until the PCM is melted totally. It is worth mentioning that the area of the different fins is constant for various dimensions, which guarantees the value of the PCM utilized in various cases to be the same. Figure 10b shows the average temperature development of the PCM for the reverse triangular fins. The temperature increases rapidly due to the conduction heat transfer and then reduces. At 1800 s, the temperature of the PCM for the longest fin is 40 • C, which is higher than the PCM temperature for cases 9, 10, 11 by 3%, 5%, and 7%, respectively. Within 3000 s, the PCM for all the cases reaches almost the same temperature.  The influence of the reverse triangular fin on the melting process is investigated in Figure 10a for cases 8 to 11. Different dimensions of the fins are considered in this part. The longest fin provides a faster charging process due to having the largest surface area, resulting in transferring more heat to the PCM domain. The effect of the fins height on the melting procedure appears at a constant rate until 2400 s (the PCM liquid fraction is 99% for case 8 and 96% for case 11). Then, the melting rate reduces until the PCM is melted totally. It is worth mentioning that the area of the different fins is constant for various dimensions, which guarantees the value of the PCM utilized in various cases to be the same. Figure 10b shows the average temperature development of the PCM for the reverse triangular fins. The temperature increases rapidly due to the conduction heat transfer and then reduces. At 1800 s, the temperature of the PCM for the longest fin is 40 °C, which is higher than the PCM temperature for cases 9, 10, 11 by 3%, 5%, and 7%, respectively. Within 3000 s, the PCM for all the cases reaches almost the same temperature.   Table 4 shows that the cases with the longest fin (base 2.5 mm) are the best cases among all the studied cases, and the case with a reverse triangular fin shows higher performance. This behavior is caused by, firstly, the longer fins, which provide a larger surface area to enhance the heat transfer to the system. Secondly, the reverse triangular fin provides a strong barrier against the circulation of the liquid PCM causing the liquid PCM to circulate within the zone rather than the whole system. Case 7 is known as the worst case due to having the highest melting time and the lowest heat storage rate. The total melting time of the PCM in case 8 is 2997 s, which is shorter than those of cases 4 and 7 by 2.1% and 7.2%, respectively. The heat storage rate of case 8 (56.1 W) is the highest among all other cases, which is greater than those of cases 4 and 7 by 1.9% and 6.5%, respectively. The outline of this section shows that the longest fin offers a higher surface area of heat transfer, resulting in a faster charging process and a higher heat storage value. The reverse triangular fin provides a better barrier against the circulation of the molten PCM, which slightly accelerates the melting process and improves the thermal storage rate. Although there is a small difference between the calculated parameters in cases 4 and 8, case 8 records slightly higher performance, and, consequently, it is selected as the reference case for the next investigations.  Figure 11 compares the liquid fraction development of the reverse triangular fin (Case 8) with that of the no-fin (case 1) and uniform-fin (case 2) cases at different times. The thermal behavior of the charging process was explained previously for case 8. Regarding case 1, the melting process initially starts at the adjacent area of the HTF channel wall. The liquid phase expands during the melting process and then it accumulates at the top region of the system, while the solid part remains at the bottom due to the Boussinesq effect. Within 2400 s, just 80% of the PCM is melted. This can be attributed to the limited heat transfer mechanism provided by the channel walls due to the absence of the fins. Considering case 2, the melting process occurs faster as the fins increase the heat transfer surface area, as well as the average thermal conductivity. Figure 12a shows the liquid fraction development of cases 1, 2, and 8. For the nofin case, the PCM melts logarithmically and needs a relatively long time to be melted. Regarding cases 2 and 8, the melting processes have almost the same behavior. It should be mentioned that case 2 benefits from higher performance compared to case 8 from t = 0 to 2350 s, of which 98% of the PCM in both cases is melted approximately. In fact, the existence of fins causes the charging process of cases 2 and 8 to occur at a higher rate compared to the no-fin case. After t = 2350 s, the melting process of case 8 occurs faster than that of case 2, and the total PCM is melted within 3000 s as illustrated in Figure 12a. Figure 12b shows the variations in the mean temperature of the PCM for cases 1, 2, and 8. Regarding the no-fin case, the mean temperature of the PCM is the lowest, and it does not reach the thermal equilibrium even within 3600 s. For cases 2 and 8, the mean temperature shows approximately the same behavior at various times. The PCM in cases 2 and 8 reach the thermal equilibrium within 3600 s and this indicates that the shape of the fin has a tiny impact on the melting time and the mean temperature distribution. Table 5 presents the total melting time and heat storage rate for cases 1, 2, and 8. It indicates that the total melting time and heat storage rate of cases 2 and 8 are almost similar, while a considerable difference between these values and those of the no-fin case is observed. Case 8 is selected to be the best case due to having the lowest total melting time and highest heat storage rate compared with the other cases. It should be mentioned that the melting times of cases 1 and 2 increase by 55.3 and 2.0% compared with that of case 8. Furthermore, the heat storage rates of cases 1 and 2 reduce by 35.4 and 1.7% compared with that of case 8.  Figure 12a shows the liquid fraction development of cases 1, 2, and 8. For the no-fin case, the PCM melts logarithmically and needs a relatively long time to be melted. Regarding cases 2 and 8, the melting processes have almost the same behavior. It should be mentioned that case 2 benefits from higher performance compared to case 8 from t = 0 to 2350 s, of which 98% of the PCM in both cases is melted approximately. In fact, the existence of fins causes the charging process of cases 2 and 8 to occur at a higher rate compared to the no-fin case. After t = 2350 s, the melting process of case 8 occurs faster than that of case 2, and the total PCM is melted within 3000 s as illustrated in Figure 12a. Figure 12b shows the variations in the mean temperature of the PCM for cases 1, 2, and 8. Regarding the no-fin case, the mean temperature of the PCM is the lowest, and it does not reach the thermal equilibrium even within 3600 s. For cases 2 and 8, the mean temperature shows approximately the same behavior at various times. The PCM in cases 2 and 8 reach the thermal equilibrium within 3600 s and this indicates that the shape of the fin has a tiny impact on the melting time and the mean temperature distribution. Table 5 presents the total melting time and heat storage rate for cases 1, 2, and 8. It indicates that the total melting time and heat storage rate of cases 2 and 8 are almost similar, while a considerable difference between these values and those of the no-fin case is observed. Case 8 is selected to be the best case due to having the lowest total melting time and highest heat storage rate compared with the other cases. It should be mentioned that the melting times of cases 1 and 2 increase by 55.3 and 2.0% compared with that of case 8. Furthermore, the heat storage rates of cases 1 and 2 reduce by 35.4 and 1.7% compared with that of case 8.    To reduce the melting time of the PCM, a flat fin is added to the bottom of the system because, as shown for the systems without the added fin, a portion of the solid part remained at the bottom of the heat storage unit. Two scenarios are considered for the RTF-base 2.5 case as the best model with an added fin. Like the first scenario, a constant value for the base of the reverse triangular fin is considered (2.5 mm) and the height of the triangular fin is then determined to have the same area of the uniform finned case with an added fin. In the second scenario, the height is supposed to be 8 mm, and the base value is calculated accordingly (the base value is obtained at 2 mm). It should be noted that the effects of adding a flat fin at the bottom of the heat exchanger are studied, which was approved in the previously published papers as a noticeable method [57]. Figure 13 shows the melting process of cases 3, 12, and 13. Case 3 stands for the uniform-fin case with an added fin at the bottom. The reversed triangular fin (2.5 mm (base) × 6.4 mm (height)) with an added fin to the bottom is case 12, while case 13 shows the reversed triangular fin (2 mm (base) × 8 mm (height)) with an added fin to the bottom of the unit. The first row of Figure 13 shows the liquid fraction contours of case 3 at various times. The PCM is melted in the areas around the fins, adjacent to the HTF channel's wall, and over the flat fin. The molten area expands gradually, and it is worth mentioning that the solid part cannot be accumulated at the bottom of the unit due to the existence of the flat fin. The second row of Figure 13 shows the melting process of case 12. Within the 1800 s, the PCM totally melts at the lower and the higher zones due to the existence of the flat fin added at the bottom and the circulation occurred through the zones via the wide distance existed between the opposite fins compared with case 13. The third row shows the liquid fraction contours of case 13, which has the longest fin and the smallest gap between the opposite fins. The solid part remains at the top, as well as the bottom zones of the system within 1800 s.

Comparing the Thermal Performance of Reverse Triangular Fins with That of No-Fin and Uniform-Fin Cases
Generally, case 13 shows higher performance, and this is clear from Figure 14, which illustrates the temperature distribution of the system at different times. Within 600 s, the average temperatures of all these cases were approximately the same because of the effects of the conduction process in all cases. The temperature increases around the fins, in the areas adjacent to the wall, and in the flat fin at the bottom. The top section of the system reaches the thermal equilibrium early, then this expands to deeper regions of the PCM. In case 13, the harmonic circulation of melted PCM can be seen in the separate zones and through the entire body system. Figure 15a illustrates the liquid fraction variations of cases 3, 12, and 13 versus time. Figure 15a shows that all the cases have approximately a similar behavior while case 13 benefits from a slight advantage over the other cases. It should be noted that case 13 has the longest fin, and, consequently, the highest surface area of the heat transfer. The average temperatures of the PCM for cases 3, 12, and 13 are also observed in Figure 15b. It is shown that case 13 reaches the thermal equilibrium earlier compared with that of the other cases. From the curves of the liquid fraction and the temperature distribution of cases 3, 12, and 13, it is clear that the thermal behaviors of the studied cases are similar, since the fins are attached to the wall and the flat fin is added to the bottom of the unit. However, small differences are observed, which can be attributed to the fins' shape.  Generally, case 13 shows higher performance, and this is clear from Figure 14, which illustrates the temperature distribution of the system at different times. Within 600 s, the average temperatures of all these cases were approximately the same because of the effects of the conduction process in all cases. The temperature increases around the fins, in the areas adjacent to the wall, and in the flat fin at the bottom. The top section of the system reaches the thermal equilibrium early, then this expands to deeper regions of the PCM. In  Figure 15a illustrates the liquid fraction variations of cases 3, 12, and 13 versus time. Figure 15a shows that all the cases have approximately a similar behavior while case 13 benefits from a slight advantage over the other cases. It should be noted that case 13 has the longest fin, and, consequently, the highest surface area of the heat transfer. The average temperatures of the PCM for cases 3, 12, and 13 are also observed in Figure 15b. It is shown that case 13 reaches the thermal equilibrium earlier compared with that of the other cases. From the curves of the liquid fraction and the temperature distribution of cases 3, 12, and 13, it is clear that the thermal behaviors of the studied cases are similar, since the fins are attached to the wall and the flat fin is added to the bottom of the unit. However, small differences are observed, which can be attributed to the fins' shape.  Table 6 shows the values of the melting time and the heat storage rate of cases 3, 12, and 13. Case 13, having the lowest melting time (1978 s) and the highest heat storage rate (81.5 W), benefits from an advantage over cases 3 and 12. It should be noted that the total melting times of cases 3 and 12 increase by 4.0 and 10.1%, respectively, compared with that of case 13. Furthermore, the heat storage rates of cases 3 and 12 reduce by 4.5 and 8.5%, respectively, compared with that of case 13.  Table 6 shows the values of the melting time and the heat storage rate of cases 3, 12, and 13. Case 13, having the lowest melting time (1978 s) and the highest heat storage rate (81.5 W), benefits from an advantage over cases 3 and 12. It should be noted that the total melting times of cases 3 and 12 increase by 4.0 and 10.1%, respectively, compared with that of case 13. Furthermore, the heat storage rates of cases 3 and 12 reduce by 4.5 and 8.5%, respectively, compared with that of case 13.

Conclusions
In this paper, the impacts of various orientations and dimensions of triangular fins on the charging process of a vertical finned triple-tube PCM heat storage system were studied. Triangular fins were installed on the inner and outer sides of the pipe inside the middle tube in two different ways: (1) the base of the triangular fins was connected to the pipe; (2) the pointed part of the triangular fins was attached to the pipe while the base of the fin was located inside the PCM domain. After finding the best case, a rectangular fin was added to the bottom of the heat storage unit to enhance the heat transfer rate and its thermal performance was compared with that of the uniform-fin case with an added fin at the bottom. The following results were achieved as follows: • The slope of the sides of the triangular fins has a determinant role in the circulation of the PCM inside the heat storage unit.

•
The case with the lowest slope (case with the highest fins height) provided a better barrier to keep the PCM inside the zone between two pairs of neighbor fins, while the case with the shortest fin height accelerated the PCM circulation through the zones.

•
The case with the highest height of the fins, and thus the lowest base, experienced the fastest phase change process due to having the largest surface area of the fin.

•
The reverse triangular fin provided a strong barrier against transferring the molten PCM through the zone, which slightly accelerated the melting process and improved the thermal storage rate.

•
The case with reverse triangular fins and the highest height showed higher performance compared with the other cases, since 98.8% of the entire zone was melted, except the region at the bottom of the unit. • Adding a rectangular fin to the bottom of the heat exchange with reverse triangular fins and the highest height was selected as the best model, resulting in the lowest melting time (1978 s) and the highest heat storage rate (81.5 W).
For future studies, the authors plan to investigate a wider range of studied parameters, the effects of the utilized materials, and the orientation of the tubes. Moreover, with the help of additive manufacturing, the authors are planning to fabricate the studied geometry and experimentally investigate the tringle twisted fin effects.