## Engineering Reference — EnergyPlus 9.6

### Engineering Reference

Water systems include a variety of components that simulate water consumption, production, and storage, including:

1. Water Use Equipment and Connections

2. Rainwater Collection

3. Groundwater Well

4. Water Storage Tanks

The input object Site:WaterMainsTemperature provides a model for the temperature of water mains that deliver water to a building via underground pipes. The water mains temperatures are an important input for plant equipment that require make-up water from the water mains, such as for domestic hot water systems.

Water mains temperatures are a function of outdoor climate conditions and vary with time of year. A correlation has been formulated to predict water mains temperatures based on two weather inputs:

• average annual outdoor air temperature (dry-bulb)

• maximum difference in monthly average outdoor air temperatures

These values can be easily calculated from annual weather data using a spreadsheet.

The correlation was developed by Craig Christensen and Jay Burch and is described in Hendron et al. (2004). The correlation is:

Tmains=(Tout,avg+6)+ratio(ΔTout,maxdiff/2)SIN(0.986(day15lag)90)

where:

TTmains is the water mains temperature (F)

TTout,avg is the average annual outdoor air temperature (F)

ΔTout,maxdiff is the maximum difference in monthly average outdoor air temperatures (F)

day is Julian day of the year (1-365)

ratio is equal to 0.4 + 0.01 * (Tamb,avg - 44)

lag is equal to 35 - 1.0 * (Tamb,avg - 44) (F).

For the Chicago-O’Hare TMY2 weather file, Tout,avg = 9.69C and ΔTout,maxdiff = 28.1C. The resulting water mains temperature profile is graphed in Figure 1.

Hendron, R., Anderson, R., Christensen, C., Eastment, M., and Reeves, P. 2004. “Development of an Energy Savings Benchmark for All Residential End-Uses”, Proceedings of SimBuild 2004, IBPSA-USA National Conference, Boulder, CO, August 4 - 6, 2004.

## Water Use Equipment and Connections[LINK]

Essential parts of most water systems are the end uses where the water is utilized for a purpose. These end uses are characterized by familiar pieces of equipment such as sinks, showers, dishwashers, clothes washers, and toilets. End uses can also include outdoor uses such as landscape irrigation and agricultural irrigation.

Water end-use modeling in EnergyPlus is accomplished via two input objects:

The WaterUse:Equipment object simulates all different types of water end uses. The WaterUse:Connections object manages multiple pieces of water equipment with a common supply and return, internally providing the services of a splitter and a mixer. The WaterUse:Connections object also allows water consuming equipment to be simulated in a closed plant loop by supplying makeup water from the water mains.

Common water equipment, such as sinks and showers, requires mixing of hot and cold water to provide a desired temperature at the tap. Water use equipment that is simulated in “unconnected” mode, i.e., without being referenced by a WaterUse:Connections object, presents a relatively simple modeling problem illustrated in Figure 2.

The variables are defined as:

˙mhot is the hot water supply mass flow rate

Thot is the hot water supply temperature

˙mcold is the cold water supply mass flow rate

Tcold is the cold water supply temperature

˙mtarget is the target mass flow rate at the tap

Ttarget is the target temperature at the tap.

The basic physics of the problem are described by two equations representing conservation of mass and conservation of energy respectively:

˙mtarget=˙mhot+˙mcold

˙mtargetTtarget=˙mhotThot+˙mcoldTcold

The hot and cold water supply temperatures, Thot and Tcold, are specified by the user with schedules (or indirectly from the water mains temperatures). The target flow rate and target temperature, ˙mtarget and Ttarget, at the tap are also specified with schedules. The target conditions can be thought of as a way to simulate occupant behavior; a person at a sink will adjust the hot and cold flow rates to get the overall water flow and temperature that is desired.

The two equations above are easily solved to calculate the necessary hot and cold flow rates needed to attain the target conditions.

˙mhot=˙mtarget(TtargetTcoldThotTcold)

˙mcold=˙mtarget˙mhot

Even though hot and cold flow rates are unlimited in “unconnected” mode, it is still possible to fail to meet the target conditions if Ttarget > Thot. In this case, the actual mixed water temperature at the tap, Tmixed, is set equal to Thot. The target flow rate is always met.

Water equipment that omits schedules for the target temperature and/or hot water suppy temperature implies that no hot water is needed. The result is all cold water at the target flow rate.

For “unconnected” water equipment, the heating rate and energy that is required to supply the hot water is calculated by the following equations.

Q=˙mhotcp(ThotTcold)

ΔE=QΔt

where:

Q is the heating rate

cp is the specific heat of water

ΔE is the heating energy

Δt is the time step interval.

All heating is assigned to “Purchased Heating”.

## Zone Heat Gain from Water Use Equipment[LINK]

Some water use equipment can be a significant source of sensible and latent heat gain when located in a zone. Showers, for example, add a noticeable latent gain as some of the hot water evaporates into the zone air (and at the same time adding some sensible gain). Other types of equipment may contribute greater sensible gain but less latent gain.

Calculating the amount of sensible and latent heat gain can be extremely complicated and depends on many variables. One set of variables is related to the initial conditions of the entering water and the undisturbed zone air, namely:

• Entering water temperature

• Entering water flow rate

• Ambient zone air temperature

• Ambient zone humidity

• Ambient zone barometric pressure

A second set of variables is tied to the specifics of the water equipment which describe the heat and mass transfer processes that occur as the water passes through the zone. These variables might include the velocity of the water in the zone air, the surface area of water exposed to the zone air, the dwell time of water in the zone, and possibly others. The complexity of these variables cannot be easily cast into a generalized EnergyPlus input object. Consider that in a shower, the velocity of water droplets is increasing as they fall, and the surface area of the spray is increasing. How does one even determine the surface area of the spray at any point?

The approach taken here is to utilize the first set of initial condition variables and forego the overly-complex second set. The initial conditions can be used to calculate a maximum possible change in sensible and latent energy between the entering water conditions and the leaving water conditions at an infinite time later. Sensible and latent energy are calculated separately and represent the maximum heat gain, respectively, that could be added to the zone air. The user simply specifies a fraction of the maximum heat gain for sensible and latent that is actually added to the zone air. The fraction is specified with a schedule to account for different modes of operation. The split between sensible and latent will vary depending on the type of equipment that is to be modeled. Typically, both fractions should be small numbers.

The maximum sensible heat gain rate is calculated by assuming that all of the water flow is allowed to cool to the zone dry-bulb temperature. The actual sensible zone heat gain rate is found by multiplying the maximum by the user fraction.

Qsens=fsens˙mtargetcp(TtargetTzone)

ΔEsens=QsensΔt

where:

Qsens is the sensible heat gain rate

fsens is the user fraction of maximum sensible heat

˙mtarget is the target mass flow rate at the tap

cp is the specific heat of water

Ttarget is the target temperature at the tap

Tzone is the zone dry-bulb air temperature

ΔEsens is the sensible heat gain

Δt is the time step interval.

The maximum latent heat gain rate is calculated by assuming that the maximum amount of water is evaporated into the zone air during the time step. The actual latent zone heat gain rate is found by multiplying the maximum by the user fraction.

mevap,zonemax=(WsatWzone)ρzoneVzone

mevap,flowmax=˙mtargetΔt

mevap,max=MIN(mevap,zonemax,mevap,flowmax)

mevap=flatmevap,max

˙mevap=mevap/Δt

Qlat=hfg˙mevap

ΔElat=QlatΔt

where:

mevap,zonemax is the mass of water that can be absorbed by the zone air in one timestep

Wsat is the fully-saturated zone humidity ratio (function of air temperature and pressure)

Wzone is the zone humidity ratio

ρzone is the zone dry-air density (function of air temperature and pressure)

Vzone is the zone air volume

mevap,flowmax is the mass of water available in the flow in one timestep

˙mtarget is the target mass flow rate at the tap

Δt is the time step interval

mevap,max is the maximum mass of water that can be evaporated

mevap is the actual mass of water evaporated into the zone air

flat is the user fraction of maximum latent heat

˙mevap is the mass rate of water evaporation into the zone air

Qlat is the latent heat gain rate

hfg is the heat of vaporization for water

ΔElat is the latent heat gain.

Sensible and latent heat gains to the zone are lagged by one zone time step.

Note that this method allows the heat gain to vary dynamically according to changing water and zone conditions instead of the fixed or scheduled value that is traditionally used. However, sensible and latent fractions may be difficult to estimate. Experimental measurement is probably the best way to determine the fractions.

A related result of any zone heat and moisture gain is that the water undergoes an equal heat and moisture loss. For a shower the water will generally leave the zone at a cooler temperature and slightly reduced flow rate when compared to the entering water conditions.

˙mdrain=˙mtarget˙mevap

Tdrain=˙mtargetcpTtargetQsensQlat˙mdraincp

where:

˙mdrain is the drainwater mass flow rate at the drain

˙mtarget is the target mass flow rate at the tap

˙mevap is the mass rate of water evaporation into the zone air

Tdrain is the water temperature at the drain

cp is the specific heat of water

Ttarget is the target temperature at the tap

Qsens is the sensible heat gain rate

Qlat is the latent heat gain rate.

The change in leaving water conditions has implications for drainwater heat recovery and storage tank conditions.

Water use equipment referenced by a WaterUse:Connections object is considered to be in “connected” mode. The WaterUse:Connections object manages multiple pieces of water equipment with a common supply and return, internally providing the services of a splitter and a mixer. The connections object also offers three internal configurations for drainwater heat recovery. External connections are available for coupling to plant loop nodes or water storage tanks for graywater reuse or reclamation.

Keep in mind that water use equipment in “connected” mode does not imply that it is connected to a plant loop. WaterUse:Equipment objects are never directly connected to a plant loop; they are connected indirectly through the WaterUse:Connections object. WaterUse:Connections can operate either within a plant loop or in stand-alone mode, i.e., without plant loop connections.

The WaterUse:Connections object overrides the hot and cold water supply temperatures of its subordinate WaterUse:Equipment objects based on the specified schedules, plant loop conditions, or water storage tank conditions. But the WaterUse:Equipment objects still maintain their individual target flow rates and target temperatures.

Solving the water subsystem described by the WaterUse:Connections object is appreciably more difficult than solving a single, unconnected piece of water use equipment. The subsystem is illustrated in Figure 4.

The nomenclature here uses uppercase letters to indicate variables related to the connections object and lowercase letters to indicate variables related to the equipment objects. The variables are defined as:

˙Mhot is the hot water supply mass flow rate

Thot is the hot water supply temperature

˙Mcold is the cold water supply mass flow rate

Tcold is the cold water supply temperature

˙Mdrain is the drainwater mass flow rate at the drain

Tdrain is the water temperature at the drain

Twaste is the wastewater temperature after the drain

˙Mreturn is the return water mass flow rate back to the plant loop

Treturn is the return temperature of makeup water back to the plant loop.

The solution method proceeds in steps:

• Initialize Connections Temperatures

• Calculate Connections Flow Rates

• Calculate Connections Drain Temperature

• Update Connections Nodes

Initializes the hot and cold supply water temperatures, Thot and Tcold, for the common supply to all the equipment. If plant-coupled, Thot is taken from the plant loop inlet node. If stand-alone, Thot is taken from the schedule.

Tcold is taken from the supply water storage tank, if specified. Otherwise Tcold is taken from the schedule. If a schedule is not specified, Tcold is taken from the water mains.

First calls on each piece of equipment to simulate itself with the given supply water temperatures Thot and Tcold. The desired hot and cold flow rates for ˙mhot,i and ˙mcold,i are calculated (as described earlier) and passed back to the connections object. The total flow rates for ˙Mhot and ˙Mcold are then calculated:

˙Mhot=i˙mhot,i

˙Mcold=i˙mcold,i

If plant-coupled, ˙Mhot is compared to the maximum flow rate available in the plant loop (after the first HVAC iteration). If ˙Mhot > ˙Mhot,max, the actual flow rate is reset to equal the maximum: ˙Mhot = ˙Mhot,max. New hot and cold flow rates, ˙mhot,i and ˙mcold,i, at the equipment level are then recalculated, decreasing the hot flow rate and increasing the cold flow rate in order to still meet the target flow rate.

α=˙Mhot,max˙Mhot

˙mhot,i=α˙mhot,i

˙mcold,i=˙mtarget,i˙mhot,i

A new mixed water temperature is also recalculated at the equipment level.

Although water storage tanks can also have a maximum flow rate on the cold side, the cold flow rate is not limited in this implementation.

At this point zone heat gains are calculated for all equipment in zones. The final drainwater temperatures and flow rates are determined for each piece of equipment by subtracting the lost heat and moisture (see above). The total drainwater temperature and flow rate for the subsystem are calculated:

˙Mdrain=i˙mdrain,i

Tdrain=i˙mdrain,itdrain,i˙Mdrain

In the case of no drainwater heat recovery, the subsystem wastewater temperature, Twaste, is equal to the drainwater temperature, Tdrain. (For drainwater heat recovery, see below.) The wastewater temperature and flow rate are propogated to the reclamation water storage tank, if specified.

Finally, if plant-coupled, the return water temperature, Treturn, is passed on to the plant outlet node. Treturn is taken from the cold water supply schedule or the water mains. The return flow rate is equal to the hot water supply flow rate in order to preserve a closed plant loop: ˙Mreturn = ˙Mhot.

For “connected” water equipment, the heating rate and energy that is required to supply the hot water for individual water equipment is calculated by the following equations.

Q=˙mhotcp(ThotTreturn)

ΔE=QΔt

where:

Q is the heating rate

cp is the specific heat of water

ΔE is the heating energy

Δt is the time interval.

The heating rate and energy for the WaterUse:Connections is the sum of the values for its constituent water equipment objects. If the WaterUse:Connections object is stand-alone, all heating is assigned to “Purchased Heating” by the individual water equipment objects. If the WaterUse:Connections object is coupled to a plant loop, all heating is supplied by the plant loop equipment.

Drainwater heat recovery adds one more complication to the water connections subsystem. A heat exchanger is added between the drain and the waste outlet to recover some heat from the water before it leaves the system. The heat is usually recovered by preheating the incoming makeup water from the mains.

The heat exchanger can be modeled as “ideal”, “counterflow”, or “crossflow”. One common type of drainwater heat exchanger is the gravity-film exchanger (GFX). The GFX can be approximated with a counterflow heat exchanger.

The destination of the preheated makeup water can be plumbed in three possible configurations: “plant”, “equipment”, “plant and equipment”. In the “plant” configuration, all preheated water flow is returned to the plant loop to match the makeup water for the hot water supply flow. In the “equipment” configuration, all preheated water flow is directed internally within the WaterUse:Connections object to provide the makeup water for the cold water supply flow for the local water use equipment. In the “plant and equipment” configuration, the preheated water is split between both of the previous configurations. This is the only configuration where the flow rate is equal on both sides of the heat exchanger.

The new variables added here are defined as:

˙Mrecovery is the makeup water mass flow rate in the heat exchanger

Trecovery is the water temperature leaving the heat exchanger

Tmakeup =is the cold makeup water temperature supplied to various flows.

Another step is inserted into the solution method:

• Initialize Connections Temperatures

• Calculate Connections Flow Rates

• Calculate Connections Drain Temperature

• Update Connections Nodes

• Calculate Connections Heat Recovery

• Update Connections Nodes

For the “equipment” and “plant and equipment” configurations, the solution requires iteration to solve because the preheated water leaving the heat exchanger is connected to the cold water supply water, creating a feedback loop.

Based on the selected configuration, the flow rate through the heat exchanger is first determined. For the “plant” configuration, ˙Mrecovery=˙Mhot. For the “equipment” configuration, ˙Mrecovery=˙Mcold. For the “plant and equipment” configuration, ˙Mrecovery=˙Mhot+˙Mcold.

The heat capacity rates on both sides of the heat exchanger are calculated and the minimum and maximum heat capacity rates, Cmin and Cmax, are determined.

Cdrain=˙Mdraincp

Crecovery=˙Mrecoverycp

Cmin=MIN(Cdrain,Crecovery)

Cmax=MAX(Cdrain,Crecovery)

Rc=CminCmax

where Rc is the heat capacity ratio.

The effectiveness ε is then calculated for the given type of heat exchanger. For the “ideal” heat exchanger, ε=1.

For the “counterflow” and “crossflow” heat exchangers, the effectiveness is calculated using the effectiveness-NTU method, where:

NTU=UACmin

For the “counterflow” heat exchanger:

ε=NTU1+NTU, if Rc = 1

For the “crossflow” heat exchanger:

ε=1EXP(NTU0.22Rc(EXP(RcNTU0.78)1))

The heat recovery rate is then calculated.

Qrecovery=εCmin(TdrainTmakeup)

The outlet water temperatures can then be calculated for both sides of the heat exchanger.

Trecovery=Tmakeup+Qrecovery˙Mrecoverycp

Twaste=TdrainQrecovery˙Mdraincp

At each iteration the difference between Trecovery and Tcold is compared and checked for convergence.

ΔTerror=TrecoveryTcold

Once the temperatures have converged, the solution is complete. Twaste is then passed on to a storage tank, if specified.