- 1Faculty of Engineering and IT, The British University in Dubai, Dubai, United Arab Emirates
- 2Department of Architecture, College of Architecture, Art and Design, Ajman University, Ajman, United Arab Emirates
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.
Therefore, the research objectives can be derived from addressing the following research question: How can a simplified multivariable model enhance control efficiency and energy performance in HVAC systems for social housing? So, the objectives of the study are:
i. Deriving a reduced-order transfer function matrix from (Yao et al., 2013) validated state-space framework, focusing on refrigerant and chilled water flow rates as inputs. A transfer function, obtained through Laplace transforms, turns time-domain dynamics into an algebraic input-output relation. For the chiller it maps refrigerant-flow changes to cooling capacity, capturing delays and response shape.
ii. Validating the model’s open-loop dynamics under disturbances (e.g., fouling, load variations) through SIMULINK software to identify control challenges.
iii. Discussing its utility as a foundation for multivariable feedback control, emphasising decoupling and stability in extreme climates.
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.
In the refrigeration cycle, the chilled water (cold carrier fluid) enters the evaporator at a temperature
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
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:
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:
Rearranging Equation 22 to solve for
where
Substituting
The transfer function matrix
Equation 26 is the transfer function matrix of the system, where each element of
where
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
It is worthy to mention that when applying advanced control techniques, the other inputs, namely,
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
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
As
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 (
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,

Figure 4. Open loop responses when applying step change in
For a different case, when the refrigerant mass flow rate (second input,

Figure 5. Open-loop responses when the refrigerant mass flow rate is increased to
For a case where a step disturbance input on inlet temperature of the condenser coolant

Figure 6. Open loop responses when applying step increase at
Based on the initial conditions of other inputs, the exit cold carrier temperature
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

Figure 7. Open loop responses when applying step increase at
The last case that will be discussed is when the disturbance step of

Figure 8. Open loop responses when applying step disturbance at
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 decrease in
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,
In Figure 5 when the refrigerant mass flow rate (second input,
Similarly, in Figure 6, a step disturbance input on the inlet temperature of the condenser coolant
In Figure 7, the exit cold carrier temperature
The cooling capacity exhibits oscillations before reaching its steady-state value of
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
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
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.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
BT: Writing – original draft, Investigation, Conceptualization, Software. MS: Resources, Supervision, Writing – review and editing, Project administration.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. The author would like to recognize the support provided by Ajman University and the Healthy and Sustainable Built Environment Research Center for the publication of this research.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that Generative AI was used in the creation of this manuscript. The entire content is developed solely by the authors. Generative AI is used to improve the language and readability of the manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Abdo-Allah, A., Iqbal, T., and Pope, K. (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). doi:10.1109/CCECE.2017.7946689
Abdo-Allah, A., Iqbal, T., and Pope, K. (2018). Modeling, analysis, and state feedback control design of a multizone HVAC system. J. Energy 2018, 1–11. doi:10.1155/2018/4303580
Addis, B. (2020). Past, current and future use of physical models in civil engineering design. Proc. Institution Civ. Eng. Civ. Eng. 174 (2), 61–70. doi:10.1680/jcien.20.00028
Afram, A., and Janabi-Sharifi, F. (2014). Review of modeling methods for HVAC systems. Appl. Therm. Eng. 67, 507–519. doi:10.1016/j.applthermaleng.2014.03.055
Ali, M. (2019). Developing a white-box modelling of residential air conditioning based on demand response algorithm. J. Build. Eng. doi:10.13140/RG.2.2.26588.13445
Asim, N., Badiei, M., Mohammad, M., Razali, H., Rajabi, A., Chin Haw, L., et al. (2022). Sustainability of heating, ventilation and air-conditioning (HVAC) systems in buildings—an overview. Int. J. Environ. Res. Public Health 19 (2), 1016. doi:10.3390/ijerph19021016
Bacher, P., and Madsen, H. (2011). Identifying suitable models for the heat dynamics of buildings. Energy Build. 43 (7), 1511–1522. doi:10.1016/j.enbuild.2011.02.005
Castleton, H. F., Stovin, V., Beck, S. B. M., and Davison, J. B. (2010). Green roofs; building energy savings and the potential for retrofit. Energy Build. 42 (10), 1582–1591. doi:10.1016/j.enbuild.2010.05.004
Goetzler, W., Guernsey, M., Young, J., Fuhrman, J., and Abdelaziz, O. (2016). The future of air conditioning for buildings. Available online at: www.osti.gov/home/.
Hassan, M. R., and Al-Samarraie, S. A. (2022). Robust nonlinear control design for the HVAC system based on adaptive sliding mode control. J. Eur. Des Systèmes Automatisés 55 (5), 593–601. doi:10.18280/jesa.550504
Homod, R. Z. (2013). Review on the HVAC system modeling types and the shortcomings of their application. J. Energy 2013, 1–10. doi:10.1155/2013/768632
Jang, H., Kang, B., Cho, K., Jang, K., and Park, S. (2019). Design and implementation of IoT-based HVAC and lighting system for energy saving. MATEC Web Conf. 260, 02012. doi:10.1051/matecconf/201926002012
Jin, G. Y., Tan, P. Y., Ding, X. D., and Koh, T. 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. doi:10.1109/ICIEA.2011.5975722
Jorissen, F., Picard, D., Six, K., and Helsen, L. (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. doi:10.3384/ecp21181315
Kathiravel, R., Zhu, S., and Feng, H. (2024). LCA of net-zero energy residential buildings with different HVAC systems across Canadian climates: a BIM-Based fuzzy approach. Energy Build. 306, 113905. doi:10.1016/J.ENBUILD.2024.113905
Lapinskiene, V., and Martinaitis, V. (2013). The framework of an optimization model for building envelope. Procedia Eng. 57, 670–677. doi:10.1016/j.proeng.2013.04.085
Li, J., Poulton, G., Platt, G., Wall, J., and James, G. (2010). Dynamic zone modelling for HVAC system control. Int. J. Model. Identif. Control 9 (1–2), 5–14. doi:10.1504/IJMIC.2010.032354
Li, P., Qiao, H., Li, Y., Seem, J. E., Winkler, J., and Li, X. (2014). Recent advances in dynamic modeling of HVAC equipment. Part 1: equipment modeling. HVAC R Res. 20 (1), 136–149. doi:10.1080/10789669.2013.836877
Li, Y., Karunathilake, D., Vilathgamuwa, D. M., Mishra, Y., Farrell, T. W., Choi, S. 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. doi:10.1109/MIE.2021.3100318
Ma, Y., Matusko, J., and Borrelli, F. (2015). Stochastic model predictive control for building HVAC systems: complexity and conservatism. IEEE Trans. Control Syst. Technol. 23 (1), 101–116. doi:10.1109/TCST.2014.2313736
Ma, Z., Cooper, P., Daly, D., and Ledo, L. (2012). Existing building retrofits: methodology and state-of-the-art. Energy Build. 55, 889–902. doi:10.1016/j.enbuild.2012.08.018
Murphy, M. D., O’sullivan, P. D., da Graça, G. C., and O’donovan, A. (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. Energies 14 (4), 871. doi:10.3390/en14040871
Palmero-marrero, A. I., and Oliveira, A. C. (2010). Effect of louver shading devices on building energy requirements. Appl. Energy 87 (6), 2040–2049. doi:10.1016/j.apenergy.2009.11.020
Pan, Y., Zhu, M., Lv, Y., Yang, Y., Liang, Y., Yin, R., et al. (2023). Building energy simulation and its application for building performance optimization: a review of methods, tools, and case studies. Adv. Appl. Energy 10, 100135. doi:10.1016/J.ADAPEN.2023.100135
Pavlak, G. S., Florita, A. R., Henze, G. P., and Rajagopalan, B. (2014). Comparison of traditional and Bayesian calibration techniques for gray-box modeling. J. Archit. Eng. 20 (2). doi:10.1061/(asce)ae.1943-5568.0000145
Silvestri, A., Lydon, G., Waibel, C., Wu, D., Schlueter, A., Silvestri, A., 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. doi:10.3929/ethz-b-000634990
Simpeh, E. K., Pillay, J. P. G., Ndihokubwayo, R., and Nalumu, D. J. (2022). Improving energy efficiency of HVAC systems in buildings: a review of best practices. Int. J. Build. Pathology Adapt. 40 (2), 165–182. doi:10.1108/IJBPA-02-2021-0019
Subbaram Naidu, D., and Rieger, C. 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. doi:10.1080/10789669.2011.540942
Touqan, B., Abdul-Ameer, A., and Salameh, M. (2022). HVAC multivariable system modelling and control. Proc. Institution Mech. Eng. Part C J. Mech. Eng. Sci. 237, 2049–2061. doi:10.1177/09544062221138836
Touqan, B., and Ameer, A. A. (2024). Enhancing HVAC multivariable system performance through hybrid modeling and direct nyquist array control. Build. Serv. Eng. Res. Technol. 45, 833–861. doi:10.1177/01436244241274107
Touqan, B., Taleb, H., and Salameh, M. (2017). Low carbon cooling approach for the residences in the UAE: a case study in Dubai. World Congr. Civ. Struct. Environ. Eng. doi:10.11159/awspt17.127
Whalley, R., and Abdul-Ameer, A. (2011). Heating, ventilation and air conditioning system modelling. Build. Environ. 46 (3), 643–656. doi:10.1016/j.buildenv.2010.09.010
Wu, Y., and Xue, W. (2024). Data-driven weather forecasting and climate modeling from the perspective of development. Multidiscip. Digit. Publ. Inst. (MDPI) 15 (Issue 6), 689. doi:10.3390/atmos15060689
Yao, Y., Huang, M., and Chen, J. (2013). State-space model for dynamic behavior of vapor compression liquid chiller. Int. J. Refrig. 36 (8), 2128–2147. doi:10.1016/j.ijrefrig.2013.05.006
Zhang, L. (2023). Challenges and opportunities of machine learning control in building operations. Springer Nature Link.
Zhao, J., Li, X., Shum, C., and McPhee, J. (2021). A review of physics-based and data-driven models for real-time control of polymer electrolyte membrane fuel cells. Energy AI 6, 100114. doi:10.1016/J.EGYAI.2021.100114
Zhou, D., Hu, Q., and Tomlin, C. 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.
Zhou, D. P., Hu, Q., and Tomlin, C. 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, USA 24-26 May 2017. doi:10.23919/ACC.2017.7963391
Zlatanović, I., Gligorević, K., Ivanović, S., and Rudonja, N. (2011). Energy-saving estimation model for hypermarket HVAC systems applications. Energy Build. 43 (12), 3353–3359. doi:10.1016/j.enbuild.2011.08.035
Glossary
Nomenclature
Abbreviations
AC Air Conditioning
HVAC Heating Ventilation, and Air Conditioning
MIMO Multiple Input, Multiple Output
SISO Single Input, Single Output
COP Coefficient of Performance COP measures chiller efficiency by comparing cooling output (kW) to energy input (kW).
FCU Fan Coil Unit
IEA International Energy Agency
OECD Organisation for Economic Co-operation and Development
ML Machine Learning
AI Artificial Intelligence
MATLAB Matrix Laboratory (Simulation Software)
Keywords: air conditioning chiller, mathematical modelling, sustainable social housing, multivariable model, transfer function matrix, open-loop responses
Citation: Touqan B and Salameh M (2025) Mathematical modeling describing the performance of open loop multivariable vapor compression chiller model: implications for sustainable HVAC in social housing development. Front. Built Environ. 11:1624022. doi: 10.3389/fbuil.2025.1624022
Received: 07 May 2025; Accepted: 10 July 2025;
Published: 25 July 2025.
Edited by:
Hasim Altan, United Arab Emirates University, United Arab EmiratesReviewed by:
Fadi Alsouda, University of Technology Sydney, AustraliaAmel Soukaina Cherif, Tunis El Manar University, Tunisia
Copyright © 2025 Touqan and Salameh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Basim Touqan, YmFzZW0udHVxYW5AYnVpZC5hYy5hZQ==