Study on the Effect of Ground Heat Storage by Solar Heat Using Numerical Simulation

Recently, energy storage techniques using renewable energy efficiently have attracted considerable attention. However, there are several problems when using renewable energy. In the case of solar energy, the energy production time is different from the consumption time, and the use of geothermal energy has high investment costs. In order to solve these problems, it is essential to develop high-efficiency systems using both solar and geothermal energy simultaneously and efficiently. Thus, in this study, the performance of underground heat storage of solar energy was examined by simulation using models of underground heat transfer and heat exchange for the development of an integrated hybrid system exploiting both geothermal and solar energy. As a result, the heat extraction performance was determined to be up to 72.75 W/m. As a result, in Kagoshima, the most southern area in Korea, a case of six hour heat storage operation achieved the highest heat exchange rate of 72.75 W/m, which is approximately 105% higher than the case of operation without heat storage.


Introduction
Recently, technologies for saving energy and greenhouse gas reduction have attracted attention from people to prevent global warming.Issues of environmentally friendly architecture and energy saving techniques have also been raised.Currently there are various regulations on energy consumption in different countries all over the world.There are still difficulties in dealing with problems associated with increasing energy demand, especially those during peak electricity demand periods.To solve these problems, technologies enabling the application of new renewable energy have emerged.Of them, geothermal heat pumping systems, has been expected to realize an effective saving in heating and cooling energy because it can exploit the geothermal power that is stable all year round.However, the comparatively high cost associated with the installation of such facilities and limited sites suitable for such installation have slowed the propagation of such facilities.To avoid these obstacles, technology that increases system performance by integrating geothermal source with another heat source system or reduces the capacity of equipment to be installed has emerged as an alternative way to exploit geothermal energy.Besides, heat (thermal energy) storage systems can resolve the issues of temporal discrepancy between the demand and the supply of thermal energy.Ground heat pump systems have the merit of utilizing the relatively stable underground temperatures, however, during long-term heating operation, the underground temperature will gradually go down due to heat extraction from the ground and, as a result, the system performance also decreases.As the solution against this problem, hybrid systems with other heat source systems are suggested and used in the system building field.The use of solar heat could be good choice for heating dominant buildings.Generally, solar heat collection systems consist of a solar heat collector, thermal storage tank, and circulation pump.The solar heat collector is the essential equipment for the system, but a thermal storage tank is additionally required when there is the discrepancy between the energy collection time and use time.Although thermal storage tanks have the advantage of easy application, they require additional space and installation costs, which sometimes are not small.Therefore, in this research, an underground heat storage system using solar heat was developed without the thermal storage tank.
Recently, many researchers have studied systems for underground heat storage [1,2].Wang and Qi [3] have conducted performance analysis of a heat pump system using underground heat storage of solar energy in residential buildings based on experimental data.The results found that the efficiency of the heat storage was over 70%.They studied on suitable ratio between the tank volume and the area of solar collectors, which was found as 20-40 L/m range.Seo et al. [4] conducted an analytical assessment of the performance factor of heat storage systems using solar heat collectors.It was found that the effect of the collector efficiency and heat loss factor on system performance was significant.Park and Kim [5] have identified the applicability of horizontal geothermal systems using heat storage in a greenhouse.Benli and Durmus [6] have used the heat storage system of a greenhouse and performed an analysis during the coldest month of a year.Their results revealed the diameter of pipe at over 1 inch (25.4 mm) with coefficient of performance (COP) of a ground source heat pump (GSHP) system ranging from 2.3 to 3.8.Jeong [7] have examined the static heat storage time of a slab heat storage air conditioning system by analyzing surveyed data and found that the optimal duration of heat storage should be about 7 h.
In addition, research on the media of heat storage systems has been conducted.Vanhoudt et al. [8] have applied an aquifer heat storage system (ATES) in a hospital with a high energy demand level.They conducted an economic analysis based on experimental data obtained from 3 years of monitoring.They found a reduction of 71% in energy consumption compared to a gas boiler installation with a payback period of 8.4 years without any subsidy.The amount of annual cost savings was € 54,000.Hesaraki et al. [9] have compared and reviewed the efficiency of each heat storage systems by examining the seasonal heat storage system literature.They found that seasonal thermal energy storage systems were more efficient for collective regional development than for single dwelling residences with a corresponding level of investment within 20%-30% per square meter of single dwelling residence.Bony and Citherlet [10] have developed a numerical model for analyzing heat storage sets in which phase change material (PCM) is used.After comparing their analysis results to data obtained from monitoring, they verified the feasibility of heat storage using PCM through an analytical and experimental approach.Miro et al. [11] showed that by-products from the potash industry are suitable for thermal energy storage materials through verification experiments.Tiari et al. [12] developed a finned heat pipe-assisted latent heat thermal storage system and analyzed the performance of the system through numerical modeling.As a result, it was confirmed that the natural convection could accelerate the charging and reduce the total charging time of approximately 30%.Furthermore, Shabgard et al. [13] developed a thermal network model in order to examine the performance of a latent heat thermal energy storage system combined with heat pipes.The model used the phase change process, the heat transfer process and the melting-solidification process.
However, few studies have examined the application of geothermal energy for underground heat storage systems.In this study, in order to analyze the performance of an underground storage system with solar heat, numerical simulation was conducted using a ground heat transfer model, a heat exchanger model and a surface heat flux model.This paper is focused on the effect of heat storage rate on the ground heat exchange rate.

Overview of the Simulation
In this study, a 3-dimensional (3D) simulation model was developed by applying a geothermal energy transference model considering phase changes (solid, liquid, gas) in the ground and heat exchanger developed in prior studies [14,15] for accurate analysis of underground heat storage using a solar-geothermal hybrid system.Computation of ground surface heat-flux considering heat exchange around the ground surface and outer air was also carried out.The ground heat transfer model used theoretical calculations for the following: ground surface long wave radiation, short & long wave radiation of the sun, convective heat transfer, and transference of the sensible & latent heat based on the standard weather data [16,17] of the past 10 year averages provided by the meteorological administration.These values were used to set boundary conditions for the model.To predict underground heat exchange rate, analytical code of FEFLOW (finite element subsurface flow) was used, which was designed to perform 3D numerical analysis based on the law of conservation of matter like mass, constitutive material, linear momentum, and energy etc. Analysis equations used in the simulation code, are as follows: the temperature of the circulating water in the U-tube was determined using the formula in Equation ( 1) which is a one-dimension advection-diffusion equation: N " 0.23Re 0.8 Pr n Here, N is 0.4 in heating mode and 0.3 in cooling mode.The heat flux delivered to the underground at the heat exchanger surface was calculated by using the following equation: In addition, the heat flux (Q) of ground surface is calculated by the heat balance in Equation ( 5): Details of the analytical model have been described in prior studies [14,15].An overview of simulation conducted in this study is illustrated in Figure 1.The system consists of a geothermal heat pump system (for heating & cooling) and a solar energy system (for underground heat storage).A simulation was carried out through the analysis model to predict the energy to be transferred from the solar collector to the underground heat storage.
Energies 2015, 8, page-page exchanger developed in prior studies [14,15] for accurate analysis of underground heat storage using a solar-geothermal hybrid system.Computation of ground surface heat-flux considering heat exchange around the ground surface and outer air was also carried out.The ground heat transfer model used theoretical calculations for the following: ground surface long wave radiation, short & long wave radiation of the sun, convective heat transfer, and transference of the sensible & latent heat based on the standard weather data [16,17] of the past 10 year averages provided by the meteorological administration.These values were used to set boundary conditions for the model.To predict underground heat exchange rate, analytical code of FEFLOW (finite element subsurface flow) was used, which was designed to perform 3D numerical analysis based on the law of conservation of matter like mass, constitutive material, linear momentum, and energy etc. Analysis equations used in the simulation code, are as follows: the temperature of the circulating water in the U-tube was determined using the formula in Equation ( 1) which is a one-dimension advection-diffusion equation: Here, N is 0.4 in heating mode and 0.3 in cooling mode.The heat flux delivered to the underground at the heat exchanger surface was calculated by using the following equation: In addition, the heat flux (Q) of ground surface is calculated by the heat balance in Equation ( 5): Details of the analytical model have been described in prior studies [14,15].An overview of simulation conducted in this study is illustrated in Figure 1.The system consists of a geothermal heat pump system (for heating & cooling) and a solar energy system (for underground heat storage).A simulation was carried out through the analysis model to predict the energy to be transferred from the solar collector to the underground heat storage.

Conditions of the Analysis
The model used in this study to predict heat exchange rate and heat storage rate is shown in Figure 2. A tetra-mesh was used to improve the prediction accuracy.The number of meshes was set to 1413 in a slice and the total mesh number was set to 12,717.A model zone of 40 m × 40 m × 100 m was set to take influential coverage of underground heat transfer into account.In the center of the zone, a borehole with 0.2 m in diameter was inserted into a heat exchanger (U-tube, 32A, outer diameter: 32 mm, inner diameter: 26 mm).The gap between the borehole and inserted heat exchanger Energies 2015, 8, 13609-13627

Conditions of the Analysis
The model used in this study to predict heat exchange rate and heat storage rate is shown in Figure 2. A tetra-mesh was used to improve the prediction accuracy.The number of meshes was set to 1413 in a slice and the total mesh number was set to 12,717.A model zone of 40 m ˆ40 m ˆ100 m was set to take influential coverage of underground heat transfer into account.In the center of the zone, a borehole with 0.2 m in diameter was inserted into a heat exchanger (U-tube, 32A, outer diameter: 32 mm, inner diameter: 26 mm).The gap between the borehole and inserted heat exchanger was filled with concrete.The ground was assumed to be granite.Values of thermal properties of the soil, grout, and U-tube are summarized in Table 1.
Energies 2015, 8, page-page 4 was filled with concrete.The ground was assumed to be granite.Values of thermal properties of the soil, grout, and U-tube are summarized in Table 1.The operation schedule is shown in Figure 3.The heat storage operation was set to collect heat and store heat.It was planned as follows: heat collection between 00:00 and 04:00, heat storage between 10:00 and 16:00, and heat collection again between 22:00 and 24:00.This was performed for one month in winter time (from December 1 to December 31).Conditions used for simulation are summarized in Table 2.  Conditions employed for this case study are summarized in Table 3.The case study was designed to examine the influence of changes (weather conditions, square measure of solar energy collector, flow rate, heat storage duration, underground heat conductivity, and grout heat conductivity) on the process of heat collection from underground solar and geothermal energy heat storage.Case 2 with a solar collector and Case 1 without a solar collector were compared to each other to find the influence of underground heat storage of solar energy on heat collection performance.To identify the effect of differences in heat storage rate on heat collection performance, Case 4 with solar energy stored for 3 h and solar collector area of 30 m 2 , Case 5 with solar energy stored for 6 h and solar collector area of 30 m 2 , and Case 6 with solar energy stored for 6 h and solar  The operation schedule is shown in Figure 3.The heat storage operation was set to collect heat and store heat.It was planned as follows: heat collection between 00:00 and 04:00, heat storage between 10:00 and 16:00, and heat collection again between 22:00 and 24:00.This was performed for one month in winter time (from December 1 to December 31).Conditions used for simulation are summarized in Table 2.
Energies 2015, 8, page-page 4 was filled with concrete.The ground was assumed to be granite.Values of thermal properties of the soil, grout, and U-tube are summarized in Table 1.The operation schedule is shown in Figure 3.The heat storage operation was set to collect heat and store heat.It was planned as follows: heat collection between 00:00 and 04:00, heat storage between 10:00 and 16:00, and heat collection again between 22:00 and 24:00.This was performed for one month in winter time (from December 1 to December 31).Conditions used for simulation are summarized in Table 2.

Selection item Heat storage mode Heat extraction mode Operation conditions
Initial ΔT 10 °C limitation temperature 5 °C Soil conditions Granite (λ: 3.5 W/mK) initial temperature 16 °C Simulation period (day) 12/1-12/31(1 month) Operation time (h) 10:00-16:00 00:00-04:00; 22:00-24:00 Conditions employed for this case study are summarized in Table 3.The case study was designed to examine the influence of changes (weather conditions, square measure of solar energy collector, flow rate, heat storage duration, underground heat conductivity, and grout heat conductivity) on the process of heat collection from underground solar and geothermal energy heat storage.Case 2 with a solar collector and Case 1 without a solar collector were compared to each other to find the influence of underground heat storage of solar energy on heat collection performance.To identify the effect of differences in heat storage rate on heat collection performance, Case 4 with solar energy stored for 3 h and solar collector area of 30 m 2 , Case 5 with solar energy stored for 6 h and solar collector area of 30 m 2 , and Case 6 with solar energy stored for 6 h and solar  Conditions employed for this case study are summarized in Table 3.The case study was designed to examine the influence of changes (weather conditions, square measure of solar energy collector, flow rate, heat storage duration, underground heat conductivity, and grout heat conductivity) on the process of heat collection from underground solar and geothermal energy heat storage.Case 2 with a solar collector and Case 1 without a solar collector were compared to each other to find the influence of underground heat storage of solar energy on heat collection performance.To identify the effect of differences in heat storage rate on heat collection performance, Case 4 with solar energy stored for 3 h and solar collector area of 30 m 2 , Case 5 with solar energy stored for 6 h and solar collector area of 30 m 2 , and Case 6 with solar energy stored for 6 h and solar collector area of 50 m 2 were compared to each other.The influence of regional weather conditions on heat collection performance was also examined by comparing regions in Korea and Japan (Seoul, Busan, and Kagoshima).To find the effect of underground heat conductivity on heat collection performance, Cases 9-11 assigned with set value of 1.0 W/mK, 1.5 W/mK, and 2.0 W/mK, respectively, were compared to Case 6 with assigned value of 1.5 W/mK.To determine the influence of grout heat conductivity on heat collection performance, Cases 12-14 with assigned value of 1.0 W/mK, 2.0 W/mK, and 3.0 W/mK, respectively, were compared to Case 6 with assigned value of 1.5 W/mK.Heat collection efficiency of solar collector for underground heat storage (assumed to be about 60%) was calculated using the Equation ( 6): where Q, A s , C, U, r, and E denoted the intensity of direct solar radiation (W/m 2 ), the area of solar collector (m 2 ), the specific heat (J/kg¨˝C), the flow rate (m/s), the radius of U-tube (m), and the efficiency of solar energy collector, respectively.By using the ∆T calculated from Equation ( 1), the input of inlet temperature (T in ) could be calculated using the Equation ( 7):

Results and Discussion
Results of the analysis on heat collection performance of each case with or without the operation of underground heat storage are shown in Figures 4-7.Case 1 was operated with underground heat collection only.For Case 2, a 20 m 2 solar collector was installed.In addition, the operation of Energies 2015, 8, 13609-13627 heat accumulation and heat collection were conducted simultaneously for 3 h.Daily average outlet temperature of circulation water during heat collection is shown in Figure 4.
Figure 4 shows the daily average outlet temperature of circulation water during heating operation in Cases 1 and 2. In Case 1, the ground temperature decreases during the heating operation, because only the heat extraction operation was conducted in this case.Therefore, the circulation water temperature seems to decrease gradually due to the drop of underground temperature.However, in Case 2, the operations of heat extraction and heat storage were conducted simultaneously, so the ground temperature increased with the underground heat storage with solar heat.Therefore, the circulation water temperature also increased with the increased underground temperature.Through this result, it is found that the underground heat storage using solar energy could improve the performance of the GSHP system.The average outlet temperatures of the circulation water calculated at each unit length and heat exchange rate in heat collection are shown in Figure 5.The temperature of circulation water from the heat exchanger in Case 2 appeared to be higher than that in Case 1, proportional to the amount of solar radiation due to the effect of the heat storage operation.Case 2 stored solar energy through underground heat storage.It exhibited a heat collection performance approximately 11% higher than that in Case 1.The changes of heat exchange rate and circulation water temperature from the operation of heat collection are shown in Figures 6 and 7. Due to the underground heat storage of solar energy, the temperature of the circulation water was kept high, which greatly contributed to the improved performance of the heat collection operation.The average outlet temperatures of the circulation water calculated at each unit length and heat exchange rate in heat collection are shown in Figure 5.The temperature of circulation water from the heat exchanger in Case 2 appeared to be higher than that in Case 1, proportional to the amount of solar radiation due to the effect of the heat storage operation.Case 2 stored solar energy through underground heat storage.It exhibited a heat collection performance approximately 11% higher than that in Case 1.The changes of heat exchange rate and circulation water temperature from the operation of heat collection are shown in Figures 6 and 7. Due to the underground heat storage of solar energy, the temperature of the circulation water was kept high, which greatly contributed to the improved performance of the heat collection operation.The average outlet temperatures of the circulation water calculated at each unit length and heat exchange rate in heat collection are shown in Figure 5.The temperature of circulation water from the heat exchanger in Case 2 appeared to be higher than that in Case 1, proportional to the amount of solar radiation due to the effect of the heat storage operation.Case 2 stored solar energy through underground heat storage.It exhibited a heat collection performance approximately 11% higher than that in Case 1.The changes of heat exchange rate and circulation water temperature from the operation of heat collection are shown in Figures 6 and 7. Due to the underground heat storage of solar energy, the temperature of the circulation water was kept high, which greatly contributed to the improved performance of the heat collection operation.Results of the comparative analysis of heat storage rate and variation of circulation water are shown in Figures 8-10.The flow rate during underground heat storage was set at 9.30 L/min.However, for heat collection, the flow rates in Cases 2 and 3 were set at 9.30 L/min and 13.95 L/min, respectively.The changes in heat storage rate and heat exchange rate in Cases 2 and 3 are shown in Figure 8.  Results of the comparative analysis of heat storage rate and variation of circulation water are shown in Figures 8-10.The flow rate during underground heat storage was set at 9.30 L/min.However, for heat collection, the flow rates in Cases 2 and 3 were set at 9.30 L/min and 13.95 L/min, respectively.The changes in heat storage rate and heat exchange rate in Cases 2 and 3 are shown in Figure 8. Results of the comparative analysis of heat storage rate and variation of circulation water are shown in Figures 8-10.The flow rate during underground heat storage was set at 9.30 L/min.However, for heat collection, the flow rates in Cases 2 and 3 were set at 9.30 L/min and 13.95 L/min, respectively.The changes in heat storage rate and heat exchange rate in Cases 2 and 3 are shown in Figure 8.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during underground heat storage are shown in Figure 9.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 10.The results obtained from the analysis showed that the difference in heat storage rates between Cases 2 and 3 was insignificant.The heat exchange rate in Case 3 appeared to be The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during underground heat storage are shown in Figure 9.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 10.The results obtained from the analysis showed that the difference in heat storage rates between Cases 2 and 3 was insignificant.The heat exchange rate in Case 3 appeared to be approximately 19% bigger than that in Case 2. The flow rate appeared to be proportional to the heat exchange rate level.However, the difference in the temperature of circulation water in heat exchangers between the two cases was only 1.33 ˝C (Case 2: 11.40 ˝C, Case 3: 10.07 ˝C).Thus, the value of flow rate in heat collection for follow-up cases was set at 0.483 m/s for the simulation.
Energies 2015, 8, page-page approximately 19% bigger than that in Case 2. The flow rate appeared to be proportional to the heat exchange rate level.However, the difference in the temperature of circulation water in heat exchangers between the two cases was only 1.33 °C (Case 2: 11.40 °C, Case 3: 10.07 °C).Thus, the value of flow rate in heat collection for follow-up cases was set at 0.483 m/s for the simulation.The comparisons between heat storage and heat collection performance according to the differences in the area of solar collector are illustrated in Figures 11 and 12.The comparisons between heat storage and heat collection performance according to the differences in the area of solar collector are illustrated in Figures 11 and 12.
Energies 2015, 8, page-page approximately 19% bigger than that in Case 2. The flow rate appeared to be proportional to the heat exchange rate level.However, the difference in the temperature of circulation water in heat exchangers between the two cases was only 1.33 °C (Case 2: 11.40 °C, Case 3: 10.07 °C).Thus, the value of flow rate in heat collection for follow-up cases was set at 0.483 m/s for the simulation.The comparisons between heat storage and heat collection performance according to the differences in the area of solar collector are illustrated in Figures 11 and 12.To identify the influence of difference in heat storage rate on heat exchange rate, the heat collection performance was examined by setting each solar collector area in Cases 5 and 6 at 30 m 2 and 50 m 2 , respectively.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during underground heat storage are shown in Figure 11.The average To identify the influence of difference in heat storage rate on heat exchange rate, the heat collection performance was examined by setting each solar collector area in Cases 5 and 6 at 30 m 2 and 50 m 2 , respectively.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during underground heat storage are shown in Figure 11.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 12. Analysis of monthly average heat storage rate and heat exchange rate showed that the heat storage rate and heat exchange rate in Case 6 equipped with a bigger solar collector area were 67% and 9% higher than those in Case 5, respectively.Based on these results, the influence of stored heat on the heat collection performance was found to be smaller than that on heat storage.This could be due to heat loss during downtime between the completion of heat storage and the heat collection operation.
Results of the analysis on heat storage and heat collection performance in accordance with varied time of heat storage operation are shown in Figures 13-15.The operation of underground heat storage was continued for 3 h or 6 h (Case 4: 10:00-13:00, Case 5: 10:00-16:00).The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during heat storage are shown in Figure 13.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 14.
Energies 2015, 8, page-page outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 12. Analysis of monthly average heat storage rate and heat exchange rate showed that the heat storage rate and heat exchange rate in Case 6 equipped with a bigger solar collector area were 67% and 9% higher than those in Case 5, respectively.Based on these results, the influence of stored heat on the heat collection performance was found to be smaller than that on heat storage.This could be due to heat loss during downtime between the completion of heat storage and the heat collection operation.
Results of the analysis on heat storage and heat collection performance in accordance with varied time of heat storage operation are shown in Figures 13-15.The operation of underground heat storage was continued for 3 h or 6 h (Case 4: 10:00-13:00, Case 5: 10:00-16:00).The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during heat storage are shown in Figure 13.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 14.The daily average heat storage rates are shown in Figure 15.The average heat storage rate calculated at each unit length in Case 4 that stored heat for 3 h appeared to be 18% higher than that in Case 5 which stored heat for 6 h.However, the total heat storage rate in Case 5 was about 1.6 times more than that in Case 4. The heat exchange rate in Case 5 also appeared to be slightly higher than that in Case 4. Such results could be attributed to the fact that less heat was lost and that more heat amount was accumulated due to higher heat storage rate in Case 5 despite the higher monthly average heat storage rate in Case 4. Thus, the time for heat storage should be determined by taking the cost of power consumption of the heat storage of solar energy into account since the difference in the heat exchange performance in accordance with varied duration of heat storage operation was found to be insignificant.
The average outer air temperature and average amount of solar radiation in each region are shown in Figures 16-19.The weather data [16,17] of the three regions in Korea and Japan (Seoul, Busan, and Kagoshima) represented in Cases 6-8 were referred to those of respective meteorological administrations in both countries.The respective regional conditions are summarized in Table 4. T a , R s , T soil , and north latitude denote the annual mean temperature of outer air, monthly average amount of solar radiation, set underground temperature, and northern latitude, respectively.The monthly average temperature of outer air of each region is shown in Figure 16.The amount of solar radiation in each region is shown in Figure 17.
Energies 2015, 8, 13609-13627 10 Busan, and Kagoshima) represented in Cases 6-8 were referred to those of respective meteorological administrations in both countries.The respective regional conditions are summarized in Table 4. Ta, Rs, Tsoil, and north latitude denote the annual mean temperature of outer air, monthly average amount of solar radiation, set underground temperature, and northern latitude, respectively.The monthly average temperature of outer air of each region is shown in Figure 16.The amount of solar radiation in each region is shown in Figure 17.    10 Busan, and Kagoshima) represented in Cases 6-8 were referred to those of respective meteorological administrations in both countries.The respective regional conditions are summarized in Table 4. Ta, Rs, Tsoil, and north latitude denote the annual mean temperature of outer air, monthly average amount of solar radiation, set underground temperature, and northern latitude, respectively.The monthly average temperature of outer air of each region is shown in Figure 16.The amount of solar radiation in each region is shown in Figure 17.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during heat storage are shown in Figure 18.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 19.During heat storage, the average inlet temperatures of the heat exchanger were 38.32 ˝C in Seoul, 49.27 ˝C in Busan, and 65.89 ˝C in Kagoshima.The heat storage rate in Kagoshima had the highest value at 191.47 W/m, which was about 119% higher than the value in Seoul at 87.37 W/m.The heat exchange rates during heating operation were 57.28 W/m in Seoul, 62.93 W/m in Busan, and 72.75 W/m in Kagoshima.The results indicate that the potential of heat exchange rate is higher in the southern region due to higher annual average temperature of outer air and higher amount of solar radiation in the south compared to the north.The difference in average heat exchange rate between Seoul and Kagoshima was about 27%, which was smaller than the value of 119% for the heat storage.The heat exchange rates during heating operation were 57.28 W/m in Seoul, 62.93 W/m in Busan, and 72.75 W/m in Kagoshima.The results indicate that the potential of heat exchange rate is higher in the southern region due to higher annual average temperature of outer air and higher amount of solar radiation in the south compared to the north.The difference in average heat exchange rate between Seoul and Kagoshima was about 27%, which was smaller than the value of 119% for the heat storage.
Results of the comparative analysis of heat collection performance according to different underground heat conductivity values are shown in Figures 20-22.In Cases 9-11, the values of underground heat conductivity were set at 1.0 W/mK, 2.0 W/mK, and 3.0 W/mK, respectively, which were comparable to the value of Case 6 at 3.5 W/mK.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during heat storage operation are shown in Figure 20.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate in the heat collection operation are illustrated in Figure 21.
The difference in heat storage rate among heat storage cases (Cases 9-11) appeared to be insignificant.For heat collection cases, the values of heat exchange rate in Cases 9-11, and 6 were 61.32 W/m, 59.61 W/m, 58.77W/m, and 57.28 W/m, respectively.

11
The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during heat storage are shown in Figure 18.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 19.During heat storage, the average inlet temperatures of the heat exchanger were 38.32 °C in Seoul, 49.27 °C in Busan, and 65.89 °C in Kagoshima.The heat storage rate in Kagoshima had the highest value at 191.47 W/m, which was about 119% higher than the value in Seoul at 87.37 W/m.The heat exchange rates during heating operation were 57.28 W/m in Seoul, 62.93 W/m in Busan, and 72.75 W/m in Kagoshima.The results indicate that the potential of heat exchange rate is higher in the southern region due to higher annual average temperature of outer air and higher amount of solar radiation in the south compared to the north.The difference in average heat exchange rate between Seoul and Kagoshima was about 27%, which was smaller than the value of 119% for the heat storage.
Results of the comparative analysis of heat collection performance according to different underground heat conductivity values are shown in Figures 20-22.In Cases 9-11, the values of underground heat conductivity were set at 1.0 W/mK, 2.0 W/mK, and 3.0 W/mK, respectively, which were comparable to the value of Case 6 at 3.5 W/mK.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during heat storage operation are shown in Figure 20.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate in the heat collection operation are illustrated in Figure 21.
The difference in heat storage rate among heat storage cases (Cases 9-11) appeared to be insignificant.For heat collection cases, the values of heat exchange rate in Cases 9-11, and 6 were 61.32 W/m, 59.61 W/m, 58.77W/m, and 57.28 W/m, respectively.Results of the comparative analysis of heat collection performance according to different underground heat conductivity values are shown in Figures 20-22.In Cases 9-11, the values of underground heat conductivity were set at 1.0 W/mK, 2.0 W/mK, and 3.0 W/mK, respectively, which were comparable to the value of Case 6 at 3.5 W/mK.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate during heat storage operation are shown in Figure 20.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate in the heat collection operation are illustrated in Figure 21.The maximum difference was about 7%, attributable to the decreased value in accordance with increased underground heat conductivity.Results obtained from the analysis for conditions of 1.0 W/mK and 3.5 W/mK of underground heat conductivity are shown in Figure 22.The temperature change around the heat exchanger situated 30 m beneath the ground was calculated right before and     The maximum difference was about 7%, attributable to the decreased value in accordance with increased underground heat conductivity.Results obtained from the analysis for conditions of 1.0 W/mK and 3.5 W/mK of underground heat conductivity are shown in Figure 22.The temperature change around the heat exchanger situated 30 m beneath the ground was calculated right before and The difference in heat storage rate among heat storage cases (Cases 9-11) appeared to be insignificant.For heat collection cases, the values of heat exchange rate in Cases 9-11, and 6 were 61.32 W/m, 59.61 W/m, 58.77W/m, and 57.28 W/m, respectively.
The maximum difference was about 7%, attributable to the decreased value in accordance with increased underground heat conductivity.Results obtained from the analysis for conditions of 1.0 W/mK and 3.5 W/mK of underground heat conductivity are shown in Figure 22.The temperature change around the heat exchanger situated 30 m beneath the ground was calculated right before and after the heat storage operation.For Case 9 (λ soil : 1.0 W/mK), the underground temperature on the boundary surface of the borehole calculated right after the heat storage operation was 19.4 ˝C, similar to the value of 19.0 ˝C in Case 6 (λ soil : 3.5 W/mK).However, the underground temperature calculated right before the heat collection operation was 11.3 ˝C, revealing a rather bigger difference of 3.4 ˝C compared to 7.9 ˝C in Case 6.From this result, it was identified that the heat accumulated through heat storage could be effectively collected through heat collection.However, the efficiency of heat collection varied greatly due to the different soil heat conductivity values.
Energies 2015, 8,[13609][13610][13611][13612][13613][13614][13615][13616][13617][13618][13619][13620][13621][13622][13623][13624][13625][13626][13627] Results of the analysis for the performance of heat storage and heat collection with different grout heat conductivities are shown in Figures 23-27.For Cases 12-14, the grout heat conductivity was 1.0 W/mK, 2.0 W/mK, and 3.0 W/mK, respectively, which was comparable to the grout heat conductivity at 1.5 W/mK in Case 6.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate in the heat storage operation are shown in Figure 23.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate in the heat collection operation are shown in Figure 24.
The differences in heat storage rate in the respective operations of heat storage among cases appeared to be insignificant.The values of inlet temperatures of each heat exchanger in the cases were calculated as 41.01 ˝C, 38.32 ˝C, 36.72 ˝C, and 35.40 ˝C, respectively.The inlet temperature of each heat exchanger in the cases appeared to be inversely proportional to the level of grout heat conductivity.However, the heat exchange rate in each case of heat collection was 51.17 W/m, 57.28 W/m, 61.16 W/m, and 65.67 W/m, respectively.And the heat collection performance appeared to be proportional to the level of grout heat conductivity.
13 after the heat storage operation.For Case 9 (λsoil: 1.0 W/mK), the underground temperature on the boundary surface of the borehole calculated right after the heat storage operation was 19.4 °C, similar to the value of 19.0 °C in Case 6 (λsoil: 3.5 W/mK).However, the underground temperature calculated right before the heat collection operation was 11.3 °C, revealing a rather bigger difference of 3.4 °C compared to 7.9 °C in Case 6.From this result, it was identified that the heat accumulated through heat storage could be effectively collected through heat collection.However, the efficiency of heat collection varied greatly due to the different soil heat conductivity values.
Results of the analysis for the performance of heat storage and heat collection with different grout heat conductivities are shown in Figures 23-27.For Cases 12-14, the grout heat conductivity was 1.0 W/mK, 2.0 W/mK, and 3.0 W/mK, respectively, which was comparable to the grout heat conductivity at 1.5 W/mK in Case 6.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate in the heat storage operation are shown in Figure 23.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate in the heat collection operation are shown in Figure 24.
The differences in heat storage rate in the respective operations of heat storage among cases appeared to be insignificant.The values of inlet temperatures of each heat exchanger in the cases were calculated as 41.01 °C, 38.32 °C, 36.72 °C, and 35.40 °C, respectively.The inlet temperature of each heat exchanger in the cases appeared to be inversely proportional to the level of grout heat conductivity.However, the heat exchange rate in each case of heat collection was 51.17 W/m, 57.28 W/m, 61.16 W/m, and 65.67 W/m, respectively.And the heat collection performance appeared to be proportional to the level of grout heat conductivity.13 after the heat storage operation.For Case 9 (λsoil: 1.0 W/mK), the underground temperature on the boundary surface of the borehole calculated right after the heat storage operation was 19.4 °C, similar to the value of 19.0 °C in Case 6 (λsoil: 3.5 W/mK).However, the underground temperature calculated right before the heat collection operation was 11.3 °C, revealing a rather bigger difference of 3.4 °C compared to 7.9 °C in Case 6.From this result, it was identified that the heat accumulated through heat storage could be effectively collected through heat collection.However, the efficiency of heat collection varied greatly due to the different soil heat conductivity values.
Results of the analysis for the performance of heat storage and heat collection with different grout heat conductivities are shown in Figures 23-27.For Cases 12-14, the grout heat conductivity was 1.0 W/mK, 2.0 W/mK, and 3.0 W/mK, respectively, which was comparable to the grout heat conductivity at 1.5 W/mK in Case 6.The average inlet temperatures of circulation water calculated at each unit length and heat storage rate in the heat storage operation are shown in Figure 23.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate in the heat collection operation are shown in Figure 24.
The differences in heat storage rate in the respective operations of heat storage among cases appeared to be insignificant.The values of inlet temperatures of each heat exchanger in the cases were calculated as 41.01 °C, 38.32 °C, 36.72 °C, and 35.40 °C, respectively.The inlet temperature of each heat exchanger in the cases appeared to be inversely proportional to the level of grout heat conductivity.However, the heat exchange rate in each case of heat collection was 51.17 W/m, 57.28 W/m, 61.16 W/m, and 65.67 W/m, respectively.And the heat collection performance appeared to be proportional to the level of grout heat conductivity.The temperature changes around the heat exchanger calculated right before and after respective operations of heat storage and heat collection for Case 12 (λgrout = 1.0 W/mK) and Case 14 (λgrout = 3.0 W/mK) were determined to identify the influence of grout heat conductivity on the performance of heat storage and heat collection.As for Case 12 (λgrout = 1.0 W/mK), the temperature change right after heat storage influenced the temperature right before heat collection by borehole during heat storage because the low grout heat conductivity caused comparatively higher temperature around the   The temperature changes around the heat exchanger calculated right before and after respective operations of heat storage and heat collection for Case 12 (λgrout = 1.0 W/mK) and Case 14 (λgrout = 3.0 W/mK) were determined to identify the influence of grout heat conductivity on the performance of heat storage and heat collection.As for Case 12 (λgrout = 1.0 W/mK), the temperature change right after heat storage influenced the temperature right before heat collection by borehole during heat storage because the low grout heat conductivity caused comparatively higher temperature around the The temperature changes around the heat exchanger calculated right before and after respective operations of heat storage and heat collection for Case 12 (λ grout = 1.0 W/mK) and Case 14 (λ grout = 3.0 W/mK) were determined to identify the influence of grout heat conductivity on the performance of heat storage and heat collection.As for Case 12 (λ grout = 1.0 W/mK), the temperature change right after heat storage influenced the temperature right before heat collection by borehole during heat storage because the low grout heat conductivity caused comparatively higher temperature around the borehole than that for Case 14 (λ grout = 3.0 W/mK).This was in consistent Energies 2015, 8, 13609-13627 with our assumption at the beginning of the simulation that an excellent heat collection performance would be facilitated by heat stored in the borehole with low grout heat conductivity.However, the result obtained from heat collection performance analysis during heat collection revealed the higher level of heat collection performance was proportional to the level of grout heat conductivity, with a maximum difference at about 28%.Contrary to the initial expectation, grout heat conductivity influenced heat collection with higher heat collection performance during the early stages of heat collection.However, it also manifested the decrease in efficiency in collecting peripheral heat, resulting in a rapid decrease in efficiency with time.Figure 27 indicates the heat exchange rate with time of the base case (Case 6).
Energies 2015, 8, page-page 15 borehole than that for Case 14 (λgrout = 3.0 W/mK).This was in consistent with our assumption at the beginning of the simulation that an excellent heat collection performance would be facilitated by heat stored in the borehole with low grout heat conductivity.However, the result obtained from heat collection performance analysis during heat collection revealed the higher level of heat collection performance was proportional to the level of grout heat conductivity, with a maximum difference at about 28%.Contrary to the initial expectation, grout heat conductivity influenced heat collection with higher heat collection performance during the early stages of heat collection.However, it also manifested the decrease in efficiency in collecting peripheral heat, resulting in a rapid decrease in efficiency with time.Figure 27 indicates the heat exchange rate with time of the base case (Case 6).

Conclusions
In this study, the performance of underground heat storage of solar energy was examined by simulation using models of underground heat transference and heat exchange for the development of integrated hybrid system exploiting geothermal and solar energy.Results obtained from the study are summarized in the Table 5.
Energies 2015, 8, 13609-13627 ‚ The effect of flow rate commonly set at 9.30 L/min for heat storage and at 9.30 L/min or 13.95 L/min for heat collection was analyzed.Our results revealed that, although the difference in heat storage amount of solar energy was insignificant, about a 19% increase in underground heat exchange rate was found.
‚ The influence of varied square measure of solar collector was also examined by employing the conditions of 30 m 2 and 50 m 2 solar collector area.Our results revealed that the difference in the amount of solar energy collection was about 67% and the difference in underground heat exchange rate was about 9%.Thus, the change in underground heat storage rate was not proportional to the change in heat collection performance.
‚ The effect of changes in the duration of heat storage operation was also analyzed by employing different conditions of the duration of underground heat storage operation for 3 h (Case 4) and for 6 h (Case 5).Result of the analysis showed the values of daily average heat storage rate were 32,645.76kWh for Case 4 and 52,235.83kWh for Case 5. Therefore, the proportionality factor was about 1.6.However, the value of the underground heat exchange rate was 51.42 W/m for Case 4 and 52.69 W/m for Case 5, indicating an insignificant effect of different durations of operation on the underground heat collection performance.Thus, the operation of heat storage should take the cost of operation into account in determining the duration of operation which could have an insignificant effect on heat collection performance.
‚ The effect of varied weather conditions on heat collection performance was also examined by simulation of heat collection using different weather conditions corresponding to three regions (Seoul and Busan in Korea and Kagoshima in Japan).Results obtained from the analysis revealed the highest level of heat exchange rate occurred in Kagoshima, at about 119% of that in Seoul.The level of annual average heat exchange rate in Kagoshima was about 27% higher than the other regions.Thus, southern regions with higher annual average outer air temperature and higher annual average amount of solar radiation might affect the performance of underground heat storage and heat collection thereof.
‚ The level of underground heat conductivity was identified to be inversely proportional to the level of heat collection performance via underground heat storage.This could be due to increased temperature of the area around the borehole where thermal diffusion was relatively retarded by the lower level of underground heat conductivity.
‚ The level of grout heat conductivity appeared to have no significant effect on heat storage performance.However, a higher level of grout heat conductivity could facilitate the enhancement of underground heat collection performance via underground heat storage.This is because low grout heat conductivity could enable efficient underground heat storage.However, low grout heat conductivity could be disadvantageous in collecting peripheral heat during heat collection.
Hence, further studies considering more factors such as borehole depth and the shape of heat exchangers should be carried out to develop an optimal design for this system.In addition, the cost of installation and operation thereof should be comprehensively evaluated.

Figure 7 .
Figure 7. Operation of heat storage and heat exchange rate (Cases 1 and 2).

Figure 8 .
Figure 8. Operation of heat storage and heat exchange rate (Cases 2 and 3).

Figure 7 .
Figure 7. Operation of heat storage and heat exchange rate (Cases 1 and 2).

Figure 7 .
Figure 7. Operation of heat storage and heat exchange rate (Cases 1 and 2).

Figure 8 .
Figure 8. Operation of heat storage and heat exchange rate (Cases 2 and 3).Figure 8. Operation of heat storage and heat exchange rate (Cases 2 and 3).

Figure 8 .
Figure 8. Operation of heat storage and heat exchange rate (Cases 2 and 3).Figure 8. Operation of heat storage and heat exchange rate (Cases 2 and 3).

Figure 8 .
Figure 8. Operation of heat storage and heat exchange rate (Cases 2 and 3).

Figure 9 .
Figure 9. Average of heat storage rate for flow rate change (Cases 2 and 3).

Figure 9 .
Figure 9. Average of heat storage rate for flow rate change (Cases 2 and 3).

Figure 10 .
Figure 10.Average of heat exchange rate for flow rate change (Cases 2 and 3).

Figure 11 .
Figure 11.Average of heat storage rate for solar collector area change (Cases 5 and 6).

Figure 10 .
Figure 10.Average of heat exchange rate for flow rate change (Cases 2 and 3).

Figure 10 .
Figure 10.Average of heat exchange rate for flow rate change (Cases 2 and 3).

Figure 11 .
Figure 11.Average of heat storage rate for solar collector area change (Cases 5 and 6).

Figure 11 .
Figure 11.Average of heat storage rate for solar collector area change (Cases 5 and 6).

Figure 11 .
Figure 11.Average of heat storage rate for solar collector area change (Cases 5 and 6).

Figure 12 .
Figure 12.Average of heat exchange rate for solar collector area change (Cases 5 and 6).

Figure 12 .
Figure 12.Average of heat exchange rate for solar collector area change (Cases 5 and 6).

Figure 13 .
Figure 13.Average of heat storage rate for storage time change (Cases 4 and 5).Figure 13.Average of heat storage rate for storage time change (Cases 4 and 5).

Figure 13 .
Figure 13.Average of heat storage rate for storage time change (Cases 4 and 5).Figure 13.Average of heat storage rate for storage time change (Cases 4 and 5).

Figure 13 .
Figure 13.Average of heat storage rate for storage time change (Cases 4 and 5).

Figure 14 .
Figure 14.Average of heat exchange rate for storage time change (Cases 4 and 5).

Figure 15 .
Figure 15.Daily average of heat storage rate for storage time change (Cases 4 and 5).

Figure 14 .
Figure 14.Average of heat exchange rate for storage time change (Cases 4 and 5).

Figure 13 .
Figure 13.Average of heat storage rate for storage time change (Cases 4 and 5).

Figure 14 .
Figure 14.Average of heat exchange rate for storage time change (Cases 4 and 5).

Figure 15 .
Figure 15.Daily average of heat storage rate for storage time change (Cases 4 and 5).Figure 15.Daily average of heat storage rate for storage time change (Cases 4 and 5).

Figure 15 .
Figure 15.Daily average of heat storage rate for storage time change (Cases 4 and 5).Figure 15.Daily average of heat storage rate for storage time change (Cases 4 and 5).

Energies 2015, 8 ,
13609-13627 11 storage rate during heat storage are shown in Figure 18.The average outlet temperatures of circulation water calculated at each unit length and heat exchange rate during heat collection are shown in Figure 19.During heat storage, the average inlet temperatures of the heat exchanger were 38.32 °C in Seoul, 49.27 °C in Busan, and 65.89 °C in Kagoshima.The heat storage rate in Kagoshima had the highest value at 191.47 W/m, which was about 119% higher than the value in Seoul at 87.37 W/m.

Figure 18 .
Figure 18.Average of heat storage rate for simulation area change (Seoul and Busan, Kagoshima).

Figure 19 .
Figure 19.Average of heat exchange rate for simulation area change (Seoul and Busan, Kagoshima).

Figure 18 .
Figure 18.Average of heat storage rate for simulation area change (Seoul and Busan, Kagoshima).

Figure 18 .
Figure 18.Average of heat storage rate for simulation area change (Seoul and Busan, Kagoshima).

Figure 19 .
Figure 19.Average of heat exchange rate for simulation area change (Seoul and Busan, Kagoshima).

Figure 19 .
Figure 19.Average of heat exchange rate for simulation area change (Seoul and Busan, Kagoshima).

Figure 20 .
Figure 20.Average of heat storage rate for soil conductivity change (Cases 6 and 9-11).Figure 20.Average of heat storage rate for soil conductivity change (Cases 6 and 9-11).

Figure 20 .
Figure 20.Average of heat storage rate for soil conductivity change (Cases 6 and 9-11).

Figure 21 .
Figure 21.Average of heat exchange rate for soil conductivity change (Cases 6 and 9-11).

Figure 21 .
Figure 21.Average of heat exchange rate for soil conductivity change (Cases 6 and 9-11).

Figure 20 .
Figure 20.Average of heat storage rate for soil conductivity change (Cases 6 and 9-11).

Figure 21 .
Figure 21.Average of heat exchange rate for soil conductivity change (Cases 6 and 9-11).

Table 4 .
Simulation area conditions.

Table 4 .
Simulation area conditions.

Table 4 .
Simulation area conditions.

T a ( ˝C) R s (W/m 2 ) T Soil ( ˝C) North latitude (
surf Upward long wave radiation from the ground surface H surf Sensible heat flux L surf Latent heat flux