# Numerical Simulation of the Heat Transfer Behavior of a Zigzag Plate Containing a Phase Change Material for Combustion Heat Recovery and Power Generation

Authors: Peilun Wang; Pengxiang Song; Yun Huang; Zhijian Peng; Yulong Ding
Journal: Journal of Combustion
 (2016)

---

## Abstract

This study presents a numerical analysis of the melting process of phase change materials (PCMs) within a latent heat thermal energy storage (LHTES) system employing zigzag plate. The numerical model used NaCl-MgCl2 mixture as PCMs and hot air as heat transfer fluid (HTF). An experimental system was built to validate the model, and the experimental data agrees reasonably well with the simulation results. The simulation results revealed the effects of the Reynolds and Stefan numbers and the surface topography of the zigzag plate on the charging process. Besides, the effect of the relationship between Reynolds and Stefan numbers on the charging process under a new boundary condition employing a fixed input power was studied. It is found that by modifying the shape of the zigzag plate surface it is feasible to enhance the heat transfer of the LHTES unit remarkably. The melting rate of PCMs increases with the value of Ste or Re numbers with only one of them changing; however, the melting rate of PCMs decreases with the increasing Ste (or decreasing Re) in a fixed input power condition.


---

## Body

## 1. Introduction

With the global development of heavy industrial sectors, an increasing consumption of fossil fuel which resulted in the risk of environmental degradation is becoming increasingly serious every day. The global energy consumption has been growing more rapidly than the economy over the past decades and in 2010 its growth reached 5.6% leading to a highest record since 1973 [1]. In 2012, China’s energy consumption grew by 11.2% overtaking the US as the world’s largest energy consumer [2]. But world primary energy consumption grew by 1.8% in 2012, well below the 10-year average of 2.6%, which benefited from the utilization of renewable energy especially the solar energy and the energy conversation. However, solar energy is susceptible thus irregular in nature and unavailable during cloudy time and nighttime. As a result, to collect and store the solar energy in the daytime for continual supply or evening use is critical. In China, the electricity generated by coal-fired power plant accounts for 77% of the total quantity of power industry. Therefore, in perspective of the energy utilization, it is important to make full use of the low-grade heat from combustion process and improve the total thermal efficiency of the whole plant. Because of its simple structure and high thermal efficiency to run in terms of the recovery of low temperature waste heat, phase change materials (PCMs) thermal heat storage technology is an ideal way to build a bridge between the energy consumption and utilization. Thermal energy can be stored in the form of latent or sensible heat or both. Compared with sensible heat thermal energy storage (SHTES) material, PCMs is particularly attractive since it can provide high storage capacity within a relatively small temperature range. Thus, the latent heat thermal energy storage (LHTES) systems have better prospect to fulfill the needs of the energy hungry society and improve the performance and reliability of the energy system [3], which is more attractive in the field of solar energy collecting [4], industrial waste heat recycling [5], and peak electricity shifting [6]. However, because of the low thermal conductivity of PCMs (usually inorganic salts or organic materials), the performance of the LHTES systems is drastically affected during energy retrieval or withdrawal, which leads to low charging and discharging rates of the LHTES systems [7].In order to overcome the limitations of the LHTES units, a large number of heat transfer enhancement techniques have been approached from different aspects. Fukai et al. [8], Mettawee and Assassa [9], and Sarı and Karaipekli [10] applied, respectively, carbon-fiber, aluminum powder and expanded graphite into the PCMs to improve the thermal conductivity of the PCMs themselves. Another key aspect is to optimize the performance of the heat exchanger, such as introducing fins into the shell-and-tube exchanger [11–13] and heat sink exchanger [14] and utilizing coil pipe [4]. In addition, other designs for enhancing the heat transfer, such as spinning tube [15], microchannel [16], and multiple tube [17], have also been studied.Because of the compact structure and large heat transfer area, the plate heat exchanger filled with PCMs is an ideal design. However, most studies on the LHTES unit employed the heat exchangers of shell-and-tube type or tube-tank type, thus the plate type LHTES unit was rarely studied. Besides, to the authors’ best knowledge, no studies have been found in the literature about the impact of the zigzag surface topography on the heat transfer of plate type LHTES unit. Therefore, this work concentrates on the impact of zigzag configuration on the heat transfer of a plate type LHTES unit. A new case of a fixed input power with various combinations of inlet velocity and inlet temperature during the charging process was also studied. In order to mitigate the heavy demand on computing resources, the heat exchange unit employing zigzag configuration was assumed to be a two-dimensional model and the calculation was focused on the melting process. Moreover, an experimental system was designed accordingly to validate the model.


## 2. Numerical Modeling of the Two-Dimensional Melting Problem

### 2.1. Physical Model

The configuration of the smooth plate (a) and zigzag plate (b) is shown in Figure1. The HTF flows between the two PCMs plates, leading to the melting or solidification of the PCMs upon absorbing/releasing heat from/to the HTF. For the smooth plate (a), the length (L) and thickness (2hl) of the computation domain are 360 and 3.75 mm; the geometries of the smooth plate are given as h1=h2=4 mm; and the thickness of the plate wall is 1 mm, as shown in Figure 1(a). For the zigzag plate, the geometries of the zigzag plate are given as H1=H2=2 mm; H3=H4=4 mm; H5=2 mm and the thickness of the plate wall is 1 mm, as shown in Figure 1(b). Since the configuration of the heat exchanger is periodic, only one segment was employed for simulation (see Figure 1). Each LHTES unit, either smooth or zigzag plates type, contains the same amount of PCMs which is salt mixture composed of sodium chloride and magnesium chloride (NaCl/MgCl2) with mass ratio 4 : 6. The physical properties of the PCMs and the HTF are tabulated in Table 1.Table 1
Thermophysical properties of the NaCl-MgCl2 salt in the validation experiment.
Name Value Density (kg/m3) 2200 Thermal conductivity (W/(m·K)) 0.6 Specific heat (J/(kg·K)) 930 Melting point (K) 713.15 Latent heat (J/kg) 280000Figure 1
The schematic of the LHTES unit.
(a)
Smooth plate (b)
Zigzag plate


### 2.2. Mathematical Model

The continuity equation is as follows:(1)∂ρ∂t+∇·ρu→=0,where ρ is density and u→ is the fluid velocity.The momentum equation is as follows:(2)∂ρu∂t+∇ρu→μ=∇μ∇·u-∂P∂x+Sx,∂ρv∂t+∇ρu→v=∇μ∇·v-∂P∂y+Sy+Sb,where u and v are the velocity inx andy directions; μ is viscosity; P is pressure;t is time; Sx and Sy are the momentum sources in the enthalpy-porosity formulation in the form of Si=Amushui1-β2/β3+ε in which Amush is mushy zone constant that measures the amplitude of the damping and ε is a small number (e.g. 0.0001) to prevent division by zero; and Sb=ρgβ(h-href)/Cp is the buoyancy source term.The conservation equation of energy is as follows:(3)∂∂tρh+∇·ρu→h=∇·α∇T+SE,where SE=∂ρΔH/∂t+∇ρu→ΔH is a phase related source term and the enthalpy of the material is the sum of the sensible enthalpy, h, and the latent heat, ΔH, in (3):(4)H=href+∫TrefTcpdT+βL,where href is the enthalpy value at the reference temperature Tref and β is the liquid fraction which is defined as follows:(5)β=1T≥TlT-TsTl-TsTl>T>Ts0T≤Ts.In (2), μ is viscosity; Sb=ρgβ(h-href)/Cp is the buoyancy source term; Sx and Sy are the momentum sources in the enthalpy-porosity formulation in the form of Si=Amushui1-β2/β3+ε in which Amush is mushy zone constant that measures the amplitude of the damping and ε is a small number (e.g. 0.0001) to prevent division by zero. In (3), u→ is the fluid velocity, α is the thermal conductivity, μ and ρ are viscosity and density, respectively, and SE=∂ρΔH/∂t+∇ρu→ΔH is a phase related source term.In order to reduce the computational cost of the simulations, the geometry of the plate is considered to be sufficiently wide and the 2-dimensional model is used to simulate the thermal behavior of the PCMs-plate. Besides, the geometry of the PCMs-plate is periodic and only one channel surrounded by two halves of the PCMs plates is simulated. The computational domain is shown in Figure1.


### 2.3. Boundary and Initial Conditions

The interface between the PCMs and the inner tube was coupled-wall boundary condition for heat transfer. The inlet was assumed as velocity inlet with uniform velocity and temperature (which was higher than the melting point of the PCMs). The outlet was set as outflow boundary. The upper boundary and lower boundary of the simulation domain are assumed to be periodic. The front and back walls are considered as adiabatic boundary conditions.The entire calculation domain was set initially at a zero velocity and a temperature ofTo (which was lower than the melting point of the PCMs).


## 2.1. Physical Model

The configuration of the smooth plate (a) and zigzag plate (b) is shown in Figure1. The HTF flows between the two PCMs plates, leading to the melting or solidification of the PCMs upon absorbing/releasing heat from/to the HTF. For the smooth plate (a), the length (L) and thickness (2hl) of the computation domain are 360 and 3.75 mm; the geometries of the smooth plate are given as h1=h2=4 mm; and the thickness of the plate wall is 1 mm, as shown in Figure 1(a). For the zigzag plate, the geometries of the zigzag plate are given as H1=H2=2 mm; H3=H4=4 mm; H5=2 mm and the thickness of the plate wall is 1 mm, as shown in Figure 1(b). Since the configuration of the heat exchanger is periodic, only one segment was employed for simulation (see Figure 1). Each LHTES unit, either smooth or zigzag plates type, contains the same amount of PCMs which is salt mixture composed of sodium chloride and magnesium chloride (NaCl/MgCl2) with mass ratio 4 : 6. The physical properties of the PCMs and the HTF are tabulated in Table 1.Table 1
Thermophysical properties of the NaCl-MgCl2 salt in the validation experiment.
Name Value Density (kg/m3) 2200 Thermal conductivity (W/(m·K)) 0.6 Specific heat (J/(kg·K)) 930 Melting point (K) 713.15 Latent heat (J/kg) 280000Figure 1
The schematic of the LHTES unit.
(a)
Smooth plate (b)
Zigzag plate


## 2.2. Mathematical Model

The continuity equation is as follows:(1)∂ρ∂t+∇·ρu→=0,where ρ is density and u→ is the fluid velocity.The momentum equation is as follows:(2)∂ρu∂t+∇ρu→μ=∇μ∇·u-∂P∂x+Sx,∂ρv∂t+∇ρu→v=∇μ∇·v-∂P∂y+Sy+Sb,where u and v are the velocity inx andy directions; μ is viscosity; P is pressure;t is time; Sx and Sy are the momentum sources in the enthalpy-porosity formulation in the form of Si=Amushui1-β2/β3+ε in which Amush is mushy zone constant that measures the amplitude of the damping and ε is a small number (e.g. 0.0001) to prevent division by zero; and Sb=ρgβ(h-href)/Cp is the buoyancy source term.The conservation equation of energy is as follows:(3)∂∂tρh+∇·ρu→h=∇·α∇T+SE,where SE=∂ρΔH/∂t+∇ρu→ΔH is a phase related source term and the enthalpy of the material is the sum of the sensible enthalpy, h, and the latent heat, ΔH, in (3):(4)H=href+∫TrefTcpdT+βL,where href is the enthalpy value at the reference temperature Tref and β is the liquid fraction which is defined as follows:(5)β=1T≥TlT-TsTl-TsTl>T>Ts0T≤Ts.In (2), μ is viscosity; Sb=ρgβ(h-href)/Cp is the buoyancy source term; Sx and Sy are the momentum sources in the enthalpy-porosity formulation in the form of Si=Amushui1-β2/β3+ε in which Amush is mushy zone constant that measures the amplitude of the damping and ε is a small number (e.g. 0.0001) to prevent division by zero. In (3), u→ is the fluid velocity, α is the thermal conductivity, μ and ρ are viscosity and density, respectively, and SE=∂ρΔH/∂t+∇ρu→ΔH is a phase related source term.In order to reduce the computational cost of the simulations, the geometry of the plate is considered to be sufficiently wide and the 2-dimensional model is used to simulate the thermal behavior of the PCMs-plate. Besides, the geometry of the PCMs-plate is periodic and only one channel surrounded by two halves of the PCMs plates is simulated. The computational domain is shown in Figure1.


## 2.3. Boundary and Initial Conditions

The interface between the PCMs and the inner tube was coupled-wall boundary condition for heat transfer. The inlet was assumed as velocity inlet with uniform velocity and temperature (which was higher than the melting point of the PCMs). The outlet was set as outflow boundary. The upper boundary and lower boundary of the simulation domain are assumed to be periodic. The front and back walls are considered as adiabatic boundary conditions.The entire calculation domain was set initially at a zero velocity and a temperature ofTo (which was lower than the melting point of the PCMs).


## 3. Experimental Validation

### 3.1. Experiment System and Procedure

The numerical model was validated by experimental work. Figure2 shows the schematics of the high temperature heat exchange system, consisting of an electric furnace, a high temperature blower, a LHTES unit, and a number of data acquisition equipment. During the charging process, the air was driven to circulate in the loop by the high temperature blower where it is heated by the electric furnace. The heat was transferred to and stored within the PCMs when the heated air flowed through the latent heat exchanger. The temperature of PCMs and outlet is measured by K-type thermocouples (the measurement range is 273.15 K~1273.15 K and error range is ± 0.2%) and the flow rate of the air is measured by a orifice meter (error range ± 2%). Upon fully charged condition (judged from the temperature of the outlet air), the charging step would be ended by switching off the electric furnace and the high temperature blower. During the discharging process, the ambient air was switched on to deliver an air flow at a preset rate to extract the stored heat from the charged LHTES unit.Figure 2
High temperature heat exchange test system.Temperature of the phase change materials was measured by thermocouples (see Figure3 for the thermocouple arrangement). As shown in Figure 3, five K-type thermocouples connecting to a PC-based temperature data acquisition system inserted into the PCMs to measure the temperature of the PCMs. The exterior of the heat exchanger is shown in Figure 4.Figure 3
Arrangement of thermocouples and structural representation of zigzag plate.Figure 4
Exterior of the heat exchanger.


### 3.2. Validation of the Simulation Models

The computational domain is represented by a 2D base grid as shown in Figure5. The finer grid is created by increasing the mesh density near both sides of the plate walls. In order to check the grid independence, grids of three sizes were checked, that is, 0.3 mm, 0.4 mm, and 0.5 mm with the cell amount of 84932, 63410, and 49163, respectively. No obvious differences in the liquid fraction were found among the three different grids. Therefore, the base grid with grid size of 0.5 mm was chosen in the following simulation, composed of 49163 cells including 19560 cells in the PCMs domain, 13770 cells in HTF domain, and the rest in the zigzag plate.Figure 5
Schematic representation of the zigzag plate on 2D base grid.An initial run was performed and compared with the experimental data to validate the computational modeling of the melting process. In the experiment, the inlet temperature and the inlet velocity were 863.15 K and 3.1 m/s and the NaCl-MgCl2 mixture with a 4 : 6 mass ratio was used, presenting the phase change temperature of 713.15 K. Other thermophysical properties of the PCMs are shown in Table 1. The specific heat and thermal conductivity of HTF (air) within range of 200~1000 K were given as follows.Specific heat (J/(kg·K))(6)Cp=1043.96081-0.327T+7.79328×10-4T2-3.5606×10-7T3.Thermal conductivity (W/(m·K))(7)λ=0.0096+5.975×10-5T.As shown in Figure 6, the simulated value agrees reasonably well with the experimental data. The measured outlet temperature as a function of time is compared with the calculated result and is found to be in good agreement. However, the deviation of thermocouple 1 between the simulation and experiment in Figure 6 could be caused by several issues. One is that the PCMs with pores in the form of powder were used in the experiment, which reduces the conductivity of PCMs and more heat is delivered to the downstream region; other reasons may include that the properties of the PCMs are correlated with the varying temperature and the heat transfer fluid flow is not identical in all different plates, which were not considered in this simulation.Figure 6
Comparison between the measured and the calculated temperature as a function of time.


## 3.1. Experiment System and Procedure

The numerical model was validated by experimental work. Figure2 shows the schematics of the high temperature heat exchange system, consisting of an electric furnace, a high temperature blower, a LHTES unit, and a number of data acquisition equipment. During the charging process, the air was driven to circulate in the loop by the high temperature blower where it is heated by the electric furnace. The heat was transferred to and stored within the PCMs when the heated air flowed through the latent heat exchanger. The temperature of PCMs and outlet is measured by K-type thermocouples (the measurement range is 273.15 K~1273.15 K and error range is ± 0.2%) and the flow rate of the air is measured by a orifice meter (error range ± 2%). Upon fully charged condition (judged from the temperature of the outlet air), the charging step would be ended by switching off the electric furnace and the high temperature blower. During the discharging process, the ambient air was switched on to deliver an air flow at a preset rate to extract the stored heat from the charged LHTES unit.Figure 2
High temperature heat exchange test system.Temperature of the phase change materials was measured by thermocouples (see Figure3 for the thermocouple arrangement). As shown in Figure 3, five K-type thermocouples connecting to a PC-based temperature data acquisition system inserted into the PCMs to measure the temperature of the PCMs. The exterior of the heat exchanger is shown in Figure 4.Figure 3
Arrangement of thermocouples and structural representation of zigzag plate.Figure 4
Exterior of the heat exchanger.


## 3.2. Validation of the Simulation Models

The computational domain is represented by a 2D base grid as shown in Figure5. The finer grid is created by increasing the mesh density near both sides of the plate walls. In order to check the grid independence, grids of three sizes were checked, that is, 0.3 mm, 0.4 mm, and 0.5 mm with the cell amount of 84932, 63410, and 49163, respectively. No obvious differences in the liquid fraction were found among the three different grids. Therefore, the base grid with grid size of 0.5 mm was chosen in the following simulation, composed of 49163 cells including 19560 cells in the PCMs domain, 13770 cells in HTF domain, and the rest in the zigzag plate.Figure 5
Schematic representation of the zigzag plate on 2D base grid.An initial run was performed and compared with the experimental data to validate the computational modeling of the melting process. In the experiment, the inlet temperature and the inlet velocity were 863.15 K and 3.1 m/s and the NaCl-MgCl2 mixture with a 4 : 6 mass ratio was used, presenting the phase change temperature of 713.15 K. Other thermophysical properties of the PCMs are shown in Table 1. The specific heat and thermal conductivity of HTF (air) within range of 200~1000 K were given as follows.Specific heat (J/(kg·K))(6)Cp=1043.96081-0.327T+7.79328×10-4T2-3.5606×10-7T3.Thermal conductivity (W/(m·K))(7)λ=0.0096+5.975×10-5T.As shown in Figure 6, the simulated value agrees reasonably well with the experimental data. The measured outlet temperature as a function of time is compared with the calculated result and is found to be in good agreement. However, the deviation of thermocouple 1 between the simulation and experiment in Figure 6 could be caused by several issues. One is that the PCMs with pores in the form of powder were used in the experiment, which reduces the conductivity of PCMs and more heat is delivered to the downstream region; other reasons may include that the properties of the PCMs are correlated with the varying temperature and the heat transfer fluid flow is not identical in all different plates, which were not considered in this simulation.Figure 6
Comparison between the measured and the calculated temperature as a function of time.


## 4. Results and Discussion

### 4.1. Comparison between the Melting Rate of the Zigzag Plate and the Smooth Plate

Under the same boundary conditions (i.e., inlet temperature is 873.15 K, velocity is 6 m/s) and initial condition (i.e., initial temperature of the PCMs is 673.15 K), the liquid fraction as a function of time was calculated for the LHTES systems employing smooth plate and zigzag plate (see Figure7), respectively. The charging process, during which the heat is absorbed and stored, can be divided into three steps. At the beginning, the heat is stored in the PCMs as sensible heat because the temperature of the PCM is below its melting temperature and the temperature difference between the HTF and the PCMs is large. Next, part of the PCMs starts to melt as its temperature increased and the heat is stored as latent heat during the phase change process, while the readily melted part and the solid part of the PCMs are absorbing and storing heat as sensible heat without phase change. At the end of the charging process, PCMs in the exchangers are melted completely, so the heat is stored in the sensible way.Figure 7
Liquid fraction as a function of time for the LHTES systems employing smooth and zigzag plates.From Figure7 it can be seen that the zigzag surface topography (with the length of 360 mm) can enhance the heat transfer of the LHTES system remarkably, since it costs about 3750 s to complete the whole melting process compared to that of about 4200 s by the system employing smooth plate of the same length with the zigzag plate. This effect can be attributed to two reasons. On one hand, the zigzag shape surface has the ability to change the flow form from laminar to turbulent by disturbing the flow direction, which can decrease the boundary layer thickness resulting in the decrease of thermal resistance at the HTF plate interface; on the other hand, the zigzag configuration presents larger heat transfer area than the smooth configuration, enhancing the heat transfer rate during the charging and discharging process. However, the smooth plate, which has the same heat transfer area with the zigzag plate, costs about 3850 s to complete the whole melting process. The zigzag surface topography does not have much effect on the heat transfer enhancement if the zigzag plate has the same heat transfer area with the smooth plate. But the zigzag topography has the ability of shortening the plate length.


### 4.2. The Effect of Tooth Shape on the Charging Process

In this section, zigzag plates of four distinct zigzag lengths, that is,M=12, 18, 24, and 30 mm, are studied and the simulation results of the liquid fraction as a function of time is shown in Figure 8. All cases are simulated under the condition of vin=6 m/s and Tin=873.15 K with the same amount of PCMs in the zigzag plates. It can be seen from Figure 8 that the four liquid fraction curves are almost coincident with each other over the entire charging process, implying that the zigzag length almost has no influence on the heat transfer capability of the LHTES system. It is well known that the heat transfer area and the flow distribution have great impact on the heat transfer capability of plate heat exchanger. With the increase of the tooth length, the heat transfer area decreases, but the adherence character of the fluid strengthens. As shown in Figure 9, the mainstream region is near the middle part of the passageway in the plates of small tooth length, but as the tooth length gets bigger, the mainstream region gets closer to the plate wall, where in the latter case the heat transfer is enhanced. Thus, the joint effect from the above two aspects leads to the tooth length almost having no impact on the melting process.Figure 8
Time variation of the liquid fraction at various zigzag lengths.Figure 9
Velocity vector for different tooth lengths.Furthermore, by changing the individual plate lengths, the impact of tooth rates on the heat transfer capability of the zigzag plate is explored. The parameter, tooth rate, is defined asU=D1/D2 along the flow direction as shown in Figure 10. Figure 11 shows that the temporal variation of liquid fraction for different tooth rates are almost coincided with each other under the condition of vin=6 m/s and Tin=923.15 K. Thus, the tooth rate has little influence on the heat transfer capability, where two aspects should be considered. On one hand, the total tooth length (D1+D2) of the five plate heat exchangers with different tooth rates are equivalent but their individual heat transfer areas are different, which are 760.3 mm, 747.984 mm, 745.032 mm, 747.984 mm, and 760.3 mm for the tooth rates of 1/5, 1/2, 1/1, 2/1, and 5/1, respectively. However, this difference is too small to produce noticeable impact on the heat transfer process. On the other hand, although the flow fields within the different plates are different with each other, they caused undistinguishable effect on the heat transfer capability of the heat exchangers. It is well known that the near-wall flow field formed near the upwind wall can enhance heat transfer more effectively than that formed near the lee side wall. The near-wall regions of the five cases showed distinct characteristic. When U<1/1, all the near-wall regions are formed near the upwind wall at the two sides of the HTF passageway. However, when U⩾1/1, the near-wall regions are formed, respectively, on the upwind wall and the lee side wall at the same side of the HTF passageway, as shown in Figure 12, which is the schematic plot of flow field near the wall. Figure 13 shows the velocity vector for different tooth ratio and, as shown in Table 2 and Figure 13, there is nearly no difference in lUlU=lU-a+lU-b between the first two cases with U<1/1, but lU of the later three cases with U⩾1/1 increases with the increase ofU. Besides, for the later three cases, lU-a (near lee side wall) increases with increasingU but lU-b (near upwind wall) decreases with increasingU. Therefore, the length of near-wall regions (lU) of the first two cases is smaller than that of the later three cases, but the near-wall regions of the first two cases locate close to the upwind wall while the two near-wall regions of the later three cases locate close to, respectively, the upwind wall and the lee side wall. Therefore, by considering the compromise between the length and the location of the near-wall regions, the change of flow filed caused by tooth rate almost has no impact on the heat exchange capability of the plates.Table 2
The values ofLU-a, LU-b, and (LU-a) + (LU-b) at different U.
U L U - a  [mm] L U - b [mm] (LU-a) + (LU-b)[mm] U = 1 / 5 5.97 6.33 12.30 U = 1 / 2 7.99 5.38 13.37 U = 1 / 1 7.71 8.37 16.08 U = 2 / 1 13.00 6.29 19.29 U = 5 / 1 17.17 1.74 18.91Figure 10
Schematic illustration of the shape of a tooth.Figure 11
The liquid fraction as a function of time for the plates of different tooth rates.Figure 12
Schematic plot of flow field near the wall.
(a)
U ≤ 1 : 1 (b)
U > 1 : 1Figure 13
Velocity vector for different tooth ratio.The heat stored by PCMs can be separated into sensible heat and latent heat due to the temperature rising and the reversible phase change between liquid and solid. In order to evaluate the ratio of latent heat (QL) in total stored heat (Q), the value of QL/Q is studied in this section and the sensible heat (QS), latent heat, and total heat can be expressed as follows:(8)QS=Cp·m·ΔT,QL=L·m·β,Q=QS+QL,where Q is the total heat stored in the PCMs and β is the liquid fraction of the PCMs.The pressure drop can be calculated as follows:(9)ΔP=∫AoPdA-∫AinPdA,where Ao and Ain are the area of the outlet and inlet of computational domain, respectively.Figure14 shows the ratio of the latent heat (QL) in the total heat (Q, composed of latent heat and sensible heat) for various zigzag lengths. As shown in Figure 11, the value of QL/Q rises quickly and gets to the first peak value at c.a. 500 s followed by a small decrease before c.a. 1000 s. This first QL/Q peak can be attributed to the large temperature difference between the HTF and the PCMs at the beginning, resulting in the melting of the PCMs at the inlet region which rapidly absorbs and stores the heat of HTF as latent heat. After certain time, once the PCMs at the inlet region are completely melted, the PCMs continue absorbing heat and store it in sensible form leading to the decrease of QL/Q. As the above process proceeds, the temperature difference between the PCMs and the HTF at the inlet region decreases and more heat transforms to the downstream region. As a result, the value of QL/Q increases with the melting of the PCMs in the downstream region, causing the second peak in QL/Q. Finally, since most of the PCMs are melted, more heat is continuously absorbed and stored in sensible form and the value of QL/Q goes down gradually.Figure 14
The proportion ofQL in Q as a function of time among the heat exchangers with different tooth lengths.In order to evaluate the heat exchange unit comprehensively, the pressure drop is also investigated for different tooth length and tooth rate plates. Figure15(a) shows the effects of zigzag tooth length on the pressure drop. The pressure drop decreases dramatically from 98 Pa to 49 Pa while the tooth length increases from 12 mm to 39 mm. Nevertheless, from Figure 15(b) it can be seen that the pressure drop decreases from 51 Pa to 32 Pa while the tooth rate increases from 1/5 to 5/1. Therefore, the increase of tooth length or tooth rate can reduce the pressure drop effectively, enhancing the comprehensive performance of the heat exchange unit.Figure 15
Pressure drop calculated for the heat exchanger unit.
(a)
Pressure drop as a function of tooth length at the same tooth rate 1 : 1 (b)
Pressure drop as a function of tooth rate at the same tooth lengthm = 30 mm


### 4.3. Effect of Different Working Conditions

In order to investigate the effect of inlet velocity on the melting process, the temporal variation of the liquid fraction at different Re numbers with the same Ste number was calculated. The Reynolds number is generally defined as(10)Re=ρvDHμ,where DH is the hydraulic diameter defined as DH=2wbc/w+bc≈2bc, while w is the plate width inside the gasket and bc is the mean spacing between plates. ρ, μ, and ν are the density, viscosity, and velocity of the HTF, respectively.Figure16 presents the liquid fraction of the PCMs as a function of time obtained by employing different Re numbers. Obviously the slope of the curves is increasing while the Re numbers are increasing from 98 to 392. This results in a shortened melting time for the relatively larger Re numbers. This variation trend is attributed to two reasons. Firstly, higher Re corresponds to faster flow velocity which prompts the heat transfer. Besides, increasing Re number is accompanied with increasingly strong turbulence, which results in the HTF temperature being more uniform in each zigzag tooth. Secondly, the inlet heat energy and the temperature of HTF in the passageway also increase with the increase of Re number, enhancing the heat transfer.Figure 16
The liquid fraction of the PCMs as a function of time during the charging process at different Re number.Figure17 presents the liquid fraction of the PCMs as a function of time obtained by employing different Ste number with the same inlet velocity. The Ste number expresses the ratio of sensible heat due to the super-heating arisen from the hot air to latent heat of heating of PCMs and the Ste number is generally defined as(11)Ste=cpTin-TmL.As definition of Ste in (11), higher Ste brings about a greater temperature difference between the phase-changing region and the HTF. Thus, it is expected that the melting rate of PCMs increases with the increase of Ste number. This speculation is verified effectively through the results shown in Figure 17. A low Ste value means a low inlet temperature. This also means a small temperature difference between the HTF and the PCMs, and thus the heat transfer is lower at the smaller Ste values. Besides, Figure 17 also indicates the change trend of the liquid fraction with increasing Ste number, which implies that the heat transfer enhancement decreases with the Ste number.Figure 17
The liquid fraction of the PCMs as a function of time during the charging process at different Ste number.Figure18 represents temporal variation of liquid fraction and outlet temperature under different conditions of the given fixed input power.Figure 18
The liquid fraction and outlet temperature under different Ste or Re numbers conditions for fixed input power.The input power can be expressed as(12)W=cp·ρvAa·Tin-Tb,Tb=12Tin+Tini.Tb is the reference temperature and the conditions of vin = 6 m/s, Tin = 873.15 K, and Tini = 663.15 K are chosen as the base conditions. Five combinations of inlet condition at an identical input power are chosen, as shown in Table 3, and the impact of Ste and Re numbers on the energy storage is studied.Table 3
The working conditions at the same inlet power.
Conditions Ste Re 1 0.737 251 2 0.823 204 3 0.881 176 4 0.957 152 5 1.06 130As shown in Figure18, the variation trend is quite similar among the five outlet temperature curves, but overall the outlet temperature rises more quickly with small Ste number (i.e., larger Re number). The outlet temperature (To) remains in low constant value in the early stage which is equal to the Tini and then rises slowly at moderate rates. When the liquid fraction approaches around 0.4, the To reaches a constant value which is close to the phase change temperature of the PCMs, but it can be seen that the onset of the constant temperature period becomes more delayed with decreasing Re numbers. At the end of the charging process, the outlet temperature rises quickly. The tendency of the outlet temperature can be explained as follows. At the beginning, the temperature of the PCMs is low and the heat is absorbed and stored in sensible form, leading to the constant To equaling Tini. At the phase change temperature of the PCMs, the material begins to melt while the heat is absorbed and stored in sensible and latent form, and due to the lower temperature of PCMs in the downstream region and the buffer effect of the melting process, Tout rises slowly but is still lower than the melting point of PCMs at this stage. In the next stage, the temperature of PCMs in the downstream region rises and To is controlled by the phase change temperature of PCMs, which leads to To equaling the melting temperature of PCMs. Because most of PCMs melt and the heat is stored mainly in sensible form in the PCMs, To rises quickly at the end of the charging process.Furthermore, the reaching time of constantTo is delayed for smaller Re numbers. This is because the temperature difference between the HTF and the PCMs increases with decreasing Re numbers, which made more thermal energy get stored in the sensible way in the inlet region, the melting speed is lower in the downstream area and the reaching time of constant value is delayed. But the time duration of constant value decreases with increasing Ste number (decreasing Re number). Also, Figure 18 shows the liquid fraction as a function of time under different conditions for fixed input power. One can see that the overall melting time decreases with increasing Re number (decreasing Ste number), but the increase tendency slows down. This is because, for a given input power, an increase in the Ste number results in an increase in the inlet temperature and a decrease in the inlet velocity. The difference between inlet temperature and initial temperature increases with decreasing Re number according to (12). This implied that more heat is stored in the inlet region as sensible heat and less heat flows to the downstream region to melt the PCMs there, which leads to a lower melting rate in the condition of smaller Re for a given inlet power.


## 4.1. Comparison between the Melting Rate of the Zigzag Plate and the Smooth Plate

Under the same boundary conditions (i.e., inlet temperature is 873.15 K, velocity is 6 m/s) and initial condition (i.e., initial temperature of the PCMs is 673.15 K), the liquid fraction as a function of time was calculated for the LHTES systems employing smooth plate and zigzag plate (see Figure7), respectively. The charging process, during which the heat is absorbed and stored, can be divided into three steps. At the beginning, the heat is stored in the PCMs as sensible heat because the temperature of the PCM is below its melting temperature and the temperature difference between the HTF and the PCMs is large. Next, part of the PCMs starts to melt as its temperature increased and the heat is stored as latent heat during the phase change process, while the readily melted part and the solid part of the PCMs are absorbing and storing heat as sensible heat without phase change. At the end of the charging process, PCMs in the exchangers are melted completely, so the heat is stored in the sensible way.Figure 7
Liquid fraction as a function of time for the LHTES systems employing smooth and zigzag plates.From Figure7 it can be seen that the zigzag surface topography (with the length of 360 mm) can enhance the heat transfer of the LHTES system remarkably, since it costs about 3750 s to complete the whole melting process compared to that of about 4200 s by the system employing smooth plate of the same length with the zigzag plate. This effect can be attributed to two reasons. On one hand, the zigzag shape surface has the ability to change the flow form from laminar to turbulent by disturbing the flow direction, which can decrease the boundary layer thickness resulting in the decrease of thermal resistance at the HTF plate interface; on the other hand, the zigzag configuration presents larger heat transfer area than the smooth configuration, enhancing the heat transfer rate during the charging and discharging process. However, the smooth plate, which has the same heat transfer area with the zigzag plate, costs about 3850 s to complete the whole melting process. The zigzag surface topography does not have much effect on the heat transfer enhancement if the zigzag plate has the same heat transfer area with the smooth plate. But the zigzag topography has the ability of shortening the plate length.


## 4.2. The Effect of Tooth Shape on the Charging Process

In this section, zigzag plates of four distinct zigzag lengths, that is,M=12, 18, 24, and 30 mm, are studied and the simulation results of the liquid fraction as a function of time is shown in Figure 8. All cases are simulated under the condition of vin=6 m/s and Tin=873.15 K with the same amount of PCMs in the zigzag plates. It can be seen from Figure 8 that the four liquid fraction curves are almost coincident with each other over the entire charging process, implying that the zigzag length almost has no influence on the heat transfer capability of the LHTES system. It is well known that the heat transfer area and the flow distribution have great impact on the heat transfer capability of plate heat exchanger. With the increase of the tooth length, the heat transfer area decreases, but the adherence character of the fluid strengthens. As shown in Figure 9, the mainstream region is near the middle part of the passageway in the plates of small tooth length, but as the tooth length gets bigger, the mainstream region gets closer to the plate wall, where in the latter case the heat transfer is enhanced. Thus, the joint effect from the above two aspects leads to the tooth length almost having no impact on the melting process.Figure 8
Time variation of the liquid fraction at various zigzag lengths.Figure 9
Velocity vector for different tooth lengths.Furthermore, by changing the individual plate lengths, the impact of tooth rates on the heat transfer capability of the zigzag plate is explored. The parameter, tooth rate, is defined asU=D1/D2 along the flow direction as shown in Figure 10. Figure 11 shows that the temporal variation of liquid fraction for different tooth rates are almost coincided with each other under the condition of vin=6 m/s and Tin=923.15 K. Thus, the tooth rate has little influence on the heat transfer capability, where two aspects should be considered. On one hand, the total tooth length (D1+D2) of the five plate heat exchangers with different tooth rates are equivalent but their individual heat transfer areas are different, which are 760.3 mm, 747.984 mm, 745.032 mm, 747.984 mm, and 760.3 mm for the tooth rates of 1/5, 1/2, 1/1, 2/1, and 5/1, respectively. However, this difference is too small to produce noticeable impact on the heat transfer process. On the other hand, although the flow fields within the different plates are different with each other, they caused undistinguishable effect on the heat transfer capability of the heat exchangers. It is well known that the near-wall flow field formed near the upwind wall can enhance heat transfer more effectively than that formed near the lee side wall. The near-wall regions of the five cases showed distinct characteristic. When U<1/1, all the near-wall regions are formed near the upwind wall at the two sides of the HTF passageway. However, when U⩾1/1, the near-wall regions are formed, respectively, on the upwind wall and the lee side wall at the same side of the HTF passageway, as shown in Figure 12, which is the schematic plot of flow field near the wall. Figure 13 shows the velocity vector for different tooth ratio and, as shown in Table 2 and Figure 13, there is nearly no difference in lUlU=lU-a+lU-b between the first two cases with U<1/1, but lU of the later three cases with U⩾1/1 increases with the increase ofU. Besides, for the later three cases, lU-a (near lee side wall) increases with increasingU but lU-b (near upwind wall) decreases with increasingU. Therefore, the length of near-wall regions (lU) of the first two cases is smaller than that of the later three cases, but the near-wall regions of the first two cases locate close to the upwind wall while the two near-wall regions of the later three cases locate close to, respectively, the upwind wall and the lee side wall. Therefore, by considering the compromise between the length and the location of the near-wall regions, the change of flow filed caused by tooth rate almost has no impact on the heat exchange capability of the plates.Table 2
The values ofLU-a, LU-b, and (LU-a) + (LU-b) at different U.
U L U - a  [mm] L U - b [mm] (LU-a) + (LU-b)[mm] U = 1 / 5 5.97 6.33 12.30 U = 1 / 2 7.99 5.38 13.37 U = 1 / 1 7.71 8.37 16.08 U = 2 / 1 13.00 6.29 19.29 U = 5 / 1 17.17 1.74 18.91Figure 10
Schematic illustration of the shape of a tooth.Figure 11
The liquid fraction as a function of time for the plates of different tooth rates.Figure 12
Schematic plot of flow field near the wall.
(a)
U ≤ 1 : 1 (b)
U > 1 : 1Figure 13
Velocity vector for different tooth ratio.The heat stored by PCMs can be separated into sensible heat and latent heat due to the temperature rising and the reversible phase change between liquid and solid. In order to evaluate the ratio of latent heat (QL) in total stored heat (Q), the value of QL/Q is studied in this section and the sensible heat (QS), latent heat, and total heat can be expressed as follows:(8)QS=Cp·m·ΔT,QL=L·m·β,Q=QS+QL,where Q is the total heat stored in the PCMs and β is the liquid fraction of the PCMs.The pressure drop can be calculated as follows:(9)ΔP=∫AoPdA-∫AinPdA,where Ao and Ain are the area of the outlet and inlet of computational domain, respectively.Figure14 shows the ratio of the latent heat (QL) in the total heat (Q, composed of latent heat and sensible heat) for various zigzag lengths. As shown in Figure 11, the value of QL/Q rises quickly and gets to the first peak value at c.a. 500 s followed by a small decrease before c.a. 1000 s. This first QL/Q peak can be attributed to the large temperature difference between the HTF and the PCMs at the beginning, resulting in the melting of the PCMs at the inlet region which rapidly absorbs and stores the heat of HTF as latent heat. After certain time, once the PCMs at the inlet region are completely melted, the PCMs continue absorbing heat and store it in sensible form leading to the decrease of QL/Q. As the above process proceeds, the temperature difference between the PCMs and the HTF at the inlet region decreases and more heat transforms to the downstream region. As a result, the value of QL/Q increases with the melting of the PCMs in the downstream region, causing the second peak in QL/Q. Finally, since most of the PCMs are melted, more heat is continuously absorbed and stored in sensible form and the value of QL/Q goes down gradually.Figure 14
The proportion ofQL in Q as a function of time among the heat exchangers with different tooth lengths.In order to evaluate the heat exchange unit comprehensively, the pressure drop is also investigated for different tooth length and tooth rate plates. Figure15(a) shows the effects of zigzag tooth length on the pressure drop. The pressure drop decreases dramatically from 98 Pa to 49 Pa while the tooth length increases from 12 mm to 39 mm. Nevertheless, from Figure 15(b) it can be seen that the pressure drop decreases from 51 Pa to 32 Pa while the tooth rate increases from 1/5 to 5/1. Therefore, the increase of tooth length or tooth rate can reduce the pressure drop effectively, enhancing the comprehensive performance of the heat exchange unit.Figure 15
Pressure drop calculated for the heat exchanger unit.
(a)
Pressure drop as a function of tooth length at the same tooth rate 1 : 1 (b)
Pressure drop as a function of tooth rate at the same tooth lengthm = 30 mm


## 4.3. Effect of Different Working Conditions

In order to investigate the effect of inlet velocity on the melting process, the temporal variation of the liquid fraction at different Re numbers with the same Ste number was calculated. The Reynolds number is generally defined as(10)Re=ρvDHμ,where DH is the hydraulic diameter defined as DH=2wbc/w+bc≈2bc, while w is the plate width inside the gasket and bc is the mean spacing between plates. ρ, μ, and ν are the density, viscosity, and velocity of the HTF, respectively.Figure16 presents the liquid fraction of the PCMs as a function of time obtained by employing different Re numbers. Obviously the slope of the curves is increasing while the Re numbers are increasing from 98 to 392. This results in a shortened melting time for the relatively larger Re numbers. This variation trend is attributed to two reasons. Firstly, higher Re corresponds to faster flow velocity which prompts the heat transfer. Besides, increasing Re number is accompanied with increasingly strong turbulence, which results in the HTF temperature being more uniform in each zigzag tooth. Secondly, the inlet heat energy and the temperature of HTF in the passageway also increase with the increase of Re number, enhancing the heat transfer.Figure 16
The liquid fraction of the PCMs as a function of time during the charging process at different Re number.Figure17 presents the liquid fraction of the PCMs as a function of time obtained by employing different Ste number with the same inlet velocity. The Ste number expresses the ratio of sensible heat due to the super-heating arisen from the hot air to latent heat of heating of PCMs and the Ste number is generally defined as(11)Ste=cpTin-TmL.As definition of Ste in (11), higher Ste brings about a greater temperature difference between the phase-changing region and the HTF. Thus, it is expected that the melting rate of PCMs increases with the increase of Ste number. This speculation is verified effectively through the results shown in Figure 17. A low Ste value means a low inlet temperature. This also means a small temperature difference between the HTF and the PCMs, and thus the heat transfer is lower at the smaller Ste values. Besides, Figure 17 also indicates the change trend of the liquid fraction with increasing Ste number, which implies that the heat transfer enhancement decreases with the Ste number.Figure 17
The liquid fraction of the PCMs as a function of time during the charging process at different Ste number.Figure18 represents temporal variation of liquid fraction and outlet temperature under different conditions of the given fixed input power.Figure 18
The liquid fraction and outlet temperature under different Ste or Re numbers conditions for fixed input power.The input power can be expressed as(12)W=cp·ρvAa·Tin-Tb,Tb=12Tin+Tini.Tb is the reference temperature and the conditions of vin = 6 m/s, Tin = 873.15 K, and Tini = 663.15 K are chosen as the base conditions. Five combinations of inlet condition at an identical input power are chosen, as shown in Table 3, and the impact of Ste and Re numbers on the energy storage is studied.Table 3
The working conditions at the same inlet power.
Conditions Ste Re 1 0.737 251 2 0.823 204 3 0.881 176 4 0.957 152 5 1.06 130As shown in Figure18, the variation trend is quite similar among the five outlet temperature curves, but overall the outlet temperature rises more quickly with small Ste number (i.e., larger Re number). The outlet temperature (To) remains in low constant value in the early stage which is equal to the Tini and then rises slowly at moderate rates. When the liquid fraction approaches around 0.4, the To reaches a constant value which is close to the phase change temperature of the PCMs, but it can be seen that the onset of the constant temperature period becomes more delayed with decreasing Re numbers. At the end of the charging process, the outlet temperature rises quickly. The tendency of the outlet temperature can be explained as follows. At the beginning, the temperature of the PCMs is low and the heat is absorbed and stored in sensible form, leading to the constant To equaling Tini. At the phase change temperature of the PCMs, the material begins to melt while the heat is absorbed and stored in sensible and latent form, and due to the lower temperature of PCMs in the downstream region and the buffer effect of the melting process, Tout rises slowly but is still lower than the melting point of PCMs at this stage. In the next stage, the temperature of PCMs in the downstream region rises and To is controlled by the phase change temperature of PCMs, which leads to To equaling the melting temperature of PCMs. Because most of PCMs melt and the heat is stored mainly in sensible form in the PCMs, To rises quickly at the end of the charging process.Furthermore, the reaching time of constantTo is delayed for smaller Re numbers. This is because the temperature difference between the HTF and the PCMs increases with decreasing Re numbers, which made more thermal energy get stored in the sensible way in the inlet region, the melting speed is lower in the downstream area and the reaching time of constant value is delayed. But the time duration of constant value decreases with increasing Ste number (decreasing Re number). Also, Figure 18 shows the liquid fraction as a function of time under different conditions for fixed input power. One can see that the overall melting time decreases with increasing Re number (decreasing Ste number), but the increase tendency slows down. This is because, for a given input power, an increase in the Ste number results in an increase in the inlet temperature and a decrease in the inlet velocity. The difference between inlet temperature and initial temperature increases with decreasing Re number according to (12). This implied that more heat is stored in the inlet region as sensible heat and less heat flows to the downstream region to melt the PCMs there, which leads to a lower melting rate in the condition of smaller Re for a given inlet power.


## 5. Conclusions

A numerical modeling of the melting process of PCMs within zigzag plate is performed using the solidification/melting model validated by experimental data. The validity of the model is investigated numerically in this study. According to this work, the following conclusions have been drawn:(i)
A finite volume method with LES model is used in this study. The modeling results showed reasonable well agreement with the experimental data.(ii)
Zigzag shape can improve the heat transfer performance of the heat exchange plate compared to the smooth shape.(iii)
The tooth length and tooth rate almost have no impact on the heat transfer capability of the unit, but they both influence the pressure drop. Pressure drop decreases with the increase of tooth rate or tooth ratio.(iv)
The melting time decreases with increasing Re and Ste numbers but decreases with increasing Re number (decreasing Ste number) when employing a fixed input power.



---

*Source: 3092508-2016-09-25.xml*