Application of Thermal Labyrinth System to Reduce Heating and Cooling Energy Consumption

To reduce the energy consumption in buildings, modern buildings are increasingly becoming airtight. In these structures, the outdoor air is supplied inside through mechanical ventilation systems, which are essential for ensuring comfortable indoor air quality. However, these systems consume a considerable amount of energy in buildings. One potential solution is using a thermal labyrinth system, which is buried underneath the building. It can pre-cool or pre-heat the outdoor air through heat transfer with the surrounding soil. In this research, a number of case studies were conducted to optimize the thermal labyrinth design. The optimized thermal labyrinth system was derived using computational fluid dynamics (CFD) simulation. In addition, operation algorithms were developed for the efficient operation of the thermal labyrinth system in buildings. The results indicated that there were five operation modes, and the thermal labyrinth could be operated for seven months of the year. The energy reduction effects of the thermal labyrinth system were analyzed and were assessed by the transient system simulation (TRNSYS) tool. A 12% reduction in the annual heating and cooling energy was achieved by applying the thermal labyrinth system.


Background
As an increasing number of air-tight buildings are being constructed in order to save energy, optimum ventilation systems that allows the outdoor air into a building have become essential [1].When the outdoor air is used to operate the heating, ventilation, and air conditioning (HVAC) system, the cooling or heating energy consumption in the building would increase in order to adjust the air temperature to a set point.Therefore, to reduce the energy consumption, an energy-saving HVAC system is required.
Energy-saving HVAC systems that use ground heat are mainly divided into tube-style (earth tube, earth tunnel, and earth-air-heat exchanger (EAHE)) and thermal labyrinth systems.Various studies have been conducted on the tube-style systems and verified its energy-saving effect [1,2].Further to their findings, research has focused on improving the EAHE heat transfer efficiency in terms of length and diameter of the tube type heat exchanger, mass flow rate, laying depth, and changes in soil [3][4][5].Previous literature has indicated that the tube-style system could apply to various climate conditions such as Italy [6], Turkey [7] and Malaysia [8].Moreover, Yu et al. [9], Ghosal et al. [10] and Bansal et al. [11] sought to maximize the system energy saving effect by combining the EAHE with another passive system or an active system, such as a solar chimney, air conditioner or an evaporative cooler.Yang et al. [12] described that the use of tube-style ventilation must be considered alongside building characteristics and categories because it cannot be quantified in isolation from the building characteristics.
Song et al. [13] analyzed the cooling, heating, and dehumidification of the thermal labyrinth system by measuring the internal temperature and humidity.Using an air current analysis program, Sohn et al. [14] analyzed the heat transfer efficiency in the ground-heat-based thermal labyrinth system based on the increase in the heat transfer area according to the change in shape of the openings which introduce outdoor air in the system.
Various studies have focused on the simulation and assessment of the tube-type ground-heatbased systems (such as EAHE, Earth Tube, and Earth Tunnel); however, few have focused on the thermal labyrinth system.In addition, few studies have investigated the potential improvements to the heat transfer efficiency of the thermal labyrinth system and its effective mode of operation.
A thermal labyrinth system introduces the outdoor air into a building through a thermal labyrinth, which is a concrete structure installed underground, and uses the ground heat produced by heat transfer to ventilate the building.According to Song et al. [13] and Sohn et al. [14], the thermal labyrinth system is easier to maintain and offers better energy efficiency than the tube type systems, because it possesses a larger contact area between the thermal labyrinth system and the outdoor air.Moreover, the thermal labyrinth system does not require an additional space, because it is installed underneath buildings.Therefore, this system is more suitable for South Korea, given the limited land availability.

Objectives
The objectives of this study were as follows: (1) To produce a thermal labyrinth design that can offer maximum heat transfer using the limited underground space of a building based on various case studies.(2) To propose an operational algorithm based on the performance analysis of the optimized thermal labyrinth system.(3) To reduce the yearly heating and cooling energy consumption of buildings by using the optimized thermal labyrinth design plan and the operational algorithm.

Description of the Sample Building and the Used Thermal Labyrinth System
In this section, various thermal labyrinth system designs that ensured maximum heat transfer efficiency were derived through different case studies by optimizing the design of the 412.68 m 2 concrete underground space.Figure 1 shows a sample building and its thermal labyrinth.The thermal labyrinth system was located three meters under the building and was constructed with concrete walls.Table 1 shows the conditions of the studied building and the thermal labyrinth system used in the case study.

Modeling
Accurately depicting the thermal conduction through soil is important for simulating the thermal labyrinth system using computational fluid dynamics (CFD) simulation.A previous study simulated the thermal conduction by modeling the cooling tube with the surrounding earth when analyzing the ground-heat-based heat transfer system, using a fluid analysis program [14] to make the simulated results as close to the actual system as possible.In this case study, the thermal labyrinth system and underlying ground were also modeled together (Figure 2).

Case Conditions
According to the previous studies, the methods that can be used to increase the ground heat-based system efficiency are: Changes in the speed of the internal air current; changes in the soil conditions; and changes in the length, section area, and heat transfer surface.Since the air current speed used in the thermal labyrinth system in the target building is determined according to the amount of ventilated air required for the target HVAC space, the air current speed is considered a constant and cannot be changed.In addition, once the site of the target building is determined, the soil conditions can no longer be changed.Therefore, in this study, changes in the surface area of the internal thermal labyrinth system were used as the case's only variable condition (Table 2).

Boundary Conditions and Solver Control
The thermal labyrinth system reduces the heating and cooling loads through heat transfer with surrounding soil which maintains the temperature below a certain temperature throughout the year.Therefore, the soil temperature is very important in analyzing the effect of the pre-heating and pre-cooling of the thermal labyrinth system.This temperature is affected by various factors in addition to climate, such as the amount of moisture in the soil, flow pattern of the ground water, or shaded areas [15].Therefore, it is difficult to describe the ground temperature data based only on the regional time and the ground depth.The ground temperatures are more influenced by weather and have greater diurnal variation near the surface.The influence of diurnal variation generally disappears about 1 m underground, and the influence of weather fades at about 6-7 m underground.Therefore, while there are significant monthly temperature changes due to solar radiation and ambient air near the ground surface, the ground temperature tends to stay constant as the depth increases, and it can be assessed with a

Case Conditions
According to the previous studies, the methods that can be used to increase the ground heat-based system efficiency are: Changes in the speed of the internal air current; changes in the soil conditions; and changes in the length, section area, and heat transfer surface.Since the air current speed used in the thermal labyrinth system in the target building is determined according to the amount of ventilated air required for the target HVAC space, the air current speed is considered a constant and cannot be changed.In addition, once the site of the target building is determined, the soil conditions can no longer be changed.Therefore, in this study, changes in the surface area of the internal thermal labyrinth system were used as the case's only variable condition (Table 2).Fin on the wall and floor (1500 mm gap)

Boundary Conditions and Solver Control
The thermal labyrinth system reduces the heating and cooling loads through heat transfer with surrounding soil which maintains the temperature below a certain temperature throughout the year.Therefore, the soil temperature is very important in analyzing the effect of the pre-heating and pre-cooling of the thermal labyrinth system.This temperature is affected by various factors in addition to climate, such as the amount of moisture in the soil, flow pattern of the ground water, or shaded areas [15].Therefore, it is difficult to describe the ground temperature data based only on the regional time and the ground depth.The ground temperatures are more influenced by weather and have greater diurnal variation near the surface.The influence of diurnal variation generally disappears about 1 m underground, and the influence of weather fades at about 6-7 m underground.Therefore, while there are significant monthly temperature changes due to solar radiation and ambient air near the ground surface, the ground temperature tends to stay constant as the depth increases, and it can be assessed with a periodic function.In this study, it was assumed that the target soil had a constant heat dissipation ratio.Using the heat dissipation coefficient, based on the heat dissipation ratio and the average temperature of The coefficients used in Equation (1) were the ground temperature coefficients in Incheon (Table 3) [17].The results of the calculation of the ground temperature using Equation (1) showed that below about 15 m, the ground temperature was 14.2 • C. Hence, the temperature of the lower part of the earth was set to be 14.2 • C in the simulation model shown in Equation (1).The heat transfer into the thermal labyrinth side could be determined by inputting the earth's thermal conductivity into the simulation (Table 4).Figure 3 is the diagram of heat transfer between soil, concrete and air.periodic function.In this study, it was assumed that the target soil had a constant heat dissipation ratio.Using the heat dissipation coefficient, based on the heat dissipation ratio and the average temperature of the surface, the distribution of the ground temperature was calculated as a function of time, based on Equation (1) [16].
, =   −   exp −  365 The coefficients used in Equation (1) were the ground temperature coefficients in Incheon (Table 3) [17].The results of the calculation of the ground temperature using Equation (1) showed that below about 15 m, the ground temperature was 14.2 °C.Hence, the temperature of the lower part of the earth was set to be 14.2 °C in the simulation model shown in Equation (1).The heat transfer into the thermal labyrinth side could be determined by inputting the earth's thermal conductivity into the simulation (Table 4).Figure 3 is the diagram of heat transfer between soil, concrete and air.The ceiling of the thermal labyrinth was adiabatic as it meets with the upper part of the structure.The part of the outer wall that touches the ground required additional arrangement as the ground temperature varied with depth.The present study estimated the temperature of this wall The ceiling of the thermal labyrinth was adiabatic as it meets with the upper part of the structure.The part of the outer wall that touches the ground required additional arrangement as the ground temperature varied with depth.The present study estimated the temperature of this wall section Energies 2018, 11, 2762 6 of 17 through linear interpolation and the basic trapezoid rule, and then uniformly applied the averaged value of the result to the wall section.
In Section 2, a CFD simulation software tool called star CCM+ was used to derive the optimal thermal labyrinth.It allows a detailed analysis of the ambient air passing through the thermal labyrinth and the thermal conductivity allows an accurate analysis of the heat transfer profile.For a more accurate analysis of the indoor thermal environment, the standard k-ε low-Reynolds-number model that allows an accurate analysis of heat transfer on the surface was applied.The low-Reynolds-number model can analyze the effect of viscosity on the surface more precisely by generating a uniform grid on the surface.Although it takes longer than the high-Reynolds-number model that makes rough calculations of the surface grid through the wall function at the surface level, it provides precision, making it suitable for assessing the thermal labyrinth system that mostly deals with heat transfer on the surface.The following are the detailed equations (Equations ( 2)-( 7)) applied to the standard k-ε low-Reynolds-number model.
Continuity Equation: Equation of Motion: Energy Equation: Energy Dissipation Equation: In the above equation, the applied values for each constant are shown in Table 5.The flow of the sublayer, which is affected by the viscosity of the wall surface, was analyzed by the low-y+ model.The value of y+ is set to <1 in the standard k-ε low-Reynolds-number model used in the present study.The total numbers of mesh elements are around 500,000-700,000 for the 7 studied cases.Detailed input conditions are shown in Table 6 and formed meshes are like Figure 4.

Assessment of the Heat Transfer Efficiency According to the Length of the Internal Thermal Labyrinth Area
The results for the comparative analysis of the heat transfer effect according to different lengths of the internal thermal labyrinth area are summarized in Table 7, and divided into cases.
Figure 5 showed the pre-cooling and pre-heating effect of Case 1-2.The simulation results showed that as the length increased, the temperature difference between the inlet and outlet and the consequent amount of the heat acquired scarcely changed.However, the pressure difference showed a rapid increase (Figure 6).This is because the increase in the length of the air current path inside the fixed space reduced the sectional area of the space where the air travelled.
As shown in Figure 7, there was a slight increase in the fan power due to the increase in the outlet pressure in the thermal labyrinth from Case 1-1 to Case 1-3.Accordingly, the coefficient of performance (COP) was shown to have been the highest in Case 1-1, where the pressure was the lowest.However, in all cases, the differences in the COP were similar; about 0.005.Thus, it was determined that there was no significant difference in the heat transfer effect unless the area that is connected to the ground changes drastically; in addition, if this area cannot be increased drastically, then it would be most efficient to minimize the increase in pressure.

Assessment of the Heat Transfer Efficiency According to the Length of the Internal Thermal Labyrinth Area
The results for the comparative analysis of the heat transfer effect according to different lengths of the internal thermal labyrinth area are summarized in Table 7, and divided into cases.
Figure 5 showed the pre-cooling and pre-heating effect of Case 1-2.The simulation results showed that as the length increased, the temperature difference between the inlet and outlet and the consequent amount of the heat acquired scarcely changed.However, the pressure difference showed a rapid increase (Figure 6).This is because the increase in the length of the air current path inside the fixed space reduced the sectional area of the space where the air travelled.
As shown in Figure 7, there was a slight increase in the fan power due to the increase in the outlet pressure in the thermal labyrinth from Case 1-1 to Case 1-3.Accordingly, the coefficient of performance (COP) was shown to have been the highest in Case 1-1, where the pressure was the lowest.However, in all cases, the differences in the COP were similar; about 0.005.Thus, it was determined that there was no significant difference in the heat transfer effect unless the area that is connected to the ground changes drastically; in addition, if this area cannot be increased drastically, then it would be most efficient to minimize the increase in pressure.

Case Length (m) Modeling
1-1 126 1-2 162 1-3 257           The comparative analysis of the heat transfer according to the change in the internal shape was categorized into cases as shown in Table 8.Owing to the internal embossed surfaces of the system, the heat transfer area between the air that passes through the interior and the ground may increase, and vortexes would be formed around these surfaces, which may result in enhanced heat delivery.
The simulation results showed that there were little changes in the temperature and the acquired heat between Case 2-1, which did not have any embossed surfaces, and Case 2-2, where there were embossed surfaces only on the walls (Figure 8).However, in Case 2-3, where there were embossed surfaces on the floor, the temperature change was relatively high, as was the subsequent heat acquired, and the pressure increased, so the overall COP was high (Figure 9).This shows that the heat transfer efficiency between the internal air and the floor improved due to the increase in the surface area of the floor.In Case 2-4, where the embossed surfaces were more tightly formed, the temperature difference between the inlet and the outlet was larger than in the other cases, thereby showing better heat transfer efficiency.
Table 8.The studied cases for analyzing the impact of the internal form change on the heat transfer efficiency.The comparative analysis of the heat transfer according to the change in the internal shape was categorized into cases as shown in Table 8.Owing to the internal embossed surfaces of the system, the heat transfer area between the air that passes through the interior and the ground may increase, and vortexes would be formed around these surfaces, which may result in enhanced heat delivery.The comparative analysis of the heat transfer according to the change in the internal shape was categorized into cases as shown in Table 8.Owing to the internal embossed surfaces of the system, the heat transfer area between the air that passes through the interior and the ground may increase, and vortexes would be formed around these surfaces, which may result in enhanced heat delivery.

Case Shape Variation Modeling
The simulation results showed that there were little changes in the temperature and the acquired heat between Case 2-1, which did not have any embossed surfaces, and Case 2-2, where there were embossed surfaces only on the walls (Figure 8).However, in Case 2-3, where there were embossed surfaces on the floor, the temperature change was relatively high, as was the subsequent heat acquired, and the pressure increased, so the overall COP was high (Figure 9).This shows that the heat transfer efficiency between the internal air and the floor improved due to the increase in the surface area of the floor.In Case 2-4, where the embossed surfaces were more tightly formed, the temperature difference between the inlet and the outlet was larger than in the other cases, thereby showing better heat transfer efficiency.
Table 8.The studied cases for analyzing the impact of the internal form change on the heat transfer efficiency.The comparative analysis of the heat transfer according to the change in the internal shape was categorized into cases as shown in Table 8.Owing to the internal embossed surfaces of the system, the heat transfer area between the air that passes through the interior and the ground may increase, and vortexes would be formed around these surfaces, which may result in enhanced heat delivery.

Case Shape Variation Modeling
The simulation results showed that there were little changes in the temperature and the acquired heat between Case 2-1, which did not have any embossed surfaces, and Case 2-2, where there were embossed surfaces only on the walls (Figure 8).However, in Case 2-3, where there were embossed surfaces on the floor, the temperature change was relatively high, as was the subsequent heat acquired, and the pressure increased, so the overall COP was high (Figure 9).This shows that the heat transfer efficiency between the internal air and the floor improved due to the increase in the surface area of the floor.In Case 2-4, where the embossed surfaces were more tightly formed, the temperature difference between the inlet and the outlet was larger than in the other cases, thereby showing better heat transfer efficiency.
Table 8.The studied cases for analyzing the impact of the internal form change on the heat transfer efficiency.

Assessment of the Heat Transfer Efficiency According to the Internal Shape
The comparative analysis of the heat transfer according to the change in the internal shape was categorized into cases as shown in Table 8.Owing to the internal embossed surfaces of the system, the heat transfer area between the air that passes through the interior and the ground may increase, and vortexes would be formed around these surfaces, which may result in enhanced heat delivery.
The simulation results showed that there were little changes in the temperature and the acquired heat between Case 2-1, which did not have any embossed surfaces, and Case 2-2, where there were embossed surfaces only on the walls (Figure 8).However, in Case 2-3, where there were embossed surfaces on the floor, the temperature change was relatively high, as was the subsequent heat acquired, and the pressure increased, so the overall COP was high (Figure 9).This shows that the heat transfer efficiency between the internal air and the floor improved due to the increase in the surface area of the floor.In Case 2-4, where the embossed surfaces were more tightly formed, the temperature difference between the inlet and the outlet was larger than in the other cases, thereby showing better heat transfer efficiency.
Table 8.The studied cases for analyzing the impact of the internal form change on the heat transfer efficiency.

Assessment of the Heat Transfer Efficiency According to the Internal Shape
The comparative analysis of the heat transfer according to the change in the internal shape was categorized into cases as shown in Table 8.Owing to the internal embossed surfaces of the system, the heat transfer area between the air that passes through the interior and the ground may increase, and vortexes would be formed around these surfaces, which may result in enhanced heat delivery.
The simulation results showed that there were little changes in the temperature and the acquired heat between Case 2-1, which did not have any embossed surfaces, and Case 2-2, where there were embossed surfaces only on the walls (Figure 8).However, in Case 2-3, where there were embossed surfaces on the floor, the temperature change was relatively high, as was the subsequent heat acquired, and the pressure increased, so the overall COP was high (Figure 9).This shows that the heat transfer efficiency between the internal air and the floor improved due to the increase in the surface area of the floor.In Case 2-4, where the embossed surfaces were more tightly formed, the temperature difference between the inlet and the outlet was larger than in the other cases, thereby showing better heat transfer efficiency.
Table 8.The studied cases for analyzing the impact of the internal form change on the heat transfer efficiency.The simulation results showed that there were little changes in the temperature and the acquired heat between Case 2-1, which did not have any embossed surfaces, and Case 2-2, where there were embossed surfaces only on the walls (Figure 8).However, in Case 2-3, where there were embossed surfaces on the floor, the temperature change was relatively high, as was the subsequent heat acquired, and the pressure increased, so the overall COP was high (Figure 9).This shows that the heat transfer efficiency between the internal air and the floor improved due to the increase in the surface area of the floor.In Case 2-4, where the embossed surfaces were more tightly formed, the temperature difference between the inlet and the outlet was larger than in the other cases, thereby showing better heat transfer efficiency.

Analysis of the Performance of the Optimized Thermal Labyrinth
Based on the previous case studies, the design methods in the cases where the COP was high were integrated to derive an optimized thermal labyrinth that can produce maximum efficiency within the given underground space.Moreover, to determine the yearly heat transfer efficiency of the optimized thermal labyrinth, a CFD simulation was performed.The monthly average temperature was set to be the thermal labyrinth inlet temperature, and the other input conditions followed those used in the preceding case studies (Table 9).
The simulation results showed that in January, when the thermal labyrinth inlet temperature was the lowest, the inlet-outlet temperature difference was the largest at 6.17 °C, and the COP was 2.49.In August, when the thermal labyrinth inlet temperature was the highest, the inlet-outlet temperature difference was 3.78 °C, and the COP was 1.53.On the other hand, in April and October (transitional periods) the inlet-outlet temperature difference was small, so the COP dropped below 1.0 (Figure 10).
Based on the annual optimized thermal labyrinth performance data and using multiple regression analysis, the relationship between the outside air and the thermal labyrinth passage temperature is expressed as follows:

Analysis of the Performance of the Optimized Thermal Labyrinth
Based on the previous case studies, the design methods in the cases where the COP was high were integrated to derive an optimized thermal labyrinth that can produce maximum efficiency within the given underground space.Moreover, to determine the yearly heat transfer efficiency of the optimized thermal labyrinth, a CFD simulation was performed.The monthly average temperature was set to be the thermal labyrinth inlet temperature, and the other input conditions followed those used in the preceding case studies (Table 9).
The simulation results showed that in January, when the thermal labyrinth inlet temperature was the lowest, the inlet-outlet temperature difference was the largest at 6.17 °C, and the COP was 2.49.In August, when the thermal labyrinth inlet temperature was the highest, the inlet-outlet temperature difference was 3.78 °C, and the COP was 1.53.On the other hand, in April and October (transitional periods) the inlet-outlet temperature difference was small, so the COP dropped below 1.0 (Figure 10).
Based on the annual optimized thermal labyrinth performance data and using multiple regression analysis, the relationship between the outside air and the thermal labyrinth passage temperature is expressed as follows:

Analysis of the Performance of the Optimized Thermal Labyrinth
Based on the previous case studies, the design methods in the cases where the COP was high were integrated to derive an optimized thermal labyrinth that can produce maximum efficiency within the given underground space.Moreover, to determine the yearly heat transfer efficiency of the optimized thermal labyrinth, a CFD simulation was performed.The monthly average temperature was set to be the thermal labyrinth inlet temperature, and the other input conditions followed those used in the preceding case studies (Table 9).
The simulation results showed that in January, when the thermal labyrinth inlet temperature was the lowest, the inlet-outlet temperature difference was the largest at 6.17 • C, and the COP was 2.49.In August, when the thermal labyrinth inlet temperature was the highest, the inlet-outlet temperature difference was 3.78 • C, and the COP was 1.53.On the other hand, in April and October (transitional periods) the inlet-outlet temperature difference was small, so the COP dropped below 1.0 (Figure 10).
Based on the annual optimized thermal labyrinth performance data and using multiple regression analysis, the relationship between the outside air and the thermal labyrinth passage temperature is expressed as follows: y = 0.0027x 2 + 0.5989x + 5.0246 (8)

Operation Mode
In this section, the results of the performance analysis of the previously discussed optimal thermal labyrinth system were used to produce the optional mode of operation.The algorithm in Figure 11 was determined based on the internal temperature, outdoor air temperature, air temperature after the outdoor air passed through the thermal labyrinth, and ground temperature.As shown in Figure 11, if the outdoor air temperature is lower than the ground temperature at 14.2 °C, then the temperature of the air after it passes through the thermal labyrinth would always be higher than the outdoor air temperature.However, if the outdoor air is higher than 14.2 °C, then the temperature of the air after it passes through the thermal labyrinth would be lower than the outdoor air temperature.Therefore, if the outdoor air temperature is lower than 14.2 °C during the heating season, the air that passes through the thermal labyrinth would be better used to heat the air handling unit (AHU) and to be supplied indoors.However, if the outdoor air temperature is

Operation Mode
In this section, the results of the performance analysis of the previously discussed optimal thermal labyrinth system were used to produce the optional mode of operation.The algorithm in Figure 11 was determined based on the internal temperature, outdoor air temperature, air temperature after the outdoor air passed through the thermal labyrinth, and ground temperature.

Operation Mode
In this section, the results of the performance analysis of the previously discussed optimal thermal labyrinth system were used to produce the optional mode of operation.The algorithm in Figure 11 was determined based on the internal temperature, outdoor air temperature, air temperature after the outdoor air passed through the thermal labyrinth, and ground temperature.As shown in Figure 11, if the outdoor air temperature is lower than the ground temperature at 14.2 °C, then the temperature of the air after it passes through the thermal labyrinth would always be higher than the outdoor air temperature.However, if the outdoor air is higher than 14.2 °C, then the temperature of the air after it passes through the thermal labyrinth would be lower than the outdoor air temperature.Therefore, if the outdoor air temperature is lower than 14.2 °C during the heating season, the air that passes through the thermal labyrinth would be better used to heat the air handling unit (AHU) and to be supplied indoors.However, if the outdoor air temperature is As shown in Figure 11, if the outdoor air temperature is lower than the ground temperature at 14.2 • C, then the temperature of the air after it passes through the thermal labyrinth would always be higher than the outdoor air temperature.However, if the outdoor air is higher than 14.2 • C, then the temperature of the air after it passes through the thermal labyrinth would be lower than the outdoor air temperature.Therefore, if the outdoor air temperature is lower than 14.2 • C during the heating season, the air that passes through the thermal labyrinth would be better used to heat the air handling unit (AHU) and to be supplied indoors.However, if the outdoor air temperature is higher than 14.2 • C, the outdoor air would be better used for air-conditioning.Moreover, when the indoor temperature is within the comfortable range of 20-26 • C, the outdoor air can be supplied indoors without heating or cooling.Therefore, the thermal labyrinth system can be used in three of the five operation modes.

Verification of the Heating and Cooling Energy Saving Effect of the Thermal Labyrinth System through TRNSYS Simulation
To determine the heating and cooling energy-saving effect when the proposed optimal thermal labyrinth system design and the operational algorithm are applied to the transient system (TRNSYS) energy simulation was performed.The yearly heating and cooling energy consumption with and without the installation of the optimized thermal labyrinth system operated according to the operational algorithm in the target building was compared.

System Description
Figure 12 shows sample building modeling and building description and Figure 13 shows the system chart of the target building.The HVAC system used for air-conditioning of the target zone contains the AHU, to which the thermal labyrinth system is connected, and the ceiling radiation air-conditioning system.This study was conducted to determine the heating and cooling energy-saving effect of the AHU with the use of the thermal labyrinth system, the TRNSYS model was set with the ceiling radiation air-conditioning system already operating.Table 10 shows input conditions for the simulation.Since none of the TRNSYS simulation components can realize the thermal labyrinth system, the temperature data after passing through it were created based on Equation ( 8), and were applied to the TRNSYS using the Type 9 Data Reader.Figure 14 shows TRNSYS modeling of the system with the building.higher than 14.2 °C, the outdoor air would be better used for air-conditioning.Moreover, when the indoor temperature is within the comfortable range of 20-26 °C, the outdoor air can be supplied indoors without heating or cooling.Therefore, the thermal labyrinth system can be used in three of the five operation modes.

Verification of the Heating and Cooling Energy Saving Effect of the Thermal Labyrinth System through TRNSYS Simulation
To determine the heating and cooling energy-saving effect when the proposed optimal thermal labyrinth system design and the operational algorithm are applied to the transient system (TRNSYS) energy simulation was performed.The yearly heating and cooling energy consumption with and without the installation of the optimized thermal labyrinth system operated according to the operational algorithm in the target building was compared.

System Description
Figure 12 shows sample building modeling and building description and Figure 13 shows the system chart of the target building.The HVAC system used for air-conditioning of the target zone contains the AHU, to which the thermal labyrinth system is connected, and the ceiling radiation air-conditioning system.This study was conducted to determine the heating and cooling energy-saving effect of the AHU with the use of the thermal labyrinth system, the TRNSYS model was set with the ceiling radiation air-conditioning system already operating.Table 10 shows input conditions for the simulation.Since none of the TRNSYS simulation components can realize the thermal labyrinth system, the temperature data after passing through it were created based on Equation ( 8), and were applied to the TRNSYS using the Type 9 Data Reader.Figure 14 shows TRNSYS modeling of the system with the building   higher than 14.2 °C, the outdoor air would be better used for air-conditioning.Moreover, when the indoor temperature is within the comfortable range of 20-26 °C, the outdoor air can be supplied indoors without heating or cooling.Therefore, the thermal labyrinth system can be used in three of the five operation modes.

Verification of the Heating and Cooling Energy Saving Effect of the Thermal Labyrinth System through TRNSYS Simulation
To determine the heating and cooling energy-saving effect when the proposed optimal thermal labyrinth system design and the operational algorithm are applied to the transient system (TRNSYS) energy simulation was performed.The yearly heating and cooling energy consumption with and without the installation of the optimized thermal labyrinth system operated according to the operational algorithm in the target building was compared.

System Description
Figure 12 shows sample building modeling and building description and Figure 13 shows the system chart of the target building.The HVAC system used for air-conditioning of the target zone contains the AHU, to which the thermal labyrinth system is connected, and the ceiling radiation air-conditioning system.This study was conducted to determine the heating and cooling energy-saving effect of the AHU with the use of the thermal labyrinth system, the TRNSYS model was set with the ceiling radiation air-conditioning system already operating.Table 10 shows input conditions for the simulation.Since none of the TRNSYS simulation components can realize the thermal labyrinth system, the temperature data after passing through it were created based on Equation ( 8), and were applied to the TRNSYS using the Type 9 Data Reader.Figure 14 shows TRNSYS modeling of the system with the building

Simulation Results
Figure 15 shows the indoor and outdoor air temperatures during heating in January.The indoor set temperature was 20 °C, and the indoor temperature during the working hours was set above 20 °C before and after the thermal labyrinth system was installed.The results showed that the heating system was working normally.It was operated only during the working hours according to the input schedule.The frequency of the heating system operation showed a difference between the two cases.Figure 16 shows that the total accumulated hours when the heater was turned on to keep the indoor temperature at 20 °C decreased by 119 h after the thermal labyrinth system was installed.Moreover, the heating energy consumption decreased by 8% from 400,754 to 369,413 kWh (Table 11).
Figure 17 shows the indoor and outdoor air temperatures during cooling in August.The set indoor temperature during the cooling season was 26 °C.The indoor temperature exceeded 26 °C immediately before and after the cooling system was turned on whether the thermal labyrinth system was installed or not; however, once the cooling system started operating during the working hours, the indoor temperature remained below 26 °C in both cases.This showed that the cooling system was working normally.The number of hours in which the cooling system operated decreased by 214 when the thermal labyrinth system was installed (Figure 18).Moreover, as the thermal labyrinth system was installed, the cooling energy consumption decreased by 83,266 (from 557,966 to 474,700 kWh), which represented 15% of the total cooling energy consumption.As shown

Simulation Results
Figure 15 shows the indoor and outdoor air temperatures during heating in January.The indoor set temperature was 20 • C, and the indoor temperature during the working hours was set above 20 • C before and after the thermal labyrinth system was installed.The results showed that the heating system was working normally.It was operated only during the working hours according to the input schedule.The frequency of the heating system operation showed a difference between the two cases.Figure 16 shows that the total accumulated hours when the heater was turned on to keep the indoor temperature at 20 • C decreased by 119 h after the thermal labyrinth system was installed.Moreover, the heating energy consumption decreased by 8% from 400,754 to 369,413 kWh (Table 11).
Figure 17 shows the indoor and outdoor air temperatures during cooling in August.The set indoor temperature during the cooling season was 26 • C. The indoor temperature exceeded 26 • C immediately before and after the cooling system was turned on whether the thermal labyrinth system was installed or not; however, once the cooling system started operating during the working hours, the indoor temperature remained below 26 • C in both cases.This showed that the cooling system was working normally.The number of hours in which the cooling system operated decreased by 214 when the thermal labyrinth system was installed (Figure 18).Moreover, as the thermal labyrinth system was installed, the cooling energy consumption decreased by 83,266 (from 557,966 to 474,700 kWh), which represented 15% of the total cooling energy consumption.As shown in

Calculation of the Payback Period
The economic feasibility of the optimized thermal labyrinth was analyzed using a simple payback period estimated from the increase in construction costs and saving in annual energy costs.The increased capital costs of the thermal labyrinth was 96,060,000 won, and the reduction cost in the annual energy use was calculated about 9,668,247 won (Table 12).Korea Electric Power Corp. data were used for unit price.Considering the increase in construction cost and the saving in energy cost, the payback period of the optimized thermal labyrinth was determined to be 9.9 years.In previous studies by Song et al. [18] and Paul et al. [19] the payback period of system was 12.1 years and 26 years respectively.This study minimized the increase of construction cost by utilizing basement floor space of existing buildings differently from existing research which required additional construction work.This could reduce the payback period.

Conclusions
In this study, a thermal labyrinth system design that can maximize the heat transfer efficiency was developed, and its optimal operation mode was estimated.The results indicated efficient operation conditions for the thermal labyrinth system.The yearly air-conditioning energy-saving effect of the optimized thermal labyrinth system was also analyzed based on the proposed operation mode.The results of this study can be summarized as follows: (1) Case studies were conducted to increase the efficiency of the thermal labyrinth system.The results showed that the internal space should be relatively wide.Moreover, embossed internal

Calculation of the Payback Period
The economic feasibility of the optimized thermal labyrinth was analyzed using a simple payback period estimated from the increase in construction costs and saving in annual energy costs.The increased capital costs of the thermal labyrinth was 96,060,000 won, and the reduction cost in the annual energy use was calculated about 9,668,247 won (Table 12).Korea Electric Power Corp. data were used for unit price.Considering the increase in construction cost and the saving in energy cost, the payback period of the optimized thermal labyrinth was determined to be 9.9 years.In previous studies by Song et al. [18] and Paul et al. [19] the payback period of system was 12.1 years and 26 years respectively.This study minimized the increase of construction cost by utilizing basement floor space of existing buildings differently from existing research which required additional construction work.This could reduce the payback period.

Conclusions
In this study, a thermal labyrinth system design that can maximize the heat transfer efficiency was developed, and its optimal operation mode was estimated.The results indicated efficient operation conditions for the thermal labyrinth system.The yearly air-conditioning energy-saving effect of the optimized thermal labyrinth system was also analyzed based on the proposed operation mode.The results of this study can be summarized as follows: (1) Case studies were conducted to increase the efficiency of the thermal labyrinth system.The results showed that the internal space should be relatively wide.Moreover, embossed internal surfaces should be formed on the walls and floor to achieve an efficient thermal labyrinth system.(2) The yearly performance analysis of the optimized thermal labyrinth system showed that the COP was highest in summer (August) and winter (January), when there was a marked difference between the outdoor air temperature and the temperature of the air that passed through the thermal labyrinth.On the other hand, in the transitional period (April or October), the outdoor air temperature was similar to that of the air that passed through the thermal labyrinth; thus, there was no active heat transfer, and the COP was low (below 1.0).Additionally, an equation that expresses the relationship between the outdoor air temperature and the temperature of the air that passed through the thermal labyrinth was (3) A yearly operation algorithm was proposed based on the indoor, outdoor air, thermal labyrinth outlet, and ground temperatures, which resulted in five operation modes.Moreover, in three of the five operation modes, air-conditioning based on the thermal labyrinth system was effective.(4) The energy-saving effect from air-conditioning using the optimized thermal labyrinth system on the target building was analyzed under the operation mode.The results showed that in 30% of the air-conditioning load produced by the thermal labyrinth system, the heating and cooling energy consumptions were reduced by 8% and 15%, respectively, which reduced the total air-conditioning energy consumption by about 12%.(5) According to the increase in construction cost and the saving in energy cost, the payback period of the optimized thermal labyrinth was determined to be 9.9 years.
The present study suggests methods for improving the efficiency of the pre-heating and pre-cooling systems using ground heat, a clean type of energy, and management plans.It suggests that the heat transfer efficiency can be greatly improved using the proposed system design when implementing a system based on ground heat.Moreover, it shows that the annual system management energy efficiency can only be increased when various methods of operation are applied according to the external environmental conditions.In addition, it suggests a methodology for analyzing these conditions according to fluid analysis simulations and energy analysis simulations.These results are very important to the future HVAC system development and in studies investigating energy efficiency.

Figure 1 .
Figure 1.Photographs and plans of sample building.(a) Sample building; (b) Building Section; (c) Thermal Labyrinth; (d) Building basement plan.

Figure 1 .
Figure 1.Photographs and plans of sample building.(a) Sample building; (b) Building Section; (c) Thermal Labyrinth; (d) Building basement plan.

Figure 2 .
Figure 2. Earth modeling for thermal labyrinth analysis.

Figure 2 .
Figure 2. Earth modeling for thermal labyrinth analysis.

Figure 6 .
Figure 6.Change in temperature, pressure, and thermal capacity with length.

Figure 6 .
Figure 6.Change in temperature, pressure, and thermal capacity with length.

Figure 6 .
Figure 6.Change in temperature, pressure, and thermal capacity with length.

Figure 6 .
Figure 6.Change in temperature, pressure, and thermal capacity with length.

Figure 6 .
Figure 6.Change in temperature, pressure, and thermal capacity with length.
wall and floor (300 mm depth, 3000 mm gap)2-4Fin on the wall and floor (300 mm depth, 1500 mm gap)

Figure 7 .
Figure 7. Air blower axial power and COP changes with length.

Figure 7 .
Figure 7. Air blower axial power and COP changes with length.

Figure 7 .
Figure 7. Air blower axial power and COP changes with length.
on the wall and floor (300 mm depth, 3000 mm gap)2-4Fin on the wall and floor (300 mm depth, 1500 mm gap)

Figure 8 .
Figure 8. Changes in temperature, pressure, and thermal capacity with the internal form change.

Figure 9 .
Figure 9. Changes in the axial power of the air blower and COP with the internal form change.

Figure 8 .
Figure 8. Changes in temperature, pressure, and thermal capacity with the internal form change.

Figure 9 .
Figure 9. Changes in the axial power of the air blower and COP with the internal form change.

Figure 9 .
Figure 9. Changes in the axial power of the air blower and COP with the internal form change.

Figure 12 .
Figure 12.Sample building modeling and building envelop description.

Figure 13 .
Figure 13.System construction in the studied building.

Figure 12 .
Figure 12.Sample building modeling and building envelop description.

Figure 12 .
Figure 12.Sample building modeling and building envelop description.

Figure 13 .
Figure 13.System construction in the studied building.Figure 13.System construction in the studied building.

Figure 13 .
Figure 13.System construction in the studied building.Figure 13.System construction in the studied building.

2 Figure 14 .
Figure 14.TRNSYS modeling of the system with the building.

Figure 14 .
Figure 14.TRNSYS modeling of the system with the building.

Figure 15 .
Figure 15.Outdoor and indoor air temperature distributions during heating.

Figure 16 .
Figure 16.Heating energy consumption and heating hours.

Figure 17 .
Figure 17.Outdoor and indoor air temperature distributions during cooling.

Figure 15 .
Figure 15.Outdoor and indoor air temperature distributions during heating.

Figure 15 .
Figure 15.Outdoor and indoor air temperature distributions during heating.

Figure 16 .
Figure 16.Heating energy consumption and heating hours.

Figure 17 .
Figure 17.Outdoor and indoor air temperature distributions during cooling.

Figure 16 .
Figure 16.Heating energy consumption and heating hours.

Figure 15 .
Figure 15.Outdoor and indoor air temperature distributions during heating.

Figure 16 .
Figure 16.Heating energy consumption and heating hours.

Figure 17 .
Figure 17.Outdoor and indoor air temperature distributions during cooling.Figure 17.Outdoor and indoor air temperature distributions during cooling.

Figure 17 .
Figure 17.Outdoor and indoor air temperature distributions during cooling.Figure 17.Outdoor and indoor air temperature distributions during cooling.

Figure 18 .
Figure 18.Cooling energy consumption and cooling hours.

Figure 18 .
Figure 18.Cooling energy consumption and cooling hours.

Table 1 .
Conditions of the studied building and thermal labyrinth system.

South Korea HVAC ** zone Information Building type Office
SystemAHU * with thermal labyrinth, radiant ceiling

Table 1 .
Conditions of the studied building and thermal labyrinth system.

South Korea HVAC ** zone Information Building type Office
System AHU * with thermal labyrinth, radiant ceiling Area 605.55 m 2 Ventilation requirement 7400 m 3 /h Minimum outdoor air requirement 2200 m 3 /h * AHU: Air handling unit ** HVAC: Heating, ventilation, and air conditioning *** COP: Coefficient of performance.

Table 2 .
Cases set considered for the simulation.

Table 2 .
Cases set considered for the simulation.

Table 3 .
Coefficients used in the calculation of the ground temperature distribution.

Table 4 .
Properties of the material applied to the simulation.

Table 3 .
Coefficients used in the calculation of the ground temperature distribution.

Table 4 .
Properties of the material applied to the simulation.
Air flow rate 2200 m 3 /h Figure 4. Scene with its mesh.

Table 7 .
The studied cases for analyzing the impact of the change in the length on the heat transfer efficiency.

Table 7 .
The studied cases for analyzing the impact of the change in the length on the heat transfer efficiency.
Air blower axial power and COP changes with length.2.3.2.Assessment of the Heat Transfer Efficiency According to the Internal Shape

Table 8 .
The studied cases for analyzing the impact of the internal form change on the heat transfer efficiency.Air blower axial power and COP changes with length.2.3.2.Assessment of the Heat Transfer Efficiency According to the Internal Shape

Table 9 .
The optimized conditions of the thermal labyrinth design.Changes in temperature, pressure, and thermal capacity with the internal form change.

Table 9 .
The optimized conditions of the thermal labyrinth design.

Table 9 .
The optimized conditions of the thermal labyrinth design.
Figure 10.Monthly average outdoor air (OA) temperature, optimized thermal labyrinth temperature, and COP.

Table 10 .
Input conditions of TRNSYS simulation.

Table 10 .
Input conditions of TRNSYS simulation.

Table 11 ,
the yearly total air conditioning energy consumption decreased by about 12% from 958,721 to 844,114 kWh.

Table 11 ,
the yearly total air conditioning energy consumption decreased by about 12% from 958,721 to 844,114 kWh.
Energies 2018, 11, x FOR PEER REVIEW 14 of 17 in Table 11, the yearly total air conditioning energy consumption decreased by about 12% from 958,721 to 844,114 kWh.
Energies 2018, 11, x FOR PEER REVIEW 14 of 17 in Table 11, the yearly total air conditioning energy consumption decreased by about 12% from 958,721 to 844,114 kWh.

Table 11 .
Energy saving ratio for air conditioning.

Table 12 .
Reduction in the annual energy costs.

Table 11 .
Energy saving ratio for air conditioning.

Table 12 .
Reduction in the annual energy costs.