Mathematical modeling describing the performance of open loop multivariable vapor compression chiller model: implications for sustainable HVAC in social housing development

Abstract

Introduction:

This study introduces a simplified multi-variable state-space model of a vapour compression chiller system, focusing on key components such as the condenser and evaporator.

Methods:

The model is based on a previously validated state-space framework and has been further developed to generate a transfer function matrix with two inputs (refrigerant and cold carrier flow rates) and two outputs (cooling capacity and outlet temperature). The simplified model offers an accurate yet computationally efficient representation of the system’s dynamics. Simulations of the open-loop system were carried out using SIMULINK/MATLAB software. The results were plotted to validate the model's ability to predict system responses accurately. These simulations encompassed various conditions, including operational changes, fouling on the evaporator, and thermal load fluctuations.

Results:

Open-loop simulation results show that the system reaches steady-state values with a response time of approximately 40 seconds for the cold carrier temperature and 70 seconds for cooling capacity. Meanwhile, the system exhibits strong coupling between inputs and outputs, with disturbances such as inlet temperature variations causing oscillations in the cooling capacity.

Discussion:

The novelty of this work lies in reducing a complex validated chiller model into a practical two-input, two-output framework, which enables simplified analysis and serves as a robust foundation for advanced multivariable control design. This contribution addresses the critical need for efficient HVAC control strategies in sustainable building applications. This research presents a foundational mathematical model for developing advanced feedback control strategies that can enhance system stability, improve energy efficiency, and support sustainable HVAC operations while maintaining thermal comfort. These advancements are especially relevant for social housing, where energy-efficient HVAC systems are essential for achieving sustainable and resilient building practices in extreme climate conditions.

1 Introduction and related scholarly work

Mathematical models are employed across various fields, including engineering, ecology, agriculture, medicine, and economics, primarily to forecast and regulate the operation of systems in real-world applications (Homod, 2013). Mathematical modelling is, for example, extensively used in construction and civil engineering (Addis, 2020), because replicating large-scale structures like buildings, bridges, or dams in a laboratory is impractical and costly, while modelling provides accurate insights into structural behaviour under various conditions. In the field of building energy management, creating precise and efficient mathematical models of a building’s components is a vital process (Pan et al., 2023). These components include the building envelope, orientation, and building services. Such models are crucial for predicting and optimising performance, especially in energy-intensive services like lighting and Heating, Ventilation, and Air Conditioning (HVAC) systems (Jang et al., 2019). Accurately predicting building service behaviour is crucial for energy-efficient operations. A reliable mathematical model supports optimal decision-making, helping to optimize energy use while ensuring performance.

HVAC systems are responsible for a significant portion of building energy use, often consuming up to 50% of a building’s total energy (Simpeh et al., 2022). They contribute substantially to greenhouse gas emissions, mainly when they operate on fossil fuels, which can harm the environment and exacerbate climate change. According to the International Energy Agency (IEA), energy consumption for air conditioning is projected to increase 4.5 times by 2050 in countries outside the Organisation for Economic Co-operation and Development (OECD), compared to a 1.3 times increase within OECD countries (Goetzler et al., 2016). This increase in demand is driven by higher incomes and better access to air conditioning, especially in developing countries. As a result, the need for HVAC system installations is likely to grow further in the coming years, making these systems even more vital for meeting future energy needs. Consequently, it is imperative to explore methods to decrease energy usage and enhance the efficiency of HVAC systems, thereby diminishing their carbon footprint and contributing to environmental conservation (Kathiravel et al., 2024).

To reduce energy use in buildings, indirect approaches such as passive design strategies are often recommended (Palmero-marrero & Oliveira, 2010). These strategies include features like improved insulation, shading devices, and green roofs, which help minimise energy demand for heating and cooling (Castleton et al., 2010). Additionally, cooling methods, including enhanced insulation, green roofs, and photovoltaic systems, have been explored to improve energy efficiency (Touqan et al., 2017). Retrofitting buildings to improve energy efficiency and sustainability is another important focus in building design (Ma et al., 2012). Despite these efforts to reduce energy consumption, HVAC systems remain essential in maintaining high levels of thermal comfort and air quality (Asim et al., 2022), which cannot be fully achieved through passive cooling strategies alone.

HVAC systems are inherently complex due to their multicomponent structure and strong interdependencies among elements such as evaporators, compressor and condensers, making their modelling and control a challenging task (Seyam, 2018; Touqan and Ameer, 2024). Inadequate control design can lead to increased energy consumption, frequent maintenance, poor performance, and user dissatisfaction (Angel, 2011). Enhancing system efficiency requires the integration of multidisciplinary knowledge through Systems Engineering, which applies structured methodologies to optimise performance (Subbaram Naidu and Rieger, 2011). Central to this process is the development of accurate and dynamic mathematical models that reliably represent system behaviour, forming the basis for effective control strategies (Li et al., 2010; Whalley and Abdul-Ameer, 2011). Such models enable the precise regulation of system outputs, thereby reducing energy waste and improving operational efficiency in real-time scenarios.

To contextualise the investigation of HVAC systems’ performance, they are vital for comfort in extreme-heat social housing, where high energy bills strain low-income residents. In the Arabian Gulf (>40°C), inefficient HVAC units worsen energy poverty and emissions; thus, despite passive design and retrofit efforts, dynamic modelling and control remain essential for maximizing efficiency without compromising comfort.

Hence, the motivation for this study is to derive a simplified and reliable mathematical model. While prior works validated a high-order state-space model (Yao et al., 2013), its complexity challenges real-time control applications. Computational efficiency is paramount for adaptive HVAC control in dynamic environments, yet most MIMO models either oversimplify coupling effects or lack experimental validation. A simplified two-input, two-output model that preserves essential dynamics can address this gap, facilitating deployment in resource-constrained contexts such as social housing.

This study aims to develop and validate a simplified two-input, two-output transfer function model of a vapour compression chiller, derived from a high-order state-space system. The model maintains key dynamic behaviour, improves computational efficiency, and supports real-time control. It facilitates advanced feedback strategies to boost stability, energy efficiency, and performance, especially in extreme climates. This model underpins cost-effective, resilient HVAC solutions aligned with energy and comfort goals in sustainable social housing.

Given the focus on energy-efficient HVAC operation and control-oriented modelling, the literature review situates the research within existing scholarly work. It outlines key modelling methods for HVAC, highlighting their advantages, limitations, and suitability for streamlined control models. This forms the basis for the methodological choices in this study.

Some researchers, especially those in (Homod, 2013; Afram and Janabi-Sharifi, 2014), have reviewed three distinct approaches to mathematical modelling. The first is the physics-based model, also known as the “white box model.” This model requires a deep understanding of the underlying processes. In the case of HVAC systems, physics-based models usually rely on equations related to mass and energy conservation, as well as principles of heat transfer. This approach necessitates comprehensive knowledge of the system being analysed. A significant benefit of physics-based modelling is its capacity to predict the system’s behaviour and estimate various outputs accurately. For instance, the research in (Ali, 2019) developed a mathematical load model for major household appliances, particularly air conditioners (ACs), based on physical and operational characteristics. Physics-based models, grounded in thermodynamic principles, offer high prediction accuracy and detailed system representation (Jorissen et al., 2021).

However, their complexity can limit real-time control applications (Silvestri et al., 2023; Zhao et al., 2021). These models are computationally demanding and require significant development effort (Y. Li et al., 2022). Moreover, while the model is built on fundamental physical relationships and validated against real data, it lacks a transparent methodology for parameter identification and uncertainty quantification, which are key requirements in white box modelling as emphasised by (Bacher and Madsen, 2011). Meanwhile, these models often encounter challenges in adapting to real-time control applications due to their inherent complexity and rigidity. Furthermore, they exhibit limited robustness when applied to systems with varying operating conditions, thereby reducing their practicality for adaptive HVAC control strategies necessary in dynamic environments such as social housing.

The second modelling approach is data-driven, often called a “black box model”. Studies such as those by (Ma et al., 2015; Zlatanović et al., 2011) have applied this method for HVAC systems, which are inherently complex with numerous parameters that make physics-based modelling difficult. Data-driven models can model complicated behaviour system (Zhou et al., 2017) and (Zhou et al., 2016). These models are beneficial for managing nonlinearities and are simpler to compute than traditional physical-based methods. However, the accuracy of these models can be significantly impacted by poor data quality, which is prevalent in large databases. Despite their advantages, data-driven methods face challenges such as poor interpretability and conservative predictions in extreme cases (Wu and Xue, 2024). In the same context, the study by (Zhang, 2023) provides a comprehensive critique of machine learning-driven building control systems, highlighting that although machine learning demonstrates potential in optimising energy consumption, considerable challenges persist, particularly concerning data quality, model interpretability, and the ability to generalise to new operational conditions. This highlights an important challenge in current data-driven methods, which often work as mysterious “black boxes” and can be difficult to validate in many real-world settings.

The third approach to modelling is the hybrid method, also known as the “grey box model,” which combines both data-driven and physics-based techniques. In this method, physics-based equations are used to develop the core model, while actual data helps to identify key parameter values. These data can be sourced from system specifications provided by manufacturers or from the system commissioning process. For example, in the study (Jin et al., 2011), a grey box model was created for cooling coil units, incorporating heat transfer equations and energy and mass balance principles. The model streamlined the calibration process by focusing on six parameters that represented the unit’s nonlinear characteristics. These parameters were derived from manufacturer data or experimental findings. The grey box approach faces challenges in parameter calibration and balancing physical insight with computational efficiency. Calibration is sensitive to data quality and quantity, with longer, diverse datasets improving model generalisation and accuracy (Murphy et al., 2021). The process requires substantial domain expertise to manage effectively. While grey-box models can provide physical interpretations of parameters, they may not fully capture parameter interactions or uncertainties (Pavlak et al., 2014). Despite these limitations, grey-box modelling can significantly reduce labour input compared to white-box approaches (Murphy et al., 2021). Overall, grey-box modelling remains the most effective modelling approach, offering a compromise between physical insight and data-driven efficiency. Still, it requires careful consideration of calibration strategies and data quality.

Traditionally, HVAC research has centred on single-input, single-output (SISO) models that isolate components such as compressors, cooling coils, air-handling units, or building envelopes (Li et al., 2014; Lapinskiene and Martinaitis, 2013). However, real HVAC systems are inherently multivariable, featuring several interacting inputs and outputs. Given the high energy consumption of these systems, it is more practical and representative to model them as multiple-input, multiple-output (MIMO) systems. This approach, where a change in one input affects multiple outputs, better captures the system’s complexity and energy demands, as discussed in (Hassan and Al-Samarraie, 2022), offering significant potential to optimize performance by manipulating various system parameters.

Recent research underscores the importance of MIMO modeling for HVAC systems, as it provides a more comprehensive representation of their complex internal interactions. MIMO modelling has been shown to provide a more comprehensive representation of the complex interactions within HVAC systems (Abdo-Allah et al., 2018; Abdo-Allah et al., 2017). State-space approaches have been effectively utilized to develop MIMO models for HVAC systems, facilitating more robust control system design and analysis (Abdo-Allah et al., 2018). Furthermore, advanced control strategies, including Direct Nyquist Array (DNA) and PID controllers, have been successfully implemented in MIMO HVAC models, leading to improved system performance and effective decoupling of system outputs (Touqan et al., 2022).

Most HVAC studies still rely on SISO models that overlook the strong inter-variable couplings of real systems. Although grey-box MIMO models offer richer dynamics, existing work is too computationally demanding for real-time control, highlighting a research gap. This research addresses this by extracting a validated high-order chiller model into a simplified grey-box 2 × 2 formulation, enabling practical multivariable control that reduces energy use, capturing the key dynamical system behaviour and enhances comfort in high-energy consumption systems.

2 Research methodology

This study builds upon the rigorously validated state-space model of a vapour compression chiller system established by (Yao et al., 2013). It will be developing a more concise, two-input, two-output multivariable model; this research accurately reflects the system’s dynamic behaviour. The model, which integrates essential components such as the condenser and evaporator based on fundamental principles of mass and energy conservation and heat transfer, aims to provide a computationally efficient yet realistic representation of chiller operations. It will indeed derive a simplified reduced-order model that still captures the same transient and steady-state behaviour of the original high-order model. Instead of presenting the detailed equations from the original work, this research focuses specifically on the sets of state equations and output equations only provided by (Yao et al., 2013), avoiding redundancy and duplication of effort.

2.1 Chiller system description

In this section, the chiller employed by (Yao et al., 2013) will be reviewed to identify the system components, inputs and outputs that will be selected for the simplified mathematical model process. It consists of four main components: the compressor, condenser, expansion device, and evaporator, as illustrated in Figure 1. Additionally, it includes several auxiliary elements, such as the chiller water pump. This figure represents the single-stage cycle, which serves as a standard model for practical refrigeration systems. Within this cycle, the refrigerant absorbs heat from the air-conditioned space via the evaporator, where its heat capacity is assumed to remain constant throughout the loop. Upon entering the compressor, the refrigerant’s pressure and temperature are significantly increased. In the condenser, the refrigerant undergoes a heat rejection process, during which its temperature decreases linearly due to heat transfer. The cycle concludes with the refrigerant passing through the expansion device, where it undergoes a pressure drop and further cooling.

FIGURE 1

In the refrigeration cycle, the chilled water (cold carrier fluid) enters the evaporator at a temperature after absorbing heat from the air-conditioned space. Inside the evaporator, the chilled water transfers this absorbed heat to the refrigerant, exiting at a reduced temperature . The cooled water then carries the cooling effect back to the indoor environment. It is subsequently routed to the Fan Coil Unit (FCU), where a blower fan distributes the cool air into the space. However, the FCU is not included in this model or analysis.

As the refrigerant moves through the compressor, its pressure and temperature rise significantly, preparing it for the condensation stage. In the condenser, the refrigerant releases the heat absorbed during the evaporation phase. This heat rejection is supported by the condenser coolant, which facilitates efficient thermal exchange. The inlet and outlet temperatures of the condenser coolant are represented by and , respectively. This heat transfer process gradually lowers the refrigerant’s temperature. The cycle is completed as the refrigerant passes through the expansion device, where its pressure drops and cools further, making it ready to re-enter the evaporator and commencing a new cycle.

2.2 Chiller mathematical model

A state-space model of the vapor compression chiller system provides a mathematical framework to describe the system’s dynamic behavior through its inputs, state variables, and outputs. This model is based on the fundamental physical principles governing the chiller, including mass and energy conservation. It is a representation in a structured matrix form, effectively capturing the interactions and dynamic relationships between the various system components.

The standard state space representation is given by the following pair of equations

Where:

is the chiller states vector and equals to

The square matrix represents the internal dynamics of the system by defining how the current state variables influence their rates of change. Its elements describe the interactions between different states, such as the effect of temperature at one stage on another within the system.

The input matrix shows how external inputs affect the rate of change of state variables. Its elements reflect the influence of inputs, such as refrigerant mass flow rate or coolant flow rate, on system states like temperatures.

The input vector that represents the external inputs to the system (e.g., external temperatures and fluid mass flow rates) and equals to .

The output vector represents the measurable outputs of the system and is equal to .

This is the output matrix that maps the state variables to the system outputs. It defines how the states directly influence the outputs. The elements of this matrix describe how the internal states of the chiller system influence measurable outputs.

This is the direct transmission matrix that maps the inputs directly to the outputs, accounting for any immediate or direct effect of inputs on outputs. The elements of this matrix describe how external inputs (e.g., compressor speed, valve position) directly and immediately affect these outputs.

Equations 3–31 present the sequential mathematical derivations that conclude the final simplified chiller multivariable input-output relationship in matrix form. Based on (Yao et al., 2013), the state space and output equations can be extracted and developed as follows:

Mapping the equations from Equations 3 to 15 into the state-space equations defined in (1) and (2) results in the matrix representation of the state-space model for the chiller system as follows:

It can be easily concluded that the state space matrices are

2.3 Transfer function matrix

To obtain the transfer function matrix, we take the Laplace transform of both the state and output equations, assuming zero initial conditions (i.e., ):

The Laplace transform of the state Equation 1 becomes:

The Laplace transform of the output Equation 2 becomes:where:

is the Laplace transform of the state vector , is the Laplace transform of the output vector , is the Laplace transform of the input vector , s is the Laplace variable.

Rearranging Equation 22 to solve for , we get:where is the identity matrix of appropriate dimensions, and is the inverse of the .

Substituting into the output Equation 23 yields:

The transfer function matrix is defined as the ratio of the output to the input in the Laplace domain:

Equation 26 is the transfer function matrix of the system, where each element of represents the transfer function from each input to each output. Substituting matrices and in Equation 26 aiming to calculate the transfer function matrix with the aid of MATLAB yieldswhere

So that input-output equation in matrix form is given by:

2.4 Chiller 2 inputs −2 outputs model

The chiller system detailed in (Yao et al., 2013), with five inputs and seven outputs, illustrates significant complexity. While such a model as open loop set up can provide detailed insights, it also poses challenges like increased computational demands and calibration difficulties when applying a multivariable control strategy. Simplifying the model by selecting fewer, strategically chosen inputs and outputs can reduce complexity while maintaining proper system functionality. The condenser and evaporator are key components in the chiller system, where heat absorption and rejection occur. Simplifying the model by focusing on just two inputs and two outputs related to these components is an effective way to reduce the system’s complexity. However, all the five system inputs can be incorporated in this configuration by considering them either as disturbances or held constant under application of a multivariable control strategy. This approach ensures that the model retains its core functionality while becoming more manageable for analysis and control, making it easier to work with while still capturing the essential thermal processes.

To simplify the model, the chosen input-output pairs are: the refrigerant mass flow rate , which impacts the chiller’s cooling capacity , and the mass flow rate of the chilled water , which affects the temperature of the chilled water . This selection is practical and follows the physical principles of chiller operation. The refrigerant mass flow rate, controlled by the voltage applied on the compressor, directly affects the cooling capacity by influencing the heat transfer efficiency in the evaporator. Similarly, the chilled water flow rate, controlled by the voltage applied on the chilled water pump, determines the outlet temperature by regulating the heat absorption rate. These outputs are measurable and critical for monitoring performance so that ensures the cooling demand is met, while indicates the cooling effect on the indoor thermal load.

It is worthy to mention that when applying advanced control techniques, the other inputs, namely, , can be considered as disturbances and can be held constant. It is also important to highlight that adding more inputs and outputs beyond a 2-inputs-2-outputs setup can make the system more complicated without providing significant performance benefits. Such a configuration would need additional sensors, actuators, and computational power, which might increase energy usage without improving the chiller’s main functions of heat rejection and absorption. Using a two-inputs, two-outputs approach concentrates on the most important variables, ensuring effective control while using fewer resources. This simple design lowers the chances of inefficiency or instability, helping to improve energy efficiency by focusing efforts where they are most beneficial.

Based on Equation 28, the two-input, two-output equations that define a multivariable chiller open loop as simplified model can be extracted as:

The multivariable input-output relationship in matrix form can be written as:

An open-loop system operates without feedback control. In this study, this helps isolate the chiller’s natural behavior before adding automated controls.

The term represents external disturbances that can impact heat transfer affecting the cooling capacity. These disturbances may arise from issues such as fouling, scaling, or dirt buildup on the evaporator surface. Additionally, also accounts for factors like a drop in compressor efficiency, which can lower the refrigerant mass flow rate and disrupt the overall heat transfer performance.

The final open loop simplified mathematical model describing the vapor compression liquid chiller in its effective reduced set is expressed in the following matrix equation:

This open loop configuration is described in the block diagram shown in Figure 2, which represents the open-loop system setup. In this system, the main manipulated inputs are and while the outputs are the chilled water outlet temperature, , and the chiller’s cooling capacity, .

FIGURE 2

As and are the main system manipulated inputs, the other inputs like and can be considered as additional disturbances that might affect both selected system outputs. can be held constant with a value of as a permanent input condition.

3 Results

This section evaluates the performance and significance of the derived open-loop chiller system model. Particular emphasis will be placed on the model’s behavior in response to perturbations in system inputs within the open-loop configuration and under various operational condition changes.

3.1 Importance of the developed 2-inputs 2-outputs open-loop model

The simplified mathematical model of the air conditioning chiller system, as shown in Equation 31, is crucial for engineers to develop control strategies in future studies that enhance system performance. As an open-loop model, it provides a straightforward view of the system’s behaviour without feedback effects. This foundational stage helps engineers understand how the system functions under various conditions and clarifies the relationship between inputs and outputs. With this understanding, they can design sophisticated control systems that aim to improve efficiency and overall performance.

Open-loop system responses are advantageous as they provide direct insights into the system’s behaviour, thereby facilitating the identification of inputs-such as refrigerant or coolant flow rates-that influence outputs like outlet temperature or cooling capacity. This model also helps pinpoint the key outputs that require monitoring and control. After understanding these aspects through open-loop modeling, engineers can design more sophisticated feedback control strategies to ensure the system runs efficiently and meets performance targets. Overall, the open-loop model serves as a strong foundation for developing responsive systems that optimise efficiency and boost performance.

3.2 Open loop responses

Open-loop responses involve analysing how the chiller’s outputs (e.g., outlet temperature and cooling capacity) react to changes in its inputs (e.g., flow rates) without any feedback control adjustments. This analysis helps in understanding the system’s natural dynamics, allowing engineers to design and optimise feedback controllers for better output regulation.

The block diagram described in Figure 2 was constructed using SIMULINK software, so when the open loop configuration is simulated with the initial conditions, including (; ; ; ; ), the output responses were as illustrated in Figure 3. Based on the initial conditions of input values, the first configured chiller output is approaching a steady state value of coupled with the second output approaching almost .

FIGURE 3

It is worth mentioning that Figure 3 will be the benchmark for analysing the performance of the system under perturbations from other inputs.

For different operational conditions, the system behaviour under alterations in system inputs can also be captured and illustrated using SIMULINK. Figure 4 shows the system behaviour when cold carrier mass flow rate (first input, ) is a step increase to , while the other inputs remain unchanged and adhere to their initial conditions. It can be noticed that there is a decrease in the exit cold carrier temperature from to which is caused by the increase of to .

FIGURE 4

For a different case, when the refrigerant mass flow rate (second input, ) is a step increase to , the resulting output responses are shown in Figure 5. Consequently, the exit cold carrier temperature is decreasing to approximately 17.3 with the coupling effect on the cooling capacity output increases to

FIGURE 5

For a case where a step disturbance input on inlet temperature of the condenser coolant is increasing to while all other inputs remain equal to their initial conditions, the output responses are depicted in Figure 6.

FIGURE 6

Based on the initial conditions of other inputs, the exit cold carrier temperature is approaching while the cooling capacity is approaching .

In a situation where a step disturbance input on the inlet temperature of the cold carrier fluid increases while all other inputs remain at their initial conditions, the output responses are shown in Figure 7. Based on the initial conditions of other inputs, the exit cold carrier temperature is approaching while the cooling capacity is approaching associated with some oscillations.

FIGURE 7

The last case that will be discussed is when the disturbance step of is applied due to dirt, scaling or foul in the evaporator while other inputs remain unchanged and as per their initial condition values. The results are illustrated in Figure 8. The responses show that no significant change in the exit cold carrier temperature ending with steady state value of and the cooling capacity is approaching steady state value of .

FIGURE 8

4 Discussions

This section provides comprehensive interpretations of the system responses illustrated in Figures 3–8, to explain the system’s behaviour within the open-loop configuration. These insights aim to augment the understanding of system engineers, thereby facilitating the proposal of suitable control strategies.

In Figure 3 where the system is subjected to initial conditions which are basically (; ; ; ; ), the first chiller output is approaching a steady state value of coupled with the second output approaching almost .

The decrease in temperature from to indicates that the evaporator is successfully absorbing heat from the cold carrier fluid, causing its temperature to drop. This temperature reduction is typical of an evaporator’s operation, where heat from the fluid is transferred to the refrigerant, cooling the cold carrier fluid as it passes through. The system reaches a stable cooling output, meaning that the heat transfer in the evaporator has balanced out with the conditions imposed by the inputs (e.g., refrigerant flow rate ). The cooling capacity of reflects the rate at which the chiller removes heat from the cold carrier fluid. This steady-state value indicates that the system has reached an equilibrium where the refrigerant is effectively absorbing heat from the cold carrier fluid at a rate of . The cooling capacity value shows that the chiller is operating at a relatively low load, typical for initial conditions with moderate refrigerant flow rates. As the cooling demand increases or the system parameters change (e.g., , would adjust accordingly.

Figure 3 will be the benchmark for analysing the performance of the system under perturbations from other inputs.

Figure 4 shows the system behaviour when cold carrier mass flow rate (first input, ) is a step increase to , while the other inputs remain unaltered and as per their initial conditions. The decrease in exit cold carrier temperature from to is caused by the increase of to . The difference is less than the previous case, when it was from to . As more fluid is passing through the evaporator, the same amount of released heat is spread across a larger mass of fluid. This means that the cold carrier fluid releases less heat per unit of mass, causing a slower temperature reduction in the cold carrier as it moves through the evaporator. In the same time, as more fluid is passing through the evaporator, the same amount of released heat is spread across a larger mass of fluid so the cooling capacity increases to approximately , up from the previous steady-state value of . The increase in does not only affect , but also impacts , demonstrating an inherent coupling between these two outputs.

In Figure 5 when the refrigerant mass flow rate (second input, ) is a step increase to the exit cold carrier temperature is decreasing to approximately 17.3 with coupling effect on the cooling capacity output increasing to . By increasing , the amount of refrigerant flowing through the evaporator increases. This allows the evaporator to absorb more heat from the cold carrier fluid (chilled water) causing a larger reduction in the temperature of the cold carrier fluid as it passes through the evaporator. The cooling capacity increases to 2.35 kW, up from the previous value of approximately 1.02 kW. The greater refrigerant flow rate allows the system to absorb more heat, which directly translates into a higher (cooling capacity). The system is now removing heat from the cold carrier at a higher rate to maintain the desired temperature decrease.

Similarly, in Figure 6, a step disturbance input on the inlet temperature of the condenser coolant is increasing to while all other inputs remain equal to their initial conditions, the exit cold carrier temperature is approaching while the cooling capacity is approaching . The increase in (inlet coolant temperature to the condenser) reduces the condenser’s ability to reject heat. This results in higher refrigerant pressures and temperatures within the system, leading to a less effective heat exchange process in the evaporator. The temperature rises to , up from its previous steady-state value of . With the evaporator unable to effectively absorb heat from the cold carrier fluid, the outlet temperature of the cold carrier fluid increases slightly, reflecting a reduced cooling effect. A negative indicates that instead of removing heat from the cold carrier, the chiller system is now adding heat. This happens because the condenser cannot reject enough heat due to the elevated , causing a thermal imbalance in the refrigerant cycle. With the refrigerant unable to condense efficiently in the condenser, the evaporator’s cooling capacity is severely compromised, leading to the system effectively transferring heat into the cold carrier rather than removing it.

In Figure 7, the exit cold carrier temperature is approaching while the cooling capacity is approaching associated with some oscillations. The increase in the inlet cold carrier temperature to imposes a significantly higher thermal load on the evaporator, which can be understood as social housing occupancy dynamics. In high-density residential environments, such as social housing in extreme climates, occupant density and intermittent usage patterns can cause significant fluctuations in thermal loads. The 40°C disturbance models scenarios like peak occupancy or severe weather events, where internal heat gains increase due to higher occupant activity. As a result, the evaporator absorbs more heat from the cold carrier fluid, but it cannot reduce the temperature as effectively due to the larger thermal input. The exit temperature stabilizes at , reflecting the system’s new equilibrium under the higher inlet temperature. The increase in raises the amount of heat the evaporator needs to absorb to maintain the cold carrier’s flow temperature. This causes the cooling capacity to increase significantly, as the evaporator works harder to remove the additional heat from the cold carrier fluid. The increase to reflects the system’s response to handle the higher thermal load introduced by the elevated inlet temperature.

The cooling capacity exhibits oscillations before reaching its steady-state value of . The evaporator’s thermal response is dynamic, and a sudden increase in creates an imbalance in the heat transfer process. This imbalance induces transient oscillations as the system attempts to reach a new equilibrium. Any overshoot in heat absorption by the evaporator leads to fluctuations in cooling capacity before stabilization. The chiller system’s physical components, such as the refrigerant cycle and heat exchanger, exhibit thermal and fluidic inertia. These inertial effects delay the system’s ability to settle immediately after the disturbance, resulting in oscillatory behavior. As increases, the heat absorbed by the evaporator rises, leading to simultaneous increases in and . This coupling is inherent in the chiller system’s design and reflects the interconnected nature of its heat transfer processes.

The results illustrated in Figure 8 due to disturbance applied on the system output show that no significant change in the exit cold carrier temperature ending with steady state value of but the cooling capacity is approaching steady state value of . Despite the disturbance, the cold carrier fluid’s exit temperature is unaffected, reducing the chilled water temperature from to . This indicates that the evaporator, though partially fouled, still manages to absorb sufficient heat to maintain the temperature difference required to cool the cold carrier fluid to . The absence of a significant temperature change could result from the reduced heat transfer being compensated by the evaporator’s thermal inertia or by the refrigerant still operating within its effective range despite the fouling. On the other hand, the fouling in the evaporator reduces its ability to transfer heat effectively, directly decreasing the cooling capacity. A reduction of corresponds to the disturbance magnitude, which accounts for the drop in to . This highlights the evaporator’s sensitivity to fouling and its direct impact on the system’s overall cooling performance. The reduced cooling capacity reflects a loss of energy efficiency, as the chiller is no longer operating at its full potential due to fouling. This could lead to higher operational costs and reduced performance in meeting cooling demands. This situation represents fouling or scaling in the evaporator, a common issue in social housing due to limited maintenance budgets and ageing infrastructure. Fouling reduces heat transfer efficiency, as evidenced by the drop in cooling capacity to without affecting outlet temperature . This scenario underscores the importance of low-maintenance design and robust control in social housing, where infrequent servicing can degrade system performance. By quantifying fouling’s impact, the model supports the development of predictive maintenance strategies or fault-tolerant controls, critical for cost-sensitive and resource-constrained environments.

The open-loop simulation results of the simplified mathematical model provide valuable insights into the dynamic behavior of the chiller system, highlighting its sensitivity to input changes and disturbances. The results show how variables like and respond passively, with phenomena such as coupling effects and oscillations, underscoring the system’s inability to maintain stability in some cases without feedback control. These findings emphasize the need for robust feedback control mechanisms to counteract disturbances, decouple interdependent outputs, and stabilize the system under varying conditions. Engineers can use this data along with the transfer function matrix in Equation 31 to design controllers, refine operational limits, and develop advanced strategies such as multivariable control techniques to enhance stability, reduce oscillations, and improve overall efficiency. By understanding these open-loop characteristics, system engineers are better equipped to optimize the chiller’s functionality and ensure reliable performance in real-world applications. The findings also highlight opportunities to integrate this model into sustainable social housing developments in the extreme climates areas, where adaptive HVAC systems can enhance energy resilience and occupant comfort.

5 Conclusion

This study presents a simplified multivariable state-space model of a vapor compression chiller system, emphasising the operation of key components such as the condenser and evaporator. The model expands upon an existing state-space mathematical framework, as detailed in (Yao et al., 2013), which has been validated through experimental data obtained from a chiller unit. While the original model in (Yao et al., 2013) involves a 5 × 7 multivariable structure with high-order dynamics, the current 2 × 2 model retains core thermal interactions while significantly reducing the state-space dimension, simplifying real-time control implementation. It delivers an accurate yet computationally straightforward mathematical model representation of the system’s behaviour. In control engineering, the mathematical model, represented through a state-space model, can be converted into a transfer function. The MIMO-related transfer function matrix plays a crucial role in the analysis of multivariable systems and controller design. It enables direct frequency-domain evaluation of input-output dynamics, such as system gain, phase lag, and coupling effects. For HVAC applications, especially vapor compression chillers, such a matrix is essential for designing decoupled or coordinated control strategies that regulate thermal outputs efficiently. The HVAC simplified mathematical model’s performance developed in this study is tested across different input scenarios using SIMULINK/MATLAB. The results are then plotted to validate the model’s ability to predict system responses. This includes scenarios such as operational changes, fouling, and thermal load variations. The open-loop simulation results demonstrate that the system reaches steady-state values after a certain response time, but the performance reveals significant interdependencies between outputs. For instance, the cooling capacity stabilizes at 1.02 kW, while the exit cold carrier temperature approaches 20.8°C under initial conditions, reflecting the system’s ability to reach equilibrium. The system reaches steady-state conditions in approximately 40 s for the exit cold carrier temperature and 70 s for the cooling capacity under nominal operating conditions. However, the simulation reveals strong input-output coupling, where a change in one input such as an increase in refrigerant mass flow rate affects both outputs simultaneously. For instance, a step increase in refrigerant flow from 0.116 kg/m3 to 0.216 kg/m3 not only lowers the cold carrier temperature to 17°C but also increases cooling capacity to 2.4 kW. Additionally, when subjected to a step disturbance in the condenser inlet temperature, the system exhibits oscillations in cooling capacity with a peak deviation of ±0.4 kW and a settling time exceeding 80 s. These dynamics highlight the system’s sensitivity to thermal and hydraulic disturbances, which would necessitate a multivariable control strategy that manages the system’s behaviour, aiming to enhance its performance while optimising its energy consumption. Therefore, implementing closed-loop multivariable control is essential. Such strategies can compensate for cross-coupling, dampen oscillations, and stabilize output responses in the presence of real-world disturbances and load variations, ensuring both performance and energy efficiency.

The novelty of this research lies in introducing a reduced-order multivariable model that retains essential dynamics while significantly lowering computational demands, making it suitable for cost-sensitive control applications. The simplified mathematical model can serve as the foundation for advanced feedback control designs that stabilise HVAC operation, decouple interacting outputs, and enhance energy efficiency, helping buildings meet sustainability targets without compromising comfort or reliability.

Future work could enhance the model’s accuracy by incorporating a secondary verification step to increase confidence. Since the model has been deemed accurate at the initial confidence level validated earlier with the Air Conditioning Experimental unit, future efforts might involve checking the simplified model’s accuracy at an additional confidence level. This can be achieved by testing the simplified mathematical model against a comparable experimental Air Conditioning unit, built based on the key parameters and components used in the model.

References

• 1

  Abdo-AllahA.IqbalT.PopeK. (2017). “Modeling and analysis of an HVAC system for the S.J. carew building at memorial university,” in Canadian Conference on Electrical and Computer Engineering, Windsor, ON, Canada, 30 April 2017 - 03 May 2017 (IEEE). 10.1109/CCECE.2017.7946689

  • CrossRef
  • Google Scholar

• 2

  Abdo-AllahA.IqbalT.PopeK. (2018). Modeling, analysis, and state feedback control design of a multizone HVAC system. J. Energy2018, 1–11. 10.1155/2018/4303580

  • CrossRef
  • Google Scholar

• 3

  AddisB. (2020). Past, current and future use of physical models in civil engineering design. Proc. Institution Civ. Eng. Civ. Eng.174 (2), 61–70. 10.1680/jcien.20.00028

  • CrossRef
  • Google Scholar

• 4

  AframA.Janabi-SharifiF. (2014). Review of modeling methods for HVAC systems. Appl. Therm. Eng.67, 507–519. 10.1016/j.applthermaleng.2014.03.055

  • CrossRef
  • Google Scholar

• 5

  AliM. (2019). Developing a white-box modelling of residential air conditioning based on demand response algorithm. J. Build. Eng.10.13140/RG.2.2.26588.13445

  • CrossRef
  • Google Scholar

• 6

  AngelW. L. (2011). HVAC design sourcebook. New York, NY, USA: McGraw Hill Professional.

  • Google Scholar

• 7

  AsimN.BadieiM.MohammadM.RazaliH.RajabiA.Chin HawL.et al (2022). Sustainability of heating, ventilation and air-conditioning (HVAC) systems in buildings—an overview. Int. J. Environ. Res. Public Health19 (2), 1016. 10.3390/ijerph19021016

  • CrossRef
  • Google Scholar

• 8

  BacherP.MadsenH. (2011). Identifying suitable models for the heat dynamics of buildings. Energy Build.43 (7), 1511–1522. 10.1016/j.enbuild.2011.02.005

  • CrossRef
  • Google Scholar

• 9

  CastletonH. F.StovinV.BeckS. B. M.DavisonJ. B. (2010). Green roofs; building energy savings and the potential for retrofit. Energy Build.42 (10), 1582–1591. 10.1016/j.enbuild.2010.05.004

  • CrossRef
  • Google Scholar

• 10

  GoetzlerW.GuernseyM.YoungJ.FuhrmanJ.AbdelazizO. (2016). The future of air conditioning for buildings. Available online at: www.osti.gov/home/.

  • Google Scholar

• 11

  HassanM. R.Al-SamarraieS. A. (2022). Robust nonlinear control design for the HVAC system based on adaptive sliding mode control. J. Eur. Des Systèmes Automatisés55 (5), 593–601. 10.18280/jesa.550504

  • CrossRef
  • Google Scholar

• 12

  HomodR. Z. (2013). Review on the HVAC system modeling types and the shortcomings of their application. J. Energy2013, 1–10. 10.1155/2013/768632

  • CrossRef
  • Google Scholar

• 13

  JangH.KangB.ChoK.JangK.ParkS. (2019). Design and implementation of IoT-based HVAC and lighting system for energy saving. MATEC Web Conf.260, 02012. 10.1051/matecconf/201926002012

  • CrossRef
  • Google Scholar

• 14

  JinG. Y.TanP. Y.DingX. D.KohT. M. (2011). “Cooling coil unit dynamic control of in HVAC system,” in Proceedings of the 2011 6th IEEE Conference on Industrial Electronics and Applications, ICIEA 2011, Beijing, China, 21-23 June 2011 (IEEE), 942–947. 10.1109/ICIEA.2011.5975722

  • CrossRef
  • Google Scholar

• 15

  JorissenF.PicardD.SixK.HelsenL. (2021). Detailed white-box non-linear model predictive control for scalable building HVAC control. Proceedings of 14th Modelica Conference 2021. Linköping, Sweden, September 20-24, 2021, 315–323. 10.3384/ecp21181315

  • CrossRef
  • Google Scholar

• 16

  KathiravelR.ZhuS.FengH. (2024). LCA of net-zero energy residential buildings with different HVAC systems across Canadian climates: a BIM-Based fuzzy approach. Energy Build.306, 113905. 10.1016/J.ENBUILD.2024.113905

  • CrossRef
  • Google Scholar

• 17

  LapinskieneV.MartinaitisV. (2013). The framework of an optimization model for building envelope. Procedia Eng.57, 670–677. 10.1016/j.proeng.2013.04.085

  • CrossRef
  • Google Scholar

• 18

  LiJ.PoultonG.PlattG.WallJ.JamesG. (2010). Dynamic zone modelling for HVAC system control. Int. J. Model. Identif. Control9 (1–2), 5–14. 10.1504/IJMIC.2010.032354

  • CrossRef
  • Google Scholar

• 19

  LiP.QiaoH.LiY.SeemJ. E.WinklerJ.LiX. (2014). Recent advances in dynamic modeling of HVAC equipment. Part 1: equipment modeling. HVAC R Res.20 (1), 136–149. 10.1080/10789669.2013.836877

  • CrossRef
  • Google Scholar

• 20

  LiY.KarunathilakeD.VilathgamuwaD. M.MishraY.FarrellT. W.ChoiS. S.et al (2022). Model order reduction techniques for physics-based lithium-ion battery management: a survey. IEEE Ind. Electron. Mag.16 (3), 36–51. 10.1109/MIE.2021.3100318

  • CrossRef
  • Google Scholar

• 21

  MaY.MatuskoJ.BorrelliF. (2015). Stochastic model predictive control for building HVAC systems: complexity and conservatism. IEEE Trans. Control Syst. Technol.23 (1), 101–116. 10.1109/TCST.2014.2313736

  • CrossRef
  • Google Scholar

• 22

  MaZ.CooperP.DalyD.LedoL. (2012). Existing building retrofits: methodology and state-of-the-art. Energy Build.55, 889–902. 10.1016/j.enbuild.2012.08.018

  • CrossRef
  • Google Scholar

• 23

  MurphyM. D.O’sullivanP. D.da GraçaG. C.O’donovanA. (2021). Development, calibration and validation of an internal air temperature model for a naturally ventilated nearly zero energy building: Comparison of model types and calibration methods. Energies14 (4), 871. 10.3390/en14040871

  • CrossRef
  • Google Scholar

• 24

  Palmero-marreroA. I.OliveiraA. C. (2010). Effect of louver shading devices on building energy requirements. Appl. Energy87 (6), 2040–2049. 10.1016/j.apenergy.2009.11.020

  • CrossRef
  • Google Scholar

• 25

  PanY.ZhuM.LvY.YangY.LiangY.YinR.et al (2023). Building energy simulation and its application for building performance optimization: a review of methods, tools, and case studies. Adv. Appl. Energy10, 100135. 10.1016/J.ADAPEN.2023.100135

  • CrossRef
  • Google Scholar

• 26

  PavlakG. S.FloritaA. R.HenzeG. P.RajagopalanB. (2014). Comparison of traditional and Bayesian calibration techniques for gray-box modeling. J. Archit. Eng.20 (2). 10.1061/(asce)ae.1943-5568.0000145

  • CrossRef
  • Google Scholar

• 27

  SeyamS. (2018). Types of HVAC systems. HVAC Syst.10.5772/intechopen.78942

  • CrossRef
  • Google Scholar

• 28

  SilvestriA.LydonG.WaibelC.WuD.SchlueterA.SilvestriA.et al (2023). Data-driven reduced order modelling using clusters of thermal dynamics conference paper data-driven reduced order modelling using clusters of thermal dynamics. Build. Simul. Conf. Proc.18. 10.3929/ethz-b-000634990

  • CrossRef
  • Google Scholar

• 29

  SimpehE. K.PillayJ. P. G.NdihokubwayoR.NalumuD. J. (2022). Improving energy efficiency of HVAC systems in buildings: a review of best practices. Int. J. Build. Pathology Adapt.40 (2), 165–182. 10.1108/IJBPA-02-2021-0019

  • CrossRef
  • Google Scholar

• 30

  Subbaram NaiduD.RiegerC. G. (2011). Advanced control strategies for heating, ventilation, air-conditioning, and refrigeration systems - an overview: part I: hard control. HVAC R Res.17 (1), 2–21. 10.1080/10789669.2011.540942

  • CrossRef
  • Google Scholar

• 31

  TouqanB.Abdul-AmeerA.SalamehM. (2022). HVAC multivariable system modelling and control. Proc. Institution Mech. Eng. Part C J. Mech. Eng. Sci.237, 2049–2061. 10.1177/09544062221138836

  • CrossRef
  • Google Scholar

• 32

  TouqanB.AmeerA. A. (2024). Enhancing HVAC multivariable system performance through hybrid modeling and direct nyquist array control. Build. Serv. Eng. Res. Technol.45, 833–861. 10.1177/01436244241274107

  • CrossRef
  • Google Scholar

• 33

  TouqanB.TalebH.SalamehM. (2017). Low carbon cooling approach for the residences in the UAE: a case study in Dubai. World Congr. Civ. Struct. Environ. Eng.10.11159/awspt17.127

  • CrossRef
  • Google Scholar

• 34

  WhalleyR.Abdul-AmeerA. (2011). Heating, ventilation and air conditioning system modelling. Build. Environ.46 (3), 643–656. 10.1016/j.buildenv.2010.09.010

  • CrossRef
  • Google Scholar

• 35

  WuY.XueW. (2024). Data-driven weather forecasting and climate modeling from the perspective of development. Multidiscip. Digit. Publ. Inst. (MDPI)15 (Issue 6), 689. 10.3390/atmos15060689

  • CrossRef
  • Google Scholar

• 36

  YaoY.HuangM.ChenJ. (2013). State-space model for dynamic behavior of vapor compression liquid chiller. Int. J. Refrig.36 (8), 2128–2147. 10.1016/j.ijrefrig.2013.05.006

  • CrossRef
  • Google Scholar

• 37

  ZhangL. (2023). Challenges and opportunities of machine learning control in building operations. Springer Nature Link.

  • Google Scholar

• 38

  ZhaoJ.LiX.ShumC.McPheeJ. (2021). A review of physics-based and data-driven models for real-time control of polymer electrolyte membrane fuel cells. Energy AI6, 100114. 10.1016/J.EGYAI.2021.100114

  • CrossRef
  • Google Scholar

• 39

  ZhouD.HuQ.TomlinC. J. (2016). Model comparison of a data-driven and a physical model for simulating HVAC systems. Available online at: http://arxiv.org/abs/1603.05951.

  • Google Scholar

• 40

  ZhouD. P.HuQ.TomlinC. J. (2017). Quantitative comparison of data-driven and physics-based models for commercial building HVAC systems. Proceedings of the American Control Conference, 2900–2906. Seattle, WA, USA24-26 May 2017. 10.23919/ACC.2017.7963391

  • CrossRef
  • Google Scholar

• 41

  ZlatanovićI.GligorevićK.IvanovićS.RudonjaN. (2011). Energy-saving estimation model for hypermarket HVAC systems applications. Energy Build.43 (12), 3353–3359. 10.1016/j.enbuild.2011.08.035

  • CrossRef
  • Google Scholar