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):
where
F represents heating area, m
2;
Pf represents the value of heating index (W/m
2 ) 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:
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/cm
2. The surface load
W of the heating wire can be expressed as follows:
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:
where
C represents the thermal storage unit specific heat, kJ/(kg × °C);
ρs represents the thermal storage unit density, kg/m
3;
V represents the heat storage unit volume, m
3 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:
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):
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];
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:
where
ε is the thermal conductivity at the average temperature of the air, kW/(m × °C);
μ is the kinematic viscosity of air at average temperature, m
2/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):
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:
where
A1 represents cross-sectional area of heat exchange tubes, m
2;
Saf1 is the air mass flow at the outlet of heat exchanger, kg/h; which can be expressed as:
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):
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:
where [
21,
22]:
where
λ represents friction coefficient in channel;
ρ2 represents medium density, kg/m
3;
ω 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 S
af (m
3/h) required for the fan-circulating system can be referred to Equation (16):
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:
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:
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]:
where
T7 represents heating element temperature, °C;
β represents radiation coefficient,
F1 represents radiation surface area of heating element, m
2.
For the convective heat
Qea between resistance element and air we can refer to Equation (20):
where
γ represents convective heat transfer coefficient,
F2 represents heating wire area, m
2.
The thermal
Qa carried by air can be written as:
where
qm,h represents the mass velocity of the air, m/s;
Vp,h represents the constant pressure specific volume of the air, m
3. Applying Equations (19)–(21) into Equation (18) gives:
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:
The thermal released
Qr also can be written as Equation (24) [
24]:
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):
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:
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:
The momentum conservation equation is:
The energy conservation equation is:
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:
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:
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 m
2. 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.