Development of Design Method for River Water Source Heat Pump System Using an Optimization Algorithm

: River water source heat pump (RWSHP) systems are being proposed to reduce the energy consumption and carbon emissions of buildings. The RWSHP system is actively applied to large-scale buildings due to its stable performance. The application of RWSHP in large-scale facilities requires an accurate capacity design with considerations of building load, heat source, and environment conditions. However, most RWSHP systems are over-designed based on peak load of buildings. These design methods, based on peak loads, are economically and environmentally disadvantageous. Therefore, this paper aims to development an optimal design method, both economically and environmentally, for the RWSHP system. To develop this optimal design method, a simulation model was created with an optimization algorithm. The economics of the RWSHP system were calculated bases on present worth of annuity factor. Moreover, CO 2 emissions were estimated using the life cycle climate performance proposed by the International Institute of Refrigeration. The total cost of the proposed RWSHP system that apply the optimum design method decreased by 24% compared to conventional RWSHP systems. Moreover, CO 2 emissions of the proposed RWSHP system reduced by 4% compared to conventional RWSHP systems.


Introduction
According to the International Energy Agency (IEA) report of 2019, energy consumption in the building sector accounts for approximately 35% of the total energy consumption [1].Carbon dioxide emissions in the construction sector accounts for 38% of the total carbon dioxide emissions [1].Furthermore, electricity consumption in building operations accounts for approximately 55% of the global electricity consumption [1].Therefore, reducing the energy consumption and carbon dioxide emissions in the construction sector is of vital importance.Renewable energy can be used to address this issue since they can lead to a reduction in fossil fuel usage, energy consumption generated from fossil fuels, and greenhouse gas (GHG) emissions.However, stable power generation from renewable energy sources cannot be guaranteed and are greatly affected by the climate.Therefore, to resolve this challenge technology is required that increases the power generation efficiency of renewable energy or ensures its stability.As of 2019, energy production from waste, biomass, and solar energy accounts for more than 80% of renewable energy production in South Korea.This implies that there is a low production of geothermal, hydrothermal, and solar heat using temperature differential energy [2].Therefore, it is necessary to investigate the performance and efficiency of temperature differential energy systems.
A natural water source heat pump (NWSHP) system uses river, sea, pond, and other water sources for renewable energy systems, is a technology that performs heating and cooling using temperature difference energy.NWSHP offers excellent performance and efficiency due to its abundant heat source, compared to the ambient temperature, which uses low and high-temperature heat sources during the cooling and heating seasons, respectively.However, a NWSHP system requires the additional installation of an intake facility and water pipe installations that can result in increased total installation costs.In addition, substantial initial investment costs may be required for large-capacity building that are based on the existing maximum load design method and incorporating the RWSHP system due to overcapacity.
NWSHP systems have been actively researched, alongside the expansion in their application, with a particular focus on system demonstration and performance prediction.Cho and Yun [3] investigated the heating and cooling performance of a heat pump (HP) using the heat energy of raw water supplied to a water treatment facility.In this study, the system performance was confirmed to decrease due to the load reductions resulting from seasonal changes.Zheng et al. [4] investigated the system efficiency of seawater source heat pump systems in areas experiencing severely cold winters.The results of the study indicated that the beach well infiltration intake systems for the seawater source heat pump system had the highest thermal performance during the entire heating season in areas of severely cold winters.Wang et al. [5] proposed a new type of groundwater-source system with a fresh air pre-conditioner.This system ensures low-grade energy stored in ground water is maximally used to reduce the energy consumption of the entire HP system and buildings.Jung et al. [6] comprehensively analyzed the NWSHP system feasibility by evaluating its performance after estimating the thermal energy of river water.The potential was measured by calculating the thermal energy of river water after measuring the water temperature and flow rate.Jung et al. [7] analyzed the NWSHP system performance, environmental impact, and economic feasibility in consideration of the following three variables: building type, water source, and intake distance.Oh et al. [8] reported on a rawwater source heat pump for a thermal storage tank in a vertical water treatment building that was dynamically simulated by TRNSYS.The condensing temperature of the HP, and heating and cooling load changes, according to the heat storage tank (HST) size, confirmed that an optimal size for the heat storage tank existed.
However, most of research focused on the HP and system performances rather than building loads and design conditions.Few studies have investigated the optimum design of an RWSHP system with real-world building applications.Therefore, the building conditions, loading, and system capacity calculation method are required when designing the system.Generally, the system capacity is determined by the building peak load with consideration to climate and standard building design conditions.This causes inefficient performance during seasonal changes and partial load conditions.In addition, conventional design methods result in system overdesign and energy overconsumption, which increases the total investment cost and CO 2 emissions.Therefore, the need for strategies to reduce both total investment costs and emissions, while actively reducing energy demands in the built environment have been emphasized.High-efficiency equipment must be used to save energy and reduce total investment costs and greenhouse gas emissions.However, the maximum energy saving effects cannot be solely achieved by using high-efficiency equipment.Therefore, an optimized design for a system considering energy production, economic, and environmental impacts is required.
Furthermore, several studies have been conducted on system optimization to prevent excessive energy consumption and reduce CO 2 emissions through optimization factors such as system design and operation.Du et al. [9] performed an area optimization for solar collectors in a combined solar heating system to minimize the unit cost of adsorption desalination.They investigated the influence of the auxiliary energy source prices, adsorption desalination performance ratios, and heat loss of solar heating on the optimization results.Mi et al. [10] performed a study on the optimization of a hot water supply system for public baths using a multi-heat source system combined with solar heat and a NWSHP.Comprehensive energy efficiency and economic evaluations were conducted on the model to reduce the total cost.Furthermore, Moon et al. [11] conducted a study on the optimal design method of ground source heat pump (GSHP) using an optimization method.Life Energies 2022, 15, 4019 3 of 19 cycle cost (LCC) analysis was performed on the initial and operating investment costs for 10 and 20 years, assuming that the optimal system design elements were HP, ground heat exchanger, and HST.Based on these analysis results, the three algorithms were confirmed to reduce the GSHP system costs in comparison with the existing design.Guo et al. [12] presented a method for designing a HST for a heating system using a NWSHP that incorporated solar heat as the auxiliary heat source.The analysis results confirmed that the HST volume could increase the energy saving rate of the system by 34.39%, and the optimum design for the HST volume was required to improve the system performance.Liu et al. [13] proposed an optimization model that could simultaneously optimize the system configuration, equipment capacity, and operating parameters for a multi-source complementary heating system.The results indicated that the LCC of the optimal mentary heating system based on air source heat pump could be reduced by 26.8%, whereas the seasonal coefficient of performance (COP) could be increased by 398.8%, compared to coal-fired boiler heating system.In addition, zero CO 2 emissions could be achieved during the operation of mentary heating system based on air source heat pump.Table 1 shows the previous studies regarding optimization of energy system.The previous studies investigating optimization focused on improving performance and economic efficiency by optimizing system variables.Recent studies have been conducted on measures to reduce CO 2 emissions or GHG, whereas studies on the reduction of these factors by optimization or the optimization of new and renewable energy systems with regards to building loads are lacking.Therefore, changes in performance, economic feasibility, and CO 2 emissions must be analyzed by optimizing the system capacity in consideration of building loads.
In this study, a design method using an optimum algorithm was developed to suitably design an HP system with river water as the heat source.The capacity of the RWSHP system applied to a large-scale building was calculated using the Hooke-Jeeves algorithm for efficient capacity design with regards to LCC.Furthermore, the environmental impact was evaluated and compared to the developed and conventional methods.For economic evaluation, the present worth of annuity factor (PWAF) analysis method was used and CO 2 emissions were calculated using the life cycle climate performance (LCCP) analysis method proposed by the International Institute of Refrigeration (IIR) [14].The results of this study could be a fundamental source for the suitable design of NWSHP systems.Moreover, this study makes a novel contribution to literature on low-cost and low-emission technology of NWSHP systems for large-scale buildings.

Natural Water Source Heat Pump System
Figure 1 shows a schematic of the NWSHP system.The NWSHP system uses heat sources such as river, sea, and pond water, and raw-water pipes.The system executes a heat exchanger using a water heat source that is more efficient than ambient temperature to heat and cool buildings.Heat exchange was performed by the heat exchanger using the water heat source from the intake pipe to conduct the heating/cooling throughout the building by an HP with a fan coil unit (FCU), radiant panel, and other equipment.

Natural Water Source Heat Pump System
Figure 1 shows a schematic of the NWSHP system.The NWSHP system uses heat sources such as river, sea, and pond water, and raw-water pipes.The system executes a heat exchanger using a water heat source that is more efficient than ambient temperature to heat and cool buildings.Heat exchange was performed by the heat exchanger using the water heat source from the intake pipe to conduct the heating/cooling throughout the building by an HP with a fan coil unit (FCU), radiant panel, and other equipment.The investigation results of international cases of NWSHP are as follows: • The Lotte World Tower in South Korea uses raw-water pipes for heating and cooling; approximately 10-20% of the heating and cooling energy is managed by the 3000-RT NWSHP system that saves approximately 700 million in energy cost annually [15].

•
Cooling and heating are conducted to office buildings and 100 houses of approximately 278,000 m 2 in a district located in Hakozaki, Japan using river water.The river water originates from the Sumida River, which is characterized by average temperatures of approximately 25 and 8 °C in summer and winter, respectively [15].

•
Cornell University in the USA developed and utilized a 20000-RT cooling and heating system using pond water with a temperature range of 5-13 °C at a depth of 76 m that reduces the energy consumption by approximately 80% [16].
These international cases confirmed that buildings have been designed to respond to the peak load or that a certain amount of heating and cooling energy can be supplied by NWSHP systems.In other words, the optimum capacity design corresponding to the building load was not achieved.

Research Method
In this study, an optimum capacity design method was proposed for the introduction of an RWSHP system.A comparative analysis was performed using existing design methodology.For comparison, an RWSHP system was built using a dynamic simulation and an optimal design was performed using an algorithm.Figure 2 shows the proposed design method simultaneously using energy simulation and the optimization algorithm.The input data included building and system operating conditions, and cost information.These were used to calculate the peak load and initial investment cost of the building according to the existing design method.The RWSHP system was created using a dynamic energy simulation, and thereafter an objective function was designated according to the designer's purpose by an algorithm to optimize the system performance, economic The investigation results of international cases of NWSHP are as follows: • The Lotte World Tower in South Korea uses raw-water pipes for heating and cooling; approximately 10-20% of the heating and cooling energy is managed by the 3000-RT NWSHP system that saves approximately 700 million in energy cost annually [15].

•
Cooling and heating are conducted to office buildings and 100 houses of approximately 278,000 m 2 in a district located in Hakozaki, Japan using river water.The river water originates from the Sumida River, which is characterized by average temperatures of approximately 25 and 8 • C in summer and winter, respectively [15].

•
Cornell University in the USA developed and utilized a 20000-RT cooling and heating system using pond water with a temperature range of 5-13 • C at a depth of 76 m that reduces the energy consumption by approximately 80% [16].
These international cases confirmed that buildings have been designed to respond to the peak load or that a certain amount of heating and cooling energy can be supplied by NWSHP systems.In other words, the optimum capacity design corresponding to the building load was not achieved.

Research Method
In this study, an optimum capacity design method was proposed for the introduction of an RWSHP system.A comparative analysis was performed using existing design methodology.For comparison, an RWSHP system was built using a dynamic simulation and an optimal design was performed using an algorithm.Figure 2 shows the proposed design method simultaneously using energy simulation and the optimization algorithm.The input data included building and system operating conditions, and cost information.These were used to calculate the peak load and initial investment cost of the building according to the existing design method.The RWSHP system was created using a dynamic energy simulation, and thereafter an objective function was designated according to the designer's purpose by an algorithm to optimize the system performance, economic feasibility, and environmental impact.The calculated results indicated the optimized capacity of system components, system performance, economics, and CO 2 emissions; the optimal system capacity could be altered according to the objective function.Figure 3 shows the flow chart of the existing and proposed design methods.The existing design method determined the capacity of the equipment required for the system by considering the safety factor after analyzing the building peak loads.If necessary, the feasibility analysis was then performed using performance or economic analyses.The existing design method was intended to respond to the building peak loads.However, this can result in excessive energy consumption and GHG emission increase, and has an added disadvantage of high initial investment costs since a system with excessive capacity operates with season changes and partial loads.
The design method proposed in this study analyzes the building peak loads similarly to the existing design method, and thereafter determines the equipment capacity required by the system by considering the safety factors.The reduction in total investment cost, increase in system performance, or reduction in CO 2 emissions were then designated as objective functions.The design was performed by calculating the optimal HP and HST capacities for the objective function by applying a suitable algorithm to the objective function.
As a result, the proposed optimal design method employed an appropriate system capacity to respond to the designer's objectives by specifying an objective function design, rather than considering the maximum load and safety factor.That is, the proposed method offsets the challenges associated with the existing design method that considers the safety factor to improve the total investment cost and CO 2 reduction.

Building Load Model
In this study, large offices in Busan, South Korea, were selected as target buildings for the RWSHP system, assuming that the entire building would be heated and cooled by the RWSHP system.The building load model was established based on the large-scale office model provided by the Pacific Northwest National Laboratory (PNNL) [17].Figure 4 shows the shape and size of the large-scale office building.

Building Load Model
In this study, large offices in Busan, South Korea, were selected as target buildings for the RWSHP system, assuming that the entire building would be heated and cooled by the RWSHP system.The building load model was established based on the large-scale office model provided by the Pacific Northwest National Laboratory (PNNL) [17].Figure 4 shows the shape and size of the large-scale office building.
The specifications of the building load model are listed in Table 2. Busan, South Korea, was assumed to be the geographical region and the associated weather data were input.The indoor setpoint temperatures for the heating and cooling periods were set to 21 and 24 • C, respectively.Infiltration was set according to the infiltration calculation method of the American Society of Heating, Refrigeration and Air-conditioning Engineers (ASHRAE) 90.1 [18].The heating season was set at eight months from January to May and October to December, and the cooling season was set at four months from June to September.The lighting, equipment, and occupancy schedules were divided into weekdays, Saturdays, and Sundays [19].
The internal heat gains of lighting and equipment were calculated based on ASHRAE 90. .The quantity of heat generated by lighting and the device were set at 11.84 and 8.07 W/m 2 , respectively [19].Furthermore, the internal heat gains of the human body were calculated based on the quantity of heat generated by the occupants while typing lightly [20]; the total number of occupants was presumed to be 192, calculated as one person per 18.5 m 2 [21].
The thermal conductivities of the walls, roofs, windows, and floors were based on energy-saving design standards for buildings in South Korea [22].
Figure 5 shows the set-point temperature boundary, ambient temperature, and monthly building loads in Busan, South Korea.The load analysis was conducted on the first floor, which had the largest heating and cooling loads.The load analysis results indicated that the heating load was dominant in comparison with the cooling load.The heating and cooling peak loads of the building were calculated as 45 and 35 kW, respectively.
U-Value (W/m Figure 5 shows the set-point temperature boundary, ambient temperature, and monthly building loads in Busan, South Korea.The load analysis was conducted on the first floor, which had the largest heating and cooling loads.The load analysis results indicated that the heating load was dominant in comparison with the cooling load.The heating and cooling peak loads of the building were calculated as 45 and 35 kW, respectively.

Dynamic Energy Simulation
Figure 6 shows a schematic of the RWSHP system designed in this study.The input and output elements, and values calculated by the simulation are indicated.The system configuration was based on the schematic.The system components consisted of a heat exchanger, water-to-water HP, HST, FCU, pump, and large-scale office building; the capacity of each system device was designed by considering the peak load and performance curve of the building.The components and values of the constructed RWSHP system are listed in Table 3.

Dynamic Energy Simulation
Figure 6 shows a schematic of the RWSHP system designed in this study.The input and output elements, and values calculated by the simulation are indicated.The system configuration was based on the schematic.The system components consisted of a heat exchanger, water-to-water HP, HST, FCU, pump, and large-scale office building; the capacity of each system device was designed by considering the peak load and performance curve of the building.The components and values of the constructed RWSHP system are listed in Table 3.The HP capacity was obtained by adding a safety factor of 20% to the building peak load.The HP power consumption was set by referring to the water-to-water HP of the climate master [23].The hot water supply load was not considered for HST, and the FCU was set by adding a safety factor of 20% to the HP capacity.The volumetric air flow rate and power consumption of the FCU were set by referring to Carrier's ceiling-type FCU [24].The pump was set within the appropriate flow rate range specified by ASHRAE [25].The validity of the numerical analysis model was confirmed in previous study by comparing the simulation and experimental results [26].

System Operation Method
Figure 7 shows the operational logic of the simulation.During the heating season, the operation was stopped when the temperature of the water intake was 2 • C or less, and when the temperature of the drain river water was 0 • C or less, the freezing of the heat exchanger due to the heat source was considered.Thereafter, the HP was set to operate to maintain the set point temperature (50 • C) of the HST that sent circulating water to the FCU for heating.In contrast, during the cooling season, the effect of the heat source was not considered.Therefore, the HP was set to operate to maintain the set point temperature (12 • C) of the HST that sent circulating water to the FCU to perform cooling.

Optimization Algorithm
In this study, optimization was performed using Hooke-Jeeves algorithms; these are frequently used for facility design due to the optimum capacity design, ability to incorporate multiple variables, and fast calculation time.Figure 8 shows the calculation method of the Hooke-Jeeves algorithms and the actual optimization process according to the objective function.

Objective function, f x = Initial investment cost + 20 years of operational costs
(1)

Optimization Conditions
The minimum, maximum, and step values of the optimization variables used in the algorithm are listed in Table 4.The maximum and minimum values of the HP heating and cooling capacity were presumed to be 50% and 200%, respectively, calculated by the existing design method.In the case of HST, 10% and 200% of the volume calculated by the existing design method were presumed.The parameter values (power consumption of the pump, flow rate, etc.) changed according to the change in capacity of HP and HST, which were altered by linear/nonlinear equations.The objective function of optimization is a function of the total investment cost.The HP cooling and heating capacities and HST volume, which significantly affect the initial investment and annual operation costs, were set as the optimization variables for the function formula.The total investment cost was calculated from the initial investment cost that consisted of the HP and HST equipment, and operational costs over 20 years.The function is expressed as: The Hooke-Jeeves algorithm is a pattern search method that solves the nonlinear optimization problem.As shown in the optimization process, this method determines the optimal value by finding the point that best fits the objective function among neighboring points from the initial value.This is repeated until a point is found to fit the objective function [11,27].

Optimization Conditions
The minimum, maximum, and step values of the optimization variables used in the algorithm are listed in Table 4.The maximum and minimum values of the HP heating and cooling capacity were presumed to be 50% and 200%, respectively, calculated by the existing design method.In the case of HST, 10% and 200% of the volume calculated by the existing design method were presumed.The parameter values (power consumption of the pump, flow rate, etc.) changed according to the change in capacity of HP and HST, which were altered by linear/nonlinear equations.For the optimization constraint during the heating/cooling season, an error was set to occur when it deviated significantly from the normal leaving water temperature (LWT) range of HP.In addition, an error was set to occur when 300 h was exceeded, based on the unmet load hour of ASHRAE 90.1-2004 [19].

Economic Analysis
The initial investment and operation costs over 20 years were calculated for the economic analysis.Table 5 lists the prices used to calculate the initial investment and annual operational costs of HP and HST.The equipment costs of HP and HST, which correspond to the initial investment cost, were set according to the Korean online e-procurement system [28].The annual operational cost was calculated by multiplying the electricity consumption by the electricity rate per kW, and the operational cost over 20 years was calculated using the PWAF method.The cost of electricity was determined from the electricity rate data obtained from the Korea Electric Power Corporation [29].However, the costs for the secondary load and facility construction, which were the same as those of the existing design method, were excluded.Based on the price information above, the calculation was performed using the following formula: C HST = C Price per cubic meter HST × Heat storage tank volume C HP = C Price per kW HP × Heat pump capacity (5)

Life Cycle Climate Performance
In this study, an LCCP analysis was performed for the RWSHP system, which aspires for economy and low CO 2 emissions.The HP CO 2 emissions calculated by the existing design method and HP calculated by the optimum design method were compared using the actual HP model of the climate master [23].The LCCP analysis method was classified into direct and indirect emissions, as provided by IIR.The direct emissions calculated the CO 2 emissions based on the annual refrigerant leakage, refrigerant leakage during disposal, and equipment lifespan.The indirect emissions were calculated based on CO 2 emissions during equipment manufacturing from the annual energy consumption, and from the refrigerant [14].

Result of Optimization
Table 6 lists the optimization results obtained in this study.Figure 9 shows a comparison of the results.The comparison between the HP capacity and the existing design method is as follow.In the optimization design method, the cooling and heating capacities decreased by 25% and 22%, respectively.Further, the volume of HST decreased by 72% compared to the existing design method.The HP capacity decreased to an appropriate capacity rather that that corresponding to the peak building load.Moreover, the calculated capacity was suitable for performing the role of a buffer tank since the HST capacity did not consider the hot water supply load.Figure 9 shows that the appropriate HP and HST capacities were calculated by the RWSHP system optimization using an algorithm; the unmet load hour was set at less than 300 h, and the feasibility of the results were verified.However, the change in climatic conditions were not considered, and thus the unmet load hour may exceed 300 h depending on the simulated climate.In addition, the simulation was conducted by assuming two constraints; however, the need to consider other variables, such as resident dissatisfaction and equipment overload, was confirmed.

Result of Energy Performance
Figure 10 shows the comparative analysis results between the energy performance of the system employing the existing and optimum capacity design methods.Since the HP and HST capacities were calculated to correspond to the peak loads in the existing design method, the existing method was confirmed to have a higher performance than the optimal capacity design method.Based on the energy performance analysis results, the HP and system COP exhibited an average decrease of 3% and 5%, respectively, compared to the existing design method.Although the optimum method had a similar COP value for the performances of HP and system during the heating season, the HP and system COP of the conventional design method were the higher during the cooling season.
The HP and system COP decrease when the equipment capacity is calculated using the optimal capacity design method.This decrease occurs since the maximum load design method reduces the operating time and frequency of operation as a result of the large equipment capacity, regardless of the system's ability to accurately respond to the building loads.Although the optimal capacity design method is designed with an appropriate device capacity, the operation time and operation frequency must be increased compared to the existing design method to cope with the same load.In addition, when separated by heating and cooling periods, there is little difference in the performance during the heating period; however, there is a difference in performance during the cooling period.This is explained by the 20 or 5 • C differences between the heat source temperature and indoor set temperature during the heating or cooling periods, respectively.In the case of a large-capacity system during the cooling period, the power consumption reduces due to the low operation frequency, and thus the performance was excellent.

Result of Economic Analysis
The economic effects of the optimum capacity design method were analyzed.The results indicated that the existing design method had the most investment in both the initial investment and operation costs over a 20-year period.The total investment costs for the original and optimization cases are shown in Figure 11.Table 7 shows the objective function of the initial and total investment costs, and the annual operation cost factors after optimization.Compared to the existing design method, the optimization case had an initial investment cost decrease of 59%, operational cost over Energies 2022, 15, 4019 15 of 19 20 years decrease of 11%, and total cost decrease of 24%.Due to the reduction in the heating and cooling capacities of HP and HST volumes, the reduction effect of both the initial investment and operational costs over 20 years was confirmed.In terms of initial investment cost, the heat storage tank had the highest price according to unit capacity.Since HST was used as a buffer tank and not for heat storage, the optimization of the HST volume was clearly observed.The volume optimization was a large influencing factor, whereas the decrease in the other factors was small.In terms of operating cost over 20 years, the system capacity and power consumption were large despite the low operating frequency of the existing design method.Although the operating cost over 20 years was considered, the operating cost of the optimal capacity design method could be lower for a period shorter than 20 years.However, the result has a limitation in that an economic evaluation was not conducted for all system components; an economic evaluation was only performed on the major factors of total investment costs.

Result of Life Cycle Climate Performance
An LCCP analysis of the optimum capacity design method was performed.The LCCP analysis was conducted by specifying the existing design method and HP of optimization case.The product was analyzed based on the water-to-water HP model of the climate master [23].
Table 8 lists the values used for LCCP.These were prepared by referring to the values in the IIR LCCP report [14].The CO 2 emission coefficient provided by the Korea Power Exchange was used for the CO 2 emissions caused by electricity usage [29].
Figure 12 shows the LCCP analysis results.Based on the CO 2 emission analysis, the existing and optimum design methods exhibited CO 2 emissions of approximately 105,259,800 and 101,956,900 kg CO 2 e, respectively.Consequently, an approximately 4% CO 2 emission reduction can be achieved by the optimal capacity design method.The optimal capacity design produced the greatest reduction in CO 2 emissions from energy use.In addition, the CO 2 emissions generated during equipment production decreased.Although all of the CO 2 emissions produced by the system could not be considered, CO 2 emissions were reduced by applying the optimal capacity design method.

Conclusions
In this study, a design method was developed to assess the system capacity using an optimum algorithm and considering LCC and LCCP.The comparison analysis was conducted between the existing and developed design methods for an RWSHP system.The validity of the optimal capacity design was verified by comparing annual simulation results.The study results are summarized as follows: 1.
The initial investment and annual operation costs were reduced by 59% and 11%, respectively, by applying the optimization design method to the RWSHP system.The optimization design method decreased the total cost by 24% compared to the existing design method.

2.
The optimization design method decreased the CO 2 emission by approximately 4% compared to the existing design method for the RWSHP system.The LCCP analysis was only evaluated for HP.However, should all system components be included in the LCCP analysis, the difference in CO 2 emission is estimated to further increase.

Figure 1 .
Figure 1.Schematic of natural water source heat pump system.

Figure 1 .
Figure 1.Schematic of natural water source heat pump system.

Figure 1 .
Figure 1.Schematic of natural water source heat pump system.

Figure 3 .
Figure 3. Flow charts of the conventional and proposed optimum design methods.

Figure 3 .
Figure 3. Flow charts of the conventional and proposed optimum design methods.
(a) Ambient air temperature in Busan, South Korea (b) Monthly loads

Figure 5 .
Figure 5. ambient air temperature and monthly building loads in a large office building; (a) ambient air temperature in Busan, South Korea and setpoint temperature boundary of large office building, (b) monthly heating and cooling loads in a large office building in Busan, South Korea.

Figure 5 .
Figure 5. Ambient air temperature and monthly building loads in a large office building; (a) ambient air temperature in Busan, South Korea and setpoint temperature boundary of large office building, (b) monthly heating and cooling loads in a large office building in Busan, South Korea.

Figure 6 .
Figure 6.Schematic of the simulation.

Figure 7 .
Figure 7. Operational logic of the RWSHP system.

Figure 7 .
Figure 7. Operational logic of the RWSHP system.

Figure 8 .
Figure 8. Hooke-Jeeves algorithm calculation method and actual optimization process; (a) calculation method of the Hooke-Jeeves algorithms, (b) actual optimization process according to the objective function.

Figure 8 .
Figure 8. Hooke-Jeeves algorithm calculation method and actual optimization process; (a) calculation method of the Hooke-Jeeves algorithms, (b) actual optimization process according to the objective function.
(a) HP Heating and cooling capacity (b) HST volume

Figure 9 .Figure 11 .
Figure 9.Comparison of optimization results; (a) comparison result between the existing design method and the optimized design method for the heating and cooling capacity of heat pumps, (b) comparison result between the existing design method and the optimized design method for the volume of the heat storage tank

Figure 9 .
Figure 9.Comparison of optimization results; (a) comparison result between the existing design method and the optimized design method for the heating and cooling capacity of heat pumps, (b) comparison result between the existing design method and the optimized design method for the volume of the heat storage tank.

Energies 2022 ,
15, x FOR PEER REVIEW 5 of 6 (a) HP Heating and cooling capacity (b) HST volume

Figure 9 .Figure 10 .
Figure 9.Comparison of optimization results; (a) comparison result between the existing design method and the optimized design method for the heating and cooling capacity of heat pumps, (b) comparison result between the existing design method and the optimized design method for the volume of the heat storage tank

Figure 10 .
Figure 10.Results of HP and system performance; (a) HP COP of the system built with the existing and optimal capacity design methods, and (b) system COP.Figure 10. Results of HP and system performance; (a) HP COP of the system built with the existing and optimal capacity design methods, and (b) system COP.

Energies 2022 ,Figure 9 .
Figure 9.Comparison of optimization results; (a) comparison result between the existing design method and the optimized design method for the heating and cooling capacity of heat pumps, (b) comparison result between the existing design method and the optimized design method for the volume of the heat storage tank

Figure 11 .
Figure 11.Result of economic analysis; (a) comparison of economic feasibility between the existing design method and the optimized design method for initial investment cost, energy cost, and total cost.(b) LCC comparison analysis result between the existing design method and the optimized design method

Figure 11 .
Figure 11.Result of economic analysis; (a) comparison of economic feasibility between the existing design method and the optimized design method for initial investment cost, energy cost, and total cost.(b) LCC comparison analysis result between the existing design method and the optimized design method.
(a) Result of economic feasibility (b) Result of LCC analysis

Figure 11 .Figure 12 .
Figure 11.Result of economic analysis; (a) comparison of economic feasibility between the existing design method and the optimized design method for initial investment cost, energy cost, and total cost.(b) LCC comparison analysis result between the existing design method and the optimized design method

Table 1 .
Previous studies regarding optimization of energy system.

Table 2 .
Specifications of building load model.

Table 3 .
Specifications of simulation components.

Table 3 .
Specifications of simulation components.

Table 5 .
Specifications of economic analysis.

Table 4 .
Specification of optimization.

Table 4 .
Specification of optimization.

Table 5 .
Specifications of economic analysis.

Table 6 .
Optimization results of capacity and volume.

Table 7 .
Optimization results of costs.

Table 8 .
Specifications of life cycle climate performance.