Air Exchange[LINK]

Air exchange and interchange between zones is treated as a convective gain. Temperature-difference-controlled or constant air mixing can be specified as a one-way or cross-zone phenomenon modeled using the ZoneMixing or ZoneCrossMixing objects. Air exchange through doorways between refrigerated spaces can be modeled using the ZoneRefrigerationMixing object.

For one-way mixing (using ZoneMixing object(s)), the mixing air flow is only used for the energy and mass balance for the receiving zone. The mass referred to in this section includes air, water vapor and CO${}_2$. The source zone energy and mass balance are not affected, although the user may choose to enter complementary pairs of one-way mixing objects. Multiple mixing flows can be defined for any receiving zone. For cross-zone mixing (using ZoneCrossMixing object(s)), the mixing air flow impacts the mass and energy balances in both the source and receiving zones. A separate ZoneCrossMixing object must be used for each of the two zones exchanging air if the mixing flow is bi-directional and based on a temperature difference greater than zero.

For refrigerated space air exchange (using ZoneRefrigerationDoorMixing object(s)), the mixing air flow impacts the mass and energy balances in both the source and receiving zones. A single object accounts for the two-way air flow with the energy and mass exchanges determined by the air density difference between the two zones.

Temperature Difference Controlled Air Exchange[LINK]

The volume of air flow into the receiving zone is specified by the user with a number of control parameters and schedules listed in the Input Output Guide. The user can turn this one-way flow on or off based on the temperature difference between the source and receiving zones, or it may be independent of the temperature difference. The density and specific heat of the air flowing into the receiving zone are determined using the average temperature and average humidity ratio in the source and receiving zones. The humidity ratio of the air flowing into the receiving zone is set equal to the humidity ratio of the source zone. The mass, moisture,and energy terms are then used as described in two previous sections, Basis for the Zone and Air System Integration, and Moisture Predictor-Corrector.

$${\rho }_{Avg}=f(\frac{T_{ReceivingZone}+T_{SourceZone}}{2},\frac{W_{ReceivingZone}+W_{SourceZone}}{2},P_{Barometric})$$

$$c_{P,Avg}=f(\frac{T_{ReceivingZone}+T_{SourceZone}}{2},\frac{W_{ReceivingZone}+W_{SourceZone}}{2})$$

$${\dot{m}}_{MixingFlowToReceivingZone}=\sum _{AllSourceZones}{\rho }_{Avg}{\dot{V}}_{Air}$$

$${\dot{Q}}_{MixingFlowToReceivingZone}=\sum _{AllSourceZones}{\rho }_{Avg}{C}_{p,Avg}{\dot{V}}_{Air}(T_{sourceZone}-T_{receivingZone})$$

$${Moisture}_{MixingFlowToReceivingZone}=\sum _{AllSourceZones}{\rho }_{Avg}{\dot{V}}_{Air}{W}_{SourceZone}$$

where

c${}_{P,Avg}$ is the average specific heat of air within the two zones (J/kg.K)

${\dot{m}}_{MixingFlowToReceivingZone}$ is the mass of moist air flowing into the receiving zone (kg${}_{air}$/s)

Moisture${}_{MixingFlowToReceivingZone}$ is the moisture mass flow rate into the receiving zone (kg${}_{H2O}$/s)

P${}_{Barometric}$ is the outside barometric pressure (Pa)

${\rho }_{Avg}$ is the average density of air within the two zones (kg/s)

${\dot{Q}}_{MixingFlowToReceivingZone}$ is the energy added to receiving zone air by mixing mass flow (W)

T${}_{ReceivingZone}$ is the temperature in the receiving zone (${}^o$C)

T${}_{SourceZone}$ is the temperature in the source zone (${}^o$C)

${\dot{V}}_{Air}$ is the volume rate of air flow defined by the user (m${}^3$/s)

W${}_{ReceivingZone}$ is the humidity ratio in the receiving zone (kg${}_{H2O}$/kg${}_{dryair}$)

W${}_{SourceZone}$ is the humidity ratio in the source zone (kg${}_{H2O}$/kg${}_{dryair}$)

For cross-mixing, the mass of moist air exchanged between the two zones is assumed to be equal. Again, the density and specific heat are based on the average conditions in the two zones. Note that the temperature and humidity ratio differences ensure that when the energy and moisture terms are used in the Moisture Predictor-Corrector, they correctly reflect a loss or gain in each zone.

$${\dot{m}}_{MixingFlowToReceivingZone}={\dot{m}}_{MixingFlowToSourceZone}={\rho }_{Avg}{\dot{V}}_{Air}$$

$${\dot{Q}}_{MixingFlowToSourceZone}={\rho }_{Avg}{C}_{p,Avg}{\dot{V}}_{Air}(T_{ReceivingZone}-T_{SourceZone})$$

$${\dot{Q}}_{MixingFlowToReceivingZone}={\rho }_{Avg}{C}_{p,Avg}{\dot{V}}_{Air}(T_{SourceZone}-T_{ReceivingZone})$$

$$Moistur{e}_{MixingFlowToSourceZone}={\rho }_{Avg}{\dot{V}}_{Air}{W}_{ReceivingZone}$$

$$Moistur{e}_{MixingFlowToReceivingZone}={\rho }_{Avg}{\dot{V}}_{Air}{W}_{SourceZone}$$

where

${\dot{m}}_{MixingFlowToSourceZone}$ is the mass of moist air flowing into the source zone (kg${}_{air}$/s)

Moisture${}_{MixingFlowToSourceZone}$ is the moisture mass flow rate into the source zone (kg${}_{H2O}$/s)

${\dot{Q}}_{MixingFlowToSourceZone}$ is the sensible energy added to source zone air by mixing mass flow (W)

${\dot{Q}}_{MixingFlowToReceivingZone}$ is the sensible energy added to receiving zone air by mixing mass flow, W

$Moistur{e}_{MixingFlowToSourceZone}$is the latent load added to source zone air by mixing mass flow (kg${}_{H2O}$/s)

$Moistur{e}_{MixingFlowToReceivingZone}$is the latent load added to receiving zone air by mixing mass flow (kg${}_{H2O}$/s)

Density Difference Controlled Air Exchange[LINK]

When closed refrigerated spaces exchange air with other closed spaces, the air flow is determined by the difference in air density between the two spaces. The fundamental assumption for this case is that the mass of dry air exchanged between the two spaces is the same.(Gosney and Olama, 1975] This assumption applies to situations where the colder of the two spaces is essentially sealed to other air flows, that is, there are no open doors or exhaust air flows. Multiple refrigeration door mixing objects can be used for the zone, but if there are multiple doors open at the same time for any significant amount of time, the model will not give results appropriate for that condition.

The sensible and latent energy loads are modeled according to the guidance specified in (ASHRAE 2006d, ASHRAE 2009, and Gosney and Olama, 1975). Equal dry air exchange is assumed, that is, the mass of dry air infiltrating into the receiving zone is assumed to equal the mass of dry air infiltrating out of the source zone.

$$h_{zoneA},{\rho }_{zoneA}=f(T_{zoneA},W_{zoneA},P_{Barometric})$$

$$h_{zoneB},{\rho }_{zoneB}=f(T_{zoneB},W_{zoneB},P_{Barometric})$$

$${\rho }_{zoneA}>{\rho }_{zoneB}$$

$${\dot{Q}}_{Mixing}={\dot{Q}}_{FullFlow}\times {\text{Schedule}}_{DoorOpen}\times F_{Flow}\times (1-F_{Protection})$$

$${\dot{Q}}_{FullFlow}=B(h_{zoneB}-h_{zoneA})$$

$$B=0.221A_{door}{\rho }_{zoneA}{F}_{Density}\sqrt{(1-\frac{{\rho }_{zoneB}}{{\rho }_{zoneA}})g{H}_{Door}}$$

$$F_{Density}={\left[\frac{2}{1+({\rho }_{zoneA}/{\rho }_{zoneB}{)}^{1/3}}\right]}^{1.5}$$

$${\dot{m}}_{DryAirZonesAB}=\frac{{\dot{Q}}_{Mixing}}{h_{zoneB}-h_{zoneA}}$$

$${\dot{m}}_{MixingFlowZoneBtoA}=\sum _{AllZoneBs}{\dot{m}}_{DryAirZonesAB}(1+W_{ZoneB})$$

$${\dot{Q}}_{MixingFlowZoneBtoA}=\sum _{AllZoneBs}{\dot{m}}_{ZoneBtoA}{C}_{p,ZoneB}(T_{zoneB}-T_{zoneA})$$

$${\text{Moisture}}_{MixingFlowZoneBtoA}=\sum _{AllZoneBs}{\dot{m}}_{ZoneBtoA}(W_{zoneB}-W_{zoneA})$$

$${\dot{m}}_{MixingFlowZoneAtoB}=\sum _{AllZoneAs}{\dot{m}}_{DryAirZonesAB}(1+W_{ZoneA})$$

$${\dot{Q}}_{MixingFlowZoneAtoB}=\sum _{AllZoneAs}{\dot{m}}_{ZoneBtoA}{C}_{p,ZoneA}(T_{zoneA}-T_{zoneB})$$

$${\text{Moisture}}_{MixingFlowZoneAtoB}=\sum _{AllZoneAs}{\dot{m}}_{ZoneBtoA}(W_{zoneA}-W_{zoneB})$$

where

A${}_{door}$ is the area of door between Zones A and B (m${}^2$)

F${}_{Flow}$ is the doorway flow factor, = 0.8 if $\Delta$T > 11${}^o$C; = 1.1 if $\Delta$T < = 11${}^o$C

F${}_{Protection}$ is the doorway protection factor, = 0 for no protection; = 0.5 for an air curtain; and 0.9 for a strip curtain (dimensionless)

g is the gravitational constant (m/s${}^2$)

h${}_{ZoneA}$ is the enthalpy of the air within Zone A (J/kg)

h${}_{ZoneB}$ is the enthalpy of the air within Zone B (J/kg)

H${}_{door}$ is the height of door between source and receiving zones (m)

Q${}_{FullFlow}$ is the sensible and latent refrigeration load (on Zone A) for fully established flow (W)

Q${}_{Mixing}$ is the sensible and latent mixing refrigeration load on Zone A for the time step (W)

m${}_{DryAirZoneAB}$ is the mass of dry air exchanged between zones A and B (kg${}_{air}$/s)

Schedule${}_{DoorOpen}$ is the value scheduled by user, fraction of time door open during time step (dimensionless)

W${}_{ZoneA}$ is the humidity ratio of the air within Zone A (kg${}_{H2O}$/kg${}_{air}$)

W${}_{ZoneB}$ is the humidity ratio of air within Zone B (kg${}_{H2O}$/kg${}_{air}$)

P${}_{ZoneA}$ is the density of air within Zone A (kg/m${}^3$)

${\rho }_{ZoneB}$ is the density of air within Zone B (kg/m${}^3$)

References[LINK]

ASHRAE. 2006d. Refrigeration Handbook, Chapter 13. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.

ASHRAE. 2009. Fundamentals Handbook, Chapter 1. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.

Gosney, W.B., Olama, G.A.-L. 1975. Heat and Enthalpy Gains through Cold Room Doorways, Proceedings of the Institute of Refrigeration, vol. 72, pp 31-41.