Engineering Reference — EnergyPlus 8.2

<< Prev | Table of Contents | Next >>

Engineering Reference

Refrigeration Equipment[LINK]

Overview[LINK]

EnergyPlus can model refrigerated case equipment consisting of a compressor rack, multiple refrigerated cases and walk-in coolers, secondary loop equipment, and optional heat reclaim air and water heating coils. The refrigerated case equipment models perform four major functions:

  • calculate the electric consumption of refrigerated cases and walk-in coolers connected to a compressor rack

  • determine the impact of refrigerated cases and walk-in coolers on zone cooling and dehumidification loads (i.e., case credits), including the effects of HVAC duct configuration

  • calculate the electric consumption and COP of the compressor rack, and the electric and water (if applicable) consumption related to cooling the compressor rack’s condenser.

  • determine the total amount of heat rejected by the compressor rack’s condenser and store this information for use by waste heat recovery models (e.g., using Desuperheater heating coil (object: Coil:Heating:Desuperheater) as an air reheat coil for high humidity control in a supermarket)

The case and walk-in models account for nearly all performance aspects of typical supermarket refrigeration equipment. Refrigerated case and walk-in performance are based on the combined effects of evaporator load, fan operation, lighting, defrost type, and anti-sweat heater operation. Optional air and water heating coils can be modeled to reclaim available waste heat (superheat) from the compressor rack.

The user has two options when describing the balance of the system. Energy used to cool the condenser is simulated in both approaches. The simplest option is to use a compressor rack object, combining the compressors and condenser into a single unit with the performance determined by the heat rejection environment and the total case load. An example schematic of a compressor rack system is shown in Figure 276 below.

A detailed refrigeration system object models compressor and condenser performance separately. The detailed refrigeration system also includes the ability to transfer refrigeration load from one system to another using subcoolers, cascade condensers, and secondary loops. An example schematic of the detailed refrigeration system is shown in Figure 277 below. Subcooler #2 is shown twice on Figure 277 because it represents a liquid suction heat exchanger. This type of subcooler uses the cool suction gas to subcool the warmer condensed liquid. Subcoolers #1 and #3 on Figure 277 represent mechanical subcoolers. These subcoolers are used to subcool the condensate on a lower-temperature system using the cold liquid refrigerant from a higher temperature system. On this example, only subcoolers #1 and #2 would be defined as a part of the refrigeration system. However, subcooler #3 would place a refrigerating load, similar to the load of a refrigerated case, on the system.

Typical Compressor Rack Equipment Schematic

Typical Compressor Rack Equipment Schematic

Typical Detailed Refrigeration System Equipment Schematic

Typical Detailed Refrigeration System Equipment Schematic

Four classes of secondary refrigeration loops can be modeled:

  • a separate water loop is used to remove heat rejected by the condenser,

  • a lower-temperature refrigeration system rejects heat to a higher-temperature refrigeration system via a cascade condenser,

  • a fluid, such as a brine or glycol solution, is cooled in a secondary evaporator and is then circulated to chill the refrigerated cases and walk-ins, and

  • a refrigerant, such as CO2, is partially evaporated in the refrigerated cases and walk-ins in a liquid-overfeed circuit, and then condensed in a secondary evaporator.

The first two classes of secondary loops are modeled using Refrigeration:System objects with Refrigeration:Condenser:WaterCooled and Refrigeration:Condenser:Cascade objects, respectively. Figure 277 shows how cascade condensers and secondary evaporators are treated as a refrigeration load on a primary detailed system. The second two classes are modeled with a Refrigeration:SecondarySystem object described later in this section.

The compressor rack, detailed and secondary refrigeration systems, refrigerated case, and other component models are described below. The optional air and water heating coils are described elsewhere in this document (Ref. objects Coil:Heating:Desuperheater and Coil:WaterHeating:Desuperheater).

Refrigeration Compressor Racks[LINK]

The refrigerated case compressor rack object works in conjunction with the refrigerated case and walk-in cooler objects (Refrigeration:Case and Refrigeration:WalkIn) to simulate the performance of a simple supermarket-type refrigeration system. This object (Refrigeration:CompressorRack) models the electric consumption of the rack compressors and the cooling of the compressor rack’s condenser. Heat removed from the refrigerated cases and walk-ins and compressor/condenser fan heat can be rejected either outdoors or to a zone. Compressor rack condenser waste heat can also be reclaimed for use by an optional air heating coil (Ref. object Coil:Heating:Desuperheater) or by a user-defined plant water loop (Ref. object Coil:WaterHeating:Desuperheater).

The performance of the compressor rack is simulated using the sum of the evaporator loads for all refrigerated cases and walk-ins connected to the rack. Whether a single refrigerated case is connected to a rack (e.g., stand-alone refrigerated case, meat cooler, or produce cooler) or several cases are connected to a rack, the rack electric consumption is calculated based on the total evaporator load for the connected cases and walk-ins and the coefficient of performance (COP) for the compressor rack. At least one refrigerated case or walk-in must be connected to the compressor rack. The model assumes the compressor rack has sufficient capacity to meet the connected refrigeration load for any simulation time step. Additionally, the model neglects compressor cycling losses at part-load conditions.

For condenser heat rejection to the outdoors, condenser cooling can be modeled as dry air cooling, wet evaporative cooling, or water loop cooling. Using evaporative cooling rather than dry air cooling will allow for more efficient condenser heat rejection based on the entering air approaching the wet-bulb temperature rather than the dry-bulb temperature. Analyses under the International Energy Agency’s (IEA) Heat Pumping Programme Annex 26 indicates that this measure can improve refrigeration system efficiency by up to 10% (IEA 2003). The use of an evaporative-cooled condenser requires a water pump and, optionally, a basin sump water heater (to protect against freezing). Makeup water will also be required to replace that lost by evaporation. In colder climates, some evaporative-cooled condensers are drained for the winter months and run as dry air units. This scenario can be modeled by using an optional evaporative condenser availability schedule.

The simulation of the evaporative cooled condenser utilizes an effective air dry-bulb temperature that is assumed to be the result of evaporation of water in the air stream (similar to object EvaporativeCooler:Direct:CelDekPad). As discussed below, this effective temperature is used by performance curves that are a function of temperature. While some designs of evaporative coolers use water film cascading across the condenser coil for evaporative cooling, the current model uses the effective temperature method as a surrogate for the more complex water film on coil calculations.

If the condenser heat rejection is specified as water cooled, an appropriate plant water loop must be defined by the user (see documentation on Plant/Condenser Loops for additional details about plant loops). This will include defining cooling supply components, such as pumps, water storage tanks, and cooling towers, as well as related branches, nodes, and connectors. The heat rejection from the refrigeration condenser is modeled as a cooling demand, which is satisfied by heat extraction devices (e.g., water tank and cooling tower) on the cooling supply side of a water loop. An example of such an arrangement is shown in Figure 278.

Example Of Condenser Heat Recovery To Water Storage Tank

Example Of Condenser Heat Recovery To Water Storage Tank

Compressor Energy Use[LINK]

Calculation of compressor rack electric power uses a simple model based on the total evaporator load (sum of the evaporator loads for all refrigerated cases and walk-ins connected to a rack) and the compressor rack operating COP which accounts for the air temperature entering the condenser:

COPoperating=COPdesign(COPfTemp)

where:

COPoperating = compressor coefficient of performance at actual operating conditions (W/W)

COPdesign = compressor coefficient of performance at design conditions (W/W)

COPfTemp = output of the normalized “Compressor Rack COP as a Function of Temperature Curve” (dimensionless)

Because the COP curve is defined only as a function of the condensing temperature, it is important that this curve definition corresponds to the lowest evaporating temperature served by the compressor rack. The air temperature used to evaluate the “Compressor Rack COP as a Function of Temperature Curve” depends on where the compressor rack’s condenser is located (Heat Rejection Location). When modeling condenser heat rejected directly to a zone (typical of a stand-alone packaged refrigerated case with integral condenser located in a building zone), the zone air dry-bulb temperature is used to calculate the change in compressor COP from the design value. If more than one refrigerated case and no walk-ins are attached to a compressor rack that rejects its condenser heat to a zone, then all cases served by this rack must reside in the same zone. When modeling a compressor rack serving at least one walk-in, OR with condenser heat rejected to outdoors, the refrigerated cases and walk-ins connected to this rack may be located in different zones. If the condenser type is specified as “Air Cooled”, the outdoor air dry-bulb temperature is used to evaluate the “Compressor Rack COP as a Function of Temperature Curve.” If the condenser type is specified as “Evap Cooled”, the air temperature leaving the condenser is related to the effectiveness of the evaporative cooling system. If the evaporative process were 100% effective, the effective temperature of air leaving the evaporative media would equal the air wet-bulb temperature. However, the efficiency of the direct evaporative process is typically less than 100%, and the effective temperature leaving the condenser is determined by:

Teffective=Towb+(1ε)[TodbTowb]

where:

Teffective = effective dry-bulb temperature of air leaving the condenser cooling coil (°C)

Towb = outdoor air wet-bulb temperature (°C)

Todb = outdoor air dry-bulb temperature (°C)

ε = evaporative condenser effectiveness.

If the user is modeling an evaporative cooled condenser and is using COPfTemp curve data (e.g., manufacturer’s data) based on wet-bulb temperature rather than dry-bulb temperature, the evaporative condenser effectiveness should be set to 1.0 for consistency.

If the condenser is water cooled, the effective temperature experienced by the condenser is based on the return water temperature from the plant loop heat rejection system (e.g., cooling tower) that is defined by the user. This return water temperature is typically related to the outdoor ambient conditions at each time step.

The electric power input to the rack compressor(s) is calculated for each simulation time step as the sum of the connected refrigerated case evaporator loads divided by the operating COP:

Prack=˙Qcase+˙QwalkinCOPoperating

where:

Prack = output variable “Refrigeration Compressor Rack Electric Power [W]”, electric power input to the rack compressor(s)

˙Qcase = evaporator load for each refrigerated case connected to the rack (W)

˙Qwalkin = refrigeration load for each walk-in connected to the rack (W)

Condenser Heat Rejection, Energy Use, and Water Use[LINK]

The compressor rack can reject heat to an air-, water-, or evaporative-cooled condenser. The condenser type determines the heat rejection temperature used for the compressor rack COP calculation. The compressor rack also allows superheat heat reclaim and heat rejection to a conditioned zone.

Condenser Fan Energy Use[LINK]

Condenser fan power for any simulation time step is calculated by multiplying the design fan power by the condenser fan power as a function of temperature curve.

PCondFan=PCondFan,design(CondFanfTemp)

where:

PCondFan = output variable “Refrigeration Compressor Rack Condenser Fan Electric Energy [W]”

PCondFan,design = design condenser fan power (W)

CondFanfTemp = output of the optional “Condenser Fan Power as a Function of Temperature Curve”

Similar to the compressor rack energy use described above, the air temperature used to evaluate the “Condenser Fan Power as a Function of Temperature Curve” depends on where the condenser rack’s condenser is located (i.e., zone air dry-bulb temperature if the condenser is located in a zone, outdoor air dry-bulb temperature if the condenser is located outdoors and is specified as air cooled, or effective temperature if the condenser is outdoors and is specified as evaporative cooled). If the sum of the evaporator loads for the refrigerated cases connected to the rack is equal to zero, the condenser fan power is set equal to zero. If the user does not provide a “Condenser Fan Power as a Function of Temperature Curve”, then the model assumes the condenser fan power is at the design power level when any of the refrigerated cases connected to this rack are operating.

If the user is modeling an evaporative cooled condenser and is using CondFanfTemp curve data based on wet-bulb temperature rather than dry-bulb temperature, the evaporative condenser effectiveness should be set to 1.0 for consistency.

For a water cooled condenser, there is no fan load at the condenser (i.e., the water/refrigerant heat exchanger). Any fan load would be related to and accounted for at the heat rejection object (e.g., cooling tower).

Superheat Reclaim Heating Coil[LINK]

EnergyPlus can simulate waste heat being reclaimed from a compressor rack for use by a refrigerant-to-air or refrigerant to water heating coil. Heat reclaimed from the compressor rack is assumed to be recovered from the superheated refrigerant gas leaving the compressor(s) and does not directly impact the performance of the compressor rack or refrigerated cases connected to the rack. The total heat rejected by the condenser (in Watts) is calculated each time step as follows:

˙Qcondenser=(˙Qcase+˙Qwalkin)(1+1COPoperating)

The heat reclaim heating coil is able to transfer a fixed percentage of this total amount of rejected energy (not to exceed 30%) and use it to heat air and water. Refer to objects Coil:Heating:Desuperheater and Coil:WaterHeating:Desuperheater for a complete description of how these coils are modeled.

NOTE: When modeling a heat reclaim coil, the heat rejection location in the Refrigeration:CompressorRack object must be “Outdoors”. If the compressor rack heat rejection location is “Zone”, the total amount of waste heat available for reclaim (e.g., by a desuperheater heating coil) is set to zero by the compressor rack object and the simulation proceeds.

Heat Rejection to Zone[LINK]

The compressor rack model can simulate condenser heat being rejected to a zone. As explained previously, if this heat rejection option is selected then all refrigerated cases connected to the rack must be located in the same zone and a superheat heat reclaim heating coil can not be modeled (Ref. Superheat Reclaim Heating Coil).

The refrigerated case and walk-in objects (Refrigeration:Case and Refrigeration:WalkIn) already calculate and report the sensible case credits which impact the zone air heat balance (Ref. Sensible Case Credits). When refrigerated cases and/or walk-ins are served by a compressor rack that rejects condenser waste heat directly to the zone (e.g., a stand-alone refrigerated case with integral compressor and condenser), this condenser waste heat also impacts the zone air heat balance and offsets some or all of the sensible case credits.

If only cases are served, the amount of condenser waste heat rejected to the zone and/or the HVAC return air (zone return air path outlet node) is calculated and reported by the refrigerated case compressor rack object as follows:

˙QZone,heating=(˙Qcase[1RAF])(˙Qcase)(˙Qcondenser+PCondFan)

˙QHVAC,heating=(˙Qcondenser+PCondFan)˙Qzone,heating

where:

˙QZone,heating = output variable “Refrigeration Compressor Rack Zone Sensible Heating Rate [W]”

RAF= return air factor for each case connected to the rack (Ref. Figure 279)

˙QHVAC,heating = output variable “Refrigeration Compressor Rack Return Air Sensible Heating Rate [W]”

If the HVAC system is off for a simulation time step (no return air mass flow), the rack condenser heat normally attributed to the HVAC return is set equal to zero and all condenser heat energy is applied to the zone air heat balance.

If, however, walk-in cooler(s) are also served by this compressor rack, no condenser heat is rejected to the HVAC return air. For walk-in cooler(s), the user must specify the zone that accepts the condenser heat rejection (because walk-ins can exchange heat with multiple zones). In that case:

˙QZone,heating=˙Qcondenser+PCondFan

Water Cooled Condenser[LINK]

If the refrigeration condenser is water cooled, a water plant loop must be defined in the input file. At a minimum, the loop must contain a pump and one or more heat sinks of sufficient capacity to remove the condenser heat load. In the system shown in Figure 278, the heat sinks are the water heater tank and the cooling tower. The water pump in the loop can be either constant (Ref. Pump:ConstantSpeed) or variable speed (Ref. Pump:VariableSpeed). A variable speed pump permits the loop flow to vary and allows for a setpoint to be established on the condenser outlet water temperature. As the refrigeration condenser heat load varies through time, the speed of the pump can be adjusted to achieve a mass flow consistent with a desired outlet water temperature according to

m=Qcondensercp(ToutTin)

where:

m = mass flow in the water loop

Qcondenser= heat rejected by the condenser

cp= specific heat of water

Tout= desired water outlet temperature

Tin= return water inlet temperature.

The desired water outlet temperature is specified using a schedule, subject to a maximum water outlet temperature (input specified). The maximum temperature is typically defined by constraints on the refrigerant loop pressures and temperatures. The desired mass flow in the water loop to meet the temperature schedule is also compared to the user-supplied maximum flow rate. If the desired mass flow is greater than the maximum allowed flow, the flow rate is set to the maximum value and the resulting water outlet temperature is determined.

The return water inlet temperature is a function of the cooling system defined by the user. A minimum return water temperature may need to be taken into consideration to prevent lowering the resulting refrigerant condensing pressure to the point that refrigerant expansion valve operation becomes impaired. When ambient conditions produce low temperature warnings based on the minimum return water temperature, an outlet temperature setpoint control may need to be placed on the water heat sink object (e.g., cooling tower) to keep the return water temperature above the minimum.

If the water loop flow is constant (i.e., driven by a constant speed pump), then the outlet water temperature will vary with the amount of heat rejected by the condenser. Using the equation above, the resulting water outlet temperature is calculated as

Tout=Qcondensercpm+Tin

Evaporative Condenser Water Pump[LINK]

If the condenser type is specified as “Evap Cooled”, a water pump is required to circulate water in the evaporative condenser. The pump power can be input directly or be autocalculated using a relationship of 0.004266 W per watt [15 W/ton] of rated total cooling capacity where the total cooling capacity is the sum of the rated total cooling capacities for the refrigeration load connected to this compressor rack. Following manufacturer’s recommendations regarding the avoidance of scaling, the water pump does not cycle when there is no cooling demand (i.e., when the compressors are not running), but rather runs continuously. However, if the evaporative condenser availability schedule is set such that evaporative cooling is not available (e.g., during very cold months to avoid freezing), then the pump power consumption will be zero during that period.

Evaporative Condenser Water Consumption[LINK]

With evaporative cooling of the condenser’s entering air, makeup water is needed to replenish the water lost due to evaporation. The quantity required is calculated as the product of the air mass flow rate and the difference between the entering and leaving air humidity ratio, divided by the density of water. The air mass flow rate is determined by multiplying the evaporative condenser air volume flow rate times the density of the entering air (i.e., at the condenser air inlet node if provided, or outdoor air conditions [e.g., no adjustment for height above ground] if the condenser air inlet node field is left blank). The volumetric air flow rate is either specified directly in the user input or is autocalculated using the relationship 0.000144 m3/s per watt of rated total cooling capacity [850 cfm/ton] where the total cooling capacity is the sum of the rated total cooling capacities for the refrigerated cases and walk-ins connected to this compressor rack (Ref. Refrigeration:Case and Refrigeration:WalkIn). The air mass flow rate is multiplied by the variable CondFanfTemp, described above, to simulate the modulation of air flow by the condenser fans (e.g., staging, multi-speed, or variable speed) as a function of temperature. Mathematically,

˙Vevaporation,makeup=˙mair(CondFanfTemp)(ωair,outletωair,inlet)ρwater

where:

˙Vevaporation,makeup = Refrigeration Compressor Rack Evaporative Condenser Water Volume Flow Rate (m3/s)

˙mair= mass flow rate of air through the evaporative condenser (kg/s)

ωair,outlet = humidity ratio of air leaving the evaporative media (kgwater/kgdry air) based on the effective dry-bulb temperature Teffective, as described above, outdoor air wet-bulb temperature, and outdoor barometric pressure

ωair,inlet = humidity ratio of inlet air (kgwater/kgdry air) based on conditions at the condenser air inlet node if provided, or outdoor air conditions (e.g., no adjustment for height above ground) if the condenser air inlet node field is left blank

ρwater = density of water evaluated at the effective air temperature (kg/m3)

The source of the makeup water may be specified as a water storage tank. If not specified, the makeup water is assumed to come from the building mains (Ref. Water Mains Temperatures).

Evaporative Condenser Basin Heater[LINK]

In cold climates, a basin heater may be needed to prevent freezing of the evaporative cooling water. This feature is included in the model whereby an electric basin heater provides heat to the sump water only when the condenser cooling system is idle (i.e., no refrigeration load) and when the outdoor air dry-bulb temperature is below a user-specified setpoint. Since heat balances and basin water temperatures are not explicitly determined, a linear loading relationship, as a function of the difference in outdoor air dry-bulb temperature and the setpoint temperature, is used calculate the power demand at a given time step by the basin heater.

Pbasinheater=Pheatercapacity(TsetpointTOutDb)

where:

Pbasinheater = electric power demand for basin heater in current time step (W)

Pheatercapacity = electric heater capacity as a function of differential temperature (W/deg K)

Tsetpoint = setpoint temperature below which the heater turns on (°C)

TOutDb = outdoor air dry-bulb temperature (°C)

A default value for the basin heater capacity of 200 W/deg K has been established based on manufacturer data.

Evaporative Condenser Availability Schedule[LINK]

Some manufacturer’s evaporative cooling systems for refrigeration condensers permit seasonal draining in the colder months and operation as an air-cooled system during that time. This optional feature is available through an availability schedule. This is important in climates subject to freezing weather in order to avoid excessive ice formation on the condenser surfaces and surroundings. (The Availability Schedule is the correct way to model the use of evaporative condensers in cold climates. However, some users may take a single input description and use it to model a building with a refrigeration system in a variety of climates. To avoid modeling the use of evaporative coolers in freezing weather, the code includes a cutout to switch to dry operation whenever the outdoor drybulb temperature drops below 4C.) During periods when evaporative cooling is not available, the outdoor condenser behaves as an air-cooled system with no water consumption or pump and basin heater loads. The effective temperature of air entering the condenser coil during this period (used to evaluate COPfTemp and CondFanfTemp) is equal to the outdoor air dry-bulb temperature at the condenser air inlet node if provided, or outdoor air conditions (e.g., no adjustment for height above ground) if the condenser air inlet node field is left blank.

Refrigerated Cases[LINK]

The refrigerated case object (Refrigration:Case) works in conjunction with the compressor rack, detailed refrigeration system, or secondary refrigeration system object (Refrigeration:CompressorRack, Refrigeration:System, or Refrigeration:SecondarySystem) to simulate the performance of a refrigerated case system. The refrigerated case model uses performance information at rated conditions along with performance curves for latent case credits and defrost heat load to determine performance at off-rated conditions. Energy use for lights, fans and anti-sweat heaters is modeled based on inputs for nominal power, schedules, and control type. The refrigerated case model accounts for the sensible and latent heat exchange with the surrounding environment (termed “case credits”) which impacts the temperature and humidity in the zone where the case is located. The simplified model described here provides the flexibility to simulate a broad range of refrigerated case types.

The total load on the refrigerated case evaporator is made up of various components:

˙Qcase=˙Qwalls+˙Qrad+˙Qinf,sens+˙Qinf,lat+˙Qlights+˙Qas+˙Qdef+˙Qfan+˙Qrestock

where:

˙Qcase = total load on the refrigerated case evaporator (W)

˙Qwalls = heat transfer through case walls due to the difference between the refrigerated case operating dry-bulb temperature and the zone air dry-bulb temperature (W)

˙Qrad = radiant heat transfer to the refrigerated case (W)

˙Qinf,sens = sensible heat transfer by air infiltration to the refrigerated case through the air curtain or via door openings (W)

˙Qinf,lat = latent heat transfer by air infiltration to the refrigerated case through the air curtain or via door openings (W)

˙Qlights = lighting heat load (W)

˙Qas = anti-sweat heater load (W)

˙Qdef = defrost heat load (W)

˙Qfan = fan heat load (W)

˙Qrestock = sensible load on the refrigerated case due to restocking of products that are at a higher temperature than the case (W)

The model assumes that these load components are known for a refrigerated case at rated ambient air conditions (typically 23.9˚C [75˚F] and 55% relative humidity) and the specified case operating temperature. A combination of user input curves and fixed correlations (defined within EnergyPlus) adjust for case performance at off-rated conditions. Several of the load components are typically provided by the case manufacturer (e.g., total rated load, fan, lighting, anti-sweat heater, and defrost loads). The remaining load components are not usually provided by the manufacturer and must be estimated (heat conduction through case walls, radiation heat transfer, sensible/latent air infiltration, and restocking).

For estimating the latent air infiltration load, the model requires that the user provide the latent heat ratio (LHR) for the refrigerated case at rated conditions. Research results are available to provide guidance in selecting this value (ASHRAE 2002, Howell 1993a, Howell 1993b). The rated LHR for refrigerated cases typically ranges from 0.1 to 0.3 depending on case configuration (e.g., glass door reach-in versus multi-deck open case) and case operating temperature.

The case loads due to wall heat conduction, radiation, and sensible air infiltration are estimated by the model as a single lumped value (sensible case credits). The sensible case credits are calculated by subtracting the known loads at rated conditions (fan, lighting, anti-sweat heater, defrost and latent case credits) from the rated total cooling capacity of the case which is provided by the case manufacturer ( ˙Qcase,rated ).

Using these assumptions and the schedule inputs provided by the user, the refrigerated case evaporator load components in Equation are determined for each simulation time step. The variation in certain loads with respect to changes in ambient air temperature and/or humidity (e.g., latent and sensible case credits, defrost load, and anti-sweat heater load) are factored into the calculation based on user-provided inputs or by the model itself.

Whenever the total heat load on the case is greater than the available evaporator capacity, such as during defrost (when the evaporator capacity is set to zero) or restocking, the load is accumulated to be met during subsequent time steps. This accounts for the energy required to bring the case back down to the rated operating temperature even though the rise in case temperature during defrost or restocking is not explicitly modeled. Following defrost, it may take multiple time steps to meet this accumulated load.

The specific calculations for case evaporator load components and electric power for these loads (as applicable) are provided below.

Case Evaporator Fan[LINK]

The refrigerated case evaporator fan electric power is calculated for each simulation time step as the product of the operating case fan power per unit length of case, the length of the refrigerated case, and the fraction of time that the case is not being defrosted. For cases with hot-gas or electric defrost (with or without temperature termination), the fan is disabled during the entire scheduled defrost drip-down time period. The evaporator fan operates continuously for off-cycle defrost or no defrost.

Pfan=Pfan,oper(Lcase)(1SCHdefrost,dripdown)

where:

Pfan = output variable “Refrigerated Case Evaporator Fan Electric Power [W]”

Pfan,oper = operating case fan power per unit length (W/m)

Lcase = case length (m)

SCHdefrost,dripdown = fraction of time case is being defrosted (0 to 1), including drip-down period (based on the defrost drip-down schedule) for hot-gas or electric defrost. For off-cycle defrost or no defrost, this value is set to zero for this calculation.

The model assumes that the evaporator fan is entirely contained within the thermal envelope of the case, and that all fan power results in a direct heat load on the case evaporator:

˙Qfan=Pfan

Case Lighting[LINK]

The refrigerated case lighting electric power is calculated for each simulation time step as the product of the installed case lighting power per unit length of case, the lighting schedule value, and the length of the refrigerated case:

Plights=Plights,installed(Lcase)(SCHlights)

where:

Plights = output variable “Refrigerated Case Lighting Electric Power [W]”

Plights,installed = installed case lighting power per unit length (W/m)

SCHlights = case lighting schedule value (0 to 1)

A maximum schedule value of 1.0 means the lights are fully on at the installed case lighting power level. Schedule values of 0.0 indicate the lights are off and 0.5 at half-power.

The user can specify the fraction of lighting energy that directly contributes to the case evaporator heat load:

˙Qlights=Plights(Fl)

where:

Fl = fraction of lighting energy to case

The remainder of the lighting energy (1 - Fl) is a heating load to the zone where the case is located, which is discussed further in section Sensible Case Credits below. This fraction (1 - Fl) can be used to represent lighting ballasts and/or bulbs located outside the air curtain of the refrigerated case.

Anti-Sweat Heater Performance[LINK]

Anti-sweat heaters warm the refrigerated case rails or doors to provide protection from moisture condensation. Different anti-sweat heater control strategies are used depending on the case temperature and the type of anti-sweat heater installed. Several types of anti-sweat heater control strategies can be simulated with this model: constant, linear variation with ambient relative humidity or dewpoint temperature, and a theoretical model that determines the minimum anti-sweat heater power required to maintain the case surface just above the temperature where condensation would occur. Additionally, anti-sweat heater performance can be disregarded if the type of refrigerated case does not warrant its use. For the control strategies described below (except “None” and “Constant Method”), the model does not allow the anti-sweat heater power to be less than the minimum power nor greater than the case anti-sweat heater power specified by the user. Each anti-sweat heater control type is described in detail below.

None[LINK]

Used for refrigerated cases that do not require an anti-sweat heater.

˙Qas=0

where:

˙Qas = anti-sweat heater load on the case evaporator (W)

Constant Method[LINK]

For refrigerated cases requiring constant anti-sweat heater output, the power use is simply calculated as the case anti-sweat heater power per unit length multiplied by the length of the case. This method is used when the manufacturer recommends that cycling of the heaters not occur.

Pas=Pas(Lcase)

where:

Pas = output variable “Refrigerated Case Anti-Sweat Heater Electric Power [W]”

Pas = case anti-sweat heater power per unit length (W)

Relative Humidity Method[LINK]

Anti-sweat heater power can be reduced at lower ambient relative humidity levels to save energy while still protecting from moisture condensation on cold surfaces. For this control type, anti-sweat heater power use is reduced linearly based on case anti-sweat heater power at the rated ambient relative humidity (typically 55% RH), the relative humidity specified by the user where no anti-sweat heater power is required, and the relative humidity of the ambient (zone) air surrounding the case.

Pas=Pas(Lcase)(1[RHratedRHairRHratedRHmin])

where:

RHair = relative humidity of the ambient (zone) air (%)

RHrated = rated ambient relative humidity (%)

RHmin = relative humidity at zero anti-sweat heater energy (%)

Dewpoint Method[LINK]

Anti-sweat heater power can also be reduced as a function of ambient air dewpoint temperature based on a similar correlation to that used by the relative humidity method. This control method varies the anti-sweat heater power linearly based on the ambient air dewpoint temperature, the case operating temperature, and the rated ambient dewpoint temperature (calculated by the model using the rated ambient temperature and rated ambient relative humidity entered by the user).

Pas=Pas(Lcase)(Tdp,airTcaseTdp,ratedTcase)

where:

Tdp,air = dewpoint temperature of the ambient (zone) air (˚C)

Tdp,rated = rated ambient dewpoint temperature (˚C)

Tcase = case operating temperature (˚C)

Heat Balance Method[LINK]

A theoretical model may also be used to simulate the performance of anti-sweat heater operation at various indoor dewpoint temperatures (Henderson and Khattar 1999). The model calculates that amount of heat required to hold the case or door surface at (or slightly above) the dewpoint temperature of the ambient air using the following simple heat balance equation:

Pas=((Tdp,airTdb,air)HcaseRair+(Tdp,airTcase)HcaseRcase)Lcase

where:

Tdp,air = dewpoint temperature of the ambient (zone) air (˚C)

Tdb,air = dry-bulb temperature of the ambient (zone) air (˚C)

Hcase = height of the case (m)

Rair = air film resistance (assumed constant at 0.3169 m2-˚C/W)

Rcase = heat transfer resistance of case (m2-˚C/W)

Tcase = case operating temperature (˚C)

Lcase = case length (m)

The model above provides a linear relationship of anti-sweat heater power with varying ambient air dewpoint temperature at constant ambient air dry-bulb and case temperatures. By assuming that the ‘nominal’ anti-sweat heater power entered by the user is required to avoid moisture condensation at rated ambient air conditions, the value of Rcase can be determined by rearranging the equation and solving as follows:

Rcase=(Tdp,ratedTcase)(PasHcase)(Tdp,ratedTdb,ratedRair)

where:

Tdb,rated = rated ambient temperature (˚C)

With Rcase known, Pas can be calculated for each simulation time step using the actual ambient (zone) air dry-bulb and dewpoint temperatures.

All Anti-Sweat Heater Control Methods[LINK]

For all control methods, the user can specify the fraction of anti-sweat heater energy that directly contributes to the case evaporator heat load:

˙Qas=Pas(Fas)

where:

Fas = fraction of anti-sweat heater energy to case

The remainder of the anti-sweat heater energy (1 - Fas) is a heating load to the zone where the case is located, which is discussed further in section Sensible Case Credits below.

Case Restocking[LINK]

The impact of restocking the refrigerated case with product that is not at the case operating temperature is modeled with the case restocking schedule. The schedule is entered as a heat gain rate per unit length of the refrigerated case (W/m). The heat load due to restocking is calculated as the scheduled load multiplied by the length of the refrigerated case. The load due to product restocking is assumed to be only sensible (temperature) heat; a latent (moisture) component is not modeled.

˙Qrestock=SCHrestock(Lcase)

where:

SCHrestock = refrigerated case restocking schedule value (W/m)

The restocking heat load is removed by the refrigerated case evaporator any time the case is not being defrosted and excess sensible cooling capacity is available. If the evaporator cooling capacity is insufficient to remove the entire restocking load, the unmet portion is carried over to the next simulation time step.

Case Defrost[LINK]

Eight refrigerated case defrost strategies can be simulated: none, off-cycle, electric, electric with temperature termination, hot-gas, hot-gas with temperature termination, hot-brine, and hot-brine with temperature termination. Some research has shown that the defrost times for cases defrosted using hot brine can be significantly shorter than defrost times for electric or hot gas.(Terrell, W. J. Jr., 1999) For each of these strategies, the refrigerated case evaporator is turned off for the required time period to allow accumulated frost to melt. Additional time can be scheduled (drip-down) to allow the water to drip from the evaporator and drain from the case.

Refrigerated cases typically require a specific number of defrost cycles per day for a pre-determined length of time. Refer to manufacturer’s recommendations for proper defrost frequency and duration. For example, a refrigerated case may have a single defrost period each day with defrost scheduled from 7:00 - 7:40 am and defrost drip-down scheduled from 7:00 - 7:55 am. Notice the drip-down schedule and the defrost schedule start at the same time, and the drip-down schedule is longer than the defrost schedule. These schedules should normally repeat for each day of the year.

For electric, hot gas, and hot brine defrost types, energy use by the defrost heater occurs during the scheduled defrost period. For defrost with temperature termination, the energy is also multiplied by the defrost ratio simulating a defrost duration shorter than the defined (maximum) period. For all non-electric defrost types, defrost electric power is set equal to zero (and is not available as an output variable). For hot gas and hot brine defrost types in cases served by a detailed system, the condenser heat rejection load is reduced by the amount of heat recovered for use in the defrost system. This condenser credit is not applied for the simple compressor rack system.

If(DefrostType=Electric)ThenPdef=Pdef(Lcase)(SCHdefrost)ElseIf(DefrostType=ElectricWithTempTermination)ThenPdef=Pdef(Lcase)(SCHdefrost)(DefrostRatio)ElsePdef=0.0EndIf

where:

Pdef = output variable “Refrigerated Case Defrost Electric Power [W]”

Pdef = case defrost power per unit length (W)

Lcase = case length (m)

SCHdefrost = case defrost schedule value (0 to 1)

DefrostRatio = fraction of maximum defrost time, used with temperature termination

Frost accumulation on the case evaporator will vary with the humidity level in the ambient air surrounding the case. Therefore, defrost heater operation can be reduced when ambient air humidity levels are low. Several methods are used to reduce unnecessary defrost heater operation, including terminating heater operation when the measured evaporator temperature indicates that the accumulated frost has been completely melted. For modeling refrigerated cases with temperature-terminated defrost, EnergyPlus allows the user to specify a defrost energy correction curve to account for variations in defrost energy as ambient air humidity levels change. The user can select from four correction curve types: None, Case Temperature Method, Relative Humidity Method, or Dewpoint Method.

None(default):DefrostRatio=1CaseTemperatureMethod:DefrostRatio=1(RHratedRHair)[a+b(Tcase)+c(Tcase)2+d(Tcase)3]RHmethod:DefrostRatio=e+f(RHair)+g(RHair)2+h(RHair)3Dewpointmethod:DefrostRatio=i+j(Tdp,air)+k(Tdp,air)2+l(Tdp,air)3

where:

RHrated = rated ambient relative humidity (%)

RHair = relative humidity of the ambient (zone) air (%)

Tcase = case operating temperature (˚C)

Tdp,air = dewpoint temperature of the ambient (zone) air (˚C)

a…l = user-defined coefficients using a cubic curve object (Curve:Cubic)

The user specifies the defrost energy correction curve type and the name of the cubic curve object (Curve:Cubic) that defines the curve coefficients. Representative curve coefficients for curve type “Case Temperature Method” are provided in Table 75.

Representative Defrost Energy Correction Curve Coefficients for Case Temperature Method
Coefficient Single-shelf horizontal display case Multi-shelf vertical display case
a 2.3632E-2 2.4598E-2
b 6.2320E-4 7.6439E-4
c 2.8320E-5 -3.8637E-5
d 4.4035E-7 7.45686E-7

: Representative Defrost Energy Correction Curve Coefficients for Case Temperature Method

Note: Coefficients derived for RHrated = 55% and a rated ambient temperature of 23.9˚C (75˚F). Source: Howell 1993b.

As mentioned above, the refrigerated case evaporator is turned off while it is being defrosted. Heat gains during defrost must be removed once the defrost period (drip-down schedule) has ended. The model assumes that heat gains due to defrost heater operation are at least partially offset by converting accumulated frost to liquid water (condensate) which drains from the case. Frost accumulation during each simulation time step is estimated by the model using the actual latent heat transfer to the refrigerated case and the heat of vaporization plus the heat of fusion for water. The model assumes that frost is not accumulated on the evaporator during the defrost drip-down time period.

Frost=Frost+˙Qcase,rated(Lcase)(RTFrated)(LHRrated)(LatentRatio)(tzn)(hf+hfg)(1SCHdefrost,dripdown)

where:

Frost = amount of accumulated frost on the case evaporator (kg)

˙Qcase,rated = case rated total cooling capacity per unit length (W/m)

Lcase = case length (m)

RTFrated = runtime fraction of the refrigerated case at rated conditions

LHRrated = latent heat ratio of the refrigerated case at rated conditions

LatentRatio = ratio of actual latent load to rated latent load on the case, based on latent case credit curve (see section Latent Case Credits below)

tzn = duration of zone simulation time step (s)

hfg = heat of vaporization of water (assumed constant at 2,498,000 J/kg)

hf = heat of fusion of water (335,000 J/kg)

SCHdefrost,dripdown = defrost drip-down schedule value (0 to 1)

During defrost (SCHdefrost), the model assumes that the hot gas, hot brine, or electric heater energy directly contributes to melting the frost (heat of fusion of water). Defrost energy not attributed to melting frost from the evaporator coil results in a heat load on the refrigerated case evaporator (˙Qdef). When the defrost drip-down time period ends, this defrost energy heat load is added to the actual case load (up to the maximum evaporator capacity) until the total defrost energy heat load is removed (which may take several simulation time steps)

If(DefrostType=ElectricorHotGasorHotBrine)Then˙Qdef=MAX(0.0,[Pdef(Lcase)(SCHdef)Frost(hf)tzn])Else˙Qdef=0.0Endif

where:

˙Qdef = defrost heat load (W)

Sensible Case Credits[LINK]

Refrigerated cases remove sensible energy from the surrounding environment (termed “sensible case credits”). In this model, the sensible case credits are composed of wall heat conduction, radiation heat transfer, and sensible heat transfer by air infiltration ( ˙Qwalls + ˙Qrad + ˙Qinf,sens in equation ). To quantify this energy transfer, the model first calculates the rated sensible case credits by subtracting the known loads at rated conditions (fan, lighting, and anti-sweat heater) from the rated sensible cooling capacity of the case. It should be noted that the lighting and fan heat discussed here are for standard-efficiency equipment. Manufacturers typically provide ratings for both standard and high-efficiency fan and lighting equipment; however, the standard equipment is used to determine rated sensible case credits. (Some manufacturers no longer include any lighting in their rated capacity values. For these cases, P’lights,std will equal zero.)

˙Qccsens,rated=[˙Qcase,rated(RTFrated)(1LHRrated)Plights,std(Fl)Pas(Fas)Pfan,std]Lcase

where:

˙Qccsens,rated = sensible case credits at rated conditions (W)

˙Qcase,rated = case rated total cooling capacity per unit length (W/m)

RTFrated = runtime fraction of the refrigerated case at rated conditions

LHRrated = latent heat ratio of the refrigerated case at rated conditions

Plights,std = standard case lighting power per unit length (W/m)

Fl = fraction of lighting energy to case

Pas = case anti-sweat heater power per unit length (W)

Fas = fraction of anti-sweat heater energy to case

Pfan,std = standard case fan power per unit length (W/m)

Lcase = case length (m)

For every simulation time step, the rated sensible case credits are then adjusted to account for variations at off-rated ambient air temperatures. The model also allows the user to define a case credit fraction using a schedule object. This case credit fraction can be useful for modeling cases that operate differently during specific time periods. For example, metal or plastic coverings may be installed on refrigerated display cases during unoccupied hours which would significantly reduce case credits (e.g., air infiltration) compared to occupied hours when the coverings are removed. If the user does not define a case credit fraction schedule, then the fraction is assumed to be 1 for the entire simulation.

˙Qccsens=˙Qccsens,rated(Tdb,airTcaseTdb,ratedTcase)(SCHcc)

where:

˙Qccsens = sensible case credits adjusted for ambient temperature and case credit fraction (W)

Tdb,air = dry-bulb temperature of the ambient (zone) air (˚C)

Tcase = case operating temperature (˚C)

Tdb,rated = rated ambient (zone) dry-bulb temperature (˚C)

SCHcc = case credit fraction (schedule value, 0 to 1)

The sensible case credits calculated above are considered heat loads on the refrigerated case evaporator. The net impact of the case credits on the surrounding zone includes adjustment for the portion of the lighting and anti-sweat heater power that does not directly contribute to the case evaporator load. Sensible case credits are negative values when heat is removed from the zone load.

˙Qccsens,NET=Plights(1Fl)+Pas(1Fas)˙Qccsens

where:

˙Qccsens,NET = net impact of the sensible case credits on the surrounding zone, negative for cooling (W)

Plights = case lighting electric power (W)

Fl = fraction of lighting energy to case

Pas = anti-sweat heater electric power (W)

Fas = fraction of anti-sweat heater energy to case

When refrigerated cases are served by a compressor rack that rejects condenser waste heat directly to the zone (e.g., a stand-alone refrigerated case with integral compressor and condenser), this condenser waste heat offsets some or all of the sensible case credits. The amount of condenser waste heat rejected to the zone is calculated and reported by the refrigerated case compressor rack object (Ref. Heat Rejection to Zone).

Latent Case Credits[LINK]

Refrigerated cases also remove latent energy (moisture) from the surrounding environment (termed “latent case credits”). In this model, the latent case credit is composed solely of the latent heat transfer by air infiltration ˙Qinf,lat in equation . The latent case credits are calculated as the product of the case length and the total cooling capacity per unit length, latent heat ratio, and runtime fraction at rated conditions. As described previously (Ref. Sensible Case Credits), a case credit fraction schedule is used to model cases that operate differently during specific time periods. The same case credit fraction is used to modify both the sensible and latent case credits. If the user does not define a case credit fraction schedule, then the fraction is assumed to be 1 for the entire simulation. The calculation of latent case credits also includes a factor (LatentRatio) that accounts for lower ambient humidity levels. Latent case credits are set to zero during the defrost-dripdown periods.

˙Qinf,lat=˙Qcclat=˙Qcase,rated(LHRrated)(RTFrated)(SCHcc)(LatentRatio)Lcase

where:

˙Qinf,lat = latent load on the refrigerated case evaporator at current ambient conditions (W)

˙Qcclat = latent case credit impact on zone load, negative for dehumidification (W)

˙Qcase,rated = case rated total cooling capacity per unit length (W/m)

LHRrated = latent heat ratio of the refrigerated case at rated conditions

RTFrated = runtime fraction of the refrigerated case at rated conditions

SCHCC = case credit fraction (schedule value, 0 to 1)

LatentRatio = ratio of actual latent load to rated latent load on the case, based on latent case credit curve

Lcase = case length (m)

Latent load on the refrigerated case evaporator will vary with ambient humidity levels. Therefore, the refrigerated case model allows the user to specify a latent case credit curve to adjust case credits based on ambient humidity, and the user can select from three curve types: Case Temperature Method, Relative Humidity Method, or Dewpoint Method.

CaseTemperatureMethod:LatentRatio=1(RHratedRHair)[m+n(Tcase)+o(Tcase)2+p(Tcase)3]RHmethod:LatentRatio=q+r(RHair)+s(RHair)2+t(RHair)3Dewpointmethod:LatentRatio=u+v(Tdp,air)+w(Tdp,air)2+x(Tdp,air)3

where:

RHrated = rated ambient relative humidity (%)

RHair = relative humidity of the ambient (zone) air (%)

Tcase = case operating temperature (˚C)

Tdp,air = dewpoint temperature of the ambient (zone) air (˚C)

m…x= user-defined coefficients using a cubic curve object (Curve:Cubic)

The user specifies the latent case credit curve type and the name of the cubic curve object (Curve:Cubic) that defines the curve coefficients. Representative curve coefficients for curve type “Case Temperature Method” are provided in Table 76.

Representative Latent Case Credit Curve Coefficients for Case Temperature Method
Coefficient Single-shelf horizontal Multi-shelf vertical
m 2.0376E-2 2.6520E-2
n 2.4378E-4 1.0780E-3
o 1.1400E-5 -6.0256E-5
p 1.8110E-7 1.2373E-6

Note: Coefficients derived for RHrated = 55% and a rated ambient temperature of 23.9˚C (75˚F). Source: Howell 1993b.

Refrigerated Case Credits With Under Case Return Air[LINK]

For certain refrigerated case types, the sensible case credits provided to the zone can create an uncomfortably cold environment in the surrounding area. For this reason, return air ducts are frequently placed behind these cases to draw this cold air under the case and direct it back to the HVAC system. This reduces localized over-cooling and improves occupant comfort.

Return Air Factor Versus Under Case HVAC Return Air Fraction

Return Air Factor Versus Under Case HVAC Return Air Fraction

Since under case return ducts reduce the temperature and humidity of the air being recirculated to the HVAC system, this can impact HVAC system performance. Figure 279 shows the relationship that is used by the refrigerated case model to determine the fraction of case credits that directly cool and dehumidify the HVAC system return air. This fraction, referred to as the Return Air Factor (RAF), is a function of the fraction of the HVAC system return air that comes from under the cases. The remaining fraction of the case credits (1-RAF) becomes part of the overall zone air energy balance. If the HVAC system is off for a simulation time step (no return air mass flow), the sensible and latent case credits normally attributed to the HVAC return are set equal to zero (even though they get calculated and reported here as non-zero values) and all case credit energy is applied to the zone air heat balance.

˙Qccsens,zone=˙Qccsens,NET(1RAF)

˙Qcclat,zone=˙Qcclat(1RAF)

˙Qccsens,HVAC=˙Qccsens,NET(RAF)

˙Qcclat,HVAC=˙Qcclat(RAF)

where:

˙Qccsens,zone = sensible case credit applied to the zone air heat balance (W)

˙Qcclat,zone = latent case credit applied to the zone air heat balance (W)

˙Qccsens,HVAC = sensible case credit applied to the HVAC return air (zone return air path outlet node) (W)

˙Qcclat,HVAC = latent case credit applied to the HVAC return air (zone return air path outlet node) (W)

RAF= return air factor (see Figure 279 above)

Variable Evaporator Temperature[LINK]

Control systems are now available that increase the evaporator temperature to improve compressor efficiency whenever the total loads on a system are less than the system capacity. To model these systems, a variable evaporator temperature is an option available with the detailed refrigeration system object (Refrigeration:System). If this option is selected, the model will compare the refrigeration load on each case to the load at rated conditions. If the case load in a particular time step is less than the rated load, an acceptable elevated evaporator temperature is determined for that case. The evaporator temperature for the whole refrigeration system is then set by the minimum evaporator temperature needed for any particular case.

LFcase=˙Qcase,actual˙Qcase,rated;0.5LFcase1.0TEvap,Allowed=TcaseLFcase(TcaseTEvap,Design)

where:

LFcase = Load factor for a particular case

Tevap= Evaporator temperature, C.

Walk-In Coolers and Freezers[LINK]

The walk-in object (Refrigeration:WalkIn) is another type of refrigeration load that can be placed on either a refrigeration compressor rack, detailed refrigeration system, or secondary refrigeration system object (Refrigeration:CompressorRack, Refrigeration:System, or Refrigeration:SecondarySystem). Walk-in coolers and freezers differ from refrigerated cases in that they may have surfaces facing more than one zone and in that they are always equipped with doors, that is, they do not have open shelves. Their sensible and latent exchange with zones is therefore calculated in a different manner than for refrigerated cases. Also, the walk-in model does not interact directly with the HVAC system, that is, the return air fraction option available in the refrigerated case model is not included.

The walk-in cooler performance is based on the ASHRAE load model, which includes infiltration through door openings and sensible loss through walls/ceilings described by the user for each zone.(ASHRAE 2006d, ASHRAE 2006e, Gosney, W.B., Olama, G.A.-L. 1975) All equipment loads (fan, light, heaters) are modeled as well. Sensible and latent exchange with multiple adjoining zones is included. A master schedule is used for the Walk In operation and additional schedules control the lights, defrost, and heater operation. Just as for cases, unmet refrigeration loads are accumulated to be met the following time step. This usually occurs during defrost and restocking.

Walk-In Sensible and Latent Heat Exchange[LINK]

A walk-in can exchange both sensible and latent energy with multiple zones. The heat transfer calculations are performed separately for each zone so that the heat transfer impact, or zone credits, can be determined. The area of all walls and ceilings facing each zone are described by the user by their thermal conductance and area. Sensible energy exchange takes place between these surfaces and the surrounding zones. Because these walls interface with conditioned zones at relatively constant temperatures, this heat exchange is modeled very simply:

QSurfacesZn = USurfacesZn x ASurfacesZn x ΔTZn

QDoorSensZn = UDoorZn x AreaDoorZn x ΔTZn

The heat transfer through the floor is similarly modeled.

QFloor = AFloor x UFloor x (TGround - TWalkIn)

Where:

AFloor = Area of the walkin floor, m2

ASurfacesZn= Area of surfaces facing Zone n, m2

QDoorSensZn = Sensible heat transfer through the closed door(s) facing Zone n, W

QsurfacesZn = Sensible heat transfer through walls and ceilings facing Zone n, W

TGround = Ground temperature, C

TWalkIn= Walk-in operating temperature, C

UFloor = Thermal conductance of floor, W/m2-K

UDoorZn = Thermal conductance of doors facing Zone n, W/m2-K

USurfacesZn= Thermal conductance of surfaces facing Zone n, W/m2-K

ΔTZn = Difference between walk-in operating temperature and Zone n drybulb temperature, C

Infiltration through doorways places both a sensible and a latent load upon the walk-in, and corresponding credits upon the adjacent zone. Two types of doors are available, nominally called ‘stock’ and ‘glass’ doors, to enable the user to model doors that differ in thermal conductance, door protection type, and frequency of opening. The sensible and latent infiltration loads are modeled according to the guidance specified in (ASHRAE 2006d, ASHRAE 2009, and Gosney and Olama, 1975). The air within the cooler is assumed to be at 90% relative humidity. Equal air exchange is assumed, that is, the mass of dry air infiltrating into the walkin is assumed to equal the mass of dry air infiltrating out of the walkin.

QInfiltration = QFullFlow x FactorDoorOpen x FactorFlow x (1 - FactorProtection)

QFullFlow=0.221*ADoor(hZoneAir-hAirWalkInAirWalkIn(1-ρZoneAirAirWalkIn)0.5(g*HDoor)0.5FactorDensity

FactorDensity = (2 /(1 + (ρAirWalkIn / ρZoneAir)0.333)) 1.5

mDryAir = QInfiltration / (hZoneAir - hAirWalkIn)

mWater = mDryAir x (WZoneAir - WAirWalkIn)

QWalkInLatentZn = mWater x ΔhIcetoVapor x (1 - SCHDefrost,DripDown)

QWalkInSensInfZn = QInfiltration - (mWater x ΔhIcetoVapor)

Where:

Adoor= Area of door facing Zone n, m2

FactorDoorOpen = Value scheduled by user, fraction of time door open during time step

FactorFlow = Doorway flow factor, = 0.8 if ΔTZn > 11C; =1.1 if ΔTZn <= 11C

FactorProtection = Doorway protection factor, = 0 for no protection; = 0.5 for an air curtain; and 0.9 for a strip curtain

g= Gravitational constant

hAirWalkIn= enthalpy of the air within the walk in, = f(TWalkIn,POudoor, 90%RH), J/kg

hZoneAir= enthalpy of the air in Zone n, J/kg

Hdoor= Height of door facing Zone n, m

QFullFlow = Sensible and latent refrigeration load for fully established flow, W

QInfiltration = Average infiltration (sensible and latent) refrigeration load for the time step, W

QWalkInLatentZn = Latent load upon the walk in facing Zone n, W

QWalkInSensInfZn = Sensible load due to infiltration upon the walkin facing Zone n, W

mDryAir = Mass of dry air infiltrating into the walk-in, kg/s

mWater = Mass of water removed from the infiltrating air, kg/s

POudoor= Outdoor air pressure, Pa

SCHDefrost,DripDown = value from 0 to 1 indicating whether the system is in the dripdown period

WAirWalkIn= Humidity ratio of the air within the walk in, = f(TWalkIn,POudoor, 90%RH), kg/kg

WZoneAir = Humidity ratio of Zone n air, kg/kg

ΔhIcetoVapor = Latent heat absorbed to change ice to vapor, J/kg

ρAirWalkIn= Density of the air within the walk in = f(TWalkIn,POudoor, 90%RH), kg/m3

ρZoneAir = Density of air in Zone n, kg/m3

The sensible load on the case and the sensible credit to the zone continue throughout the defrost and dripdown periods. However, to be consistent with the treatment of refrigerated cases, there is no latent credit to the zone or latent load upon the cooler during the dripdown period. Latent load and latent credit are both based on reducing the infiltrating vapor to ice. The sensible heat exchange between the walk in and the zone is then the total of the heat transfer through the doors and surfaces and the infiltration sensible load. The latent load upon the walkin is converted to the amount of frost added to the coils during each time step. This accumulating value is used later to determine the load placed upon the walkin during the defrost cycle.

QWalkInSensZn = QWalkInSensInfZn+QDoorZn + QsurfacesZn

QZoneLatent = - QWalkInLatentZn

QZoneSens = - QWalkInSensZn

ΔFrostZn = (mWater *Δtime)* (1- SCHDefrost,DripDown)

Where:

QWalkInSensZn = Total sensible heat exchange between the walkin and Zone n, W

QZoneLatent = Latent load upon the Zone n, W

QZoneSens = Sensible load upon Zone n , W

ΔFrostZn= Change in frost inventory, kg

Δtime= Length of time step, s

After the heat exchange with each zone is calculated, the total load on the walkin is calculated:

QWalkInLatentTot = ∑QWalkInLatentZn

QWalkInSensTot= ∑QWalkInSensZn~ + QLight + QFan + QHeater + QDefrost + QStocking + QFloor

QWalkInTotal = QWalkInLatentTot+ QWalkInSensTot

ΔFrostTot= ∑ΔFrostZn

Where QLight, QFan, QHeater, QStocking , and QDefrost are described below.

Walk-In Fans, Heaters, Lighting, and Restocking[LINK]

Sensible heat loads are placed on a walk-in by fans, heaters, and lighting. Unlike refrigerated cases, there is no option to allocate any portion of these heat loads to the surrounding zone(s). Larger walk-ins will have separate fans at the cooling coil and for general circulation. The general circulation fan is assumed to run at all times. The cooling coil fan is assumed to be off for Hot-Fluid and Electric defrost. Lighting, heating, and restocking are modeled according to the schedule values entered by the user. For lighting and heating, the maximum power is entered along with a scheduled ratio (between 0 and 1) to be applied for any point in time. The heating power includes all heaters except those used for defrost purposes. The heater power should include anti-sweat, door, floor, and drain-pan heaters. For restocking, the total sensible load in Watts is scheduled for each point in time (the restocking latent load is assumed to be zero).

QLight = RatedQLighting * SCHLighting

QFan = PowerCircFan + PowerCoilFan * ( 1 - SCHDripDown )

QHeater = PowerHeater * SCHHeater

QStocking = SCHStocking

Where:

QLight = Refrigeration load due to lighting during current time step, W

RatedQLighting = Maximum lighting load specified for the walk-in, W

SCHLighting = Scheduled value between 0 and 1 for the current time step

QFan = Refrigeration load due to fan power during the current time step, W

PowerCircFan = Rated circulating fan power, W

PowerCoilFan = Rated coil fan power, W

SCHDripDown = Scheduled value between 0 and 1 for the current time step

QHeater = Refrigeration load due to heaters during current time step, W

PowerHeater= Rated total heater(s) power (including anti-sweat, floor, door, etc.) , W

SCHHeater = Scheduled value between 0 and 1 for the current time step

QStocking = Refrigeration load due to stocking during the time step, W

SCHStocking = Scheduled value of load due to stocking, W

Defrost[LINK]

The defrost types available for the walk-in model include none, off-cycle, electric, and hot-fluid. Defrosts are started according to scheduled times and can be ended either by schedule or by temperature termination. Dripdown schedules are used to keep the cooling coil off long enough to drain any condensate from the system.

For defrost types none and off-cycle, the refrigeration load on the walk-in due to defrost is zero. For off-cycle, the walk-in refrigeration capacity is set to zero during the drip-down scheduled time.

The energy required for hot-fluid defrost is assumed to be reclaimed from the compressor exhaust (for detailed systems, this energy appears as a credit against the heat rejection needed at the condenser). The energy used by electric defrost is available as an output variable.

If the defrost cycle is controlled by the schedule, the refrigeration load placed upon the walk-in is calculated as the product of the defrost capacity and the defrost schedule. The load is then reduced according to the amount of accumulated ice melted during that time step.

QDefrost = CapacityDefrost*SCHDefrost - Δfrost x ΔhIceMelt / Δtime

Where:

QDefrost = Refrigeration load imposed by defrost heat, W

CapacityDefrost = Rated defrost power, W

SCHDefrost = Scheduled value between 0 and 1 for the current time step

Δfrost = amount of frost melted during time step, kg

ΔhIceMelt = heat of fusion for ice, J/kg

Δtime = time in time step, s

If the defrost is controlled by temperature termination, the defrost cycle is assumed to end when all the ice is melted. However, we need to recognize not all defrost heat goes to melt ice. Some of the defrost heat goes to raising the temperature of the coil mass to greater than 0C, and some is transferred to the walk-in environment as some of the coils are defrosted before others. The user enters a ‘defrost energy fraction’ to specify the portion of the defrost energy that goes directly to melting ice. The default for defrost energy fraction is 0.7 for electric defrost and 0.3 for warm fluid defrost.( Baxter, V. D., Mei, V.C., 2002) For this type of defrost control, the model calculates the amount of energy available to melt the ice in each time step. The accumulated amount of ice is then reduced accordingly. When all the ice is melted, the defrost schedule value is set to zero and no further defrost load is placed upon the walk-in cooler. If the defrost schedule ends before the ice is melted, the schedule is used and the ice continues to accumulate until the next defrost cycle. The refrigeration capacity is kept at zero until the end of the drip-down schedule. Until the accumulated ice is melted, the defrost heat load upon the walk-in is:

QDefrost = CapacityDefrost x SCHDefrost x (1- FractionDefrostEnergy)

Air Chillers and Air Chiller Sets[LINK]

The Air Chiller object (Refrigeration:AirChiller) is another type of refrigeration load that can be placed on either a refrigeration compressor rack, detailed refrigeration system, or secondary refrigeration system object (Refrigeration:CompressorRack, Refrigeration:System, or Refrigeration:SecondarySystem). Air chillers are used to model the type of equipment typically used in refrigerated warehouses. For that reason, there is a major difference between the air chiller model and those for refrigerated cases or walk-ins. For cases and walk-ins, a portion of the model is directed toward calculating the amount of refrigeration needed to maintain the refrigerated volume at the desired temperature due to heat exchange with the surrounding zone, and that zone is conditioned to a nearly constant temperature. In a refrigerated warehouse, the refrigeration load is caused by heat exchange with a variable external environment. For that reason, the loads for these zones are calculated by the usual EnergyPlus zone heat balance. The amount of refrigeration needed to maintain the specified temperature set points is then passed to the air chiller model, in a similar fashion to the load passed to a window air conditioner model. The air chillers are therefore solved using the system time step, not the zone time step used for cases and walk-ins.

The air chiller performance is based on three types of manufacturers ratings, Unit Load Factor, Total Capacity Map, or a set of European standards. Correction factors for material and refrigerant are applied to all of these ratings.

Unit Load Factor Capacity[LINK]

Bruce Nelson has provided a useful description of the Unit Load Factor approach.(Nelson, B.I., 2010)

“One well-known method used to calculate the sensible cooling capacity of evaporators is the effectiveness method.(Kays,* *W.M., A.L. London, 1964) Heat exchanger effectiveness is defined as the ratio of the actual amount of heat transferred to the maximum possible amount of heat that could be transferred with an infinite area. This method is extremely useful because cooling capacity can be calculated directly knowing only the dimensional characteristics of the coil and the initial temperature difference (entering air temperature minus the evaporating temperature). This initial temperature difference is referred to as”DT1" … in the refrigeration industry. Sensible cooling capacity is calculated as follows:

qsens=˙m×cp×ε×(TcoilinletTevap)=˙m×cp×ε×DT1

For a given size of coil operating with constant airflow rate, the effectiveness can be considered constant over the small operating temperature ranges typical of refrigeration applications, and therefore, capacity can be considered to be proportional to the ratio of DT1. Hence, if evaporator coil sensible capacity is known for a given DT1, then capacity at a new initial temperature difference, DT1␣, can be found by multiplying the original capacity by the ratio DT1␣/DT1."

Where:

qsens = Cooling capacity (sensible only), W

˙m = mass flow rate of air, kg/s

cp = specific heat capacity of moist air, J/kg-C

ε = effectiveness (=(Tcoil inlet - Tcoil exit)/(Tcoil inlet - Tevap)

Tcoil inlet = dry-bulb air temperature entering the coil, C

Tevap = average refrigeratnt evaporating temperature, C

Tcoil exit= dry-bulb air temperature leaving the coil, C

DT1 = initial temperature difference, C

Using this approach, the manufacturer specifies the Unit Load Factor in terms of sensible capacity per degree of temperature difference.

ULF=CapacityRated,Sensible/DT1Rated

The total capacity is the sum of the sensible and latent capacity. The sensible heat ratio (SHR) is the sensible heat transfer divided by the total (sensible plus latent) heat transfer. Again, from Nelson, (Nelson, B.I., 2010)

The mass transfer process is much more “thermally effective” than the sensible heat transfer process, that is, the heat flux through the evaporator surfaces during the mass transfer process is extremely high.(AHRI, 2001) Consequently, if the surface effectiveness of the coil were to remain constant, the increase in the evaporator cooling capacity during combined sensible and latent cooling would be equal to the sensible cooling capacity divided by the SHR… However, the increase in heat flux through the fin surfaces has the effect of decreasing fin efficiency and overall surface effectiveness due to an increase in the fin surface temperature gradient.7 The result is a slightly lower total cooling capacity.

Qideal=qsensSHR;SHR=qsensQtotal

Where:

Qideal = Cooling capacity (total) if fin efficiency and total effectiveness were constant, W

QTotal = Cooling capacity (total), actual

The total capacity is therefore a function of the sensible heat ratio, which is a function of the total capacity, and they are both, of course a function of the psychometrics of the air flowing through the chiller. This is handled with a two step estimation process.

ΔT=Minimum(ΔTmax,(TCoilinletTevap))qsens,max=ULFΔT×(1SCHDefrost,DripDown)×SCHCoilTCoilexitestimate=TCoilinletqsens,max˙mDryAir×cp,CoilInletDryAirhCoilexitestimate=f(TCoilexitestimate,PBarometric)ataRelativeHumidityof1.0QTotalestimate=(hCoilInlethCoilexitestimate)×˙mmaxSHR=qsens,maxQTotalestimateCorrection=f(SHR);Functioninputbyuser,linearorquadraticcurveQTotal=Correction×qsens,max

Where:

ΔT = Temperature difference between the inlet air and the average evaporating temperature, C

ΔTMax = Maximum temperature difference specified by the user, C

SCHCoil = Coil availability schedule

hCoil exit = Enthalpy of air at the coil exit

hCoil inlet = Enthalpy of air at the coil inlet

PBarometric = Barometric air pressure,Pa

The “Correction” function must be obtained from the chiller manufacturer. Some curves typical of ammonia chillers have been published (see, for example, Fig. 2 in (Nelson, B.I., 2010)). A default linear approximation of this curve is provided as an input option.

European Standard Ratings[LINK]

Five standard rating conditions have been defined in a European rating system. The capacity is reported at the rating condition as either the “Nominal” or “Standard” capacity. The “Nominal” capacity includes both latent and sensible loads and the “Standard” capacity includes sensible loads only. “Wet Coil Factors” are provided with the ratings to translate between the two, along with a chart giving the impact of Air Inlet Temperature on the Wet Coil Factor. The user identifies the rating condition used and whether the capacity input is “Nominal” or “Standard”. These rating factors, along with the air inlet temperature and evaporating temperature are used to calculate the actual cooling capacity.

QTotal=QNominal×WetCoilFactor(TCoilinlet)WetCoilFactor(StandardCondition)×ΔTΔTRated

Total Capacity Map[LINK]

Some manufacturers are beginning to provide more comprehensive performance information. For these air chillers, the manufacturers specify a Rated Total Capacity at a given inlet air relative humidity. A table or set of curves is then provided to calculate the total capacity QTotal, as a function of the inlet air temperature and relative humidity, and the average evaporating temperature.

Sensible and Latent Capacity[LINK]

The sensible and latent loads served are then calculated as:

hCoilexit=hCoilInletQTotal˙VAir,Max×ρCoilInletTCoilexit=f(hCoilexit)ataRelativeHumidityof1.0HRCoilexit=f(TCoilexit,hCoilexit)˙mWater=˙mdryair,max×(HRCoilexitHRCoilinlet)qlatent=˙mwater×hicetovaporqsens=QTotalqlate