Melting of NEPCM within a Cylindrical Tube: Numerical Study Using the Lattice Boltzmann Method

This article presents a numerical investigation on melting of phase change material (PCM) enhanced by nanoparticles inside a cylindrical tube using the lattice Boltzmann method. Water (ice) and copper particles are chosen as the base fluid (PCM) and nanoparticles, respectively. Results show that the melting rate is the same for all regions of the cylinder for a low Rayleigh number, while it intensifies at the top half of the cylinder for a moderate Rayleigh number. Also, existence of strong unstable flow in the bottom portion of the cylinder at a Rayleigh number of 10⁶ causes the melting rate to keen after a definite time. Nanoparticles have no significant effect on the melting rate at the beginning of melting, where the conduction mode of heat transfer dominates between the hot wall and solid PCM, while full melting of PCM occurs earlier by the increase of solid concentration.     

Investigation of heat transfer during melting of a PCM by a U‐shaped heat source

An experimental study was conducted to investigate the melting process of a phase change material (PCM) and the associated convection heat transfer due to a U‐shaped heat source embedded in the PCM. The experiments were conducted at four input heat fluxes that varied from 3450 to 5840 W/m². The results showed that the heat transfer behavior, interface movement, and the heat transfer coefficients differed both axially and vertically inside the chamber. The local convective heat transfer coefficient in the inner region, enclosed by the U‐tube, was found to be about 35% higher than that in the outer region over the input heat flux range, resulting in faster melting in the inner region than in the outer region. As melted domain grew vertically from 15% to 100%, it was observed that the overall h in the inner region increased by 40–55% from the lowest to highest heat flux. The melting rate was also found comparatively high up until 65–70% of the total PCM volume melted because of the higher contribution from the inner region. It was also observed that the Rayleigh number increased by approximately 23% in the inner region and 18% in the whole domain as the heat flux increased from 3450 to 5840 W/m². A new Nusselt–Rayleigh number correlation is proposed for the heat transfer during the melting process due to a U‐shaped heat source. Copyright © 2017 John Wiley & Sons, Ltd.     

Investigation on the thermal performance of phase change material/porous medium‐based battery thermal management in pore scale

With the latent heat, the phase change material (PCM) is widely used in battery thermal management (BTM) to control the temperature. In this paper, the porous medium is employed to enhance the heat transfer of PCM. The lattice Boltzmann model for PCM/porous medium in pore scale is considered, where the mesh system with porous medium (fixed point) is generated by quartet structure generation set (QSGS) method. The effects of the Rayleigh number and porosity on the heat transfer process in BTM are investigated. The results show that decreasing the porosity will accelerate the melting rate. When the porosities are 0.9, 0.8, 0.7, and 0.6, the total melting times are decreased by 23.7%, 43.3%, 58.0%, and 75.4%, compared with pure PCM. The heat is transferred through the high‐conductivity framework. The natural convection in the porous medium is weak, and the conduction is the dominated heat transfer. As a result, the area of solid–liquid interface will be increased, and the heat‐transferred rate is accelerated. However, when the Rayleigh number is raised to 10⁵, applying the porous medium with porosity of 0.9 will increase the total melting time, resulted from the stronger natural convection of PCM. The present study is helpful for design of PCM/porous medium‐based BTM.     

Heat transfer enhancement of fatty acids when used as PCMs in thermal energy storage

Phase change materials (PCM) used in latent heat storage systems usually have very low thermal conductivities. This is a major drawback in maintaining the required heat exchange rate between PCM and heat transfer fluid. This paper investigates the enhancement of the heat transfer between PCM and heat transfer fluid, using high thermal conductivity as additives like stainless steel pieces, copper pieces and graphite–PCM composite material. In the experiments, palmitic–lauric acid (80:20) (PL) and stearic–myristic acid (80:20) (SM) were used as PCMs. Test results show that heat transfer enhancement of copper pieces was better at 0.05 Ls^?1 flow rate compared to 0.025 Ls^?1. Using copper as an additive increased the heat transfer rate 1.7 times for melting and 3.8 times for freezing when flow rate was 0.050 Ls^?1. Decreasing the flow rate from 0.050 to 0.025 Ls^?1, increased the melting times 1.3 times and freezing times 1.8 times, decreasing heat transfer rates accordingly. The best result of heat transfer enhancement was observed for the PCM–graphite composite. However, changing the flow rate did not affect the heat transfer rate when graphite was used as additive. Copyright © 2007 John Wiley & Sons, Ltd.     

Accelerated melting of PCM in a multitube annulus‐type thermal storage unit using lattice Boltzmann simulation

Melting of the ice/water as the phase change material in a horizontal single‐tube annulus is sluggish when the stable stratification exists at the bottom of the configuration. To obviate this problem, three heat transfer enhancement techniques could be implemented using the enthalpy‐based lattice Boltzmann method with the double distribution function method to accelerate the process. The multifarious arrangements of the tubes in this horizontal annulus are investigated to expand the region affected by the natural convection. Also, the dispersion of the Cu nanoparticles in the base PCM could boost the thermal conductivity and melting rate. Finally, the metallic porous matrix made of nickel–steel alloys and saturated with the base PCM could be used to enhance the thermal conductivity of the base PCM. The solid–liquid phase change process is defined as the constrained melting of ice‐water in the tube heating mode. There is a thermal equilibrium between ice/water and the nickel–steel porous matrix and the Cu nanoparticles. The Prandtl number, Stefan number, Rayleigh number, and Darcy number are 6.2, 1, 10⁴–10⁵, and 10^?3, respectively. The volumetric concentric of the nanoparticles is between 0 and 0.02 and the porosity ranging from 1 to 0.9 in the representative elementary volume scale.     

Melting in a side heated tall enclosure by a uniformly dissipating heat source

Melting of an organic phase change material (PCM) n-triacontane (C₃₀H₆₂) in a side heated tall enclosure of aspect ratio 10, by a uniformly dissipating heat source has been studied computationally and experimentally. While heat transfer data for melting in enclosures under isothermal wall boundary condition are available in the literature, other boundary conditions, such as constant heat flux often arise in applications of PCM for transient thermal management of electronics. An implicit enthalpy–porosity approach was utilized for computational modeling of the melting process. Experimental visualization of melt front locations was performed. Comparisons between experimental and computational heat transfer data and melt interface locations were good. Fluid flow and heat transfer characteristics during melting suggested that natural convection plays a dominant role during initial stages of melting. At later times, the strength of natural convection diminishes as melting is completed. Correlations of heat transfer rate and melt fraction with time were obtained.     

Crystallization of supercooled spherical nodules in a flow

The general context of this work concerns the latent heat storage (cold storage) by encapsulated spherical nodules of 7.7 cm in diameter containing a Phase Change Material (PCM) and filling a cylindrical tank (about 2500 nodules per m³). During the storage process, a cold heat transfer liquid flows through the tank to crystallize the PCM inside the nodules. The stored energy is released when a hotter heat transfer liquid flows through the tank to induce the melting of the PCM. As the velocity in the tank is relatively small (a few mm·s⁻¹), natural convection is expected on the coolant when the latent heat is released at the crystallization of the PCM. The present study concerns a single nodule surrounded by a flowing heat transfer liquid at a temperature T_C lower than the melting temperature T_F. This text presents an enthalpic modelling of the phase change inside the nodule coupled with a CFD simulation of the external flow to describe the mutual influence of the natural convection and the kinetics of crystallization. This study is afterwards extended to the case of two superposed nodules to investigate the influence of the crystallization of the lower nodule on the upper one.     

Fluid flow and heat transfer in PCM panels arranged vertically and horizontally for application in heating systems

The rectangular panel is within the most common geometries in phase change materials (PCM). Nevertheless, there is a lack of knowledge regarding how the arrangement (vertical or horizontal) and thickness affect heat flux and phase change duration. Such knowledge would be very helpful both in the design of PCM heat exchangers. This paper studies the behavior of the RT60 paraffin, Pr ≈ 330, phase change temperature between 53 and 61 °C; and the fatty acid Palmitic Acid, Pr ≈ 110, phase change temperature of 65 °C. Parametric studies analyzing the influence of: temperature of the walls (5 × 10⁵ < Ra ≤ 2 × 10⁶), phase change processes (melting and solidification), panel position and PCM thickness (aspect ratio1/20 and 3/20) are performed. The flow behavior is analyzed over time using velocity plots and volumetric liquid fraction contours. The formation of Bénard cells and their evolution is described. Free convection dependence of dimensionless string function and Rayleigh number is discussed. Free convection during melting for horizontal panels becomes very important and the mean heat fluxes increase up to twofold compared to vertical panels. For solidification, however, conduction becomes more relevant. The importance of both mechanisms is highlighted by calculating heat transfer rates.     

Melting of nano-enhanced phase change material in a cavity heated sinusoidal from below: Numerical study using lattice Boltzmann method

Natural convection and melting of ice as a phase change material dispersed with copper nanoparticles are numerically investigated. Square cavity filled with nano-mixture (Cu−ice) subjected to sinusoidal temperature distributions from the hot bottom boundary. The phase change process and heat transfer are formulated and solved using the enthalpy-based lattice Boltzmann method. Home-built numerical code is developed and validated. The effect of Rayleigh number (Ra = 10⁴, 10⁵, and 10⁶) and copper nanoparticle concentration (ϕ = 0%, 1%, 3%, and 5%) on the flow characteristics and thermal performance of NePCM during the melting process is examined. According to the numerical results, the melting and charging times decrease by increasing the Rayleigh number. It is also observed that increasing the volume fraction of nanoparticle decrease melting time by up to 10%.     

Thermal and heat transfer characteristics in a latent heat storage system using lauric acid

The thermal and heat transfer characteristics of lauric acid during the melting and solidification processes were determined experimentally in a vertical double pipe energy storage system. In this study, three important subjects were addressed. The first one is temperature distributions and temporal temperature variations in the radial and axial distances in the phase change material (PCM) during phase change processes. The second one is the thermal characteristics of the lauric acid, which include total melting and total solidification times, the nature of heat transfer in melted and solidified PCM and the effect of Reynolds and Stefan numbers as inlet heat transfer fluid (HTF) conditions on the phase transition parameters. The final one is to calculate the heat transfer coefficient and the heat flow rate and also discuss the role of Reynolds and Stefan numbers on the heat transfer parameters. The experimental results proved that the PCM melts and solidifies congruently, and the melting and solidification front moved from the outer wall of the HTF pipe (HTFP) to the inner wall of the PCM container in radial distances as the melting front moved from the top to the bottom of the PCM container in axial distances. However, it was difficult to establish the solidification proceeding at the axial distances in the PCM. Though natural convection in the liquid phase played a dominant role during the melting process due to buoyancy effects, the solidification process was controlled by conduction heat transfer, and it was slowed by the conduction thermal resistance through the solidified layer. The results also indicated that the average heat transfer coefficient and the heat flow rate were affected by varying the Reynolds and Stefan numbers more during the melting process than during the solidification process due to the natural convection effect during the melting process.     

Conjugate natural convection in a square cavity with finite thickness horizontal walls

The effect of conduction of horizontal walls on natural convection heat transfer in a square cavity is numerically investigated. The vertical walls of the cavity are at different constant temperatures while the outer surfaces of horizontal walls are insulated. A code based on vorticity–stream function is written to solve the governing equations simultaneously over the entire computational domain. The dimensionless wall thickness of cavity is taken as 0.1. The steady state results are obtained for wide ranges of Rayleigh number (10³ < Ra < 10⁶) and thermal conductivity ratio (0 < K < 50). The variation of heat transfer rate through the cavity and horizontal walls with Rayleigh number and conductivity ratio is analyzed. It is found that although the horizontal walls do not directly reduce temperature difference between the vertical walls of cavity, they decrease heat transfer rate across the cavity particularly for high values of Rayleigh number and thermal conductivity ratio. Heatline visualization technique is a useful application for conjugate heat transfer problems as shown in this study.     

Experimental study of a compact PCM solar collector

The performance of a compact phase change material (PCM) solar collector based on latent heat storage was investigated. In this collector, the absorber plate–container unit performs the function of both absorbing the solar energy and storing PCM. The solar energy was stored in paraffin wax, which was used as a PCM, and was discharged to cold water flowing in pipes located inside the wax. The collector's effective area was assumed to be 1 m² and its total volume was divided into 5 sectors. The experimental apparatus was designed to simulate one of the collector's sectors, with an apparatus-absorber effective area of 0.2 m². Outdoor experiments were carried out to demonstrate the applicability of using a compact solar collector for water heating. The time-wise temperatures of the PCM were recorded during the processes of charging and discharging. The solar intensity was recorded during the charging process. Experiments were conducted for different water flow rates of 8.3–21.7 kg/h. The effect of the water flow rate on the useful heat gain (Qu) was studied. The heat transfer coefficients were calculated for the charging process. The propagation of the melting and freezing front was also studied during the charging and discharging processes. The experimental results showed that in the charging process, the average heat transfer coefficient increases sharply with increasing the molten layer thickness, as the natural convection grows strong. In the discharge process, the useful heat gain was found to increase as the water mass flow rate increases.     

Natural convection heat transfer of microemulsion phase-change-material slurry in rectangular cavities heated from below and cooled from above

An experimental study has been conducted in dealing with natural convection heat transfer characteristics of microemulsion slurry in rectangular enclosures. The microemulsion slurry used in the present experiment was composed of water, surfactant, and fine particles of phase-change-material (PCM). The PCM mass concentration of the microemulsion slurry was varied from a maximum 30 mass% to a diluted minimum 5 mass%, and the experiments have been done separately in three subdivided temperature ranges of the dispersed PCM particles in a solid phase, two phases (coexistence of solid and liquid) and a liquid phase. The results showed that the Nusselt number increased slightly with the PCM mass concentration for the slurry in solid phase. In the phase change temperature range, the Nusselt number increased with an increase in PCM mass concentration of the slurry at low Rayleigh numbers, while it decreased with increasing PCM mass concentration of the slurry at high Rayleigh numbers. There was not much difference in natural heat transfer characteristics of the PCM slurry with low PCM concentrations (<10 mass%), however, the difference was getting greater with increasing the PCM concentration, especially for the enclosure at a lower aspect ratio (width/height of the rectangular enclosure). The enclosure height was varied from 5.5 to 24.6 mm under a fixed width and depth of 120 mm. Hence, the experiments were performed for a wide range of modified Rayleigh number from 3 × 10² to 1.0 × 10⁷. The correlation generalized for the PCM slurry in a single phase was derived in the form of Nu=0.22(1−C₁C_me^−C₂AR)Ra^1/(3n+1), where C₁ and C₂ were the optimum fitting constants obtained by the least square method. While the PCM was in a phase changing region, the correlation could be expressed as Nu=0.22(1−C₁C_me^−C₂AR)Ra^1/(3n+1)Ste^−0.25, where the Ste was the modified Stefan number.     

Thermal Performance of a Phase Change Material‐Based Latent Heat Thermal Storage Unit

An experimental analysis is presented to establish the thermal performance of a latent heat thermal storage (LHTS) unit. Paraffin is used as the phase change material (PCM) on the shell side of the shell and tube‐type LHTS unit while water is used as the heat transfer fluid (HTF) flowing through the inner tube. The fluid inlet temperature and the mass flow rate of HTF are varied and the temperature distribution of paraffin in the shell side is measured along the radial and axial direction during melting and solidification process. The total melting time is established for different mass flow rates and fluid inlet temperature of HTF. The motion of the solid–liquid interface of the PCM with time along axial and radial direction of the test unit is critically evaluated. The experimental results indicate that the melting front moves from top to bottom along the axial direction while the solidification front moves only in the radial direction. The total melting time of PCM increases as the mass flow rate and inlet temperature of HTF decreases. A correlation is proposed for the dimensionless melting time in terms of Reynolds number and Stefan number of HTF. © 2013 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com/journal/htj ). DOI 10.1002/htj.21120     

Heat transfer by natural convection in a vertical porous layer

Heat transfer by a natural convection in a vertical porous layer heated from below and cooled from above is studied analytically. In the case of linear theory, the normal mode technique is used to find the criteria for the onset of convection and it is shown that convection sets in when the cirtical Reyleigh number exceeds π². The nonlinear theory is investigated using normal mode technique combined with the orthonormal sequences which determines the amplitudes and hence the heat transfer. It is shown that uni-cellular pattern exist and the corresponding heat transfer increases with Rayleigh number.     

Exergy analysis of two phase change materials storage system for solar thermal power with finite-time thermodynamics

A mathematical model for the overall exergetic efficiency of two phase change materials named PCM1 and PCM2 storage system with a concentrating collector for solar thermal power based on finite-time thermodynamics is developed. The model takes into consideration the effects of melting temperatures and number of heat transfer unit of PCM1 and PCM2 on the overall exergetic efficiency. The analysis is based on a lumped model for the PCMs which assumes that a PCM is a thermal reservoir with a constant temperature of its melting point and a distributed model for the air which assumes that the temperature of the air varies in its flow path. The results show that the overall exergetic efficiency can be improved by 19.0-53.8% using two PCMs compared with a single PCM. It is found that melting temperatures of PCM1 and PCM2 have different influences on the overall exergetic efficiency, and the overall exergetic efficiency decreases with increasing the melting temperature of PCM1, increases with increasing the melting temperature of PCM2. It is also found that for PCM1, increasing its number of heat transfer unit can increase the overall exergetic efficiency, however, for PCM2, only when the melting temperature of PCM1 is less than 1150 K and the melting temperature of PCM2 is more than 750 K, increasing the number of heat transfer unit of PCM2 can increase the overall exergetic efficiency. Considering actual application of solar thermal power, we suggest that the optimum melting temperature range of PCM1 is 1000-1150 K and that of PCM2 is 750-900 K. The present analysis provides theoretical guidance for applications of two PCMs storage system for solar thermal power.     

Lattice Boltzmann simulation of tree-shaped fins enhanced melting heat transfer

Abstract

A thermal lattice Boltzmann model is developed to simulate the melting process with natural convection in a cavity filled with tree-shaped solid fins, in which the velocity field and temperature field distribution functions are considered. The present model incorporates the total enthalpy and a free parameter in the equilibrium distribution function to handle conjugate heat transfer. The results indicate that natural convection of liquid phase change material (PCM) plays a significant role in the melting heat transfer of PCM. Increasing the number of branching levels leads to a more rapid melting process, and selecting appropriate bifurcation angle has more efficient heat transfer performance.     

Numerical Studies on Natural Convection in a Trapezoidal Enclosure With Discrete Heating

Abstract

A steady state laminar natural convection flow in a trapezoidal enclosure with discretely heated bottom wall, adiabatic top wall, and constant temperature cold inclined walls is performed. The finite volume based commercial code “ANSYS-FLUENT” is used to investigate the influence of discrete heating on natural convection flows in a trapezoidal cavity. The numerical solution of the problem covers various Rayleigh numbers ranging from 10³ to 10⁶, non-dimensional heating length ranging from 0.2 to 0.8 and Prandtl number is 0.7. The performance of the present numerical approach is represented in the form of streamfunction, temperature profile and Nusselt number. Heat transfer increases with increase of Rayleigh numbers at the corners of the cavity for same heating length from center of the bottom wall. However, the heat transfer rate is less and almost constant for the Rayleigh numbers considered. It is found that the average Nusselt number monotonically increases with increase of Rayleigh number and length of heat source. The variation of local and average Nusselt numbers is more significant for larger length of heating than smaller one. The heat transfer correlations useful for practical design problems have been predicted.     

Novel ventilation cooling system for reducing air conditioning in buildings.: Part I: testing and theoretical modelling

A latent heat storage unit incorporating heat pipes embedded in phase change material (PCM) is developed and tested for a novel application in low energy cooling of buildings. A one-dimensional mathematical model of the heat transfer from air to PCM is presented to allow sizing of a test unit. Details of the construction and testing of one heat pipe/PCM unit in a controlled environment are described, and measurements of heat transfer rate and melting times are presented. When the difference between air and PCM temperature was 5°C, the heat transfer rate was approximately 40 W over a melt period of 19 h. The heat transfer rate could be improved, and the phase change time shortened, with an alternative design for finning of the heat pipe inside the PCM.     

Numerical simulations on performance enhancement of a cross-flow latent thermal energy storage heat exchanger

基于列管式换热器具有传热面积大、结构紧凑、操作弹性大等优点,使其在相变储能领域具有广阔的应用前景。本文建立一种新型列管式相变蓄热器模型,在不考虑自然对流的情况下,利用Fluent软件对相变蓄热器进行二维储热过程的数值模拟。本文主要研究斯蒂芬数、雷诺数、列管排列方式、肋片数以及相变材料的导热系数对熔化过程的影响,并对熔化过程中固液分界面的移动规律进行了分析。模拟结果表明,内肋片强化换热效果明显,特别是对应用低导热系数相变材料[导热系数小于1 W/(m·K)]的列管式蓄热器,相对于无肋片结构,加入肋片（N_fn=2）可缩短熔化时间52.6%。