Next Article in Journal
Numerical and Physical Simulation of Heat Transfer Enhancement Using Oval Dimple Vortex Generators—Review and Recommendations
Next Article in Special Issue
The Effect of Experimental Conditions on Methane (95%)–Propane (5%) Hydrate Formation
Previous Article in Journal
Look Ahead Based Control Strategy for Hydro-Static Drive Wind Turbine Using Dynamic Programming
Previous Article in Special Issue
Study on Effects of Inlet Resistance on the Efficiency of Scroll Expander in Micro-Compressed Air Energy Storage System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermal Calculation and Experimental Investigation of Electric Heating and Solid Thermal Storage System

1
School of Electrical Engineering, Shenyang University of Technology, Shenyang 110870, China
2
Shenyang Lanhao New Energy Technology Company, Shenyang 110006, China
*
Author to whom correspondence should be addressed.
Energies 2020, 13(20), 5241; https://doi.org/10.3390/en13205241
Submission received: 7 August 2020 / Revised: 28 September 2020 / Accepted: 2 October 2020 / Published: 9 October 2020
(This article belongs to the Special Issue Energy Storage Systems and Conversion Processes)

Abstract

:
Electric heating and solid thermal storage systems (EHSTSSs) are widely used in clean district heating and to flexibly adjust combined heat and power (CHP) units. They represent an effective way to utilize renewable energy. Aiming at the thermal design calculation and experimental verification of EHSTSS, the thermal calculation and the heat transfer characteristics of an EHSTSS are investigated in this paper. Firstly, a thermal calculation method for the EHSTSS is proposed. The calculation flow and calculation method for key parameters of the heating system, heat storage system, heat exchange system and fan-circulating system in the EHSTSS are studied. Then, the instantaneous heat transfer characteristics of the thermal storage system (TSS) in the EHSTSS are analyzed, and the heat transfer process of ESS is simulated by the FLUENT 15 software. The uniform temperature distribution in the heat storage and release process of the TSS verifies the good heat transfer characteristics of the EHSTSS. Finally, an EHSTSS test verification platform is built and the historical operation data of the EHSTSS is analyzed. During the heating and release thermal process, the maximum temperature standard deviation of each temperature measurement point is 28.3 °C and 59 °C, respectively. The correctness of the thermal calculation of the EHSTSS is thus verified.

1. Introduction

In recent years, renewable energy utilization has developed rapidly in China. However, the abandonment of the renewable electricity generation is serious, especially during the wintertime heating period [1]. In the northern parts of China, CHP units are used for central heating in winter. The power output of CHP unit is greatly constrained by the heating load according to the principle of “power determined by heat” [2]. This has become one of the main reasons for the consumption of renewable energy in winter. The decoupling of the heat and power control of the CHP unit can be realized by adding a large amount of electric heating and heat storage unit to the CHP unit, which can effectively improve the adjustment capacity of the CHP unit and the renewable consumption capacity of the grid [3,4]. In addition, the traditional coal-fired boiler heating system for urban and factory thermal energy supply can also be replaced by the EHSTSS. It has become one of the effective ways to control air pollution. In the multi-energy system of combined cooling, heating and power, the thermal storage system can also be used as an energy storage system and a thermal energy supply system [5]. It can improve the diversity of thermal energy supply in a multi-energy system. In various thermal energy storage system, sensible heat storage is relatively popular because of its simple technology application and low cost [6]. More research on sensible heat storage heating systems mainly include the EHSTSS and the electrode boilers. Compared with the electrode boiler, the EHSTSS has advantages in heat storage density, safety, area covered and electrothermal response speed. It has become a research hotspot in recent years. The thermal calculation and the heat transfer characteristics analysis of the EHSTSS is a key link in the design and manufacture of the equipment. Similar to the thermal calculation of a traditional coal-fired boiler, the thermal calculation of the EHSTSS is mainly to calculate the key parameters of the equipment when the rated parameters are known [7].
Thermal energy storage (TES) is a technology under investigation since the early 1970s [8]. At present, the research on EHSTSS mainly focuses on the thermal storage medium, heat transfer characteristics of the system and structural optimization of the system. However, less attention has been paid to the thermal calculation of EHSTSS. Zhu et al. [9] report a helical-coil tube heat exchanger. The temperature difference of thermal energy storage is analyzed with and without considering radiative heat transfer. Gasia et al. [10] proposed a heat transfer enhancement technique that adds cheap and commercially available metallic wool. The thermodynamic performance of heat exchanger is analyzed from heat transfer process and enhancement of heat transfer characteristics in the above literature. However, only a single component is analyzed for thermodynamic performance, and the calculation method and process of heat exchanger thermal parameters are not investigated. Mousavi, et al. [11] investigated the melting process of phase change materials in an internal melt ice-on-coil thermal storage system. The effects of different operating parameters such as the inlet temperature and flow rate of the heat transfer fluid are analyzed, but there is no design calculation and heat transfer performance analysis of the heat storage module. Fadl et al. [12] present a latent heat thermal energy storage system. The influence of different heat transfer fluid inlet temperatures and volume flow rates of the system is evaluated by the experimental investigation. However, it lacks the heat exchanger parameter thermal design calculation process, and the optimization parameter curve of the overall system is not very obvious. In the aspect of TSS design, the existing research mainly focuses on the thermal calculation of parts of the TSS, such as separate heat exchangers, heat storage modules, etc. Raczka et al. [13] and Fujii et al. [14] introduced the design and calculation method of flue gas/waste water-heat exchanger and indirect heat exchanger, respectively. Du et al. [15] also reported a design method for heat storage units. However, the abovementioned research lacks any thermal design calculation of heat storage equipment from the perspective of a complete system and further systematic verification of equipment performance. It is necessary for the whole system as the research object to carry out thermal design calculation and systematic verification for EHSTSS.
Therefore, a systematic thermal calculation method for the EHSTSS is presented by this work. According to the structural characteristics of the system, thermal calculation of the system mainly calculates the parameters of the heating element, the TSS, the fan-circulating system and heat exchange system. Meanwhile, the fluid-solid coupling characteristics of the heat transfer in the TSS are analyzed. The numerical simulation of temperature distribution is carried out in the condition of heat storage and heat release. An experimental correlation for the EHSTSS is derived in order to verify correctness of the thermal calculation.
The two main contributions of this paper are summarized as follows:
  • A thermal calculation process specifically designed for the EHSTSS is proposed.
  • Systematic verification of the rationality and correctness of the EHSTSS from three aspects: case design, simulation analysis and experimental verification.
  • The multi angle and systematic verification results can provide the basis for the optimal design of energy storage system.
The rest of this paper is organized as follows: Section 2 presents the thermal calculation method and process of the electric heating and heating storage system. Section 3 analyses the heat transfer characteristics of the TSS. Experimental verification is presented in Section 4. Section 5 concludes the paper.

2. Thermal Calculation Flow and Method of the EHSTSS

2.1. Thermal Calculation Flow of the EHSTSS

The EHSTSS is composed of the TSS, including a thermal storage module and embedded heating element, the heat exchanger, the frequency converter fan, the thermal insulation layer, an external controller, etc. A structural diagram is shown in Figure 1. The overall dimensions of EHSTSS, i.e., length, width and height are 1560 mm, 720 mm and 1100 mm, respectively. In the EHSTSS, the thermal storage unit uses a solid sensible thermal storage medium such as magnesium oxide. The length, width and height of the thermal storage unit are 240 mm, 115 mm and 53 mm, respectively. The heating elements using NiCr or FeCrAl materials are embedded in the thermal storage module [16]. In the thermal storage process of EHSTSS, the thermal is generated by the heating element and stored in the thermal storage module. The thermal is extracted from the thermal storage module by the frequency conversion fan, and the thermal is exchanged through the heat exchanger to meet the heat load demand, when the heat load is in demand [17]. The heat transfer process mentioned above, the heat flow of the EHSTSS is mainly composed of heat generation, thermal storage, heat extraction and heat exchange. Therefore, taking heat flow and system structure into account, the EHSTSS is divided into heating system, thermal storage system, fan-circulating system and heat exchange system.
The four subsystems of the EHSTSS are interrelated. The thermal calculation of the EHSTSS should clarify the parameter relationship between the subsystems. The parametric relationship among the subsystems of the EHSTSS is shown in Figure 2.
Thermal calculation of EHSTSS usually start with thermal storage system. Before that, the key parameters of the EHSTSS need to be determined, mainly including thermal storage capacity, heating power, initial and final temperature of thermal storage module, temperature of supply and return water, heating time, heating voltage, etc.
Thermal calculation of the thermal storage system mainly determines the number of the thermal storage units and their arrangement. When the arrangement of thermal storage units is determined, the number of thermal storage channels and structure parameters of the TSS can be obtained. The number of thermal storage channels and structure parameters of the TSS are the vital input parameters of the heating system. The length and surface load of the heating element are calculated according to heating voltage, heating power, number of the thermal storage channel and other parameters. The parameter of the heat exchange system is affected by the maximum heat load, the temperature of upper and lower air ducts determined by the thermal storage system. The parameter of fan-circulating system is mainly affected by the air flow and the flow resistance which is determined by heat exchanger system and the thermal storage system.
The thermal calculation process for the four subsystem of the EHSTSS discussed above is shown in Figure 3. The calculation methods of key parameters in each subsystem are investigated below.

2.2. Thermal Calculation of Key Parameters of EHSTSS

Thermal calculation is mainly aimed at calculating the key parameters in the four subsystems of the EHSTSS. The following is a detailed investigation to the calculation of key parameters in the subsystem.

2.2.1. Thermal Calculation of Key Parameters in Heating System

The heating power, the length of heating wire and the surface load of the heating element are the key parameters in the heating system. The heating power is the basic parameter of the heating system, which is mainly related to the building area, the thermal index of building heating and the heating time. The heating power Ph (kW) of the heating system can be determined by Equation (1):
P h = 24 P f F 1000 t 1 η
where F represents heating area, m2; Pf represents the value of heating index (W/m2 ) shown in Table 1, t1 represents heating time, h; η thermal efficiency of the system, thermal efficiency is the ratio of the thermal Q1 generated by the heating element to the effective thermal Q released by the system, i.e., η = Q1/Q.
When the physical parameters and diameter of the heated wire are determined, the length of the heating wire used for a single heating element is the key parameter related to the parameters such as heating voltage and heating power. The length L of the heating wire can be expressed as follows:
L = 3 π U 2 D 2 4 k ρ Ω N P h
where U represents the rated voltage, V; ρΩ represents the resistivity, μΩ × m, k represents a temperature coefficient, N represents the heating element number, D represents the diameter of the heated wire, mm.
The surface load of the heating wire is the parameter that affects the service life of the heating element, so the rationality of the selection of the heating element can be judged according to the surface load. As the surface load increases, the temperature of heating element rises. In general, in high temperature heating applications, the surface load of heating elements can be controlled at 3~8 W/cm2. The surface load W of the heating wire can be expressed as follows:
W = P h 3 π L D

2.2.2. Thermal Calculation of Key Parameters in Thermal Storage System

In a thermal storage system, the key thermal calculation parameters are mainly the number and arrangement of thermal storage units. At the end of the preheating of the TSS, the average temperature is recorded as T0 °C. When the thermal storage module is heated to a set temperature, the average temperature is T °C. When the physical and structural parameters of the thermal storage unit are known, the number n of thermal storage units can be expressed as follows:
n = m P h t 1 ρ s C V ( T T 0 )
where C represents the thermal storage unit specific heat, kJ/(kg × °C); ρs represents the thermal storage unit density, kg/m3; V represents the heat storage unit volume, m3 and m represents the margin of thermal storage capacity.
According to the number of thermal storage units, the arrangement of TSS is calculated. The horizontal row number a and the height row number d of the thermal storage unit are determined, and the longitudinal row number e can be expressed as follows:
e = [ Q a d ρ s C V ( T T 0 ) + 0.5 ]
where Q represents thermal storage capacity, kWh; Q = mP1t1.

2.2.3. Thermal Calculation of Key Parameters in Heat Exchange System

The key parameters of the heat exchange system mainly include the heat transfer area, the number of heat exchange tubes and the air flow rate in the heat exchange tube.
In the heat exchanger, the heat transfer area of the heat exchanger is related to the heat transfer coefficient, the temperature difference of the fluid and the maximum heat load. The calculation of heat transfer area A can refer to Equation (6):
A = Q t 2 α × Δ T
where α represents heat transfer coefficient of heat exchanger; t2 represents the fastest heat release time, h; Δt represents logarithmic temperature difference, °C; which can be referred to Equation (7) [19,20];
Δ T = ( T 2 T 3 ) ( T 1 T 4 ) ln [ ( T 2 T 3 ) / T 1 T 4 ]
where T1, T2 represents air temperature at inlet and outlet of heat exchanger, °C; T3, T4 represents the supply and return water temperature of the heat exchanger, °C.
In Equation (6), the heat transfer coefficient α of heat exchanger can be expressed as:
α = 0.023 ε D dl ( ω D dl μ ) 0.8 P r 0.4 c t c l ψ
where ε is the thermal conductivity at the average temperature of the air, kW/(m × °C);μ is the kinematic viscosity of air at average temperature, m2/s; Ddl is the outer diameter of the heat exchange tube, m; Pr is the Prandtl number of air at average temperature; ct and cl are pipeline correction coefficients; ψ is thermal efficiency coefficient.
The number of heat exchanger tubes is related to the heat transfer area, the length and diameter of the tubes. The length of the heat exchanger tube is suitable to the diameter, and the ratio of the tube length (L1) to the diameter (D1) is about 4~6 [19]. For the thermal calculation of the number of heat exchanger tubes we can refer to Equation (9):
N 1 = A π D 1 L 1
where N1 represents the number of heat exchanger tubes.
The flow resistance of heat exchanger accounts for a large proportion of the flow resistance of thermal storage system. The velocity of air in the tube is the key parameter affecting flow resistance.
In the heat exchange tube, the velocity of air under standard state can be expressed as follows:
ω 0 = S af 1 A 1 × 3600 × 1.293
where A1 represents cross-sectional area of heat exchange tubes, m2; Saf1 is the air mass flow at the outlet of heat exchanger, kg/h; which can be expressed as:
S af 1 = P h × 3600 C 1 T 1 C 2 T 2
where C1, C2 the specific heat of the air at the outlet and inlet of the heat exchanger, kJ/(kg × K).
The actual velocity of air in heat exchanger tube can be calculated as (12):
ω = ω 0 [ ( T 1 + T 2 ) + 273 ] / 273

2.2.4. Thermal Calculation of Key Parameters in Fan-circulating System

In a fan-circulating system, the motor power of the frequency conversion fan is an extremely important parameter. The motor power is determined by the required flow resistance and flow. In the EHSTSS, flow resistance consists of two parts. One part is the high temperature channel flow resistance ∆Pf1(Pa) and the back flow resistance ∆Pf2 (Pa) of heat exchange system, the other part is low temperature channel flow resistance∆Pf3(Pa), flow resistance ∆Pf4 (Pa)of the heat exchanger and the thermal storage channel flow resistance∆Pf5(Pa). The total resistance of the EHSTSS can be expressed as follows:
Δ P = Δ P f 1 + Δ P f 2 + Δ P f 3 + Δ P f 4 + Δ P f 5
where [21,22]:
Δ P f 1 , 2 = λ I ω 2 ρ a 2 d e ( 2 ( T w / T a v ) 0.5 + 1 ) 2
Δ P f 3 , 4 , 5 = λ I ω 2 ρ a 2 d e
where λ represents friction coefficient in channel; ρ2 represents medium density, kg/m3; ω represents air velocity, m/s; I represents length of channel, m; de represents channel section equivalent diameter, m; Tw represents thermal storage unit surface average temperature, °C; Tav represents air average temperature, °C.
The air volume flow Saf (m3/h) required for the fan-circulating system can be referred to Equation (16):
S af = 10.20 S af 1 ( T 2 + 273 ) 273 × 1.293
When the air pressure and air flow required for the fan-circulating system are determined, the motor power Pm(kW) can be calculated as follows:
P m = 1.1 Δ P × S af 3600 × 1020 δ
where δ represents motor efficiency.
In order to verify the rationality and correctness of the thermal calculation process and method of the EHSTSS, the case design is based on the system with heating power of 100 kW, daily thermal storage of 1MWh and heat transfer power of 125 kW. The case design is shown in Appendix A.

3. Analysis of Heat Transfer Characteristics of EHSTSS

The heat transfer process of EHSTSS is a complex process which includes heat conduction, convective heat-transfer and radiation heat-transfer [23]. There are different heat transfer modes between structures of the system, as shown in Figure 4.
The operation process of EHSTSS is an alternating process between thermal storage and thermal release. The temperature of the TSS is a key index to judge whether design parameters of the system are suitable. Therefore, the heat transfer characteristics analysis of the system is carried out to explore the temperature variation rule of the heat storage device.
In the thermal storage process of an EHSTSS, heat transfer methods include convective heat-transfer between heating elements and cold air, radiation heat-transfer between heating elements and thermal storage module and heat conduction inside thermal storage module. In the thermal release process of the system, heat transfer methods include heat conduction inside thermal storage module and convective heat-transfer between cold air and thermal storage module. The heat transfer process in the system involves two regions of solid and fluid. Therefore, the analysis of the heat transfer characteristics between the two regions is of great significance to improve the heat storage efficiency of the heat storage module.

3.1. Heat Transfer Analysis in Thermal Storage Process

The thermal in the thermal storage module is derived from the heating element. The thermal generated by the heating element interacts with the heat storage module through radiation heat exchange and convection heat exchange. The net heat Qs in the thermal storage module can be expressed as follows:
Q s = Q es + Q ea Q a = C M s ( T 6 T 5 )
where Qes represents radiant heat between heating element and thermal storage module, W; Qea represents convective heat between resistance element and air, W; Qa represents heat carried by air, W; Ms is the weight of the heat storage module, kg; T5 and T6 represent initial temperature and temperature at any time of thermal storage module during heating, °C.
The radiant heat Qes between the resistance element and the thermal storage module can be expressed as follows [24]:
Q es = β [ ( T 7 100 ) 4 ( T 5 100 ) 4 ] F 1
where T7 represents heating element temperature, °C; β represents radiation coefficient, F1 represents radiation surface area of heating element, m2.
For the convective heat Qea between resistance element and air we can refer to Equation (20):
Q ea = γ ( T 7 T 1 ) F 2
where γ represents convective heat transfer coefficient, F2 represents heating wire area, m2.
The thermal Qa carried by air can be written as:
Q a = q m , h V p , h ( T 6 T 2 )
where qm,h represents the mass velocity of the air, m/s; Vp,h represents the constant pressure specific volume of the air, m3. Applying Equations (19)–(21) into Equation (18) gives:
T 6 ( C M s + q m , n V p , h ) = β [ ( T 7 100 ) 4 ( T 5 100 ) 4 ] + γ F 2 T 7 γ F 2 T 1 + q m , h V p , h T 2 + C M s T 5
The temperature T6 of the thermal storage module at a certain moment can be calculated by Equation (22).

3.2. Heat Transfer Analysis in Thermal Release Process

In the thermal release process, the thermal source of the system is mainly the thermal storage module. The thermal released Qr by the thermal storage module is equal to the sum of the effective thermal Qc and the lost thermal of the thermal storage system. The thermal released Qr can be written as:
Q r = Q c / η
The thermal released Qr also can be written as Equation (24) [24]:
Q r = C M s ( T 8 T 9 )
where T8 and T9 represent the initial temperature and the temperature at the moment of the thermal storage module during heat release, °C. The effective heat released by the system can be expressed as Equation (25):
Q c = M h C p . h ( T 3 T 4 )
where Mh represents the mass flow rate, kg/s; cp.h represents the specific heat, kJ/ (kg× K).
Applying Equations (24) and (25) into Equation (23) gives:
T 9 = T 8 M h C p . h ( T 3 T 4 ) C M c η
The temperature T9 of the thermal storage module at a certain moment can be solved by Equation (26).

3.3. Simulation Results and Discussion

3.3.1. Parameters of Geometry Model

According to the thermal calculation process and method of the EHSTSS proposed in the paper, the physical model of the TSS is built. The physical model parameters are set as shown in Table 2.
In order to measure the thermal storage module and fluid temperature in the system, six sets thermocouples are arranged to monitor the temperature of the thermal storage module and air in the channel as shown in Figure 5.

3.3.2. Mathematical Model

According to the previous analysis, the heat storage system is mainly divided into the fluid region formed by the incompressible high-temperature air and the solid region formed by the heat storage module. In the fluid region, the fluid heat transfer process can be described by three equations: the mass conservation equation, the momentum conservation equation and the energy conservation equation [25].
The mass conservation equation is:
ρ a t + div ( ρ a u j ) = 0
The momentum conservation equation is:
ρ a d v d t = F p + ξ 2 v
The energy conservation equation is:
ρ a C p . h d T a d t = 2 ( λ a T a ) + ϕ + S i
where Ta is the air temperature in the thermal storage channel, K; v is the air velocity in the thermal storage channel, m·s−1; λa is the thermal conductivity of air, W·m−1·K−1; F is the mass force on the fluid, N; p is the fluid pressure, Pa; ζ is the dynamic viscosity of air, kg·m−1·s−1; ϕ is the loss equation; t is the time, s; Si is the energy source term. In this study, the energy source term is not considered, so Si is set as 0.
The heat storage module transfers energy by heat conduction. The heat conduction problem of the heat storage module can be simplified as a steady-state heat conduction problem. and its governing equation can be expressed as:
ρ s C T s t = λ s ( 2 T s 2 x + 2 T s 2 y + 2 T s 2 z ) + S j
where λs is the thermal conductivity of the thermal storage unit, W·m−1·K−1; Ts is the temperature of the thermal storage module, K; x, y and z are the coordinate values, m; Sj is the heating power per unit volume of the thermal storage module, W.
On the thermal exchange interface between the high-temperature air and the thermal storage channel, the directly coupled fluid-solid interface must meet the energy continuity condition, that is, the temperature and heat flux density of the air and the heat storage channel are equal. It can be expressed as:
{ T a = T s q a = λ a ( T a n ) = λ s ( T s n ) = q s
where qa is the heat flux density on the fluid side, J·m−2·s−1; qs is the heat flux density on the solid side, J·m−2·s−1; n is the fluid-solid interface normal vector.

3.3.3. Mesh Generation and Boundary Conditions

The model was calculated using the commercial simulation software FLUENT15. The tetrahedral structured grid is used to divide the heat storage system. The grid division of heat storage system is shown in the Figure 6. The mesh irrelevance test is implemented under the 4 types of grid numbers. After the simulation process reaches the steady state, the results are shown in the Table 3. Compared with the high number of grids, there is a certain difference in the temperature simulation results when the number of grids is less than 1.06 million, and the simulation results change little when the number of grids exceeds 1.06 million. Therefore, the calculation speed and time is considered comprehensively, and the number of meshes is determined to be 1.06 million.
In the numerical simulation, time step is 10. The air flowing in the heat storage channel and its front and rear cavities is a turbulent flow, and the k-ε turbulence model is applied. The unsteady Reynolds time-averaged N-S equation is used in simulation. The coupling between velocity and pressure is realized by SIMPLIC algorithm. The finite volume method is used to discretize the governing equations. The convection term difference scheme adopts the second-order upwind scheme. Dimensionless residuals of continuous equations reduced to less than 1 × 10−3 are the convergence criteria.
Thermal storage: Heating power is 270 kW. Absorption coefficient is 0.7. Cold air flow rate is 1 m/s. The outlet is set to the pressure outlet and the boundary condition of the wall is adiabatic and no slip. The interface between the fluid region and solid region is set as coupling interface. The total thermal storage time is 26.190 s.
Thermal release: The heating power is set to 0 kW. The inlet of cold air speed is 28 m/s and air temperature is given according to field test data. Thermal release time is 41.130 s.

3.3.4. Simulation Results Analysis

In the thermal storage process, the temperature variation of the TSS is shown in Figure 6. The temperature near the heating element is higher, and the temperature of the outlet section is more uniform. At the initial moment of heating, as shown in Figure 7a–c, the temperature at the inlet section emerges a peak, and the direction of the peak points to the inside of the thermal storage channel. It is caused by the difference of air velocity due to the accelerated flow of cold air inside the device and the uneven temperature distribution at the inlet of the heat storage channel. Because of the turbulence at the outlet of the thermal fluid, the temperature of the thermal fluid at the outlet of is lower. When the heating time reaches the set value, as shown in Figure 7d, the temperature of the thermal storage module is around 900 K, and it reaches the preset temperature. Meanwhile, the heating element temperature is approximately 1100 K, which is lower than the maximum temperature of the element. In conclusion, the temperature distribution of the thermal storage module is relatively uniform. The design and arrangement of the module and heating elements are reasonable, which verifies correctness of the thermal calculation of the system.
In the heat release process, the temperature distribution of the thermal storage module is shown in Figure 8. At the end of the thermal release, the temperature of module is uniformly distributed around 380 K, which is consistent with the preset temperature. The temperature near the inlet is lower than other locations, which is caused by the rapid heat release due to the low air temperature and fast flow rate at the entrance. In conclusion, the TSS has a good match with air flow and velocity provided by the fan-circulating system, which makes the heat storage module complete the heat release within the setting time. The correctness and rationality of the thermal calculation and guidance for the EHSTSS is proved.

4. Experimental Results and Discussion

In the paper, the operating data of the EHSTSS with a rated power of 125 kW in Anshan City, (Liaoning Province, China) are extracted and processed. The external and internal structure of the experimental equipment is shown in Figure 9. The heating area of the equipment is 1000 m2. The average outdoor temperature in winter at the location of the equipment is −0.2 °C and the heating time is 114 days. In order to make the temperature of the measuring point more accurate at the testing time, the temperature at the sampling time and the temperature at 20 s before and after the sampling are measured, and the average of the three values is taken as the temperature at this time. The control parameters of the device are adjusted to stabilize the various operating data of the equipment. The operation data of the device under stable conditions was used.
In the experiment, thermocouples are used to measure the temperature at different positions of the thermal storage module, and the positions the arrangement of thermocouples refer to in Figure 5.

4.1. Analysis of Experimental Results Under Thermal Storage

In the thermal storage process, temperature variation at each measuring point is shown in Figure 10. At the end of thermal storage, the temperature of measuring point 2 is the highest, and the temperature value is 715 °C. The lowest temperature was measured at point 4, with a temperature value of 683 °C. The maximum temperature difference of the measuring point is 33 °C. The temperature standard deviation of each temperature measurement point at each moment is calculated, as shown by the dotted line in Figure 10a. In the middle of thermal storage, the temperature deviation at each point is relatively large, the maximum is 28.3 °C. At the beginning and end of thermal storage, the temperature deviation of the each point is small, only about 10 °C. Therefore, the temperature of the thermal storage module is relatively uniform, and the thermal stress of the thermal storage unit is smaller. The TSS design is reasonable.
At the same time, the test value and simulation value of the temperature at measuring points 1, 3, and 6 are compared and analyzed, as shown in Figure 10b. The standard deviation statistics of the simulation value and the test value error rate of each measuring point are investigated, and the error deviations of the simulation value and the test value of the measuring points 1, 3, and 6 in the whole thermal storage process are respectively 0.796%, 0.925% and 0.805%. The simulation value is close to the experimental value, so it can be seen that the EHSTSS designed in the paper has good heat transfer performance.
The average temperature curve of the thermal storage module in different typical working days is shown in Figure 11. The average temperature of the thermal storage module rises from 130 °C to about 700 °C during the heating time of 10 h. The average temperature variation of the typical day has a small fluctuation, the temperature difference of the thermal storage module at the same time is small. The maximum temperature deviation is only 33.2 °C. It can be seen that the EHSTSS has good operational stability, and it can meet the design requirements.
The water supply temperature is constant as the control target. The average temperature of the TSS and frequency curve of the fan are shown in Figure 12. With the average temperature rising gradually in the preset range, the frequency of the fan declines, and the water supply temperature is maintained at the setting temperature about 50 °C. It can be obtained that the parameters of heating system, thermal storage system, heat exchange system and fan-circulating system have good adaptability and compatibility. They ensure that the temperature change of the thermal storage module and the supply water can be stabilized in the set range.

4.2. Analysis of Experimental Results Under Thermal Release

In the thermal release process, temperature at each measuring point of the module is shown in Figure 13. At the end of thermal release, the temperature of measuring point 5 is the highest, and the temperature value is 121 °C. The lowest temperature was measured at point 3, with a temperature value of 106 °C. The maximum temperature difference is 15 °C. The temperature standard deviation of each measuring point of the system at different times is analyzed. The maximum standard deviation of the system is 19.2 °C, and the temperature standard deviation of the system is maintained at about 10 °C for most of the time. The thermal storage module of the system is proved to have good temperature uniformity, and the structure of thermal storage module is reasonable. The coordination between air circulation system and thermal storage system is good.
The test value and simulation value of the temperature at measuring points 1, 3, and 6 are compared and analyzed in the thermal release process, as shown in Figure 13b. The standard deviation statistics of the simulation value and the test value error rate of each measuring point are analyzed, and the error deviations of the simulation value and the test value of the measuring points 1, 3, and 6 in the whole thermal storage process are respectively 8.6%, 6.7% and 3.34%. The simulated value has a certain deviation from the experimental value. This is mainly caused by the unavoidable heat leakage in the thermal insulation layer design during the experiment of the EHSTSS. Therefore, more attention should be paid to thermal insulation design in the system.
The average temperature variation curve of the thermal storage module on different typical working days is shown in Figure 14. The average temperature of thermal storage module on the typical daily varies from 700℃ to 130 °C during the thermal release time of 14 h. At the same time, the standard temperature deviation changes relatively large between each day, and the maximum temperature deviation reaches 59 °C. Therefore, the operating fluctuation of the system during the heat release process is relatively large, and the system can be optimized from the angle of thermal release.
The water supply temperature is constant as the control target. The average temperature of the TSS and frequency curve of the fan are shown in Figure 15. The average temperature gradually decreases in the preset range. Meanwhile, the frequency of the blow rises. The temperature of the supply water still remains at a preset temperature of about 50 °C. It can be concluded that the parameters of subsystem in the EHSTSS are reasonable, and each part has good compatibility with each other. Therefore, supply water temperature can be maintained stability by adjusting frequency of the fan to control the heat release.

4.3. EHSTSS Electric Heating Characteristics and Pressure Drop Analysis

The electrothermal conversion process of EHSTSS under operation condition is analyzed. The average temperature of each measuring point for the heat storage module and the power of EHSTSS are shown in Figure 16 within 72 h. As shown in Figure 16 before the temperature of the heat storage module reaches the set value, the EHSTSS operates at full power to ensure that the system can store heat while outputting heat. After reaching the set temperature, in order to ensure that the temperature of the system is maintained at the set temperature before pure heat release, EHSTSS operates at half power. When the system is only exothermic, the power output of EHSTSS is stopped. The power output of EHSTSS changes periodically with the temperature of the heat storage module, which verifies that the system has good heat exchange and heat transfer performance.
The pressure drop in the heat storage system under different wind speeds is tested by adjusting the speed of the hot air fan. The change curve is shown in Figure 17.
As shown in Figure 17, the pressure drop of the heat storage system is 42.4 Pa at the rated wind speed of the fan, which is close to the calculated value 38.7 Pa. At the same time, the calculated value of the pressure drop of the heat storage system at other wind speeds is also close to the test value. The rationality of fan selection and channel design is verified.

5. Conclusions

In the paper, the thermal calculation method and heat transfer characteristics of an EHSTSS is studied. The key parameters of the subsystems in the EHSTSS were investigated, and the feasibility of the proposed method was verified by a case design. The following conclusions are drawn through the actual case design, simulation analysis and experimental verification of the EHSTSS.
Systematic verification shows that within the rated heating and thermal release time, the internal thermal of the thermal storage module can be fully stored and released. The maximum temperature of the thermal storage module is controlled at about 700 °C, and the outlet water temperature is stabilized at about 50 °C. The design system can fully meet the design requirements.
During the heating process, the temperature gradient of each measuring point of the system is small. The maximum temperature standard deviation of each temperature measurement point at each moment is 28.3 °C. The simulation value is consistent with the test value, and error deviations of the temperature is only 0.796%, 0.925% and 0.805%. The good heat transfer performance of the system is proved.
In the thermal release process, the temperature distribution of the system is uniform, and the temperature standard deviation of the system in each measuring point is maintained at about 10 °C for most of the time. However, the standard temperature deviation changes relatively large between typical day in each month, and the maximum temperature deviation reaches 59 °C. This phenomenon indicates that there are some defects in the long-term operation stability of the system, which needs to be further improved. Corresponding parameter optimization can be made in terms of the number of channels of the heat storage module, the channel structure, and the location of the air inlet and outlet of the device.
From the perspective of case design, simulation analysis and experimental verification, the thermal calculation method, process and heat transfer characteristics analysis proposed in this paper are suitable for the design of an EHSTSS and it can provide a reference for the design and verification method of EHSTSS.

Author Contributions

All authors conceived and designed the study. H.Z. performed the experiment verification, the analysis of results and wrote the manuscript with guidance from N.Y., Z.X. and L.J., L.C. performed the simulation of the study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by “The National Key Research and Development Program of China” (2017YFB0902100).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Parameter design results of EHSTSS case.
Table A1. Parameter design results of EHSTSS case.
NameValueNameValue
Basic parameter calculationConfiguration of heating power/kW100Enthalpy value of return water/kJ·kg−1189
Heating time/h10Enthalpy value of supply water/kJ·kg−1231
Maximum heat release time/h8Pipe diameter/mm50
Inlet temperature of water/℃45Maximum heating load/ kW125
Outlet temperature of water/℃55Total thermal storage capacity/kW·h1000
Design power/kW100Number of single phase heating elements11
Calculation of heating element parameterThree-phase voltage/V400Single-phase voltage/ V231
Total number of heating elements33Resistance of heating element (single)/Ω0.57
Diameter of heating elements/ mm3Heating power(single)/ kW0.74
Temperature coefficient1.08Length of heating element/ mm2474
Resistivity/Ω·mm−2·m−14.79Surface load of heating element/ W·cm−23.18
Calculation of thermal storage module structure designThermal storage margin1.1Total thermal storage capacity /kW·h1100
Specific heat capacity of thermal storage unit/ kJ·kg−1·℃−11.064Volume of thermal storage unit/ m³0.0041
Final average temperature of thermal storage module/℃700Weight of thermal storage unit/ kg11.55
Initial average temperature of thermal storage module/℃150Design quantity of thermal storage unit716
Horizontal row number of thermal storage unit6Longitudinal row number7
Height row number of thermal storage unit20Actual number of thermal storage unit724
Thermal storage unit density/kg·m−32800Weight of thermal storage module/ kg8357
Thermal storage capacity of thermal storage unit/ kW·h1.71Real thermal storage capacity/ kW·h1111
Calculation of heat exchanger selectionDesign heat transfer load/ kW125Outlet air flow of heat exchanger/ m³·h−1605
Inlet air temperature of heat exchanger/℃650Inlet air flow of heat exchanger/ m³·h−11967
Outlet air temperature of heat exchanger/℃105Heat transfer area/m2115
Heat transfer coefficient/kW·m−2·k−11.1Heat transfer area of heat exchanger pipe/m20.503
Length of heat exchanger tube/mm1000Number of heat exchanger tubes5
Diameter of heat exchanger tube/mm25Real air velocity in heat exchanger tube/ m·s−127.26
Calculation of frequency conversion fan selectionFlow resistance of thermal storage channel/Pa6Total flow resistance/ Pa20177
Flow resistance of heat exchanger/Pa19688Outlet air flow of circulating fan/ m³·h−11969
Flow resistance of low temperature duct/Pa11Motor efficiency0.75
Flow resistance of high temperature duct/Pa22Calculation power of circulating fan/kW6.5
Flow resistance after heat exchanger/Pa450Real power of circulating fan/kW7.5

References

  1. Wang, H.X.; Yang, J.Y.; Chen, Z.; Li, G.; Liang, J. Optimal dispatch based on prediction of distributed electric heating storages in combined electricity and heat networks. Appl. Energy 2020, 267, 1–10. [Google Scholar] [CrossRef]
  2. Dai, Y.; Chen, L.; Min, Y.; Chen, Q.; Hu, K.; Hao, J.; Zhang, Y.; Xu, F. Dispatch model of combined heat and power plant considering heat transfer process. IEEE Trans. Sustain. Energy 2017, 8, 1225–1236. [Google Scholar] [CrossRef]
  3. Takahiro, N.; Tomohiro, A. High-temperature latent heat storage technology to utilize energy of solar heat and industrial exhaust heat. Int. J. Energy Res. 2016, 41, 240–251. [Google Scholar]
  4. Dai, Y.; Chen, L.; Min, Y.; Mancarella, P.; Chen, Q.; Hao, J.; Hu, K.; Xu, F. Integrated dispatch model for combined heat and power plant with phase-change thermal energy storage considering heat transfer process. IEEE Trans. Sustain. Energy 2018, 9, 1234–1243. [Google Scholar] [CrossRef]
  5. Liu, B.; Meng, K.; Dong, Z.Y.; Wei, W. Optimal dispatch of coupled electricity and heat system with independent thermal energy storage. IEEE Trans. Power Syst. 2019, 34, 3250–3263. [Google Scholar] [CrossRef]
  6. Liu, X.; Wu, J.; Jenkins, N.; Bagdanavicius, A. Combined analysis of electricity and heat networks. Appl. Energy 2016, 162, 1238–1250. [Google Scholar] [CrossRef] [Green Version]
  7. Dai, Y.H.; Chen, L.; Min, Y.; Mancarella, P. A general model for thermal energy storage in combined heat and power dispatch considering heat transfer constraints. IEEE Trans. Sustain. Energy 2018, 9, 1518–1528. [Google Scholar] [CrossRef]
  8. Fernandez, A.I.; Martinez, M.; Segarra, M.; Martorell, I.; Cabeza, L.F. Selection of materials with potential in sensible thermal energy storage. Sol. Energy Mater. Sol. Cells 2010, 94, 1723–1729. [Google Scholar] [CrossRef]
  9. Zhu, Y.; Yuan, Y.; Zhang, C.; Xie, M.; Tan, H. Numerical study on heat transfer enhancement of thermal energy storage systems considering radiation of molten salt. Sol. Energy 2019, 183, 337–344. [Google Scholar] [CrossRef]
  10. Mousavi, A.S.S.; Poncet, S.; Sedighi, K.; Mojtaba, A.D. Numerical modeling of the melting process in a shell and coil tube ice storage system for air-conditioning application. Appl. Sci. 2019, 9, 2726. [Google Scholar] [CrossRef] [Green Version]
  11. Fadl, M.; Eames, P.C. An experimental investigation of the heat transfer and energy storage characteristics of a compact latent heat thermal energy storage system for domestic hot water applications. Energy 2019, 188, 1–13. [Google Scholar] [CrossRef]
  12. Gasia, J.; Maldonado, J.M.; Galati, F.; Marilena, D.S. Experimental evaluation of the use of fins and metal wool as heat transfer enhancement techniques in a latent heat thermal energy storage system. Energy Convers. Manag. 2019, 184, 530–538. [Google Scholar] [CrossRef]
  13. Raczka, P.; Kazimierz, W. Methods of thermal calculations for a condensing waste-heat exchanger. Chem. Process Eng. 2014, 35, 447–461. [Google Scholar] [CrossRef]
  14. Fujii, S.; Horie, N.; Nakaibayashi, K.; Kanematsu, Y.; Kikuchi, Y.; Nakagaki, T. Design of zeolite boiler in thermochemical energy storage and transport system utilizing unused heat from sugar mill. Appl. Energy 2019, 238, 561–571. [Google Scholar] [CrossRef]
  15. Du, Y.X.; Xiao, G.M.; Gui, Y.W.; Lei, L.; Zhang, L.N.; Yu, M.X. Study on heat transfer characteristics of phase-change energy storage unit for thermal management. Int. J. Thermophys. 2014, 35, 1577–1589. [Google Scholar] [CrossRef]
  16. Barrios, G.; Huelsz, G.; Rechtman, R. Heat transfer and flow transitions of a thermal plume generated by a heating element on the enclosure bottom wall. Eur. J. Mech. B/Fluids 2019, 77, 17–24. [Google Scholar] [CrossRef]
  17. Yu, G.J.; Yu, B.; Liang, Y.T.; Wang, M.; Joshi, Y.; Sun, D.L. A new general model for phase-change heat transfer of waxy crude oil during the ambient-induced cooling process. Numer. Heat Transf. Part A Appl. 2017, 71, 511–527. [Google Scholar] [CrossRef]
  18. Design Code for City Heating Network: Industry standard of China CJJ 34-2010; China Building Industry Press: Beijing, China, 2010; No. CJJ 34-2010.
  19. Kim, M.H.; Lee, Y.T.; Gim, J.; Awasthi, A.; Chung, J.D. New effective thermal conductivity model for the analysis of whole thermal storage tank. Int. J. Heat Mass Transf. 2019, 131, 1109–1116. [Google Scholar] [CrossRef]
  20. Tumsa, T.Z.; Mun, T.Y.; Lee, U.; Yang, W. Effects of coal characteristics to performance of a highly efficient thermal power generation system based on pressurized oxy-fuel combustion. Int. J. Energy Res. 2016, 41, 127–138. [Google Scholar] [CrossRef]
  21. Gautam, A.; Saini, R.P. Experimental investigation of heat transfer and fluid flow behavior of packed bed solar thermal energy storage system having spheres as packing element with pores. Sol. Energy 2020, 204, 530–541. [Google Scholar] [CrossRef]
  22. Xu, H.J.; Zhao, C.Y. Analytical considerations on optimization of cascaded heat transfer process for thermal storage system with principles of thermodynamics. Renew. Energy 2019, 132, 826–845. [Google Scholar] [CrossRef]
  23. Guo, H.; Xu, Y.J.; Guo, C.; Chen, H.S.; Wang, Y.F.; Zheng, Y.; Ye, H. Thermodynamic analysis of packed bed thermal energy storage system. J. Therm. Sci. 2020, 29, 445–456. [Google Scholar] [CrossRef]
  24. Mosaffaa, A.H.; Talati, F.; Basirat Tabrizib, H.; Rosenc, M.A. Analytical modeling of PCM solidification in a shell and tube finned thermal storage for air conditioning systems. Energy Build. 2012, 49, 356–361. [Google Scholar] [CrossRef]
  25. Guo, D.C.; Liu, M.; Xie, L.Y.; Wang, J. Optimization in plate-fin safety structure of heat exchanger using genetic and Monte Carlo algorithm. Appl. Eng. 2014, 70, 341–349. [Google Scholar] [CrossRef]
Figure 1. Structure diagram of the EHSTSS.
Figure 1. Structure diagram of the EHSTSS.
Energies 13 05241 g001
Figure 2. Thermal calculation relation of energy storage system.
Figure 2. Thermal calculation relation of energy storage system.
Energies 13 05241 g002
Figure 3. Thermal calculation flow diagram of solid thermal storage system.
Figure 3. Thermal calculation flow diagram of solid thermal storage system.
Energies 13 05241 g003
Figure 4. Flow-solid coupling diagram of solid thermal storage system.
Figure 4. Flow-solid coupling diagram of solid thermal storage system.
Energies 13 05241 g004
Figure 5. Thermocouple arrangement diagram (left view).
Figure 5. Thermocouple arrangement diagram (left view).
Energies 13 05241 g005
Figure 6. Grid division of heat storage system.
Figure 6. Grid division of heat storage system.
Energies 13 05241 g006
Figure 7. Cloud profile of temperature distribution at section during heating. (a) heat for 2 h; (b) heat for 4 h; (c) heat for 6 h; (d) heat for 7.2 h.
Figure 7. Cloud profile of temperature distribution at section during heating. (a) heat for 2 h; (b) heat for 4 h; (c) heat for 6 h; (d) heat for 7.2 h.
Energies 13 05241 g007
Figure 8. The temperature distribution of the system after the end of the thermal release.
Figure 8. The temperature distribution of the system after the end of the thermal release.
Energies 13 05241 g008
Figure 9. Experimental facility.
Figure 9. Experimental facility.
Energies 13 05241 g009
Figure 10. Temperature variation curves of measuring points of thermal storage module in the thermal storage process (a) experimental value and deviation of temperature at measuring point; (b) curve of test point temperature test value and simulation value.
Figure 10. Temperature variation curves of measuring points of thermal storage module in the thermal storage process (a) experimental value and deviation of temperature at measuring point; (b) curve of test point temperature test value and simulation value.
Energies 13 05241 g010
Figure 11. Temperature curve of the heat storage module under the heating.
Figure 11. Temperature curve of the heat storage module under the heating.
Energies 13 05241 g011
Figure 12. Supply water temperature, thermal storage module average temperature and fan frequency variation curve under thermal storage.
Figure 12. Supply water temperature, thermal storage module average temperature and fan frequency variation curve under thermal storage.
Energies 13 05241 g012
Figure 13. Temperature variation curves of measuring points of thermal storage module in the thermal release process (a) experimental value and deviation of temperature at measuring point; (b) curve of test point temperature test value and simulation value.
Figure 13. Temperature variation curves of measuring points of thermal storage module in the thermal release process (a) experimental value and deviation of temperature at measuring point; (b) curve of test point temperature test value and simulation value.
Energies 13 05241 g013
Figure 14. Average temperature curve of the module under thermal release process.
Figure 14. Average temperature curve of the module under thermal release process.
Energies 13 05241 g014
Figure 15. Supply water temperature, thermal storage module average temperature and fan frequency variation curve under thermal release.
Figure 15. Supply water temperature, thermal storage module average temperature and fan frequency variation curve under thermal release.
Energies 13 05241 g015
Figure 16. EHSTSS power-temperature curve in 72 h.
Figure 16. EHSTSS power-temperature curve in 72 h.
Energies 13 05241 g016
Figure 17. Variation curve of aerodynamic loss in the thermal storage system.
Figure 17. Variation curve of aerodynamic loss in the thermal storage system.
Energies 13 05241 g017
Table 1. Recommended value of heating index [18].
Table 1. Recommended value of heating index [18].
TypeEnergy SavingNon Energy SavingTypeEnergy SavingNon Energy Saving
Residence40~4558~64Office50~7060~80
Hospital55~7065~80Hostel50~6060~70
Store55~7065~80Canteen100~130115~140
Cinema80~10595~115Gym100~150115~165
Table 2. Simulation parameters of the TSS.
Table 2. Simulation parameters of the TSS.
Parameter TypeValue
Thermal storage unit length, width and height240 mm × 115 mm × 53 mm
Thermal storage capacity (MgO)1000 kWh
Thermal storage unit specific heat960 J·kg−1·K−1
Bulk density2900 kg·m−3
Heat transfer coefficient2.7 W·m−1·K−1
Initial temperature of thermal storage module380 K
Final temperature of heat storage module 950 K
Table 3. Simulation results of the grid independence verification.
Table 3. Simulation results of the grid independence verification.
PorosityGrid TypeNumber of Grids/ MillionOutlet Temperature of Regenerator/KCore Temperature of Regenerator/K
15%Structured hexagonal grid0.9648688
1.06667702
1.17669704
1.28670705

Share and Cite

MDPI and ACS Style

Zhao, H.; Yan, N.; Xing, Z.; Chen, L.; Jiang, L. Thermal Calculation and Experimental Investigation of Electric Heating and Solid Thermal Storage System. Energies 2020, 13, 5241. https://doi.org/10.3390/en13205241

AMA Style

Zhao H, Yan N, Xing Z, Chen L, Jiang L. Thermal Calculation and Experimental Investigation of Electric Heating and Solid Thermal Storage System. Energies. 2020; 13(20):5241. https://doi.org/10.3390/en13205241

Chicago/Turabian Style

Zhao, Haichuan, Ning Yan, Zuoxia Xing, Lei Chen, and Libing Jiang. 2020. "Thermal Calculation and Experimental Investigation of Electric Heating and Solid Thermal Storage System" Energies 13, no. 20: 5241. https://doi.org/10.3390/en13205241

APA Style

Zhao, H., Yan, N., Xing, Z., Chen, L., & Jiang, L. (2020). Thermal Calculation and Experimental Investigation of Electric Heating and Solid Thermal Storage System. Energies, 13(20), 5241. https://doi.org/10.3390/en13205241

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop