A New Model Predictive Control Method for Eliminating Hydraulic Oscillation and Dynamic Hydraulic Imbalance in a Complex Chilled Water System

: To enhance the energy performance of a central air-conditioning system, an effective control method for the chilled water system is always essential. However, it is a real challenge to distribute exact cooling energy to multiple terminal units in different ﬂoors via a complex chilled water network. To mitigate hydraulic imbalance in a complex chilled water system, many throttle valves and variable-speed pumps are installed, which are usually regulated by PID-based controllers. Due to the severe hydraulic coupling among the valves and pumps, the hydraulic oscillation phenomena often occur while using those feedback-based controllers. Based on a data-calibrated water distribution model which can accurately predict the hydraulic behaviors of a chilled water system, a new Model Predictive Control (MPC) method is proposed in this study. The proposed method is validated by a real-life chilled water system in a 22-ﬂoor hotel. By the proposed method, the valves and pumps can be regulated safely without any hydraulic oscillations. Simultaneously, the hydraulic imbalance among different ﬂoors is also eliminated, which can save 23.3% electricity consumption of the pumps.


Introduction
The energy consumptions in the buildings sector account for almost 40% of the total final energy consumption [1,2]. Previous studies showed that heating, ventilation, and air conditioning (HVAC) systems take up almost 50% of energy use in buildings [3]. To enhance energy savings in HVAC systems, many solutions have been applied in the last few decades. Obviously, the control methods are effective approaches with minimal additional cost. With the decreased costs of Data Collecting System (DCS), Cloud Service (CS), and Remote Control (RC), etc. in recent years, the design and implementation of more complex control techniques have become feasible [4][5][6].
Classical on/off and PID control methods are still widely used in many HVAC systems due to simplicity and lower investment, causing inconsistent performance among these systems. With the help of the modern automation and information technologies, many advanced control methods can be adopted to improve the safety and efficiency of the HVAC systems. Model Predictive Control (MPC) is one of the promising control approaches for future low energy building. Research on the MPC method has attracted a lot of attention in the past few years due to its many advantages, such as: performing anticipatory control instead of corrective control [7] and ability to deal with constraints and uncertainties [8]. Recently, Abdul and Farrokh [9] presented a comprehensive literature review of control methods with an emphasis on the theory and applications of MPC for HVAC systems. Recently, Abdul and Farrokh [9] presented a comprehensive literature review of control methods, focusing on the theory and applications of MPC for HVAC systems. Compared with the PID, Robust and Fuzzy control methods, MPC control strategies chillers to achieve high efficiency. Recently, Gianni et al. [35] proposed a specialized MPC approach to control the HVAC system and the storage devices under thermal comfort and technological limitations.
In this paper, the contributions are mainly in two aspects: (1) A new MPC method is proposed to achieve the dynamic hydraulic balance for the chilled water distribution system. In this procedure, a system model is developed with considering the hydraulic characteristics of the chillers, terminal units, throttle valves, pumps and pipe network in detail. An optimization model using a GA algorithm is also presented to provide the calibrated parameters for the system model; (2) two types of MPC-based controllers at the local level are elaborately designed for adjustments of throttle valve and variable-speed pump, respectively. Instead of those PID-based controllers, the local MPC-based controllers can work more smoothly with the supervisory MPC platform. In addition, to validate the proposed MPC method, a round of regulation is demonstrated step by step in a real-life case study. With the help of the automation system and information technologies, the dynamic hydraulic balance was achieved without any hydraulic oscillations in the tested chilled water system.

System Description and Two Control Mechanism
For a complex chilled water system, how to distribute the exact cooling energy to multiple terminal units depends on accurate load prediction and an effective control method. There are already many previous studies dedicated to air-conditioning load prediction [36,37]. In this paper, we focus on the safety and efficiency of the new control method.

Chilled Water Distribution System
A typical chilled water distribution system is composed of chillers, terminal units (such as AHU or Coil), pumps, valves and numerous pipes, etc. A diagram of a typical configuration of chilled water distribution system in a high-rise building is shown in Figure 1.
The chillers provide cooling energy to circulating water. There are also heat rejection units for the chillers, which are not shown in Figure 1. The primary pumps that serve those chillers usually operate at a constant speed. The supplementary pump and water tank are used to provide a stable pressure head for the inlet of the primary pumps. At the outlet of the chillers, a common leg is employed to guarantee the chillers with enough and stable water flow. The secondary pumps are often operated at variable speed according to the cooling demands. The throttle valves are installed at the entrance of each floor, which are adjusted by the controllers to achieve a hydraulic balance between different floors. There are also many on/off valves installed at the entrance of the terminal units, which can be controlled by a local user or an upper supervisory controller.

Feedback Control Method
In a conventional control method, the valves are usually tuned by PID-based controllers. A schematic control diagram of typical PID control method is shown in Figure 2.
In a typical PID-based control method, the openness of the valves and rotation speed of the pumps are adjusted by those PID controllers. The set point for a controller is vitally important. The set-points come from the upper supervisory platform or onsite user. Usually, one or several indoor temperature sensors can be employed to provide the current indoor temperature (T indoor ). The temperature errors between current points and set points are used as feedback signals for the controllers to adjust valves. The return water temperature (T r ) can be used to regulate pumps as well. The feedback mechanism can usually work well when there are not too many controllers with coupling control aims. Unfortunately, there is severe nonlinear hydraulic coupling between the flow rates of the valves and pumps in a complex chilled water distribution system. For instance, when one of the valves is adjusted to a new set point, the flow rates of the other valves would be influenced instantly. Then, the PID controllers must act on their valves to keep their own set points. The adjustments of the valves would thereafter influence each other continuously, which would cause hydraulic oscillations in the pipe network. This phenomenon can significantly deteriorate the control performance. When the variable-speed pumps are also involved in this control process, the hydraulic oscillations would be aggravated. Sometimes, hydraulic oscillations may cause the water hammer in the pipe network, which can lead to leaks or the rupturing of the pipes. To mitigate the hydraulic oscillations, a trade-off between control quality and control stability has to be made for the feedback control methods.

Model Predictive Control Method
A new MPC method is introduced here, which is intrinsically immune to hydraulic oscillation while regulating valves and pumps in a complex water distribution system. A schematic control diagram of the proposed MPC method is shown in Figure 3.
In a typical PID-based control method, the openness of the valves and rotatio of the pumps are adjusted by those PID controllers. The set point for a controller important. The set-points come from the upper supervisory platform or onsite us ally, one or several indoor temperature sensors can be employed to provide the indoor temperature (Tindoor). The temperature errors between current points and s are used as feedback signals for the controllers to adjust valves. The return water ature (Tr) can be used to regulate pumps as well. The feedback mechanism can work well when there are not too many controllers with coupling control aims. U nately, there is severe nonlinear hydraulic coupling between the flow rates of th and pumps in a complex chilled water distribution system. For instance, when o valves is adjusted to a new set point, the flow rates of the other valves would be in instantly. Then, the PID controllers must act on their valves to keep their own se The adjustments of the valves would thereafter influence each other continuousl would cause hydraulic oscillations in the pipe network. This phenomenon can cantly deteriorate the control performance. When the variable-speed pumps are volved in this control process, the hydraulic oscillations would be aggravated. Som hydraulic oscillations may cause the water hammer in the pipe network, which to leaks or the rupturing of the pipes. To mitigate the hydraulic oscillations, a between control quality and control stability has to be made for the feedback contr ods.

Model Predictive Control Method
A new MPC method is introduced here, which is intrinsically immune to h oscillation while regulating valves and pumps in a complex water distribution sy schematic control diagram of the proposed MPC method is shown in Figure 3. In the proposed method, the valves and pumps are adjusted by so-called M trollers. There are no set-points for those MPC controllers. Instead, the simulatio from the chilled water distribution model are used for the MPC controllers to   In the proposed method, the valves and pumps are adjusted by so-called MPC controllers. There are no set-points for those MPC controllers. Instead, the simulation results from the chilled water distribution model are used for the MPC controllers to produce control signals to adjust the openness of valves or rotation speed of pumps. Without a feedback mechanism for the controllers, there is no possibility to produce hydraulic oscillation. To improve the accuracy of the simulation results, the parameters of the system model, such as hydraulic resistance of the pipes, chillers, and terminal units, should be calibrated regularly. Obviously, the proposed MPC control method is an intrinsically open loop control method. Therefore, the proposed method can be immune to any hydraulic oscillation. However, the effectiveness of the proposed method highly depends on an elaborate chilled water distribution model. A dedicated model for complex chilled water distribution system is presented in the next section. The details of the simulation results and process of the proposed method would be introduced in Section 5.

Chilled Water Distribution Model
The aim of the simulation is to obtain the required water pressure, flow rate, and temperature distribution in the pipe network. Then, the pumps and valves can be adjusted accordingly to supply proper water flow to the terminal units as demanded. The modeling for chilled water distribution system can be developed from two aspects: the hydraulic model and thermal model. As the water pressure waves in a pipe can propagate in sound speed, a chilled water system would reach hydraulic stability within a few minutes following one round of regulation. The calibration and simulation process can be performed adequately based on a static hydraulic model. On the other side, the supply water temperatures of the chillers are usually from 5 to 10 • C. The return water temperatures also change in a small range. Therefore, the adjustments of the valves and pumps depend mainly on the simulation results of the hydraulic model than that of the thermal model. The chilled water distribution modeling is focused on the hydraulic model with the following assumptions that can reduce the computational burden without significant decline of accuracy:

•
The friction factors between chilled water and pipe inner surface are temperature independent.

•
The chilled water is incompressible.

•
There is no leakage along the pipes.

•
The thermal losses from pipe insulation layer to ambient are neglected.
During hydraulic modeling, a chilled water distribution system can be treated analogous to an electric circuit [38]. Besides the pipe junctions, the chillers and terminal units which connect the supply and return pipes are regarded as the "nodes". The pipes with valves and fittings are regarded as the "branches".

Model Predictive Control Method
From the view of the hydronic distribution, the hydraulic behaviors of the chillers can be simplified as that of local hydraulic resistances. The relationship between a chiller's resistance R CH and its water flow rate . V CH can be regressed by using the measured data of onsite meters. Then, the hydraulic characteristic of a chiller can be written as Equation (1): where r 0 , r 1 , and r 2 are the fitting coefficients which can be calibrated with a few pairs of . V CH , R CH . Considering that the coefficients r 0 and r 1 have far less influence on resistance than r 2 , the coefficients r 0 and r 1 can often be neglected to reduce computation burden.
The supply water temperature of a chiller can be treated as an already known value from onsite temperature sensors or set temperature value of the chiller in operation, which can be written as Equation (2): where T s is the supply water temperature of a chiller; and T set CH is the set value of the chiller in operation.
The primary pumps are often operated at constant speed. Then, the hydraulic characteristic of a primary pump can be given as Equation (3): where ∆H P1 is the water head of the pump; . V P1 is the water flow rate of the pump; k 0 , k 1 and k 2 are the fitting coefficients which are usually provided by the pump manufacturer; Subscript "P1" indicates a primary pump.
At the inlet of the primary pump, there are often a supplementary pump and a water tank to provide a designed static pressure point in case of water leakage or vaporization. Then, the return water pressure P r at the inlet of the primary pumps can be given as Equation (4): where P set is the designed static pressure. The secondary pumps are usually operated at variable speed. The hydraulic characteristic of a variable-speed pump can be given as Equation (5): where n P2 is the current rotation speed of the pump; n 0 P2 is the rated rotation speed of the pump; Fr P2 is the operation frequency corresponding to n P2 ; Fr 0 P2 is the rated operation frequency corresponding to n 0 P2 ; k 0 , k 1 , and k 2 are still the fitting coefficients of the pump.

Terminal Unit and Valve
The terminal units in a chilled water system can also be treated as local resistances as that of the chillers. In the calibration process, r 0 , r 1 , and r 2 are still the fitting coefficients for hydraulic resistance of a terminal unit. However, the flow rate of a terminal unit should be regulated to fulfill its customers' cooling energy demand. In the simulation process, the terminal units should be treated as nodes with required mass flow rates . m TU , which is given as Equation (6): where Q TU is the predicted cooling load of the terminal unit; T s TU is the supply water temperature of the terminal unit; and T r TU is the return water temperature of the terminal unit. Q TU is usually estimated from a cooling load prediction software tool. Considering that the insulation layers of the chilled water pipes are often effective to prevent the heat transferring to ambient, T s TU is approximately equal to the supply water temperature of the chillers. T r TU depends on the onsite conditions and effectiveness of the terminal unit. For instance, a coil in good condition can reach the designed ∆T syndrome with required mass flow rates . m TU . The designed ∆T is normally provided by the manufactory, such as 5.0 • C. Then, T r TU can be selected to be T s TU − 5 . When the required mass flow rate of a terminal unit had been determined by Equation (6), the corresponding pressure drop of this terminal unit can be obtained by the simulation process.
The pressure drop of a throttle valve ∆P V can be estimated from its hydraulic characteristics as Equation (7): Energies 2021, 14, 3608 where . V V is the water flow rate of the valve; ρ v is the water density at 16 • C; ρ is the water density by the valve; K v is the flow coefficient of the valve, which can be given as Equation (8): K v = f (ϕ) (8) where ϕ is the openness of the valve; and the relation between K v and ϕ is provided by the manufacturer.
When the required valve's pressure drop and flow rate have been obtained by hydraulic simulation, the required openness of the valve can be obtained by Equation (7) and Equation (8).

Pipe Network
The hydraulic model of a chilled water pipe network can also be treated analogous to that of an electricity circuit [38]. Specifically, the pressure drop, flow, and pipeline resistance characteristic coefficient can be regarded as the voltage, current, and electric resistance. Considering the topology of a pipe network layout, k is the total branch number, and (n + 1) is its total node number. Then, the associated matrix A and the basic circuit matrix B of the pipe network can be obtained from graph theory.
According to the Kirchhoff's current law, it can be written as Equation (9): where . V b is the flow rate column vector of each branch. From the Kirchhoff's voltage law, it can be written as Equation (10), where R b is the resistance column vector of each branch; H p is the pump head column vector of each branch; and Z b is the height difference column vector for each branch. The resistance of a pipe can be estimated by the equation of Darcy [39] as Equation (11): where R is the resistance of the pipe; l is the length of the pipe; d is the inner diameter of the pipe; . V is the volume flow rate of the pipe. f is friction factor of the pipe, which can be estimated by the equation of Colebrook-White [40] as Equation (12), where ε is the absolute roughness of the inner surface of the pipe; Re is the Reynold number.

Calibration Model
The resistances of the pipes, chillers, and terminal units, which are necessary parameters for hydraulic simulations, must be calibrated firstly. A valid calibration model can find a set of optimal resistance values by the measurements. Based on the calibrated parameters, the simulation model can reach an acceptable match between the measured and the estimated data. Considering the non-analytical and non-smooth characteristics of the decision variables, the calibration model in the paper is also developed with an optimization method, as shown in Table 1.
In the objective function, Equation (13)

Genetic Algorithm Solution
The genetic algorithm (GA) is employed to get a global optimization result considering the non-analytical and non-smooth properties of a pipe network hydraulic model. During the GA processing, the decision variables will be discretized in their range. For instance, the roughness of the pipes ε k can be discretized into 2 10 different values between its upper limit (i.e., 1 mm) and lower limit (i.e., 0.01 mm). For instance, we made the roughness of the pipes ε k discrete into 2 10 different values between its upper limit (i.e., 1 mm) and lower limit (i.e., 0.01 mm). The resolution of ∆ε k can be obtained by Equation (19): Such a resolution ∆ε k may be precise enough to meet the calibration requirement of most chilled water systems. The other two decision variables R i , R j can be discretized similarly. For instance, R i,min and R i,max are set to be 1 kPa and 80 kPa, respectively. R j,min and R j,max are set to be 1 kPa and 50 kPa, respectively. After R i , R j are calibrated, the fitting parameters r 0 , r 1 , and r 2 can be regressed as Equation (1). Figure 4 shows the scheme of the GA algorithm for the calibration model. The number of generations are defined as variable "Gen". The GA algorithm would iterate until the value of "Gen" reaches the maximum number "MaxGen". Due to the sophisticated mechanisms of the GA algorithm, this paper does not discuss the impact of the population size, probability of crossover, and probability of mutation on the algorithm.

Design and Operation
Based on the advanced automation system of the chilled water system, the design of the local MPC controllers and the onsite operation approach are presented in detail in this section.

Local MPC Controller Design
The proposed MPC controllers are used to adjust the openness of the throttle valves and variable-speed pumps. With regard to an MPC controller of a valve, the openness can be calculated by its pressure drop and flow rate. According to Equation (7), it is rewritten as Equation (20): Combining Equation (20) with Equation (8), it is given as Equation (21): where f −1 is the reverse function of Equation (8); the relation between ϕ and K v depends on the hydraulic characteristic of the valve. It is often expressed as a nonlinear function as Equation (22), where a 1 and a 0 are the fitting coefficients provided by the manufacturer. When the simulation results of ∆P V and . V V are obtained by the proposed hydraulic model, the MPC controller can calculate the expected openness of the valve instantly by Equation (21) or Equation (22). Then, the openness of the valve can be adjusted by an accurate signal from its actuator.
With regard to an MPC controller of a variable-speed pump, the rotation speed can also be calculated by its water head and flow rate. According to Equation (5), it is given as Equation (23), n P2 n 0 where a P , b P , and c P are the factors as follows: Since the pump characteristic curve (∆H P2 -. V P2 ) is an inverted parabola, k 0 is a negative number, and so is a P . Therefore, the ± sign in Equation (23) means that the positive sign is for all points to the right of the curve maxima, while the negative sign applies to points to the left of it. In some situations, there are two possible rotation speeds corresponding to the same water head of the pump. It can lead to physical instabilities for the pump to alternate between two rotation speeds. The MPC controller of pump should choose the value closer to the current rotation speed.
Similarly, when the simulation results of ∆H P2 and . V P2 are obtained by the proposed hydraulic model, the MPC controller can calculate the expected operation frequency of the pump by Equations (23) and (5). Then, the operation frequency of the pump can be adjusted by an accurate signal from its invertor.

Real-Time Regulation Approach
With the help of automation and information technologies, such as the Building Information Modeling (BIM) system, Building Automation (BA) system, Cloud Service (CS) and Remote Control (RC), the proposed MPC method can be implemented smoothly. The real-time control process is usually following the hourly cooling load prediction. Therefore, the period of the real-time regulations on the valves and pumps in a chilled water system is also 1 h accordingly. To realize the proposed method, there are four successive steps in a round of regulation, as shown in Figure 5: (1) Real-time data collection The current hydraulic state of the chilled water system should be collected by an effective data access system, such as the Building Automation (BA) system. The network layout, location of chillers and terminal units, integrity of pipes and pumps, etc. can be recorded by a database system, such as BIM or CAD-based tools. The valid data are treated as the initial and boundary conditions, which are the foundation for the hydraulic calibration and simulation.
(2) Hydronic distribution system calibration The onsite meters are installed at selected positions of the water pipe network. For instance, the entrance of each floor, the inlet and outlet of the chillers, valves, pumps, and some key terminal units are usually selected to install pressure meters, flow rate meters, or temperature sensors. Based on the measured data, the resistance of pipes, chillers, and terminal units can be calibrated by the proposed calibration method. Generally, the accuracy of a calibration model is positively correlated with the amount of valid data. The calibration process can be performed in Cloud Service or the local server. CS is more adapted to serve many buildings simultaneously. It is more economic and efficient in recent years to deploy all algorithms and operation data in the public cloud.

(3) Hydraulic simulation
With the calibrated parameters in step (2), an hourly hydraulic simulation on the chilled water system can be carried out based on the hourly cooling load prediction. The pressure drops and flow rates of the throttle valves and the water heads and flow rates of the variable-speed pumps are calculated to achieve hydraulic balance for terminal users in different floors. The simulation process can also be performed in CS. Then, the simulation results can be transferred to the MPC controllers via internet or local network.

(4) Pump and valve adjustment
With the values of pressure and flow rate from the hydraulic simulation, the MPC controller can produce the corresponding control signals to adjust valves or pumps. In order to avoid the hydraulic oscillation, all valves and pumps should be adjusted to the latest values simultaneously. To fulfill this task, an RC system can be used to adjust the pumps and valves in different locations.
After a round of regulation, water flow in pipes would stabilize in one or two minutes. Then, the control effect would be evaluated instantly. The next round of calibration process would use the collected pressures and flow rates as new input data to obtain more accurate calibration parameters. Because the errors in the calibration and simulation process would be improved in the next round of regulation, sometimes, the onsite regulation had to be performed for 2-3 rounds to meet control expectation as close as possible. Generally, the errors between observed and expected flow rates are good criteria to determine whether the control expectation has been met. The criterion for the onsite regulation is given as Equation (24): whereV e n is the expected flow rate;V o n is the observed flow rate of onsite meter; max is the allowable relative error; n is the relative error of flow meter n. Subscript n is the index of the flow meter; N is the total assembly of the flow meters.

Case Study
To validate the proposed MPC method, the control system of the chilled water system in a hotel in Shanghai was retrofitted in April of 2019. The topology diagram of the chilled water distribution system is shown in Figure 6. There are 22 floors in this high-rise building with several Air Handling Units (AHUs) in each floor. To overcome the severe hydraulic imbalance among different floors, 22 throttle valves with MPC controllers are installed in the supply water pipe at the entrance of each floor. There are two identical chillers in operation. The rated cooling capacity of the chiller is 865.6 kW, and the rated power consumption is 245.2 kW. The designed temperature of chilled water supply/return is 7/12 • C. The rated volume flow rate of the primary pump is 90.0 m 3 /h, and the rated water head is 27.9 mH 2 O. There are six identical secondary pumps with a rated volume flow rate as 55.0 m 3 /h and a rated water head as 30.7 mH 2 O. There are six identical secondary  The lengths and diameters of pipes and locations of the chillers, valves, pumps, and terminal units had been obtained from the CAD drawings of the project. To simulate the water distribution system, two chillers, 22 throttle valves, three primary pumps, six secondary pumps, 202 terminal units, and 946 pipes are involved in the simulation model. The adjustments of valves and pumps are performed hourly according to the hourly cooling load prediction. The cooling load prediction was provided from a third-party software based on the well-known tool "Energy-plus", which considered local climate and occupants (1) Step 1.
The real-time observed data of 27 pressure meters and 24 flow meters were collected every 30 s by wireless signals. Then, 120 groups of observed values of pressure and flow rate can be collected in an hour. However, the hydraulic conditions between two rounds of adjustments (in an hour) were usually similar. Then, the collected data each 15 min were averaged as only one sample for calibration. In this case, the observed date in recent 72 h (288 samples) was used to calibrate the hydraulic resistances of pipes, chillers, and AHUs. The primary pumps are running at a constant speed (50 Hz), while the secondary pumps are running at variable speed according to the control strategy. The hydraulic characteristic curves of the primary and secondary pumps are shown in Figures 7 and 8, respectively.  (2) Step 2.
In step 2, the hydraulic resistances of the chillers and terminal units, and the roughness of the inner surface of the pipes were calibrated by the proposed method. According to Equation (13), W = 288, M = 27, and N = 24. The initial values of all the roughness of the pipes were set as 100 (µm). The initial values of resistance coefficients of all chillers (r 0 , r 1 , r 2 ) were set as (0, 0, 0.80). That of all terminal units was set as (0, 0, 8.0). Coefficients r 0 and r 1 were neglected for simplification. Using the GA algorithm, the decision variables were coded into 100 chromosomes for the initial population. In addition, the "MaxGen" was set to 100 in the case. The calibrated roughness values of all pipes are shown in Figure 9. The average roughness of all pipes is 141.1 (µm). The calibrated resistance coefficient values of all terminal units are shown in Figure 10. The average resistance coefficient r 2 of all terminal units is 25.8. The calibrated local resistance coefficient r 2 of 2 chillers are 0.940 and 0.958, respectively.  With the latest calibrated parameters, the simulation step can be performed based on the current cooling load prediction. The required flow rates of terminal units can be obtained by Equation (6). Then, the water pressure and flow rate distribution in the pipe network can be simulated accordingly. The simulation program is developed by the authors of this paper on the Microsoft Visual studio 2017 platform. One simulation took less than 2 min. The simulation results can be illustrated clearly as color-coded maps according to the construction CAD drawings. For simplicity, not all 24 simulation results from 12:00 a.m. to 11:00 p.m. are illustrated there. The simulation results of the water pressure and flow rate distribution at 12:00 a.m. are shown in Figures 11 and 12, respectively. The required pressure drops and flow rates of the throttle valves are shown in Table 2. The simulation results of the valves would be transferred to the MPC controllers via wireless signal either. The required pressure head and flow rate of each secondary pump is 31.05 mH 2 O and 201.5 m 3 /h, respectively. To save electricity as much as possible, the pressure head of the secondary pump should be set at a low value for which at least one of the 22 throttle valves consumes a water pressure drop of no more than 2.0 kPa. The lowest pressure head of the secondary pump can be obtained in 2-3 tries.
When the MPC controllers of the 22 valves and six pumps obtained the latest simulation results, the openness of valves and frequency of the pumps can be adjusted by Equations (22) and (23), respectively. The openness of the valves at 12 o'clock is also shown in Table 2. The frequency of the secondary pump is 46.846 Hz accordingly.
From above steps 1 to 4, a round of regulation at 12 a.m. was accomplished. After 3-5 min of this regulation finished, the observed data of the flow meters were used to determine whether the control expectation has been met according to Equation (24). As shown in Figure 13, the maximal relative error between observed and expected flow rates is less than 1.0%. It indicated that the control aim had been achieved after this round of regulation. The next round of regulation would be performed according to the latest cooling load prediction at 1:00 p.m. During the regulations by the proposed MPC method, the hydraulic oscillations in the chilled water system had been eliminated. Considering the large thermal inertia of the building, the hourly regulations are usually enough for the chilled water system in such a 22-floor hotel. Furthermore, the dynamic hydraulic imbalance is almost eliminated after the proposed regulations. The hydraulic balance can significantly improve the system efficiency and hotel guests' comfort. In a conventional PID control method, the flow rate of the chilled water is regulated to fulfill the cooling demand of the least favorable circulating loop. The chilled water systems are often involved in a large flow rate and degrading ∆T syndrome in most PID-based control systems. Then, the terminal units in most floors are in an over-flow state, which causes excessive thermal waste and electricity consumption of the pumps. To evaluate the energy saving effects of the proposed (MPC) method, the hourly electricity consumed by the circulation pumps are simulated in conditions of the proposed method and the conventional method, respectively. The hourly simulation results in a typical summer day are shown in Figure 14. In this case, the average electricity saving of the proposed MPC method is as high as 23.3% in 24 h. It is noted that the electricity consumed by the chillers or pumps cannot be measured separately from other electricity devices by the onsite meters. Therefore, the measured electricity data of the pumps have to be absent in Figure 14. On the other hand, the actual energy consumptions of the terminal units in an over-flow state are not easy to estimate accurately in the simulation scenario by the conventional method. Therefore, the energy consumption of the chillers cannot be estimated accurately as well. As a rule of thumb, the real-life chillers can save more than 10% electricity consumption when the large flow rate and degrading ∆T syndrome are eliminated.

Conclusions
This paper proposed a safe and highly efficient MPC method for the chilled water system with the help of advanced automatic and information technologies including BIM, BA, CS, and RC. In a hotel of 22 floors, the proposed method is illustrated in detail about how to adjust valves and pumps step by step. The main conclusions of this paper are as follows: (1) In a complex chilled water system with throttle valves and variable-speed pumps, a new MPC method is proposed for real-time operation to achieve dynamic hydraulic balance and immune hydraulic oscillation. By the proposed method, the MPC controller is designed instead of the conventional PID controller. A round of regulation encompasses four steps, which are (a) Real-time data collection; (b) Hydronic distribution system calibration; (c) Hydraulic simulation; and (d) Pump and valve adjustment. After a round of regulation, the control effect would be evaluated instantly. It can be improved to meet the control expectation as close as possible in the next round of regulation. (2) A chilled water distribution model is developed to simulate a complex chilled water system in a high-rise building. The CAD drawings or BIM diagram of a real-life project can be used to produce the topology of the pipe network layout. The required water pressures and flow rates of the valves and variable-speed pumps can be obtained according to the cooling load demands by the proposed hydraulic model. (3) A calibration method is presented with a dedicated optimization model. The genetic algorithm is used to search the global optimal values of the hydraulic resistances (R i , R j , ε k ) with non-analytical and non-smooth characteristics. The data collection system is used to provide the calibration model with real-time measured pressure and flow rate data. (4) The proposed method got validated in a real chilled water system in a 22-floor hotel in Shanghai. The 22 throttle valves and six variable-speed pumps of the chilled water system were adjusted by the proposed MPC controllers in a typical summer day. The results indicated that the dynamic hydraulic balance was achieved without any hydraulic oscillations. The hourly electricity consumed by the circulation pumps are also compared between that of the proposed method and the conventional method. The average electricity savings of the proposed method is as high as 23.3% in a typical summer day.
Author Contributions: Y.Y. and N.Z. wrote this paper and developed the method. H.Z. and H.W. were responsible for data collecting and the validation of the proposed method. All authors have read and agreed to the published version of the manuscript.