Engineering Reference — EnergyPlus 9.3

<< Prev | Table of Contents | Next >>

Engineering Reference

Variable Refrigerant Flow Heat Pumps[LINK]

A Variable Refrigerant Flow (VRF, or Variable Refrigerant Volume) system is an air-conditioning system that varies the refrigerant flow rate using variable speed compressor(s) in the outdoor unit, and the electronic expansion valves (EEVs) located in each indoor unit. The system meets the space cooling or heating load requirements by maintaining the zone air temperature at the setpoint. The ability to control the refrigerant mass flow rate according to the cooling and/or heating load enables the use of as many as 60 or more indoor units with differing capacities in conjunction with one single outdoor unit. This unlocks the possibility of having individualized comfort control, simultaneous heating and cooling in different zones, and heat recovery from one zone to another. It may also lead to more efficient operations during part-load conditions. VRF can be divided into two common types:

  • Heat Pump (HP) type: the most general type that can be used for either cooling or heating, but not simultaneously.

  • Heat Recovery (HR) type: can deliver simultaneous heating and cooling to different zones by transferring heat between the cooling and heating indoor units. This generally occurs in the winter season in medium-sized to large-sized commercial buildings with a substantial core such as computer rooms.

There are two alternative VRF models available in EnergyPlus to simulate the energy performance of VRF:

  1. System Curve based Model (VRF-SysCurve). In this model, a number of system level curves are used to describe the VRF system performance. It uses the AirConditioner:VariableRefrigerantFlow object to describe both the HP and HR systems.

  2. Physics based Model (VRF-FluidTCtrl). This model is able to consider the dynamics of more operational parameters and is applicable for fluid temperature control. It uses the AirConditioner:VariableRefrigerantFlow:FluidTemperatureControl:HP object for the HP system and emphAirConditioner:VariableRefrigerantFlow:FluidTemperatureControl:HR for the HR system.

System Curve based VRF Model (VRF-SysCurve)[LINK]

Overview[LINK]

This section describes the System Curve based VRF Model (VRF-SysCurve), which can address both the heat pump and heat recovery operating modes and support air-, evaporatively-, or water-cooled condenser equipment ( throughout this section, the term “condenser” refers to the outdoor unit where the compressor is located. )

The figure below schematically depicts the AirConditioner:VariableRefrigerantFlow (VRF AC) system. The outdoor unit is connected directly to the zone terminal units using a zone terminal unit list (ref: ZoneTerminalUnitList). The VRF AC system conditions multiple zones and is controlled by thermostats located in each zone. Zone terminal units operate to meet the zone sensible cooling or sensible heating requirements as determined by the zone thermostat schedule.

When the heat pump does not operate to reclaim waste heat, the VRF AC system can only operate in either cooling or heating mode. Based on the master thermostat priority control selection, the operating mode is determined by polling the appropriate zone(s) served by the VRF HP system. When the system is operating in cooling mode, the cooling coils will be enabled only in the terminal units where zone cooling is required. When the system is operating in heating mode, the heating coils will be enabled only in the terminal units where zone heating is required. Supply air fans will continue to operate if the zone terminal unit’s fan operating mode is set to continuous fan.

When the heat pump does operate to reclaim waste heat, the VRF AC system can simultaneously cool and heat multiple zones. The heat pump will select an operating mode according to the dominant load as reported by the zone thermostat(s). The calculation of the dominant load is based on the master thermostat priority control selection and may either be based on individual zone loads, the number of zones requiring cooling or heating, the master thermostat zone load, or an operating mode schedule. The heat pump will operate in cooling mode, and provide waste heat to zones with a heating load, when the dominant load among all zone terminal units is cooling. The heat pump will operate in heating mode, and absorb heat from zones with a cooling load, when the dominant load among all zone terminal units is heating.

The figure below shows the VRF AC terminal units with draw through fan placement. Blow through fan placement can also be modeled by connecting the supply air fan inlet node to the outside air mixer’s mixed air node if an outdoor air mixer is used or to the zone terminal unit inlet node if an outdoor air mixer is not used. The variable refrigerant flow heat pump coordinates the operation of these components and is modeled as a type of zone equipment where the zone terminal units are specified in a zone equipment list (Ref. ZoneHVAC: EquipmentList and ZoneHVAC:EquipmentConnections). The AirConditioner: VariableRefrigerantFlow object is not specified in an air primary loop or a zone equipment list object.

Variable Refrigerant Flow Heat Pump (draw through fan placement) [fig:variable-refrigerant-flow-heat-pump-draw]

Variable Refrigerant Flow Heat Pump (draw through fan placement) [fig:variable-refrigerant-flow-heat-pump-draw]

Energyplus object connections for VRF-SysCurve model [fig:energyplus-object-connections-for-vrf]

Energyplus object connections for VRF-SysCurve model [fig:energyplus-object-connections-for-vrf]

The terminal unit is able to model supply air fan operation in two modes: cycling fan – cycling coil (i.e., AUTO fan) and continuous fan – cycling coil (i.e., fan ON). Fan:Simple:OnOff must be used to model AUTO fan, while Fan:Simple:OnOff or Fan:Simple:ConstVolume can be used to model fan ON.

Output variables reported by the VRF AC object include the heat pump’s operating capacity (cooling or heating), electric consumption, operating COP, defrost electric consumption, part-load ratio (ratio of actual capacity to available capacity), runtime fraction (when cycling), cycling ratio (actual compressor cycling rate), crankcase heater power, and other report variables associated with an evaporative condenser. Report variables are also available to indicate the maximum available terminal unit cooling or heating capacity (i.e., when the condensers available capacity is insufficient to supply the capacity requested by all terminal units). Reporting of other variables of interest for the VRF AC (fan electric power, DX coil cooling rate, DX coil heating rate, terminal unit total cooling rate, etc.) is done by the individual system components (e.g., fan, DX cooling coil, DX heating coil, and zone terminal unit). For convenience, the total cooling and heating coil capacities (i.e., the coil loads) is also reported.

Model Description[LINK]

As described previously, the VRF AC system conditions multiple zones and is controlled by one or more zone thermostats (ZoneControl:Thermostatic). Each simulation time step, EnergyPlus performs a zone air heat balance to determine if cooling or heating is required to meet each zone’s thermostat set point. When heat pump mode is selected (i.e., heat recovery is not selected), an operating mode is chosen based on the master thermostat priority control and all terminal units are operated in that specific mode. When heat recovery mode is selected, the dominant load is determined based on the master thermostat priority control and the outdoor unit operates in this mode. The indoor terminal units operate in either cooling or heating mode based on the individual coil loads. The heat pump condenser is then modeled to determine any impact the condensing unit might have on the zone terminal units (i.e., capacity limitations due to over-sized terminal units, operating limitations due to allowable operating temperature range, impacts of defrost operation, changes in performance when heat recovery mode is active, etc.). The following sections describe the performance calculations for cooling and heating.

Cooling Operation[LINK]

The operating capacity of the heat pump is calculated based on the user defined input for the heat pump’s rated cooling capacity, the ratio of indoor terminal unit capacity to outdoor unit capacity (combination ratio), the actual operating conditions, and whether or not the system’s heat recovery mode is active. Performance correction factors are used to correct for off-design performance as follows:

Cooling Combination Ratio Correction Factor[LINK]

The cooling combination ratio correction factor is defined as the total terminal unit rated cooling capacity divided by the heat pump’s rated cooling capacity. The combination ratio is important when multiple terminal units (or DX coils) are attached to a single variable-speed condensing unit. If the combination ratio is less than 1, the condenser’s rated capacity is assumed to be able to meet the indoor terminal unit’s demand request. However, if the combination ratio is determined to be greater than 1, the maximum cooling capacity available from the heat pump’s condenser may be higher than the user specified rated cooling capacity. The combination ratio capacity correction factor is based on a linear, quadratic or cubic curve and provides a multiplier > 1 to correct for combination ratio’s greater than 1. For example, a combination ratio of 1.3 (130%) may result in a combination ratio capacity correction factor of 1.06 (outdoor condenser can actually provide 106% of rated capacity) . The cooling combination ratio correction factor is applied to cooling performance calculations. If the cooling combination ratio correction curve name is not specified by the user, the cooling combination ratio correction factor (CRcooling,correction) in the following equation is assumed to be equal to 1.

CRcooling,rated=n1˙Qcoil(i),cooling,rated˙Qcooling,total,rated

CRcooling,correction=a+b(CRcooling,rated)+c(CRcooling,rated)2+d(CRcooling,rated)3

CRcooling,correction=MAX(1.0,CRcooling,correction)

where:

˙Qcoil(i),cooling,rated is the rated total (sensible + latent) cooling capacity in zone i (W)

˙Qcooling,total,rated is the rated total cooling capacity of heat pump (W)

CRcooling,correction is the Cooling Combination Ratio capacity correction factor at rated conditions (this value is reported in the eio file)

ad are the equation coefficients for cooling combination ratio correction factor

CRcooling,rated is the cooling combination ratio defined as the total indoor terminal unit’s rated total cooling capacity divided by the rated total cooling capacity of the heat pump condenser.

The sum of the individual zone total cooling requirements is used to calculate the performance of the outdoor condensing unit. The operating capacities of the indoor cooling coils are calculated based on the indoor cooling coil’s rated cooling capacity and the actual operating conditions. The operating capacity of the heat pump condenser is calculated in a similar fashion using a load-weighted average indoor wet-bulb temperature of all operating cooling coils.

Cooling Capacity Ratio Modifier Function of Low Temperature[LINK]

The heat pump’s cooling capacity correction factor (function of temperature) is determined by averaging the zone cooling coil inlet air wet-bulb temperature for all operating cooling coils. A load-weighted average inlet air wet-bulb temperature is used in the calculations. The weighted average cooling coil inlet air wet-bulb temperature and the outdoor condenser entering air dry-bulb temperature are then used to calculate the temperature correction factor in cooling mode for the heat pump condenser. If the Condenser Type is WaterCooled, then the cooling capacity curve will be function of weighted average cooling coil inlet air wet-bulb temperature and outdoor condenser entering water temperature. The actual zone air wet-bulb temperature is used to calculate the operating capacity of the individual zone cooling coils.

˙Qzone,total=i1˙Qzone(i)

Twb,avg=i1(Twb,i)˙Qzone(i)˙Qzone,total

The operating capacity of the heat pump is calculated using a bi-quadratic equation using a load-weighted average indoor wet-bulb temperature and outdoor dry-bulb temperature as the independent variables.

CAPFTHP,cooling=a+b(Twb,avg)+c(Twb,avg)2+d(Tc)+e(Tc)2+f(Twb,avg)(Tc)

The operating capacity of the terminal unit’s DX cooling coil is calculated using either a linear, quadratic, or cubic curve. If more information is available to more accurately model the DX cooling coil, a bi-quadratic curve may be used (i.e., if the performance of the DX cooling coil is also a function of outdoor dry-bulb temperature).

CAPFTcoil,cooling=a+b(Twb,i)+c(Twb,i)2+d(Twb,i)3

or

CAPFTcoil,cooling=a+b(Twb,i)+c(Twb,i)2+d(Tc)+e(Tc)2+f(Twb,i)(Tc)

where:

Twb,i is the wet-bulb temperature of the air entering the cooling coil in zone i (C)

Twb,avg is the load-weighted average wet-bulb temperature of the air entering all operating cooling coils (C)

˙Qzone(i) is the total (sensible + latent) cooling load in zone i (W)

˙Qzone,total is the total (sensible + latent) cooling load in all zones (W)

CAPFTHP,cooling is the heat pump Cooling Capacity Ratio Modifier (function of temperature)

CAPFTcoil,cooling is the zone coil Cooling Capacity Ratio Modifier (function of temperature)

af are the equation coefficients for Cooling Capacity Ratio Modifier

Tc is the temperature of the air entering an air-cooled or evaporatively-cooled condenser (C).

Using multiple curves to define Cooling Capacity Ratio Modifier[LINK]

The cooling capacity ratio modifier determines the change in total (sensible + latent) capacity with respect to the heat pump rated cooling capacity. This modifier corrects for off-design performance and provides the operating total (sensible + latent) cooling capacity for the heat pump condenser. The performance data for VRF AC systems may be specified using a single curve object, however, if the performance data does not reflect a smooth change in performance as outdoor conditions vary, the overall performance of the system may be described using two cooling capacity ratio modifier curves (i.e., a separate curve at low and high outdoor temperatures) and a boundary curve (i.e., the curve defining the separation at specific outdoor temperatures). If the Condenser Type is WaterCooled, the cooling capacity and EIR modifier curves may be represented using the low temperature curve input field only.

Application of Dual Performance Curves[LINK]

Two additional inputs are available to more accurately model the cooling performance of a VRF AC system. The first additional input is a boundary curve which is used to distinguish differences in performance at “low” and “high” outdoor temperature regions and the second input is a cooling capacity performance curve for the “high” outdoor temperature region. The boundary curve is a linear, quadratic or cubic curve defining the outdoor temperature as a function of average indoor wet-bulb temperature. The cooling capacity ratio function of high temperature curve is a biquadratic curve using average indoor wet-bulb temperature and outdoor dry-bulb temperature as the two independent variables. These additional curves are used to define the performance over two separate performance regions as shown in Figure 4. The red circles in the figure identify the points at which the performance changes abruptly (i.e., slopes or shapes are significantly different) and are used to create the boundary curve. The performance on either side of the boundary curve can be a smooth curve changing with outdoor condition or a flat plateau where system controls maintain a constant operating condition. This type of performance cannot be accurately simulated using a single performance curve object. For this reason, dual capacity performance curves may be used.

Note: If the cooling performance can be described using a single performance curve object, the boundary and high temperature performance curve objects are not required.

Manufacturers may also provide this information in a tabular format. In this case the data should first be graphically displayed to identify any non-linearities and to also identify where an abrupt change in performance occurs so that a boundary curve can be created.

The cooling capacity ratio boundary curve object is used to differentiate between dual cooling capacity ratio performance curves. This curve defines the outdoor temperature below which the cooling capacity ratio function of low temperature curve object is used, otherwise, the cooling capacity ratio function of high temperature curve object is used.

TOA,DB=a+b(TI,WB)+c(TI,WB)2+d(TI,WB)3

where:

TOA,DB is the outside air dry-bulb temperature (C)

a-d are the coefficients for cooling capacity ratio boundary curve object

TI,WB is the weighted average indoor wet-bulb temperature (C).

For the performance data shown below in Figure 3, the boundary curve would be created by identifying the outdoor dry-bulb temperature at which the performance changed abruptly for each of the indoor wet-bulb temperature curves shown in the figure. The following example shows the data interpreted from the figure used to regress the cooling capacity ratio boundary curve coefficients. For this example, the regression coefficients were found to be: A1 = 29.87396, B1 = -0.6928, C1 = 0.01928, D1 = -0.000532.

Performance Data for Variable Refrigerant Flow Air Conditioner Model [fig:performance-data-for-variable-refrigerant-flow-air-conditioner-model]

Performance Data for Variable Refrigerant Flow Air Conditioner Model [fig:performance-data-for-variable-refrigerant-flow-air-conditioner-model]

Although the capacity and energy performance curves each have an independent boundary curve input, in this example the same boundary curve may be used for both the capacity and energy input ratio curves. When the “low” and “high” predicted performance data do not line up at the boundary curve points, the boundary curve may have to be modified slightly to allow a smooth transition in performance from the “low” to “high” temperature region. In this case, the boundary curves for capacity and energy may be different.

Non-Linear Performance of VRF Heat Pump in Cooling Mode [fig:non-linear-performance-of-vrf-heat-pump-in-cooling-mode]

Non-Linear Performance of VRF Heat Pump in Cooling Mode [fig:non-linear-performance-of-vrf-heat-pump-in-cooling-mode]

A regression is then performed on the data to the left of the boundary curve to calculate the low temperature cooling performance curve coefficients. A regression is also performed on the data to the right of the boundary curve to create the high temperature cooling performance curve coefficients. The model then uses the boundary curve to determine which performance curve (low or high) to use during the simulation. For example, given a particular average indoor wet-bulb temperature, if the boundary curve object calculates an outdoor dry-bulb temperature that is above the actual outdoor dry-bulb temperature then the cooling capacity ratio function of low temperature performance curve is used to determine AC system performance for that specific simulation time step. When creating the boundary curve, be careful to make sure the low and high performance curves meet, as closely as possible, at the boundary curve points (i.e., that discontinuities do not exist or are minimized to the extent possible). Tabular data may also be used to specify performance and will usually eliminate the need for dual performance curves.

Determining Cooling Coil Sensible Capacity[LINK]

The zone terminal unit’s cooling coil is controlled by a zone thermostat and must meet a zone “sensible” load. The cooling coil’s sensible capacity is defined as the cooling coil’s total capacity multiplied by the coil’s sensible heat ratio (SHR) at the current operation conditions. Since the SHR of the variable refrigerant flow cooling coil changes as the inlet air wet-bulb temperature and the operating part-load ratio change, an iterative solution technique must be used to determine the coil’s sensible capacity. The target solution is found when the sensible capacity of the zone cooling coil equals the zone sensible cooling load (if sufficient capacity is available). The iterative solution converges when the difference (error) between the zone sensible cooling load and the terminal unit sensible cooling capacity is within tolerance (0.001) or the terminal unit has insufficient capacity to meet the zone sensible cooling load. Refer to the description of the VRF cooling coil model for further details (ref: Variable Refrigerant Flow Cooling Coil).

˙Qcoil(i),cooling,total=˙Qcoil(i),cooling,rated(CAPFTcoil,cooling)

PLRi=1

Iteration Loop (PLRi)    SHRPLR=f(Twb,i,PLRi,˙mi)    ˙Qcoil(i),cooling,sensible=˙Qcoil(i),cooling,total(SHRPLR)(PLRi)    Error=(˙Qzone(i)˙Qcoil(i),cooling,sensible)/˙Qzone(i)End Iteration Loop

where:

˙Qcoil(i),cooling,total is the zone terminal unit total (sensible + latent) cooling capacity (W), report variable “Zone VRF Air Terminal Total Cooling Rate"

˙Qcoil(i),cooling,sensible is the zone terminal unit sensible cooling capacity (W), report variable “Zone VRF Air Terminal Sensible Cooling Rate"

PLRi is the cooling coil sensible part-load ratio in zone i

SHRPLR is the cooling coil sensible heat ratio (function of PLR, inlet air wet-bulb temperature, and cooling coil inlet air mass flow rate)

˙mi is the cooling coil inlet air mass flow rate (m3/s).

The terminal unit total cooling (sensible + latent) capacity is then summed to provide the total terminal unit cooling requirement. The cooling coils capacity includes the impacts of fan heat and any outdoor air provided to the zone.

˙Qcooling,TerminalUnits=i1(˙Qcoil(i),cooling,total)(PLR(i))

where ˙Qcooling,TerminalUnits is the total terminal unit cooling requirement (sensible + latent) in all zones (W).

The piping correction factor is then used to adjust the total zone cooling requirement to account for piping losses in the air conditioner’s refrigeration piping system.

Cooling Piping Correction Factor in cooling mode[LINK]

The cooling piping correction factor is calculated using either one or two independent variables. Both the equivalent piping length and the combination ratio are used together to determine piping losses, or the equivalent piping length itself is used. The vertical height specified by the user is also added to the result calculated by either. If a single independent variable is used, a linear, quadratic, or cubic equation is used to determine the equivalent length correction factor. If two independent variables are used, a biquadratic equation is used to determine the equivalent length correction factor. The program will automatically use the correct equation based on the performance curve type. In either case, a single coefficient adjusts the piping correction factor based on the difference in height from the highest to lowest terminal unit. Coefficients a-d (e) in the following equation(s) are inputs in the piping correction factor for length in cooling mode performance curve object. Coefficient f is a direct input to the VRF model as the piping correction factor for height in cooling mode coefficient. The equivalent piping length in cooling mode (PEQ,cooling) is also a direct input in the VRF model. The cooling combination ratio (CRcooling) is automatically calculated by the program. The vertical height, the difference between the highest and lowest terminal unit (e.g., 12 m higher than condenser – 3 m lower than condenser = 9 m height), is a common input for both cooling and heating. When all terminal units are above or below the outdoor unit, an average positive or negative value is used. The limits on the calculated cooling piping correction factor are 0.5 < Pcorrection,cooling < 1.

CRcooling=n1˙Qcoil(i),cooling,rated˙Qcooling,total,rated

Pcorrection,cooling=a+b(PEQ,cooling)+c(PEQ,cooling)2+d(PEQ,cooling)3+f(PH)

or

Pcorrection,cooling=a+b(PEQ,cooling)+c(PEQ,cooling)2+d(CRcooling)+e(CRcooling)2+f(PEQ,cooling)(CRcooling)+g(PH)

where:

Pcorrection,cooling is the Piping Correction Factor in Cooling Mode

a – f are equation coefficients for piping correction factor in cooling mode

g is the user specified piping correction factor for height in cooling mode coefficient

PEQ,cooling is the user specified equivalent piping length in cooling mode (m)

CRcooling is the combination ratio in cooling mode (total rated indoor terminal unit capacity divided by the rated condenser cooling capacity) (reported to eio file)

PH is the user specified vertical height used for piping correction factor calculation (m).

An example piping correction factor chart is shown in the following figure. The height selected for use in the equation above is selected to minimize the piping correction factor and serves to identify the worst case piping losses.

Typical Piping Factor Correction Chart [fig:typical-piping-factor-correction-chart]

Typical Piping Factor Correction Chart [fig:typical-piping-factor-correction-chart]

The total demand on the heat pump condenser is then calculated as the quotient of the total terminal unit cooling capacity and the cooling piping correction factor. The piping losses are fixed throughout the simulation (i.e., these losses are not based on which terminal units are requesting refrigerant).

˙Qcooling,total=˙Qcooling,TerminalUnitsPcorrection,cooling

The heat pump’s total available cooling capacity is then calculated as:

˙QHP,cooling,total=˙Qcooling,total,rated(CAPFTHP,cooling)(CRcooling,correction)

where:

˙Qcooling,total is the total heat pump condenser cooling load (W)

˙QHP,cooling,total is the heat pump total available cooling capacity (W).

Heat-Recovery-Cooling-Capacity-Modifier[LINK]

When operating in heat recovery mode, the heat pump’s available cooling capacity is typically different than the available capacity when operating in cooling only mode. This modifier is used to adjust the available cooling capacity using a fraction when heat recovery is active. This fraction is based on a bi-quadratic equation with indoor and outdoor temperatures used as the independent terms.

HRCapModHP,cooling=a+b(Tent)+c(Tent)2+d(Tcond)+e(Tcond)2+f(Tent)(Tcond)

This equation can be used to provide a constant fractional difference for available cooling capacity in heat recovery mode (i.e., only a is non-zero) or a fractional term that varies with indoor and outdoor conditions. With very limited performance data available at this time, it is recommended that only the constant (a) term be used at this time. When the VRF system is not operating in heat recovery mode, this fraction is set to 1. The available cooling capacity in heat recovery mode is then:

˙QHR,cooling,total=˙QHP,cooling,total(HRCapModHP,cooling)

where ˙QHR,cooling,total = heat recovery total available cooling capacity (W).

The figure below shows VRF system laboratory data for cooling only mode (solid characters) and heat recovery mode (dotted characters). Using the limited laboratory data, the available cooling capacity fraction used to model heat recovery mode is approximately 0.91 and the cooling energy fraction is approximately 1.14. This is the only data available at this time to estimate the impact of heat recovery mode on performance. In the bi-quadratic equation, only coefficient a should be used until more complete data sets exist. Laboratory testing will eventually provide more data and better estimates of performance in heat recovery mode.

Comparison of cooling only and heat recovery mode operation [fig:comparison-of-cooling-only-and-heat-recovery-mode-operation]

Comparison of cooling only and heat recovery mode operation [fig:comparison-of-cooling-only-and-heat-recovery-mode-operation]

Transition from Cooling Only mode to Heat Recovery mode[LINK]

When the VRF system transitions from cooling only operation to heat recovery operation, this transition takes some finite amount of time. During the transition period the available cooling capacity can change significantly, The following figure illustrates the transition between cooling only mode and heat recovery mode. For this test, the VRF system was turned on and allowed to reach steady-state operation. Three of the four indoor terminal units were operating in cooling mode, When the fourth terminal unit was enabled in heating mode, the transition from cooling only mode to heat recovery mode took approximately 45 minutes. During this time, the available cooling is significantly reduced and recovers over time. When the system again reaches steady-state operation, the available cooling capacity and power consumption are markedly different. Although computer models do not typically simulate this type of transitional performance, efforts to model this aspect of performance were included in the VRF AC heat recovery model. The initial heat recovery cooling capacity fraction and heat recovery cooling capacity time constant are used to model this transition period. The initial heat recovery cooling capacity fraction identifies the fraction of available heat recovery mode cooling capacity at the start of the transition period, the heat recovery cooling capacity time constant identifies the time needed to recover to 99% of the steady-state value. This exponential model used for simulating the transition period can be turned off by setting the initial heat recovery cooling capacity fraction to 1.

Laboratory test identifies performance changes during transition period [fig:laboratory-test-identifies-performance-changes-during-transition-period]

Laboratory test identifies performance changes during transition period [fig:laboratory-test-identifies-performance-changes-during-transition-period]

Heat Recovery Cooling Capacity Fraction and Time Constant[LINK]

When the heat pump changes operating modes (i.e., from cooling only to heat recovery mode), the transition does not happen immediately. There is a time delay, and a period of time where the available cooling capacity is well below the steady-state capacity available after the system has had time to adjust. When this type of transition is modeled, an exponential decay curve is used. At the start of the transition period, only a fraction (kcool) of the steady-state capacity in heat recovery mode is available. The remaining capacity is recovered over a period of 5 time constants (tc) and is modeled with the following equation. This equation was used to replicate both the cooling capacity and condenser power curves in the previous figure.

˙QHR,avail,cooling=kcool(˙QHR,cooling,total)+(1kcool)(˙QHR,cooling,total)(1e(t/τc))

where:

kcool is the fraction of steady-state capacity at beginning of heat recovery mode

tc is the time constant to reach steady-state operation, 5 time constants equals 99% of steady-state capacity (hr)

QHR,avail,cooling is the available cooling capacity in heat recovery mode (W)

(1e(t/tc,cool)) is the multiplier used during the transition period from cooling mode to heat recovery mode. This multiplier can be viewed using the report variable “VRF Heat Pump Heat Recovery Status Change Multiplier” (ref: object Output:Variable).

The operating part-load ratio of the VRF system can then be calculated:

PLR=˙Qcooling,total˙QHR,avail,cooling

where PLR is the heat pump part-load ratio, report variable “VRF Heat Pump Part Load Ratio”.

Note: for calculation purposes ˙QHR,avail,cooling is equivalent to ˙QHR,cooling,total when heat recovery mode is not active.

The heat pump total available cooling capacity must be greater than or equal to the total cooling capacity requested by the zone terminal units. When the total operating capacity of all terminal unit’s will be greater than the available operating capacity of the heat pump condenser, one or more of the terminal unit’s operating capacity must be reduced to the point where the sum of the indoor terminal unit demand request plus piping losses is equal to the total available cooling capacity of the outdoor condenser. At this point, the part-load ratio of the heat pump condenser will be equal to 1.

A maximum terminal unit cooling capacity limit is used to restrict the cooling capacity of each indoor terminal unit. The capacity limit is equivalent to a maximum allowed operating capacity for all zone terminal units. This limit is used to conserve energy between multiple indoor terminal units and a single outdoor condensing unit. Assuming no terminal unit can provide more capacity than can be delivered by the outdoor condenser, the terminal unit total cooling capacity calculation previously discussed is revised as follows:

˙Qcoil(i),cooling=i1MIN((˙Qcoil(i),cooling,total)(PLR(i)),CoolingCapacityLimit)

where CoolingCapacityLimit(CLcooling) is the report variable (ref: Output:Variable) describing the “VRF Heat Pump Maximum Capacity Cooling Rate” (W).

When multiple terminal units are operating, the terminal units near their maximum capacity are more likely to be capacity limited than those terminal units operating well below their available capacity. The assumption here is that terminal units that are not capacity limited can provide more refrigerant to meet the same load. When the model finds that there is no terminal unit capacity limit, this variable will report 1E+20 indicating that no limit exists. The figure below demonstrates the application of the capacity limit factor for the case where all indoor terminal units are operating at different capacities. A solution is reached when the sum of the indoor terminal unit’s cooling capacities (accounting for piping losses) is equal to the heat pump condenser’s available cooling capacity.

Example of Cooling Capacity Limit [fig:example-of-cooling-capacity-limit]

Example of Cooling Capacity Limit [fig:example-of-cooling-capacity-limit]

When the heat pump’s part-load ratio is less than 1 (i.e., the total capacity of all terminal unit’s is less than the actual operating capacity of the heat pump condenser), the heat pump’s part-load ratio is compared to the minimum heat pump part-load ratio. If the heat pump’s part-load ratio is less than the minimum heat pump part-load ratio, the heat pump will cycle on and off to meet the cooling load. A correction factor is used to account for startup losses of the compression system.

Cooling Part-Load Fraction Correlation (function of cycling ratio)

The cooling part-load fraction correlation (function of heat pump cycling ratio) is a linear, quadratic or cubic curve with the independent variable being cycling ratio (part-load ratio / minimum part-load ratio). The output of this curve is used in combination with the rated EIR and EIR modifier curves to give the “effective” EIR for a given simulation time step.

CyclingRatioFrac=a+b(CyclingRatio)+c(CyclingRatio)2+d(CyclingRatio)3

where:

CyclingRatio=PLRPLRmin

a-c (d) are the coefficients for the quadratic (cubic) cycling ratio correlation curve equation.

The cycling ratio can be viewed using the report variable “VRF Heat Pump Cycling Ratio” (ref: Output:Variable).

The cooling part-load ratio correlation should be normalized to a value of 1.0 when the cycling ratio equals 1.0 (i.e., no efficiency losses when the heat pump runs continuously [PLR PLRmin] for the simulation time step). For cycling ratio values between 0 and 1 (0 CyclingRatio < 1), the following rules apply:

CyclingRatioFrac0.7 and CyclingRatioFracCyclingRatio

If CyclingRatioFrac < 0.7 a warning message is issued, the program resets the CyclingRatioFrac value to 0.7, and the simulation proceeds. The runtime fraction of the heat pump is defined as CyclingRatio/CyclingRatioFrac. If CyclingRatioFrac < CyclingRatio, the runtime fraction will exceed 1. In this case a warning message is issued and the runtime fraction of the heat pump is limited to 1.0.

Heat Pump Runtime Fraction (HPRTF)=CyclingRatioCyclingRatioFrac

The heat pump runtime fraction can be viewed using the report variable “VRF Heat Pump Runtime Fraction” (ref: Output:Variable).

Since manufacturers data do not typically provide information defining the cycling losses of VRF AC systems, a typical part-load fraction correlation for a conventional DX refrigeration system (e.g., residential or small commercial unit) may be substituted here as:

CyclingRatioFrac=0.85+0.15(CyclingRatio)

The electrical power consumed by the heat pump condenser is then calculated based on the heat pump’s rated electric power consumption (including condenser fan power), the actual operating conditions, and the cycling ratio of the refrigeration system. Performance correction factors are used here to correct for off-design performance as follows:

Cooling Energy Input Ratio Modifier Function of Low Temperature[LINK]

As described previously (Ref. Application of Dual Performance Curves), the cooling energy input ratio modifier curve can either be applied as a single curve object as shown below, or applied using dual performance curve objects. If the Condenser Type is WaterCooled, then the cooling energy input ratio modifier curve will be function of weighted average indoor air wet-bulb temperature and outdoor condenser entering water temperature. The equation form used for a single or dual performance curve application is shown here.

EIRFTcooling=a+b(Twb,avg)+c(Twb,avg)2+d(Tc)+e(Tc)2+f(Twb,avg)(Tc)

where EIRFTcooling is the cooling energy input ratio modifier (function of temperature).

Cooling Energy Input Ratio Modifier Function of Part-Load Ratio[LINK]

The cooling energy input ratio calculation is based on 2 EIR curves. One is used when the heat pump part-load ratio is less than or equal to 1, and a second curve is used to describe the variation of cooling energy input ratio when the part-load ratio is greater than 1. The part-load ratio curve when PLR>1 is not a required input and when not used, the energy use is assumed to be constant when PLR > 1.

IF(PLR1)THEN    EIRFPLRcooling=a+b(PLR)+c(PLR)2+d(PLR)3ELSE    EIRFPLRcooling=e+f(PLR)+g(PLR)2+h(PLR)3

where:

EIRFPLRcooling is the cooling energy input ratio correction factor (function of part-load ratio)

a-d are the coefficient for cooling energy input ratio correlation when part-load ratio 1

e-h are the coefficient for cooling energy input ratio correlation when part-load ratio > 1.

The total power consumed by the heat pump condenser in cooling mode is based on the user specified coefficient of performance (COP) and calculated as:

CoolingPower=˙Qcooling,total,ratedCAPFTcoolingCOPcooling,reference(EIRFTcooling)(EIRFPLRcooling)(HPRTF)

The cooling COP is then calculated as the ratio of the outdoor unit cooling capacity divided by the total electric consumption rate of all components associated with the outdoor unit. Although included in the equiation, defrost power is typically 0 during cooling mode operation.

COPcooling=˙Qcooling,totalCoolingPower+PCrankCaseHeater+PEvapCoolerPower+Pdefrost

where:

COPcooling is the operating coefficient of performance, report variable “VRF Heat Pump Cooling COP”

CoolingPower is the operating electric consumption rate (W), report variable “VRF Heat Pump Cooling Electric Power”

PCrankcaseHeater is the report variable for electric consumption rate of crankcase heater (W)

PEvapCoolerPump is the report variable for electric consumption rate of evaporatively-cooled condenser water pump (W)

Pdefrost is the report variable for average defrost power for the simulation time step (W).

When operating in heat recovery mode, equations similar to those used for available cooling capacity are used to model heating electric consumption rate.

HREIRModHP,cooling=a+b(Tent)+c(Tent)2+d(Tcond)+e(Tcond)2+f(Tent)(Tcond)

This equation can be used to provide a constant fractional difference for cooling electric consumption rate in heat recovery mode (i.e., only a is non-zero) or a fractional term that varies with indoor and outdoor conditions. With very limited performance data available at this time, it is recommended that only the constant (a) term be used at this time. When the VRF system is not operating in heat recovery mode, this fraction is set to 1. The cooling electric consumption rate in heat recovery mode is then:

CoolingPowerHR=CoolingPower(HREIRModHP,cooling)

PowerHR=kEIR,cool(CoolingPowerHR)+(1kEIR,cool)(CoolingPowerHR)(1e(t/τc,EIR,cool))

where:

kEIR,cool is the fraction of steady-state cooling electric consumption rate at beginning of heat recovery mode

τc,EIR,cool is the time constant to reach steady-state operation, 5 time constants equals 99% of steady-state capacity (hr)

CoolingPowerHR is the cooling electric consumption rate in heat recovery mode (W)

(1e(t/τc,EIR,cool)) is the multiplier used during the transition period from cooling mode to heat recovery mode. This multiplier can be viewed using the report variable “VRF Heat Pump Heat Recovery Status Change Multiplier” (ref: Output:Variable).

Heating Operation[LINK]

Calculations of the heat pump’s heating performance is nearly identical to the calculations described above for cooling operation.

The sum of the individual zone total heating requirements are used to calculate the performance of the outdoor condensing unit. The operating capacity of the indoor heating coils are calculated based on the indoor heating coil’s nominal heating capacity and the actual operating conditions. A capacity correction factor is used to correct for off-design performance as follows:

Heating Combination Ratio Correction Factor (function of capacity ratio)

The heating combination ratio correction factor is defined as the total terminal unit rated heating capacity divided by the heat pump’s rated heating capacity. In some instances, the cooling combination ratio may be used to define the correction factor used to modify the heat pump’s rated heating capacity and is manufacturer specific. The combination ratio is important when multiple terminal units (or DX coils) are attached to a single variable-speed condensing unit. If the combination ratio is less than 1, the condenser has sufficient capacity to meet the indoor terminal unit’s capacity request. However, if the combination ratio is determined to be greater than 1, the maximum heating capacity available from the heat pump’s condenser may be higher than the user specified rated heating capacity. The combination ratio capacity correction factor is based on a linear, quadratic or cubic curve and provides a multiplier > 1 to correct for combination ratio’s greater than 1. For example, a combination ratio of 1.3 (130%) may result in a combination ratio capacity correction factor of 1.06 (outdoor condenser can actually provide 106% of rated capacity) . The heating combination ratio correction factor is applied to heating performance calculations. If the heating combination ratio correction curve name is not specified by the user, the heating combination ratio correction factor (CRheating,correction) in the following equation is assumed to be equal to 1.

CRheating=n1˙Qcoil(i),heating,rated˙Qheating,total,rated

CRheating,correction=a+b(CRheating)+c(CRheating)2+d(CRheating)3

CRheating,correction=MAX(1.0,CRheating,correction)

where:

˙Qcoil(i),heating,rated is the rated total heating capacity in zone i (W)

˙Qheating,total,rated is the rated total heating capacity of heat pump (W)

CRheating,correction is the heating Combination Ratio capacity correction factor at rated conditions

ad are the cubic equation coefficients for heating combination ratio correction factor

CRheating is the combination ratio in heating mode (total rated indoor terminal unit capacity divided by the rated condenser heating capacity) (reported to eio file).

Heating Capacity Ratio Modifier (function of temperature)[LINK]

As described previously (Ref. Application of Dual Performance Curves), the heating capacity ratio modifier curve can either be applied as a single curve object as shown below, or applied using dual performance curve objects. The single curve object application is described here.

The heat pump’s heating capacity ratio modifier (function of temperature) is determined by averaging the zone heating coil inlet air dry-bulb temperature for all operating heating coils. A zone load-weighted average is used when calculating the performance of the heat pump condenser. This weighted-average heating coil inlet air dry-bulb temperature and the outdoor condenser entering air wet-bulb temperature are then used to calculate the heat pump’s temperature correction factor in heating mode.

Note that some manufacturers do not provide performance data as a function of outdoor wet-bulb temperature. In this case, substituting outdoor dry-bulb temperature is permitted and the Heating Performance Curve Outdoor Temperature Type input should be specified as DryBulbTemperature. This also means that performance curve coefficients for both capacity (CAPFT) and energy (EIRFT) should be calculated using indoor and outdoor dry-bulb temperature.

Tdb,avg=i1(Tdb,i)˙Qzone(i)˙Qzone,total

CAPFTHP,heating=a+b(Tdb,avg)+c(Tdb,avg)2+d(Tc)+e(Tc)2+f(Tdb,avg)(Tc)

CAPFTcoil,heating=a+b(Tdb,i)+c(Tdb,i)2+d(Tc)+e(Tc)2+f(Tdb,i)(Tc)

where:

Tdb,i is the dry-bulb temperature of the air entering the heating coil in zone i (C

Tdb,avg is the weighted-average dry-bulb temperature of the air entering all operating heating coils (C

CAPFTheating is the heating Capacity Correction Factor (function of temperature)

af are the bi-quadratic equation coefficients

Tc is the wet-bulb temperature of the air entering an air-cooled condenser (C.

The total terminal unit heating capacity required is the simple sum of the terminal unit capacity. If the heat pump is off or there is no zone heating requirement, the terminal unit total heating requirement will be zero. The heating coils capacity includes the impacts of fan heat and any outdoor air provided to the zone.

˙Qheating,TerminalUnits=i1(˙Qcoil(i),heating,total)(PLR(i))

where:

˙Qheating,TerminalUnits is the heat pump’s outdoor condenser total zone heating requirement (W)

˙Qcoil(i),heating,total is the zone terminal unit total heating capacity (W), report variable “Zone VRF Air Terminal Total Heating Rate“.

The piping correction factor is then used to adjust the zone heating requirement to account for piping losses in the AC refrigeration system.

Piping Correction Factor in heating mode[LINK]

The piping correction factor in heating mode is calculated based on the length of the farthest terminal unit and the difference in height from the highest to lowest terminal unit. The piping losses are fixed throughout the simulation (i.e., these losses are not based on which terminal units are requesting refrigerant). Coefficients a-f in the following equations are inputs in the piping correction factor for length in heating mode curve object. The curve may use either one or two independent variables. Coefficient g is a direct input to the VRF model as piping correction factor for height in heating mode coefficient. The equivalent piping length in heating mode (PEQ,heating) is also a direct input in the VRF model. The vertical height (PH), the difference between the highest and lowest terminal unit (e.g., 12 m higher than condenser – 3 m lower than condenser = 9 m height), is a common input for both cooling and heating.

Pcorrection,heating=a+b(PEQ,heating)+c(PEQ,heating)2+d(PEQ,heating)3+g(PH)

or

Pcorrection,heating=a+b(PEQ,heating)+c(PEQ,heating)2+d(CRheating)+e(CRheating)2+f(PEQ,heating)(CRheating)+g(PH)

where:

Pcorrection,heating is the piping correction factor in heating mode

a-f are the equation coefficients for piping correction factor in heating mode

g is the user specified piping correction factor for height in heating mode coefficient

PEQ,heating is the equivalent piping length for heating specified by the user (m)

PH is the user specified vertical height used for piping correction factor (m).

The heat pump’s total available heating capacity is then determined using the previously described modifiers.

˙QHP,cooling,total=˙Qcooling,total,rated(CAPFTHP,cooling)(CRcooling,correction)

Heat Recovery Heating Capacity Modifier[LINK]

When operating in heat recovery mode, the heat pump’s available heating capacity is typically different than the available capacity when operating in heating only mode. This modifier is used to adjust the available heating capacity using a fraction when heat recovery is active. This fraction is based on a bi-quadratic equation with indoor and outdoor temperatures used as the independent terms.

HRCapModHP,heating=a+b(Tent)+c(Tent)2+d(Tcond)+e(Tcond)2+f(Tent)(Tcond)

This equation can be used to provide a constant fractional difference for available heating capacity in heat recovery mode (i.e., only a is non-zero) or a fractional term that varies with indoor and outdoor conditions. With very limited performance data available at this time, it is recommended that only the constant (a) term be used at this time. The available heating capacity in heat recovery mode is then:

˙QHR,heating,total=˙QHP,heating,total(HRCapModHP,heating)

Transition from Heating Only mode to Heat Recovery mode[LINK]

When the VRF system transitions from heating only operation to heat recovery operation, this transition takes some finite amount of time. During the transition period the available heating capacity can change significantly, As described for the transition from cooling only operating to heat recovery mode, the initial heat recovery heating capacity fraction and heat recovery heating capacity time constant are used to model the transition from heating only mode to heat recovery mode. The initial heat recovery heating capacity fraction identifies the fraction of available heat recovery mode heating capacity at the start of the transition period, the heat recovery heating capacity time constant identifies the time needed to recover to 99% of the steady-state value. This exponential model used for modeling the transition period can be turned off by setting the initial heat recovery heating capacity fraction to 1.0.

Heat Recovery Heating Capacity Fraction and Time Constant[LINK]

When the heat pump changes operating modes (i.e., from heating only to heat recovery mode), the transition does not happen immediately. There is a time delay, and a period of time where the available heating capacity is different from the steady-state capacity available after the system has had time to adjust. When this type of transition is modeled, an exponential decay curve is used. At the start of the transition period, only a fraction (kheat) of the steady-state capacity in heat recovery mode is available. The remaining capacity is recovered over a period of 5 time constants (tc,heat) and is modeled with the following equation. This equation was used to replicate both the cooling capacity and condenser power curves in the previous figure.

˙QHR,avail,heating=kheat(˙QHR,heating,total)+(1kheat)(˙QHR,heating,total)(1e(t/τc,heat))

where:

kheat is the fraction of steady-state heating capacity at beginning of heat recovery mode

τc,heat is the time constant to reach steady-state operation, 5 time constants equals 99% of steady-state capacity (hr)

QHR,avail,heating is the available heating capacity in heat recovery mode (W)

(1e(t/τc,heat)) is the multiplier used during the transition period from heating mode to heat recovery mode. This multiplier can be viewed using the report variable “VRF Heat Pump Heat Recovery Status Change Multiplier” (ref: Output:Variable).

Defrost Energy Input Ratio Modifier (function of temperature)[LINK]

The defrost energy input ratio (EIR) modifier curve (function of temperature) is a bi-quadratic curve with two independent variables: heating coil average entering air wet-bulb temperature and outdoor air dry-bulb temperature. The output of this curve is multiplied by the heating coil capacity, the fractional defrost time period and the runtime fraction of the heating coil to give the defrost power at the specific temperatures at which the coil is operating. The use of outdoor wet-bulb temperature is explicit in this curve object and the defrost equation below and is independent of the selection for Heating Performance Curve Outdoor Temperature Type.

DefrostEIRTempModFac=a+b(Twb,avg)+c(Twb,avg)2+d(Tc)+e(Tc)2+f(Twb,avg)(Tc)

The heat pump condenser’s total heating load is then calculated as the quotient of the total terminal unit capacity and the piping correction factor. Additional load due to defrost is also included (see following section).

˙Qheating,total=˙Qheating,TerminalUnitsPcorrection,heating+˙Qdefrost

The heat pump’s total (gross) heating capacity is then calculated based on the capacity correction factor as a function of temperatures. The impact of defrost on total heat pump heating capacity is also accounted for (see following section). The part-load ratio of the heat pump condenser can then be calculated.

˙QHP,total=˙Qtotal,rated(CAPFTHP,heating)(CRheating,correction)(HCAPdefrost)

PLR=˙Qheating,total˙QHP,total

where:

˙QHP,total is the heat pump total available heating capacity (W)

PLR is the heat pump part-load ratio

HCAPdefrost is the heating capacity correction factor for defrost mode.

As described for cooling operation, the available heating capacity of the heat pump is compared to the requested heating capacity of all indoor terminal units. If the requested heating capacity of all indoor terminal units is greater than the available heating capacity of the heat pump, one or more of the indoor terminal unit’s capacity is limited such that the sum of the zone terminal unit heating demand plus piping losses is equal to the available heat pump heating capacity (including the impact due to defrost).

The electrical power consumed by the heat pump condenser is calculated based on the heat pump’s nominal electric power consumption in heating mode and the actual operating conditions. A performance correction factor is used to correct for off-design performance as follows:

Heating Part-Load Fraction Correlation (function of heat pump cycling ratio)

The part-load fraction correlation (function of heat pump cycling ratio) is a linear, quadratic or a cubic curve with the independent variable being cycling ratio (part-load ratio / minimum part-load ratio). The output of this curve is used in combination with the rated EIR and EIR modifier curves to give the “effective” EIR for a given simulation time step. The cycling ratio correlation accounts for startup losses of the heat pump’s compression system.

CyclingRatioFrac=a+b(CyclingRatio)+c(CyclingRatio)2

or

CyclingRatioFrac=a+b(CyclingRatio)+c(CyclingRatio)2+d(CyclingRatio)3

where

CyclingRatio=PLRPLRmin

The cycling ratio can be viewed using the report variable “VRF Heat Pump Cycling Ratio” (ref: Output:Variable).

The cycling ratio correlation should be normalized to a value of 1.0 when the cycling ratio equals 1.0 (i.e., no efficiency losses when the compressor(s) run continuously for the simulation time step). For cycling ratio values between 0 and 1 (0 CyclingRatio < 1), the following rules apply:

CyclingRatioFrac0.7 and CyclingRatioFracCyclingRatio

If CyclingRatioFrac < 0.7 a warning message is issued, the program resets the CyclingRatioFrac value to 0.7, and the simulation proceeds. The runtime fraction of the heat pump is defined as CyclingRatio/CyclingRatioFrac. If CyclingRatioFrac < CyclingRatio, then a warning message is issued and the runtime fraction of the heat pump is limited to 1.0.

Heat Pump Runtime Fraction(HPRTF)=CyclingRatioCyclingRatioFrac

The heat pump runtime fraction can be viewed using the report variable “VRF Heat Pump Runtime Fraction” (ref: Output:Variable).

Since manufacturers data do not typically provide information defining the cycling losses of VRF AC systems, a typical part-load fraction correlation for a conventional DX refrigeration system (e.g., residential or small commercial unit) may be substituted here as:

CyclingRatioFrac=0.85+0.15(CyclingRatio)

Heating Energy Input Ratio Modifier (function of temperature)

As described previously (Ref. Application of Dual Performance Curves), the heating energy input ratio modifier curve can either be applied as a single curve object as shown below, or applied using dual performance curve objects. The single curve object application is discussed here.

EIRFTheating=a+b(Tdb,avg)+c(Tdb,avg)2+d(Tc)+e(Tc)2+f(Tdb,avg)(Tc)

where EIRFTheating = heating energy input ratio correction factor (function of temperature) (0-1).

Heating Energy Input Ratio Modifier (function of part-load ratio)

IF(PLR1)THEN    EIRFPLRheating=a+b(PLR)+c(PLR)2+d(PLR)3ELSE    EIRFPLRheating=e+f(PLR)+g(PLR)2+h(PLR)3

where:

EIRFPLRheating is the heating energy input ratio modifier (function of part-load ratio)

a-d are the coefficients for the heating energy input ratio modifier when part-load ratio 1

e-h are the coefficients for the heating energy input ratio modifier when part-load ratio > 1

The total power consumed by the heat pump condenser in heating mode is then calculated. A correction for power consumed due to defrost is also included (see following section).

The total power consumed by the heat pump condenser in heating mode is based on the user specified coefficient of performance (COP). A correction for power consumed due to defrost is also included (see following section).

HeatingPower=(˙Qheating,total,ratedCAPFTheatingCOPheating,reference)(EIRFTheating)(EIRFPLRheating)(HPRTF)(EIRdefrost)

The heating COP is then calculated as the ratio of the outdoor unit heating capacity divided by the total electric consumption rate of all components associated with the outdoor unit. Evaporatively-cooled condenser pump power is typically 0 during heating mode operation.

COPheating=˙Qheating,total/(HeatingPower+PCrankcaseHeater+PEvapCoolerPower+Pdefrost)

where:

COPheating is the operating coefficient of performance, report variable “VRF Heat Pump Heating COP”

HeatingPower is the operating electric consumption rate (W), report variable “VRF Heat Pump Heating Electric Power”.

When operating in heat recovery mode, equations similar to those used for available heating capacity are used to model operating power.

HREIRModHP,heating=a+b(Tent)+c(Tent)2+d(Tcond)+e(Tcond)2+f(Tent)(Tcond)

This equation can be used to provide a constant fractional difference for heating electric consumption rate in heat recovery mode (i.e., only *a* is non-zero) or a fractional term that varies with indoor and outdoor conditions. With very limited performance data available at this time, it is recommended that only the constant (a) term be used at this time. When the VRF system is not operating in heat recovery mode, this fraction is set to 1.The available heating electric consumption rate in heat recovery mode is then calculated as:

HeatingPowerHR=HeatingPower(HREIRModHP,heating)

PowerHR=kEIR,heat(HeatingPowerHR)+(1kEIR,heat)(HeatingPowerHR)(1e(t/τc,EIR,heat))

where:

kEIR,heat is the fraction of steady-state heating electric consumption rate at beginning of heat recovery mode

τc,EIR,heat is the time constant to reach steady-state operation, 5 time constants equals 99% of steady-state capacity (hr)

HeatingPowerHR is the heating electric consumption rate in heat recovery mode (W)

(1e(t/τc,EIR,heat)) = The multiplier used during the transition period from heating mode to heat recovery mode. This multiplier can be viewed using the report variable “VRF Heat Pump Heat Recovery Status Change Multiplier” (ref: Output:Variable).

Operating Coefficient of Performance[LINK]

Similar to the cooling and heating COP report variables, a report variable is included to identify the overall COP of the system. The numerator represents the total cooling and heating coil capacities (or loads) where piping losses have been accounted for. If heat recovery is not used only one of these terms is non-zero. When heat recovery is used, one or both of these terms can be non-zero, therefore, the operating COP includes recovered energy. The denominator includes the electric consumption rates of all system components. For water-cooled VRF AC systems, the plant pump power is not included.

COPoperating=(˙Qcooling,TerminalUnits+˙Qheating,TerminalUnits)(OutdoorUnitPower+PCrankcaseHeater+PEvapCoolerPower+Pdefrost+PTU,fan+PTU,parasitic)

where:

OutdoorUnitPower is the cooling or heating electric consumption rate of outdoor unit (W)

PTU,fan is the electric consumption rate of all terminal unit fans (W)

PTU,parasitic is the electric consumption rate of all terminal unit parasitic electric (W).

Defrost Adjustment Factors[LINK]

Frost formation on the outdoor coil, and the need to periodically defrost this coil, has a significant impact on heating capacity and energy use by the DX heating system. This model uses a timed reverse-cycle defrost model. If the outdoor air dry-bulb temperature is below the specified maximum temperature for defrost operation, then the model calculates adjustment factors for heating capacity and input power due to frost formation This method of accounting for the impacts of frosting/defrost was taken from the model used in DOE-2.1E (ESTSC 2001, Miller and Jaster 1985).

The model first estimates the outdoor coil temperature according to a linear empirical relationship with outdoor air dry-bulb temperature as the independent variable.

Tcoil,out=0.82Tdb,o8.589

The difference between the outdoor air humidity ratio (from the weather file) and the saturated air humidity ratio at the estimated outdoor coil temperature is then calculated, and this value is used as an indication of frost formation on the outdoor coil.

Δωcoil,out=MAX[1.0E6,ωoutdoorωsat(Tcoil,out,OutBaroPress)]

Frost formation on the outdoor coil must be periodically removed. The fraction of compressor runtime when the coil is being defrosted is entered by the user. Adjustment factors to total heating coil capacity and input power due to frost formation on the outdoor coil are calculated by an empirical model with Δωcoil,out as the independent variable as shown below.

HCAPdefrost=0.909107.33(Δωcoil,out)

EIRdefrost=0.936.45(Δωcoil,out)

If the outdoor air dry-bulb temperature is above the specified maximum temperature for defrost operation, the fractional defrost time period is set to zero and the heating capacity/input power multipliers are set to unity (1).

Defrost Operation[LINK]

If the fractional defrost time period is greater than zero for the simulation time step, then the model calculates the electrical power used during defrost. The additional heating load due to defrost (indoor cooling during defrost) is also calculated so that it may be added to the existing heating load when calculating input power for the compressor(s) and outdoor coil fan(s).

Qdefrost=0.01(tfrac,defrost)(7.222Tdb,o)(Qtotal,rated1.01667)

Pdefrost=DefrostEIRTempModFac(Qtotal,rated1.01667)(tfrac,defrost)(RTF)

where:

Qdefrost is the additional indoor heating load due to reverse-cycle defrost (W)

Qtotal,rated is the total full-load heating capacity of the coil at rated conditions (W)

Pdefrost is the report variable for average defrost power for the simulation time step (W)

DefrostEIRTempModFac is the energy input ratio modifier curve applicable during defrost

tfrac,defrost is the defrost time period fraction specified by user

RTF=(PLRPartLoadFrac)=runtime fraction of the heating coil

Physics based VRF Heat Pump Model (VRF-FluidTCtrl-HP)[LINK]

Overview[LINK]

This section describes the Physics based VRF Model (VRF-FluidTCtrl) to simulate the energy performance of Variable Refrigerant Flow (VRF) systems in the heat pump (HP) operation mode, i.e., either cooling or heating is provided but not simultaneously.

Compared with the empirical system curves based VRF-SysCurve model, the VRF-FluidTCtrl model adheres to a physics-based development providing the ability to simulate the refrigerant loop performance and consider the dynamics of more operational parameters. This is essential for the representation of many enhanced control logics, e.g., the adjustment of superheating degrees during low load operations. Furthermore, this model implements component-level curves rather than the system-level curves, and thus requires much fewer curves as model inputs. The definition of VRF performance data for simulation as part of ASHRAE Standard 205 will ensure VRF manufacturers provide adequate data to create the required performance curves.

The main features of the VRF-FluidTCtrl model include:

  • Introduction of separate curves for capacities and power inputs of indoor and outdoor units instead of overall curves for the entire system.

  • Allowing of variable evaporating and condensing temperatures in the indoor and outdoor units under various operational conditions.

  • Allowing of variable fan speed based on the temperature and zone load in the indoor unit.

  • Allowing of further modifications of operational parameters (e.g., evaporating temperature, superheating degrees, and supply air flow rate) during low load conditions.

  • Development of a physics-based model to calculate thermal loss in the refrigerant piping network (varies with refrigerant flow rate, operational conditions, pipe length, and pipe and insulation materials) instead of a constant correction factor.

These features enhance the simulation accuracy of the VRF system performance in both heating and cooling modes, especially during low load operations. In addition, these features enable the modeling of a multi-stage control strategy to adjust the system capacity during low load conditions, as shown in Figure VRF-FluidTCtrl-1. Moreover, the VRF-FluidTCtrl model enables the potential simulation of demand response of VRF systems by directly slowing down the speed of compressors in the outdoor units with inverter technology.

Multi-stage control strategy to adjust the system capacity during low load conditions [fig:multi-stage-control-strategy-to-adjust-the-system-capacity-during-low-load-conditions]

Multi-stage control strategy to adjust the system capacity during low load conditions [fig:multi-stage-control-strategy-to-adjust-the-system-capacity-during-low-load-conditions]

Note that a number of calculation steps are coupled together in the VRF-FluidTCtrl model, for instance, the piping loss calculation and the system performance calculation. More specifically, the piping loss changes the operating conditions of the system, which may lead to a different control strategy and thus in reverse affect the amount of piping loss. This makes it difficult to obtain an analytical solution for a number of operational parameters (e.g., enthalpy of refrigerant entering the indoor unit), and therefore numerical iterations are employed to address this problem (refer to Figure VRF-FluidTCtrl-3 for more details). Therefore, the VRF-FluidTCtrl model can have a longer execution time to perform the simulation than the VRF-SysCurve model.

The object connections for VRF-FluidTCtrl model is similar to those for VRF-SysCurve model. The difference lies in the object used to describe a specific components. For example, VRF-SysCurve model uses AirConditioner:VariableRefrigerantFlow object to describe the VRF outdoor unit performance, while in VRF-FluidTCtrl model the AirConditioner:VariableRefrigerantFlow:FluidTemperatureControl object is used.

EnergyPlus object connections for VRF-FluidTCtrl model [fig:energyplus-object-connections-for-vrf-fluidtctrl-model]

EnergyPlus object connections for VRF-FluidTCtrl model [fig:energyplus-object-connections-for-vrf-fluidtctrl-model]

Model Description[LINK]

Figures 11, 12, 13, and 14 illustrate the overall logic of the VRF-FluidTCtrl algorithms. Detailed calculation procedures and explanations of equations are described in the following sections.

Schematic Pressure-Enthalpy Diagram for VRF Operation – Heating Mode [fig:vrf-fluid-t-ctrl-2a]

Schematic Pressure-Enthalpy Diagram for VRF Operation – Heating Mode [fig:vrf-fluid-t-ctrl-2a]

Schematic Pressure-Enthalpy Diagram for VRF Operation – Cooling Mode [fig:vrf-fluid-t-ctrl-2b]

Schematic Pressure-Enthalpy Diagram for VRF Operation – Cooling Mode [fig:vrf-fluid-t-ctrl-2b]

Flow Chart of the VRF-FluidTCtrl Model – Heating Mode [fig:vrf-fluid-t-ctrl-3a]

Flow Chart of the VRF-FluidTCtrl Model – Heating Mode [fig:vrf-fluid-t-ctrl-3a]

Flow Chart of the VRF-FluidTCtrl Model – Cooling Mode [fig:vrf-fluid-t-ctrl-3b]

Flow Chart of the VRF-FluidTCtrl Model – Cooling Mode [fig:vrf-fluid-t-ctrl-3b]

The VRF-FluidTCtrl algorithms are described in details below, including the models of both the indoor and outdoor unit(s) of the VRF system. Because of the different control algorithms used in the outdoor unit for cooling and heating modes, the heating and cooling modes are described separately.

Modeling of the indoor unit (I/U) – Part I[LINK]

Step 1.1: Obtaining zonal load/condition information

Obtain the following information for each zone from the zone modules within EnergyPlus:

  • zone sensible loads Qin,sensible

  • zone total loads Qin,total

  • indoor air temperature Tin

  • indoor air humidity ratio Win

The operation mode of the VRF system can therefore be determined from Qin,sensible.

Step 1.2: Calculate I/U required evaporator/condenser temperature

Evaluate the required coil surface air temperature Tfs and the required evaporator/condenser refrigerant temperature Te,req/Tc,req for each indoor unit.

Assuming the coil air flow rate is at the maximum (Ga,rate), the entering air temperature after the indoor coil fan can be calculated for each indoor unit:

Tcoil,in=Tin+Qfan/(Ga,rate×cp×ρin)

where:

cp is the air heat capacity (kJ/(kgK))

Ga,rate is the volumetric flow rate of the air through the indoor unit, at the rated conditions (m3/s)

Qfan is the sensible heat released by fan (W)

Tcoil,in is the temperature of air entering the indoor coil (°C)

ρin is the density of indoor air, fρ(Tin,Win) (kg/m3).

The supply air temperature Tout can be calculated by:

Tout=Tcoil,inQin,sensible/(Ga,rate×cp×ρin)

where Tout is the temperature of the air leaving the indoor unit (supply air) (C).

The coil surface air temperature Tfs can be calculated, given the bypass factor BF:

Tfs=Tcoil,in(Tcoil,inTout)/(1BF)

where:

BF is the bypass factor (different values for cooling mode and heating mode)

Tfs is the indoor unit coil surface air temperature (C).

In the cooling mode, the superheating degree SH is given and thus the required evaporating temperature Tereq for the associated indoor unit can be obtained:

Tereq=Tfs[AcSH2+BcSH+Cc]

Where

Ac, Bc, Cc are coefficients

SH is the superheating degrees for the indoor unit (C)

SHref is the reference superheating degrees for the indoor unit (e.g., 3) (C)

Tereq is the required evaporating temperature for a specific indoor unit (C).

In the heating mode, the subcooling degree SC is given and thus the required condensing temperature Tcreq for the associated indoor unit can be obtained:

Tcreq=Tfs+[AhSC2+BhSC+Ch]

where:

Ah, Bh, Ch are coefficients

SC is the subcooling degrees for the indoor unit (C)

SCref is the reference subcooling degrees for the indoor unit (e.g., 5) (C)

Tcreq is the required condensing temperature for a specific indoor unit (C).

Step 1.3: Calculate effective evaporating /condensing temperature

There are two refrigerant temperature control strategies for the indoor unit: (1) ConstantTemp, (2) VariableTemp.

  • In the ConstantTemp strategy, Te and Tc are kept at constant values provided by the user.

  • In the VariableTemp strategy, Te and Tc are determined using the equations below:

In the cooling mode:

Te=min{max[min(Tereq,i),Temin],Temax}

where:

Tereq,i is the required evaporating temperature for the ith indoor unit (C)

Temin is the minimum evaporating temperature for the indoor unit (e.g. 3) (C)

Temax is the maximum evaporating temperature for the indoor unit (e.g. 13) (C).

In the heating mode:

Tc=min{max[min(Tcreq,i),Tcmin],Tcmax}

where:

Tcreq,i is the required condensing temperature for the ith indoor unit (C)

Tcmin is the minimum condensing temperature for the indoor unit (e.g. 42) (C)

Tcmax is the maximum condensing temperature for the indoor unit (e.g. 46) (C).

Modeling of the outdoor unit (O/U) - Cooling Mode[LINK]

Step 2c.1: Piping loss calculation in the cooling mode

Piping loss occurs when the refrigerant flows through the main pipe connecting the indoor and outdoor units. It includes both the refrigerant pressure drop and heat loss. It affects the VRF system operation in several ways. First, the heat loss creates an extra load to the system and leads to higher energy consumptions. Second, the pressure drop and heat loss changes the operational conditions of the compressor (i.e., compressor suction pressure and compressor suction temperature) and thus affects its operational performance. Therefore, the piping loss calculation is critical for the analysis of the whole system.

Note that the change of compressor operational conditions may lead to different control strategies of the system, which in reverse affects the amount of piping loss. So the piping loss analysis and system performance analysis are coupled together. Numerical iterations are designed to address the coupling effect, as described below.

Calculate the refrigerant flow rate for each indoor unit using assumed hHexin:

Gi=Qi/(hHexout,ihHexin)

where:

Gi is the refrigerant flow rate for the ith indoor unit (kg/s)

Qi is the total cooling load of the ith zone (W)

hHexin is the enthalpy of the refrigerant entering the indoor unit (kJ/kg)

hHexout,i is the enthalpy of the refrigerant leaving a specific indoor unit (kJ/kg)

hHexout,i is a function of Pe, Te, and SHi. It can be calculated using refrigerant thermodynamic property equations fg_h(Pe,Te+SHi).

The enthalpy of the refrigerant in the main loop hHexout can be obtained by calculating the weighted average of the enthalpy of the refrigerant in each indoor unit:

hHexout=(hHexout,i×Gi/Gtot)

Gtot=Gi

where:

Gtot is the refrigerant flow rate in the main loop (kg/s)

hHexout is the enthalpy of the main loop refrigerant leaving indoor units (kJ/kg).

The viscosity of the refrigerant within the pipe can be calculated via:

kv,1=Pe/4926000

kv,2=hHexout/383.55

kv,3=(THexout+273.15)/344.39

μ=106×(4.302×kv,1+0.81622×kv,12120.98×kv,2+139.17×kv,22+118.76×kv,3+81.04×kv,32+5.7858×kv,1×kv,28.3817×kv,1×kv,3218.48×kv,2×kv,3+21.58)

where:

μ is the viscosity of the refrigerant within the pipe (Pa-s)

kv,i is the coefficients to calculate the refrigerant viscosity

THexout is the temperature of main loop refrigerant leaving indoor units, fg_t(Pe,hHexout) (C)

fg_t are functions calculating the temperature of superheating refrigerant.

Given hHexout and Pe=fs_p(Te), the following dimensionless quantities describing the refrigerant flow state can be obtained:

Re=Gtot/3600/(0.25×π×D2)×D/μ

Pr=μ×fg_Cp(Pe,hHexout)×0.001/fg_λ(Pe,hHexout)

Nu=0.023×Re0.8×Pr0.3

St=Nu/Re/Pr

where:

Pe is the evaporating pressure (Pa)

Re is the Reynolds number

Pr is the Prandtl number

Nu is the Nusselt number

St is the Stanton number

μ is the viscosity of the refrigerant within the pipe (Pa-s)

fs_p are functions calculating the pressure of saturated refrigerant

fg_Cp are functions calculating the specific heat of superheating refrigerant

fg_λ are functions calculating the conductivity of superheating refrigerant.

Then the piping pressure loss ΔPpipe can be obtained using the above dimensionless quantities:

ΔPpipe=8×St×Pr2/3×L/D×fg_ρ(Pe,hHexout)×V2/2+H×fg_ρ(Pe,hHexout)×9.80665

Ps=PeΔPpipe

where:

fg_ρ are functions calculating the density of superheating refrigerant

D is the main pipe diameter (m)

H is the height difference between the outdoor unit node and indoor unit node of the main pipe (m)

L is the main pipe length (outdoor unit to the first branch joint) (m)

Ps is the compressor suction pressure (Pa)

ΔPpipe is the pressure drop in the pipe (Pa)

V is the refrigerant flow velocity (m/s).

The suction saturated temperature Te (i.e., saturated vapor temperature corresponding to the compressor suction pressure) can be obtained via refrigerant thermodynamic property equations:

Te=fs_t(PeΔPpipe)

where:

fs_t are functions calculating the temperature of saturated refrigerant

Te is the suction saturated temperature at the compressor inlet (C).

The heat loss through the pipe can be obtained via:

k1=Nu×fg_λ(Pe,hHexout)

k2=2×λi/Ln(1+2×Wi/D)

k3=h×(D+2×Wi)

Qpipe=(π×L)×(TaTHexout)/(1/k1+1/k2+1/k3)

where:

fg_λ are functions calculating the conductivity of superheating refrigerant

Qpipe is the heat loss through the pipe (W)

Ta is the average of outdoor air temperature and indoor temperature (C)

THexout is the temperature of main pipe refrigerant leaving indoor units (C)

wi is the pipe insulation thickness (m)

ki are coefficients for the piping loss calculation.

Note that Qpipe is calculated using an assumed hHexin. Its value affects the compressor operation calculations as shown in Steps 2c.2 - 2c.6 and may change the value of condensing temperature Tc. This leads to an updated hHexin=f(Pc,TcSC). If the difference between the calculated hHexin and the assumed hHexin is higher than the assigned tolerance (5%), a new round of iteration is performed using the calculated hHexin.

Step 2c.2: Initialize O/U operation conditions

The compressor operation is coupled with a number operational parameters, which makes it difficult to obtain an analytical solution. In the algorithm, an iterative approach is designed to obtain the energy consumption of the compressor (Step 2c.2 – 2c.6).

For the first iteration,

  • Initialize outdoor unit SC with the reference value (from IDF input, e.g., 5C)

  • Initialize the compressor power Ncomp with the value calculated from the reference COP (e.g., 3.5):

Ncomp=iQin,total,iCOP

where:

Ncomp is the compressor power (W)

Ncomp is the assumed compressor power for the first iteration (W)

Qin,total,i is the total cooling load for zone i (W).

For the following iterations,

  • Initialize SC with the calculated value in the previous iteration

  • Initialize the compressor power Ncomp with the calculated value in the previous iteration

The heat rate released from the outdoor unit can be calculated by:

Qout=iQin,total,i+Ncomp+Qpipe

where:

Qout is the heat rate released from the outdoor unit (W)

Qpipe is the heat loss through the pipe (W).

Step 2c.3: Calculate O/U effective condensing temperature

(1) Calculate the required coil surface air temperature Tfs for the outdoor unit.

The temperature of the air leaving the outdoor unit Tout can be calculated by:

Tout=Tin+Qout/(Ga,rate×cp×ρin)

where:

Ga,rate is the volumetric flow rate of the air through the outdoor unit, at the rating conditions (m3/s)

Tout is the temperature of the air leaving the outdoor unit (C)

Tin is the temperature of the air entering the outdoor unit, i.e., outdoor air (C)

ρin is the air density at outdoor conditions, fρ(Tin,Win) (kg/m3).

The condenser surface temperature Tfs can be calculated by:

Tfs=Tin+(ToutTin)/(1BFc)

where BFc is the bypass factor for the outdoor unit.

(2) Calculate required condensing temperature for the outdoor unit Tcreq and then the effective condensing temperature Tc (between 42-46C)

Tcreq=Tfs+[AhSC2+BhSC+Ch]

Tc=Tcreq

where:

Ah, Bh, Ch are coefficients

SC is the subcooling degrees for the outdoor unit (C)

SCref is the reference subcooling degrees for the outdoor unit (C)

Tcreq is the required condensing temperature for the outdoor unit (C)

Tc is the effective condensing temperature (C).

Step 2c.4: Calculate required Loading Index

Loading Index is used to describe the operational mode of the compressor system. The first index represents minimal capacity operation, while the last index represents full capacity operation.

The required Loading Index can be calculated by the following procedures.

  1. Calculate the evaporative capacity at a variety of Loading Index:

Mcap=(r1+r2Tc+r3Te+r4T2c+r5TcTe+r6Te2)

Qrps=Ccap,system×Mcap×Qref

where:

Ccap,system is the evaporative capacity correction factor, describing the possible system configuration difference between test bed and real system (a constant value obtained from manufacturer data)

r1, …, r6 are the empirical coefficients corresponding to rps

rps is the compressor speed (r/s)

Te is the suction saturated temperature at the compressor inlet (C)

Tc is the effective condensing temperature (C)

Mcap is the multiplier for the evaporative capacity calculation (–)

Qrps is the evaporative capacity corresponding to rps (W)

Qref is the rated evaporative capacity (W).

An example of resulting capacity for different rps (at Tc=36C and Te=9C) is presented in Table 1.

Evaporative capacity at different Loading Index
Loading Index (rps) Evaporative Capacity (kW), Tc=36C, Te=9C Loading Index (rps) Evaporative Capacity (kW), Tc=36C, Te=9C
1 (18) 3.951 2 (25) 5.400
3 (30) 6.435 4 (36) 8.516
5 (45) 10.756 6 (60) 14.445
7 (75) 18.412 8 (90) 22.045
9 (108) 26.162
  1. Find the rps range that covers the required evaporative capacity Qrps,modify .

Qrps,modify=Ccap,operation×(Qin,total+Qpipe)

Ccap,operation=Ccap,density×Ccap,enthalpy

Ccap,density=ρtest/ρreal

Ccap,enthalpy=hEvapout,testhEvapin,testhCompin,realhEvapin,real

hCompin,real=hHexout,real+Qpipe/Gtot

where:

Ccap,operation is the evaporative capacity correction factor, describing the operational difference between test cases and real cases (i.e., SH and SC )

Ccap,density is the evaporative capacity correction factor, describing the variations of refrigerant density at test conditions and real operational conditions

Ccap,enthalpy is the evaporative capacity correction factor, describing the variations of refrigerant enthalpy at test conditions and real operational conditions

Gtot is the refrigerant flow rate in the main loop (kg/s)

hEvapin,real is the enthalpy of refrigerant entering the evaporators (IU) at real conditions (kJ/kg)

hEvapout,real is the average enthalpy of refrigerant leaving the evaporators (IU) at real conditions (kJ/kg)

hEvapin,test is the enthalpy of refrigerant entering the evaporator at test conditions (i.e., SC at test conditions (e.g., 5C) and Tc) (kJ/kg)

hEvapout,test is the enthalpy of refrigerant leaving the evaporator at test conditions (i.e., SH at test conditions (e.g., 8C) and Te) (kJ/kg)

hCompin is the enthalpy of refrigerant entering the compressor (kJ/kg)

Qpipe is the heat loss through the pipe (W).

For example, if the required capacity is 8 kW, the rps range is 30 to 36.

  1. Calculate the rps that meets the required capacity by interpolation. In the above example, the resulting rps is 34.5 rps.

  2. If the calculated rps is lower than the minimum rps (e.g. 18 rps ), go to Step 2c.5, otherwise skip Step 2c.5 and directly go to Step 2c.6.

Step 2c.5: Modify evaporating temperature to further reduce outdoor unit capacity

If the calculated rps is lower than the minimum rps (e.g. 18 rps ), it means that the zone cooling load is even lower than the system evaporative capacity corresponding to the minimum compressor speed. In this situation, the evaporating temperature Te as well as the superheating degree SH is modified to further reduce the outdoor unit capacity. More specifically:

  1. Set rps at its minimum value (e.g., 18 rps).

  2. Update Te to meet the required evaporative capacity, using equations described in Step 2c.4a.

  3. Update Te to meet the updated Te . Note that due to the Te updates, the refrigerant state and flow rate are changed and thus the piping loss analysis should also be repeated (Step 2c.1). So is the calculation of Ccap,operation (Step 2c.2 - 2c.3).

  4. SH can be updated based on the updated Te, using the equations shown in Step 1.2.

Step 2c.6: Calculate the compressor power

  1. Calculate the compressor power by the following procedures.
<!-- -->
  1. Calculate the compressor power at a variety of loading index using the following equation. The resulting table from the same example used above is shown below in Table 2.

Mcomp=c1+c2Tc+c3Te+c4T2c+c5TcTe+c6Te2

Ncomp,rps=Mcomp×Ncomp,ref

where:

c1, …, c6 are the empirical coefficients corresponding to rps

Te is the suction saturated temperature at the compressor inlet (C)

Tc is the effective condensing temperature (C)

Mcomp is the multiplier for the compressor power calculation (–)

Ncomp,ref is the rated compressor power (W)

Ncomp,rps is the compressor power corresponding to rps (W).

Outdoor unit compressor power at different Loading Index
Loading Index (rps) Evaporative Capacity (kW), Tc=36C, Te=9C
1 (18) 3.951
2 (25) 5.400
3 (30) 6.435
4 (36) 8.516
5 (45) 10.756
6 (60) 14.445
7 (75) 18.412
8 (90) 22.045
9 (108) 26.162
  1. According to the rps range determined, calculate the compressor power Ncomp by interpolation. In the above example, the compressor power is 1.155 kW.
<!-- -->
  1. Compare the calculated Ncomp above with the initialized Ncomp in Step 2c.2:
  • If NcompNcomp>δ then go to Step 2c.2 for a new round of iteration.

  • Else, end the iteration and go to Step 2c.7.

Step 2c.7: Total power consumption of the outdoor unit[step-2c.7-total-power-consumption-of-the-outdoor-unit]

Calculate the total electric power consumption by the outdoor unit:

Nout=Nfan+Ncomp/einv

where:

einv is the efficiency of the inverter of compressor

Nfan is the electric power consumption by the outdoor fan (W)

Nout is the total electric power consumption by the outdoor unit (W).

Modeling of the outdoor unit (O/U) - Heating Mode[LINK]

Step 2h.1: Piping loss calculation in the heating mode

Similarly to that in the cooling mode, the piping loss analysis in the heating mode needs to address the coupling effect between the piping loss and system operations. However, due to the control strategy difference between the two modes, the piping loss algorithm in the heating mode requires different calculation steps and more numerical iterations, as shown below.

Calculate the refrigerant flow rate for each indoor unit using assumed hHexin:

Gi=Qi/(hHexinhHexout,i)

where:

Gi is the refrigerant flow rate for the ith indoor unit (kg/s)

Qi is the total heating load of the ith zone (W)

hHexin is the enthalpy of the refrigerant entering the indoor unit (kJ/kg)

hHexout,i is the enthalpy of the refrigerant leaving a specific indoor unit (kJ/kg)

hHexout,i is a function of Pc, Tc, and SCi. It can be calculated using refrigerant thermodynamic property equations fgh(Pc,TcSCi).

The refrigerant flow rate in the main loop can be obtained by summing up the flow rate in each indoor unit:

Gtot=Gi

where Gtot refrigerant flow rate in the main loop (kg/s).

Calculate the average refrigerant pressure and enthalpy within the pipes, using assumed piping pressure loss ΔPpipe and hTdi :

Pave=Pc+ΔPpipe/2

have=(hHexin+hTdi)/2

where:

have is the average refrigerant enthalpy within the main pipe (kJ/kg)

hHexin is the enthalpy of the refrigerant entering the indoor unit (kJ/kg)

hTdi is the enthalpy of the refrigerant leaving the compressor (kJ/kg)

Pave is the average refrigerant pressure within the main pipe (Pa)

Pc is the condensing pressure (Pa)

ΔPpipe is the pressure drop in the main pipe (Pa).

The viscosity of the refrigerant within the pipe can be calculated given Pave and have, using the following equations:

kv,1=Pave/4926000

kv,2=have/383.55

kv,3=(Tave+273.15)/344.39

μ=106×(4.302×kv,1+0.81622×kv,12120.98×kv,2+139.17×kv,22+118.76×kv,3+81.04×kv,32+5.7858×kv,1×kv,28.3817×kv,1×kv,3218.48×kv,2×kv,3+21.58)

where:

μ is the viscosity of the refrigerant within the pipe (Pa-s)

kv,i are coefficients to calculate the refrigerant viscosity

have is the average refrigerant enthalpy within the pipes (kJ/kg)

Pave is the average refrigerant pressure within the pipes (Pa)

Tave is the average temperature of refrigerant leaving indoor units, which corresponds to Pave and have (C)

Given Pave and have, the following dimensionless quantities describing the refrigerant flow state can be obtained:

Re=Gtot/3600/(0.25×π×D2)×D/μ

Pr=μ×fgCp(Pave,have)×0.001/fgλ(Pave,have)

Nu=0.023×Re0.8×Pr0.4

St=Nu/Re/Pr

where:

Re is the Reynolds number

Pr is the Prandtl number

Nu is the Nusselt number

St is the Stanton number

μ is the viscosity of the refrigerant within the pipe (Pa-s)

fgCp are functions calculating the specific heat of superheating refrigerant

fgλ are functions calculating the conductivity of superheating refrigerant.

Then, the piping pressure loss ΔPpipe can be obtained using the above dimensionless quantities:

ΔPpipe=8×St×Pr2/3×L/D×fgρ(Pave,have)×V2/2H×fgρ(Pave,have)×9.80665

where:

fgρ are functions calculating the density of superheating refrigerant

D is the main pipe diameter (m)

H is the height difference between the outdoor unit node and indoor unit node of the main pipe (m)

L is the main pipe length (outdoor unit to the first branch joint) (m)

St is the Stanton number

Re is the Reynolds number

Pr is the Prandtl number

ΔPpipe is the pressure drop in the pipe (Pa)

V is the refrigerant flow velocity (m/s).

Compare the calculated ΔPpipe and the one assumed above. If the difference is higher than the assigned tolerance (5%), a new round of iteration is performed using the calculated ΔPpipe.

The compressor discharge saturated temperature Tc (i.e., saturated vapor temperature corresponding to compressor discharge pressure) can be obtained via:

Tc=fst(Pc+ΔPpipe)

where:

fst are functions calculating the temperature of saturated refrigerant

Tc is the discharge saturated temperature at the compressor outlet (C)

Pc is the condensing pressure (Pa)

ΔPpipe is the pressure drop in the pipe (Pa).

The heat loss through the pipe can be obtained via:

k1=Nu×fgλ(Pave,have)

k2=2×λi/Ln(1+2×Wi/D)

k3=h×(D+2×Wi)

Qpipe=(π×L)×(TaTHexin)/(1/k1+1/k2+1/k3)

where:

fgλ are functions calculating the conductivity of superheating refrigerant

Qpipe is the heat loss through the pipe (W)

Ta is the average of outdoor air temperature and indoor temperature (C)

THexin is the average of entering indoor units (C)

wi is the pipe insulation thickness (m)

ki is the coefficients for the piping loss calculation.

The enthalpy of the refrigerant entering the indoor unit can be updated using the calculated Qpipe:

hTdi=hHexin+Qpipe/Gtot

where:

hHexin is the enthalpy of the refrigerant entering the indoor unit (kJ/kg)

hTdi is the enthalpy of the refrigerant leaving the compressor (kJ/kg)

Gtot is the refrigerant flow rate in the main loop (kg/s).

Compare the calculated hTdi and the one assumed above. If the difference is higher than the assigned tolerance (5%), a new round of iteration is performed using the calculated hTdi.

Note that Qpipe is calculated using an assumed hHexin at the beginning of the piping loss algorithm. Its value affects the compressor operation calculations as shown in Step 2h.2~2h.6 and may change the value of evaporating temperature Te. This leads to an updated hHexin=f(Pe,Te+SH,Ncomp,Gtot). If the difference between the calculated result and the assumed value is higher than the assigned tolerance (5%), a new round of iteration is performed using the calculated hHexin.

Step 2h.2: Initialize O/U operation conditions

Similar to that in cooling mode, an iteration approach is designed to determine the energy consumption of the compressor (Step 2h.2 to Step 2h.6).

For the first iteration,

  • Initialize outdoor unit SH with the reference value (from IDF input, e.g., 1.5C)

  • Iinitialize the compressor power Ncomp with the value calculated from the reference COP (e.g., 3.5):

Ncomp=iQin,total,iCOP

where:

Ncomp is the compressor power (W)

Ncomp is the assumed compressor power for the first iteration (W)

Qin,total,i is the total heating load for zone i (W).

For the following iterations:

  • Initialize SH with the calculated value in the previous iteration

  • Initialize the compressor power Ncomp with the calculated value in the previous iteration

Calculate the heat rate extracted by the outdoor unit by:

Qout=iQin,total,i+QpipeNcomp

where:

Qout is the heat rate extracted by the outdoor unit (W)

Qpipe is the heat loss through the pipe (W).

Step 2h.3: Calculate O/U effective evaporating temperature

  1. Calculate the required coil surface air temperature Tfs for the outdoor unit.

The enthalpy of the air leaving the outdoor unit can be calculated by:

Hfs=HinQout/(Ga,rate×ρo)/(1BF)

The coil surface air temperature Tfs can be calculated from the Hfs:

if Hfs<H98%,Wo, then Tfs=f(Hfs,98%)

if HfsH98%,Wo, then Tfs=f(Hfs,Wo)

where:

BF is the bypass factor for the outdoor unit

Ga,rate is the volumetric flow rate of the air through the outdoor unit, at the rating conditions (m3/s)

Hfs is the enthalpy of the air leaving the outdoor unit (kJ/kg)

Hin is the enthalpy of the air entering the outdoor unit, i.e., outdoor air (kJ/kg)

Wo is the humidity ratio of the outdoor air (kg/kg)

ρo is the density of the outdoor air (kg/m3).

  1. Calculate required evaporating temperature for the outdoor unit Tereq and then the effective evaporating temperature Te.

Tereq=Tfs[AcSH2+BcSH+Cc]

Te=Tereq

where:

Ac,Bc,Cc are coefficients

SH is the superheating degrees for the outdoor unit (C)

SHref is the reference superheating degrees for the outdoor unit (C)

Tereq is the required evaporating temperature for the outdoor unit (C)

Te is the effective evaporating temperature (C).

Step 2h.4: Calculate required compressor Loading Index

Calculate the required compressor Loading Index by the following procedures.

  1. Calculate the evaporative capacity at a variety of speeds:

Mcap=Ccap,system×(r1+r2Tc+r3Te+r4Tc2+r5TcTe+r6T2e)

Qrps=Ccap,system×Mcap×Qref

where:

Ccap,system is the evaporative capacity correction factor, describing the possible system configuration difference between test bed and real system (obtained from manufacturer data)

r1, …, r6 is the empirical coefficients corresponding to rps

rps is the compressor speed (r/s)

Te is the effective evaporating temperature (C)

Mcap is the multiplier for the evaporative capacity calculation (–)

Tc is the discharge saturated temperature at the compressor outlet (C)

Qrps is the evaporative capacity corresponding to rps (W)

Qref is the rated evaporative capacity (W).

An example of resulting capacity for different rps (at Tc = 36C and Te = 9C) is presented in Table 3.

Table 3– evaporative capacity at different Loading Index

  1. Find the rps range that covers the required evaporative capacity Qrps,modify.

Qrps,modify=Ccap,operation×(Qin,total+QpipeNcomp)

Ccap,operation=Ccap,density×Ccap,enthalpy

Ccap,density=ρtest/ρreal

Ccap,enthalpy=hEvapout,testhCondout,testhCompin,realhCondout,real

where:

Ccap,operation is the evaporative capacity correction factor, describing the operational difference between test cases and real cases (i.e., SH and SC)

Ccap,density is the evaporative capacity correction factor, describing the variations of refrigerant density at test conditions and real operational conditions

Ccap,enthalpy is the evaporative capacity correction factor, describing the variations of refrigerant enthalpy at test conditions and real operational conditions

hCondout,test is the enthalpy of refrigerant leaving the condensers (IU) at test conditions which corresponds to SH at test conditions (e.g., 8C) and Te) (kJ/kg)

hCondout,real is the average enthalpy of refrigerant leaving the condensers (IU) at real conditions (kJ/kg)

hCompin,real is the enthalpy of refrigerant entering the compressor at real conditions which corresponds to SH and Te at real conditions (kJ/kg)

Qpipe is the heat loss through the pipe (W).

For example, if the required capacity is 8 kW, the rps range is 30 to 36 rps based on Table 3.

  1. Calculate the rps that meets the need by interpolation. In the above example, the resulting rps is 34.5 rps.

  2. If the calculated rps is lower than the minimum rps (e.g. 18rps), go to Step 2h.5; otherwise, skip Step 2h.5 and go to Step 2h.6.

Step 2h.5: Modify evaporating temperature to further reduce outdoor unit capacity

If the calculated rps is lower than the minimum rps (e.g. 18 rps), it means that the zone heating load (indoor unit side) is so low that it leads to an evaporative capacity (outdoor unit side) which is even lower than the system evaporative capacity corresponding to the minimum compressor speed. In this situation, the evaporating temperature Te as well as the superheating degree SH is modified to further reduce the outdoor unit capacity. More specifically:

  1. Set rps at its minimum value (e.g., 18 rps).

  2. Update Te to meet the required evaporative capacity, using equations described in Step 2h.4a.

  3. SH for each indoor unit can be updated using the equations shown in Step 1.2.

It should be noted that, different from the corresponding step in the cooling mode (Step 2c.5), the Te and SH updates in the heating mode do not affect the refrigerant state and flow rate calculations (as shown in Step 2h.1). Therefore, the piping loss analysis does not need to be repeated here.

Step 2h.6: Calculate the compressor power

  1. Calculate the compressor power by the following procedures.
<!-- -->
  1. Calculate the compressor power at a variety of Loading Index using the following equation. The resulting table (Table [table:atmospheric-variables-at-two-different]) from the same example is shown below.

Mcomp=c1+c2Tc+c3Te+c4Tc2+c5TcTe+c6T2e

Ncomp,rps=Mcomp×Ncomp,ref

where:

c1, …, c6 empirical coefficients corresponding to rps

Te is the effective evaporating temperature (C)

Te,ref is the reference evaporating temperature (C)

Tc is the discharge saturated temperature at the compressor outlet (C)

Tc,ref is the reference condensing temperature (C)

Mcomp is the multiplier for the compressor power calculation (–)

Ncomp,ref is the rated compressor power (W)

Ncomp,rps is the compressor power corresponding to rps (W).

Table [table:atmospheric-variables-at-two-different] – Outdoor unit compressor power at different Loading Index

  1. According to the rps range determined, calculate the compressor power Ncomp by interpolation. In the above example, the compressor power is 1.155 kW.
<!-- -->
  1. Compare the calculated Ncomp above with the initialized Ncomp in Step 2h.2:
  • If NcompNcomp>δ then go to Step 2h.2 for a new round of iteration.

  • Else, end the iteration and go to Step 2h.7.

Step 2h.7: Total power consumption of the outdoor unit

Same as that in the cooling mode (Step 2c.7)

Modeling of the indoor unit (I/U) - Part II - Cooling Mode[LINK]

Step 3c.1: Update air flow rate for each indoor unit

  1. Calculate coil surface temperature for each indoor unit:

Tfs=Te+[AcSH2+BcSH+Cc]

where:

Te is the evaporating temperature decided in the outdoor unit calculations (C)

Tfs is the coil surface temperature (C)

SH is the superheating degrees decided in the outdoor unit calculations (C).

  1. Calculate the enthalpy of the air at the coil surface, Hfs:

ifTfs<T98%,Win, thenHfs=f(Tfs,98%)

ifTfsT98%,Win, thenHfs=f(Tfs,Win)

where:

T98%,Win is the dew point temperature of the indoor air (C)

Hfs is the enthalpy of the air at the coil surface (kJ/kg).

  1. Calculate the required air flow rate Ga for each indoor unit:

Ga=Qin,total/[(Hcoil,inHfs)×(1BF)×ρin]

where:

Qin,total is the total cooling load for the zone (W)

Hcoil,in is the enthalpy of the entering air of the indoor unit (kJ/kg)

ρin is the density of indoor air,fρ(Tin,Win) (kg/m3).

  1. Decide whether to modify SH for further indoor unit capacity reduction.
  • If Ga < Ga,min (e.g., 0.7× Ga,rate ) go to Step 3c.2

  • Else, directly go to Step 3c.3

Step 3c.2: Modify SH to adjust the indoor unit capacity

Set Ga at its minimum value (e.g., 0.7× Ga,rate).

The required coil surface temperature Tfs can be determined as described in Step 1.2.

Given Tfs and Te, SH can be determined using the equation shown in Step 1.3.

  • If SH > 15C, Set SH as 15C. On/Off control strategy may be implemented when needed.

Step 3c.3: Total power consumption of the indoor unit

The power consumption of the indoor unit comes from the fan operations. This can be calculated using the existing VAV fan model in EnergyPlus. Please refer to the current EnergyPlus Engineering Reference for more details.

Modeling of the indoor unit (I/U) - Part II - Heating Mode[LINK]

Step 3h.1: Update air flow rate for each indoor unit.

  1. Calculate coil surface temperature for all the indoor units:

Tfs=Tc[AhSC2+BhSC+Ch]

where:

Tc is the condensing temperature decided in the outdoor unit calculations (C)

Tfs is the coil surface temperature (C)

SC is the subcooling degrees decided in the outdoor unit calculations (C).

  1. Calculate the required air flow rate Ga for each indoor unit:

Ga=Qin,total/[(TfsTcoil,in)×(1BF)×ρin]

where:

Qin,total is the total heating load for the zone (W)

Tcoil,in is the temperature of the entering air of the indoor unit (C)

ρin is the density of indoor air,fρ(Tin,Win) (kg/m3).

  1. Decide whether to modify SH for further indoor unit capacity reduction
  • If Ga < Ga,min (e.g., 0.7× Ga,rate ) go to Step 3h.2

  • Else, directly go to Step 3h.3

Step 3h.2: Modify SC to modify the indoor unit capacity

Set Ga at its minimum value (0.7× Ga,rate).

The required coil surface temperature Tfs can be determined as described in Step 2.

Given Tfs and Te, SC can be determined using the equation shown in Step 1.3.

  • If SC > 20C, Set SC as 20C. On/Off control strategy may be implemented when needed.

Step 3h.3: Total power consumption of the indoor unit

Calculate electric power consumption by the indoor unit using the existing VAV fan model in EnergyPlus. Please refer to the current EnergyPlus Engineering Reference for more details.

Additional energy consumption by defrost and crankcase heater[LINK]

There may be additional energy consumption due to the defrost operation and crankcase heater operation. These components have no impact on the heat pump operations. The calculation methods in the VRF-FluidTCtrl model are the same as those in VRF-SysCurve model. Please refer to the VRF-SysCurve section for more details.

Physics based VRF Heat Recovery Model (VRF-FluidTCtrl-HR)[LINK]

Overview[LINK]

Introduction of VRF-HR Systems[LINK]

This section introduces the model to simulate the energy performance of Variable Refrigerant Flow (VRF) systems with heat recovery (HR) configurations, which is capable of achieving heat recovery from cooling zones to heating zones and providing simultaneous zone cooling and heating operations. The model is developed for the 3-pipe VRF-HR systems, which are the dominant system configuration in the current VRF-HR market.

Compared with the VRF-HP system, VRF-HR is more complicated in terms of system configuration and operational controls. Figure 15 shows the schematic chart of a typical 3-pipe VRF-HR system. As can be seen in the figure, it has dedicated refrigerant pipes for suction gas, liquid and discharge gas. To enable simultaneous cooling and heating, complex refrigerant management loop and more system components are implemented, including one additional heat exchanger in the outdoor unit (OU) and multiple Branch Selector (BS) Units. The two heat exchangers in the outdoor unit can work at different evaporator/condenser combinations to enable specific operational modes to handle diverse and changing indoor heating/cooling load requirements. The Four-Way Directional Valves (FWV) and BS units enable the system to provide separate refrigerant piping connections for different operational modes. This leads to varying refrigerant flow directions and different control logics for various operational modes, and therefore specific algorithm is needed for different major operational modes.

Schematic chart of a 3-pipe VRF-HR system [fig:VRF-HR-Chart-Schematic]

Schematic chart of a 3-pipe VRF-HR system [fig:VRF-HR-Chart-Schematic]

VRF-HR operational modes[LINK]

Depending on the indoor cooling/heating requirements and the outdoor unit operational states, the operations of the VRF-HR system can be categorized into six modes:

  • Mode 1: Cooling load only. No heating load. Both OU heat exchangers operate as condensers.

  • Mode 2: Simultaneous heating and cooling. The sum of the zone cooling loads and compressor heat is much larger than the sum of the zone heating loads. Both OU heat exchangers operate as condensers.

  • Mode 3: Simultaneous heating and cooling. The sum of the zone cooling loads and compressor heat is slightly larger than the sum of the zone heating loads. One OU heat exchanger operates as a condenser while the other as an evaporator.

  • Mode 4: Simultaneous heating and cooling. The sum of the zone cooling loads and compressor heat is slightly smaller than the sum of the zone heating loads. One OU heat exchanger operates as a condenser while the other as an evaporator.

  • Mode 5: Simultaneous heating and cooling. The sum of the zone cooling loads and compressor heat is much smaller than the sum of the zone heating loads. Both OU heat exchangers operate as evaporators.

  • Mode 6: Heating load only. No cooling load. Both OU heat exchangers operate as evaporators.

System-level Heat Balance Diagram for the Six VRF-HR Operational Modes [fig:VRF-HR-Chart-HeatBalance]

System-level Heat Balance Diagram for the Six VRF-HR Operational Modes [fig:VRF-HR-Chart-HeatBalance]

The system-level heat balance diagram for all the six operational modes are shown in Figure 16. Note that the heat recovery loss (HR loss) only exists in Mode 3 and Mode 4, in which the OU evaporator and condenser run simultaneously. In these modes, the following two items are at similar levels: (a) sum of IU heating loads and IU condenser side piping loss, and (b) sum of IU cooling loads, IU evaporator side piping loss and heat released by the compressor. Take Mode 3 for example, item (b) is higher than item (a), and therefore the system requires the operation of OU condenser to release the extra heat to ensure the system-level heat balance. However, the extra heat is at a relatively low level so that the system needs to release more heat than required via OU condenser and meanwhile runs OU evaporator to ensure the heat balance as well as system reliability. This leads to the presence of HR loss.

System operation at various modes[LINK]

With the help of FWV and BS units, every operational mode has its own refrigerant piping connections to achieve different refrigerant flow directions, as shown in Figure 17. This leads to different refrigerant operations and piping loss situations, as shown in the Pressure-Enthalpy Diagrams in Figure 18.

Take Mode 3 for instance, the system utilizes one OU heat exchanger as an evaporator and the other as a condenser. The superheated refrigerant coming out of the compressor is divided into two parts. One part passes through the OU condenser within which the refrigerant is condensed and subcooled (segment 2-4’), while the other part passes through the IU condenser to satisfy the indoor heating requirements (segment 3-4). Note that before entering the IU condenser, the second part of the refrigerant flow passes through the main discharge gas pipe connecting the indoor and outdoor units, where the refrigerant pressure drop and heat loss occurs. This leads to extra heating load and a lower condensing temperature level at IU condenser than that at OU condenser. The two parts of refrigerant leaving condensers are then expanded through expansion valves following segment 4“-5/5’ and pass through the IU evaporator (segment 5-6) and OU evaporator (segment 5’-1’). Similarly to the refrigerant passing through the main discharge gas pipe, the refrigerant leaving IU evaporators passes through the main suction gas pipe where the refrigerant pressure drop and heat loss occurs (segment 6-1”). This leads to extra cooling load and a higher evaporating temperature level at IU evaporator than that at OU evaporator. Finally, the refrigerant flow leaving the OU evaporator and that leaving the main discharge gas pipe are combined (segment 1’/1“-1) and enter the compressors for the next cycle. In this regard, the system provides simultaneous heating and cooling to the indoor environment.

A number of operational parameters are controlled to ensure system heat balance and stable and reliable system operation, including evaporating temperature levels, condensing temperature levels, superheating and sub-cooling degrees, and the refrigerant flow rates at various components. The operational control logics for various modes are different and therefore particular algorithm needs to be designed for different operational modes in the new VRF-HR model.

Piping diagrams for the six VRF-HR operational modes [fig:VRF-HR-Chart-Piping]

Piping diagrams for the six VRF-HR operational modes [fig:VRF-HR-Chart-Piping]

Illustration of the six VRF-HR operational modes on Pressure-Enthalpy diagrams [fig:VRF-HR-Chart-EnthalpyPressure]

Illustration of the six VRF-HR operational modes on Pressure-Enthalpy diagrams [fig:VRF-HR-Chart-EnthalpyPressure]

Model Description[LINK]

The algorithm of the VRF-HR Cooling Only Mode and Heating Only Mode are the same as those in the VRF-HP Cooling Mode (Figure 14) and Heating Mode (Figure 13), respectively. The holistic logics of the VRF-HR Simultaneous Heating and Cooling Mode are illustrated in Figure 19. In general, the effective evaporating temperature and/or condensing temperature for the indoor units is first determined based on the load requirements and indoor unit configurations. Then the pressure and heat losses through the main pipe are calculated at the given operational conditions. After that, the effective condensing temperature and/or evaporating temperature of the outdoor unit is calculated, taking into account the piping loss. With the above calculations, the compressor speed can be determined using the performance curves describing the evaporative capacity of the system at specific condensing/evaporating temperature combinations at various compressor speeds. Finally, the compressor power at that speed is calculated and the total electric power consumption by the outdoor unit is obtained. Note that particular control logics are designed for various operational modes, and therefore lead to different refrigerant temperature adjustments. Also, note that a number of calculation steps are coupled together, for instance, the piping loss calculation and the system performance calculation. More specifically, the piping loss changes the operating conditions of the system, which may lead to a different control strategy and thus affect the amount of piping loss. This makes it difficult to obtain an analytical solution for a number of operational parameters (e.g., enthalpy of refrigerant entering the indoor unit), and therefore numerical iterations are employed to address this problem.

Schematic chart of the new VRF-HR algorithm: Simultaneous Heating and Cooling Mode [fig:VRF-HR-AlgorithmOverview-Mode2-5]

Schematic chart of the new VRF-HR algorithm: Simultaneous Heating and Cooling Mode [fig:VRF-HR-AlgorithmOverview-Mode2-5]

Detailed calculation procedures for the VRF-HR Simultaneous Heating and Cooling Mode are summarized as follows.

Step 1: Obtaining zonal load/condition information

The following information for each zone from the zone modules within EnergyPlus are obtained:

  • zone sensible loads Qin,sensible

  • zone total loads Qin,total

  • indoor air temperature Tin

  • indoor air humidity ratio Win

If there are only zone cooling loads and no zone heating loads, go to the VRF-HR Cooling Only Mode, the algorithms of which is the same as those for the VRF-HP Cooling Mode. If there are only zone heating loads and no zone cooling loads, go to the VRF-HR Heating Only Mode, the algorithms of which is the same as those for the VRF-HP Heating Mode. Otherwise, go to the VRF-HR Simultaneous Heating and Cooling Mode as described below.

Step 2: Calculate IU required evaporating temperature and/or the condensing temperature Evaluate the required coil surface air temperature Tfs and then the required evaporator refrigerant temperature Te,req for each indoor unit with cooling requirements. Likewise, evaluate the required condenser refrigerant temperature Tc,req for each indoor unit with heating requirements. (Refer to Step 1.2 in the VRF-FluidTCtrl-HP model for more details.)

Step 3: Calculate IU effective evaporating temperature Te and/or the condensing temperature Tc There are two refrigerant temperature control strategies for the indoor unit: (1) ConstantTemp, (2) VariableTemp.

  • In the ConstantTemp strategy, Te and Tc are kept at constant values provided by the user.

  • In the VariableTemp strategy, Te and Tc are determined using the required evaporating/condensing temperature calculated in Step 2.

(Refer to Step 1.3 in the VRF-FluidTCtrl-HP model for more details.)

Step 4: IU condenser side piping loss calculations

This section calculates the IU condenser side piping loss, which occurs at the High and Low Pressure Gas Pipe where the refrigerant flowing from the OU compressor outlets to the IU condensers. It includes both the refrigerant pressure drop ΔPpipe and heat loss Qpipe. Note that the change of compressor operational conditions may lead to different control strategies of the system, which in reverse affects the amount of piping loss. So the piping loss analysis and system performance analysis are coupled together. Numerical iterations are designed to address the coupling effect.

At this step, the compressor discharge saturated temperature Tc (i.e., saturated vapor temperature corresponding to compressor discharge pressure) can be obtained using the calculated refrigerant pressure drop ΔPpipe.

(Refer to Step 2h.1 in the VRF-FluidTCtrl-HP model for more details.)

Step 5: IU evaporator side piping loss calculations

This section calculates the IU evaporator side piping loss, which occurs at Suction Gas Pipe where the refrigerant flowing from the IU evaporators to the OU compressor inlets. Similarly to the IU condenser side piping loss, it includes both the refrigerant pressure drop ΔPpipe and heat loss Qpipe.

At this step, the compressor suction saturated temperature Te (i.e., saturated vapor temperature corresponding to compressor suction pressure) can be obtained using the calculated refrigerant pressure drop ΔPpipe.

Note that one key input of the IU evaporator side piping loss calculations is the enthalpy of the refrigerant at IU evaporator inlets. It is assumed to be equal to the average enthalpy of the refrigerant at IU condenser outlets, which is obtained in the IU condenser side piping loss calculations.

(Refer to Step 2c.1 in the VRF-FluidTCtrl-HP model for more details.)

Step 6: Determine the operational mode for simultaneous heating and cooling operations

As noted earlier, simultaneous heating and cooling operations include the following modes:

  • Mode 2: Cooling dominant w/o HR loss

  • Mode 3: Cooling dominant w/ HR loss

  • Mode 4: Heating dominant w/ HR loss

  • Mode 5: Heating dominant w/o HR loss

This section determines the operational mode based on the load requirements and operational conditions:

a. Calculate the Loading Index LI1 satisfying IU cooling load (Refer to Step 2c.4 in the VRF-FluidTCtrl-HP model for more details.)

b. Calculate the Loading Index LI2 satisfying IU heating load (Refer to Step 2h.4 in the VRF-FluidTCtrl-HP model for more details.)

c. If LI1 LI2, the system operates at Mode 5

d. If LI1 > LI2 and Te < To - Tdiff, the system operates at Mode 2 (To: outdoor air dry-bulb temperature; Tdiff: a constant value representing the difference between OU evaporating temperature and outdoor air temperature during simultaneous heating and cooling)

e. If LI1 > LI2 and Te To - Tdiff, the system operates at Mode 3 or 4 (these two modes can be handled by one set of algorithms)

Step 7: OU operation analysis at Mode 5

If Te < To - Tdiff, perform the following procedures:

a. Select the compressor speed corresponding to LI2

b. Calculate the compressor power corresponding to LI2 and the previously obtained Tc and Te (Refer to Step 2c.6 in the VRF-FluidTCtrl-HP model for more details.)

c. Calculate the evaporative capacity Captot,evap provided by the compressor at LI2 and the previously obtained Tc and Te (Refer to Step 2c.4 in the VRF-FluidTCtrl-HP model for more details.)

d. Calculate the OU evaporator load Capou,evap based on system-level heat balance

e. Obtain the OU fan flow rate mairevap corresponding to Capou,evap, and thus the fan power (Refer to Step 2c.3 in the VRF-FluidTCtrl-HP model for more details.)

If Te To - Tdiff, perform the following procedures:

a. Select the compressor speed corresponding to LI1

b. Perform iterations between step b-i to identify the compressor Loading Index and power consumption.

c. Initialized compressor power Ncompini

c.1 For the 1st iteration step, calculate Ncomp = fpow,comp( Tc, To – 5, LI2)

c.2 For the following iteration steps, update Ncomp = (Ncompini + Ncompnew)/2

d. Calculate the OU evaporator load Capou,evap based on system-level heat balance

e. Obtain the OU evaporating temperature Te level using Capou,evap and the rated air flow rate (Refer to Step 2c.3 in the VRF-FluidTCtrl-HP model for more details.)

f. Update Te level and IU evaporator side piping loss, corresponding to Te update

g. Identify the compressor Loading Index LInew to provide sufficient evaporative capacity Captot,evap at updated Te level (Refer to Step 2c.4 in the VRF-FluidTCtrl-HP model for more details.)

h. Calculate the compressor power Ncompnew corresponding to LInew and the updated Te (Refer to Step 2c.6 in the VRF-FluidTCtrl-HP model for more details.)

i. Compare Ncompnew and Ncompini. Start a new round of iteration if the difference is greater than the tolerance.

Step 8: OU operation analysis at Mode 2

a. Select the compressor speed corresponding to LI1

b. Calculate the compressor power corresponding to LI1 and the previously obtained Tc and Te (Refer to Step 2c.6 in the VRF-FluidTCtrl-HP model for more details.)

c. Calculate the evaporative capacity Captot,evap provided by the compressor at LI1 and the previously obtained Tc and Te (Refer to Step 2c.4 in the VRF-FluidTCtrl-HP model for more details.)

d. Calculate the OU condenser load Capou,cond based on system-level heat balance

e. Obtain the OU fan flow rate mair,cond corresponding to Capou,cond, and thus the fan power

Step 9: OU operation analysis at Mode 3 or 4

a. Select the compressor speed corresponding to LI1

b. Perform iterations between step b-e to identify the updated Te level within the range of To-5 and the original Te.

c. Calculate the evaporative capacity Captot,evap provided by the compressor at LI1 and the previously obtained Tc and assumed Te (Refer to Step 2c.4 in the VRF-FluidTCtrl-HP model for more details.)

d. Calculate the OU evaporator load Captot,evap at assumed Te level and rated fan flow rate (Refer to Step 2h.3 in the VRF-FluidTCtrl-HP model for more details.)

e. Perform iterations to identify the updated Te level to ensure the heat balance for the b. and c. calculations

f. Update Te level and IU evaporator side piping loss, corresponding to Te update

g. Calculate the compressor power corresponding to LI1 and the updated Te (Refer to Step 2c.6 in the VRF-FluidTCtrl-HP model for more details.)

h. Calculate the OU condenser load Capou,cond based on system-level heat balance

i. Obtain the OU fan flow rate mair,cond corresponding to Capou,cond, and thus the fan power

Step 10: Modify IU operational parameters for capacity adjustments

The air flow rate and SH/SC value of each indoor unit can be manipulated to adjust the cooling/heating capacity. (Refer to Step 3c and 3h in the VRF-FluidTCtrl-HP model for more details.)

Zone Terminal Unit List[LINK]

The zone terminal unit list identifies the terminal units that are connected to a single variable refrigerant flow heat pump. The zone terminal unit list is used exclusively in the variable refrigerant flow (VRF) heat pump object (ref: AirConditioner:VariableRefrigerantFlow) and VRF zone terminal units (ref: ZoneHVAC: TerminalUnit:VariableRefrigerantFlow). Up to 20 terminal units may be connected to a single VRF outdoor condensing unit. This list is extensible if additional indoor terminal units are required. The following figure shows the connection scheme between the zone terminal units, the zone terminal unit list, and finally the VRF AC system. The zone terminal units are connected to the zone through zone inlet and outlet zone nodes. Each zone terminal unit is entered in a list which represents all terminal units connected to a single VRF AC system. And finally, the zone terminal unit list name is entered in the corresponding VRF AC object.

Zone Terminal List connections in EnergyPlus objects [fig:zone-terminal-list-connections-in-energyplus]

Zone Terminal List connections in EnergyPlus objects [fig:zone-terminal-list-connections-in-energyplus]

References[LINK]

R. Raustad. A variable refrigerant flow heat pump computer model in EnergyPlus, ASHRAE Transactions (2013), 119 (1):1-9.

T. Hong, K. Sun, R. Zhang, R. Hinokuma, S. Kasahara, Y. Yura. Development and validation of a new variable refrigerant flow system model in EnergyPlus, Energy Build. 117 (2015): 399–411.

T. Hong, R. Zhang, K. Sun. Development, Validation and Performance Evaluation of VRF Heat Recovery Systems. LBNL Report, 2016.