General

This is the multi-zone office hydronic simple emulator model of BOPTEST, emulating a 2-zone building. Each zone is equipped with independent air handling units (AHUs) for ventilation and circuits connected to fan-coil units (also known as ventiloconvectors) as the emission system. Hot and cold water serving the fan-coil units is produced by means of an air-source heat pump and an air-cooled chiller.

Building Design and Use

Architecture

The building test case is modeled as a rectangular building of 40x25 m and 15 m of height. The height of each floor is considered to be 3 m, and at each floor internal floors are added to account for the thermal mass. Hence, the total floor area to be conditioned is 5000 m2. The building is divided in two zones/spaces of equal floor area, with their main façades oriented towards north (NZ) and south (SZ) respectively, as shown in the figure below. Each zone has a window-to-wall ratio of 50%.

\"Floor

Constructions

Exterior walls

The walls are modelled using IDEAS.Buildings.Components.OuterWall and consist of the following layers:

Name

Thickness [m]

Thermal Conductivity [W/(m·K)]

Specific Heat Capacity [J/(kg·K)]

Density [kg/m3]

Layer 1 (external masonry)

0.08

0.89

800

1920

Layer 2 (glasswool)

0.10

0.035

800

60

Layer 3 (interior masonry)

0.14

0.3

880

850

Layer 4 (gypsum)

0.015

0.38

840

1120

The average U-value is 0.28 W/(m2K)

Floor

The floor is modelled using IDEAS.Buildings.Components.SlabOnGround IDEAS.Buildings.Components.SlabOnGround and consists of the following layers:

Name

Thickness [m]

Thermal Conductivity [W/(m·K)]

Specific Heat Capacity [J/(kg·K)]

Density [kg/m3]

Layer 1 (concrete)

0.1

1.4

900

2240

Layer 2 (polyurethane foam)

0.07

0.025

1500

40

Layer 3 (screed)

0.05

1.3

1000

2000

Layer 4 (tile)

0.01

1.4

840

2100

The average U-value is 0.32 W/(m2K)

Roof

The roof is modelled using IDEAS.Buildings.Components.OuterWall and consist of the following layers:

Name

Thickness [m]

Thermal Conductivity [W/(m·K)]

Specific Heat Capacity [J/(kg·K)]

Density [kg/m3]

Layer 1 (screed)

0.02

1.3

1000

2000

Layer 2 (glasswool)

0.10

0.035

800

60

Layer 3 (plywood)

0.02

0.15

1880

540

The average U-value is 0.32 W/(m2K)

Internal walls

The internal walls that separate the zones in the building are modelled using IDEAS.Buildings.Components.InternalWall and consist of the following layers:

Name

Thickness [m]

Thermal Conductivity [W/(m·K)]

Specific Heat Capacity [J/(kg·K)]

Density [kg/m3]

Layer 1 (gypsum)

0.015

0.38

840

1120

Layer 2 (interior masonry)

0.07

0.3

880

850

Layer 3 (gypsum)

0.015

0.38

840

1120

The average U-value is 1.78 W/(m2K)

Internal floors

The internal floors that separate each zone floors in the building are modelled to add thermal mass using IDEAS.Buildings.Components.InternalWall and consist of the following layers:

Name

Thickness [m]

Thermal Conductivity [W/(m·K)]

Specific Heat Capacity [J/(kg·K)]

Density [kg/m3]

Layer 1 (tile)

0.01

1.4

840

2100

Layer 2 (concrete)

0.25

1.4

900

2240

Layer 3 (tile)

0.01

1.4

840

2100

The average U-value is 2.26 W/(m2K)

Windows

The windows are modelled using IDEAS.Buildings.Components.Window and the glazing information of IDEAS.Buildings.Data.Glazing.Ins2Ar2020, with an U-value of 1.028 W/(m2K) and a g-value of 0.551. The window model assumes that the frame occupies 15% of the area and is made of insulated aluminium, with a U-value of 2.8 W/(m2K).

Occupancy schedules

The design occupancy density is one occupant per 15 m2. The number of occupants present in each zone is equally divided and takes the following normalized profile for one day.

\"Normalized

The occupied time for the HVAC system is between 7 AM and 7 PM during the weekdays. The unoccupied time is outside of this period. During summer months (July and August), the occupancy is reduced by half since it is assumed people take holiday and the work load is reduced. A bank holiday schedule is also implemented according to the Belgian calendar for the following days:

Internal loads and schedules

Internal gains from occupancy are taken as 45 W of latent heat and 73 W of sensible heat for typical office work, with a convective-radiative split of 40/60% respectively.

Internal gains from lighting are taken for office lighting requirements of 500 lx (based on standard EN 12464-1) and LEDs with an efficiency of 150 lm/W and with a convective-radiative split of 65/35% respectively.

Internal gains from appliances are assumed to have a convective-radiative split of 70/30% respectively and are distributed as follows:

these values are taken from ASHRAE appliances data in ANSI/ASHRAE Fundamentals 2017.

Climate data

The weather data is from TMY3 for Uccle, Brussels, Belgium, between the years 2007 and 2021. The weather file is hosted in IDEAS.

HVAC System Design

Primary and secondary system designs

The emission system (ventiloconvectors) nominal heating and cooling powers are determined by using the following conventional rules for design.

As such, the nominal powers of the emission system are 75 and 100 kW for heating and cooling respectively for the north zone and 62.5 and 112.5 kW for heating and cooling respectively for the south zone. The mass flow rates of the emission loops are determined using the relation:

Qdesign = mflow cp ΔT

where Qnominal is the nominal heat flow, mflow is the nominal mass flow rate, cp is the specific heat capacity of the media and ΔT is the nominal temperature difference. This nominal temperature difference is a design parameter and it is assumed to be 20 K and 5 K for the water stream for heating and cooling respectively, and 10 K for the air stream (both for the heating and cooling coil) in the ventiloconvector. The heat exchange coefficients (UA) of the heat exchangers in the ventiloconvectors are computed by the model from these nominal values using the NTU method.

The ventilation system is composed of one air handling unit (AHU) per zone. These AHUs are equipped with a double mechanical flux system (i.e., a supply and a return fan), a heat recovery system and a heating and cooling coils. No humidifier/dehumidifer system is installed. No variable air volume (VAV) boxes are equipped in the ventilation ducts. However, the fans of the AHU can modulate their pressure head to achieve the required ventilation rate. The nominal supply and return air flow rates are computed based on the volume of the zones to be conditioned, and they are 3 kg/s and 2.5 kg/s of air respectively. Based on these nominal air flows, the nominal powers and nominal mass flow rates of the heating and cooling coils are calculated using the same relations as for the emission system. The nominal temperature differences in this case are 20 K and 5 K for the water streams for heating and cooling respectively, and 30 K and 15 K for the air stream for heating and cooling respectively. As a result, for each individual AHU the nominal heating power is 90 kW whereas the nominal cooling power is 45 kW. The heat recovery exchanger is assumed to have a constant effectiveness of 84%.

The production system is composed of an air-source heat pump for heating and an air-cooled chiller for cooling. The nominal powers of the production system are determined by the sum of the nominal powers of each emission and ventilation coil loop, which are 75+62.5+90+90=317.5 kW and 100+112.5+45+45=302.5 kW for heating and cooling respectively. The nominal mass flow rates are calculated from the resulting heating and cooling nominal powers.

The production system is connected to the emission and ventilation systems through the distribution system shown in the schematic below. The production system supplies hot and cold water to supply manifolds or collectors. From here, the different circuit loops distribute the water to the different heat exchangers in the emission system or the ventilation system. Each circuit loop is equipped with a pump that allows for activating the circuit and a mixing 3-way valve placed downstream of the pump, which allows for regulating the supply temperature in the distribution loop.

In addition, each heat exchanger is internally equipped with a 3-way valve that allows for regulating the heating or cooling load. However, this component should be seen as an internal part of the heat exchanger and not the distribution system.

\"Building

Equipment specifications and performance maps

Heat pump

The coefficient of performance (COP) of the heat pump is dependant on the supply temperature imposed and the ambient temperature, and takes the following bi-linear relation:

COP = -15.11 - 0.05 Tsupply + 0.125 Tambient

where Tsupply and Tambient are in Kelvin. In essence, the heat pump nominal COP (at 35 °C supply and 7 °C ambient temperature) is 4.5.

Chiller

The energy efficiency ratio (EER) of the chiller is dependant on the supply temperature imposed and the ambient temperature, and takes the following bi-linear relation:

EER = -68.5 + 0.4 Tsupply - 4/30 Tambient

where Tsupply and Tambient are in Kelvin. In essence, the chiller nominal EER (at 7 °C supply and 35 °C ambient temperature) is 2.5.

Fluid movers

The supply and return fans of the air handling units have a nominal pressure drop of 250 Pa and 180 Pa respectively, and follow the pressure curve given by IDEAS.Examples.PPD12.Data.FanCurvePP12 with adapted motor effiency such that their power at nominal speed is 3.50 and 2.15 kW respectively.

The ventiloconvector fans follow the generic curve of the model and have their nominal pressure drop is calibrated such that at nominal mass flow rate their power use is 2.5 kW, based on data from JAGA manufacturer.

The pumps of the distribution and production system follow the generic curve of the model and have the following nominal pressure heads and powers:

Rule-based or local-loop controllers (if included)

The model is equipped with a BMS baseline controller that attempts to keep the temperature and IEQ requirements, also depicted in the hydraulic schematic. Internal controls of the emission and ventilation system are not communicated/exported to the BMS, hence cannot be overridden.

A minimum and maximum zone set-points can be defined by the user, which are 21/25 °C when the building is occupied. When the building is empty, a night setback of ±5 °C is applied. These set-points are input to the ventiloconvectors, whose internal controls (not connected to the BMS) will actuate on their fans and internal 3-way valves to keep the building within the defined boundaries. The distribution pumps connected to the emission system are activated following a schedule control: the heating/cooling pumps are activated every working day from 6AM (to allow pre-conditioning of the space) and de-activated at 7PM, except from Mondays where the activation time is 4AM to compensate for the weekend.

The AHUs are controlled based on the building occupancy. Whenever the building is occupied, the AHUs supply and extract air at their nominal flow rate. The baseline controller does not have any feedback from the CO2 measurements in the zones. The distribution pumps connected to the ventilation system are activated following the same schedule. A constant air supply set-point of 21 °C is set by the baseline controller. The internal controls of the AHUs (not connected to the BMS) manage the internal 3-way valves in the heating/cooling coils and the recovery by-pass to achieve the desired air supply set-point. The by-pass control of the AHUs is connected to the BMS, unlocking the possibility of free cooling during nights. However, this feature is not implemented in the baseline controller.

The mixing 3-way valves in the distribution system modulate to track a given supply set-point. For the ventiloconvectors, the set-points follow a heating/cooling curve based only in the outdoor temperature (no room compensation is included).

For the ventilation coils, these set-points are set constant to: The modulation is achieved by means of PI controllers. All PIs have a constant gain of 0.2 and an integral gain of 30s. The model allows to override either the set-point to be tracked by the PI or the valve position directly.

The production system ideally tracks a given supply temperature set-point. The heating production system follows the maximum of the heating distribution set-points. The cooling production system follows the minimum of the cooling distribution set-points. The pumps in the production system are activated whenever at least one of the distribution pumps is activated.

Model IO's

Inputs

The model inputs are:

Outputs

The model outputs are:

Forecasts

The model forecasts are:

Additional System Design

Lighting

Lighting heat gain is included in the zone model whenever there are occupants in the building and it is not controllable.

Shading

There is no shading on this building.

Onsite Generation and Storage

There is no onsite generation or storage on this building site.

Model Implementation Details

Moist vs. dry air

A moist air model is used. Relative humidity is tracked based on latent heat gain from occupants, outside air relative humidity, and cooling coil models in the fan coil units and the AHU that includes condensation. Since the AHU does not include a humidifier/dehumidifer system, this parameter is not subjected to control.

Pressure-flow models

The duct airflows and pipe waterflows are modeled using a pressure-flow network. Air exchange between zones is not modeled.

Infiltration models

Airflow due to infiltration is calculated using the IDEAS.Buildings.Components.InterzonalAirFlow.n50FixedPressure model and a n50 value of 5.

CO2 models

CO2 generation in the zones is calculated using the IDEAS.Buildings.Components.Occupants and IDEAS.Buildings.Components.OccupancyType.OfficeWork models. Outside air CO2 concentration is 400 ppm.

Scenario Information

Time Periods

The Peak Heat Day (specifier for /scenario API is 'peak_heat_day') period is:

The Typical Heat Day (specifier for /scenario API is 'typical_heat_day') period is:

The Peak Cool Day (specifier for /scenario API is 'peak_cool_day') period is:

The Typical Cool Day (specifier for /scenario API is 'typical_cool_day') period is:

The Mix Day (specifier for /scenario API is 'mix_day') period is:

Energy Pricing

The gas price is considered constant at a value of 0.0464 EUR/kWh, corresponding to the average prices for Belgian buildings with an energy use between 20 and 500 MWh/year for the year 2015. (Eurostat. Gas prices for industrial consumers - bi-annual data (from 2007 onwards). Tech. rep., 2016. and Eurostat. Electricity prices components for industrial consumers - annual data (from 2007 onwards). Tech. rep., 2017.)

The Constant Electricity Price (specifier for /scenario API is 'constant') profile is:

The Dynamic Electricity Price (specifier for /scenario API is 'dynamic') profile is:

The Highly Dynamic Electricity Price (specifier for /scenario API is 'highly_dynamic') profile is:

Emission Factors

It is used a constant emission factor for gas of 0.22 kgCO2/kWh as reported by the Centre for Environmental Economics and Environmental Management (CEEM - Centre for Environmental Economics and Environmental Management. Universiteit Gent. https://ceem.ugent.be/en/index.htm.).

It is used a hourly variable emission factor for electricity, extracted from Electricity Maps for the year 2019.