1 s2.0 S1359431116335037 Main
1 s2.0 S1359431116335037 Main
1 s2.0 S1359431116335037 Main
Research Paper
Dynamic target high pressure control of a VRF system for heating energy
savings
Geun Young Yun a,⇑, Je Hyeon Lee b, Inhan Kim c
a
Department of Architectural Engineering, Kyung Hee University, 1732, Deogyeong-daero, Giheung-gu, Yongin si, Gyeonggi-do 17104, South Korea
b
Department of Digital Appliances R&D Team, Samsung Electronics, 129, Samsung-ro, Yeongtong-gu, Suwon-si, Gyeonggi-do 443-742, South Korea
c
Department of Architecture, Kyung Hee University, 1732, Deogyeong-daero, Giheung-gu, Yongin si, Gyeonggi-do 17104, South Korea
h i g h l i g h t s
a r t i c l e i n f o a b s t r a c t
Article history: Variable refrigerant flow (VRF) systems are widely used because of their ability to provide individualized
Received 20 July 2016 comfort control with energy-saving potential. This study develops load responsive high pressure control
Revised 18 November 2016 of a VRF system with the aim of reducing the heating energy consumption of a VRF system under part
Accepted 20 November 2016
load conditions. The developed control consists of two parts: one part that determines the level of heating
Available online 22 November 2016
load, and the other that assigns a target high pressure based on the level of heating load. In this way, the
compressor speed can be accurately matched to heating load, which improves the energy performance of
Keywords:
the VRF system. A series of multicalorimeter experiments revealed that the heating capacity of the VRF
Building energy
Heat pump
system varied by 45% by modulating the target high pressure and that its efficiency was enhanced by 21%
Energy efficiency by changing the high pressure from 30 kgf/cm2 to 25 kgf/cm2. To evaluate the annual heating energy per-
Simulation formance of the VRF system with the developed control, a custom computer code was developed to
EnergyPlus implement the developed control using a programming language called EnergyPlus Runtime Language.
Variable refrigerant flow system Simulation outcomes showed that the annual heating energy consumption of a medium-size office build-
ing was reduced by 22% when the developed control was applied.
Ó 2016 Elsevier Ltd. All rights reserved.
http://dx.doi.org/10.1016/j.applthermaleng.2016.11.148
1359-4311/Ó 2016 Elsevier Ltd. All rights reserved.
G.Y. Yun et al. / Applied Thermal Engineering 113 (2017) 1386–1395 1387
the improvement of their energy efficiency is critical to decrease 2. Thermal load based control of high pressure in a VRF system
electricity demands in winter.
One feature of the VRF system is that it enables concurrent The developed control in this study consists of two parts: one
heating and cooling modes of different indoor units through mode part that determines the level of heating load, and the other that
change units or heat recovery units [9,10]. A key technology is to assigns a target high pressure based on the level of heating load.
vary refrigerant flow rates between outdoor units and multiple A unique aspect of the control is the direct modulation of the set-
indoor units [11,12]. Individual zone controls, along with seamless point of high pressure to match the heating loads of indoor units.
switchover between heating and cooling operation modes are main The existing studies [26–28] focused on the control of the com-
reasons for the wide use of VRF systems in various buildings such pressor speed to change the heating capacity for a fixed high tem-
as offices, schools, retail stores, and hotels [13]. perature. This study utilizes the high pressure as a control variable,
Prior research on VRF systems has focused on reducing electric- so that the heating capacity can be adjusted with high accuracy to
ity consumption in cooling mode. Simulation and experimental varied heating demands.
results indicated that the water-cooled VRF system considerably The first part of the developed control determines the capacity-
lowered the amount of energy used by the compressor in the cool- weighted part load condition of indoor units in the single refriger-
ing mode [14,15]. Several studies [16–19] have confirmed the ant circuit (Fig. 1). This part to decide the thermal load of indoor
importance of control of electronic expansion value (EEV) opening units is based on the method developed by Yun et al. [24]. The
and compressor control of the part load cooling efficiency of VRF algorithm is applied to a part load condition when the deviation
systems. Tu et al. [16] found that the average suction pressure of (DTrs) of room temperature from the setpoint temperature is less
refrigerant was a reliable variable to control the cooling capacity than 3 K. When the deviation is within 3 K, the algorithm evaluates
of the compressor. Choi and Kim [17] experimentally showed that that the heating load is low enough to apply the developed control.
cooling capacity can be effectively controlled by modulation of EEV The algorithm estimates the two factors of the current indoor unit:
opening to meet the cooling load of two indoor units with a super- Ld1 and Ld2. Ld1 is determined based on the change in the temper-
heat temperature of 4 °C. Control algorithms for multi-evaporator ature of the room over the past five minutes (DTrm) and the crite-
air conditioning systems that regulate the cooling load and EEV rion for the room temperature change (Trc):
openings to determine compressor speed have also been developed
T rc
[18,19]. Ld1 ¼ ð1Þ
There are few experimental studies of the performance of VRF T rm
systems in heating mode. Joo et al. [20] found that the control of Trc can be determined from Table 1 [29]. In the next step, the
the compressor speed and the EEV openings was a key factor in runtime fraction of the indoor unit (Ld2), which is the ratio of
regulating heating capacity. Jeong and Choi [21] investigated the the time that current indoor unit is on over the past thirty minutes
effects of EEV opening and compressor speed on the heating per- is calculated. Then, the smaller value between Ld1 and Ld2
formance of the heat pump and found that heating capacity becomes the heating load coefficient of the current indoor unit
increased with an increase in EEV opening and compressor speed. (Ld(n)). Then, the heating load coefficients for all the indoor units
The setpoint of evaporating temperature remains unchanged in in the refrigerant circuit is calculated. Finally, the capacity
general VRF controls [22], although a rise in evaporating tempera- weighted average heating load coefficient (Ld) is calculated using
ture results in the efficiency of the VRF system [23]. Recent studies Eq. (2), as follows:
[24,25] revealed the large potential of controlling evaporating tem- Pm
perature to improve the cooling efficiency of VRF systems espe- n¼1 ðLdðnÞ CapacityðnÞ Þ
Ld ¼ Pm ð2Þ
cially for part load conditions because the cooling capacity could n¼1 CapacityðnÞ
be accurately adjusted to cooling load conditions.
Here, m is the total number of indoor units in the refrigerant circuit
The heating capacity of the VRF system needs to be adjusted in
and Capacity(n) is the heating capacity of the indoor unit, n.
order to match the varied heating loads of indoor units in the VRF
The second part of the algorithm decides the dynamic target
system. The common way of modulating the heating capacity of
high pressure as a function of the capacity-weighted average heat-
the VRF system is to change its compressor speed [9]. Awan [26]
ing load coefficient (Ld), as shown in Table 2. A series of experi-
experimentally investigated the change in heating capacity by
ments were carried out to determine the range of target values.
adjusting the compressor speed from 30 Hz to 90 Hz. Kim et al.
Both the stability of the heating capacity with various loads and
[27] proposed the design of the variable speed heat pump that
the pressure range of the compressor with the refrigerant,
could response to direct load control signal. Tu et al. [28] devel-
R410A, were considered to determine the target values. As the
oped a control model for the VRF system, which determined the
heating load decreased, the high pressure of the compressor in
optimal range of compressor frequency. However, the high pres-
heating mode decreased from 30 kgf/cm2 to 25 kgf/cm2. The rela-
sure is fixed in the heating capacity control of the VRF system in
tionship between the discharge refrigerant temperature (Tds in
a normal operation [22], although the change in the high pressures
°C) in Table 1 and high pressure (Ph in kgf/cm2) is given by Eq. (3).
has a direct influence on the heating capacity. Thus, in this study
the high pressure was explored as a control variable in the adjust- T ds ¼ 0:0173ðPh Þ2 þ 2:4443ðPh Þ 7:9824 ð3Þ
ment of the heating capacity to different heating demands.
In order to decrease the energy consumption of the VRF system
in heating mode, this study developed a control strategy, which 3. Methods
regulates high pressure in the VRF system as a function of the heat-
ing load of indoor units. We also investigated the energy implica- 3.1. Experimental methods
tions of the proposed control model when applied to the VRF
system of an office building. After describing the developed control The effects of changes in target high pressure on the heating
in this study, we elucidate the experimental and simulation meth- capacity and power input of the VRF system were measured in a
ods to examine the energy efficiency and annual heating energy multicalorimeter, which consisted of two indoor unit chambers
consumption with the new control developed in this study. Finally, and one outdoor unit camber. Outdoor unit chamber represented
we present results from a series of experiment and simulations and outdoor conditions while indoor conditions were simulated with
discuss main outcomes. indoor unit chambers. The range of the heat dissipation capacity
1388 G.Y. Yun et al. / Applied Thermal Engineering 113 (2017) 1386–1395
Fig. 1. The process for determining the setpoint for high pressure.
Table 1
Value for a change in room temperature. the capacity of the outdoor unit. The accuracy error of the mul-
ticalorimeter is within ±2% and enables constant temperature
Difference between air and discharge Criterion for air
refrigerant temperatures (DTre, K) temperature change (Trc, K)
and humidity control. The temperature range of the outdoor unit
chamber was 30 to 60 °C, while the temperature range of 0 to
DTre 6 18 4
18 < DTre 6 23 5
50 °C was simulated in the indoor unit chamber. Thus, the mul-
23 < DTre 6 27 6 ticalorimeter used in this study was adequate for representing var-
27 < DTre 6 32 7 ious indoor and outdoor conditions.
DTre > 32 8 A code tester made in accordance with the calorimeter air
enthalpy method [30] was used to determine the air side heating
capacity of the indoor unit chamber. Air temperature, humidity,
Table 2 and airflow rates of the indoor and outdoor unit chambers were
Target high pressures as a function of capacity-weighted average heating load
recorded by the monitoring system. Dry-bulb and wet-bulb tem-
coefficient.
peratures were recorded at the inlet and outlet of each indoor unit
Level Capacity-weighted average Target high chamber. The monitoring system also measured and saved the
heating load coefficient pressure (kgf/cm2)
temperature and pressure of the refrigerant and the heating capac-
1 Ld = 1.0 30 ity and electricity input. The ratios of EEV openings were also
2 0.9 6 Ld 6 1.0 29
recorded.
3 0.8 6 Ld < 0.9 28
4 0.7 6 Ld < 0.8 27 The outdoor unit chamber was air-conditioned by one 24 HP
5 0.6 6 Ld < 0.7 26 outdoor unit (AM240FXVGHH1) while each indoor unit chamber
6 Ld < 0.6 25 enclosed three 4 HP indoor units (ND11104HXB1). The total piping
length from the outdoor unit to the branch joint was 35 m and
each indoor unit was 15 m from the branch joint. Two compressors
with a power of 12.78 kW each were used in the outdoor unit. The
of the indoor unit chamber was from 5.58 to 55.82 kW, while that power of each fan was 400 W. R410A (14.8 kg) was used as the
of the outdoor unit chamber was from 10.47 to 104.67 kW. The refrigerant. The minimum outdoor temperature that the outdoor
heating capacity of the two indoor unit chambers was equal to unit could operate at was 20 °C and heat dissipation capacity
G.Y. Yun et al. / Applied Thermal Engineering 113 (2017) 1386–1395 1389
was 61 kW at the outdoor temperature of 15 °C. Each indoor unit to test the significance of individual regression coefficients. There
had a heating capacity of 12.8 kW and the power consumption of was a clear difference in the VRF system performance in low and
each unit was 78 W. The specification of the indoor unit is given high outdoor temperature regions, and thus separate curves for
in Table 3. In the experiment, the refrigerant condensing tempera- those two distinct areas were required. The heating capacity ratio
ture was controlled by the target high pressure. The range of the boundary curve was used to separate the low and high tempera-
target high pressure was 25–30 kgf/cm2 at 1 kgf/cm2 intervals. ture areas with different performances. The cubic equation form
The target high pressure was changed by regulating the compres- was used to generate performance curve coefficients of the heating
sor speed. The heating capacity and power consumption of the sys- capacity ratio boundary curve in this study (Eq. (4)):
tem were first measured at the reference high pressure of 27 kgf/
cm2. Indoor dry-bulb temperature ranged from 16 to 24 at 2 °C T OA;WB ¼ a þ bðT I;DB Þ þ cðT I;DB Þ2 þ dðT I;DB Þ3 ð4Þ
intervals while outdoor dry-bulb temperature ranged from Here, TOA,WB is the outdoor wet-bulb temperature, TI,DB is the indoor
20 °C to 15 °C, as shown in Table 4. Heating capacity and electric- dry-bulb temperature, and a, b, c, and are performance curve coef-
ity consumption at other target high pressures were measured ficients. Heating capacity modifier performance curves in the VRF
under standard conditions (indoor dry-bulb temperature of 20 °C, model identified changes in heating capacity as the indoor and out-
outdoor wet-bulb temperature of 6 °C, and outdoor dry-bulb tem- door temperatures were varied. Heating capacity correction factor
perature of 7 °C). (CAPFTheating), which is used to create performance curve coeffi-
cients, is dependent on the dry-bulb temperature of the air entering
the heating coil (Tdb) and the wet-bulb temperature of the air enter-
3.2. Simulation
ing a condenser (Tc). CAPFTheating is the ratio of the actual heating
capacity to the rated heating capacity. The bi-quadratic function
A dynamic building energy simulation program, EnergyPlus
show in Eq. (5) was used.
[31], was used in this study. EnergyPlus calculates thermal loads
of a building and the energy demands of heating, ventilating, and CAPFT heating ¼ a þ bðT db Þ þ cðT db Þ2 þ dðT c Þ þ eðT c Þ2 þ f ðT db ÞðT db Þ
air-conditioning (HVAC) systems using the heat balance method
ð5Þ
[32]. A key feature of EnergyPlus is its rigorous validation, which
guarantees the accuracy and reliability of the simulation outcomes Here, a, b, c, d, e, and f are performance curve coefficients. The heat-
[32]. The VRF heat pump air conditioner model in EnergyPlus [33] ing energy input ratio modifier describes changes in heating power
has been validated against experimental data [34]. input as a function of indoor and outdoor temperatures. The bi-
Fig. 2 shows the conceptual diagram of the VRF system mod- quadratic regression equation for the performance curve of the
elled in this study. A system model consisted of six energy recovery heating energy input ratio modifier consists of the heating energy
ventilators (ERV) linked to a VRF system. A 24 HP outdoor unit input ratio correction factor (EIRFTheating), the dry-bulb temperature
(AM240FXVGHH1) provided heating to a typical office floor. Both of the air entering the heating coil (Tdb), and the wet-bulb temper-
full-load and part-load performance data available from manufac- ature of air entering the condenser (Tc) as shown in Eq. (6).
turers are required as inputs for a VRF model in EnergyPlus. The
full-load performance data identifies variations in power input EIRFT heating ¼ a þ bðT db Þ þ cðT db Þ2 þ dðT c Þ þ eðT c Þ2 þ f ðT db ÞðT db Þ
and heating capacity as a function of outdoor and indoor condi- ð6Þ
tions, while part-load performance data reflects capacity and
power input characteristics as the compressor speed varies. Perfor- To derive EIRFT, the ratio of the actual power input to the rated
mance curves to simulate variations in the power input and capac- power input was first calculated and then divided by CAPFT. Heat-
ity with changes in both indoor and outdoor conditions and ing energy input ratio modifier function of the part-load ratio
compressor speeds are key inputs for VRF models in EnergyPlus. defines the power input of the VRF system as a function of part-
The setpoint of high pressure in the developed control varies load ratio (PLR). In this study, we used the cubic function of the
with the heating load and the performance of the VRF system at part-load ratio (Eq. (7)):
full-load and part-load conditions changes with the varied high EIRFPLROA;WB ¼ a þ bðPLRÞ þ cðPLRÞ2 þ dðPLRÞ3 ð7Þ
pressure. Thus, this study generated a series of performance curves
at the target high pressures from 25 kgf/cm2 to 30 kgf/cm2 to where EIRFPLR is the ratio of operating power to the rated power of
model the VRF system in EnergyPlus. The heating capacity and the VRF system. EIRFPLR was calculated according to the combina-
electricity consumption data measured at the indoor and outdoor tion ratio, which is the ratio of the total capacity of the indoor units
conditions specified in Table 4 were used to create the five sets to the capacity of the outdoor unit.
of performance curves at the setpoints of high pressures from Heating performance changes of the VRF system due to varia-
25 kgf/cm2 to 30 kgf/cm2 at 1 kgf/cm2 intervals. Raustad [35] tion in target high pressure were modelled with a custom com-
developed a method to create the performance curves of the out- puter code to implement the developed control. A programming
door unit for the modelling of the VRF system in EnergyPlus using language called EnergyPlus runtime language (Erl) was used for
experimental data and this study used the method by Raustad [35]. the development of the custom computer code [36]. The custom
Regression analysis was used to generate performance curve computer code first determined the capacity-weighted average
coefficients. The F-test was used to evaluate the statistical signifi- heating load coefficient (Ld) in load responsive control of high
cance of the overall regression model while the T-test was applied pressure at each simulation time step. The custom code obtained
Table 3
Specification of the indoor unit.
Table 4
Indoor and outdoor conditions for the experiment.
Fig. 2. Conceptual diagram of the VRF system in the case study building.
changes in zone air temperatures, setpoint temperature, heating Korea [38]. The thermal transmittances of the external walls and
coil runtime, and operation status of an indoor unit in order to cal- the windows were 0.270 W/(m2 K) and 2.100 W/(m2 K), respec-
culate Ld. The code then evaluated the level of the indoor heating tively. Meteonorm [39] was used to get the weather data of Suwon,
loads. Eventually, the code replaced the heating performance South Korea for EnergyPlus simulations. This study followed the
curves of the VRF system in accordance with the level of the indoor specification of the US DOE medium model for modelling both
heating loads. internal heat gains from people, lighting, equipment, and sched-
We selected the prototype of the medium office model devel- ules, including occupancy and HVAC operations.
oped by the US Department of Energy (US DOE) [37], to evaluate
the energy saving potentials of thermal load based control of the
VRF system. Thee medium office model was divided into six 4. Results
perimeters and two internal zones with a floor area of 1661 m2.
The external walls and the windows was defined according to Fig. 3 illustrates the effect of the target high pressure on the
the specification of the energy saving standard for buildings in electricity consumption and heating capacity of the VRF system
G.Y. Yun et al. / Applied Thermal Engineering 113 (2017) 1386–1395 1391
Fig. 4. Performance of the VRF system as a function of high pressure at different outdoor conditions: (a) heating capacity, (b) electricity consumption, (c) energy efficiency.
G.Y. Yun et al. / Applied Thermal Engineering 113 (2017) 1386–1395 1393
Fig. 5. Heating electricity consumption with and without the target high pressure control.
SD = 4.9 °C) in January to 6.6 °C (SD = 5.2 °C) in February (Fig. 6). energy delivered to the building on January 3 was decreased by
The mean monthly outdoor temperature during November to 25% when the developed control was used. The heat delivered
February was less than 10 °C. This confirms the crucial impact of was simulated with the code in the simulation model, which chan-
outdoor conditions on the heating energy use of a building. ged the sets of the performance curves according the level of
The use of the heat load based control of the VRF system indoor heating demands. The performance curves included the
reduced the annual heating electricity consumption by 22%. The heating capacity correction factor (CAPFTheating) and the heating
annual heating electricity consumption reduced from 5325 kW h energy input ratio corrector factor (EIRFTheating) as explained in
to 4160 kW h as the developed control was used to operate the Section 3.2. The heating energy saving predicted by the simulation
VRF system. The energy savings by the developed control are sig- model was smaller than a value from the experiment results. This
nificant compared to the advanced control [40] that shows the was because the setpoint of high pressure was not always lower
energy reduction by through the employment of artificial neural than 30 kgf/cm2. There was a variation in the rates of reduction
network models. The heat delivered to the building also decreased in monthly heating electricity consumption when the developed
when the developed control was applied. For example, the heating control was applied. The use of the developed control reduced
1394 G.Y. Yun et al. / Applied Thermal Engineering 113 (2017) 1386–1395
[31] US Department of Energy, Building Technologies Office: EnergyPlus Energy [36] US DOE, EnergyPlus Documentation: Application Guide for EMS, US
Simulation Software, <http://apps1.eere.energy.gov/buildings/energyplus/> Department of Energy, 2015.
(accessed June 18, 2016). [37] Office of Energy Efficiency & Renewable Energy, Commercial Reference
[32] R.H. Henninger, M.J. Witte, EnergyPlus testing with global energy balance test, Buildings|Department of Energy, <http://energy.gov/eere/buildings/
Gard Analytics, 2013. commercial-reference-buildings> (accessed June 18, 2016).
[33] R. Raustad, Computer modeling VRF heat pumps in commercial buildings [38] Korean Energy Agency, Energy saving standard for buildings, 2015, Ministry of
using EnergyPlus, ASHRAE Trans. 119 (2013) 10. Land, Infrastructure and Transport.
[34] Y.P. Zhou, J.Y. Wu, R.Z. Wang, S. Shiochi, Y.M. Li, Simulation and experimental [39] Meteotest. Meteonorm – Global Meteorological Database, Meteotest, 2015.
validation of the variable-refrigerant-volume (VRV) air-conditioning system in [40] J.W. Moon, S.K. Jung, Development of a thermal control algorithm using
EnergyPlus, Energy Build. 40 (6) (2008) 1041–1047, http://dx.doi.org/10.1016/ artificial neural network models for improved thermal comfort and energy
j.enbuild.2007.04.025. efficiency in accommodation buildings, Appl. Therm. Eng. 103 (2016) 1135–
[35] R. Raustad, Verification of a VRF heat pump computer model in EnergyPlus, 1144, http://dx.doi.org/10.1016/j.applthermaleng.2016.05.002.
ASHRAE Trans. 119 (2013) 101.