Search

With increasing global concerns regarding indoor air quality (IAQ) and air pollution, concerns about regularly replacing ventilation devices, particularly high-efficiency particulate air (HEPA) filters, have increased. However, users cannot easily determine when to replace filters. This paper proposes models to estimate the dust loading levels of HEPA filters for an energy-recovery ventilation system that performs air purification. The models utilize filter pressure drops, the revolutions per minute (RPM) of supply fans, and rated airflow modes as variables for regression equations. The obtained results demonstrated that the filter dust loading level could be estimated once the filter pressure drops and RPM, and voltage for the rated airflow were input in the models, with a root mean square error of 5.1–12.9%. Despite current methods using fewer experimental datasets than the proposed models, our findings indicate that these models could be efficiently used in the development of filter replacement alarms to help users decide when to replace their filters.

Semi-transparent thin-film solar cells with ultra-thin Cu(In,Ga)Se2 (CIGS) absorber layers were fabricated on fluorine-doped tin oxide coated glass substrates using a single-stage co-evaporation process. The effects of the deposition conditions, such as the absorber thickness, deposition temperature, and post-deposition treatment using NaF on the material properties of the CIGS films were evaluated. In addition, the power conversion efficiency (PCE) and average visible transmission (AVT) in the wavelength range of 420–720 nm for semi-transparent CIGS solar cells using these absorber layers were systematically investigated. The AVT of the solar cells was directly dependent on the thickness of the CIGS absorber films. When the film thickness of the absorber layers increased from 300 nm to 950 nm, the AVT decreased gradually from 9.04% to 0.30%, while their PCE values did not show a clear dependence on the absorber thickness. Narrow columnar grains were observed in the microstructure of the ultra-thin CIGS absorber layers, irrespective of the film thickness and deposition temperature. Doping the ultra-thin CIGS absorber layers with Na after film deposition enhanced solar cell performance (the open circuit voltage (VOC) and fill factor (FF) were improved). The enhanced cell performance was attributed to the lack of a resistive GaOx film at the CIGS/FTO interface and thermal damage to the FTO back contact. The PCE values of semi-transparent solar cells with ultra-thin CIGS absorber layers prepared under the optimal deposition conditions varied from 6.46% to 9.78% when the CIGS absorber thickness was varied from 200 nm to 300 nm, respectively, while the corresponding AVT values changed from 18.59% to 9.04%.

Lead-free double perovskite halides are emerging optoelectronic materials that are alternatives to lead-based perovskite halides. Recently, single-crystalline double perovskite halides were synthesized, and their intriguing functional properties were demonstrated. Despite such pioneering works, lead-free double perovskite halides with better crystallinity are still in demand for applications to novel optoelectronic devices. Here, we realized highly crystalline Cs2AgBiBr6 single crystals with a well-defined atomic ordering on the microscopic scale. We avoided the formation of Ag vacancies and the subsequent secondary Cs3Bi2Br9 by manipulating the initial chemical environments in hydrothermal synthesis. The suppression of Ag vacancies allows us to reduce the trap density in the as-grown crystals and to enhance the carrier mobility further. Our design strategy is applicable for fabricating other lead-free halide materials with high crystallinity.

The thermal environment experienced by human is significantly affected by the radiant heat. Information regarding the effect of radiant temperature on the occupants in a building is necessary for a comprehensive evaluation of indoor thermal comfort and for a corresponding control of the building's passive and active systems. This paper presents a real-time estimation method of MRT distributions in indoor spaces that does not require a direct installation of measurement devices at occupant locations. To consider the effect of longwave radiation, a pan-tilt IR camera was applied to remotely measure interior surface temperatures. For shortwave radiation, a single pyranometer installed outside was applied to estimate diffuse and direct irradiance at the occupant locations. The real-time MRT monitoring process and the Pan-tilt IR scanning system⁺ are derived. An experiment was performed to validate the applicability of the derived MRT monitoring process. The MRT at three target locations were measured by a globe thermometer and also estimated by the derived method. The measured and estimated MRT values were compared, and the results showed high goodness of fits as the R2 values were higher than 0.80 for all target locations. The monitoring results confirm that the MRT value can differ up to 15 °C according to occupant locations because of solar radiation. This paper discusses the applicability in various fields and possible future works of the proposed MRT estimation process.

Solar radiation is a key element that has a large impact on the indoor thermal environment. It is therefore important to incorporate the impact of indoor solar irradiance when evaluating the thermal environment within a real building. This paper proposes a real-time estimation method for solar irradiance reaching multiple occupant locations within a building from global irradiance data measured at a single point outside. First, Watanabe's decomposition model and Perez's transposition model are selected and applied to calculate direct and diffuse irradiance incident on a tilted facade. Then, the sky view factor (SVF), horizontal view factor (HVF), and unshaded fraction (Fu) are applied to calculate reductions in diffuse solar irradiance incident on the window area. The sunlit factor (Fsunlit) is applied to determine whether the direct solar beam incident on the occupant location is blocked or not. To calculate these four factors, a method of generating a simple sky model and a building model is presented for application to existing buildings. An experiment is performed to verify the accuracy of the calculation method. The indoor irradiance at a target location was measured by a pyranometer and also estimated by the calculation method. The measured and estimated indoor irradiance values were compared, and the results show a high estimation performance of the calculation method with a coefficient of determination of 0.95. This paper discusses the applicability of the proposed calculation method to the real-time evaluation of an indoor environment and to the control of building passive systems.

Airtightness measurement standards using the fan pressurization method are widely used for measuring residential units such as detached houses and apartment units. The measurement of airtightness for small buildings can be readily conducted through existing fan pressurization methods. However, it can be difficult to achieve accurate measurement results for large buildings, as the height and volume of the buildings are considerably increased. This paper empirically investigates factors that are to be considered in regard to airflow rate (Q) and pressure difference (ΔP), which are measurement values required for deriving more accurate airtightness values for large buildings. The specifications of the standards ISO 9972, EN 13829, ASTM E779 , JIS A 2201, ATTMA TS 1, and CAN/CGSB 149.15 relating to Q and △P are compared. Test conditions with respect to Q and △P are discussed for large building airtightness measurements. In measuring airtightness, a method is proposed in which the pressure difference at each floor of the building is measured and an average pressure difference value is used. Experiments were conducted on two multi-story office buildings using the method based on measuring pressure difference according to existing standards in tandem with the proposed pressure difference measurement method. The results of the experimental study show that the vertical pressure distribution across the overall building envelope can be significantly different when a building is pressurized. The airtightness values may also be different according to the measurement methods. This paper demonstrates how important it is to derive accurate △P values in measuring large building airtightness.

External movable shading devices are mainly used to improve the thermal and lighting environment of an interior space by preventing overheating during the cooling season as well as by preventing glare. An external movable shading device can be controlled using various methods, where parameters related to the lighting environment, such as daylighting and glare, are commonly adopted as major variables associated with the control of the shading device. In this paper, an algorithm is provided with which an external movable shading device can be operated in any one of three control modes as selected according to the operation purpose preferred by the occupant. The algorithm can be applied to external movable shading devices of various shapes. The shaded fraction is adopted as a major control element, and the control purpose of the external shading device is configured to include energy conservation, illuminance satisfaction, and glare protection. The algorithm includes an hourly evaluation procedure for obtaining the range of allowable shading device movement that satisfies each of the three operation purposes. Real-scale mockup tests for the three control modes were conducted during different time periods in summer to verify the performance of the control strategies. The test results regarding the energy and environment performance for each control mode show that the control algorithm can be selected to provide satisfactory compromises between energy and visual comfort factors. The result of this study can provide users with options for controlling movable shading devices on a building according to their preferred operation purposes.

In this study, we developed a virtual air flow sensor using an in-situ damper performance curve to secure the stability of control of the variable air volume (VAV) terminal unit, and also established the in-situ measurement procedure. The minimum air flow rate of the VAV terminal unit was related to the energy consumption, it was important to determine the minimum air flow rate suitable for the situation of each room in terms of energy saving. However, it was difficult to set the minimum air flow rate low setpoint due to the low accuracy of the air flow sensor of the VAV terminal unit. This paper suggested a virtual air flow sensing method using an in-situ damper performance curve in the VAV terminal unit. The input factors of the virtual sensor were developed on the basis of the supply fan speed and damper opening ratio, which could be easily obtained from the existing control system. In addition, the in-situ measurement procedure of the virtual air flow sensor was developed by dividing the procedure into five steps. Finally, reliability of the virtual air flow sensor was evaluated through uncertainty analysis of input variables and relative error analysis, in comparison with the conventional air flow rate measurement method. The developed virtual air flow sensor was found to have an uncertainty of up to 8.8%, and it was also found that the closer to the maximum the values of the input variables, the lower the uncertainty. In addition, verification of relative error with respect to the measured values by the hot-wire anemometer was conducted by varying operation conditions to a total of 12 cases, and as a result, relative error was found to be up to 5.6%. In addition, the results of long-term experiment showed that relative error was within about 9.5%, and thus, the feasibility for field application and control was confirmed.

The mean radiant temperature (MRT) is an indicator for evaluating the radiant heat environment near occupants and is determined by the radiant heat exchange between the occupants and their surroundings. To control various heating and cooling systems according to the occupants’ thermal comfort, it is essential to consider MRTs in the real-time evaluation of thermal environment. This study proposes a pan–tilt infrared (IR) scanning method to estimate the MRTs at multiple occupant locations in real buildings. The angle factor was calculated by defining the specific classification criteria for dividing the entire indoor surface into sub-surfaces. The coupling IR camera and pan–tilt motor were applied to enable storing data pairs of IR thermal image frame (IR image frame) and pan–tilt angle so each surface area taken by the IR camera can have its direction information. The measurement method of the mean surface temperature using the pan–tilt IR system was presented. The pan–tilt IR system hardware and MRT monitoring software were developed. An experiment was performed to verify the applicability of the proposed pan–tilt IR scanning method. By comparing the surface temperatures measured using a contact thermometer and the proposed IR system, the contact thermometer could cause inaccurate measurement of surfaces with a non-uniform distribution of temperature. The difference between surface temperatures increased by up to 15 °C and, accordingly, the MRT distributions differed by up to 6 °C within the same space. The proposed IR scanning method showed good applicability in various aspects. This paper reports that the MRT has a significant effect on the occupants’ thermal comfort and also suggests considering MRTs in the real-time evaluation of thermal environment to control various heating and cooling systems appropriately.

A movable shading device installed on an adaptive facade is able to change the positions of the shades according to the various control objectives (i.e. daylighting, thermal comfort, glare, energy, etc.) of the occupant. Since the key variables and the control method may change according to the control objectives, it is important to determine the hourly states of the shades in consideration of all of these factors. This study proposes a multi-purposes control strategy of adaptive facades for heating, cooling, lighting energy and glare. To depict the movements of the movable shading devices, the operable range is represented by an n number of fixed shade states having constant intervals. The control strategy entails deriving allowable shade states for glare prevention and energy conservation independently and then deriving the intersection range between the two shade state ranges. A case study was performed to verify the changes in hourly shade states and the energy performance of the envelope obtained by applying the proposed control strategy to an adaptive facade having a movable folding type shading device. The results of the case study show that hourly shade states were determined to satisfy the objectives of glare prevention and energy conservation simultaneously based on weather conditions and indoor requirement settings. From the case study, it can be seen that the proposed method can block the entry of direct solar radiation by lowering the unshaded fraction according to whether or not there would be discomfort glare while increasing the view factor in certain sections to promote the entry of diffuse solar radiation into the indoor space during the cooling and heating periods. The results of the case study also show that there can be large differences in the heating, cooling, and lighting energy performance of the building envelope for each orientation depending on the order of priority given to the objectives of the shading control, i.e. the objective of energy conservation only and the objective of simultaneous glare prevention and energy conservation. However, it was possible to effectively reduce the differences in energy proportions by adopting translucent shades.

In this study, a model equation is derived that uses a statistical analysis based on empirical models to predict the airtightness of reinforced concrete apartment buildings popular in Asian regions. Airtightness data from 486 units personally measured by the authors in the past eight years are used. As major variables used in the prediction model, two groups of variables are configured for the geometric components of the envelope, which is a major path of airflow in a building and is where air infiltration and leakage occur. The two groups of variables represent (1) the areas of the individual components forming the envelope and (2) the connection lengths between different components of the envelope. For the effective prediction of airtightness, correlation analysis and multiple regression analysis were applied step by step in this study. The results of the correlation analysis indicated that the areas of the slab and the window are the area variables that present the greatest impact, whereas the perimeter length of the window is the connection length variable that presents the greatest impact. Using a multiple linear regression analysis method, airtightness prediction model equations can be derived, and it is found that the model with variables for area is able to predict airtightness more accurately compared to the two models derived from variables for connection length and all variables for area and connection length. Although the statistical approach in this study shows a limitation in that the prediction results may vary depending on the attributes and type of data collected by countries, the methodology and procedure in this study contribute to similar studies for making prediction models and finding the influence of variables in the future with high applicability and feasibility.

Accurate measurement of air flow rate is essential in automatic building control using the variable air volume (VAV) system. In order to solve the problems of the existing air flow measurement method and improve the accuracy of air flow control, this study developed a data-based multiple regression air flow prediction model. The independent variables used in the development of the predictive model were selected as the factors used for control and monitoring when operating with variable air flow rate in the existing air conditioning system. Data collection and correlation between independent variables and air flow rate of the terminal unit were analyzed. Using the IBM SPSS statistics version 25, an air flow rate prediction model was developed using multiple regression analysis. Reliability of model was evaluated by comparing the measured airflow. The relative error of &minus, 9.3% to 10.4% is shown when comparing the estimated air flow rate by the developed model with the measured air flow rate.

Since the thermal environment of large space buildings such as stadiums can vary depending on the location of the stands, it is important to divide them into different zones and evaluate their thermal environment separately. The thermal environment can be evaluated using physical values measured with the sensors, but the occupant density of the stadium stands is high, which limits the locations available to install the sensors. As a method to resolve the limitations of installing the sensors, we propose a method to predict the thermal environment of each zone in a large space. We set six key thermal factors affecting the thermal environment in a large space to be predicted factors (indoor air temperature, mean radiant temperature, and clothing) and the fixed factors (air velocity, metabolic rate, and relative humidity). Using artificial neural network (ANN) models and the outdoor air temperature and the surface temperature of the interior walls around the stands as input data, we developed a method to predict the three thermal factors. Learning and verification datasets were established using STAR CCM+ (2016.10, Siemens PLM software, Plano, TX, USA). An analysis of each model’s prediction results showed that the prediction accuracy increased with the number of learning data points. The thermal environment evaluation process developed in this study can be used to control heating, ventilation, and air conditioning (HVAC) facilities in each zone in a large space building with sufficient learning by ANN models at the building testing or the evaluation stage.

Exterior venetian blinds (EVBs) are widely used in commercial and residential buildings, and show a high performance level for both shading and lighting purposes. However, an EVB installed over an open window also influences natural ventilation rates. In this study, applying the pressure loss rate of an EVB to wind-driven natural ventilation rates was proposed. A mock-up building was set that the EVB was installed on the front opening. Wind generator was installed toward the front opening in order to describe wind-driven cross ventilation in a single zone. To investigate wind pressure changes through the EVB plane, 16 pressure measuring taps were installed on each of the front and back side of the EVB plane. Various wind speeds were induced toward the front opening so as to derive the pressure loss rate for the four EVB slat angle cases (0°, 45°, 90°, and no shading). The pressure loss rates for each case were derived from the field measurements. The results show that the pressure loss rates have a range of 0.22 to 0.90. In addition, the measurement results indicated that an EVB can change the air velocity by about 50% based on the slat angle. Therefore, when an EVB is installed on a window opening, the effect of the EVB on wind-driven cross ventilation rate should be taken into account.

Considerable efforts have been made to reduce the energy consumption of buildings due to the energy crisis, and, the Korean government has supported the use of renewable energy through various grants. Among the possible renewable energy sources, geothermal energy can be used regardless of the outside weather. Therefore, energy consumption can be reduced considerably in summer and winter. Despite the increasing use of renewable energy, the use of renewables has not been operating appropriately. Therefore, this study examined some of the problems of the operation of renewable energy and some possible improvements. The aim of the study is to evaluate a building containing an actual installed multi-geothermal heat pump (Multi-GHP) system, in terms of the energy efficiency. In addition, this study evaluated the present control system and the method of complex operation regarding existing heat sources systems and GHP systems through a simulation. The results can be regarded as the result of a hybrid operation method for the improvement of an existing operation. Therefore, the Multi-GHP system energy use of a hybrid operation condition of the Multi-GHP systems and the absorption (ABS) chiller-heater system was reduced compared to the operation condition of the Multi-GHP system, and the total energy consumption of the heat source equipment was reduced. The proposed operation plan was evaluated after applying the system to a building. These results showed that the efficient operation of a Multi-GHP hybrid operation method is possible. As a result, the GHP energy use of Multi-GHP systems and the ABS chiller-heater system was reduced by 30% compared to existing operation and the total energy consumption of heat source equipment was reduced by 78%.

The energy performance of buildings depends on how effectively the building envelope responds to climate. Architects, therefore, need to design building envelopes with the consideration of local climate characteristics in the early design stage. Simplified formulas were used that evaluate the heating and cooling energy demand of building envelopes, which were applied to a model building with envelope and climate properties according to eight climate zones. Two climate indices, P and S , were developed. P enables the comparison of the heating and cooling energy demand of building envelopes, and S is for comparing the solar heat gain during heating and cooling seasons to review the feasibility of installing shading devices. The physical properties of envelopes were set differently according to the requirements in each climate zone proposed by American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) 90.1. Using local climate data, the P and S of 24 cities over eight climate zones in the United States were derived, which can be used to evaluate the heating and cooling energy characteristics of envelopes. The indices not only enable users to understand the characteristics of the local climate conditions in a simple manner, but also to carry out quantitative assessments on whether shading devices are feasible and, if so, what type is recommended.

Exterior dynamic shading devices, installed on “kinetic façades”, generate shaded areas of various shapes on windows according to the shape of the shading elements and the direction of their movement. The calculation of the shaded area is vitally important because it is directly related to the solar heat gain calculation process in building energy assessments. This paper dis-cusses a dynamic calculation method for deriving shaded fractions in consideration of the irregular shapes and unique movements of the shading elements in kinetic façades. The planar-polygon method was adopted for calculating accurate shaded areas on a window generated by irregularly shaped shade elements. To account for movements of the shading elements, the range of movement directions (i.e., rotating, sliding, etc.) was divided into steps of equivalent intervals. Applying these two methods, a shaded area calculation tool for the kinetic façade was developed. Three movement directions of shading devices were chosen for calculating shaded area, and the values of shaded fractions for six kinetic façade types were derived for different façade orientations during the summer and winter solstices. Lastly, to simplify the detailed calculation method, estimation equations for two types of kinetic façade were derived from a trend analysis of the shaded fraction values. This study deals with both detailed and simplified methods (estimation equation) for deriving the shaded fraction. The detailed method can be a more accurate solution in deriving the shaded fractions generated by complex exterior movable shading devices. However, the simplified method can be adopted in the early design stages to review various shading devices within a brief duration of time.

Terminal boxes are critical components of variable air volume (VAV) systems. The minimum airflow of terminal boxes is a key factor for comfort, indoor air quality (IAQ) and energy cost. If the minimum airflow is higher than needed, terminal boxes will have significantly more simultaneous heating and cooling, and airhandling units (AHUs) will consume more fan power. Buildings will have IAQ problems if the minimum airflow is less than required. A dynamic building simulation and field experiment were carried out to evaluate the impact of minimum airflow on the energy consumption of single duct VAV terminal boxes. A comparison of the simulation data with the existing constant setpoint revealed 20% annual energy savings in the improved minimum airflow setpoint. The experiment showed that the room air temperature can be maintained stably and the vertical difference in the room air temperature is lower than the comfort value when the improved minimum airflow setpoint is applied. Measurements of the CO 2 levels also revealed no indoor air quality problems.

In large-space buildings such as stadiums or auditoriums, the mean radiant temperature (MRT) at an occupant location can be highly influenced by the outdoor thermal environment. Therefore, it is important to measure a wide range of MRT values to understand thermal comfort conditions. This study reports on the importance of MRT measurements when evaluating the thermal environments of large spaces. We present an MRT measurement method using an infrared thermal imaging camera (IR camera). First, MRT measurement methods using a globe thermometer and an IR camera with respect to large spaces were reviewed to assess their limitations and applicability. A method utilizing the IR camera is proposed to enable multiple-location measurements with a single installed device. In our proposed method, MRT values are calculated using angle factors (AF) and mean indoor surface temperatures. A solid-angle algorithm is adopted to calculate AFs at occupant locations using a simplified model of the indoor space. To confirm the applicability of the proposed method in a large space, experiments are conducted using the globe-thermometer method and the IR method. Results show that the MRT values estimated from both methods have similar trends with differences less than 1 °C. The measured surface temperatures using an IR camera show that indoor surface temperatures can be significantly different according to location and orientation. We conclude that the proposed method is capable of measuring MRT values at multiple locations.

This study focuses on the control method of an adaptive façade based on the combination of energy conservation and glare protection strategies. Two evaluation methods are applied. The first method is to determine the appropriate range of the shading states which allow a minimum amount of solar heat gain and lighting energy requirements in the passive zone. The second method is to determine the appropriate range of the shading states which protect discomfort glare. Finally, the hourly optimal shading state can be determined through the two evaluation methods. In the application, the hourly shading state of a dynamic folding system, which can fully close and open in the two directions, is derived and evaluated. The results of the evaluation showed that the optimal control method suggested in this study can deliver effective energy performance while avoiding discomfort glare.

Modular technology has become a major issue of the construction industries to enhance their productivity. Modular bridge construction generally requires the collaboration between the contractor, designer, fabricator and constructor. Therefore, a readily accessible information model based on BIM technology should be provided for their communication during a construction project life-cycle. In this study, BIM based 3D information modeling was carried out for the modular bridge pier. First, the product breakdown structure (PBS) and level of detail (LOD) of the pier were defined. Based on them, 3D models were created by using parametric modeling method. In addition, database was constructed for the exchange of geometry and property data of 3D models. Finally, application areas of 3D information model were suggested, including the quantity estimation and the 4D simulation.

In buildings, windows have a major influence on space heating demand and indoor environment both with respect to climate and daylight. To reduce the window coefficient of the overall heat transmission, we use aerogel. Aerogels have a high surface area, low density, open pore structure, and excellent insulation properties. We mixed pressure sensitive adhesive and aerogel (10, 15, and 20 mass%) using a homogenizer. A mixture of the adhesives and silica aerogels attached film can reduce thermal conductivity. Silica aerogels are characterized by a surface area analyzer (BET), a Fourier transform infrared spectrometer, a thermogravimetry (TG) analyzer, and probe tack method. Thermal conductivity was measured by a TCi thermal conductivity analyzer.

An air-handler fan pressurization method can be used to estimate envelope leakage in large or tall buildings. However, because the size of a building's air handling system is designed only to meet heating and cooling demands, there may be many unanticipated problems on the test site. For example, one may not be able to perform the envelope leakage test using the air-handler fan(s) with relatively high pressure differences across the building envelope (e.g. 50–75 Pa). In order to overcome the insufficient pressurization problem faced during the air handler fan pressurization test, a combination of an air-handler and a blower door is proposed and verified by performing a pressurization test in an office building. It was found that the proposed method is feasible for evaluating the overall envelope leakage in large or tall buildings.

Reliable airtightness data is needed to estimate the level of air infiltration and the thermal loads for improving building energy efficiency and indoor comfort. While useful information on air lea...

In order to evaluate the thermal and lighting energy performance of a kinetic facade using external movable shading devices, it is important to consider the operation of the shading devices since it can influence the performance significantly. This study proposes a calculation methodology which assesses the performance of the movable shading devices with the consideration of the movements of the shading devices. Calculation methods were derived by which solar heat gain, lighting energy requirement, and the primary energy equivalent to heating and cooling energy requirement can be obtained. Using the calculation methods, the optimal operation scenario for the movable shading devices was presented which can minimize the solar heat gain and lighting energy requirement. A comparison case study was carried out to demonstrate the use of the control strategies for a building where a drop awning was installed. The results showed that the optimal operation scenario suggested in this study can potentially deliver effective energy performance. The limitation and applicability of the suggested method were also discussed.

The Korean government is fostering a renewable energy industry as a means of handling the energy crisis. Among the renewable energy systems available, geothermal energy has been highlighted as highly efficient, safely operable and relatively unaffected by outdoors air conditions. Despite the increasing use of renewable energy, the devices using renewables may not be operating appropriately. This study examined current problems in the operation of a geothermal heat pump (GHP) system. The efficiency of a geothermal heat pump system to studied to maximize the operation plan. Our study modelled the target building and analyzed the energy using TRNSYS, which is a dynamic energy simulation tool, to apply the coefficient of performance (COP) and evaluate the operation method. As a result, the GHP total energy consumption from the COP control method was reduced by 46% compared to the current operation. The proposed control method was evaluated after applying the system to a building. The results showed that efficient operation of a geothermal heat pump system is possible.