In this paper, energy consumption analysis and a process to identify appropriate models based on heat dynamics for large structures are presented. The analysis uses data from heating, ventilation, and air-conditioning (HVAC) system sensors, as well as data from the indoor climate and energy software (IDA Indoor Climate and Energy (IDA-ICE) 4.7 simulation program). Energy consumption data (e.g., power and hot water usage) agrees well with the new models. The model is applicable in a variety of applications, such as forecasting energy consumption and controlling indoor climate. In the study, both data-derived models and a grey-box model are tested, producing a complex building model with high accuracy. Also, a case study of the S. J. Carew building at Memorial University, St. John’s, Newfoundland, is presented.
Ministry of Higher Education of Libyan Government1. Introduction
Heating, ventilation, and air-conditioning (HVAC) systems are crucial for indoor climate management and air quality. These systems are also a key factor in overall operational costs. For industrial buildings, nearly one-third of the energy usage depends on HVAC system operation [1–3]. The recent rapid industrialization of the world’s developing nations has led to an increase in energy demand, followed by a rapid rise in pollution levels. As a result, researchers are investigating ways to mitigate or prevent further environmental damage through a combination of conservation methods and wide-scale adoption of renewable energy systems [4].
Ideally, HVAC (heating, ventilation, and air-conditioning) systems are developed to form an interior environment that provides user-comfort with operational cost-efficiency. To maintain consistent user-comfort and affordability amidst changing variables, a suitable control system is needed. Several options have been modeled. One popular method uses data to create a mathematical-based HVAC system that considers input and output variables to find and set system parameters. Data-driven HVAC can readily identify strategies for system refinement and enhancement. These types of model determination are termed system identification (SI) in the literature (ASHRAE, 2005) [5].
In previous studies, researchers categorized modeling approaches into two main types, namely, black box and grey box. For the black box method, no prior information is required, but for the grey-box strategy, there must be a reservoir of preexisting knowledge. Due to these constraints, the black box modeling approach is generally better represented in the literature. Examples of black box models applied to HVAC systems are polynomial forms such as ARX, ARMAX, BJ, and OE. Despite its popularity, the black box model strategy overlooks the physical features of a system, leading to issues around the practical application in real-world designs.
For example, Chi-Man Yiu et al. [6] looked at black box strategies for air-conditioning systems. The researchers contrasted two ARMAX models: the first, a single-input/single-output system, and the second, a multi-input/multi-output (MIMO) system. For the MIMO system, Chi-Man Yiu et al. [6] employed parameters derived from the recursive extended least squares method. Mustafaraj et al. [7] investigated temperature and humidity models (ARX, ARMAX, BJ, and OE) for office environments, using a black box approach. The same researchers [8] continued their work by applying nonlinear autoregressive models with NARX inputs to gauge temperature and humidity levels while comparing and contrasting the outcomes for these models with those of linear ARX models. Additionally, Mustafaraj et al. [8] examined CO2 concentrations’ effect on model performance, considering that occupancy levels in a building are directly correlated to CO2. Rabl [9] provided a summary of approaches applied for dynamic analysis of power usage by modeling heat dynamics. These models were applied in studies by Sonderegger [10] and Boyer et al. [11], using differential equations. For dynamic models, parameter estimation and system identification are essentially the same processes.
This paper uses the grey-box method for modeling dynamic systems. This approach, which is accurate and comprehensive, enables the collection of information on a structure’s thermal properties [12–14]. Grey-box models employ discrete-time measurement equations and continuous time stochastic differential equations. An HVAC system’s yearly power usage can be predicted using energy performance analysis tools, such as SIMBAD, EnergyPlus, eQUEST, HVACSIM+, IDA-ICE, and TRNSYS at set time frames (hourly or less) by a set of equations describing a building’s thermal performance. Calculations comparing various design options are usually made for part-load and full-load performance [15–18].
This paper simulates a whole building (the S. J. Carew building in St. John’s, Newfoundland) using the IDA Indoor Climate and Energy (IDA-ICE) 4.7 simulation program. In addition to examining the modeled structure’s power use, the study investigates a 3D model, a heat model (with variable parameters), and an IDA-ICE model library. The IDA-ICE was developed to investigate different thermal climate zones occurring in indoor environments [19].
By actual details and logged data of a large building in a cold climate, we used the IDA-ICE software and building logged data to model the building in IDA-ICE software. We propose 12 inputs and 12 outputs dynamic model for the system. The dynamic model is required to design and test system controllers before actual implementation. To determine a state-space system model, we use MATLAB system identification toolbox. For the model determination, we used data from IDA-ICE software. Contributions of this paper are building data, proposed system dynamic model, a method to determine the system model, and developed system dynamic model parameters. The primary objectives of this paper are as follows:
Apply the IDA-ICE program to model the S. J. Carew building (Memorial University, St. John’s, Newfoundland) using all real dimensions and building materials information from Department of Facilities Management and the Honeywell office, which is responsible for running and monitoring the system.
By using the IDA-ICE software we can divide the single valve of hot water coming from the main room to four distinct units, enabling each air-handling unit (AHU) to have an individual valve for control the mechanical hot water flow for each zone as another input of the system are supply fan speed and fresh air dampers position.
Compare the data from the IDA-ICE program with the building logged data for validating the power use outcomes.
Determine the potential of applying the system identification approach to reduce the time needed to simulate the building and use the system model simulation results in identifying the dynamics related to a building’s climate control.
2. The Building for This Case Study
A case study on the S. J. Carew building, with an interior size of 25,142 m2, is conducted. The building is located on the campus of Memorial University, St. John’s, Newfoundland, and includes several teaching rooms and research labs for the Memorial’s Faculty of Engineering and Applied Science. The building also features a large cafeteria. There are four individual air-handling units (AHUs) in 300 zones within the building. Figure 1 illustrates a 3D model for the structure, applying the IDA-ICE program mentioned in the previous section, while Table 1 provides an energy report.
Energy report for the building.
Input data Report
Building envelope
Windows
Project
Building
Building
Area [m2]
U [W/(m2 K)]
U∗A [W/K]
% of total
Side
Area [m2]
U Glass [W/(m2 K)]
U Frame [W/(m2 K)]
U Total [W/(m2 K)]
U∗A [W/K]
Customer
Model floor area
25141.7 m2
Walls above ground
7450.19
0.54
4002.24
44.85
N
174.36
2.86
2
2.77
483.51
Created by
Almahdi Abdo-Allah
Model volume
128952.9 m3
Roof
10529.8
0.17
1811.13
20.3
E
131.08
2.86
2
2.77
363.49
Location
Newfoundland (St. John’s Airport)_718010
Model ground area
10544.5 m2
Windows
713.73
2.77
1979.24
22.18
S
213.24
2.86
2
2.77
591.35
Climate file
CAN_NF_St.Johns.718010_CWEC
Model envelope area
29440.0 m2
Doors
201.72
0.59
118.38
1.33
W
195.05
2.86
2
2.77
540.89
Case
building2017_AHU8
Average U-value
0.3031 W/(m2 K)
Walls below ground
0
0
0
0
Total
713.73
2.86
2
2.77
1979.24
Simulated
1/19/2017 22:04
Envelope area per Volume
0.2283 m2/m3
Floor towards ground
10544.5
0.09
991.82
11.11
Total
18567.2
0.41
7689.82
100
3D model for the structure.
3. Simulation Tool
The S. J. Carew building is modeled by employing IDA-ICE as a dynamic thermal simulation tool. This program is selected because it is widely accepted as a viable thermal building performance simulator towards the study of power usage and indoor thermal climate of whole buildings [19]. The IDA-ICE program uses symbolic equations framed in a modeling language and a variable time-step differential-algebraic (DAE solver). The models can be expressed through the Neutral Model Format (NMF)/Modelica and act as both computer code and readable document, which are applicable to various simulation environments [20, 21].
The simulation tool IDA-ICE 4.7 is employed to predict the power usage and interior climate of the S. J. Carew building. The IDA-ICE 4.7 tool is ideal for modeling of multiple-zone HVAC systems as in the S. J. Carew building. IDA-ICE 4.7 is able to determine the general thermal comfort level of the building by measuring the internal air quality (IAQ) and performing dynamic simulations. The heat exchanger uses controllers to maintain zonal temperatures, which can be set as fixed points by modulating control valves. Meanwhile, in the real system (as shown in Figure 2), a hot water valve collects relevant data for water heated by a heating coil. Although the building’s system features a single valve for its hot water production, the IDA-ICE software divides the single valve into four distinct units, enabling each air-handling unit (AHU) to have an individual valve [22, 23].
AHU2 for S. J. Carew building.
4. Building Model
IFC files were used to develop a simulation model from the building information model (BIM). As shown in Figure 3, the third floor of the building features an AHU and all floors in the building have their own individual AHU.
The third floor of the building.
Data from the Facilities Management Department at Memorial University are used to source construction information regarding building dimensions and elements such as windows, doors, and walls. The data then are inputted to the IDA-ICE program.
Figure 4 depicts eight different kinds of windows used in the Carew building, while Figure 5 illustrates the building’s south elevation. Heating system information (e.g., radiator type and position) is presented in Figure 6, and Figure 7 depicts the main room’s ventilation system for AHU1, AHU2, and AHU3. Each of the areas has individual internals loads (i.e., occupancy and light based on floor type and usage) that have been determined by applying national building code monthly values. The ventilation design determines the supply and exhaust air flows, with standard commercial building pressure coefficients applied. There are doors dividing nearly all of the areas in the structure, enabling bidirectional air flow (even through closed doors). Air tightness is measured as an n50-parameter, while infiltration and exfiltration are simulated using the IDA-ICE air flow network feature.
Details of windows used in the Carew building.
South elevation of the building.
Fourth-floor layout.
Main room’s ventilation system.
5. Simulation Results
The IDA-ICE 4.7 program was used to analyze energy use in the S. J. Carew building at Memorial University, St. John’s, Newfoundland. The analysis involved a number of factors, such as weather data, infiltration, external/internal heat gain, and overall heat capacity. The simulation was done for the course of one full year (January 1, 2016, to December 31, 2016). The space heating and total energy consumption analysis results for the building in five points as following.
5.1. AHU Results
(i) AHUs Temperature. Figure 8 shows AHU1 supply air and return air, as well as the outdoor air temperature. As can be seen, the air temperature represents a mixture of temperatures from individual zones, while the air supply temperature represents the temperature of the air terminal zone following any alterations made to the duct or fan systems. The supply air temperature set-point refers to the air temperature prior to these alterations.
AHU1 supply air, return air, and outdoor air temperature.
(ii) AHU Airflow. The flows represent the total flow from every zone impacted by AHU1 (Figure 9) and have been multiplied according to weight (i.e., how many zones are the same type). This value is included in every zone. The flow is volumetric and assessed for actual temperatures. Therefore, it can differ from the set-point flow for each zone. Set-point flows are determined from mass flows based on variations in density.
AHU1 supply and return air flow.
(iii) Heating and Cooling AHUs Coils. Figure 10 depicts the central cooling and energy for AHU1. In circuits that are water-based, the energy consumption can be measured in the circuit according to temperature changes and mass flow (i.e., supply/return). Therefore, heat can be calculated after generation losses are included, but prior to the calculation of emission and distribution losses. Cooling energy is included as a positive quantity in the calculations.
Central heating and cooling power of AHU1.
5.2. Heat Balance
Figure 11 depicts the zones’ sensible (dry) and full latent (moist) heat balance. To find the sensible heat balance only, the details for the zone’s power report need to be logged. In this setup, the control volume indicates air-wetted surface area located at the room unit zone-side, backed by an air gap. Heat balance contributions can be allocated as shown in Figure 11.
Results for building heat balance.
(i) Equipment Heat. This type of heat emanates out of equipment like printers and computers as a form of radiant or convective heat.
(ii) Floor and Wall Heat. As the control volume is positioned directly below the surface of floors and walls, any measurement of heat indicates the presence of conductive heat passing through the structural element. This type of thermal energy can include net transmission and heat storage, in addition to internal heat functions (e.g., in-floor radiant heat). The thermal energy that has been stored as a function of room masses, such as in furniture, also is in this category.
(iii) Daylight Heat. This type of heat describes sunlight streaming into open doors or through windows, taking into consideration any short-wave radiation which exists. Solar radiation that has been absorbed and retransmitted is excluded from this category.
(iv) Heating/Cooling Room Heat. This type of heat is represented by controlled room units (e.g., radiators or chilled beams). In hydronic systems, there is an automatic calculation to account for the radiation and convection aspects, as described in the manual. Floor heating is excluded from this category.
(v) Window Heat. This type of heat describes heat emitted from window surfaces, such as through retransmitted absorbed solar radiation or through conduction. Long-wave radiation entering via openings such as opened doors falls under this category of heat. Solar radiation can have two major impacts on a room’s heat:
It can be absorbed by the window covering or pane and then emanate through the room as a radiative or convective process.
It can be directly transmitted as short-wave radiation and be reflected by room surfaces until finally absorbed by internal room masses (furniture, equipment, people, etc.).
(vi) Air Flow Heat. This describes all air flows, including infiltration, flowing from other zones and mechanical ventilation.
5.3. Energy Delivered
The report on the energy delivered provides a general overview of the total energy purchased or generated in the S. J. Carew building, as shown in Figure 12. The reported items are matched to the energy meters. The report also shows the primary form of energy employed, as well as the cost and estimates for CO2 emitted. These are presented according to the structure’s floor area and with regard to absolute values. Conversion factors from the meter energy to other measures are given as energy meters.
Energy delivered.
5.4. Results: Energy from Systems
The results also provide a general review for HVAC system energy flow. As shown in Table 2, the review is in three categories: use energy, utilized free energy, and auxiliary energy. The results provide a means to validate real data and then apply this data for system identification to obtain the Carew building’s HVAC system’s state-space model.
Energy use distribution.
Used energy
Utilized free energy
Auxiliary energy
Months
Zone heating (kWh)
Zone cooling (kWh)
AHU heating (kWh)
AHU cooling (kWh)
AHU heat recovery (kWh)
AHU cold recovery (kWh)
Fans (kWh)
Pumps (kWh)
1
244041
108.9
104998
0
639269
0
54672
119.2
2
250651
100.8
108788
0
643065
0
51008
137
3
133363
844.5
99535
0
620252
0
54752
91.4
4
122152
1679
69885
0
527862
0
53215
53.9
5
110590
13874
7648
12537
404205
0
55451
31.4
6
9588.9
37040
5861.1
112465
215982
-91.8
54228
180.1
7
791.4
59348
56780
189447
101874
-196.3
56426
402.5
8
482.2
44485
45673
205557
94640
-0.2
56442
187.9
9
6229.5
27183
75142.4
66266
198208
-0.1
54267
47.6
10
90232
7316
113167
845.2
369420
0
55551
14.4
11
113625
689.4
122842
0
502646
0
53290
42.4
12
240978
40
129368
0
617644
0
54769
100.8
5.5. Results: Energy from Zones
The results provide information on the sensible heat balance in the Carew building’s zones. Information on total (i.e., dry and wet) heat balance, are presented in the heat balance Figure 6. The data are provided for monthly and one-year (simulation period) basis. Figure 13 provides details on envelope transmission losses, with control volume being positioned at the surface of the floor as well as on the ceiling and inside walls. For slab (embedded) cooling and heating processes, the control volume also involves activated layers and thus includes large thermal masses.
Details on envelope transmission losses.
6. Simulation Validation for IDA-ICE
A viable model must provide accurate results and also meet the required specifications. In the present work, the building data used was provided by Memorial University’s Department of Facilities Management and the Honeywell Office. The data provided in Tables 3 and 4 (energy and hot water consumption) for the S. J. Carew building were collected between April 2012 and May 2017, inclusive. These were used to compare, and contrast power consumption derived from real data with power consumption derived from the IDA-ICE program data.
Energy from hot water consumption.
Y2012
Y2013
Y2014
Y2015
Y2016
Y2017
MONTH
MMBTU
MMBTU
MMBTU
MMBTU
MMBTU
MMBTU
APR
1,646
1,838
1,932
2,333
2,109
1,640
MAY
1,318
1,446
1,806
1,956
1,745
1,985
JUN
1,164
1,251
1,273
1,615
1,391
0
JUL
774
1,067
856
1,464
1,150
0
AUG
525
932
974
1,019
1,071
0
SEP
391
958
1,017
1,289
1,368
0
OCT
1,438
1,519
1,510
1,793
1,906
0
NOV
1,659
1,740
1,803
2,117
1,946
0
DEC
1,915
2,313
2,143
2,073
2,732
0
JAN
2,389
2,228
2,638
2,546
2,898
0
FEB
2,051
2,175
2,512
2,276
2,473
0
MAR
1,939
2,405
2,615
2,502
2,605
0
TOT
17,210
19,871
21,079
22,983
23,394
3,625
LITRES
565,029
633,115
666,618
719,885
762,559
116,177
Electrical power consumption.
Y2012
Y2013
Y2014
Y2015
Y2016
Y2017
MONTH
kWh
kWh
kWh
kWh
kWh
kWh
APR
366,300
375,458
462,454
526,557
485,348
444,139
MAY
393,773
425,824
471,612
503,663
398,352
434,982
JUN
434,982
357,143
407,509
503,663
489,927
0
JUL
407,509
508,242
563,187
425,824
467,033
0
AUG
366,300
366,300
471,612
462,454
494,505
0
SEP
402,930
393,773
425,824
430,403
370,879
0
OCT
439,560
407,509
425,824
526,557
434,982
0
NOV
412,088
476,190
512,821
439,560
489,927
0
DEC
393,773
444,139
526,557
434,982
508,242
0
JAN
439,560
462,454
476,190
448,718
476,190
0
FEB
398,352
439,560
494,505
444,139
430,403
0
MAR
402,930
416,667
467,033
508,242
476,190
0
TOT
4,858,059
5,073,260
5,705,128
5,654,762
5,521,978
879,121
The first step for comparison was to verify design details for the Carew building. These details included aspects such as building materials, location, dimensions, total area, etc. The second step was to make a comparison using the file for outdoor air temperature/weather as represented in the IDA-ICE program (based on readings from St. John’s Airport, [ASHRAE, 2013]) and the building’s actual outdoor air temperature obtained from the Honeywell software data. Figure 3 depicts the sensor (-1.4° OA); the average of the temperature readings from 2016 in one-hour time samples for both data was the same.
A viable model needs to have both accurate results and the ability to satisfy any required specifications. Figure 14 shows the IDA-ICE model of hot water usage from January to December 2016. The energy consumption for hot water was more than 800,000 kWh in Jan and Dec. Also, it was almost 300,000 in the summer time (Jul and Aug). Furthermore, although the actual data for hot water usage measured slightly low in some months and slightly high in others, it compared well to the IDA simulated data. Regarding overall energy consumption, the modeled data are only somewhat different than the actual data. Figure 15 shows the actual (measured) data as moderately higher than the simulation data, but these slight differences could be due to discrepancies in the lab readings due to miscalibrated equipment.
Energy use from hot water.
Electrical power use.
7. System Identification
Our study used the IDA-ICE 4.7 simulation software for measuring the interior environment as well as the overall energy performance. This software is able to simulate and model multiple-zoned HVAC systems and is also gauge interior air quality (IAQ), energy requirements, and thermal comfort levels. To model the S. J. Carew building, the zonal inputs and outputs must be included in the identification data. There are three main steps in system identification [24–27]:
(i) Collecting the data towards model identification.
(ii) Choosing an appropriate model structure.
(iii) Building a model that provides the best functionality (i.e., satisfies specifications and gives accurate results).
During these steps, the focus is on optimizing the chosen model to suit a real-life system. In this study, a structure is used that features four AHUs as a means for identifying the state-space model of the system. The data used for system identification were collected during the winter months, which means that the cold-water valve was not operating. Additionally, because the S. J. Carew building has four floors, the system features twelve inputs and twelve outputs overall, calculated from three inputs (U) and three outputs (Y) per floor. Figure 16 illustrates the details.
Zonal Temperature °C (Tz) (y1,y4,y7,y10). These data are derived from the IDA-ICE software. Although the actual system features sensors in every room, the temperature on each floor still needs to be measured. The data from the IDA-ICE software are used to control the hot water valve. Figure 17 illustrates the outputs.
Hot Water Valve for Heating Coil/Zones u1,u4,u7,u10. These data are also derived from the IDA-ICE software. In the actual system (as shown in Figure 16), a hot water valve collects data on hot water use. Note that this system only has one valve for hot water production, whereas in the IDA-ICE software there are four valves, which enables every floor to have a separate valve. Figure 18 shows these inputs as percentage of opening and closing operation of the hot water valves.
Fresh Air Dampers u3,u6,u9,u12. As shown in Figure 18, the fresh air sensors positioned in AHUs are able to gauge, in percentage, the amount of fresh air entering the building. The sample time Ts used in these calculations is used in all data.
CO2 Levels (CO2) (y3,y6,y9,y12). This data is obtained from the sensors for return air flow ducts for individual AHUs. Figure 19 depicts CO2 levels occurring in AHUs. These outputs can be applied in moderating fresh air dampers.
Static Air Pressure (Ps) (y2,y5,y8,y11). This data comes from two sensors: one for hot ducts and one for cold ducts. As illustrated in Figure 20, these outputs can be applied to the control of supply fan speed.
Supply Fan Speed u2,u5,u8,u11. This data is derived from AHUs, the sensor. Figure 18 depicts the sensors measuring the fan speed of AHUs. The sample time Ts for gathering the data is one hour, and the input signals are obtained in percentages.
Inputs and outputs of the system.
Zones temperature.
All inputs of the system.
A CO2 level of AHUs.
The static air pressure of AHUs.
System disturbances (w) can occur with changes to wind speed/direction and outside temperature. These changes are recorded in the IDA-ICE software and the input/output signal data series organized through MATLAB. The ordering of the data is imperative before moving onto the next stage, which is system identification using the System Identification (SI) Toolbox. Every individual data set is cut in two: one-half represents estimation data, while the other half represents validation data.
Data from the first half (Data 1) are organized as input/output sets using MATLAB. The inputs are arranged into 12 columns, with every column relating to a specific input signal. Note that the number of rows is equal to the number of simulation period hours. Similarly, the outputs are also arranged into 12 columns, with every column relating to a specific output signal and the number of rows being equal to the number of input arrays. In these calculations, both the estimation and validation data have a 90-day time frame. Figures 21, 22, and 23 show some of the time plots of this data for estimation. The CO2 level of AHU3, static air pressure of AHU2 and zone temperature of AHU1 is shown. The output is the upper plot and the input is the lower plot.
Input and output signals of CO2 level of AHU3.
Input and output signals of the static air pressure of AHU2.
Input and output signals of zone temperature of AHU1.
Every state-space model is estimated in the SI Toolbox. The models undergo a comparison based on the degree of accuracy between the validation data sets’ estimated and measured (i.e., real) outputs. In comparison, the estimated and real outputs are plotted for every model, after which a numerical value is allotted regarding the model’s “fit.” Using the SI Tool, the estimated outputs for numerous models are able to be plotted quickly and at the same time, with the model showing the highest value (i.e., the best “fit”) deemed to have the greatest reliability. As an outcome of the process, we can obtain state-space models for Data 1 data groups and compare models across different seasons. For detrending the data, there are no alterations made to any relative differences among inputs and outputs.
Figure 24 shows the output of the model that follows the temperature of a zone in AHU1 with the same output as the real system. The agreement between these graphs can be seen as a percentage of the error. Ideally, this result is 57%. Also, Figure 25 shows that the system performance percentage for the estimated model and the actual system of CO2 level of zone 3 in AHU3 was 75%. Also, the percentage of the estimation model of the static air pressure the simulated or predicted model output is shown together with the measured validation data in Figure 26.
Validation of real measurements and outputs of CO2 level model fit 75%.
Validation of real measurements and outputs of zone temperature model fit 57%.
Validation of real measurements and outputs of the static air pressure model fit 49%.
The part of the system data that the model could not describe is called the residuals [28–30]. They contain important information about the quality of the estimated model. The cross-correlation between residuals and the correct model does not exceed the confidence level [28]. If this is the case, the original model has captured the underlying properties of the system. The remaining autocorrelation indicates whether the error pattern is accurate. Standard process models do not evaluate the error model and unknown interference is not in the original model, thus the remaining runtime is not used for model verification. Figures 27, 28, and 29 show plots of the autocorrelation and cross-correlation of system responses to inputs. It is clear from the cross-correlation diagram of these figures that the estimated model is very similar to the responses of the system to the inputs; the correlation curves lie between the dashed lines.
Autocorrelation y1 and cross-correlation of system responses u1 & y1.
Autocorrelation y5 and cross-correlation of system responses u5 & y5.
Autocorrelation y12 and cross-correlation of system responses u12 & y12.
To determine model settings for the system, a linear parametric model can be estimated from a state-space structure. In general, the state-space model discrete-time settings generally feature the following structure:(1)Xt+Ts=Axt+But+KetYt=Cxt+Dut+etwhere the xt represents the states of the system and yt,utande(t) represent the output, input and error.
The A, B, C, DandK matrices contain the model parameters, and Ts is the sampling time of the system.
For modeling multi-input/multi-output (MIMO) systems, state-space models have been proven to be the most popular option, likely because the state-space method is relatively simple. For the system used in the present work (twelve inputs and twelve outputs), the discrete-time state-space model for order 12 (sampled as Ts = 3600 s) and the A, C, and K matrices are as follows: (2)A=1.0185960.024987-0.09841-0.018090.0298260.0115510.122706-0.09652-0.02257-0.033340.0211340.044136-0.048020.953931-0.18782-0.109130.0179030.0127050.022660.044524-0.00967-0.02612-0.007550.0130650.0067350.1614721.001750.153348-0.17807-0.0057-0.074450.089848-0.18638-0.04874-0.0658-0.090620.1025720.1061590.0150120.876272-0.056670.0400260.0348990.0638950.046081-0.03180.0580220.151734-0.010120.1193880.0320580.1938050.8757730.1050070.06194-0.11106-0.276650.118682-0.04009-0.123570.0341820.0019810.218893-0.064990.0036830.698464-0.156370.16687-0.063750.068068-0.13830.073127-0.06088-0.037180.0182760.0266440.006168-0.004840.7000560.1346080.025221-0.058180.065294-0.061740.128313-0.04206-0.25685-0.114380.260537-0.275070.0008510.6324010.220702-0.200810.2536060.220525-0.014090.0718270.152001-0.061850.0795890.4492750.2203420.1269010.6791720.0924660.1459620.1147910.029159-0.04508-0.122380.0735430.029797-0.28537-0.131110.07257-0.047010.726139-0.029710.1366160.147829-0.16104-0.06659-0.355560.222772-0.21727-0.0814-0.011910.0911890.1175580.9004310.2623430.0183220.1249130.198312-0.061540.143552-0.172690.228241-0.131410.238771-0.21135-0.018980.833612C=-32.38193.0325220.8503611.3204611.6540510.228377-1.254750.711672-1.122640.011717-0.989780.435679-4.962814.8258742.8347661.0769082.428522-0.775371.6508160.8299841.838859-1.4162-1.793840.42955-789.668315.0027-248.299258.1054147.2952-161.542-225.011100.9488120.897383.13542-194.598-118.225-32.5823.7390061.1892150.4568382.628477-0.976740.321862-0.63062-1.013030.477141-0.117920.126173-7.894853.6020450.5342730.073331-0.296110.2605422.0931790.5935470.695332-0.010572.583058-3.38128-596.62338.8736128.9263-51.9804-268.546-222.936-531.204-263.77958.86106194.4376142.4278-27.9924-32.52223.3754690.211781-0.906222.133119-0.3489-0.90599-0.03896-1.003580.530528-0.06908-0.21883-5.45871.7588071.5881130.6614211.924165-2.585041.5629192.2868611.189903-0.49241-0.055580.912613-928.097287.0009-266.691428.3191393.2388-234.874-209.79590.12368-223.088-365.80631.59467-280.605-32.4433.893139-0.188643.3452913.703252-0.71339-0.366980.158930.455992-0.003490.5980771.910353-6.37360.3621661.906773.0732532.246129-1.032950.390075-1.420841.8386770.77659-1.909331.183602-736.847466.1483-328.71973.65842118.4599-142.761-464.838183.0777-43.4737-86.4524-281.692-65.8117K=-0.022470.002415-0.000170.006293-0.001730-0.00404-0.00078-0.000120.0063620.003690.000173-0.037660.0088260.0002010.0194320.01211300.006975-0.0074300.012768-0.015390.0002420.0335930.01661100.0047040.0028910.000146-0.01335-0.00477-0.00877-0.025820.004979-0.000430.03453-0.003960.000281-0.016390.0029480.000201-0.059720.0036920.0001460.0404230.0015160.003516-0.038920.023262-0.000180.037843-0.00404-0.000280.00469-0.013570-0.017510.01406600.0518350.000717-0.00015-0.03922-0.00076-0.00837-0.00337-0.011720-0.008870.006588-0.00045-0.025820.0052040.0004550.0304620.009431-0.000260.0119970.0027810.0001450.00114-0.02152-0.0003-0.01015-0.013340-0.049760.014821-0.000840.0019320.01672700.036388-0.024710-0.014450.0181990.000477-0.027530.0067310.0001560.001845-0.00527-0.000220.0307580.0105230.000299-0.03034-0.009410.0004460.0106680.00654800.0215020.01097-0.00015-0.01222-0.0214-0.00045-0.01857-0.026410-0.045460.014413-0.000110.011490.0138940.0001320.034855-0.01938-0.000530.012010.001892-0.000490.008802-0.024050-0.066580.00417-0.000370.0673780.021370.000465
8. Conclusions
In this paper, the S. J. Carew building was modeled using the IDA-ICE program using all details of HVAC system and instructions of the building. This model provides good approximations comparing the consumption of hot water and electricity with the measured data for a full year (2016). It also compares the average of the outside temperature of the weather file of IDA-ICE program and the measured data. All system inputs and outputs were selected, and a linear state-space model was identified describing the thermal response of the system. The dynamic model is required to design and test system controllers before actual implementation. The model was derived using MATLAB’s System Identification (SI) Toolbox. The model has twelve state variables, twelve inputs, and twelve outputs. The model responses when compared with actual data are within the allowed range. Validation data and autocorrelation function for the residuals as well as the cross-correlation function between input and residuals are computed and presented.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work is funded by the Ministry of Higher Education of Libyan Government, which is managed by CBIE in Canada. Also, the authors thank the Honeywell Office and Department of Facilities Management at Memorial University for providing them with structure details of the building and access to Honeywell software of the HVAC system.
AbelE.Per-Erik NilssonL. E.Per FahlénL. J.Achieving the desired indoor climate-energy efficiency aspects of system design2003AwbiH.2003Spon PressGoodfellowH.2001ElsevierAmerican Society of Heating Refrigerating and Air-Conditioning Engineers2009Abdo-AllahA.IqbalT.PopeK.Modeling, analysis, and state feedback control design of a multizone HVAC system2018201811430358010.1155/2018/4303580YiuJ. C.-M.WangS.Multiple ARMAX modeling scheme for forecasting air conditioning system performance2007488227622852-s2.0-3425088466410.1016/j.enconman.2007.03.018MustafarajG.ChenJ.LowryG.Development of room temperature and relative humidity linear parametric models for an open office using BMS data20104233483562-s2.0-7264909090010.1016/j.enbuild.2009.10.001MustafarajG.LowryG.ChenJ.Prediction of room temperature and relative humidity by autoregressive linear and nonlinear neural network models for an open office2011436145214602-s2.0-7995374928710.1016/j.enbuild.2011.02.007RablA.Parameter estimation in buildings: Methods for dynamic analysis of measured energy use1988110152662-s2.0-002396036410.1115/1.3268237SondereggerR. C.1978NJ, USAPrinceton Univ.BoyerH.ChabriatJ. P.Grondin-PerezB.TourrandC.BrauJ.Thermal building simulation and computer generation of nodal models199631320721410.1016/0360-1323(96)00001-72-s2.0-0030143790AndersenK. K.MadsenH.HansenL. H.Modelling the heat dynamics of a building using stochastic differential equations200031113242-s2.0-003389063210.1016/S0378-7788(98)00069-3JiménezM. J.MadsenH.BloemJ. J.DammannB.Estimation of non-linear continuous time models for the heat exchange dynamics of building integrated photovoltaic modules20084021571672-s2.0-3634898672110.1016/j.enbuild.2007.02.026MadsenH.HolstJ.Estimation of continuous-time models for the heat dynamics of a building199522167792-s2.0-002927300210.1016/0378-7788(94)00904-XMeiL.InfieldD.EickerU.FuxV.Thermal modelling of a building with an integrated ventilated PV façade20033566056172-s2.0-003773889610.1016/S0378-7788(02)00168-8JeongY.HaghighatF.Modelling of a hybrid-ventilated building using ESP-r20161212713910.1080/14733315.2002.11683629JokisaloJ.KurnitskiJ.KorpiM.KalameesT.VinhaJ.Building leakage, infiltration, and energy performance analyses for Finnish detached houses20094423773872-s2.0-5574909166710.1016/j.buildenv.2008.03.014StadlerM.FirestoneR.CurtilD.MarnayC.On-site generation simulation with EnergyPlus for commercial buildings, ACEEE summer study on energy efficiency in buildings2006California, USAPacific GroveEQUA Simulation AB2010http://www.equa.se/en/ida-iceSahlinP.ErikssonL.GrozmanP.JohnssonH.ShapovalovA.VuolleM.Whole-building simulation with symbolic DAE equations and general purpose solvers20043989499582-s2.0-3374843450510.1016/j.buildenv.2004.01.019SahlinP.ErikssonL.GrozmanP.JohnssonH.ShapovalovA.VuolleM.Equation-based building simulation make it? Experiences from the introduction of IDA indoor climate and energy39Proceedings of the 8th International (IBPSA) International Building Performance Simulation Association ConferenceAugust 2003Eindhoven, The Netherlands94995810.1016/j.buildenv.2004.01.019Abdo-AllahA.IqbalT.PopeK.Modeling and analysis of an HVAC system for the S.J. carew building at memorial universityProceedings of the 30th IEEE Canadian Conference on Electrical and Computer Engineering, CCECE 2017May 2017Windsor, CanadaIEEE1410.1109/CCECE.2017.79466892-s2.0-85021814536Abdo-AllahA.IqbalT.PopeK.Modeling, analysis, and design of a fuzzy logic controller for an ahu in the s.j. carew building at memorial university2018201811454038710.1155/2018/4540387LjungL.2006USAPrentice-HallMR3444790RaillonL.GhiausC.Study of error propagation in the transformations of dynamic thermal models of buildings2017201715563614510.1155/2017/5636145MR3699151ChintalaR. H.RasmussenB. P.Automated multi-zone linear parametric black box modeling approach for building HVAC systemsProceedings of the ASME 2015 Dynamic Systems and Control Conference2015USAAmerican Society of Mechanical EngineersV002T29A00410.1115/DSCC2015-9933AfrozZ.ShafiullahG. M.UrmeeT.HigginsG.Modeling techniques used in building HVAC control systems: A review20188364842-s2.0-8503773213610.1016/j.rser.2017.10.044LjungL.ZhangQ.LindskogP.IouditskiA.SinghR.An integrated system identification toolbox for linear and non-linear models200639193193610.3182/20060329-3-AU-2901.00148ScottonF.HuangL.AhmadiS. A.WahlbergB.Physics-based modeling and identification for HVAC systems?Proceedings of the 12th European Control Conference (ECC '13)July 2013IEEE140414092-s2.0-84893323157BarkleyA.SantiE.Improved online identification of a DC-DC converter and its control loop gain using cross-correlation methods2009248202120312-s2.0-6944908682910.1109/TPEL.2009.2020588