Modeling , Analysis , and State Feedback Control Design of a Multizone HVAC System

A HVAC system is modeled by applying a state space MIMO (multi-input/multioutput) system method for control system design and analysis. Thermal models are developed using the simulation program IDA Indoor Climate and Energy. The building has four floors in total, with separate air-handling units (AHUs) on each floor. The system’s eight main input data are hot water and the energy usage for each AHU, while the eight main outputs are return airflow temperature and CO 2 levels for AHUs. The factors of wind direction and velocity are also applied as disturbances. By comparing usage data on simulated power consumption versus measured data for the three months of October, November, and December 2016, good agreement was achieved with simulated data. The main aim is to develop a state feedback controller and then apply it toward optimal functionality of a control system. After utilizing the MATLAB identification toolbox, a MIMO system-based state space model is developed.


Introduction
A main aim in developing an optimal HVAC (heating, ventilation, and air-conditioning) system is to create a comfortable environment for occupants with reduced energy inputs [1].However, heating and cooling loads typically change according to the exterior environment, as well as with the specific needs of the users.HVAC systems require a control system to keep the comfort level and air quality relatively constant with variable conditions.Furthermore, power usage can be greatly decreased if the system is suitably controlled.
A method that is based solely on measured data is one possible means for obtaining a mathematical interpretation of the system.The model can be used to determine system parameters in cases where input and output variables are already available.This modeling method is useful if the system is constructed and data related to performance can be readily obtained.Compared to forward models, those that are data-driven can identify system approaches that can prove to be both easier to use and better performance predictors.Modeling methods that primarily use data are classified as "system identification" (ASHRAE 2005).
A wide range of research over the past few decades supports the suitability of applying the system identification approach in energy simulation and in determining and analyzing the moisture, cooling, and heating environments in buildings.Applying data gathered from a building management system, Lowry and Lee (2004) studied the outcome of using a data-driven model to gauge thermal response [2].A few years earlier, Madsen and Holst (1995) utilized nearly the same system identification approach to determine a structure's heating dynamics obtained from data on discrete time performance [3].Cunningham (2001) applied system identification methods to find moisture release rates in a structure based on psychometric data [4], and Mechaqrane and Zouak (2004) incorporated system identification in their investigations on the prediction of interior air temperature in residential structures [5].A few studies have also focused on comparing test approach models with theoretical predictions using simulation software like TRNSYS.Peippo et al. (1991) determined the dynamics of a structure according to discrete time results obtained from simulation software using a onehour time step [6].Despite the approach's clear potential, there are numerous obstacles to the simulation of a structure's dynamic response when applying simulation software.Currently, simulation software can be significantly time-consuming and the results are often missing important information on the fast dynamic behavior of a structure (i.e., in the order of seconds), as the majority of available software utilizes a discrete time step that is typically set at one hour.Problems arise when, for instance, data is needed on fast dynamic behavior for a control strategy (e.g., on/off), but it cannot be obtained because it is situated within the time step.
This paper presents a simulation of an entire structure, applying the simulation program IDA Indoor Climate and Energy 4.7.The energy usage (including heating and cooling) of the S. J Carew Building of Memorial University of Newfoundland is investigated with a 3D numerical model using IDA ICE software.Established in 1998, the IDA Indoor Climate and Energy program can investigate separate thermal climate zones [7].There are four main objectives of the present work: (1) Use real data to create entire structure and then validate the results.
(2) Determine the viability of using system identification in decreasing the time required for calculating the simulation of complex structures.
(3) Determine the viability of using system identification for identifying the dynamics of climate control design in a structure when using discrete time data of onehour sample time.
(4) Use state feedback (classical) control by applying a state space model.  1 provides a description of the structure's amenities and Figure 1 illustrates the structure.Initially, a section of the building is investigated using the simulation software IDA ICE [7].Data is applied on the weather for St. John's and also construction details (e.g., windows, doors and walls, and radiators, fans, and pumps) to obtain energy consumption information for the simulation.In addition, we apply AutoCAD files drawings for determining the dimensions of the building such as the height of the building, size and positions of the windows and doors, and information on ventilation and heating systems.

IDA ICE Simulation
Validation.Determining whether the model satisfies the requirements and whether the results are accurate is critical to developing a feasible model.In this study, data is obtained from the Department of Facilities Management at Memorial University and compared with the structure's hot water and power usage.Figure 2 illustrates the power usage and Figure 3 illustrates the hot water for the entire structure of the S. J. Carew building from April 2012 to May 2017.Table 2 and Figure 4 illustrate the IDA ICE data for three months, October to December 2016, including energy consumption in the building, from lighting, HVAC, and equipment, such as computers, TV, cookers, and fridge-freezer.
The output variables of exterior temperature, hot water power usage, and power usage are compared with simulation data (three months' worth) from IDA ICE software, as follows.

Exterior Temperature.
Exterior temperature from the IDA ICE program and data weather files were compared with data from October to December 2016, using one-hour time samples, as shown in Figure 5.

Power Usage.
The modeled data differs slightly from measured data for power usage.As illustrated in Figure 6 the measured data slightly exceeds the modeled data, which is likely due to differences in laboratory equipment.This information is needed for the identification data and reference model for modeling the building.

System Identification
Three distinct stages can be defined when using system identification [8]: (i) Data gathering for identifying the model.
(ii) Choosing a suitable model structure.
(iii) Developing a model which offers optimal system functionality.The S. J. Carew building has four AHUs which are needed to identify the system's state space model.The system has eight inputs; 1, 2, 3, and 4 represent the hot water power usage (kWh) and 5, 6, 7, and 8 represent the power usage (kWh) inputs of the system.The outputs, 1, 2, 3, and 4 represent the level of CO 2 (PPM), while the 5, 6, 7, and 8 outputs represent the zone temperature ( ∘ C), as shown in Figure 8.

Input Signals.
A rise in airflow temperatures and radiator heat of the zonal temperature and hot water flow (minimum and maximum temperatures) is highly plausible.These form the zones' sole control variables; thus identification system signals can be used as input signals.The system's power usage (PU) and hot water power usage (WPU) as a time function, for each AHU, are illustrated in Figure 9.

Output Signals.
The system outputs are defined as the temperature and CO 2 level in the return air flow.Returned air CO 2 levels are shown in Figure 10 and returned air temperature variations are shown in Figure 11.

Choosing a Model Structure.
A model structure is selected from range of structures which are roughly categorized as being either linear or nonlinear.However, because we are using a nonlinear system in this work, we choose the ARMAX model.

Identifying the Model.
For the identification decision process used in preprocessing the data, the decision process can be categorized into five stages [9]: (a) Deciding the optimal model structure (e.g., ARX, ARMAX, and process models) for our purposes.
(b) Deciding the model order.
(c) Deciding the optimal estimation approach.
(d) Launching the identification process.
(e) Reviewing and validating the results.

State Space Model.
The total quantity of the system's independent components is assumed equal to the total number of the state variables, .The system's total power and state time derivative variables set the system's power change rate, and the system's state variables, at time, , give enough data for calculating the values of the system's variables for that time [10].The multizone HVAC system was developed as single system with 42 states and nine outputs [11,12].
Matrix  8 * 8 and matrix  8 * 8 form the system's properties.Hence, the output variables impact the output equation matrices ( 8 * 8 and  8 * 8 ).The matrices are calculated by applying the MATLAB tool box for system identification.Additionally, using MATLAB, we can formulate the system's dynamic behavior for an arbitrary input and lsim (sys, u, t, xo) function simulation.The system uses eight inputs, each depicted as , while the time samples are given as  vector and   indicates the system's starting values.The system  outputs response, in Figure 12, shows the zonal temperatures, while Figure 13 indicates the zonal level of CO 2 .
As illustrated in Figure 13, the system responses fall within the correct range if the starting transient is neglected.

Control Strategies
HVAC systems control interior environmental factors such as room humidity and temperature in commercial or residential structures, with the overall aim of giving users a comfortable working and living environment.HVAC systems are so widely used that they comprise at least 50% of global power usage [13][14][15][16].In addition to creating a comfortable working and living environment, HVAC controls typically need to also consider maximizing energy efficiency.System identification and state space model for one AHU of the S. J. Carew Building was presented without control [17].In earlier studies, researchers investigated mainly humidity and temperature levels in modeling HVAC systems [18][19][20][21][22].A nonlinear HVAC model was introduced [18,19] that involves a temperature/humidity ratio and an observer for determining approximate moisture and thermal loads.An adaptive fuzzy output feedback controller was developed [20] that can be premised on an HVAC system observer.Researchers proposed applying both a decentralized nonlinear adaptive controller [21].and a back-stepping controller [22] on a model.CO 2 concentrations are generally seen as having a notable impact on room comfort levels [23,24].Accordingly, some researchers are suggesting developing a hybrid HVAC system that can develop temperature continuous states, as well as CO 2 concentration as a discrete state [25,26].As these states can be highly interrelated, a more viable approach would be to integrate the presented discrete and continuous dynamics to form a model that assumes CO 2 concentration and temperature as states.
A system's possible future development can be predicted by the state of a dynamical system, the latter is essentially a group of variables.The control system is used to bring the nonlinear system into a stable state, while achieving the control targets.For one AHU unit system design feedback linearization technique has been applied [27]; the control of the cross-water heat exchanger [28] and feedback controller achieves global input-output linearization of greenhouse environments [29].In creating system dynamics by referring to state feedback, the multi-input system can be controlled as a linear state model.Hence, the feedback control can be devised through a step-by-step process that is based on putting closed loop eigenvalues in specific places.Figure 14 illustrates a control system that utilizes state feedback.The system which features linear process dynamics, processes disturbances, , reference input, , and controller elements, , is feedback gain and   is input gain.The feedback controller's main purpose is regulating the system output, , until it tracks the reference input during process dynamics uncertainty and disturbances.Performance specifications are critical to control design and the most important feature is stability.The aim is for the system's equilibrium point to stabilize asymptotically, but increasingly complex specifications can include obtaining the preferred properties from the system's step or frequency responses.Such desired specifications include rise times, overshoot, and settling times of step responses.In optimizing functionality, the system's disturbance rejection properties can be analyzed to find the best way to handle disturbance inputs, , while maintaining output, , at the required value specifications.
The system is represented by the following linear differential matrix equation: It is assumed that  = 0 while neglecting the disturbance signal, .The aim is to give the output a value, .In this scenario, state vector components must first be calculated.Thus, as the state at time, , already has the required data for predicting the system's behavior, the following time invariant control law becomes both a state function and reference input: The equation can be formulated as follows, if the feedback is maintained linear: In this equation, the reference value, , has been deemed constant, which refers back to Figure 14.Here, the negative sign works as an indicator of negative feedback being standard procedure.The closed loop system which we created after the feedback (2) has been applied to the system (4) can be calculated as To find feedback gain, , to set the characteristic polynomial of the closed loop system: This formulation is known as the "pole placement" or eigenvalue assignment problem.In the formulation,   input gain has no effect on system stability (the latter is instead created by − eigenvalues) but does have an impact on the steady state solution.Thus, we can set the equilibrium point and steady state output as In this formulation,   is the best choice, which then gives   = , which is the required value.Furthermore, because   is scalar: The variable   represents the opposite of the closed looped system's zero-frequency gain.Therefore, by applying   input gain and  which is (8 × 8) feedback gain matrix, the dynamics of a closed loop system can be modified until the required specifications are achieved.In this paper, the state feedback gain matrix is

Open Loop System.
The system has eight inputs and eight outputs.In this part, the system is presented by 64 transfer functions each relating one input to one output.(i) Input gain, addition of an input gain to the state control law.
(ii) Integral, including an overall follow-up measure.

Input Gain.
For the input gain   , presented in (8), the input gain controller is introduced to eliminate the stationary error associated with the complete status feedback controller for each constant input.The control law is designed such that the output () follows the reference input ().The controller maintains tracking in steady state only if the reference inputs are stage inputs.The gain, , is outside the feedback loop and makes the overall system sensitive to noise and disturbances.However, this is not "robust" since any change to the system parameter causes a nonzero error.

Integral Action.
The technique of integral control is another type of technique of placing the poles.It is also known as a tracking controller when it needs an output to track the input control signal.The output feedback is transferred to the controlled system via the integrator.The integrator, also called integral action, is used to increase the system type and reduce finite errors to become zero [30].
The configuration of the integrated control technique is shown in Figure 16.The introduction of the input integrator causes the controller to have a pole at  = 0, which helps to eliminate the constant reference and improve the robustness of the system.However, Figure 16 shows the block diagram of a system with state feedback and integral control using MATLAB Simulink.Simulations are performed for a controller structure where a unit's step input is [465 459 453 471 23.9 22.9 21.9 24.9] and signals are used as the reference signal.To accomplish one of the design requirements, the output signals should follow the given reference signals.Through simulation, mathematical modeling for the system is verified and the performances for the controller structures are analyzed.Also, the initial state,  0 , of the system for concentrations and indoor temperatures is taken from measured data.
Simulation results of the system with a measured initial condition of CO 2 level are  0 = [446.4440.6 435.44 453.4], and change set points to investigate the system's responses with state feedback control and integral action.The responses of the CO 2 level for the zones are illustrated in Figure 17.
Figure 18 shows the responses of the zones temperature with measured initial condition of the zone temperatures,  0 = [21.21 20.35 19.39 22.13].
For full state feedback controller, , with input gain,   , and integral action, the steady state error is zero.Figure 19 illustrates the controller action (steady state error) of CO 2 level responses (CO 2 -Er-Zones).
Figure 20 shows controller action (steady state error) of zones temperature responses (T-Er-Zones).

Conclusion
In this paper, the HVAC system of the S. J. Carew building was modeled using the IDE ICE program.This model provides    In addition, the outside temperature for the program and the measured data are compared for three months as the first part of the process.In the second part, the system identification toolbox was used to obtain the state space model of the multi-input and multioutput system MIMO.The model has eight status variables, eight inputs, and eight outputs, and model responses are within the permissible range.In the third part, a novel HVAC system model was developed that considers temperature and CO 2 concentration as the quantitative indices of comfort in a building.In applying an input-output feedback linearization method to linearize the HVAC system, one type of linear controller, pole placement controller with input gain, and integral action were able to regulate the linearized HVAC system at the desired set point without steady state error.Simulation results validated the proposed HVAC model, demonstrating its effectiveness in maintaining comfortable conditions.

Figure 1 :
Figure 1: 3D model of the S. J. Carew building using IDA Indoor Climate and Energy program.

Figure 4 :
Figure 4: IDA ICE data for three months (October to December 2016).

Figure 5 :
Figure 5: Difference between outdoor temperatures of IDA ICE and actual data.

Figure 7 :
Figure 7: Difference between hot water power usage of IDA ICE and actual data.

Figure 8 :
Figure 8: System with eight inputs and eight outputs.

Figure 9 :Figure 10 :
Figure 9: Actual data of power usage (PU) and water power usage (WPU) of the building, for each AHU.

Figure 11 :Figure 12 :
Figure 11: Actual data of return air temperature of the building for each AHU.

Figure 14 :
Figure 14: Block diagram of a system with state feedback () and input   controllers.

Figure 15 :Figure 16 :
Figure 15: Steps response of the open loop system with sampling time of 3600 s.

Figure 18 : 2 - 2 Figure 19 :
Figure 18: The outputs response of the zones temperature T-Zones with steps references T-Refs.

Figure 20 :
Figure 20: Controller action (steady state error) of temperature responses for each zone T-Er-Zones.

Table 1 :
Details of the building.
Water Power Usage.As illustrated in Figure7, in October to December 2016 IDA ICE model, hot water power usage measured 6,583 MMBTU.Hot water usage in October for actual data was lower than IDA data, but in December was opposite and November was almost the same.
2.3.Simulation Model.The IDA Indoor Climate and Energy 4.7 (IDA ICE) simulation tool is used to make an assessment on both the energy performance and interior climate.The tool can model multiple-zoned structures housing HVAC and can also be used for assessing thermal comfort, interior air Figure 6: Difference between power usage of IDA ICE and actual data.

Table 2 :
IDA ICE data for three months (October to December 2016).
Figure 17: Outputs response of the CO 2 level for the zones CO 2 -Zones with steps references CO 2 -Refs.
good approximation results in which hot water consumption and electricity consumption are compared with actual data.