Numerical study on thermal energy storage performance of phase change material under non-steady-state inlet boundary
Y.B. Tao,Y.L. He
DOI: https://doi.org/10.1016/j.apenergy.2011.04.039
IF: 11.2
2011-01-01
Applied Energy
Abstract:Highlights ► Based on the unstable solar radiation, a model was established for phase change process under unsteady boundary. ► The PCM melting time decreases with the initial inlet temperature increase under the same average inlet temperature. ► The melting time reduces about 51.9% with the initial inlet temperature increase from 30°C to 90°C. ► The melting time decreases with the initial inlet mass flow rate increase under the same average inlet mass flow rate. ► The melting time reduces about 36.5% with the initial inlet mass flow rate increase from 2.0 × 10 −4 kg/s to 8.0 × 10 −4 kg/s. Abstract Due to the solar radiation intensity variation over time, the outlet temperature or mass flow rate of heat transfer fluid (HTF) presents non-steady-state characteristics for solar collector. So, in the phase change thermal energy storage (PCTES) unit which is connected to solar collector, the phase change process occurs under the non-steady-state inlet boundary condition. In present paper, regarding the non-steady-state boundary, based on enthalpy method, a two dimensional physical and mathematical model for a shell-and-tube PCTES unit was established and the simulation code was self-developed. The effects of the non-steady-state inlet condition of HTF on the thermal performance of the PCTES unit were numerically analyzed. The results show that when the average HTF inlet temperature in an hour is fixed at a constant value, the melting time (time required for PCM completely melting) decreases with the increase of initial inlet temperature. When the initial inlet temperature increases from 30 °C to 90 °C, the melting time will decrease from 42.75 min to 20.58 min. However, the total TES capacity in an hour reduces from 338.9 kJ/kg to 211.5 kJ/kg. When the average inlet mass flow rate in an hour is fixed at a constant value, with the initial HTF inlet mass flow rate increasing, the melting time of PCM decreases. The initial inlet mass flow rate increasing from 2.0 × 10 −4 kg/s to 8.0 × 10 −4 kg/s will lead to the melting time decreasing from 37.42 min to 23.75 min and the TES capacity of PCM increasing from 265.8 kJ/kg to 273.8 kJ/kg. Under all the studied cases, the heat flux on the tube surface increases at first, until it reaches a maximum then it decreases over time. And the larger the initial inlet temperature or mass flow rate, the earlier the maximum value appearance and the larger the maximum value. Keywords Enthalpy method Phase change thermal energy storage Non-steady-state boundary Nomenclature a coefficient in Eq. (5) a ′ coefficient in Eq. (6) b coefficient in Eq. (5) b ′ coefficient in Eq. (6) c p specific heat, J kg −1 K −1 f liquid fraction h heat transfer coefficient, W m −2 K −1 k thermal conductivity, W m −1 K −1 L length of the PCM unit, m m ˙ mass flow rate, kg s −1 Pr Prandtl number Q thermal storage capacity, kJ kg −1 r radial coordinate, m R i inner radius of the tube, m R o inner radius of the shell side, m Ra Rayleigh number T temperature, °C t time, s x axial coordinate, m Greek symbols δ melt layer thickness, m μ dynamic viscosity, Pa s θ relative temperature ( T − T m ), K ρ density, kg m −3 υ kinetic viscosity, m 2 s −1 Δ H enthalpy, kJ kg −1 Φ heat flux, W Superscripts ∗ last time layer value Subscripts f heat transfer fluid i initial state in inlet boundary l liquid m melting point out outlet boundary p phase change material s solid 1 Introduction A number of researches on the utilization of solar energy have been carried out due to the globe energy crisis and environmental pollution. However, solar energy has a serious shortcoming that it is unstable and discontinuous with different weathers, times and seasons. So, in order to ensure the solar energy system continuous and stable operation with high efficiency, thermal energy storage (TES) unit becomes a necessary component in the solar thermal utilization systems, such as solar buildings [1,2] , solar water heating systems [3,4] and solar energy generation systems [5,6] . Due to the high energy storage density and constant phase change temperature, phase change thermal energy storage (PCTES) has gradually become the preferred TES pattern. Thermal energy storage with phase change materials has been a main topic in research for the last 20 years. Zalba et al. [7] performed a review on the history of solid–liquid phase change thermal energy storage on materials, heat transfer and applications. Sharma et al. [8] summarized the investigation and analysis of the available thermal energy storage systems incorporating PCMs for different applications. Agyenim et al. [9] carried out a review of materials, heat transfer and phase change problem formulation for latent heat thermal energy storage system over the last three decades. In recent years, a lot of numerical and experimental studies on the characteristics of PCTES process have been performed. In numerical studies, Sharma et al. [10] numerically studied the heat transfer performance of a latent heat storage system with fatty acids as PCMs. The effects of PCMs physical properties, heat exchanger materials and patterns on the PCTES performance were analyzed. A mathematical model regarding the conjugated problem of transient forced convection and solid–liquid phase change heat transfer based on the enthalpy formulation was presented by Trp et al. [11] . The transient heat transfer phenomenon during charging and discharging of the shell-and-tube latent thermal energy storage system was analyzed. Fang and Chen [12] investigated the effects of different multiple PCMs on the melted fraction, stored thermal energy and fluid outlet temperature of the shell-and-tube latent thermal energy storage unit. Guo and Zhang [13] numerically studied the effects of geometry parameters and boundary conditions on the performance of a new type high temperature latent heat thermal energy storage system. Qarnia [14] developed a theoretical model to predict the thermal performance of a solar latent heat storage unit and three kinds of PCM ( n -octadecane, Paraffin wax and Stearic acid) were examined to find the optimum design for a given climatic conditions. Adine and Qarnia [15] numerically studied a latent heat storage unit consisting of a shell-and-tube filled with P116 and n -octadecane. Wu et al. [16] performed the simulations of a packed bed cool thermal energy storage system with n -tetradecane as phase change material. In experimental studies, Trp [17] experimentally analyzed the transient heat transfer characteristics during phase change material melting and solidification. Akgun et al. [18,19] analyzed the latent thermal energy storage system of the shell-and-tube type with three kinds of paraffin as PCMs. A novel tube-in-shell storage geometry was introduced and the effects of the Reynolds number and Stefan number on the melting and solidification behaviors were examined. Long [20] investigated heat transfer performance of a triplex concentric tube thermal energy storage unit. Wang et al. [21] adopted β -Aluminum nitride as additive to enhance the thermal conductivity and thermal performance of phase change materials. Mddrano et al. [22] experimentally evaluated the performance of commercial heat exchanger used as PCM thermal storage systems. The foregoing literature review shows that a lot of studies have been performed on the shell-and-tube PCTES unit. But most of them were performed under the constant HTF inlet temperature and mass flow rate. However, for the solar thermal utilization system, the HTF comes from the solar collector, the HTF inlet conditions for the PCTES unit depends on the HTF outlet conditions at the solar collector. Due to the unstable and discontinuous characteristics of the solar radiation intensity and the complicated heat transfer and fluid flow process in the collector [23–25] , the HTF outlet condition presents obvious non-steady-state. So, the inlet conditions (HTF temperature or mass flow rate) for the PCTES unit inevitably present non-steady-state characteristics, which will greatly affect the melting rate and stored thermal energy of PCM and the outlet temperature of HTF. In present study, regarding the non-steady-state characteristics of HTF temperature or mass flow rate at the inlet of the PCTES unit, based on enthalpy method, a two dimensional physical and mathematical model for the phase change process in water/ n -octadecane shell-and-tube PCTES unit was established and the simulation code was self-developed. Then the numerical study on the performance of the PCTES unit was performed and the effects of the non-steady-state HTF inlet conditions on the PCTES performance were analyzed. 2 Physical model and governing equations 2.1 Physical model The physical model for PCTES unit is shown in Fig. 1 , which is a shell-and-tube configuration. Water is selected as HTF and flows in inner tube. The shell side is full of PCM ( n -octadecane). In the charging process, thermal energy transfers from HTF to PCM through the tube surface and is stored in PCM, and in discharging process, thermal energy stored in PCM releases to HTF. The length for the computation domain ( L ) is 1.0 m, the inner radius for the tube ( R i ) is 6.35 mm, inner radius for shell side ( R o ) is 11.35 mm. The thickness of tube wall is neglected. The thermophysical properties for HTF and PCM are shown in Table 1 . In order to simplify the physical and mathematical model, the following assumptions are adopted [15] . (1) The axial heat conduction and viscous dissipation in the HTF is negligible. The flow of HTF is treated as one dimensional fluid flow. The phase change material region is treated as an axisymmetric model. (2) The effect of natural convection of PCM during melting is taken into account with an effective thermal conductivity of the liquid phase of PCM, as shown in: (1) k k l = CRa n δ R - R i m where m = 0.8, n = 0.25, C dependents on the inlet temperature of the HTF [14,15] . (2) C = 0.24 θ f,in ⩽ 10 0.18 10 < θ f,in ⩽ 20 0.16 20 < θ f,in 2.2 Governing equations The enthalpy method is adopted to deal with the moving boundary problem in solid–liquid phase change process. The corresponding governing equations are shown as follows [15] . For the HTF (water) (3) ∂ θ f ∂ t = - A ∂ θ f ∂ x - B ( θ f - θ ∗ ) where θ f = T f − T m ; A = m ˙ f ρ f π R i 2 , B = 2 h ( ρ c p ) f R i . For the PCM ( n -octadecane) (4) ( ρ c p ) p ∂ θ ∂ t = ∂ ∂ x k p ∂ θ ∂ x + 1 r ∂ ∂ r rk p ∂ θ ∂ r - ρ p Δ H ∂ f ∂ t where θ = T − T m , f is the melting fraction of the PCM. The energy Eq. (4) is formulated by the enthalpy method. And the melting fraction is determined as: f = 0 , θ < 0 0 < f < 0 , θ = 0 f = 1 , θ > 0 2.3 Initial conditions and boundary conditions The initial conditions are For HTF: θ f ( x , t = 0 ) = T i - T m , For PCM: θ ( x , r , t = 0 ) = T i - T m . The boundary conditions are For PCM: inlet and outlet boundaries ∂ θ ( x = 0 , r , t ) ∂ x = ∂ θ ( x = L , r , t ) ∂ x = 0 , inner surface boundary - k p ∂ θ ( x , r = R i , t ) ∂ r = U ( θ f ( x , t ) - θ ( x , r = R i , t ) ) , outer surface boundary ∂ θ ( x , r = R o , t ) ∂ r = 0 . For HTF: inlet boundary θ f ( x = 0 , r , t ) = θ f,in = f ( t ) , outer surface boundary - k p ∂ θ ( x , r = R i , t ) ∂ r = U ( θ f ( x , t ) - θ ( x , r = R i , t ) ) . 2.4 Model validation The finite volume method (FVM) is used to discrete the governing equations. Due to the coupled energy transfer process between the HTF and PCM, the energy equations both for HTF and PCM are integrated solving in the whole computational domain. After the grid checking, a grid number 100( x ) × 40( r ) is adopted as the final computational grid. In the solving process, the time step is set to 5 s. The convergence criterion is that the maximum relative heat balance deviation between HTF and PCM is less than 1.0 × 10 −3 . In order to validate the reliability of the physical model and simulation code, the comparisons between the present numerical predictions and the experimental data [28] were performed under the same geometric parameters and operating conditions. In the experiment, two copper-constantan thermocouples are installed inside the PCM at locations 1 ( x = 0.95 m, r = 0.001 m) and 2 ( x = 0.51 m, r = 0.002 m). The comparative results are shown in Fig. 2 . The predicted temperatures for the two reference points are good accordant with the experimental results in Ref. [28] both under θ f in = 10 K and θ f in = 20 K conditions. The good agreements show that the physical model and simulation code in the present paper are correct and reliable. In solar thermal utilization system, affected by the unstable solar radiation intensity, the energy gained by HTF in solar collector will present unsteady characteristics. So, if the HTF mass flow rate is fixed, the HTF outlet temperature for solar collector (inlet temperature for PCTES unit) must be unsteady. If we want to control the temperature at a constant value, the mass flow rate of HTF must be unsteady. In the following sections, the numerical analysis will be performed both for the above two conditions: temperature variation with time and mass flow rate kept constant; mass flow rate variation with time and temperature kept constant. 3 Effects of non-steady-state HTF inlet temperature In order to analyze the effect of the non-steady-state HTF inlet temperature condition on the thermal energy storage performance of PCM, the inlet temperature is assumed as linearly changing with time. So, the inlet temperature should satisfy the following linear equation. (5) T f,in = b + t / a where T f,in is the inlet temperature of the HTF, unit °C; t is time, unit seconds. a and b are the coefficients for the linear equation. a represents the variation rate, b presents the initial inlet temperature value. In order to satisfy the average inlet temperature in an hour fixed at 60 °C, the detailed values for a and b are shown in Table 2 . Different a and b values present different temperature variation tendencies. For case 1–3, the HTF inlet temperature increases over time and the increasing tendency weakens with value b augment. For case 4–6, the HTF inlet temperature decreases over time and the decreasing tendency augments with value b increase. In the following sections, using different initial inlet temperatures presents different inlet temperature variation tendencies. The mass flow rate for the HTF is fixed at 5.0 × 10 −4 kg/s and the initial temperature both for PCM and HTF is 25 °C. 3.1 Effects of b on melting time and melting fraction Fig. 3 shows the effects of initial HTF inlet temperature ( b ) on PCM melting time. With the increase of initial inlet temperature, in the early, the temperature difference between HTF and PCM augments, which leads to the heat transfer rate increasing. So, the time required for the PCM complete melting is reduced. When the average inlet temperature is fixed at 60 °C and the initial inlet temperature ( b ) increasing from 30 °C to 90 °C, the melting time will decrease from 42.75 min to 20.58 min, which reduces about 51.9%. The effect of initial HTF inlet temperature on melting fraction is shown in Fig. 4 . The melting fraction quickly increases over time until it reaches 1 which means the PCM completely melted. And the higher initial inlet temperature, the larger slope of the melting fraction over time, which presents the faster melting rate of the PCM. So, the melting time as shown in Fig. 3 will quickly decrease with initial inlet temperature increase. 3.2 Effects of b on TES capacity and solid–liquid interface Fig. 5 presents the effect of initial inlet temperature on the total thermal energy storage capacity in an hour. Although with the increase of b , the melting rate increases, the thermal storage capacity in an hour is reduced. This is because that when b is larger, the increasing degree of inlet temperature over time will be weakened for case 1–3. And for case 4–6, the inlet temperature will quickly decrease over time, which even leads to the presentation of PCM releasing energy to HTF in the late. So when the initial inlet temperature increases from 30 °C to 90 °C, the thermal storage capacity of PCM in an hour reduces from 338.9 kJ/kg to211.5 kJ/kg, which reduces about 37.6%. The effects of initial inlet temperature on solid–liquid interface is shown in Fig. 6 at t = 12 min. With the increase of initial inlet temperature, the melting rate increases and the melting fraction at the same time step augments. So, at the same time step, the location of solid–liquid interface moves to the outer surface of shell side with b increase. 3.3 Effects of b on surface heat flux and HTF outlet temperature Fig. 7 shows the effect of initial inlet temperature on the total surface heat flux of PCTES unit. In all the studied cases, the heat flux increases over time at first, until it reaches a maximum value, then it will decrease. And the larger the initial inlet temperature, the earlier the maximum value presentation and the larger the maximum value. In the early, only the PCM near the inlet region tube surface begins to melt and the thermal energy storage rate is slow. So the heat flux on the total tube surface is lower. Then the PCM near to the whole tube surface starts to melt, the thermal energy storage rate increases which leads to the heat flux quickly augment and gradually tending to a maximum. After that, with the solid–liquid interface moving to the inner of PCM unit, the heat conduction distance augments and thermal resistance increases, which results in the heat storage rate slowing down. And the heat flux decreases again. For case 4–6 ( b = 70, 80, 90 °C), the inlet temperature of HTF decreases over time. In the late, the HTF temperature is lower than PCM temperature. The heat flux at the tube surface becomes to negative value, which means the PCM releases thermal energy to the HTF and presents a discharging process. So, the total thermal energy storage capacity in an hour will reduces with the increase of initial inlet temperature as shown in Fig. 5 . The effects of initial HTF inlet temperature on outlet temperature is shown in Fig. 8 . In the early, only a little of PCM melts and the heat storage capacity is smaller. The HTF outlet temperature quickly increases over time. Then with the melting fraction increase, the heat storage capacity augments. Affected by the inlet temperature variation patterns, the outlet temperature presents different variation tendencies. For the case 4–6, the inlet temperature decreases over time. The outlet temperature quickly decreases at first. Then with the solid–liquid interface moving to the inner of PCM, the heat conduction resistance increases and heat storage capacity rate decreases. The HTF outlet temperature augments in intermediate time. At last, with the more reduction of inlet temperature, the outlet temperature quickly decreases. But for case 1–3, the inlet temperature increases over time. So, although the heat storage capacity augments over time, the outlet temperature only has a little decrease for case 3 and slowly increases for case 1 and case 2. Then with the solid–liquid interface moving to the inner of PCM, the heat storage rate decreases. At the same time, the inlet temperature linearly increases over time. So the outlet temperature quickly augments. 4 Effects of non-steady-state HTF inlet mass flow rate In the following sections, the inlet mass flow rate of HTF is assumed as linearly changing with time. So, the inlet mass flow rate should satisfy the following linear equation. (6) m ˙ f = ( b ′ + t / a ′ ) × 10 - 4 where m ˙ f is the inlet mass flow rate of the HTF, unit kg/s; t is time, unit seconds. a ′ and b ′ are the coefficients for the linear equation. a ′ represents the variation rate and b ′ presents the initial value. In order to satisfy the average inlet mass flow rate in an hour fixed at 5.0 × 10 −4 kg/s, the detailed values a ′ and b ′ are shown in Table 2 . Different a ′ and b ′ values present different mass flow rate variation tendencies. For case 1–3, the HTF inlet mass flow rate increases over time and the larger value b ′, the weaker increasing tendency. For case 4–6, the HTF inlet mass flow rate decreases over time and the larger value b ′, the quicker decreasing tendency. In the following sections, using different initial mass flow rate ( b ′) presents different mass flow rate variation tendencies. The inlet temperature for the HTF is fixed at 60 °C and the initial temperature is 25 °C both for PCM and HTF. 4.1 Effects of b ′ on melting time and melting fraction Fig. 9 shows the effects of initial HTF inlet mass flow rate on PCM melting time. With the increase of initial inlet mass flow rate, in the early, the thermal energy taken by HTF augments. So, the decreasing degree of HTF temperature in the PCTES unit slows down and the temperature difference between HTF and PCM augments. At the same time, the convection heat transfer coefficient increases. As a results, the heat transfer rate between PCM and HTF increases, which leads to the time required for PCM completely melting reduction. When the average inlet mass flow rate is fixed at 5.0 × 10 −4 kg/s, the initial inlet mass flow rate ( b ′) increasing from 2.0 × 10 −4 kg/s to 8.0 × 10 −4 kg/s will lead to the melting time decreasing from 37.42 min to 23.75 min, which reduces about 36.5%. The effect of initial HTF inlet mass flow rate on melting fraction is shown in Fig. 10 . The melting fraction quickly increases over time. And the higher the initial inlet mass flow rate, the larger slope the melting fraction over time, which presents the faster PCM melting rate. So, the melting time as shown in Fig. 9 will quickly decrease with the increase of initial HTF inlet mass flow rate. 4.2 Effects of b ′ on TES capacity and solid–liquid interface Fig. 11 presents the effect of initial inlet mass flow rate on thermal energy storage capacity of PCM in an hour. With the increase of initial inlet mass flow rate, the TES capacity in an hour is augmented which is different from the effect of initial inlet temperature. This is because the inlet temperature for HTF is constant at different inlet mass flow rate conditions. And the temperature for PCM increases over time. So the temperature difference between PCM and HTF decreases over time. When b ′ is larger, the HTF mass flow rate in the early is larger and the HTF temperature is slowly decreased along the tube. The heat transfer coefficient and temperature difference between PCM and HTF are larger and the heat transfer rate is larger. This means that in the early, the larger initial inlet mass flow rate, the larger TES capacity. In the late, with the increase of PCM temperature, the temperature difference decreases. For case 4–6 the mass flow rate decreases which leads to the TES rate quickly decreasing. And for case 1–3, the mass flow rate increases, which leads to the TES rate slowly decreasing. So, in totally, the TES capacity augments with b ′ increase, but the augmentation degree is limited. When the initial inlet mass flow rate increases from 2.0 × 10 −4 kg/s to 8.0 × 10 −4 kg/s, the TES capacity of PCM in an hour augments from 265.8 kJ/kg to 273.8 kJ/kg, which augments only about 3%. The effects of initial inlet mass flow rate on solid–liquid interface is shown in Fig. 12 at t = 12 min. With the increase of initial inlet mass flow rate, the melting rate increases and the melting fraction at the same time step augments. The location of solid–liquid interface moves to the outer surface of shell side. 4.3 Effects of b ′ on surface heat flux and HTF outlet temperature Fig. 13 shows the effects of initial inlet mass flow rate on the total surface heat flux of PCTES unit. In all the studied cases, the heat flux increases over time at first, until it reaches a maximum value, then it will decrease. And the larger the initial inlet mass flow rate, the earlier the maximum value presentation and the larger the maximum value. In the early, only the PCM near the inlet region tube surface begins to melt and the TES rate is slow. So the heat flux on the total tube surface is less. Then the PCM near to the total tube surface starts to melt, the TES rate increases which leads to the heat flux quickly augment and gradually tending to a maximum value. In the late of the melting process, with the solid–liquid interface moving to the inner of PCM, the heat conduction distance augments and thermal resistance increases, which results in the TES rate slowing down. So the heat flux quickly decreases. Due to the HTF inlet temperature is fixed at 60 °C, which is always higher than PCM, so the surface heat flux is always a plus value. But at the last stage, with the completely melting of PCM, the PCM temperature quickly increases. The temperature difference between PCM and HTF dramatically decreases, so the surface heat flux nears to zero. The effect of initial HTF inlet mass flow rate on the outlet temperature is shown in Fig. 14 . In the early, only a little of PCM melts and the heat storage capacity is small. The HTF outlet temperature quickly increases over time and reaches a local maximum. And the larger the initial inlet mass flow rate, the earlier the maximum value presentation and the larger the maximum value. Then with the melting fraction increase, the heat storage capacity augments. The outlet temperature of HTF slowly decreases. In the late of the melting process, with the solid–liquid interface moving to the inner of PCM, the heat conduction resistance increases, which leads to heat storage capacity decreasing. So, the HTF outlet temperature augments again. At last with the completely melting of PCM, the temperature for PCM quickly increases and the temperature difference between PCM and HTF decreases to zero, the HTF outlet let temperature nears to a constant value which equals the inlet temperature, 60 °C. From the figure, it also can be seen that with the increase of initial inlet mass flow rate, the outlet temperature for HTF increases. 5 Conclusions In present paper, regarding that in solar energy utilization process, the temperature or mass flow rate of HTF presents non-steady-state characteristics at the inlet of the PCTES unit, a two dimensional physical and mathematical model for the phase change process in water/ n -octadecane shell-and-tube latent thermal energy storage unit was established based on enthalpy method. The effects of the non-steady-state inlet conditions of HTF on melting time, melting fraction, TES capacity, solid–liquid interface, heat flux on tube surface and HTF outlet temperature were analyzed. The following conclusions can be derived. 1. The PCM melting time decreases with the increase of initial HTF inlet temperature, when the average inlet temperature in an hour is fixed. In the studied conditions, the initial inlet temperature increasing from 30 °C to 90 °C will lead to the melting time decreasing from 42.75 min to 20.58 min, which reduces about 51.9%. However, the total TES capacity in an hour reduces from 338.9 kJ/kg to 211.5 kJ/kg, which reduces about 37.6%. 2. The PCM melting time decreases with the increase of initial HTF inlet mass flow rate, when the average inlet mass flow rate in an hour is fixed. In the studied conditions, the initial inlet mass flow rate increasing from 2.0 × 10 −4 kg/s to 8.0 × 10 −4 kg/s will lead to the melting time decreasing from 37.42 min to 23.75 min, which reduces about 36.5%. At the same time, the TES capacity of PCM in an hour augments from 265.8 kJ/kg to 273.8 kJ/kg, which augments about 3%. 3. The heat flux on the tube surface increases at first, until it reaches a maximum then it decreases over time under all the studied cases. And the larger the initial inlet temperature or mass flow rate, the earlier the maximum value appearance and the larger the maximum value. Acknowledgments The present work is supported by the National Natural Science Foundation of China (Nos. 50906070 and 50736005) and the National Basic Research Program of China (973 Program) (2010CB227102). References [1] W. Xiao X. Wang Y.P. Zhang Analytical optimization of interior PCM for energy storage in a lightweight passive solar room Appl Energy 86 2009 2013 2018 [2] F. Kuznik J. Virgone Experimental assessment of a phase change material for wall building use Appl Energy 86 2009 2038 2046 [3] C. Garnier J. Currie T. Muneer Integrated collector storage solar water heater: temperature stratification Appl Energy 86 2009 1465 1469 [4] K. Sutthivirode P. Namprakai N. Roonprasang A new version of a solar water heating system coupled with a solar water pump Appl Energy 86 2009 1423 1430 [5] D. Brosseau J.W. Kelton D. Ray M. Edgar Testing of thermocline filler materials and molten-salt heat transfer fluids for thermal energy storage systems in parabolic trough power plants J Solar Energy – Trans ASME 127 2005 109 116 [6] Z. Yang S.V. Garimella Molten-salt thermal energy storage in thermoclines under different environmental boundary conditions Appl Energy 87 2010 3322 3329 [7] B. Zalba J.M. Marin L.F. Cabeza H. Mehling Review on thermal energy storage with phase change: materials heat transfer analysis and application Appl Therm Eng 23 2003 251 283 [8] A. Sharma V.V. Tyagi C.R. Chen D. Buddhi Review on thermal energy storage with phase change materials and applications Renew Sustain Energy Rev 13 2009 318 345 [9] F. Agyenim N. Hewitt P. Eames M. Smyth A review of materials, heat transfer and phase change problem formulation for latent heat thermal energy storage systems Renew Sustain Energy Rev 14 2010 615 628 [10] A. Sharma L.D. Won D. Buddhi J.U. Park Numerical heat transfer studies of the fatty acids for different heat exchanger materials on the performance of a latent heat storage system Renew Energy 30 2005 2179 2187 [11] A. Trp K. Lenic B. Frankovic Analysis of the influence of operating condition and geometric parameters on heat transfer in water-paraffin shell-and-tube latent thermal energy storage unit Appl Therm Eng 26 2006 1830 1839 [12] M. Fang G.M. Chen Effects of different multiple PCMs on the performance of a latent thermal energy storage system Appl Therm Eng 27 2007 994 1000 [13] C.X. Guo W.J. Zhang Numerical simulation and parametric study on new type of high temperature latent heat thermal energy storage system Energy Convers Manage 49 2008 919 927 [14] H.El. Qarnia Numerical analysis of a coupled solar collector latent heat storage unit using various phase change materials for heating the water Energy Convers Manage 50 2009 247 254 [15] H.A. Adine H.El. Qarnia Numerical analysis of the thermal behaviour of a shell-and-tube heat storage unit using phase change materials Appl Math Model 33 2009 2134 2144 [16] S. Wu G. Fang X. Liu Thermal performance simulations of a packed bed cool thermal energy storage system using n -tetradecane as phase change material Int J Therm Sci 49 2010 1752 1762 [17] A. Trp An experimental and numerical investigation of heat transfer during technical grade paraffin melting and solidification in a shell-and-tube latent thermal energy storage unit Solar Energy 79 2005 648 660 [18] M. Akgun O. Aydin K. Kaygusuz Experimental study on melting/solidification characteristics of a paraffin as PCM Energy Convers Manage 48 2007 669 678 [19] M. Akgun O. Aydin K. Kaygusuz Thermal energy storage performance of paraffin in a novel tube-in-shell system Appl Therm Eng 28 2008 405 413 [20] J. Long Numerical and experimental investigation for heat transfer in triplex concentric tube with phase change material for thermal energy storage Solar Energy 82 2008 977 985 [21] W. Wang X. Yang Y. Fang J. Ding J. Yan Enhanced thermal conductivity and thermal performance of form-stable composite phase change materials by using β -aluminum nitride Appl Energy 86 2009 1196 1200 [22] M. Mddrano M.O. Yilmza M. Nogues I. Martorell J. Roca L.F. Cabeza Experimental evaluation of commercial heat exchangers for use as PCM thermal storage systems Appl Energy 86 2009 2047 2055 [23] Y.B. Tao Y.L. He Numerical study on coupled fluid flow and heat transfer process in parabolic trough solar collector tube Solar Energy 84 2010 1863 1872 [24] Y.L. He J. Xiao Z.D. Cheng Y.B. Tao A MCRT and FVM coupled simulation method for energy conversion process in parabolic trough solar collector Renew Energy 36 2011 976 985 [25] Z.D. Cheng Y.L. He J. Xiao Y.B. Tao R.J. Xu Three-dimensional numerical study of heat transfer characteristics in the receiver tube of parabolic trough solar collector Int Commun Heat Mass Transfer 37 2010 782 787 [26] F.P. Incropera D.P. DeWitt Fundamentals of heat and mass transfer 4th ed. 1996 Wiley New York [27] H.E.S. Fath Heat exchanger performance for latent heat thermal energy storage system Energy Convers Manage 31 1991 149 155 [28] M. Lacroix Numerical simulation of a shell and tube latent thermal energy storage unit Solar Energy 50 1993 357 367