Utilization of Buildings’ Foundations for a Seasonal Thermal Energy Storage Medium to Meet Space and Water Heat Demands

Seasonal thermal energy storage (STES) can be utilized to cover a portion or meet the whole space and water heat demands in residential and commercial buildings. With the scarcity of fuel resources, scientists have become aware of the importance of the utilization of thermal energy. An STES system can be charged using solar collectors to heat a storage medium when solar radiation is available. Solar irradiance is seasonal so this type of system can compensate for the shortage of energy in the winter by storing surplus solar energy in the summer. This study offers an innovative design of an STES system that takes advantage of the backfill space in a building’s underground foundation to install the thermal storage medium on the excavated surface among the support pillars’ base. To study the feasibility and the feasibility of such STES system, the design was simulated using the TRNSYS® tool where the STES medium was charged by a solar thermal system. The chosen building type was a four-story hotel located in three cities each with its unique climate: Potsdam (Germany), Zarqa (Jordan), and Doha (Qatar). The results showed the coverage rate of the building’s annual heat demand was 56%, 82%, and 84%, and the payback periods were 3.4, 4.4, and 29 years for Potsdam, Zarqa, and Doha, respectively. Zarqa is the most efficient and feasible for STES compared to the other two cities.


Introduction
Buildings are one of the biggest energy consumers with 32% of the total final energy consumption in the world according to the International Energy Agency. The fossil energy systems can be replaced by solar energy systems especially for locations with high annual solar radiation, which lead to prevent large greenhouse gas emissions and reduce the negative impact on the environment as well. Solar energy, while clean and renewable, is diurnal and seasonal in nature, hence, the importance of implementing techniques to store this energy for later use when solar radiation is low or entirely unavailable. Seasonal thermal energy storage (STES) is implemented widely by using different technologies. Depending on the STES medium, surplus thermal energy from the summer can be stored and then used in the winter at various scales ranging from single-family homes, hotels, districts, towns, or regions.
Humans have been searching for means to store energy; therefore, energy storage is considered an old practice. In ancient Greece and Romans times, architects used solar energy inside their buildings to heat water. They designed their bathhouses' openings to acquire solar heating to support their existing bath fire, thus reducing the fuel needed for water heating loads [1]. Native Americans inhabited caves that best utilize solar energy in winter and summer at the Grand Canyon of Arizona. The low azimuth angle of the sun in the winter provided heat while the shading of the natural cliff blocked the summers' high azimuth angle.
Nowadays, there is a continuous increase in energy demand globally. Climate change significantly affects cooling and heating loads as well as the annual energy demands of for a building [2]. Solar thermal energy could be stored for months and seasons solar thermal systems' efficiency, sustainability, and productivity can be improved by implementing an STES criteria. Solar thermal energy can be stored by changing the internal energy of different materials in different phases. Materials that could be utilized are sensible, latent, and thermochemical as a medium of an STES. Furthermore, materials could be combined to achieve the required thermal properties. Due to the rising temperature of the thermal storage medium, the atoms and molecules of its materials vibrate faster and collide with each other to release the thermal energy. This allows the thermal energy to be stored as a change of the internal energy of the material in many forms. Sensible heat, latent heat, and thermochemical heat are forms of energy storage; a combination of these forms is also possible for storage [3].
The sun provides enough energy to meet global demand; however, many variables affect the accumulated energy from the sun. Therefore, some sort of energy storage system must be introduced to maximize the benefit of solar energy. Solar thermal systems use heat collectors to capture solar energy to be utilized according to the demand load. Energy can be stored for short-term and extended periods (seasonal) according to the storage system's size. Seasonal energy storage is evaluated to be 100 to 1000 times larger in capacity than short-term storage size [4].
A reference approach of STES was designed, modeled, and simulated at four buildings in different cities under various climatic zones in Europe (Stockholm, Zurich, Barcelona, and Madrid). A solar thermal system was linked with the STES to discharge the thermal energy to the buildings' heating applications. The simulation was implemented by using TRNSYS® software. The authors of the study found that a STES with 1 m 3 volume that was charged by 20 m 2 flat plate collectors and supported by 10 kW auxiliary heater can cover the buildings' annual heating demand [5].
The performance of underground fabricated sand bed STES in six different typical single-family homes' models in Richmond, Virginia, USA, was tested by using TRNSYS®, where the simulation was performed to calculate the stored thermal energy in 15 m 3 STES that gained from flat plate solar thermal collectors to cover the home's space heat demand. As a result, the overall system efficiency ranged from 50 to 70% compared with total useful solar energy. And the maximum temperature of the STES reached in the 5-year simulation was 96.1°C [6], while other researchers studied the performance of STES under the same simulation's parameters for a five-story student housing complex at Virginia Commonwealth University. The simulation was performed by using TRNSYS® on 15 thermal energy storage locations to make optimization to choose the most efficient STES system. As a result, the STES can cover more than 91% of the building's space heat demand after the fifth year of its operation [7].
Multiple authors investigated the use of solid materials as a sensible heat storage medium. A study of storing of the thermal solar energy in the soil in the summer season in Al-kreaat quarter in Baghdad used the quasi-threedimensional model and finite line source model in their computations. The results found that 38% of the total collected energy can be stored in the soil which indicated that the amount of stored energy depends on the borehole length [8]. An investigation on covering the heat demand of four buildings within a complex by using STES placed in the ground. The capacity of STES was 5,700 m 3 that was charged by a solar thermal system with collectors' area 4,700 m 2 . Several scenarios were analyzed based on the buildings' heat demand to optimize the best size of the STES, and as a result, the optimized STES could give around 53% solar fraction factor where that could decrease heat consumption by 22-100% throughout the year [9].
A TRNSYS® simulation was performed on a typical STES system for a typical 120 m 2 building in Thessaloniki, Greece. The simulation performed to calculate the covered building's space and water heat demand by 36 m 3 hot water STES. The operation of the simulated STES system was examined for two years. As a result, the hot water STES covered more than 52.3% of the building's heat demand, and the maximum temperature of the STES reached in the second year simulation was in the range 84-90°C [10].
Sand-bed thermal energy storage (TES) system experiment was conducted in a region with an extended freezing period at a two-story house located in Palmer, Alaska, using TRANSYS® simulation, and the simulated temperature of the sand-bed was compared to the measured one with a relative error of 15%. The results indicate that such types of solar thermal storage systems are viable options for very cold climates [11]. Conversely, in the United Arab Emirates (UAE), a study demonstrated that the desert sand at UAE could be used in concentrated solar power (CSP) facilities to store thermal energy up to 1000°C. The research project, Sand Stock, used sand particles as the heat collector, heat transfer, and TES media. The study concluded that desert sand can be used not only as a TES material but also as a direct solar absorber under concentrated solar flux [12].
The feasibility of STES using sensible, latent, and thermochemical heat storage technologies through a residential sector in eight cities in the United Kingdom (UK) was investigated, and a calculation modeling framework was used to estimate the critical solar collector areas and the required capacity of a TES [13].
In the case of using solid material as a medium of STES, a study examined the performance of a novel combined receiver storage system for beam down solar concentrator and tested different mixing ratios of sand and basalt to be used as the thermal storage medium. The finite element heat transfer model was used to solve the unsteady governing equations for the system. Accordingly, the achieved results showed that the utilization factor is in the range of 0.53-0.73 depending on the conductive rods configuration and the sand-basalt mixture ratio. Furthermore, the utilization factor increases by increasing the basalt percentage [14]. Others experimented with the utilization of the basalt stone as a material in sensible heat storage and found that it has a high energy density of 4:2 × 103 kJ m −3 K −1 [15].
The performance of underground water STES in a cold climatic zone in Ottawa, Canada, was tested where the storage volume was 36 m 3 and charged by 41.6 m 2 evacuated tube collectors for covering the heat demand of a single-family home. The simulation tools used were TRNSYS® and ESPr; as a result, the covered home's heat demand was in the range of 87-98% using this system [16].
A study assessed using STES instead of using an air-source heat pump to cover the space heat demand of a 240 m 2 house in Hangzhou, China. A TRNSYS simulation was performed to develop and optimize short-term solar STES to reach the 2 International Journal of Photoenergy lowest annual cost over 20 years. The electricity demand reduction was 74% by the optimized STES [17]. All researches emphasized that there is an urgent need to develop and apply smart and flexible heat systems by integrating thermal energy storage. In this case, thermal energy storage can provide high-demand response, which will need smart operating system and control techniques.
This study investigates a novel STES systems by utilizing silica sand, basalt stone, and copper mesh, in placed of the backfilled area in the foundations' layer underneath a building. Using the available spaces for the underground STES could reduce the cost of the required excavation and also providing a method for taking advantage of an otherwise underutilized area. The STES system will be simulated in three different climatic zones to assess its performance based on the variations in the climatic features of the selected locations: (i) Zarqa-Jordan, (ii) Doha-Qatar, and (iii) Potsdam-Germany. The authors chose the sensible STES type instead of latent/thermochemical because it is the safest, the most common, the simplest, and the cheapest one [1].
This research investigated the possibility and the feasibility of using such STES to meet space and water heat demands of a hotel building in the three climatic zones using TRNSYS® simulation software.

Methodology
The heat demand for space heating (SH) and domestic hot water (DHW) of the buildings will be calculated utilizing the degree days method. A proper solid material will be utilized in a sensible STES to cover the building's heat demand. This study investigates the thermal energy application of an STES system under a proper boundary and initial conditions that will be discussed in detail in the next sections. Investigation of the reliability of the designed STES system relied on mainly three energy pillars: energy storage, renewable production, and STES system sustainability [18]. The authors represented the three energy pillars in the results. Thermal energy transfers throughout the whole system are illustrated in Figure 1 as a simple schematic diagram. A simulation study will be performed for the STES system in the selected locations by using TRNSYS® software. An economical evaluation for capital and operation cost of the STES system compared with the heat demand electricity cost will be done in the nest sections. Finally, the simulation and the economical evaluation results will be analyzed to assess the STES performance and feasibility in the selected locations. The steps and procedures to be followed are presented in a chart shown in Figure 2

Design and Simulation
Three different climatic zones were selected to test the performance of the STES and its feasibility for different weather conditions, Zarqa in Jordan, Doha in Qatar, and Potsdam in Germany. The climatic and solar properties for the three cities are shown in Table 1 according to the National Climatic Data Center [19] and Photovoltaic Geographical Information System [20], based on the different climate features, including temperature, humidity, and cloud cover. The authors will assess the best performance of the STES that will be relevant to the best solar energy features for one of the selected locations. Jordan is characterized by having available solar energy along with more than 300 sunny days/year, where it is one of the sun-belt states and most of its area is exposed to a very high solar irradiance [21].

Architectural
Design of the Building. The four-story building, hotel with floor area 400 m 2 , was chosen for this study's design and simulation. The building was designed with the same dimensions and orientation (north orientation) in the three cities but with different insulation materials due to the different weather conditions. The hotel's occupation capacity is 80 persons including the staff, while the capacity of each bedroom hosts two persons. The hotel building was drawn using SketchUp® in a three-dimensional view as shown in Figure 3(a) and a two-dimensional floorplan for the ground floor and the one above it as shown in Figure 3(b). Three different buildings' envelope constructions were designed following the codes as, listed in Table 2, to reduce the heat loss/gain due to the different levels of weather conditions. The authors used Polyurethane insulation panels between outside stone and internal concrete hollow block as a wall assembly in different thicknesses based on the climatic features of each selected location [23].  International Journal of Photoenergy demand. Space heating (SH) loads of the three buildings were calculated based on the degree days method [25]. The designation of heating degree days for the three different climatic zones provides a tempo-spatial map of heating demand. The knowledge of total hotel's space heating demand was used to identify the accumulated heat demand curve, appropriate for seasonal thermal energy storage system modelling.

Hotel's Heat
The conduction and infiltration loads of the hotel buildings were calculated using equations ((1)-(3)) with 21°C indoor design temperature.
where Q c is the heat gain/loss due to conduction through walls, windows, and roof (kWh); Q i is the heat gain/loss due to infiltration (kWh); U is the overall heat transfer coefficient (W m -2 ·K -1 ); A is the surface area for walls, windows, or roofs (m 2 ); HDD is the heating degree days (K day); ρ is the density of air (1.225 kg m -3 ); C p is the specific heat of air (1000 J kg -1 K -1 ); V is the volume changed per hour based on the Air Change Method (m 3 h -1 ).
Two schedules of the hotel's internal load that generated by occupancy and electrical equipment were estimated by the authors to calculate the amount of heat gain.
The resulting annual space heating load is shown in Figure 4. It can be noticed that the annual space heating load for Potsdam was the highest at 101 kWh m -2 year -1 , which is comparable to the 100 kWh m -2 year -1 for the hotel sector in Germany [26]. On the other hand, the annual space heating load for Zarqa was 57 kWh m -2 year -1 ; this result is a bit lower than what was previously recorded at around 95 kWh m -2 year -1 [25]. This could be attributed to the fact that this study used a building envelope's construction with U values lower by half than the one that was used in [25]. Finally, the annual space heating load for Doha city was the lowest one at 14 kWh m -2 year -1 due to its warm climate throughout the year. This result is comparable to 15 kWh m -2 year -1 previously recorded for Doha [27].

3.2.2.
Hotel's Domestic Hot Water (DHW) Load. It was assumed that the percentage of the annual occupation of the hotel is 60% (around 45 residents). The daily DHW consumed per person was estimated to be equal to 30 liters person -1 day -1 at a set temperature for the hot water at 55°C according to the Standard Assessment Procedure (SAP) [28]. And the cold water temperatures were generated from NOAA website for the three cities, respectively [19]. The hotel's DHW load was calculated based on equations (4)-(5) [29].  [19,20].      International Journal of Photoenergy where Q w is the required energy for DHW load (kWh), M d is the volume of daily consumed hot water (liter day -1 ), C water is the specific heat capacity of water (1.16 kWh m -3 K -1 ), ΔT is the temperature difference between cold and hot water (K), and ρ water is the density of water (0.988 kg m -3 ). The resulting annual DHW load is shown in Figure 5(a); it is observed that it varies depending on the climatic zone, whereas the annual DHW load of Zarqa was 17,583 kWh year -1 , the lowest annual DHW load was for Doha at 12,597 kWh year -1 , and the highest load was that of Potsdam at 19,095 kWh year -1 . All DHW loads are comparable to previously recorded results [29]. Consequently, the resulting total annual heating demand of the three cities is shown below in Figure 5(b).

Solar Thermal Combisystem
3.3.1. Solar Thermal Collector. Solar thermal collectors are the main source of the thermal energy that will supply to the building's thermal applications to cover its heat demand. The solar thermal collectors charge the STES for storing the thermal energy from the summer season to be used in the winter season. Such systems that consist of solar thermal collectors charging DHW and space heating (SH) applications is called a solar thermal combisystem.
Evacuated tube collectors are more efficient than flat plate collectors, and they have a longer life cycle with lower required maintenance [30,31]. The used solar thermal collectors in this study were the Compound Parabolic Collector (CPC 18) due to their high efficiency and sustainability. The collectors were oriented to the south direction, and the simulation was performed with different collector tilts due to the incident solar irradiance in the three cities. The required area of collectors for three locations was calculated using equation (6) [32], where the solar fraction was estimated by authors to equal 0.8 [33]. The compound parabolic collectors (CPC) are typical collectors that are able to collect both beam and diffuse   [34]. The technical details of the used solar thermal collector are shown in Table 3.

A c =
No:of days × Q water × solar fraction yearly solar irradince × Avg system efficiency : Two working fluid were used to be circulated in two loops that cross the heat exchanger. Ethylene glycol has far superior heat transfer properties more than propylene glycol. Accordingly, the ethylene glycol was circulated in the primary loop between the collectors and the heat exchanger because it provides effective freeze protection as long as the proper antifreeze concentration is maintained. Ethylene glycol would be changed every 3-5 years because the antifreeze fluids degrade over time [35]. The circulated hot water will pass through the second loop driven by pump power, and a part of that circulated hot water will be stored temporarily in diurnal thermal energy storage (DTES) for daily usage. The working fluid loops were illustrated as shown in Figure 1. Consequently, the volume of the DTES was calculated by using equation (7) where HWD p is the average hot water demand per person (liters) and HWD r is the average hot water demand for the hotels' kitchens (liters) and the annual occupancy is 60% for the three cities.

Solar Thermal Heat Exchanger.
Transferring the thermal energy from the solar collectors' circulation to the STES needs an effective heat exchanger working on a liquid to solid transferring approach as shown in Figure 1. According to the heat atlas [37], the shell and tube heat exchanger was the best choice to use in this STES system. Shell and tube heat exchanger has high effectiveness due to its large available heat transfer area. It was assumed that the effectiveness of the used heat exchanger is up to 90%.

Seasonal Thermal Energy Storage (STES).
An STES system was designed to fit in the filler area above the foundations' layer underneath the building. STES medium was solid materials like silica sand and basalt stone to be utilized as sensible heat storage (SHS) combined with a copper mesh to promote the distribution and diffusion of heat throughout the STES as shown in Figure 6. The whole parameters that illustrated in Figure 6; they were inserted in the STES component's models in TRNSYS® simulation software to clarify the STES thermal variables and its boundary conditions. Three types of solid materials were investigated in this research to be utilized in the STES for storing the produced thermal energy by the solar thermal system throughout the year. The first one is the basalt stone; basalt is a very hard and dark-colored stone that is formed by volcanic activity, and it is available in many countries [38,39]. The second type was silica sand; it is a mineral resource that can be used for different applications due to its good thermal properties and availability in both Europe and Jordan [40]. The third type was copper; a copper mesh was used to distribute the thermal energy within the STES and then to transfer it to the building's thermal applications through the heat exchanger when it is required. Copper is characterized by its high hardness and its high transmissivity of thermal energy due to its high melting point and its high resistance to corrosion [41]. The thermal specifications of the STES medium's materials are presented in Table 4.
The required area of the buildings' foundations in the three cities is dependent on the type of land and live and dead  International Journal of Photoenergy loads of the buildings. Based on that, the available area of the STES can be calculated by subtracting the required area of the building foundations from the overall building area. With the high variety of that characteristic of the three cities, it was assumed that the ratio of building's foundation area is 20% of the building's floor area depending on best practice records. Therefore, the available volume under the building is 300m 3 . The overall STES-specific heat was calculated using equation (8) and in accordance with the STES materials thermal specifications listed in Table 4. The mass percentage was divided due to the highest specific heat (C p ) of the STES materials. The percentage ratios are 65% of basalt stone, 30% of silica sand, and 5% of copper mesh.
where m silica sand , m basalt , and m copper are the selected mass (kg) of the silica sand, basalt, and copper, respectively; m total is the total mass of the STES medium (kg); C silica sand , C basalt , and C copper are the specific heat (kJ kg -1 K -1 ) of the STES materials silica sand, basalt, and copper, respectively. The STES was insulated with 0.2 m of polyurethane foam and 0.2 m of concrete stratification to avoid the heat loss between the STES and the underground environment and to avoid the leakage from the surrounding groundwater as well [44].

Simulation Parameters.
The simulation was constructed based on the deck file of the task 32 [5] by following the simulation running process steps inherent in TRNSYS®. The simulation was performed for a five-year span starting from June where the STES can operate without an auxiliary heater. The TRNSYS® simulation parameters are shown in Table 5.
Based on different heat demands of the buildings, the required collector area was calculated using equation (6) The collectors' tilt angles were fixed to a certain angle according to each city's coordinates. The authors did not use a solar tracking system according to its high cost, where it needs periodic maintenance as well [45]. The initial STES temperatures were equal to the ground temperatures for the three cities, respectively, where it is 21.5°C in Zarqa [46], 28.9°C in Doha [47], and 12.7°C in Potsdam [48].
Based on the high building's heat demand and the solid materials into the available STES volume, the temperature of the STES was set to the maximum 150°C for the three cities to avoid the heat dissipation and losses. Then, the collector boiling temperature was calibrated to 155°C [49]. The STESspecific heat was calculated using equation (8).
The rate of the covered heat demand (Q HD ) by the STES for each month will be calculated by using equation (9).

Model Validation
The stored thermal energy for the long term under certain climate conditions is evaluated by the capacitance of the STES to cover the building's heat demand. Therefore, the STES capacitance is the most critical metric tool that will be taken into consideration in any STES design process. TRNSYS® simulation tool can provide an accurate illustration of the thermal performance of a STES.
In order to validate this study, another study, [16], was replicated through a TRNSYS® simulation using the same STES parameters. Accordingly, the maximum relative error for the monthly stored thermal energy between the validation and the referenced study was 1%, while the average relative error was 0.6% as shown in Figure 7.

Results
The utilization of STES allows the storage of thermal energy from the summer season to be recalled during the winter season. The design of the STES system in three different climatic zones based upon a simulation on TRNSYS® was performed, while the hotel's heat demand was calculated based on the illustrated methods and then compared with the resulted stored thermal energy from TRNSYS® simulation process. For more details about the used deck file in TRNSYS® simulation software (see Appendix (II)).    Figure 10(a) along with the percentage of seasonal and annual covered heat demand simultaneously for the three cities in Figure 10(b For more details, see Appendix (I) that contains the simulation results of the thermal energy transfer in the STES for the three locations.

Comparison with Literature: STES Pros and Cons.
The results of this study were compared with the results of the study used for the simulation's validation in reference [16] and are listed in Table 6.

Economic Evaluation
The capital cost of the solar thermal combisystem combined with STES was calculated to be compared with electricity   (7)). 13 International Journal of Photoenergy consumption for covering the building's heat demand. The calculating approach of the initial cost for each component in the solar combisystem was calculated according to reference [50], while the initial cost of the STES components was calculated though best practices due to the difficulty of its pricing in the market.
6.1. The Capital Cost of the STES System for the Three Cities. The components' initial costs are listed in Table 7 for the three cities along with the replacement period of each component. The life span relevant to each component is different rendering it difficult to have a fixed life span for the whole system.
6.2. The Annual Electricity Cost and Savings for the Three Cities. The electricity consumption for covering the building's heat demand was calculated using equation (10) for the three cities in order to determine the electricity cost for each city; the results are shown in Table 8.
The savings from the electricity bills were calculated using equation (11) for the three cities.
Annual Electricity Savings $ kWh ! = Q HD kWh ½ × Covered Q HD % × Electricity Cost $ kWh ! : 6.3. The Cash Flow Diagram and STES System's Payback Periods for the Three Cities. Based on the calculated system's capital cost and the savings from the electricity bills for the three cities, the cash flow diagrams of the three cities were generated as shown in Figure 11. It can be observed from Figure 11 that the net annual savings from electricity costs are around 22,533 $ kWh -1 , 12,963 $ kWh -1 , and 858 $ kWh -1 for Potsdam, Zarqa, and Doha, respectively. When the capital cost of each city, shown in Table 7, is divided by the net annual savings from electricity cost, the yield payback periods are 3.4 years, 4.4 years, and 29.4 years for Potsdam, Zarqa, and Doha, respectively.

Conclusion
All buildings are built on foundation pillars that are based deep into the ground. To increase the soil consistency between these foundation's pillars, the space is backfilled with durable basaltic materials. The underground space can be utilized for an innovative solution for a seasonal energy storage system that could be buried under the building by using the same excavation for the building's foundation, replacing regular basaltic backfill materials with solid materials that have good thermal properties. This study simulated a hotel building with 400 m 2 floor area that could benefit from STES to store thermal energy during the summer and autumn seasons to take advantage of this stored energy to cover the hotel's heat demand from domestic hot water and space heating loads in the winter. The selected solid materials were silt materials that are available in the global market like silica sand and basalt stone combined with a copper mesh for heat distribution through the STES. The STES system was charged by a solar thermal combisystem for a five-year span with different collector areas in three different climate zones in Potsdam (Germany), Zarqa (Jordan), and Doha (Qatar).