The Inﬂuence of the Shape of Granite on the Heat Storage Process in a Rock Bed

: Granite is one of those materials that due to its thermal parameters is used as a ﬁlling for storage beds, including high-temperature ones. The article analyzes local material that was extracted in Strzegom, Poland. The purpose of the paper is the assessment of storage material with regard to its cooperation with a heat source that is available for a short time, e.g., a solar installation. Three di ﬀ erent shapes of granite material were tested: rock, cube and sphere. Each shape has its advantages and disadvantages, which are associated with economic and strength aspects. The article presents experimental tests of the material, which were conducted in order to determine the e ﬃ ciency of the charging process. The results show that rock-shaped granite ﬁlling elements are characterized with the best parameters during the charging process, and that they obtained the highest ﬁrst- and second-law e ﬃ ciency in the entire tested range of inlet air temperature and ﬂow rate. The e ﬃ ciency of the cube-shaped granite was lower than the sphere-shaped granite. This means that the e ﬃ ciency does not directly depend on the coe ﬃ cient of sphericity of the elements that ﬁll the storage bed. The determination of the second law e ﬃ ciency showed that the highest use of energy supplied with hot air occurs after 1 h of charging the accumulator in the case of all the analyzed geometries. At the end of the paper, the inﬂuence of the obtained results on the process of modelling the charging of a storage bed ﬁlled with elements of non-spherical geometry is also discussed.


Introduction
Accumulation of heat is nowadays becoming more and more important due to the development of systems that use renewable energy sources, and also due to treatments that aim to improve the energy efficiency of enterprises e.g., by waste heat recovery. Among the various accumulation techniques, the use of materials such as rocks [1], preferably those that are locally available, is characterized by its simplicity, as well as many other advantages. It is worth noting that this material adapts well to a heat source and operates stably within a wide temperature range. This is an advantage when compared to accumulators that use a phase change material. It is especially important in cooperation with many renewable energy sources installations. Such a source of heat can be, e.g., a flat [2] or concentrating [3] air solar collector, which, due to the weather conditions in the climatic zone of Central Europe, works very irregularly (apart from in the summer period). The rock storage beds that are described in literature cooperate with both low-temperature [4] and high-temperature installations [5]. The chemical stability of this material and its resistance to corrosion are also important, as they affect the operational safety of rock storage beds [6].
The choice of storage bed filling material plays a significant role. In addition to the assessment of parameters such as the specific heat of the material, other parameters that affect the accumulation process must also be evaluated. The purpose of the storage bed also plays an important role in this process [7]. The aim of this work is the assessment of three different shapes of storage material i.e., locally available Strzegom granite. The storage bed is meant to be used for permanent cooperation with a heat source with variable time availability, such as a solar concentrating collector that uses direct solar radiation. The research of other scientists often refers to the analyses and studies concerning filling materials in the form of spheres [1,8]. Yang B et al. dealt with ceramic spheres randomly poured into a cylindrical tank [9].
The maximum temperature during the tests was close to 600 • C and the material diameter was equal to 25 mm [9]. Yang J et al. in [10] analyzed the thermal performance of a storage unit with steel particles (diameters of 6 and 12 mm). In paper [11], Yang J et al. studied three shapes: spherical, flat ellipsoid and long ellipsoid. The results show that the Nusselt number was almost equal for the spheres and flat ellipsoid and was higher for the long ellipsoid. The overall heat transfer performance was the highest for the long ellipsoid. Materials with an irregular shape, such as rocks intended for solar power plants [11], as well as those with a cylindrical shape [12], and large materials such as bricks [13,14], were also investigated. The shape influences the heat storage process and the drops of pressure in a storage bed, which Allen et al. examined in [7] while investigating the impact of eight different material shapes on the drop of pressure in the system.
The efficiency of the heat accumulation process plays an important role in the assessment of the material of the storage bed. It can be determined by, e.g., the thermal efficiency of the charging process [15]. Such optimizations also use the minimization of entropy generation [16,17] and exergy efficiency [5], as well as the second-law efficiency [15,17]. The literature lacks a broader analysis of the influence of shape on the heat transfer process, and the performed analyses are carried out using various optimization methods. It is therefore difficult to compare the results of various researchers.
In the further part of the paper, the tested bed filling material is described in detail, along with the test stand and the experiment plan. Afterwards, the calculation algorithm for determining the optimization parameters is presented. The third chapter contains the results of the experimental tests and calculations. The last chapters are devoted to a discussion regarding the obtained results and the drawn conclusions from the conducted research.

Materials and Methods
Strzegom granite, which has high technical properties (compressive strength of 110-170 MPa [18]), high durability, and is locally available in Poland, was used in the tests. The specific heat of this material is equal to 780 J/(kg·K), which makes it attractive in terms of heat storage. It is also resistant to working in high temperatures of even several hundred degrees Celsius [19].

Filling of the Storage Bed
The three shapes of the tested filling material are shown in Figure 1. The cheapest of the considered types of material is rock granite, while the production of spheres is the most expensive. Granite cubes are characterized by the greatest stability during installation and operation, while rock and sphere granites require a special stabilizing structure. Granite cubes have the largest contact surface, followed by rock granite, with sphere granite having the smallest. This affects the heat transfer between the filling elements.
For each type of storage bed filling, a high degree of similarity regarding the arrangement of the material inside the heat accumulator was ensured. The material was selected carefully, and 230 elements with a unit weight of approximately 0.174 kg were selected for each shape. The total weight of the material in the storage bed was equal to 40.02 kg. The same mass and density allow for the assumption that the characteristic dimension, which is in the form of the substitute diameter for individual elements, was the same. The way of arranging the material inside ensured the same height of the stack, as well as the same bed filling factor of 0.65. An experiment plan was designed for the material and the tests were carried out on an experimental stand. For each type of storage bed filling, a high degree of similarity regarding the arrangement of the material inside the heat accumulator was ensured. The material was selected carefully, and 230 elements with a unit weight of approximately 0.174 kg were selected for each shape. The total weight of the material in the storage bed was equal to 40.02 kg. The same mass and density allow for the assumption that the characteristic dimension, which is in the form of the substitute diameter for individual elements, was the same. The way of arranging the material inside ensured the same height of the stack, as well as the same bed filling factor of 0.65. An experiment plan was designed for the material and the tests were carried out on an experimental stand. Figure 2 shows a diagram and photo of the stand that was used for the experimental tests of the heat storage process in the granite. The used equipment enables the air temperature at the inlet to the storage bed and the flow rate of the air sucked in from the environment to be adjusted. The measured parameters include air temperature at the inlet to the storage bed, air outlet temperature, ambient temperature, and also air flow rate.   Figure 2 shows a diagram and photo of the stand that was used for the experimental tests of the heat storage process in the granite. The used equipment enables the air temperature at the inlet to the storage bed and the flow rate of the air sucked in from the environment to be adjusted. The measured parameters include air temperature at the inlet to the storage bed, air outlet temperature, ambient temperature, and also air flow rate. For each type of storage bed filling, a high degree of similarity regarding the arrangement of the material inside the heat accumulator was ensured. The material was selected carefully, and 230 elements with a unit weight of approximately 0.174 kg were selected for each shape. The total weight of the material in the storage bed was equal to 40.02 kg. The same mass and density allow for the assumption that the characteristic dimension, which is in the form of the substitute diameter for individual elements, was the same. The way of arranging the material inside ensured the same height of the stack, as well as the same bed filling factor of 0.65. An experiment plan was designed for the material and the tests were carried out on an experimental stand. Figure 2 shows a diagram and photo of the stand that was used for the experimental tests of the heat storage process in the granite. The used equipment enables the air temperature at the inlet to the storage bed and the flow rate of the air sucked in from the environment to be adjusted. The measured parameters include air temperature at the inlet to the storage bed, air outlet temperature, ambient temperature, and also air flow rate.  The experimental stand is made of a metal sheet casing insulated with mineral wool. Inside, there is a chamber that can be filled with the tested material. Air is the working medium that heats the storage material. Details concerning the test stand are provided in Table 1.

Test Stand and Experiment Plan
Due to the multitude of results that are needed to perform a broad analysis of the operation of a storage bed with various fillings, it was decided to use an experiment plan with a shortened testing time. The input factors, which are important for the heat storage process, include the air temperature at the inlet to the bed T air_inlet , and its volumetric flow rate . V air . The design of the Hartley experiment, which was described in detail in [19,20], was selected for the investigation. It enabled the tests to be carried out for three input factors at three levels of variation. The sphericity of the examined granite elements, which was marked as Ψ and determined on the basis of [21], was defined as the third factor. Parameter Ψ for the rock granite was equal to 0.6, for the cube granite to 0.8, and for the sphere granite to 1.0. Table 2 summarizes the input factors with three levels of variability, which were determined according to the principles described in [20].   Table 3 shows the matrix of the Hartley design for the performed experiment. In accordance with these settings, the test stand arrangements and bed filling material were prepared for each of the 11 experiments, and then the tests were carried out. The use of the experimental plan during the tests enables the final equation for the selected parameter (characterizing the heat accumulation process) to be determined. This facilitates the analysis of the selected shape of the material in the given temperature and flow ranges. During the experimental tests for the given autotransformer settings, the volumetric flow rate and the level of the inlet temperature to the bed were monitored. Data from thermocouples were recorded at one-minute increments.
In order to carry out comparative analyses for the indicated three shapes of granite, it was required to first calculate the average temperature of the material inside the bed. For this purpose, it was necessary to perform the thermal balance of the accumulator, which is described in the next section.

Thermal Balance of the Storage Bed
The temperature of the storage bed's filling material is used in the calculations of entropy generation. This temperature cannot be measured due to the temperature gradient within the bed, which is a result of the heat exchanged along the accumulator from the inlet to the outlet, and the fact that air flows unevenly through the bed. Therefore, for further calculations, it was necessary to determine the average temperature of the accumulator's filling material. Knowing the parameters of the air at the inlet and outlet of the bed, and the material parameters of the casing, it was possible to make a thermal balance of the amount of heat supplied with hot air . Q input and the heat losses from the surface of the accumulator . Q loss_all . This allowed the amount of stored energy . Q accu (1) and the average temperature of the storage bed's filling material T mat to be determined from (2). .
The flux of the supplied heat . Q input is described by Equation (3).
where the specific heat of air at a constant pressure c p_air is calculated as in the case of a semi-perfect gas for which the temperature is the average value T air_avg of the inlet and outlet air temperatures. The mass flow rate . m air is derived from (4) while assuming that the air is treated as a semi-perfect gas, where density ρ air is a function of the inlet temperature T air_inlet . This is because the volumetric flow rate . V air is measured at the air inlet of the accumulator. .
The tests were carried out in a laboratory, and therefore it is assumed that there was convective heat exchange between the walls and the environment outside the accumulator. The occurrence of radiation heat exchange was also taken into account. The accumulator was placed on wheels that are used for transport, and thus there was also a free convection on the bottom surface. Due to the fact that the level of filling in the accumulator was equal to 0.65, the speed of the air flowing through the accumulator was small. As a result, the performed heat loss calculations for the forced convection inside the accumulator gave a lower value than for the free convection. Due to the fact that a higher value is always considered in loss calculations, only the algorithm for determining heat losses is presented below. The algorithm also uses the equations for free convection on the inner side of each wall of the storage bed.
In the heat loss calculations that take into account the occurrence of free convection on both sides of the accumulator's casing, there are unknowns T in , Gr in , Nu in , h c_in , T out , Gr out , Nu out , h c_out , h r_out , and . Q loss . They are determined for each surface of the accumulator. This means that there is a need to determine 30 unknowns that describe the process of the heat losses from the accumulator. For this purpose, a set of equations was created for the process of heat transfer to the environment through all the surfaces of the accumulator.
The Grashof number is determined from Equation (6).
In the heat exchange inside the accumulator, the characteristic dimension is: for the side of the bed, it is its height b; and for the top and bottom surfaces, it is c h [22], which is given by (7).
where a is the width and length of the bed. The temperature difference ∆T inside the accumulator: for the side surface is T air_avg − T in_side , for the top surface is T air_avg − T in_top , and for the bottom surface is T air_avg − T in_bot . The linear coefficient of expansion β is calculated for the average value of the temperature of the layer adjacent to the wall for each of the surfaces of the accumulator. It is similar to the kinematic viscosity ν, which is a function of the average temperature of the layer adjacent to the wall for the assessed surface. The Nusselt number Nu in is determined from Equation (8) [22].
Pr wall 1 4 , where the Prandtl number Pr avg is determined for the average fluid temperature T air_avg , and Pr wall is determined for the wall temperature T in . The convective heat transfer coefficient is described by Equation (9). For the top surface, it is increased by 30%, and for the bottom surface, it is decreased by 30%.
For the outer side surface of the accumulator, the characteristic dimension is the height that is equal to b + δ ins_top + δ ins_bot , while for the top and bottom surfaces it is d h , which is expressed by (10).
The Nusselt number Nu out for the external surfaces is described by general Equation (11): where the constants are as follows: The Grashof and Prandtl numbers are determined analogously to those determined for the inside of the accumulator. Only the characteristic dimensions and the selected temperature differences are changed. The Prandtl number and the physical properties that occur in the Grashof number are determined for the average temperature of the layer adjacent to the wall. Moreover, radiative heat exchange occurs outside the accumulator and is described by the radiative heat transfer coefficient h r_out , which is determined from Equation (12).
From Equation (6), the Grashof number for the heat transfer inside and outside the accumulator was determined. The Nusselt numbers were determined from Equations (8) and (11), and the convective heat transfer coefficient for both sides of the partition were determined from Equation (9). The radiative coefficient of heat transfer to the outside of the accumulator was calculated from Equation (12). Therefore, seven equations describing the process of heat transfer through the partition are obtained. The last three equations for determining all ten unknowns for each of the surfaces are general equations for the heat transfer through the surface of the accumulator . Q loss . The first Equation (13) describes the heat transfer on the inside of the accumulator. .
The heat flux that enters the partition is equal to the heat flux received on the other side, which for the analyzed accumulator is described by Equation (13).
The third equation takes into account the heat flux conducted through the wall. It is equal to the heat flux determined from Equations (13) and (14) and is expressed by formula (15).
In Equation (15), the surface area is the average value of the partition's surface area. Due to the fact that it varies linearly from A in inside the accumulator to A out outside the accumulator, it can be written as (16): The ten equations described above allowed for the determination of ten unknowns that describe the heat exchange process between the air flowing inside the accumulator and the environment for one surface. The total heat losses from the bed . Q loss_all , which are described by Equation (5), take into account three surfaces. Therefore, these equations were written for each of the surfaces, resulting in a system of 30 equations. In order to determine the searched average bed temperature, Equations (1)-(3) were added to the thermal balance, and the system of 33 equations and 33 numerically determined unknowns was created. The average granite temperature, determined in this way, was used in the optimization analysis, which is discussed in the next part of the paper.

Calculation of the First-Law and Second-Law Efficiencies
The collected experimental data and the results of the heat balance allowed the thermal efficiency of the process and the entropy generation to be specified in order to determine the second-law efficiency.
There are many definitions of the first law of thermodynamics, also known as the first-law efficiency [23]. However, they mean the same and result from the energy balance. Generally, the first-law efficiency can be defined as the quotient of the energy effect obtained for the system in question, and the energy supplied to this system. It is written as (17).
Energy e f f ect Energy supplied (17) The first-law efficiency for heat accumulators was written, e.g., by Bejan [17] as (18).
For the analyzed case, it can be written as (19): Then, the thermal efficiency of the heat accumulator's charging process, also called the first-law efficiency, was derived in [15] and presented as (20).
The introduced dimensionless time θ is expressed by (21).
m air ·c p_air n·m mat ·c p_mat ·t The heat exchanger's parameter y is expressed as (22).
The number of heat transfer units N tu is determined from dependencies (23) [17].
Calculations related to the concept of entropy and exergy were carried out. Cengel and Boles [23] define the second-law efficiency as a measure of approximation to a reversible process. Therefore, such efficiency for any system can be written as (24).
According to the investigation conducted by Bejan [17], it can be expressed by (25). For the case of the storage processes, Equation (26) was used to estimate the value of the entropy generation rate . S gen [17].
. S gen = . m air ·c p_air ·ln T amb T air_inlet + . Q out T amb + n·m mat ·c mat · dT mat dt (26) where . Q 0ut is the heat transfer to the environment, which was calculated from (27).
Moreover, the total supplied exergy was estimated using Equation (29) on the basis of [17].
m air ·c p_air ·t· T air_inlet − T amb − T amb ·ln T air_inlet T amb The results of measurements and calculations are presented graphically in the next chapter.

Results
The analysis of the influence of the shape of granite on the charging process of the heat accumulator was based on the well-known thermal efficiency of the process. The change in efficiency over time was shown for different temperatures of air entering the accumulator and for different flow rates in order to determine the influence of these parameters on the analyzed geometries of the granite. Thermal efficiency does not define the charging time of an accumulator, and therefore it is worth checking the second-law efficiency, which provides information regarding this matter.

Experimental Research Results
The ambient temperature was around 22 • C during all the experiments. Figure 3 shows the results of the measurements of the air flowing out of the bed, with the measuring points being designated every 10 min for better orientation. The various inlet temperatures are marked with colors, while the shape of the filling material and the volumetric air flow rate are marked with symbols.

Experimental Research Results
The ambient temperature was around 22 °C during all the experiments. Figure 3 shows the results of the measurements of the air flowing out of the bed, with the measuring points being designated every 10 min for better orientation. The various inlet temperatures are marked with colors, while the shape of the filling material and the volumetric air flow rate are marked with symbols. The temperature characteristics shown in the graph tend to flatten. Such flattening took place almost completely in the subsequent hours of the charging process. Due to the purpose of the comparative research carried out in the article, only the first three hours of the process were analyzed, which are very important with regard to the cooperation between the rock bed and the heat source, which is available for a short time. In the first minutes of charging, the curves for rock granite are almost linear in all the tested cases. The lowest achieved outlet temperature may be due to the intense heat transfer at this time. The less streamlined shape of the rock-shaped granite can cause a more turbulent flow, while at the same time increase the heat transfer.
During the tests, drops in pressure in the bed were also measured, and in the tested range they were not dependent on the shape of the material, but instead only on the volumetric flow rate. For the lowest flow of 0.006 m 3 /s, the pressure drop in the bed was equal to 3 Pa; for 0.008 m 3 /s, the pressure drop in the bed was equal to 5 Pa; and for 0.010 m 3 /s, the pressure drop was equal to 7 Pa. This means that the amount of energy lost to overcome hydraulic resistances is many times smaller than the transferred thermal energy. Therefore, further analyses omit the influence of hydraulic drops on the second-law efficiency.

Thermal Efficiency of the Charging Process
The used Hartley design serves as a means to generate a quadratic equation describing the unknown parameter yexp with the use of three independent variables. In the analyzed case, a series of eleven experiments, performed in accordance with the design, allow the following polynomial to be The temperature characteristics shown in the graph tend to flatten. Such flattening took place almost completely in the subsequent hours of the charging process. Due to the purpose of the comparative research carried out in the article, only the first three hours of the process were analyzed, which are very important with regard to the cooperation between the rock bed and the heat source, which is available for a short time. In the first minutes of charging, the curves for rock granite are almost linear in all the tested cases. The lowest achieved outlet temperature may be due to the intense heat transfer at this time. The less streamlined shape of the rock-shaped granite can cause a more turbulent flow, while at the same time increase the heat transfer.
During the tests, drops in pressure in the bed were also measured, and in the tested range they were not dependent on the shape of the material, but instead only on the volumetric flow rate. For the lowest flow of 0.006 m 3 /s, the pressure drop in the bed was equal to 3 Pa; for 0.008 m 3 /s, the pressure drop in the bed was equal to 5 Pa; and for 0.010 m 3 /s, the pressure drop was equal to 7 Pa. This means that the amount of energy lost to overcome hydraulic resistances is many times smaller than the transferred thermal energy. Therefore, further analyses omit the influence of hydraulic drops on the second-law efficiency.

Thermal Efficiency of the Charging Process
The used Hartley design serves as a means to generate a quadratic equation describing the unknown parameter y exp with the use of three independent variables. In the analyzed case, a series of eleven experiments, performed in accordance with the design, allow the following polynomial to be obtained (30).
where b 0 , b 1 , . . . , b 10 are regression coefficients. Due to the use of the experiment plan, the equations were determined. Thanks to this, the spatial characteristics of the change in the given parameters within the range of the assumed boundary conditions, which are defined in Table 2, were plotted. In the first step, Equation (30) was used to determine the first-law efficiency, which resulted in the creation of the 3D graphs as shown in Figure 4. The characteristics shown in Figure 4 present the change in thermal efficiency over time within the examined temperature range and mass flow rate.
Energies 2020, 13, x FOR PEER REVIEW  10 of 16 where b0, b1, …, b10 are regression coefficients. Due to the use of the experiment plan, the equations were determined. Thanks to this, the spatial characteristics of the change in the given parameters within the range of the assumed boundary conditions, which are defined in Table 2, were plotted. In the first step, Equation (30) was used to determine the first-law efficiency, which resulted in the creation of the 3D graphs as shown in Figure 4. The characteristics shown in Figure 4 present the change in thermal efficiency over time within the examined temperature range and mass flow rate. As can be seen from Figure 4, the highest thermal efficiency in all the analyzed cases is achieved for rock granite. Lower values are obtained for sphere granite, and the lowest values for cube granite. The biggest differences are between the 60th and the 100th minute of charging and are especially visible for the highest flow and lowest inlet temperature. However, the influence of temperature and flow rate on the thermal efficiency for each of the analyzed geometries is not large. The charging time, i.e., the changing temperature of the air and granite, has a much greater influence. Therefore, it is worth analyzing the accumulator's charging process with regard to the second-law efficiency.

The Generation of Entropy-The Second-Law Efficiency
By performing calculations of the generation of entropy, the total exergy content of the air drawn from the hot air supply, and by applying the experiment plan, it was possible to plot the change in the second-law efficiency over time for the examined area of the air inlet temperature and the flow rate ( Figure 5). The calculations were also conducted with the use of the method described in Section 3.2 and Equation (30), but this time in order to determine the second-law efficiency.
The determined second-law efficiency achieves the highest values for rock granite in almost the entire tested range, and the lowest values for cube granite. As can be seen from Figure 5, the efficiency of the charging process for rock granite reaches the highest values after 60 min for the highest inlet air temperature. For the sphere and cube granites, the second-law efficiency also reaches the highest values after 1 h. However, after this time, the efficiency for rock granite drops very quickly, and reaches the lowest values at the end of the measurements for the entire tested range of the flow rate. For the charging time shorter than 2 h, the use of the energy potential for the cube and sphere granites is smaller than for the rock granite. For these geometries, the second-law efficiency decreases slower over time. As can be seen from Figure 4, the highest thermal efficiency in all the analyzed cases is achieved for rock granite. Lower values are obtained for sphere granite, and the lowest values for cube granite. The biggest differences are between the 60th and the 100th minute of charging and are especially visible for the highest flow and lowest inlet temperature. However, the influence of temperature and flow rate on the thermal efficiency for each of the analyzed geometries is not large. The charging time, i.e., the changing temperature of the air and granite, has a much greater influence. Therefore, it is worth analyzing the accumulator's charging process with regard to the second-law efficiency.

The Generation of Entropy-The Second-Law Efficiency
By performing calculations of the generation of entropy, the total exergy content of the air drawn from the hot air supply, and by applying the experiment plan, it was possible to plot the change in the second-law efficiency over time for the examined area of the air inlet temperature and the flow rate ( Figure 5). The calculations were also conducted with the use of the method described in Section 3.2 and Equation (30), but this time in order to determine the second-law efficiency.

The Generation of Entropy-The Second-Law Efficiency
By performing calculations of the generation of entropy, the total exergy content of the air drawn from the hot air supply, and by applying the experiment plan, it was possible to plot the change in the second-law efficiency over time for the examined area of the air inlet temperature and the flow rate ( Figure 5). The calculations were also conducted with the use of the method described in Section 3.2 and Equation (30), but this time in order to determine the second-law efficiency.
The determined second-law efficiency achieves the highest values for rock granite in almost the entire tested range, and the lowest values for cube granite. As can be seen from Figure 5, the efficiency of the charging process for rock granite reaches the highest values after 60 min for the highest inlet air temperature. For the sphere and cube granites, the second-law efficiency also reaches the highest values after 1 h. However, after this time, the efficiency for rock granite drops very quickly, and reaches the lowest values at the end of the measurements for the entire tested range of the flow rate. For the charging time shorter than 2 h, the use of the energy potential for the cube and sphere granites is smaller than for the rock granite. For these geometries, the second-law efficiency decreases slower over time.

Discussion and Conclusions
The literature review showed a negligible amount of papers concerning the influence of the geometry of heat accumulator filling elements on the efficiency of the storage process. This especially applies to materials with a non-spherical geometry. Therefore, the authors prepared appropriate research materials and performed experiments based on the experimental plan. The results of different outlet air temperatures were obtained for different granite geometries and the same initial conditions. It was shown that the geometry of the accumulator's filling elements affects the thermal efficiency of the charging process. In the entire tested range, the highest efficiency was obtained for rock granite, which is characterized by the lowest coefficient of sphericity. The process of charging an accumulator filled with sphere granite has a lower efficiency, and the lowest efficiency is obtained for an accumulator filled with cube granite. This means that the coefficient of sphericity is not a parameter that determines the charging efficiency.
The determination of the second-law efficiency enabled the degree of using the energy supplied with hot air to be determined. The change in this efficiency over time showed that the best use of the supplied energy is obtained after 1 h. The efficiency then reaches the highest values for all the tested geometries. However, the charging process is most intensive for rock granite and the least intensive The determined second-law efficiency achieves the highest values for rock granite in almost the entire tested range, and the lowest values for cube granite. As can be seen from Figure 5, the efficiency of the charging process for rock granite reaches the highest values after 60 min for the highest inlet air temperature. For the sphere and cube granites, the second-law efficiency also reaches the highest values after 1 h. However, after this time, the efficiency for rock granite drops very quickly, and reaches the lowest values at the end of the measurements for the entire tested range of the flow rate. For the charging time shorter than 2 h, the use of the energy potential for the cube and sphere granites is smaller than for the rock granite. For these geometries, the second-law efficiency decreases slower over time.

Discussion and Conclusions
The literature review showed a negligible amount of papers concerning the influence of the geometry of heat accumulator filling elements on the efficiency of the storage process. This especially applies to materials with a non-spherical geometry. Therefore, the authors prepared appropriate research materials and performed experiments based on the experimental plan. The results of different outlet air temperatures were obtained for different granite geometries and the same initial conditions. It was shown that the geometry of the accumulator's filling elements affects the thermal efficiency of the charging process. In the entire tested range, the highest efficiency was obtained for rock granite, which is characterized by the lowest coefficient of sphericity. The process of charging an accumulator filled with sphere granite has a lower efficiency, and the lowest efficiency is obtained for an accumulator filled with cube granite. This means that the coefficient of sphericity is not a parameter that determines the charging efficiency.
The determination of the second-law efficiency enabled the degree of using the energy supplied with hot air to be determined. The change in this efficiency over time showed that the best use of the supplied energy is obtained after 1 h. The efficiency then reaches the highest values for all the tested geometries. However, the charging process is most intensive for rock granite and the least intensive for cube granite. As a result, rock granite, apart from its attractive price, also has the best geometry with regard to thermal analysis due to its non-spherical surface.
Finally, it is worth noting that during the model calculations, the flux of the heat transferred between hot air and the material filling the accumulator . Q accu can be determined from Equation (31). This means that this flux for the same starting temperatures of air and granite depends only on the heat transfer coefficient h and the heat transfer surface area A mat .
Literature provides formulas for the Nusselt number for the determination of the heat transfer coefficient, which were developed by Ranz and Marchall [24], Kacnelson and Timofiejewa [25], Ketterning [26], and Delton [27]. The characteristic dimension is the diameter of the material, or the equivalent diameter, which is calculated from the volume of the sphere. This means that when modelling an accumulator filled with elements with the investigated geometry, the heat transfer flux . Q accu should assume the same values. However, as the research shows, the differences between the various geometries are visible. This is mainly due to the different heat exchange surface, the contact surface between the elements of the accumulator's filling, or the different flow around non-spherical bodies. Therefore, according to the authors, the modelling of the heat storage in materials with complex geometry requires deeper analysis and should be subjected to further research.