Exploratory Research to Improve Energy-E ﬃ ciency of a Ground-Coupled Heat Pump Utilizing an Automatic Control Device of Circulation Pump Speed

: Ground-coupled heat pumps (GCHPs) are an e ﬃ cient thermal energy production system that can satisfy the gap between heating and air-conditioning. Be that as it may, exploratory research on GCHPs is still lacking. The ﬁrst objective of this article is to describe a utilitarian energy-e ﬃ ciency improvement device for a vertical GCHP system that includes a bu ﬀ er tank (BT) between the heat pump unit and the fan coil units and user supply, utilizing the quantitative regulation of water ﬂow rate with a variable-speed circulation pump. At that point, the investigative estimations are utilized to test the performances of the GCHP system in various operating modes. Fundamental e ﬃ ciency parameters (coe ﬃ cient of performance (COP) and CO 2 emission) are achieved for one month of running utilizing two control strategies of the GCHP—standard and optimized regulation of the water pump speed—and a benchmarking of these parameters is achieved. Exploratory research has indicated higher e ﬃ ciency of the system for the ﬂow regulation solution utilizing a BT and programmed control equipment for the circulation pump speed compared with the standard regulation solution (COP sys with 7–8% higher and CO 2 emission level 7.5–8% lower). The second objective is to elaborate a simulation model of the necessary heat / cold in heating and air-conditioning periods, utilizing the Transient Systems Simulation (TRNSYS) program. Finally, the simulation, acquired utilizing the TRNSYS program, is analyzed and compared with experimental information, leading to a good agreement and, along these lines, the simulation model is approved.


Introduction
The European Green Deal is an ambitious growth strategy aiming to transform Europe, setting out to achieve net zero emissions by 2050 [1,2]. A priority area for change is the building sector that accounts for up to 40% of energy use [3]. A large amount of energy used for space heating and cooling in buildings comes from fossil fuels [4], which, with their high carbon dioxide (CO 2 ) emissions, considerably contribute to greenhouse emissions. To limit the impact of fossil fuels, the renewable energy technologies are increasingly used. Among the various technologies to heat and cool buildings, the heat pumps (HPs) are certainly one of the most efficient [5]. The HP moves heat from a thermal source to buildings and rejects the thermal load from the building to the heat sink in the heating and cooling modes, respectively. The use of ground-coupled heat pump (GCHP) systems for heating, ventilating and air conditioning (HVAC) in buildings has been pushed in the last few years as a promising good solution for reducing CO 2 emissions and operating costs [6]. Various GCHPs have been utilized in residential and commercial structures globally in light of their recognizable high 2 of 19 proficiency and environmental friendliness. Several analyses of GCHP systems are reported in the literature by a few researchers [6][7][8][9][10][11][12][13][14].
For GCHP, the coupling between the HP unit and the ground is obtained by means of heat exchangers [15] where the heat-carrier fluid flows, thus connecting the ground with the evaporator or the condenser of the HP unit. Ground heat exchangers (GHEs) can be vertically or horizontally oriented, but vertical configurations (borehole heat exchangers, BHEs) are more frequently used in order to limit the use of land area [16], especially in commercial buildings. When outdoor air is used as a thermal source/sink, the energy efficiency of the HP is affected by the variation in its temperature throughout the year. Ground temperature is more stable during the year and also closer to the indoor air temperature of the building than the outdoor air temperature; consequently the energy efficiency of the GCHP is higher than that of air-source HPs (ASHPs) [17].
Esen et al. [36] compared a GSHP with an ASHP in heating mode from techno-economic point of view, concluding that GSHP was economically convenient due to lower operating costs involved by the higher energy efficiency. Additionally, Lu et al. [37] found similar results in the Melbourne climate. Man et al. [38] made an in situ running efficiency trial of a GCHP for heating and air-conditioning solutions in a temperate zone. The trial demonstrated that the GCHP efficiency was influenced by its irregular or persistent operating modes.
Petit and Meyer [24] analyzed the thermal performances of a GCHP and an ASHP system, concluding that the GCHP, with horizontal or vertical GHEs, was more favorable regarding financial feasibility. Esen and Inalli [39] utilized an in situ thermal response test for a ground-source HP system in Turkey, finding that the thermal resistance of the ground fluctuated marginally with depth. Del Col et al. [40] effectuated a test of the ideal operation mode of a ground-source HP system with variable-speed compressor, water pumps and fans in the loops. They presented the modeling of this HP system type and the model obtained was compared with trial information taken during a heating season. Zeng et al. [41] additionally performed an exploratory performance analysis of GCHP for heating and air-conditioning operating modes in the karsts area. Following the results of the study, they concluded that the GCHP efficiency was affected by its operating modes.
This article focuses on the energy and environmental analysis of an existing vertical GCHP system mounted in an experimental laboratory located in Timisoara, Romania, in a continental temperate climate. The system consisted of a simple GCHP with a single U-tube borehole heat exchanger that had been operational since 2012. Exploratory measurements was utilized to determine the GCHP performances at various operating modes. Generally, the circulation pump of GCHP system operates at fixed speed and supplies a constant flow rate. That is why a variable flow rate control method of the circulation pump is proposed to improve the efficiency of the system and save energy during the heating/cooling (air-conditioning) operation. Thus, the development and implementation of a system energy optimization equipment utilizing quantitative regulation with a buffer tank (BT) and a variable-speed circulation pump is one of the principal creative contributions of this survey. A comparative analysis of primary efficiency parameters (e.g., coefficient of performance (COP), CO 2 emission) achieved for one month of working, utilizing the two control strategies-the standard and optimized regulation of the water pump speed in various operating modes is accomplished. In addition, a model of numerical simulation of heat/cold consumption in heating/air conditioning Energies 2020, 13, 5016 3 of 19 operating modes utilizing Transient Systems Simulation (TRNSYS) is built. At last, the simulation acquired utilizing TRNSYS program is investigated and compared with test estimations.

Testing Laboratory Description
Exploratory analysis of GCHP efficiency was led in a laboratory (Figure 1) situated on the ground floor of the Civil Engineering Faculty building of Timisoara, Romania. The City of Timisoara is located in the West of Romania, and its latitude and longitude are 45 • 47 N and 21 • 17 E, respectively. The climate is continental temperate with four unique seasons. The heating season in Timisoara runs from 1 October to 30 April, and the air-conditioning season runs from 1 May to 30 September.
is located in the West of Romania, and its latitude and longitude are 45°47′ N and 21°17′ E, respectively. The climate is continental temperate with four unique seasons. The heating season in Timisoara runs from 1 October to 30 April, and the air-conditioning season runs from 1 May to 30 September.
The lab room has an area of 47 m 2 , and its height is 3.70 m. The exterior walls are made of a 200 mm permeable block with a 100-mm thermal protecting layer and a 20 mm lime mortar. The U-values are as per the following: walls, 0.345 W/m 2 K, and twofold windows, 2.22 W/m 2 K. The window area is 16 m 2 , and the interior door area is 2.1 m 2 . The indoor air set temperature is 20 °C for the cold season (winter) and 26 °C for the hot season (summer). The outside air set temperature is −15 °C for the cold season and 32.6 °C for the hot season.
The GCHP was mounted in this exploratory lab room heated and cooled through fan coil units. With the input information referred to above, a heating demand of 3.11 kW and an air-conditioning demand of 2.15 kW were determined [42]. The lab room area was assimilated with a three-man apartment area in Timisoara. Taking into account the average daily consumption of domestic hot water (DHW) of 50 l/person, a tank hot-water temperature of 45 °C and a cold water temperature of 20 °C, a DHW load of 4.36 kW was calculated. Figure 2 delineates the month to month energy requirement for heating (positive values) and air-conditioning (negative values) of the research space.  The lab room has an area of 47 m 2 , and its height is 3.70 m. The exterior walls are made of a 200 mm permeable block with a 100-mm thermal protecting layer and a 20 mm lime mortar. The U-values are as per the following: walls, 0.345 W/m 2 K, and twofold windows, 2.22 W/m 2 K. The window area is 16 m 2 , and the interior door area is 2.1 m 2 . The indoor air set temperature is 20 • C for the cold season (winter) and 26 • C for the hot season (summer). The outside air set temperature is −15 • C for the cold season and 32.6 • C for the hot season.
The GCHP was mounted in this exploratory lab room heated and cooled through fan coil units. With the input information referred to above, a heating demand of 3.11 kW and an air-conditioning demand of 2.15 kW were determined [42]. The lab room area was assimilated with a three-man apartment area in Timisoara. Taking into account the average daily consumption of domestic hot water (DHW) of 50 l/person, a tank hot-water temperature of 45 • C and a cold water temperature of 20 • C, a DHW load of 4.36 kW was calculated. Figure 2 delineates the month to month energy requirement for heating (positive values) and air-conditioning (negative values) of the research space.

Description of Experimental Set-Up
The GCHP exploratory system comprised of a BHE, HP unit, BT, water circulation pumps, fan coil units, sink, information acquisition tools, and auxiliary elements, as illustrated in Figure 3.
The heat carrier fluid can be conveyed towards two fan coils units in two discharge regulation types [42]:  The heat carrier fluid can be conveyed towards two fan coils units in two discharge regulation types [42]:

1.
By a distribution pump connected within the HP unit of the GCHP (standard case); 2.
By a constant-speed flow pump connected with a BT. The GCHP computerization can adjust the circulation pump connected to the BT by ON/OFF exchanging. This equipment improves the functioning of whole framework. The BT permits diminishing the GCHP ON/OFF exchanging due to its thermal inertia, and in this way, the energy performance increments. The solution for heat carrier fluid discharge regulation was advanced, utilizing a programmed control equipment of water pump speed [43]. The principal elements of the programmed equipment for pump speed adjustment are presented in Figure 4 [44].  • For heating, a range of outdoor air temperature was considered between −20 °C and 22 °C, taking into account that the indoor air temperature is constant (ti = 22 °C), and the outdoor air temperature varies depending on the day-night and seasonal alternation; • For cooling, the range of the outdoor air temperature was considered to be between 26 °C and 42 °C, taking into account that the indoor air temperature is steady (ti = 26 °C), and the outdoor air temperature fluctuates, relying upon the day-night and seasonal alternation. For the proposed setting, various estimations of the heat demand were determined, utilizing the Romanian standard SR 1907 and of the cold demand, utilizing the Romanian standard SR 6648, contingent upon the temperature difference ∆t among interior and outside, as follows: • For heating, a range of outdoor air temperature was considered between −20 • C and 22 • C, taking into account that the indoor air temperature is constant (t i = 22 • C), and the outdoor air temperature varies depending on the day-night and seasonal alternation; • For cooling, the range of the outdoor air temperature was considered to be between 26 • C and 42 • C, taking into account that the indoor air temperature is steady (t i = 26 • C), and the outdoor air temperature fluctuates, relying upon the day-night and seasonal alternation.
Each estimation of the heat and cold demand corresponds to a flow rate of heat carrier fluid that must be pumped in the system by the circulation pump at a certain speed ν of its engine and, in this way, the speed control of the circulation pump joined to the BT is controlled.
Following the calculations performed, the variation curve-by points of the frequency ν, in Hz-of the frequency converter was obtained depending on the temperature difference ∆t in • C. Applying the geometric interpolation by the numerical method of the smallest squares [45], the analytical expression of this curve was determined, under the forms:

•
For heating: Figures 5 and 6 graphically illustrate both the curve of the function given by points for heating and cooling, respectively, and the corresponding geometric interpolation curve analytically expressed by Equation (1) or Equation (2), which shows a good concordance between them.   . Figure 6. Graph of the function ν(Δt) for cooling operation.   The temperature difference between within and outside of the heated/air-conditioned room is estimated by temperature sensors TS1 and TS2, coupled with a programmable logic controller (PLC). In the interior memory of the PLC, a calculation algorithm, dependent on Equations (1) and (2), is actualized for ascertaining the frequency of the frequency converter, which relies upon the measured temperature difference. The PLC sends to the frequency converter the comparing frquency incentive to guarantee the fluid discharge, as indicated by the heating/air-conditioning requirement at that time. Additionally, the PLC permits the ON/OFF exchanging control of the The temperature difference between within and outside of the heated/air-conditioned room is estimated by temperature sensors TS1 and TS2, coupled with a programmable logic controller (PLC). In the interior memory of the PLC, a calculation algorithm, dependent on Equations (1) and (2), is actualized for ascertaining the frequency of the frequency converter, which relies upon the measured temperature difference. The PLC sends to the frequency converter the comparing frquency incentive to guarantee the fluid discharge, as indicated by the heating/air-conditioning requirement at that time. Additionally, the PLC permits the ON/OFF exchanging control of the circulation pump. For the concurrent ON/OFF exchanging of the circulation pump and GCHP, common temperature values of this procedure were determined.

Borehole Heat Exchanger
In the GCHP application, the BHE ensure that the heat extraction or rejection from the ground. The running the BHE instigates concurrent heat and moisture flows in the soil in which the loop is installed. The heat exchange among BHE and soil is fundamentally practiced by thermal conduction and, in a specific way, by moisture movement. In this manner, soil types and temperature, similar to moisture gradients, have critical impacts on the heat exchange process [41]. Thus, the properties of the soil at the trial site are fundamental.
The straightforward vertical borehole of this exploratory GCHP has a depth of 80 m. Antifreeze solution (glycol 30%) flows in a sole polyethylene U-tube of 32 mm interior diameter, with a 60-mm partition between the return and supply tubes, placed in the borehole. The borehole with a total diameter of 110 mm was loaded up with sand and finished with a bentonite layer at the top to stay away from the interruption of contaminations in the aquifers. The mean temperature over the full borehole depth was 15.2 • C.
The ground proprieties are based on estimations acquired from the Banat Water Resources Management Agency [46]. The mean thermal conductivity and thermal diffusivity of the ground over the full borehole depth were 1.91 W/(mK) and 0.80 × 10 −6 m 2 /s, respectively. The boreholes were backfilled with grout, blended in with penetrating mud, cement, and sand in explicit extents. The thermal conductivity and diffusivity of the grout tried by the maker were 2.31 W/(mK) and 0.92 × 10 −6 m 2 /s, respectively.

Heat Pump Unit Equipment
The HP unit utilized in this study is a ground-to-water HP with a scroll air-tight compressor that works with the refrigerant R410A. The nominal heating and air-conditioning limits are 6.5 kW and 3.8 kW, respectively. The HP unit is an electrically-powered HP, having an interior refrigerating system and a DHW tank with a 175-L capacity. Running the HP was administered by an electronic regulator, which, contingent upon the water return temperature, turned the HP compressor on or off. The circulation pump of heat source was adjusted by the HP regulator that starts the heat source pump 30 s before compressor actuation.
The energy efficiency of a reversible GCHP is expressed by COP hp , which can be calculated utilizing the following formulas: • For heating operating mode: in which E t is the useful heat and E el is the electricity consumption of HP.

•
For cooling operating mode: in which the energy efficiency ratio EER hp , in Btu/(Wh), describes the cooling performance, and the value 3.412 is the transformation factor from Watt to Btu/h (British Thermal Units per hour). Another way of expressing COP hp is dependent on the temperature of the hot environment (T h ) and that of the cold environment (T c ) [10]: The COP of the overall system (COP sys ) is calculated with Equation (3), in which E el is the electricity consumption of the GCHP, including the electricity consumption of the HP compressor, circulation pumps, fan coil units, frequency converter, and PLC.
The carbon dioxide (CO 2 ) emissions C CO 2 during the heating system operation are evaluated as follows: where g el = 0.547 kg CO 2 /kWh is the specific CO 2 emission factor for electrical energy [47]. The thermal power E t for heating/cooling operating modes and electrical consumption E el , utilized by the HP unit or the GCHP system, must be measured to obtain the COP hp or COP sys and CO 2 emissions.

Water Circulation Pumps
The water circulation circuits of the GCHP comprised of a GCHP-BT circuit and BT-fan coil unit circuit. Two centrifugal pumps with the rated flow rate of 2.8 and, respectively, 5.5 m 3 /h were used for the two water circuits. The primary circulation pump (constant-speed pump coupled to HP unit) was adjusted by the GCHP computerization, and the subsequent pump (variable-speed pump coupled to BT) was adjusted by a programmed control equipment.

Fan Coil Units
As terminal units of the GCHP, two parallel coupled fan coil units, whose global thermal power was 3.2 kW, were utilized.

Information Acquisition System
To more readily comprehend the behavior of the GCHP, an information acquisition system was created to screen the most relevant parameters of the system. The GCHP information procurement system comprises the inside and outside temperatures, supply/return temperatures, heat source temperature (outlet BHE temperature), DHW temperature, relative air humidity, and primary working parameters of the system elements.

Measurement Instrumentations
A network of sensors was created to permit the registration of the most important parameters of the system, which are described in detail in [42], together with the other measuring instruments, including their accuracy.

Methodology
Various experiments were performed on the GCHP system under steady-state operation regimes to assess its COP and CO 2 emission [42]. The tests were partitioned into heating and cooling tests with or without the production of DHW. The system was examined for 2 years. The exploratory estimations were accomplished for two strategies of the flow rate regulation of the heat carrier fluid (water) in the system: solution (1)-standard regulation, by constant-speed circulation pump coupled with HP unit, and solution (2)

Heating Operation
Variation in the indoor air temperature (t i ) and outside temperature (t e ), registered for a 1-month time frame, is illustrated in Figure 8. It ought to be noticed that, in the solution (2), a diminished indoor temperature was acquired, around the design temperature of 22 • C, which leads to all the more suitable comfort in the heated space. Furthermore, the diminished wavering of the antifreeze solution temperature prompts a lower heat source requirement. The heat source temperature in the cold season is up to 12-13 • C higher than the outside temperature, which raises the capacity and improves the performance of a GCHP. Table 1 presents the summarization of the average estimation of the temperatures (t i , t e , t s ), power consumption (E el ), helpful heat for heating (E t ), COP of the system (COP sys ) and the HP unit (COP hp ), and CO 2 emission (C CO 2 ).

Heating Operation
Variation in the indoor air temperature (ti) and outside temperature (te), registered for a 1-month time frame, is illustrated in Figure 8. It ought to be noticed that, in the solution (2), a diminished indoor temperature was acquired, around the design temperature of 22 °C, which leads to all the more suitable comfort in the heated space. Furthermore, the diminished wavering of the antifreeze solution temperature prompts a lower heat source requirement. The heat source temperature in the cold season is up to 12-13 °C higher than the outside temperature, which raises the capacity and improves the performance of a GCHP. Table 1 presents the summarization of the average estimation of the temperatures (ti, te, ts), power consumption (Eel), helpful heat for heating (Et), COP of the system (COPsys) and the HP unit (COPhp), and CO2 emission (   The COPhp estimations of the HP unit for standard and improved solution are 4.81 and 5.07, respectively. For the optimized solution, the COPsys = 4.40 is 7.5% higher and the CO2 emissions level is 7% lower than on the classic solution. Because of the characteristics of the climate, ground, and the site where the estimations were performed and a higher underground water flow rate, the heat source temperature is increased, and the COPhp and COPsys are remarkable values for the two cases. Consequently, the comparison of these results with others obtained in similar analyses [6,48], but in different local conditions, showed that the efficiency obtained here are significantly improved.   The COP hp estimations of the HP unit for standard and improved solution are 4.81 and 5.07, respectively. For the optimized solution, the COP sys = 4.40 is 7.5% higher and the CO 2 emissions level is 7% lower than on the classic solution. Because of the characteristics of the climate, ground, and the site where the estimations were performed and a higher underground water flow rate, the heat source temperature is increased, and the COP hp and COP sys are remarkable values for the two cases. Consequently, the comparison of these results with others obtained in similar analyses [6,48], but in different local conditions, showed that the efficiency obtained here are significantly improved. Figures 9 and 10 show the variation in the time of the characteristic temperatures (t e , t i , and t s ) recorded during the experiments in the cooling operating mode. It should be noted that, in case (2), a more decreased indoor air temperature was acquired, around the design temperature of 26 • C, prompting higher comfort in the heated space. Table 2 summarizes the average estimations of the temperatures (t i , t e , t s ), power consumption (E el ), useful heat for air-conditioning (E t ), system EER (EER sys ) and CO 2 emissions (C CO 2 ). The energy efficiency ratio EER sys of the GCHP system was 8% higher, and the CO 2 emissions was 8% lower in the optimized solution than in the classic solution.

Cooling Operation
Energies 2020, 13 (2), a more decreased indoor air temperature was acquired, around the design temperature of 26 °C, prompting higher comfort in the heated space. Table 2 summarizes the average estimations of the temperatures (ti, te, ts), power consumption (Eel), useful heat for air-conditioning (Et), system EER (EERsys) and CO2 emissions ( 2 CO C ). The energy efficiency ratio EERsys of the GCHP system was 8% higher, and the CO2 emissions was 8% lower in the optimized solution than in the classic solution.   The correlation of the exploratory GCHP efficiency in the heating and air-conditioning period (Tables 1 and 2) demonstrate that the efficiency in heating and air-conditioning mode was practically equivalent. This is because the heating burden was higher than the air-conditioning demand, and additionally, the power consumption in the heating mode was higher than the power consumption  The correlation of the exploratory GCHP efficiency in the heating and air-conditioning period (Tables 1 and 2) demonstrate that the efficiency in heating and air-conditioning mode was practically equivalent. This is because the heating burden was higher than the air-conditioning demand, and additionally, the power consumption in the heating mode was higher than the power consumption in air-conditioning operation.
The COP estimations of the system were compared with the current COP values recorded in GCHPs investigation. The exploration of Man et al. [38] uncovered COPs of GCHPs of 4.18-4.57 in the winter and 3.9-4.54 in the summer. Likewise, the investigation in the hot season of Michopoulos et al. [31], which utilized a vertical GHE at a depth of 80 m, obtained a COP sys of 4.4-4.5. It is seen that the efficiency esteems achieved here are comparative.

Heating and DHW Provision
For the comparative analysis of the two flow control strategies, an equivalent DHW volume was utilized, V dhw = 1.21 m 3 . The variation in the indoor air (t i ), outside air (t e ), DHW (t dhw ) and heat source (t s ) temperatures-registered for seven days' time span for every one of the two investigated solutions, from 15 January 2013 to 21 January 2013 and 10 January 2014 to 16 January 2014-is illustrated in Figures 11 and 12. The heat source temperature in the cold season is roughly 15 • C higher than the outside temperature.   Table 3 summarizes the average estimations of the characteristic temperatures (ti, te, tdhw, ts), DHW volume (Vdhw), electrical energy consumption (Eel), helpful heat (Et), COP of the system (COPsys) and of the HP unit (COPhp), and CO2 emission ( 2 CO C ).
Although the COPsys recorded equivalent values for the two solutions, the results of the experiments show that, when utilizing programmed control equipment for the circulation pump speed, an electrical energy saving of 3% and a CO2 emissions reduction of 3% were obtained.
Examining the test information (Tables 1 and 3) results demonstrate that the COPsys of the system working in heating and DHW provision, comparative with the heating operating mode, diminished fundamentally in the domain 20.6-23.9%, in contrast with the two solutions, from 4.06-4.40 to 3.22-3.36, respectively. The COPhp estimations of the HP unit for solutions (1) and (2) are  Table 3 summarizes the average estimations of the characteristic temperatures (t i , t e , t dhw , t s ), DHW volume (V dhw ), electrical energy consumption (E el ), helpful heat (E t ), COP of the system (COP sys ) and of the HP unit (COP hp ), and CO 2 emission (C CO 2 ). Table 3. GCHP efficiency for heating and DHW operation. Although the COP sys recorded equivalent values for the two solutions, the results of the experiments show that, when utilizing programmed control equipment for the circulation pump speed, an electrical energy saving of 3% and a CO 2 emissions reduction of 3% were obtained.

Flow Regulation
Examining the test information (Tables 1 and 3) results demonstrate that the COP sys of the system working in heating and DHW provision, comparative with the heating operating mode, diminished fundamentally in the domain 20.6-23.9%, in contrast with the two solutions, from 4.06-4.40 to 3.22-3.36, respectively. The COP hp estimations of the HP unit for solutions (1) and (2) are 3.80 and 3.96, respectively.

Cooling and DHW Provision
To establish the GCHP efficiency in summer, test estimations over 1 week for every one of the two analyzed control strategies was accomplished from 26 June 2013 to 2 July 2013 and 29 June 2014 to 5 July 2014. During the experiments, both the cooling demand and the DHW demand for a family that utilizes a DHW volume V dhw = 1.35 m 3 were provided. Figure 13 outlines the evolution of indoor air (t i ) and outside air (t e ) temperatures, and Figure 14 shows the variation in the DHW (t dhw ) and heat source (t s ) temperatures. The test information (Tables 2 and 4) on the COPsys of the GCHP in air-conditioning and DHW operation, comparative with air-conditioning operating, show that there is an abatement in the domain of 3.2-8.6% compared with the two solutions, from 4.03-4.43 to 3.90-4.05, respectively.   Table 4 includes a summarization of the average estimations of the characteristic temperatures (t i , t e , t dhw , t s ), DHW volume (V dhw ), power consumption (E el ), helpful heat (E t ), COP of the system (COP sys ) and CO 2 emissions (C CO 2 ). Table 4. GCHP efficiency for air-conditioning and DHW operation.     Though the COP sys recorded values about equal for the two solutions, the trial results show that, when utilizing the programmed control equipment for the circulation pump speed, an electrical energy saving of 5% and a CO 2 emissions reduction of 5% were obtained.

Flow Regulation
The test information (Tables 2 and 4

Uncertainty Analysis
Uncertainty analysis (the analysis of vulnerabilities in test estimation and results) is important to assess the exploratory information [49,50]. The methodology utilized to carry out uncertainty analysis was described by Holman [49]. An outcome Z is a function given by the independent variables x1, x2, x3…xn. On the off chance that the errors in the independent variables w1, w2, w3…wn are all given with same chance, then uncertainty in the outcome wZ with this chance is determined utilizing the following formula [49]: In the current survey, the temperatures, heat, and electricity were estimated with suitable instruments, already clarified. Error analysis for assessing the maximum uncertainty in the trial results was conducted utilizing Equation (7). It was discovered that the most extreme uncertainty in

Uncertainty Analysis
Uncertainty analysis (the analysis of vulnerabilities in test estimation and results) is important to assess the exploratory information [49,50]. The methodology utilized to carry out uncertainty analysis was described by Holman [49]. An outcome Z is a function given by the independent variables x 1 , x 2 , x 3 . . . x n . On the off chance that the errors in the independent variables w 1 , w 2 , w 3 . . . w n are all given with same chance, then uncertainty in the outcome w Z with this chance is determined utilizing the following formula [49]: In the current survey, the temperatures, heat, and electricity were estimated with suitable instruments, already clarified. Error analysis for assessing the maximum uncertainty in the trial results was conducted utilizing Equation (7). It was discovered that the most extreme uncertainty in the outcomes is in the COP sys , with an acceptable uncertainty interval of 1.31-1.69% in the heating operating mode and of 2.29-3.38% in the air-conditioning operating mode. The accuracy of the COP hp was evaluated in the range of 1.82-2.33% in the heating mode.

Numerical Simulation of Useful Thermal Energy Using TRNSYS Software
The essential rule of the TRNSYS software [51] is the implementation of algebraic and first-order conventional differential equations depicting physical components into programming subroutines (types) with a standard interface. The STEC library is based on steady-state energy preservation formulated in thermodynamic amounts (temperature, pressure, and enthalpy). One of the principal advantages of TRNSYS program for the modeling and design of HP systems is that it incorporates parts for the computation of building thermal loads, specific components for HVAC systems, HP systems, circulation pumps, modules for vertical GHEs and storing heat, similar to climatic information files, which make it an entirely appropriate instrument to model a GCHP system to provide heating and air-conditioning to a building. The METEONORM program [52] has a database of over 7000 worldwide stations that can generate data on monthly, daily, or hourly time scales.
For the validation of the numerical simulation model, four statistical assessment criteria (indices) were utilized-the coefficient of multiple determinations (R 2 ), the root mean square error (RMSE), the coefficient of variation (c v ), and relative error (e r )-as defined below [53,54].
The coefficient R 2 presents the general concordance among measured (real) and simulated (predicted) time series, and ranges from 0, for a poor model, to 1 for a perfect model, and is expressed as: y sim,j − y mea,j 2 n j=1 y 2 mea,j (8) in which n is the number of data points; y mea,j is the measured value of data point j; y sim,j is the predicted value.
The RMSE is a measure of overall performance over the entire range of the data set expressed by the formula: Greater or equal to 0, the RMSE indicates an ideal model fit when it equals 0. The c v (in %) can be interpreted as an order of magnitude of the repeatability relative standard uncertainty and it is defined by: where y mea,j is the average value of all measured data points. The relative error e r , in %, can be calculate as follows: e r = y sim, j − y mea, j y mea,j 100 (11)

Definition of the Operational Plan
To simulate the heat energy utilized for the heating/air-conditioning of the exploratory laboratory space, the operational connections were set up between the construction and all inner and outer factors. Figure 17 shows the operation plan made in TRNSYS, where the construction thermal behavior was modeled utilizing a "Type 56" subroutine elaborated with the TRNBuild interface. Climate information for the Timisoara was obtained from the METEONORM information base [52], and the climate information readers, "Type 109" and "Type 89d", were utilized to change over the information in a form intelligible from TRNSYS. The simulation model considered the outside air penetration, heat/cold source type, and internal gains. Likewise, light thresholds, air-conditioning, and shading were characterized for a decent way to deal with the model. The results were extracted with an online plotter ("Type 65").
where j y , mea is the average value of all measured data points.
The relative error er, in %, can be calculate as follows:

Definition of the Operational Plan
To simulate the heat energy utilized for the heating/air-conditioning of the exploratory laboratory space, the operational connections were set up between the construction and all inner and outer factors. Figure 17 shows the operation plan made in TRNSYS, where the construction thermal behavior was modeled utilizing a "Type 56" subroutine elaborated with the TRNBuild interface. Climate information for the Timisoara was obtained from the METEONORM information base [52], and the climate information readers, "Type 109" and "Type 89d", were utilized to change over the information in a form intelligible from TRNSYS. The simulation model considered the outside air penetration, heat/cold source type, and internal gains. Likewise, light thresholds, air-conditioning, and shading were characterized for a decent way to deal with the model. The results were extracted with an online plotter ("Type 65").  Figure 17. Simulation scheme of the thermal energy consumption in TRNSYS software.

Comparison of Simulated Results with Exploratory Information
The simulated values of heat energy utilized for space heating and conditioning were obtained from the simulations performed for a period of 1 year (8760 h) and included the values measured in Table 5. The statistical criteria (RMSE, c v , and R 2 ) for the assessment of simulation model accuracy are presented in Table 6, considering the operation of GCHP system in various modes, and the relative errors e r , which are included in Table 5.
The maximum relative error e r of the simulated values in TRNSYS, compared to those measured, was about 1.60% during heating and around 1.65% during air-conditioning, which is very satisfactory. The values of the RMSE and c v performance indices in the heating period are approximately 2.7222 and 0.0141, respectively, and in the air-conditioning period they are approximately 3.080 and 0.0238, respectively. The R 2 -values in the two operating modes are approximately 0.9998, which is very acceptable. Accordingly, the simulation model was validated by the experimental measurements.

Conclusions
The utilization of HPs in present-day buildings with improved thermal protection and the diminished heat demands is a good option, in contrast to the old-style heating/air-conditioning and DHW heating solutions. This article focused on the assessment of the GCHP system efficiency, giving heating/air-conditioning and DHW a laboratory trial. The performance of the GCHP system, in both optimized and standard mode, was assessed via both experiments and computer simulations with TRNSYS. The principal conclusions of this experimental survey are as follows: 1.
The accomplished trial investigations showed higher efficiency of the GCHP for the flow rate regulation solution utilizing a BT and programmed control device for circulation pump speed, comparative with an old-style regulation solutions (COP sys with 7-8% increase, and CO 2 emissions with 7.5-8% decrease).

2.
The overall COP of the GCHP in the two solutions was COP sys > 4 for the heating or air-conditioning period, and 3 < COP sys < 4 for heating or air-conditioning and DHW provision.

3.
In the standard and optimized flow regulation solutions, the experimental evaluations of the COP for HP unit (COP hp ) during the heating and DHW period were 3.80 and 3.96, respectively, and during the heating period they were 4.81 and 5.07, respectively. 4.
The GCHP system utilized for cooling and DHW provision has higher performance (COP sys = 3.90 for the classic case and 4.05 for optimized case) than that for heating and DHW production (COP sys = 3.22 for the classic case and 3.36 for optimized case) because the heat extracted from the soil for heating and DHW production is higher than the heat injected into the soil for air-conditioning and DHW heating.

5.
When utilizing the variable-speed water pump, electrical energy saving and a decrease in the CO 2 emissions of 3% for air heating in the laboratory and 5% for air cooling in the space were obtained simultaneously with DHW generation.

6.
The proposed numerical simulation model can be utilized as an instrument to establish the GCHP performance in various operating modes in order to optimize the energy efficiency of the system and guarantee the occupant's comfort consistently. 7.
A limitation of this study is that the optimized solution for regulating the flow rate of heat-carrying fluid using a buffer tank and an automatic control device of the circulation pump speed is applicable only to residential buildings. 8.
In order to possibly improve the energy efficiency of the GCHP system, in-depth research is needed in the future, aimed at integrating solar photovoltaic collectors into the system, which will produce electricity to drive the circulation pump in the water pumping process.