Mathematical Modelling and Operational Analysis of Combined Vertical–Horizontal Heat Exchanger for Shallow Geothermal Energy Application in Cooling Mode

: Geothermal heat exchangers (GHEs) represent a buried pipe system, which can be utilised to harness renewable thermal energy stored in the ground to improve the e ﬃ ciency of heating and cooling systems. Two basic arrangements of GHEs have been widely used: vertical and horizontal. Vertical GHEs generally have a better performance in comparison with the horizontal arrangement, and these systems are particularly suitable for conﬁned spaces. Nevertheless, the main technical challenge associated with GHEs, for either the vertical or the horizontal arrangement, is the performance deterioration associated with an increase in the operation times during summer or winter seasons. In this paper, a combined horizontal-vertical GHE arrangement is proposed to address the current challenges. The combined GHE arrangement can be operated in ﬁve di ﬀ erent modes, corresponding to di ﬀ erent thermal loading conditions. These ﬁve operation modes of the combined GHE are analysed based on the transient ﬁnite di ﬀ erence models previously developed for the horizontal and vertical arrangements. The simulation results reveal that for the single operation mode (horizontal or vertical only), the vertical GHE performs better than the horizontal GHE due to relatively stable ground temperature deep down. While, for the combined operation mode, the series operations (horizontal to vertical or vertical to horizontal) of the GHE are superior to the split mode. It is found that the e ﬀ ect of the ﬂuid mass ﬂow rate ratio is trivial on the heat dissipation of the split mode GHE. The highest heat transfer rate in the split ﬂow operational mode is rendered by the ratio of the mass ﬂow rate of 40% horizontal and 60% vertical. In addition, the climate condition has more e ﬀ ect on GHE’s performance and the increase of the ﬂuid ﬂow rate it can enhance the amount of energy released by the GHE. as the contact area, where heat was exchanged with the surrounding soil, increased. It is seen that the series operations could release 40.8% and 39.6% more energy than the vertical mode, for the horizontal to vertical and vice versa modes, respectively. In the split ﬂow mode, the ﬂuid is split to ﬂow with a ratio of 50% in the horizontal GHE and 50% in the vertical GHE. The results demonstrate that the split ﬂow mode released 2.8% less energy than the horizontal to vertical mode. Even though, the amount of energy released in the intermittent operation regime was less than that generated by the continuous operation due to the total operation period in the intermittent condition was less compared with the continuous condition. However, the average heat transfer rates increased. They were 60.1% for the horizontal mode, 68.5% for the vertical mode, 54% for the horizontal to vertical mode, 53.5% for the vertical to horizontal mode, and 52.6% for the split ﬂow mode. In the intermittent operation condition, the deterioration of the ground temperature during the operation hours is possible to recover during the time when the system is switched o ﬀ . As a result, it increases the heat transfer rate of the GHE and produces lower outlet ﬂuid temperature over the next day’s operation as reﬂected by an example of ﬂuid temperature generated by the vertical to horizontal mode (see Figure 11).


Introduction
Increasing energy demands and crude oil depletion have forced the world to explore new energy resources, those that are not only sustainable but also have a minor environmental footprint. Shallow geothermal energy is one of the sustainable energy resources that could be applied to improve the efficiency of current and future heating and cooling systems. This renewable thermal energy is usually exploited in a geothermal heat exchanger (GHE); a system of pipes buried in the ground, in which the heat carrier fluid is circulated. GHEs are normally integrated with heat pumps and air Studies on hybrid ground source heat pump systems with a cooling tower as a supplemental heat rejecter were presented by a number of researchers [9][10][11][12][13]. Park et al. [9] proposed a new parallel system comprising of a hybrid ground source heat pump-cooling tower. The new parallel system enables the GHE to be switched off during the recovery period of the soil thermal condition. The performance of the heat pump was investigated at different flow rates of the fluid in the primary flow loop (the heat pump), GHE, and cooling tower. The results indicated that the parallel system of the ground source heat pump-cooling tower generates 21% COP more than that produced by the conventional ground source heat pump system. Man et al. [10] provided the technical and economic analysis of a hybrid cooling tower-ground source heat pump, based on the hourly load of a two storey residential building located in Hong Kong. The hybrid system in this study not only solves the thermal degradation problem, but also reduces the operation and capital costs of the air conditioning system. A study of the operation strategy of a hybrid cooling tower-ground source heat pump system was presented by Wang et al. [11]. The operation strategy consists of a fixed cooling set point, outside air reset, wet bulb reset, and load reset. Fan et al. [13] presented a theoretical design of a hybrid cooling tower-ground source heat pump, which takes into account the effect of borehole distance, borehole depth, and thermal properties of the grout. In this study, a combined strategy of operation was introduced. The results showed that the lowest energy consumption was obtained when the control strategies of the entering water temperature and wet-bulb temperature differences are combined.
Canelli et al. [14], presented an analysis of the energy, economic, and environmental performances of three different hybrid ground source heat pump systems including (1) a hybrid boiler-chiller-ground source heat pump, (2) a hybrid boiler-chiller-ground source heat pump and fuel cell, and (3) a hybrid boiler-chiller-ground source heat pump and photovoltaic thermal system. The system was optimised to meet the heating and cooling conditions of both residential and commercial buildings, which are in a sharing load. The results indicated that the hybrid system with the fuel cells and photovoltaic thermal system has a definite advantage in terms of energy savings, operational costs, and carbon emission reductions.
A detailed numerical study was conducted by Zhu et al. [15] to investigate the performance and economic characteristics of a combined vertical ground source heat pump and phase change material cooling storage system. The optimal performance of the hybrid system was achieved by varying the ratio of the phase change material cooling storage system to the total cooling load of the system. The obtained optimal cooling storage ratio was close to 40%.
A hybrid system of the ground source electrical heat pump and ground source absorption heat pump was proposed by Wu et al. [16]. The motivation for this study was to combine the features of both heat pumps, as the ground source electrical heat pump has higher energy efficiency in the cooling mode. Conversely, the ground source absorption heat pump has higher energy efficiency in the heating mode. This study also presented the effects of supply ratios on thermal imbalance ratios, annual primary energy efficiency, and cost-efficiency characteristics.
A theoretical study on a hybrid air source heat compensator-vertical ground source heat pump system was conducted by You et al. [17], using a simulation tool, TRNSYS. Four operation strategies were analysed including an air source heat compensator for direct heat compensation, a combined air source heat compensator-ground source heat pump for heat compensation, a combined air source heat compensator-ground source heat pump for space heating, and an air source heat compensator-ground source heat pump for domestic hot water. The results showed that the hybrid air source heat compensator-ground source heat pump reduces energy consumption by 23.86% compared with the boiler-split air conditioner system. Additionally, the operational costs are reduced by 50%. Another study has done by Xianting Li et al. [18] by developing an energy-efficient heat pump system using the TRNSYS simulation tool in a combination of a hybrid source and three types of hybrid source heat pumps to calculate the efficiency of year-round operation. The results show that the hybrid heat pump system can maintain heat reliability during year-round operation including winter. This hybrid source heat pump system saves energy up to 15% and the payback period is approximately five years. In addition, Gaoyang Hou et al. [19] also conducted a simulation using the TRNSYS tool to analyse the work system of heat pumps sourced from hybrid soil with optimal control strategies. His research combines a horizontal ground loop and a liquid dry cooler with a short and long term simulation process on TRNSYS to analyse soil thermal conditions and energy variations. Gaoyang Hou et al. also performed a simulation using the TRNSYS tool to analyse the heat pump working system sourced from hybrid soil with optimal control strategies. Their research combines a horizontal ground loop and a liquid dry cooler with a short and long term simulation process on TRNSYS to analyse soil thermal conditions and energy variations. The simulation results show that the overall performance in the short-term simulation is influenced mainly by the temperature of the diverter heating set rather than by cooling. Diverter heating set was recommended about 8-10 • C in climatic zones similar to the Birmingham area by combining long-term coefficient of performance (COP) values and soil thermal variation.
A hybrid active air source regeneration-ground source heat pump system was studied by Allaerts et al. [20]. In this study, the borehole area is divided into two different regions, namely warm and cold regions. These two different regions were proposed to balance the extraction/rejection of heat during heating and cooling periods. In addition, a supplementary dry cooler was used to capture heat/cold during summer/winter to recover the degradation of the soil thermal condition. According to Allaerts et al. the proposed hybrid system can significantly reduce the size of the borehole area by up to 47% in the cost-optimal configuration.
The brief review of the literature  above shows that most researchers proposed a hybrid cooling tower-ground source heat pump for a cooling load dominated regime and a hybrid solar system-ground source heat pump for heating load dominated conditions. Some other studies proposed a combined system of a ground source heat pump with different additional heat rejecter/absorber systems including a boiler, chiller, fuel cell, photovoltaic, and phase change material cooling storage systems. Additionally, a ground source heat pump system with two different regions of BHEs including warm and cold was recently suggested and analysed. However, it seems a comprehensive study of a combined GHE arrangement with different operation modes has yet to be presented. This paper is focused on the operational analysis of a combined horizontal-vertical GHE to address various demands and thermal loading conditions. The performance of the combined GHE is studied based on the results of a transient finite difference model developed in this paper. The effects of continuous and intermittent operation conditions, climate condition and fluid mass flow rate on the GHE's performance are investigated with this new model. Based on the analysis of the outcomes of numerical simulations, the recommendations for the optimum operation of combined GHEs are summarised in the conclusion.

Physical Model
The physical model is described first including a number of simplifications, which are usually needed for quantitative analysis of complex systems or phenomena. A horizontal GHE, as shown in Figure 1a, represents a multi-U-shaped pipe buried in a trench at a specific depth h close to the ground surface. A horizontal GHE has a total length L, diameter d, wall thickness t, and pipe spacing L s . The heat transfer process at the horizontal GHE occurs when the fluid at a temperature T fi enters the GHE and exchanges the heat energy with the pipe's inner surface by convection. The heat then conducts through the pipe wall to the surrounding soil. The ground exchanges the heat at the ground surface with the atmosphere through diverse heat transfer mechanisms including reflection, convection, radiation, and evaporation. A vertical GHE, as shown in Figure 1b, could be modelled as a single U-shaped pipe, which is buried inside a borehole having a diameter D b and depth z. The pipe is described by diameter d, length L, thickness t, and shank spacing L s . The borehole is usually filled with a material, which has a relatively higher thermal conductivity than the surrounding soil, in order to enhance the GHE's thermal performance. The borehole could be confined by different types of soil/rock layers that have different thermal conductivities. The heat transfer at GHE occurs when the fluid is operating at a temperature of T fi and the mass flow rate m is circulated through a U-shaped pipe, exchanges its heat with the pipe's wall by convection. The heat then flows through the pipe's wall, grout, and surrounding soil by conduction sequentially. The convection occurs when the soil surface exchanges the heat with the atmosphere. As mentioned above, a vertical GHE is appropriate for installation in areas with limited land access. It may benefit from the relatively stable ground temperature at the deeper ground depth, and it leads to an improvement in performance. However, this performance decreases with the increase in the operation time, and it may take a season to recover fully. The poor soil conductivity and relatively deep depth of burring are the main reasons for the slow recovery of the soil temperature.
Energies 2020, 13, x FOR PEER REVIEW 5 of 21 when the fluid is operating at a temperature of Tfi and the mass flow rate m is circulated through a U-shaped pipe, exchanges its heat with the pipe's wall by convection. The heat then flows through the pipe's wall, grout, and surrounding soil by conduction sequentially. The convection occurs when the soil surface exchanges the heat with the atmosphere. As mentioned above, a vertical GHE is appropriate for installation in areas with limited land access. It may benefit from the relatively stable ground temperature at the deeper ground depth, and it leads to an improvement in performance. However, this performance decreases with the increase in the operation time, and it may take a season to recover fully. The poor soil conductivity and relatively deep depth of burring are the main reasons for the slow recovery of the soil temperature. To take advantage of different arrangements and negate their disadvantages, a combined, structured, horizontal-vertical GHE was proposed, see  To take advantage of different arrangements and negate their disadvantages, a combined, structured, horizontal-vertical GHE was proposed, see

Mathematical Model
The performance of the combined GHE was studied based on the previously developed mathematical models for the horizontal [21] and vertical [22] GHE, respectively. The key equations in these models are briefly summarised below.

Governing Equation of the Soil
The thermal model of the soil domain around the horizontal GHE is developed based on the Cartesian coordinates as it is easy to handle when applying the soil's internal source term, in which the value varies with the increase in the soil depth. The governing equation of the soil domain around the horizontal GHE is given as: The thermal model of the soil domain around the vertical GHE is developed by considering the cylindrical coordinates. The cylindrical coordinates are selected because these represent the shape of the borehole, thus the boundary conditions can be easily defined and prescribed. As a result, an accurate result can be obtained. The governing equation of the soil domain around the vertical GHE is given as An internal heat source term was used to take into account the effect of seasonal changes in soil temperature [21,23], expressed by the following equation:

Mathematical Model
The performance of the combined GHE was studied based on the previously developed mathematical models for the horizontal [21] and vertical [22] GHE, respectively. The key equations in these models are briefly summarised below.

Governing Equation of the Soil
The thermal model of the soil domain around the horizontal GHE is developed based on the Cartesian coordinates as it is easy to handle when applying the soil's internal source term, in which the value varies with the increase in the soil depth. The governing equation of the soil domain around the horizontal GHE is given as: The thermal model of the soil domain around the vertical GHE is developed by considering the cylindrical coordinates. The cylindrical coordinates are selected because these represent the shape of the borehole, thus the boundary conditions can be easily defined and prescribed. As a result, an accurate result can be obtained. The governing equation of the soil domain around the vertical GHE is given as An internal heat source term was used to take into account the effect of seasonal changes in soil temperature [21,23], expressed by the following equation: The direct measurement or analytical approach can be used to determine deviations in soil temperature during the process of seasonal change. The experimental results are preferred as they represent the actual soil temperature. However, it is sometimes hard to obtain the measurement Energies 2020, 13, 6598 7 of 20 data, especially for a specific location and a certain depth. Thus, the analytical equation presented by Baggs [24] was used in this work:

Governing Equation of the Temperature Exchange in the Grout for the Vertical GHE
The governing equation of the grout is given by a similar equation to the one above:

Governing Equation of the Pipe
The energy balance equation for the pipe is The area of the pipe for the vertical GHE was determined based on an equivalent pipe diameter, which is given later by Equation (8).

Governing Equation of the Fluid
For the working fluid, the energy balance equation states Equation (6) can be applied directly to the evaluation of performance of horizontal GHE. While for the vertical GHE, the internal pipe diameter d in was replaced with an equivalent pipe diameter D eq since the single U shaped pipe is represented as a pipe, which has an equivalent diameter (refer to Figure 3). The equation to calculate the equivalent diameter has been developed and validated by Gu and Oneal [25] as: Energies 2020, 13, x FOR PEER REVIEW 7 of 21 The direct measurement or analytical approach can be used to determine deviations in soil temperature during the process of seasonal change. The experimental results are preferred as they represent the actual soil temperature. However, it is sometimes hard to obtain the measurement data, especially for a specific location and a certain depth. Thus, the analytical equation presented by Baggs [24] was used in this work:

Governing Equation of the Temperature Exchange in the Grout for the Vertical GHE
The governing equation of the grout is given by a similar equation to the one above:

Governing Equation of the Pipe
The energy balance equation for the pipe is The area of the pipe for the vertical GHE was determined based on an equivalent pipe diameter, which is given later by Equation (8).

Governing Equation of the Fluid
For the working fluid, the energy balance equation states Equation (6) can be applied directly to the evaluation of performance of horizontal GHE. While for the vertical GHE, the internal pipe diameter was replaced with an equivalent pipe diameter since the single U shaped pipe is represented as a pipe, which has an equivalent diameter (refer to Figure 3). The equation to calculate the equivalent diameter has been developed and validated by Gu and Oneal [25] as: The equation is valid if ≤ ≤ Energies 2020, 13, 6598 8 of 20 The boundary conditions of the computational domain of the horizontal GHE are: • The boundary conditions of the computational domain of the vertical GHE are: The system of differential Equations (1), (2), and (4)-(6) were solved by using an explicit finite difference scheme, which provided the accurate solution for the temperature distribution in the soil, grout, pipe, and fluid. As an example, Figure 4 shows the computational domain of the vertical GHE. The working fluid exchanges the heat energy with the pipe's surface, which then flows through the pipe wall in a 1D way. The heat is then transferred through the grout and surrounding soil in a 2D manner. Two different boundary conditions are applied at the boundary of the soil domain. The convection heat transfer is considered at the ground surface, while the adiabatic conditions are applied at the bottom and lateral edges of soil domain boundaries.  To achieve the stability of the explicit numerical calculations, the time step ∆ must be within the Courant-Friedrichs-Lewy stability range [26], which is given by: Total of heat exchange by the GHE was calculated in accordance with the following equation: To achieve the stability of the explicit numerical calculations, the time step ∆t must be within the Courant-Friedrichs-Lewy stability range [26], which is given by: Total of heat exchange by the GHE was calculated in accordance with the following equation: For the split flow mode, the final temperature of mixed fluid from both the horizontal and the vertical GHE was found as: The described GHE's models were validated against the experimental data [21][22][23] and will be used to simulate the performance of a combined GHE.

A Hypothetical Case Study
In this study, each arrangement of the GHEs (horizontal or vertical) was set to have the same parameters including the pipe length, pipe diameter, and fluid flow rates. The soil domain around the GHEs was assumed to be a single soil layer, which was homogeneous and isotropic. Table 1 presents the parameters of the reference case used in the simulation to be described next. The initial temperature of the ground was estimated using the analytical equation (4) suggested by Baggs [24]. The Baggs' equation was determined based on the input parameters of the climate conditions data. The amplitude data for annual air temperature and the average air temperature was obtained from the Australian Bureau of Meteorology [29]. The variable ground temperature for the local site was based on the data given by Baggs [24]. Table 2 summarises the parameters of the reference used to estimate the soil temperature conditions. Figure 5 shows the typical changes of soil temperatures in Adelaide and Brisbane at the end of winter (in August) and summer (in February). As observed from the figure that the soil temperature in Adelaide was lower than in the Brisbane. This tendency could be affected by mild, and generally a warm and temperate climate of Adelaide compared to humid subtropical climate of Brisbane as represented by the parameter of average annual air temperature (T m ) in Equation (11). Besides the local site variable of the ground temperature (∆T m ), which varied with geographic location, it also affects the soil temperature. As an example, Figure 6 shows the ambient temperature in Adelaide and Brisbane on three consecutive summer days. These were used in the case study. In general, it can be seen from the figure that the pattern of air temperature was almost the same. The profile of air temperature fluctuated during day and night. For these three randomly selected summer days it was found that the air temperature in Adelaide was slightly higher than that in Brisbane.    The effect of seasonal changes on the soil temperature is modelled by incorporating the internal source term concept into the GHE model [21]. The value of the internal source term varied with the soil depths. It was higher in a shallow region and lower in a deeper zone. At a depth of 12 m below the ground surface, the value of the internal source term was assumed to be zero (this corresponded to the depth at which the effect of the ambient temperature on the soil temperature was negligible). The effect of seasonal changes on the soil temperature is modelled by incorporating the internal source term concept into the GHE model [21]. The value of the internal source term varied with the soil depths. It was higher in a shallow region and lower in a deeper zone. At a depth of 12 m below the ground surface, the value of the internal source term was assumed to be zero (this corresponded to the depth at which the effect of the ambient temperature on the soil temperature was negligible). Figure 7 shows the absolute value of the internal source term at various depths for two reference locations, namely Adelaide and Brisbane. As is shown in Figure 7, at the upper layer (0-3.5 m depth), the internal source term of the soil in Brisbane was relatively higher than in Adelaide. However, at a deeper layer, the condition was in the opposite. This phenomenon could be affected by meteorological, terrain, and subsurface conditions. The value of the internal heat source was positive when the ground temperature was warm (August to February) and conversely it was negative when the ground temperature was cool (February to August), as is shown in Figure 5.
Energies 2020, 13, x FOR PEER REVIEW 12 of 21 Figure 7 shows the absolute value of the internal source term at various depths for two reference locations, namely Adelaide and Brisbane. As is shown in Figure 7, at the upper layer (0-3.5 m depth), the internal source term of the soil in Brisbane was relatively higher than in Adelaide. However, at a deeper layer, the condition was in the opposite. This phenomenon could be affected by meteorological, terrain, and subsurface conditions. The value of the internal heat source was positive when the ground temperature was warm (August to February) and conversely it was negative when the ground temperature was cool (February to August), as is shown in Figure 5.

Results and Discussion
This section shows the simulation outcomes of the GHE subjected to five operational modes, as described above. The analysis and discussion on continuous operation, intermittent operation, split flow operation, climate condition, and variations of the fluid mass flow rate will be discussed in the following section. Figure 8 shows the simulation results of the GHE, which was operated on three consecutive

Results and Discussion
This section shows the simulation outcomes of the GHE subjected to five operational modes, as described above. The analysis and discussion on continuous operation, intermittent operation, split Energies 2020, 13, 6598 12 of 20 flow operation, climate condition, and variations of the fluid mass flow rate will be discussed in the following section. Figure 8 shows the simulation results of the GHE, which was operated on three consecutive summer days, with respect to the geological and climate conditions corresponding to Adelaide (South Australia). The amount of energy released by the GHE was calculated using Equation (21). It is observed that the lowest value of the outlet fluid temperature could be attained so the highest amount of energy can be released. As an example, Figure 9 provides the profile of the fluid temperature generated by operation of the horizontal, vertical, and the horizontal to the vertical GHE. This figure demonstrated that the outlet fluid temperature increased with the increase in the operation period. The accumulation of heat in the surrounding soil during the operation of the GHE led to a reduction in the heat transfer rate. As is shown in the Figure 9, the profile of the outlet fluid temperature of the vertical GHE was relatively stable when compared with those generated by the horizontal and the series operation modes. This phenomenon occurred because the ambient temperature, which fluctuated diurnally, did not have a significant effect on the performance of the vertical GHE. From Figure 8, it followed that the lowest energy was released when operating the horizontal GHE only. The horizontal GHE released 1002 kW less energy than that released in 3 days by the vertical GHE. The relatively stable temperature of the ground at a deeper layer might enhance the heat transfer capacity of the vertical GHE and led to better thermal performance. The amount of energy released could be increased by using the combined operation mode of the GHE, including the split flow and series operation modes. It can be seen that the series operation mode, from horizontal to vertical, could release slightly higher energy, namely 22 MJ more than the opposite operation mode. While the split flow mode, with a ratio of fluid flow at 50% in the horizontal GHE and 50% in the vertical GHE, released less energy than the series operation mode. The average heat transfer rate was calculated by dividing the amount of energy released with the operation period, namely 72 h. The single arrangement of the horizontal or the vertical GHE generated an average heat transfer rate of 6.3 kW and 10.2 kW, respectively. The combined arrangements produced a relatively higher heat transfer rate namely: 15.7 kW, 15.6 kW, and 15.4 kW for operation modes from the horizontal to the vertical, from the vertical to the horizontal GHE, and by splitting the fluid flow, respectively. Therefore, the horizontal or vertical GHE may be operated when the loading load was relatively low. At peak loads or when heating/cooling demands were relatively high, the combined GHE's operation could have a significant advantage, see Figure 8. respectively. Therefore, the horizontal or vertical GHE may be operated when the loading load was relatively low. At peak loads or when heating/cooling demands were relatively high, the combined GHE's operation could have a significant advantage, see Figure 8.

Intermittent Operation
The GHE can also be operated in the intermittent condition to cope with cyclic load conditions. Figure 10 shows the amount of energy released by the GHE operating intermittently in three summer days, see Figures 5 and 6. The GHE ran for 8 h (during the working hours) and was off for 16 h daily. The results display that during 24 h of operation, the horizontal GHE released 41% less

Intermittent Operation
The GHE can also be operated in the intermittent condition to cope with cyclic load conditions. Figure 10 shows the amount of energy released by the GHE operating intermittently in three summer days, see Figures 5 and 6. The GHE ran for 8 h (during the working hours) and was off for 16 h daily. The results display that during 24 h of operation, the horizontal GHE released 41% less energy than the vertical GHE. The combined operation modes of the GHE could increase the amount of energy released, as the contact area, where heat was exchanged with the surrounding soil, increased. It is seen that the series operations could release 40.8% and 39.6% more energy than the vertical mode, for the horizontal to vertical and vice versa modes, respectively. In the split flow mode, the fluid is split to flow with a ratio of 50% in the horizontal GHE and 50% in the vertical GHE. The results demonstrate that the split flow mode released 2.8% less energy than the horizontal to vertical mode. Even though, the amount of energy released in the intermittent operation regime was less than that generated by the continuous operation due to the total operation period in the intermittent condition was less compared with the continuous condition. However, the average heat transfer rates increased. They were 60.1% for the horizontal mode, 68.5% for the vertical mode, 54% for the horizontal to vertical mode, 53.5% for the vertical to horizontal mode, and 52.6% for the split flow mode. In the intermittent operation condition, the deterioration of the ground temperature during the operation hours is possible to recover during the time when the system is switched off. As a result, it increases the heat transfer rate of the GHE and produces lower outlet fluid temperature over the next day's operation as reflected by an example of fluid temperature generated by the vertical to horizontal mode (see Figure 11). flow mode. In the intermittent operation condition, the deterioration of the ground temperature during the operation hours is possible to recover during the time when the system is switched off. As a result, it increases the heat transfer rate of the GHE and produces lower outlet fluid temperature over the next day's operation as reflected by an example of fluid temperature generated by the vertical to horizontal mode (see Figure 11).

Split Flow Operation
In this section, the effect of the fluid mass flow rate ratio in the split flow operation mode was investigated. The ratio of the fluid mass flow rate varied as follows: 30%:70%; 40%:60%; 50%:50%; 60%:40%; and 70%:30% for the horizontal and the vertical GHE, respectively. From the current numerical simulations, it was found that the ratio of fluid mass flow rate in the split flow mode did not significantly affect the amount of energy released by the GHE, as shown in Figure 12. As an example, Figure 13, presents the outlet fluid temperature of the GHE at three different flow rate ratios, namely: 30%:70%; 50%:50%; and 70%:30% for the horizontal and the vertical GHE, respectively. It is found that the difference in the outlet fluid temperature was relatively small. The highest fluid outlet temperature was yielded by the GHE with a flow rate ratio of 70%:30%. It can be seen that the GHE operated with a flow rate ratio of 70% horizontal and 30% vertical released the lowest amount of energy across three consecutive days of operation, namely, 3% less than the one with a flow rate ratio of 60% horizontal and 40% vertical. The amount of energy released gradually increased with the increase of the ratio of the fluid mass flow rate in the vertical GHE and declined at the ratio of 30% horizontal and 70% vertical GHE namely, 0.7% less than that with a ratio of 40% horizontal and 60% vertical GHE. This tendency was due to a significant reduction in the fluid mass flow rate in the horizontal GHE (only 30% of the total fluid flow rate). At a lower mass flow rate, the thermal resistance of the horizontal GHE increased as a reduction in fluid mass flow rate was directly proportional to the decrease in the convective heat transfer coefficient between the working fluid and the inner pipe surface. As a result, it decreased the heat transfer rate of the horizontal GHE. In addition, with 70% of mass flow rate in the vertical GHE, the surrounding soil temperature deteriorated quickly due to heat accumulation, leading to degradation in the heat transfer capacity of the vertical GHE. From the Figure 12, it can be seen that the highest heat transfer rate was produced by the GHE with a ratio of flow rate of 40% horizontal and 60% vertical namely, 15,490 kW.

Climate Condition
The performance of the GHE installed in a temperate climate corresponding to Adelaide was compared with that installed in a subtropical climate, such as near Brisbane. Figure 14 summarises the performance of the GHE, which was installed in two different regions, with different climate conditions. It was found that during 72 h (3 days: 10-12 January) of the continuous operation, the GHE installed in Adelaide's temperate climate could release 34.3% and 26.9% more energy compared with that installed in Brisbane, for the horizontal and vertical mode, respectively. While

Climate Condition
The performance of the GHE installed in a temperate climate corresponding to Adelaide was compared with that installed in a subtropical climate, such as near Brisbane. Figure 14 summarises the performance of the GHE, which was installed in two different regions, with different climate conditions. It was found that during 72 h (3 days: 10-12 January) of the continuous operation, the GHE installed in Adelaide's temperate climate could release 34.3% and 26.9% more energy compared with that installed in Brisbane, for the horizontal and vertical mode, respectively. While

Climate Condition
The performance of the GHE installed in a temperate climate corresponding to Adelaide was compared with that installed in a subtropical climate, such as near Brisbane. Figure 14 summarises the performance of the GHE, which was installed in two different regions, with different climate conditions. It was found that during 72 h (3 days: 10-12 January) of the continuous operation, the GHE installed in Adelaide's temperate climate could release 34.3% and 26.9% more energy compared with that installed in Brisbane, for the horizontal and vertical mode, respectively. While for the combined modes including horizontal to vertical, vertical to horizontal, and split flow were 31.7%, 31.8%, and 31% higher, respectively. It is observed that operating the vertical GHE only, in the temperate climate, Adelaide, could release an amount of energy that was almost the same as that released by the combined operation mode, in the subtropical climate. This tendency occurred because the climate condition influenced the initial soil temperature in both regions. The initial soil temperature in Adelaide was lower than that in Brisbane, as seen in Figure 5. In addition, the difference in soil types in both locations contributed to the heat transfer capacity of the GHE. As a result, the GHE installed in a temperate climate, Adelaide, produced a lower outlet fluid temperature as demonstrated in Figure 15.        Figure 16 shows the effect of the fluid mass flow rate on the amount of energy released by the GHE in 3 summer days in Adelaide. The performance of the GHE under two different fluid mass flow rates namely, 0.6 kg/s and 1 kg/s, was compared. The results of the simulation indicate that the amount of energy released by the GHE increased as the fluid flow rate increased. This tendency occurred because the mass flow rate affected the convective heat transfer coefficient of the fluid inside the GHE. The higher the mass flow rate, the higher the coefficient convective heat transfer was attained. As a result, it increased the heat transfer capacity and the amount of energy released by the GHE. As an example, Figure 17 shows the profile of the fluid temperature generated from the horizontal, vertical, and horizontal to vertical modes at two different mass flow rates namely, 0.6 kg/s and 1 kg/s. This figure shows that the outlet fluid temperature increased with the increase of fluid mass flow rate. This tendency occurred because at a higher mass flow rate, the time period during which the fluid makes contact with the pipe was shorter when compared to a relatively low flow rate. In addition, at a higher mass flow rate, the GHE released more heat into the surrounding soil and led to a quick increase in soil temperature, resulting in increased outlet fluid temperature. Varying the mass flow rate from 0.6 to 1 kg/s increased the amount of energy released by 2.4% for the horizontal mode and 3.3% for the vertical mode. The rate increases were 4.9%, 5.1%, and 5.6% for the horizontal to vertical, vertical to horizontal, and split flow mode, respectively.

Variations of the Fluid Mass Flow Rate
Energies 2020, 13, x FOR PEER REVIEW 18 of 21 attained. As a result, it increased the heat transfer capacity and the amount of energy released by the GHE. As an example, Figure 17 shows the profile of the fluid temperature generated from the horizontal, vertical, and horizontal to vertical modes at two different mass flow rates namely, 0.6 kg/s and 1 kg/s. This figure shows that the outlet fluid temperature increased with the increase of fluid mass flow rate. This tendency occurred because at a higher mass flow rate, the time period during which the fluid makes contact with the pipe was shorter when compared to a relatively low flow rate. In addition, at a higher mass flow rate, the GHE released more heat into the surrounding soil and led to a quick increase in soil temperature, resulting in increased outlet fluid temperature.
Varying the mass flow rate from 0.6 to 1 kg/s increased the amount of energy released by 2.4% for the horizontal mode and 3.3% for the vertical mode. The rate increases were 4.9%, 5.1%, and 5.6% for the horizontal to vertical, vertical to horizontal, and split flow mode, respectively.   Energies 2020, 13, x FOR PEER REVIEW 18 of 21 attained. As a result, it increased the heat transfer capacity and the amount of energy released by the GHE. As an example, Figure 17 shows the profile of the fluid temperature generated from the horizontal, vertical, and horizontal to vertical modes at two different mass flow rates namely, 0.6 kg/s and 1 kg/s. This figure shows that the outlet fluid temperature increased with the increase of fluid mass flow rate. This tendency occurred because at a higher mass flow rate, the time period during which the fluid makes contact with the pipe was shorter when compared to a relatively low flow rate. In addition, at a higher mass flow rate, the GHE released more heat into the surrounding soil and led to a quick increase in soil temperature, resulting in increased outlet fluid temperature.
Varying the mass flow rate from 0.6 to 1 kg/s increased the amount of energy released by 2.4% for the horizontal mode and 3.3% for the vertical mode. The rate increases were 4.9%, 5.1%, and 5.6% for the horizontal to vertical, vertical to horizontal, and split flow mode, respectively.