Next Article in Journal
Comment on ‘Repurposing Hydrocarbon Wells for Geothermal Use in the UK: The Onshore Fields with the Greatest Potential. Watson et al. (2020)’
Next Article in Special Issue
Controlling the Thermal Power of a Wall Heating Panel with Heat Pipes by Changing the Mass Flowrate and Temperature of Supplying Water—Experimental Investigations
Previous Article in Journal
Effect of Demagnetization on a Consequent Pole IPM Synchronous Generator
Previous Article in Special Issue
Model-Based Condenser Fan Speed Optimization of Vapor Compression Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Exact Algorithm for the Optimal Chiller Loading Problem and Its Application to the Optimal Chiller Sequencing Problem

by
Federica Acerbi
1,
Mirco Rampazzo
2,* and
Giuseppe De Nicolao
1
1
Department of Industrial and Information Engineering, University of Pavia, 27100 Pavia, Italy
2
Department of Information Engineering, University of Padova, 35131 Padova, Italy
*
Author to whom correspondence should be addressed.
Energies 2020, 13(23), 6372; https://doi.org/10.3390/en13236372
Submission received: 25 October 2020 / Revised: 18 November 2020 / Accepted: 26 November 2020 / Published: 2 December 2020

Abstract

:
The optimal management of multiple chiller systems calls for the solution of the so-called optimal chiller loading (OCL) problem. Due to the interplay of continuous and logical constraints, OCL is an NP-hard problem, so that a variety of heuristic algorithms have been proposed in the literature. Herein, an algorithm for its exact solution, named X-OCL, is developed under the assumption that the chillers’ power consumption curves are quadratic. The proposed method hinges on a decomposition of the solution space so that the overall OCL problem is decomposed to a set of equality constrained quadratic programming problems that can be solved in closed form. By applying the new X-OCL solver to well known case studies, we assess and compare the performances of several literature algorithms, highlighting also some errors in the published results. Moreover, X-OCL is used to design a greedy optimal chiller sequencing (OCS) solver, called X-OCS. The X-OCS is tested on two literature benchmarks and on the model of the heating, ventilation and air-conditioning (HVAC) system of a semiconductor plant, over a two-year period. The performances of X-OCS are remarkably close to the theoretical optimal performance.

Graphical Abstract

1. Introduction

Heating, ventilation and air-conditioning (HVAC) systems typically utilize a large percentage of the total building energy consumption, amounting to approximately 25–30% for dwellings [1] and up to more than 50% in industries relying on clean-rooms [2]. Between 40% and 60% of this requirement is due to chillers, which are therefore responsible for a significant fraction of total energy use. Such numbers indicate that the efficiency of HVAC systems is closely dependent on the efficiency of the chiller unit. Since a multiple chiller system typically employs machines, there are usually several combinations of chillers’ part loads that are able to satisfy the load demand. The problem of determining the load fraction that each chiller has to deliver in order to minimize the system power consumption is known as optimal chiller loading (OCL) problem. In the last decade, several methods have been proposed to solve the OCL problem [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]. Chang et al. [3] assume that the chillers’ efficiency is a quadratic function of the partialization and maximize the sum of chillers’ coefficients of performance (COP) using the Lagrangian method. Of course, maximizing the sum of the efficiencies is not equivalent to the “canonical” OCL problem, where system power consumption is instead minimized. In a second paper, Chang [4] has addressed the canonical OCL problem again by the Langrangian multipliers methods, assuming that all the chiller power consumption curves are convex and cubic in the part load ratio. The poor convergence properties of the lambda iteration at low cooling load demand were obviated by a suitable gradient method.
Nevertheless, it is well-known that the chillers’ power consumption could be a concave function, which makes OCL an NP-hard mixed-integer problem and poses a further challenge to the search of the optimal solution as global convergence of iterative methods cannot be guaranteed. This has motivated the development of several heuristic methods which do not guarantee optimality, but are effective to obtain a fair solution in a viable execution time.
Assuming quadratic models of power consumption, Geem [5] resorts to the generalized reduced gradient (GRG) method. Salari et al. [6] show that the mixed integer problem can be solved using the general algebraic modeling system (GAMS). In parallel, many nature-inspired heuristic algorithms have been proposed, namely, genetic algorithm (GA) [7,8], simulated annealing (SA) [9,10], particle swarm optimization (PSO) [11,12], evolution strategy (ES) [13], differential evolution (DE) [14], cuckoo search algorithm using differential operator (DCSA) [15], differential search (DS) [16], improved firefly algorithm (IFA) [17], teaching-learning-based optimization (TLBO) [18], improved invasive weed optimization (EIWO) [19], distributed chaotic estimation of distribution algorithm (DCEDA) [20], imperialist competitive algorithm (ICA) [21] and wild goat’s algorithm (WGA) [22].
Many of these algorithms have been tested on the Hsinchu benchmark, a widely used case study consisting of a six-chiller system installed in a semiconductor factory located in the Hsinchu Scientic Garden (Taiwan). Hence, it is therefore natural to use the performances achieved on this benchmark to compare and rank alternative solution methods. A selection consisting of the six best algorithms is reported in Table 1. In order of publication, they are IFA (2013), DCSA (2014), GAMS (2015), TLBO (2017), EIWO (2018), DCEDA (2020). Two recent algorithms, ICA (2020) and WGA (2020), claim to have improved over the existing literature. However, an inspection of the chillers’ cooling capacities used in these two papers shows that they are different from those of the standard benchmark so that the claim of improvement is void. It is worth noting that, since an exact solver of the Hsinchu benchmark is not available, the final word on the optimality of solutions found by heuristic algorithms still needs to be had.
In the practical management of a real chiller system, solving the OCL is not sufficient, because there exist further dynamic constraints, namely minimum up-time and down-time requirements on chillers’ operation. When these constraints are accounted for, the power consumption minimization problem goes under the name of optimal chiller sequencing (OCS). Again, this is a problem that is hardly tractable without resorting to some heuristics. In particular, the knowledge of all future cooling loads is required, which raises the problem of forecasting it with reasonable accuracy, a task that can be successfully addressed only on a finite prediction horizon. In the literature, OCS solvers resulting from the combination of heuristic OCL algorithms with dynamic programming schemes have been proposed [8,23,24,25].
The present paper addresses both the OCL and OCS problems. Concerning the former one, two main issues are investigated. First of all, we derive an exact algorithm when the chillers’ power consumption is a quadratic function of the partial load ratio (PLR), as it happens for the Hsinchu benchmark. Our X-OCL algorithm allows one to have the final word on the existing heuristic methods, highlighting also some erroneous results reported in the literature. The second issue has to do with the practical applicability of X-OCL to real-world plants. In particular, the execution time is compared with a state-of-art mixed integer solver and the adequacy of the quadratic power consumption model is discussed.
Concerning the OCS problem, we exploit the X-OCL algorithm to derive a lower bound on the minimum power consumption achievable by any OCS solver. Second, a greedy OCS algorithm leveraging on X-OCL is proposed and compared with [23]. Finally, the lower bound is used to quantitatively assess the degree of suboptimality ensuing from the lack of preview implicit in the greedy approach.
Literature benchmarks, i.e., the Hsinchu one and two OCS benchmarks are used for test and comparison. Moreover, an extensive dataset collected during two years in a semiconductor fab is used to build and demonstrate a comprehensive solution of both OCL and OCS, including the data-based estimation of the chillers’ power consumption models. In particular, the potential energy saving with respect with the current HVAC energy management is assessed.

2. The Optimal Chiller Loading Problem

Assuming n chillers operated in parallel, let Q i , i = 1 , , n , denote the cooling power delivered by the i-th chiller and P i = P i ( Q i ) , i = 1 , , n the associated power consumption. For a prescribed overall cooling load demand Q l o a d , the goal of the OCL problem is finding the cooling powers Q i , i = 1 , , n that the n chillers have to deliver in order to minimize the system total energy consumption P t o t :
min Q 1 , , Q n i = 1 n P i s . t . Q l o a d = i = 1 n Q i
For each chiller, PLR (part load ratio) denotes the cooling load fraction, given by
P L R i = Q i Q 100 % , i
where Q 100 % , i is the maximum power supplied under full capacity operation. The vector of all PLRs is denoted by
PLR = [ P L R 1 , , P L R n ] T
When the i-th chiller is turned on, it should not operate under a minimum P L R , denoted by P L R m i n , i , 1 . For the subsequent derivations, it is convenient to introduce a binary variable δ i that indicates the status of the i-th chiller and a real-valued variable x i , such that:
P L R i = δ i · x i , δ i { 0 , 1 } , x i [ P L R m i n , i , 1 ]
The power consumption P i of the i-th chiller is assumed to depend mainly on P L R i and the condenser inlet water temperature T i , that is P i = P i ( P L R i , T i ) .
The consumption surface P i ( P L R i , T i ) is either obtained from laboratory experiments or field data collected during operation.
Note that the problem of finding the optimal part load ratios PLR * can be formulated as a mixed-integer nonlinear program (MINLP):
PLR * : = arg min PLR i = 1 n P i ( P L R i )
subject   to i = 1 n P L R i · Q 100 % , i = Q l o a d
P L R i = δ i · x i
P L R m i n x i 1
δ i { 0 , 1 }
In the above problem, two types of constraints are present: the cooling demand constraint, Equation (1b) and a set of operational constraints Equation (1c–e) regarding the admissible operating regions of the chillers.
In the following, it is implicitly assumed that the cooling demand constraint is such that the admissible solutions set is nonempty. Then, given that the constraints define a closed set of admissible solution, the cost function admits a minimum.

2.1. Quadratic Power Consumption Model

With the exception of particular structures, mixed-integer programming problems are classified as NP-hard, which means that in the worst case, the solution time grows at least exponentially with the problem size. Although its combinatorial nature might suggest the use of heuristics, we will show that a significant subclass of industrial OCL problems, characterized by a quadratic power consumption model quadratic in P L R i , may still be successfully attacked by a carefully designed exact method. It is worth asking whether the quadratic assumption of the power consumption is more or less realistic. We will discuss this point later. Even if the consumption is not quadratic, a good quadratic approximation is likely computable. Moreover, in the case in which it is not easy to workout a good quadratic approximation, it can still be addressed by means piece-wise quadratic model, as suggested in Section 6.
Assumption 1.
The power consumption P i of the i-th chiller obeys the following model:
P i ( P L R i , T i ) = β 0 , i + β 1 , i P L R i + β 2 , i P L R i 2 + f ( T i ) , if P L R m i n P L R i 1 0 , otherwise
where β p , i are the model parameters and f ( · ) is a suitable function. Moreover, it is assumed that β 2 , i 0 .
When the OCL problem is solved in a given time slot, the condenser inlet water temperature T i can be assumed to be known. Then, for a given T i , in the interval [ P L R m i n , 1 ] , the consumption surface is a quadratic function of P L R i alone:
P i ( P L R i , T i ) = a i + c i P L R i + q i P L R i 2 , if P L R m i n P L R i 1 0 , otherwise
with a i = β 0 , i + f ( T i ) , c i = β 1 , i , and  q i = β 2 , i .
In view of the quadratic nature of the cost function, the system total energy consumption can be expressed in matrix form as follows:
P t o t = i = 1 n P i = 1 2 x T Q x + c T x + a
where
x = { P L R i } , x R n × 1 a = { a i } , a R n × 1 c = { c i } , c R n × 1 Q = d i a g { 2 q i } , Q R n × n
Moreover, the equality constraint (1b) can be rewritten as
E x = d
where
E = { Q 100 % , i } , E R 1 × n d = Q l o a d , d R
Equation (4a,b) represent the standard OCL problem with quadratic power consumption, that so far has been attacked by either mixed integer programming solvers or a variety of heuristic algorithms. Differently from these approaches, below we derive an algorithm that computes the exact solution.

2.2. Partition of the Solution Space

In view of (4), the admissible set F for P L R i is
P L R i F , F = s = 0 3 F σ
F 0 = { 0 } , F 1 = { P L R m i n } ,
F 2 = { 1 } , F 3 = ( P L R m i n , 1 )
where each subset F σ is associated with one of the following operating conditions: switched off ( σ = 0 ), minimum part load ( σ = 1 ), maximum part load ( σ = 2 ), intermediate part load ( σ = 3 ). In the sequel, σ i { 0 , 1 , 2 , 3 } will denote the state of the i - th chiller.
In order to satisfy the operational constraints of a multiple chiller system, the solution of the OCL problem must be searched within the cartesian product of the chillers’ admissible sets, i.e., 
PLR S , S = F n
Since the admissible set F of a single chiller can be partitioned in four subsets, the overall admissible set S can be partitioned in 4 n subsets S j , j = 1 , , 4 n , each of which is in a one-to-one correspondence with the n digit multichiller code
s j = s j 1 s j n , s j i { 0 , 1 , 2 , 3 }
formed by the state codes σ i , i = 1 , , n , of the n chillers.
To make an example, consider the case of n = 3 chillers. Then, the possible 4 n = 64 subsets S j are associated to the multi-chiller codes as follows:
S 1 = F 0 × F 0 × F 0 s 1 = [ 0 , 0 , 0 ] S 2 = F 0 × F 0 × F 1 s 2 = [ 0 , 0 , 1 ] S 3 = F 0 × F 0 × F 2 s 3 = [ 0 , 0 , 2 ] S 4 = F 0 × F 0 × F 3 s 4 = [ 0 , 0 , 3 ] S 5 = F 0 × F 1 × F 0 s 5 = [ 0 , 1 , 0 ] S 6 = F 0 × F 1 × F 1 s 6 = [ 0 , 1 , 1 ] S 7 = F 0 × F 1 × F 2 s 7 = [ 0 , 1 , 2 ] S 8 = F 0 × F 1 × F 3 s 8 = [ 0 , 1 , 3 ] S 9 = F 0 × F 2 × F 0 s 9 = [ 0 , 2 , 0 ] = S 64 = F 3 × F 3 × F 3 s 64 = [ 3 , 3 , 3 ]
The n elements of the set S j are in a one-to-one correspondence with the chiller number i , 1 i n . Given the chiller number i, the notation S j [ i ] will denote the operating condition, either a point or a range, of the i-th chiller. For instance, within  S 8 , whose multi-chiller code is s 8 = [ 0 , 1 , 3 ] , the three chillers i = 1 , 2 , 3 will operate at the following conditions:    
S 8 [ 1 ] = 0 S 8 [ 2 ] = P L R m i n S 8 [ 3 ] = ( P L R m i n , 1 )

2.3. Divide and Conquer Strategy

For j = 1 , , 4 n , let QP ( j ) indicate the OCL problem (4) restricted to the subset S j .
Problem 1.
QP ( j )
min PLR S j i = 1 n P i ( P L R i )
subject   to i = 1 n P L R i · Q 100 % , i = Q l o a d
The 4 n problems QP ( j ) can be partitioned in two subsets, C and C ¯ :
-
C = { j s j i 3 , i } : all the partial load ratios are fixed so that S j has cardinality one;
-
C ¯ , when there is at least one chiller operating at intermediate part load ( i : s j i = 3 ).
Consider, for example, S 7 , whose multi-chiller code is s 7 = [ 0 , 1 , 2 ] : chiller #1 is switched off, chiller #2 is operating at minimum part load, and chiller #3 is operating at maximum part load. Within this subset, no optimization is actually needed because all the chillers’ P L R s are fixed and
S 7 = 0 P L R m i n 1 T
has cardinality one. Therefore, only a feasibility check is required: if constraint (8b) is satisfied, then S 7 is the optimal solution of QP ( 7 ) . Otherwise, QP ( 7 ) does not admit a solution. It is easy to see that the number of elements of the subset C is 3 n .
As a second example, consider S 8 , whose multichiller code is s 8 = [ 0 , 1 , 3 ] : here, the chiller #3 is operating at intermediate part load, i.e.
S 8 = PLR | PLR = 0 P L R m i n x 3 T
with x 3 F 3 . In this case, QP ( 8 ) is a constrained quadratic programming problem in the unknown x 3 .
As will be shown later, for  j C ¯ , the  QP ( j ) problems enjoy a remarkable property: their optimal solution, if it exists, is a critical point and no more than one critical point exists. In the following, PLR * ( j ) will denote:
  • the optimal solution of QP ( j ) , if  j C ;
  • the unique feasible critical point of QP ( j ) , if  j C ¯ .
Let A denote the set of integers j s.t. PLR * ( j ) exists. The associated value of the cost function will be denoted by
P t o t * ( j ) : = i = 1 n P i ( P L R i * ( j ) )
Then, the key idea is to exploit the partition of the solution space by a two-step procedure:
  • Solve QP ( j ) , obtaining PLR * ( j ) and P t o t * ( j ) for j = 1 , , 4 n ;
  • Letting
    j * = arg min j A P t o t * ( j )
    where obtain the globally optimal part load vector as PLR * = PLR * ( j * ) .
In the next subsection, it is shown how to reduce the inequality-constrained problems associated with j C ¯ to equality-constrained quadratic problems (EQP), for which a closed form solution is available.

2.4. Reduction to Equality-Constrained Problems

For a given j C ¯ , we denote by V j = { i { 1 , 2 , , n } s j i = 3 } the set of chillers operating at intermediate part loads, i.e., between P L R m i n and 1. We also let κ ( j ) denote the cardinality of V j . It is then possible to rewrite QP ( j ) as a reduced-order quadratic problem in the κ ( j ) n unknowns P L R i , i V j .
min P L R i , i V j i V j P i ( P L R i )
subject   to i V j P L R i · Q 100 % , i = Q ˜ l o a d
P L R m i n < P L R i < 1 , i V j
above, Q ˜ L o a d is
Q ˜ L o a d = Q L o a d i V j ¯ P L R i · Q 100 % , i
is the cooling load that must be supplied by the chillers operating at intermediate part loads.
Now, we associate to each QP ( j ) the corresponding equality-constrained quadratic problem EQP ( j ) , that is obtained by removing the inequality constraints (10c).
Problem 2.
EQP ( j )
min P L R i , i V j i V j P i ( P L R i )
subject   to i V j P L R i · Q 100 % , i = Q ˜ l o a d
It is convenient to rewrite EQP ( j ) in matrix form. For this purpose, we introduce a selection matrix M ( V j ) R κ ( j ) × n that selects κ ( j ) elements out of n.
Note that, being integers, the elements of V j admit an obvious ordering. Then,
M ( V j ) = [ M ] k , i = 1 if i = k - th element of V j 0 otherwise
In short, M ( V j ) will be denoted by M j . The reduced-order matrices are thus given by:
Q ˜ = M j Q M j T , Q ˜ R κ × κ E ˜ = E M j T , E ˜ R 1 × κ c ˜ = M j c , c ˜ R κ × 1 d ˜ = Q ˜ L o a d
Then, EQP ( j ) can be restated as:
x ˜ * ( j ) : = arg min x ˜ 1 2 x ˜ T Q ˜ x ˜ + c ˜ T x ˜ , x ˜ R κ × 1
s . t . E ˜ x ˜ = d ˜
By applying the first-order Karush–Kuhn–Tucker (KKT) necessary condition to the EQP ( j ) problem, the following linear system is obtained:
Q ˜ E ˜ T E ˜ 0 x ˜ * λ * = c ˜ d ˜
where x ˜ * identifies a critical point, either maximum, minimum or saddle, and  λ * R is the associated Lagrange multiplier. In order to guarantee the existence of the solution, a technical assumption is introduced.
Assumption 2.
E ˜ Q ˜ 1 E ˜ T 0 .
Note that Assumption 2 is immediately satisfied if q i > 0 , i , that is when the quadratic power consumption curves (3) are all convex, although this is not necessary.
Theorem 1.
Under Assumption 1, EQP ( j ) admits a unique critical point x * , given by
λ * = d ˜ + E ˜ Q ˜ 1 c ˜ E ˜ Q ˜ 1 E ˜ T
x ˜ * = Q ˜ 1 ( c ˜ + E ˜ T λ * )
Proof of Theorem 1.
In view of Assumption 1, det ( Q ) 0 , because  q i 0 , i . Moreover, Assumption 2 guarantees that E ˜ Q ˜ 1 E ˜ T 0 . Then, it is immediate to see that the KKT condition (15) admits (16)–(17) as a unique solution. □
Given the critical point x ˜ * , it is easy to obtain a critical point for QP ( j ) , as well. For this purpose, it is convenient to introduce an auxiliary vector x ¯ R n × 1 , such that:
x i ¯ = S j [ i ] , if i V j ¯ 0 , otherwise
The candidate critical point for QP ( j ) is thus given by:
PLR ^ * ( j ) = x ¯ + M j T x ˜ *
The keystone of the solution procedure is the connection between the critical points of QP ( j ) and EQP ( j ) , as stated in the following theorem.
Theorem 2.
The critical point PLR * ( j ) exists if PLR ^ * ( j ) S j . In such a case, PLR * ( j ) = PLR ^ * ( j ) .
Proof of Theorem 2.
Sufficiency. Assume that the critical point PLR ^ * for EQP(j) belongs to S j . Observe that EQP(j) has less contraints than QP(j). Therefore, if a critical point for EQP(j) is feasible for QP(j), it is ipso facto a critical point for QP(j), as well.
Necessity. Assume that QP(j) admits a critical point, say PLR * ( j ) . Given that EQP(j) and QP(j) differ only for strict inequality constraints, any critical point for QP(j) is critical also for EQP(j). Since EQP(j) admits at most one critical point, necessity is proven. □

2.5. Summary of the X-OCL Algorithm

We are now in a position to summarize the steps of the proposed algorithm, hereafter named X-OCL, for the exact solution of the OCL problem.
The partition of the solution set in the 4 n subsets S j , j = 1 , 4 n allows one to divide the MINLP problem into 4 n sub-problems QP ( j ) (8). For each subset S j , two situations can occur:
  • all the chillers part load ratios are fixed;
  • there is at least one chiller working at intermediate part load.
In the former case, only a feasibility check is required to decide whether the PLR configuration is to be kept as a candidate solution. In the latter case, the solution of the QP ( j ) problem is reduced to the solution of an EQP ( j ) problem that admits a unique critical point, easily computable in closed form. As stated in Theorem 2, if the critical point of EQP ( j ) belongs to S j , it coincides with the critical point of QP ( j ) . Otherwise, QP ( j ) does not admit a critical point.
Once all the QP ( j ) have been processed, the corresponding set of critical points { PLR * ( j ) , j = 1 , , 4 n } include the optimal solution of the overall OCL problem, which can be found just by comparing the associated system power consumptions P t o t * ( j ) . A pseudo-code summary is reported in Algorithm 1.
Algorithm 1: X-OCL.
  1:
Input: a , c , q R n × 1 , s R 4 n × n , Q l o a d R
  2:
Output: PLR * , P t o t *
  3:
  
  4:
for j = 1 , , 4 n do
  5:
   V j = { i { 1 , 2 , , n } : s j i = 3 } , κ = | V j |  
  6:
  if V j = then
  7:
                                                                                                                       { feasibility check}
  8:
   all the chillers part load ratios are fixed
  9:
    x S j
10:
   if i = 1 n x i · Q 100 % , i = Q l o a d then
11:
     PLR * ( j ) = x
12:
     P t o t * ( j ) = i = 1 n P i ( P L R i * ( j ) )
13:
   else
14:
     PLR * ( j ) =
15:
     P t o t * ( j ) =  
16:
   end if
17:
  else
18:
                                                                                   { solve the EQP ( j ) associated to QP ( j ) }
19:
    V j Q ˜ , E ˜ , c ˜ , d ˜  
20:
    λ * = d ˜ + E ˜ Q ˜ 1 c ˜ E ˜ Q ˜ 1 E ˜ T , λ * R  
21:
    x ˜ * = Q ˜ 1 ( c ˜ + E ˜ T λ * )  
22:
    x ¯ = S j [ i ] , if i V j ¯ 0 , otherwise  
23:
    PLR ^ * ( j ) = x ¯ + M j T x ˜ *  
24:
   if PLR ^ * ( j ) S j then
25:
     PLR * ( j ) = PLR ^ ( j )
26:
     P t o t * ( j ) = i = 1 n P i ( P L R i * ( j ) ) ,
27:
   else
28:
     PLR * ( j ) =
29:
     P t o t * ( j ) =  
30:
   end if
31:
  end if
32:
end for
33:
  
34:
j * = min j P t o t ( j )  
35:
PLR * = PLR * ( j * )
36:
P t o t * = P t o t ( j * )
Summarizing the difference between OCL and OCS is that the former one is a static problem, whereas the second one is intrinsically a dynamic problem. An OCS problem stripped of its dynamic constraints could be solved just as a sequence of OCL problems.

3. Test on Hsinchu Benchmark Model

3.1. Hsinchu Cooling Plant Model

The Hsinchu chiller system, originally described in [7], has become a widely used benchmark for the testing and comparison of OCL algorithms [6,15,17,18,19]. The case study involves six chillers installed in a semiconductor factory located in Hsinchu Scientific Garden (Taiwan) with a 7620 kW total cooling capacity. Quadratic models of the chillers’ energy consumption were obtained and validated from data collected every 5 min over a 5-month period [7]. The benchmark problem assumes that the condenser inlet water temperature is 24.5 ° C. The coefficients of the six chillers’ energy consumption models (3) are reported in [7] (Table 3) and the corresponding P-PLR curves are displayed in Figure 1. It is asked to solve the OCL problem for five different cooling loads, ranging from 70% to 90% of the system total cooling capacity. It is also required that the partial load ratio of each chiller never goes below 0.3. According to our notation, the following parameter settings are used:
-
a i , c i , and q i from [7];
-
Q L o a d = 90 % , 85 % , 80 % , 75 % , and 70 % of the chillers’ maximum capacity ( i = 1 n Q n o m );
-
P L R m i n , i = 0.3 , i .

3.2. OCL Benchmark: Results

The comparison between the X-OCL optimal solution and those obtained by the six literature methods in Table 1 sheds light on some errors in the published results.
For the three highest cooling load demands ( Q L o a d = { 6858 , 6477 , 6096 } [kW]), the X-OCL optimal solution coincides with the common solution provided by the six algorithms, thus confirming that they had reached the optimum. In the remaining two cases ( Q L o a d = 5717 and 5334 [kW]), the solution computed by X-OCL coincides with those of GAMS, EIWO and DCEDA, which, however, are apparently outperformed by IFA and DCSA. The worst performance is that of TLBO.
However, concerning the solution provided by IFA and DCSA, four wrong values of power consumption (identified by asterisks in Table 1) were published in [15,17]:
-
P t o t , I F A at load 75% and 70%;
-
P t o t , D C S A at load 75% and 70%.
In fact, Table 1 reports P t o t , I F A = P t o t , D C S A 3507.3 [kW] at 70% cooling load, which is inconsistent with the P L R i reported in the same table. Such inconsistency is easily verified by plugging the P L R i ’s into the chillers’ power consumption Equation (3) in order to obtain the individual chiller power consumptions of the P i [kW] column (the papers [15,17] describing algorithms IFA and DCSA do not report these individual consumptions). Apparently, the source of the error is the mechanical use of the chiller’s quadratic power consumption model out of its operational range, that is, in correspondence with the null partial load ratio. Obviously, the associated power consumption is null as well, but if the quadratic model has a negative constant term, as is the case for chillers 3–5, the formula will return a negative power consumption as if a turned-off chiller could generate free power.
In Table 2, the published chillers’ power consumption P p , i and the corresponding corrected values P c , i are reported for both IFA and DCSA. The values highlighted in red are the unfeasible chillers’ power consumption, of which the correction is reported in column P c , i .
Once the errors have been corrected, IFA and DCSA are no more optimal at 75 % and 70 % loads. At load 75 % P t o t , I F A rises from 3840.0690 to 3960.5709 [kW] ( + 3.14 % ), while at 70 % , it rises from 3507.2848 to 3627.7578 [kW] ( + 3.43 % ). Analogously, for DCSA, P t o t , D C S A goes from 3840.0545 to 3960.5580 [kW] ( + 3.14 % ) at 75 % and from 3507.2848 to 3627.7578 [kW] ( + 3.43 % ) at 70 % . The corrected results for all the six literature methods are reported in Table 3.
In conclusion, on the Hsinchu benchmark, GAMS, EIWO and DCEDA prove to be the best heuristic algorithms, as their solutions coincide with the optimal ones computed by X-OCL, for all considered loads.

4. The Optimal Chiller Sequencing Problem

The cooling load demand of a building can be subject to significant variations during the day. Consequently, solving the OCL problem in each time step t (for example, 20 min), just ignoring minimum up/down time constraints on the chillers, could lead to frequent switchings (chiller startups and shutdowns). In order to preserve chillers from excessive mechanical stress and increase their operating life, each machine should not be switched off before a minimal up-time is reached. Analogously, it should not be switched on too quickly. To comply with these requirements, minimum up/down-time constraints must be enforced in the formulation of the so-called optimal chiller sequencing (OCS) problem. In its full formulation, OCS is a dynamical problem, because, in order to minimize the cumulative power consumption, the current decision should also take into account future constraints. As a consequence, the solution approaches proposed in the literature range from dynamic programming to heuristic methods were designed to alleviate the complexity of the problem.

4.1. A Lower Bound to the OCS Problem

Any solution to the OCS problem must face some level of approximation. Even when dynamic programming is used, there is the necessity of forecasting future loads, which introduces a suboptimality margin with respect to the ideal solution based on perfect knowledge of the future load profile. The availability of an easy-to-compute limit of performance against which the results of heuristic methods can be benchmarked is therefore of interest. In order to derive such a limit, one can consider the relaxed OCS problem (R-OCS), that is, an OCS problem without up- and down-time constraints. The relaxed OCS problem boils down to a sequence of independent OCL problems to be solved at each step in correspondence with the associated load. The availability of an exact OCL solver, such as X-OCS, makes it possible to compute the exact solution of the relaxed OCS as well. Notably, in view of the independence of the OCL problems, what matters is not the load sequence but the load distribution, so that the R-OCS bound could be easily derived based on statistical distributions reflecting different production and weather scenarios.
Such a bound can be used to quantitatively assess the existing margin of improvement for a given heuristic OCS solver. In fact, if the achieved power consumption is close enough to the bound, there is no scope for the search of further improvements. Along this direction, in the following section, the X-OCL solver is used to derive a greedy OCS algorithm, whose performance is then assessed against the R-OCS bound.

4.2. X-OCS, a Greedy Approach to OCS

X-OCS is a greedy algorithm that reduces OCS to a sequence of OCL problems, solvable through X-OCL. The approach is greedy because at each time step, future constraints are ignored, and the optimal OCL solution, compatible with the current minimum up/down time constraint, is searched for. In the mathematical form, the greedy OCS problem can be written as a mixed-integer quadratic problem with linear constraints:
Problem 3.
greedy OCS
PLR * ( t ) : = arg min PLR i = 1 n P i ( P L R i ( t ) , T i ( t ) )
subject   to i = 1 n P L R i ( t ) · Q 100 % , i = Q l o a d ( t )
P L R i ( t ) = δ i ( t ) · x i ( t )
P L R m i n x i ( t ) 1
( T i O N ( t 1 ) M U T i ) ( δ i ( t ) δ i ( t 1 ) ) 0
( T i O F F ( t 1 ) M D T i ) ( δ i ( t ) δ i ( t 1 ) ) 0
δ i { 0 , 1 }
where M U T i and M D T i are the i-th chiller’s minimum up-time and minimum down-time limits, expressed in number of time steps. The time counters T i O N ( t ) and T i O F F ( t ) are expressed as:
T i O N ( t ) = ( 1 + T i O N ( t 1 ) ) δ i ( t )
T i O F F ( t ) = ( 1 + T i O F F ( t 1 ) ) ( 1 δ i ( t ) )
It is easy to observe that the greedy OCS problem is an OCL problem with the two additional constraints (19e–f) which force some chillers to be online/offline depending on their previous states δ ( τ ) , τ = t 1 , t 2 , .
At each time step t, two situations can occur:
-
all the chillers’ states δ i ( t ) are free, i.e., all the chillers have been online/offline for more time steps than those prescribed by MUT/MDT;
-
there is at least one chiller, say the i-th one, whose state δ i ( t ) is constrained to be online or offline ( δ i ( t ) = 1 or δ i ( t ) = 0 ) by the MUT or MDT.
Concerning the first case, the minimum up/down-time constraints are not active, therefore the step of the greedy OCS boils down to an OCL problem and its optimal solution can be found by the X-OCL algorithm. In the second case, instead, at each time t, the optimal solution is found by applying the X-OCL algorithm to a suitable subset of the OCL solution space S . At each time step t, the solution space of the greedy OCS is therefore given by:
P ( t ) = { j { 1 , , 4 n } s j , i 0 , i B O N ( t ) , s j , i = 0 , i B O F F ( t ) }
where
B O N ( t ) = { i { 1 , , n } ( T i O N ( t 1 ) M U T i ) < 0 } B O F F ( t ) = { i { 1 , , n } ( T i O F F ( t 1 ) M D T i ) < 0 }
are the sets of chillers which must be on and off, respectively.
The idea is to exploit the partitions of the solution space P j by the typical two-step procedure of X-OCL:
  • Solve QP ( j ) , obtaining PLR * ( j ) and P t o t * ( j ) for j A ( t ) P ( t ) ;
  • Letting
    j * = arg min j A ( t ) P ( t ) P t o t * ( j )
    obtain the globally optimal part load vector as PLR * = PLR * ( j * ) .
Herein, A ( t ) , which denotes that the set of integers j s.t. PLR * ( j ) exists, is a function of t because the load changes with time.

4.3. OCS Simulated Example

In this section, the performance of the X-OCS method is assessed on two case studies taken from the literature [23]. Case study 1 involves a hotel in Taipei with two 450 refrigeration tons (RT) chillers and two 1000 RT chillers, while case study 2, again from the Hsinchu Science Industrial District, features nine 1250 RT chillers. In both cases, chillers are described by their COP-PLR curves, expressed by a second-order polynomial model:
C O P i = α i + β i P L R i + γ i P L R i 2
where α i , β i and γ i are the chiller’s coefficients, reported in Table A1.
The first step was the identification of quadratic power consumption models (3) using data sampled from the COP-PLR curves. Further details regarding the identification procedure are reported in Appendix A. The coefficients of the quadratic P-PLR curves in Figure 2 are reported in Table A2.
The aim of both benchmarks is to compute the sequence of chillers’ partializations over one day, assuming 20-min stages, so as to minimize the cumulative power consumption, while satisfying the cooling demand constraint at each stage. The load demand profiles are displayed in Figure 3.
The parameters’ settings were as follows:
-
P L R m i n = 0.5 ;
-
M U T i = 3 i ;
-
M D T i = 1 i .
Recall that M U T i = 3 means that the i-th chiller must be on at least 3 consecutive steps before being turned off. Likewise, M D T i = 1 indicates that the i-th chiller, once turned off, must remain off at least 1 step.
The results obtained via X-OCS were compared with those obtained by dynamic programming, as reported in [23]. Where optimization was carried out under the ideal condition that all future loads are known in advance.
For sake of comparability with DP, the chillers power consumption associated with X-OCS was evaluated by plugging the PLRs computed by X-OCS into the original benchmark’s COP model [23] and not via the approximate quadratic power consumption curves used by X-OCS.

4.4. OCS Benchmark: Results

The results are shown in Table 4 and Table 5. In case study 1, the power consumption obtained by the X-OCS method was, at each stage, lower than or equal to that obtained by the DP method. From stage 1 to stage 29, the DP and X-OCS methods selected the same chillers. Since the MUT and MDT constraints were not active, the solutions coincided with the R-OCS ones. At stage 30, the R-OCS solution had chillers 4 and 3 switched on. However, the MUT constraint forced the DP and X-OCS to leave chillers 1 and 2 switched on. The same goes for stage 31. At stage 32, 44 and 59 the X-OCS method performed better than the DP one. Although the chillers were not constrained by minimum up/down time limits, apparently the DP method, as implemented in [23], could not find the global minimum. Notably, the greedy X-OCS method achieves a total electric power consumption (46,495.91 [kW]) which is less than 0.1 % greater than the minimum achievable bound R-OCS (64,432.56 [kW]), meaning that there is no practical margin of improvement left).
In the second case study, the cooling load demand varied slowly over time, so that the MUT and MDT constraints were never active. Therefore, the X-OCS, notwithstanding its greedy nature, attains the best achievable performance bound R-OCS. For the majority of the cooling loads, DP and X-OCS gave the same results, the only exceptions being at stages 22, 48–62, 68–72, where X-OCS performed slightly better. For both the case studies, the cumulative daily power consumptions obtained by DP (marked by asterisks) had been reported incorrectly in [23]:
-
Case study 1: P d a y = 645,220.08 [kW] instead of 64,883.84 [kW]
-
Case study 2: P d a y = 298,425.69 [kW] instead of 289,525.25 [kW]
As a matter of fact, the daily power consumption values reported in the paper did not match with the sums of the power consumptions at each step, which we used for the correction.
So far, the performances of alternative methods have been compared on OCL and OCS benchmark problems whose quadratic consumption models were taken from the literature. Moreover, a limited number of loads were considered.
In this section, the feasibility of HVAC efficient management based on the exact solution of the OCL problem is validated against a real-world scenario that includes the estimation of the chillers’ consumption models from the field data.

4.5. Field Data

The study was based on experimental data from a large HVAC system installed in a semiconductor fab located in Austria. The chiller unit is composed by five water-cooled centrifugal chillers, see Table 6, connected in parallel, in a constant primary flow chilled water system. Quarter-hourly data were recorded at different working conditions over a period of almost two years, from February 2017 to January 2019. Collected data include: temperatures, PLR, power consumption.
The time series of the cooling load demand is shown in Figure 4.
The chillers were subject to the following operating constraints:
-
P L R m i n = 0.2 ;
-
M U T i = 4 (one hour), i ;
-
M D T i = 2 (half an hour), i .

4.6. Chiller Energy Consumption Models

The chillers’ power consumption models were estimated using the evaporator cooling capacity Q e v a p [kW] and the condenser inlet water temperature T [ ° C], as covariates, and the compressor power consumption data P [kW], as target.
For each chiller, the dataset was randomly partitioned in two subsets: 70% for training and 30% for testing, respectively. The parameters β p , i of the model (2) were estimated via least squares fitting of the training data, discarding data with P L R < P L R m i n = 0.2 . The values of the estimated parameters β p , i are reported in Table 6, together with the percentage coefficient of variation, defined as CV % = 100 × SE ( β p , i ) / | β p , i | . The cooling capacity Q n o m are the nominal values provided by the manufacturer.
In Figure 5, the fitted surfaces P ( P L R , T ) are displayed against the validation data for each of the five chillers. It is seen that the the quadratic model, in spite of a few outliers that are not uncommon when data are collected in industrial frameworks, predicts the consumptions at different operating conditions well, as also confirmed by the goodness-of-fit (GOF) plots in Figure 6.
The 3-D histogram of the covariates ( Q e v a p , T ) sketched in Figure 7 shows that T, i.e., the condenser inlet water temperature depicted in Figure 8, is mainly concentrated in a narrow range centered around the setpoint, namely 21.5 ° C. Following what is usually done in the literature benchmarks, one could neglect temperature variations around the set point and solve the OCL and OCS problems using the chillers consumption curves at 21.5 ° C, displayed in Figure 9. However, the inspection of Figure 5 shows that for chiller 4, the power consumption is significantly affected by the temperature, especially in summer. Therefore, differently from other literature studies, OCL and OCS solutions, we computed this based on the complete model P ( P L R , T ) .

4.7. Real HVAC System: Assessment of Potential Savings

The field data were used to perform a retrospective analysis of the efficiency of the HVAC system management. More precisely, starting from the historical decisions and the associated cumulative power consumption, two comparisons were performed for the 2-year OCS problem. First of all, the lower bound R-OCS on the best achievable consumption was computed in order to quantitatively assess the potential improvement margin. Being a lower bound, R-OCS may be overly optimistic, so that it is important to evaluate the performance that can be obtained in practice. This was done by running the X-OCS solver, whose energy consumption could then be compared with the (ideal) R-OCS bound and the historically recorded power consumption.
The cumulative energy consumption recorded during the 2-year monitoring was 1.758 × 10 6 [kWh]. This figure can be compared with the R-OCS lower bound, equal to 1.600 × 10 6 [kWh]. This means that the potential margin of improvement is not larger than 8.97 % .
When the X-OCS algorithm was applied, the cumulative energy power consumption was 1.601 × 10 6 [kWh]. As a matter of fact, for this HVAC system, the loss of performance due to the suboptimality of the greedy algorithm is definitely negligible (it is less than 0.1 % ). On average, the power saving is 0.157 × 10 6 / ( 2 × 365 × 24 ) = 8.96 [kW]. In Figure 10, the difference between the consumption achieved by X-OCS and the lower bound is plotted on a weekly basis. The black and blue lines overlap almost perfectly, so that they cannot be distinguishable.

5. Execution Time

The X-OCL and X-OCS algorithms, coded in Matlab®, were executed on a standard laptop (Intel(R) i7-7500U dual-core with hyperthreading, RAM 16GB, 2.7 GHz). While no explicit parallelization of the algorithm was implemented, the solution of the 4 n QP problems was formulated as a unique algebraic computation using sparse matrices. This means that the algorithm may have benefited from some optimization automatically enforced by the Matlab® compiler. In order to compare the numerical efficiency, the Hsinchu benchmark was also solved using the CPLEX solver for costrained mixed integer problems as implemented in the GAMS environment, used by [6]. The results concerning the computational times are summarized in Table 7.
Concerning the Hsinchu OCL benchmark, the computational cost of X-OCL was negligible ( 0.11 s) and, more importantly, definitely lower than the time ( 1.991 s) spent by GAMS.
Concerning the three OCS benchmarks, in all cases, the X-OCS computation time was larger than that for R-OCS, because the greedy algorithm performs the additional task of looking for solutions that satisfy the up-/down-time constraints.
In the OCS examples, the larger number of loads makes it possible to assess the average computation time per load. For a small number of chillers, the average computation time per load of X-OCS is remarkably small: 0.368 / 72 = 5 × 10 3 s for OCS benchmark 1 (4 chillers) and 573.96 / 68,110 = 8.4 × 10 3 s for the OCS on field data (5 chillers). In the latter case, X-OCS took less than 6 min to solve the OCS problem over 2-year data with quarter-hour sampling.
As expected, in view of the exponential growth of the number of QP problems, the maximal average computation time per load is found in correspondence with the OCS benchmark 2, where 9 chillers are present. Nevertheless, the average time per load amounts to 228.29 / 72 = 3.17 s, which is totally affordable.
Finally, a simulated benchmark was set up in order to better assess the average computation time per load as a function of the number of chillers. In particular, nine OCL simulated experiments, characterized by a variable number of chillers, ranging from two to ten, were considered. The ten-chiller system was composed by the six Hsinchu chillers, whose models are reported in [7] (Table 3), with the first four replicated. For each experiment, the OCL problem was solved for five different cooling loads ranging from 70% to 90% of the chiller’s total cooling capacity. In Figure 11, the average computation time per load is displayed as a function of the number of chillers. The plot highlights the exponential growth but, at the same time, it shows that, even for a large HVAC system made of 10 chillers, the computation time per load (10 s) is all but prohibitive.

6. Discussion

The main purpose of the paper was the derivation of an exact algorithm for the solution to the OCL problem. This goal was successfully completed by a decomposition approach that exploits a suitable partition of the solution space. The availability of an exact method has been immediately exploited along two directions. First, it has become possible to give a definitive assessment on the performances of different literature methods that had been applied to some consolidated benchmark problems. By the way, the comparison with the exact solution revealed that some unrealistic performances had been declared in the literature, due to erroneous extrapolations of the power consumption curve. The exact method has been exploited also in relation with the optimal chiller sequencing problem. For a given cooling demand profile, if the dynamic up/down-time constraints are neglected, a sequence of OCL problems can be exactly solved to yield a lower bound limit, called R-OCS, to the performance of any feasible solution.
Concerning the practical applicability of the exact solution in place of heuristic approaches discuss in the literature, two main objections may be raised: (i) the explosion of the computational cost with the number of chillers, (ii) the need to make overly restrictive assumptions on the shape of the power consumption curves.
On the first side, it is indeed true that the OCL problem as formulated in (4) is NP-hard, as confirmed by Figure 11, where the exponential growth of the computations time is apparent. At the same time, the figure shows that, even for a medium/large-sized chiller system, e.g., 6–10 chillers, the computation time per load is less than 2 s. This suggests that, even for chiller systems used in large semiconductor factories, the cost of the exact solution is all but prohibitive. The numerical efficiency is a direct consequence of the tiny number of computations required to solve the elementary EQP problems, see (12). Moreover, the partition strategy underlying X-OCL implies that it is totally parallelizable into 4 n threads, a feature that has not been explicitly exploited and that could further speed up the solution.
The second objection to the general applicability of the exact solution has to do with the quadratic assumption made on the power consumption curve. Even if the majority of benchmarks satisfy this assumption, there is no stringent reason to rule out other functional descriptions. Nevertheless, when confronted with real data (see the field data benchmark), we found that a quadratic power consumption fitted well the recorded data, see Figure 6. Moreover, even when a single quadratic function is not adequate, it is still possible to switch to a piece-wise quadratic description. In that case, it would be immediate to generalize the exact algorithm by just increasing the number of partitions S j , in such a way that the problem still boils down to a set of easy-to-solve EQP problems.
In view of its numerical efficiency, the application of X-OCS to the solution of the OCS problem also appears very promising. Indeed, in the OCS benchmarks taken from the literature and in the field data OCS benchmark, the performance of X-OCS is very close to the lower bound, implying that there is no scope for the use of more sophisticated algorithms.

Author Contributions

Conceptualization, F.A., M.R., G.D.N.; methodology, F.A.; software, F.A.; validation, F.A., M.R., G.D.N.; formal analysis, F.A.; data curation, F.A.; writing—original draft preparation, F.A.; writing—review and editing, M.R., G.D.N.; project administration, G.D.N.; funding acquisition, G.D.N.; supervision, G.D.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially supported by the Italian Ministry for Research in the framework of the 2017 Program for Research Projects of National Interest (PRIN), Grant no. 2017YKXYXJ. A part of the work has been performed in the project Power Semiconductor and Electronics Manufacturing 4.0 (SemI40), under grant agreement No 692466. The project is co-funded by grants from Austria, Germany, Italy, France, Portugal and -Electronic Component Systems for European Leadership Joint Undertaking (ECSEL JU).

Conflicts of Interest

The authors declare no conflict of interest.

Acronyms

P i Electric power consumption of the i-th chiller [kW]
P t o t Multi-chiller system total power consumption [kW]
Q i Cooling power delivered by the i-th chiller [kW]
Q 100 % , i Maximum cooling power supplied under full capacity by the i-th chiller [kW]
Q l o a d Cooling load demand [kW]
P L R i Part load ratio of the i-th chiller
P L R m i n , i Minimum part load ratio of the i-th chiller
T i Condenser inlet water temperature on the i-th chiller [ ° C]
O C L Optimal Chiller Loading problem
O C S Optimal Chiller Sequencing problem
X - O C S Exact Optimal Chiller Sequencing algorithm
X - O C L Exact Optimal Chiller Loading algorithm
R - O C S Relaxed Optimal Chiller Sequencing problem
H V A C Heating, Ventilation and Air-Conditioning

Appendix A. Derivation of Power Consumption Models from the COP Ones

Assuming that the available COP-PLR curves represent the true model of chillers’ efficiency, quadratic approximated power consumption models can be easily derived as follows:
  • Sample N data points from the COP-PLR curve of the i-th chiller, uniformly in the range [ P L R m i n , i ,1] obtaining the pairs { P L R i ( k ) , C O P i ( k ) } , k = 1 , N ;
  • Compute the chiller’s power consumptions as
    P i ( k ) = P L R i ( k ) · Q 100 % , i C O P i ( k ) , k = 1 , N ;
  • Using the training set made of input-output paired samples { P L R i ( k ) , P i ( k ) } , k = 1 , N estimate the parameters of the quadratic model via ordinary least squares (OLS).
The above procedure was applied to the case studies presented in Section 4.3. For both the case studies, 50 data points were randomly sampled (uniformly in [0.5,1]) from the COP-PLR curves reported in Table A1 to obtain the training datasets represented in the left panels of Figure A1 and Figure A2 as red dots. The estimated P-kW curves are shown in the right panels of Figure A1 and Figure A2, and their parameters are reported in Table A2.
Table A1. Coefficients of chillers’ COP-PLR curves.
Table A1. Coefficients of chillers’ COP-PLR curves.
SystemsChiller α i β i γ i Q nom
Case 110.15613.9023−2.5909450
20.90001.8432−1.4188450
30.29323.0419−2.00541000
40.14153.6376−2.24691000
Case 211.56521.8094−0.98031250
21.05194.1471−2.41731250
30.57033.1602−2.09121250
40.32572.3513−1.42651250
50.54381.8668−1.23611250
61.52711.0634−0.72381250
70.78651.8473−1.16331250
80.84993.7768−2.28591250
91.11911.0228−0.75421250
Table A2. Estimated parameters of chillers’ P-PLR curves.
Table A2. Estimated parameters of chillers’ P-PLR curves.
SystemsChiller a i c i q i Q nom
Case 11243.58−398.01504.00450
2130.81−103.53309.65450
3417.51−444.57771.991000
4383.79−347.84611.641000
Case 2195.54321.92103.601250
2170.6841.43235.461250
3371.09−307.99693.761250
4477.85−217.09733.351250
5433.17−186–20810.891250
6104.21358.47205.501250
7272.33116.22457.961250
8218.68−20–94333.721250
9191.69276.56429.511250
Figure A1. (Left): COP-PLR curves of chillers in case study 1. Red dots are the 50 samples used to identify the quadratic approximate power consumption model. (Right): Identified quadratic P-PLR curves of chillers in case study 1. Red dots are the energy consumption data used for the training.
Figure A1. (Left): COP-PLR curves of chillers in case study 1. Red dots are the 50 samples used to identify the quadratic approximate power consumption model. (Right): Identified quadratic P-PLR curves of chillers in case study 1. Red dots are the energy consumption data used for the training.
Energies 13 06372 g0a1
Figure A2. (Left): COP-PLR curves of chillers in case study 2. Red dots are the efficiency sampled data. (Right): P-PLR curves of chillers in case study 2. Red dots are the energy consumption data used for the training.
Figure A2. (Left): COP-PLR curves of chillers in case study 2. Red dots are the efficiency sampled data. (Right): P-PLR curves of chillers in case study 2. Red dots are the energy consumption data used for the training.
Energies 13 06372 g0a2

References

  1. Pérez-Lombard, L.; Ortiz, J.; Pout, C. A review on buildings energy consumption information. Energy Build. 2008, 40, 394–398. [Google Scholar] [CrossRef]
  2. Tschudi, W.; Benschine, K.; Fok, S.; Rumsey, P. Cleanroom Energy Benchmarking in High-Tech and Biotech Industries; Technical Report; Lawrence Berkeley National Lab.: Berkeley, CA, USA, 2001. [Google Scholar]
  3. Chang, Y.C. A novel energy conservation method—optimal chiller loading. Electr. Power Syst. Res. 2004, 69, 221–226. [Google Scholar] [CrossRef]
  4. Chang, Y.C.; Chan, T.S.; Lee, W.S. Economic dispatch of chiller plant by gradient method for saving energy. Appl. Energy 2010, 87, 1096–1101. [Google Scholar] [CrossRef]
  5. Geem, Z.W. Solution quality improvement in chiller loading optimization. Appl. Therm. Eng. 2011, 31, 1848–1851. [Google Scholar] [CrossRef]
  6. Salari, E.; Askarzadeh, A. A new solution for loading optimization of multi-chiller systems by general algebraic modeling system. Appl. Therm. Eng. 2015, 84, 429–436. [Google Scholar] [CrossRef]
  7. Chang, Y.C. Genetic algorithm based optimal chiller loading for energy conservation. Appl. Therm. Eng. 2005, 25, 2800–2815. [Google Scholar] [CrossRef]
  8. Chang, Y.C.; Lin, F.A.; Lin, C.H. Optimal chiller sequencing by branch and bound method for saving energy. Energy Convers. Manag. 2005, 46, 2158–2172. [Google Scholar] [CrossRef]
  9. Chang, Y.C.; Chen, W.H.; Lee, C.Y.; Huang, C.N. Simulated annealing based optimal chiller loading for saving energy. Energy Convers. Manag. 2006, 47, 2044–2058. [Google Scholar] [CrossRef]
  10. Chang, Y.C. An innovative approach for demand side management—optimal chiller loading by simulated annealing. Energy 2006, 31, 1883–1896. [Google Scholar] [CrossRef]
  11. Ardakani, A.J.; Ardakani, F.F.; Hosseinian, S. A novel approach for optimal chiller loading using particle swarm optimization. Energy Build. 2008, 40, 2177–2187. [Google Scholar] [CrossRef]
  12. Lee, W.S.; Lin, L.C. Optimal chiller loading by particle swarm algorithm for reducing energy consumption. Appl. Therm. Eng. 2009, 29, 1730–1734. [Google Scholar] [CrossRef]
  13. Chang, Y.C.; Lee, C.Y.; Chen, C.R.; Chou, C.J.; Chen, W.H.; Chen, W.H. Evolution strategy based optimal chiller loading for saving energy. Energy Convers. Manag. 2009, 50, 132–139. [Google Scholar] [CrossRef]
  14. Lee, W.S.; Chen, Y.T.; Kao, Y. Optimal chiller loading by differential evolution algorithm for reducing energy consumption. Energy Build. 2011, 43, 599–604. [Google Scholar] [CrossRef]
  15. Dos Santos Coelho, L.; Klein, C.E.; Sabat, S.L.; Mariani, V.C. Optimal chiller loading for energy conservation using a new differential cuckoo search approach. Energy 2014, 75, 237–243. [Google Scholar] [CrossRef]
  16. Sulaiman, M.H.; Ibrahim, H.; Daniyal, H.; Mohamed, M.R. A new swarm intelligence approach for optimal chiller loading for energy conservation. Procedia-Soc. Behav. Sci. 2014, 129, 483–488. [Google Scholar] [CrossRef] [Green Version]
  17. dos Santos Coelho, L.; Mariani, V.C. Improved firefly algorithm approach applied to chiller loading for energy conservation. Energy Build. 2013, 59, 273–278. [Google Scholar] [CrossRef]
  18. Duan, P.; Li, J.; Wang, Y.; Sang, H.; Jia, B. Solving chiller loading optimization problems using an improved teaching-learning-based optimization algorithm. Optim. Control Appl. Methods 2018, 39, 65–77. [Google Scholar] [CrossRef]
  19. Zheng, Z.; Li, J. Optimal chiller loading by improved invasive weed optimization algorithm for reducing energy consumption. Energy Build. 2018, 161, 80–88. [Google Scholar] [CrossRef]
  20. Yu, J.; Liu, Q.; Zhao, A.; Qian, X.; Zhang, R. Optimal chiller loading in HVAC System Using a Novel Algorithm Based on the distributed framework. J. Build. Eng. 2020, 28, 101044. [Google Scholar] [CrossRef]
  21. Jabari, F.; Mohammadpourfard, M.; Mohammadi-ivatloo, B. Energy efficient hourly scheduling of multi-chiller systems using imperialistic competitive algorithm. Comput. Electr. Eng. 2020, 82, 106550. [Google Scholar] [CrossRef]
  22. Jabari, F.; Dibavar, A.A.; Mohammadi-Ivatloo, B. Economic Dispatch of Multiple-Chiller Plants Using Wild Goats Algorithm. In Optimization of Power System Problems; Springer: Berlin/Heidelberg, Germany, 2020; pp. 171–194. [Google Scholar]
  23. Chang, Y.C. An outstanding method for saving energy-optimal chiller operation. IEEE Trans. Energy Convers. 2006, 21, 527–532. [Google Scholar] [CrossRef]
  24. Huang, S.; Zuo, W.; Sohn, M.D. Amelioration of the cooling load based chiller sequencing control. Appl. Energy 2016, 168, 204–215. [Google Scholar] [CrossRef] [Green Version]
  25. Beghi, A.; Cecchinato, L.; Cosi, G.; Rampazzo, M. A PSO-based algorithm for optimal multiple chiller systems operation. Appl. Therm. Eng. 2012, 32, 31–40. [Google Scholar] [CrossRef]
Figure 1. Hsinchu benchmark: P-PLR curves at T = 24.5 ° C.
Figure 1. Hsinchu benchmark: P-PLR curves at T = 24.5 ° C.
Energies 13 06372 g001
Figure 2. Estimated chillers’ P-PLR curves on case study 1 (left) and case study 2 (right).
Figure 2. Estimated chillers’ P-PLR curves on case study 1 (left) and case study 2 (right).
Energies 13 06372 g002
Figure 3. Cooling load demand profile for case study 1 (left) and case study 2 (right).
Figure 3. Cooling load demand profile for case study 1 (left) and case study 2 (right).
Energies 13 06372 g003
Figure 4. Cooling load demand time series.
Figure 4. Cooling load demand time series.
Energies 13 06372 g004
Figure 5. Three-dimensional representation of power consumption models. Red dots represent the experimental data.
Figure 5. Three-dimensional representation of power consumption models. Red dots represent the experimental data.
Energies 13 06372 g005
Figure 6. Goodness-of-Fit (GOF) plots of the quadratic power consumption models. Blue dots are the validation data.
Figure 6. Goodness-of-Fit (GOF) plots of the quadratic power consumption models. Blue dots are the validation data.
Energies 13 06372 g006
Figure 7. 3-D histogram of the pairs ( Q l o a d , T ) recorded from February 2017 to January 2019.
Figure 7. 3-D histogram of the pairs ( Q l o a d , T ) recorded from February 2017 to January 2019.
Energies 13 06372 g007
Figure 8. Condenser inlet water temperature time series.
Figure 8. Condenser inlet water temperature time series.
Energies 13 06372 g008
Figure 9. P-PLR curves of Field data benchmark at T = 21.5 ° C.
Figure 9. P-PLR curves of Field data benchmark at T = 21.5 ° C.
Energies 13 06372 g009
Figure 10. Weekly power consumption time series: comparison bewteen actual recorded data (red), the lower bound R-OCL (black) and the consumption associated with X-OCS (blue).
Figure 10. Weekly power consumption time series: comparison bewteen actual recorded data (red), the lower bound R-OCL (black) and the consumption associated with X-OCS (blue).
Energies 13 06372 g010
Figure 11. Simulated X-OCL experiment: average computation time per load.
Figure 11. Simulated X-OCL experiment: average computation time per load.
Energies 13 06372 g011
Table 1. Comparison of the results of IFA, DCSA, GAMS, TLBO, EIWO and DCEDA on Hsinchu benchmark.
Table 1. Comparison of the results of IFA, DCSA, GAMS, TLBO, EIWO and DCEDA on Hsinchu benchmark.
Q load [kW]iIFA (2013)DCSA (2014)GAMS (2015)TLBO (2017)EIWO (2018)DCEDA (2020)
PLR i P i [kW] PLR i P i [kW] PLR i P i [kW] PLR i P i [kW] PLR i P i [kW] PLR i P i [kW]
6858 ( 90 % ) 10.812774-0.812726-0.8127808.97360.8186-0.8127-0.8126-
20.749527-0.749619-0.7496740.72750.7523-0.7492-0.7489-
31.000000-1.000000-1.0000903.34501.0000-1.0000-1.0000-
41.000000-1.000000-1.0000781.48901.0000-1.0000-1.0000-
51.000000-1.000000-1.0000755.20101.0000-1.0000-1.0000-
60.838603-0.838559-0.8386748.83920.8297-0.8390-0.8395-
4738.576 4738.575 4738.5753 4738.54 4738.575 4738.58
6477 ( 85 % ) 10.727803-0.727731-0.7277718.50400.727731-0.7275-0.7280-
20.656174-0.656132-0.6561641.19600.656132-0.6563-0.6564-
31.000000-1.000000-1.0000903.34501.000000-1.0000-1.0000-
41.000000-1.000000-1.0000781.48901.000000-1.0000-1.0000-
51.000000-1.000000-1.0000755.20101.000000-1.0000-1.0000-
60.716408-0.716524-0.7165621.91350.716524-0.7166-0.7160-
4421.649 4421.649 4421.6486 4421.65 4421.649 4421.65
6096 ( 80 % ) 10.642725-0.642735-0.6427639.14110.6431-0.6427-0.6431-
20.562642-0.562645-0.5626553.89550.5621-0.5628-0.5622-
31.000000-1.000000-1.0000903.34501.000000-1.0000-1.0000-
41.000000-1.000000-1.0000781.48901.000000-1.0000-1.0000-
51.000000-1.000000-1.0000755.20101.000000-1.0000-1.0000-
60.594504-0.594490-0.5945510.63470.5946-0.5944-0.5946-
4143.706 4143.706 4143.7064 4143.64 4143.706 4143.71
5717 ( 75 % ) 10.842218-0.843697-0.00000.00000.55765-0.0000-0.0000-
20.781365-0.783794-0.7150702.48090.46918-0.7151-0.7144-
30.000002-0.000001-1.0000903.34500.99995-1.0000-1.0000-
40.999995-1.000000-1.0000781.48901.00000-1.0000-1.0000-
51.000000-1.000000-1.0000755.20101.00000-1.0000-1.0000-
60.887053-0.883049-0.7934700.03730.47250-0.7933-0.7941-
3840.063   * 3840.055   * 3842.5532 3904.70 3842.553 3843.07
5334 ( 70 % ) 10.759350-0.749969-0.00000.00000.64179-0.0000-0.0000-
20.691121-0.682477-0.5835572.30740.66219-0.5834-0.5831-
30.000021-0.000012-1.0000903.34500.33009-1.0000-1.0000-
41.000000-1.000000-1.0000781.48900.99059-1.0000-1.0000-
51.000000-1.000000-1.0000755.20100.99900-1.0000-1.0000-
60.757897-0.776363-0.6217534.09500.58047-0.6218-0.6221-
3507.286   * 3507.270   * 3546.4375 3642.51 3546.438 3546.48
Table 2. Published chillers’ power consumptions and the corresponding corrected values.
Table 2. Published chillers’ power consumptions and the corresponding corrected values.
Q load [kW]iIFADCSA
PLR i [ ] P p , i [kW] P c , i [kW] PLR i [ ] P p , i [kW] P c , i [kW]
5717 ( 75 % ) 10.842218843.0026843.00260.843697844.7431844.7431
20.781365777.3214777.32140.783794780.1789780.1789
30.000002−120.50190.00000.000001−120.50350.0000
40.999995781.4855781.48551.000000781.4890781.4890
51.000000755.2010755.20101.000000755.2010755.2010
60.887053803.5604803.56040.883049798.9460798.9460
3840.0690 3960.5709 3840.0545 3960.5580
5334 ( 70 % ) 10.759350750.8680750.86800.749969741.1047741.1047
20.691121677.0164677.01640.682477668.0074668.0074
30.000021−120.47300.00000.000012−120.48670.0000
41.000000781.4890781.48901.000000781.4890781.4890
51.000000755.2010755.20101.000000755.2010755.2010
60.757897663.1834663.18340.776363682.1826682.1826
3507.2848 3627.7578 3507.4980 3627.9847
Table 3. Comparison of the results of IFA, DCSA, GAMS, TLBO, EIWO, DCEDA and X-OCL on Hsinchu benchmark.
Table 3. Comparison of the results of IFA, DCSA, GAMS, TLBO, EIWO, DCEDA and X-OCL on Hsinchu benchmark.
Q load [kW]iIFA (2013)DCSA (2014)GAMS (2015)TLBO (2017)EIWO (2018)DCEDA (2020)X-OCL
PLR i P i [kW] PLR i P i [kW] PLR i P i [kW] PLR i P i [kW] PLR i P i [kW] PLR i P i [kW] PLR i P i [kW]
6858 ( 90 % ) 10.812774809.05410.812726808.99990.8127808.97360.8186815.66540.8127808.97050.8126808.85950.8127809.0002
20.749527740.63130.749619740.73530.7496740.72750.7523743.77050.7492740.26180.7489739.92300.7496740.7352
31.000000903.34501.000000903.34501.0000903.34501.0000903.34501.0000903.34501.0000903.34501.0000903.3450
41.000000781.48901.000000781.48901.0000781.48901.0000781.48901.0000781.48901.0000781.48901.0000781.4890
51.000000755.20101.000000755.20101.0000755.20101.0000755.20101.0000755.20101.0000755.20101.0000755.2010
60.838603748.85200.838559748.80350.8386748.83920.8297739.06660.8390749.29030.8395749.84250.8386748.8030
4738.576 4738.575 4738.5753 4738.54 4738.575 4738.66 4738.5733
6477 ( 85 % ) 10.727803718.57420.727731718.50230.7277718.50400.727731718.50230.7275718.27150.7280718.77310.7277718.5010
20.656174641.23790.656132641.19590.6561641.19600.656132641.19590.6563641.36380.6564641.46370.6561641.1960
31.000000903.34501.000000903.34501.0000903.34501.000000903.34501.0000903.34501.0000903.34501.0000903.3450
41.000000781.48901.000000781.48901.0000781.48901.000000781.48901.0000781.48901.0000781.48901.0000781.4890
51.000000755.20101.000000755.20101.0000755.20101.000000755.20101.0000755.20101.0000755.20101.0000755.2010
60.716408621.80000.716524621.91320.7165621.91350.716524621.91320.7166621.98740.7160621.40200.7165621.9136
4421.649 4421.649 4421.6486 4421.65 4421.649 4421.6739 4421.6466
6096 ( 80 % ) 10.642725639.12690.642735639.13560.6427639.14110.6431639.45260.6427639.10520.6431639.45460.6427639.1359
20.562642553.89800.562645553.90060.5626553.89550.5621553.42760.5628554.03520.5622553.51440.5626553.9010
31.000000903.34501.000000903.34501.0000903.34501.000000903.34501.0000903.34501.0000903.34501.0000903.3450
41.000000781.48901.000000781.48901.0000781.48901.000000781.48901.0000781.48901.0000781.48901.0000781.4890
51.000000755.20101.000000755.20101.0000755.20101.000000755.20101.0000755.20101.0000755.20101.0000755.2010
60.594504510.64420.594490510.63240.5945510.63470.5946510.72570.5944510.55610.5946510.72570.5945510.6325
4143.706 4143.706 4143.7064 4143.64 4143.706 4143.7297 4143.3704
5717 ( 75 % ) 10.842218843.00260.843697844.74310.00000.00000.55765570.83540.00000.00000.00000.00000.00000.0000
20.781365777.32140.783794780.17800.7150702.48090.46918478.86570.7151702.55570.5831701.79870.7150702.4808
30.0000020.00000.0000010.00001.0000903.34500.99995903.31891.0000903.34501.0000903.34501.0000903.3450
40.999995781.48551.000000781.48901.0000781.48901.00000781.48901.0000781.48901.0000781.48901.0000781.4890
51.000000755.20101.000000755.20101.0000755.20101.00000755.20101.0000755.20101.0000755.20101.0000755.2010
60.887053803.56040.883049798.94600.7934700.03730.47250414.99160.7933699.92290.7941700.76820.7934700.0374
3960.579 3960.558 3842.5532 3904.70 3842.553 3842.602 3842.5532
5334 ( 70 % ) 10.759350750.86800.749969741.10470.00000.00000.64179638.31580.00000.00000.0000.00000.00000.0000
20.691121677.01640.682477668.00740.5835572.30740.66219647.27510.5834572.22360.5831571.95450.5835572.3068
30.0000210.00000.0000120.00001.0000903.34500.33009328.49621.0000903.34501.0000903.34501.0000903.3450
41.000000781.48901.000000781.48901.0000781.48900.99059774.87021.0000781.48901.0000781.48901.0000781.4890
51.000000781.20101.000000755.20101.0000755.20100.99900754.70261.0000755.20101.0000755.20101.0000755.2010
60.757897663.18340.776363682.18260.6217534.09500.58047489.84740.6218534.18030.6221534.44330.6217534.0957
3627.758 3627.98 3546.4375 3642.51 3546.438 3546.4327 3546.6437
Table 4. Case study 1.
Table 4. Case study 1.
Stage Load R-OCSDPX-OCS
kW Chiller kW Chiller kW Chiller [ PLR 1 , PLR 2 , PLR 3 , PLR 4 ]
1700441.134441.154441.134[0.00, 0.00, 0.00, 0.70]
20890556.534556.554556.534[0.00, 0.00, 0.00, 0.89]
211100695.324 2695.564 2695.324 2[0.00, 0.63, 0.00, 0.82]
221250801.664 2801.784 2801.664 2[0.00, 0.73, 0.00, 0.92]
231100695.324 2695.564 2695.324 2[0.00, 0.63, 0,00, 0.82]
241000634.424 2636.734 2634.424 2[0.00, 0.50, 0,00, 0.77]
28920579.964579.984579.964[0.00, 0.00, 0.00, 0.92]
2916101048.594 2 11048.434 2 11048.594 2 1[0.75, 0.75, 0.00, 0.94]
3016501081.634 31082.074 2 11082.074 2 1[0.76, 0.77, 0.00, 0.96]
3116701096.434 31099.694 2 11099.434 2 1[0.77, 0.78, 0.00, 0.97]
3217001118.964 31122.124 2 31118.964 3[0.00, 0.00, 0.78, 0.91]
3619001246.744 3 21247.614 2 31246.744 3 2[0.00, 0.67, 0.74, 0.86]
4317501150.684 3 21151.604 2 31150.684 2 3[0.00, 0.61, 0.69, 0.79]
4415751020.824 2 11030.564 31020.824 2 1[0.74, 0.73, 0.00, 0.92]
4519001246.814 3 21250.324 3 11249.424 3 1[0.70, 0.00, 0.73, 0.85]
4617901175.084 3 21180.754 3 11179.564 3 1[0.67, 0.00, 0.69, 0.79]
4717251135.324 3 21142.764 3 11135.324 3 2[0.00, 0.50, 0.70, 0.80]
4817751165.894 3 21166.784 3 21165.894 3 2[0.00, 0.62, 0.69, 0.80]
5719001246.754 3 21247.614 3 21246.754 3 2[0.00, 0.67, 0.74, 0.86]
5821001385.034 3 2 11386.084 3 2 11385.034 3 2 1[0.68, 0.63, 0.70, 0.81]
5920001319.664 3 11326.204 3 2 11319.664 3 1[0.73, 0.00, 0.77, 0.90]
6021001385.034 3 2 11386.084 3 2 11385.034 3 2 1[0.68, 0.63, 0.70, 0.81]
6119001246.754 3 21247.614 3 21246.754 3 2[0.00, 0.67, 0.74, 0.86]
6717501150.684 3 21151.604 3 21150.684 3 2[0.00, 0.61, 0.69, 0.79]
681475948.044 2 1948.084 2 1948.044 2 1[0.71, 0.67, 0.00, 0.85]
691200763.594 2765.204 1765.454 1[0.72, 0.00, 0.00, 0.88]
701050664.474 2671.504 1670.464 1[0.65, 0.00, 0.00, 0.76]
711100695.324 2700.104 1695.324 2[0.00, 0.63, 0.00, 0.82]
72800495.804495.824539.234 2[0.00, 0.50, 0.00, 0.57]
P d a y 64,432.56 kW64,883.84 * kW64,495.91 kW
Table 5. Case study 3.
Table 5. Case study 3.
Stage Load X-OCLDPX-OCS
kW chiller kW Chiller kW Chiller [ PLR 1 , PLR 2 , PLR 3 , PLR 4 , PLR 5 , PLR 6 , PLR 7 , PLR 8 , PLR 9 ]
162102899.688 6 3 2 12899.698 6 3 2 12899.688 6 3 2 1[1.00 1.00 0.97 0.00 0.00 1.00 0.00 1.00 0.00]
2161502848.508 6 3 2 12852.348 6 3 2 12848.508 6 3 2 1[1.00 1.00 0.92 0.00 0.00 1.00 0.00 1.00 0.00]
2262802974.848 6 3 2 1 93039.868 6 3 2 1 72974.848 6 3 2 1 9[1.00 1.00 0.72 0.00 0.00 0.80 0.00 1.00 0.50]
2773553658.938 6 3 2 1 73672.698 6 3 2 1 73658.938 6 3 2 1 7[1.00 1.00 0.93 0.00 0.00 1.00 0.95 1.00 0.00]
2875203768.288 6 3 2 1 7 93847.818 6 3 2 1 7 93768.288 6 3 2 1 7 9[1.00 1.00 0.79 0.00 0.00 1.00 0.73 1.00 0.50]
3487204659.678 6 3 2 1 7 94659.658 6 3 2 1 7 94659.678 6 3 2 1 7 9[1.00 1.00 0.98 0.00 0.00 1.00 1.00 1.00 1.00]
3589104791.008 6 3 2 1 7 9 44891.138 6 3 2 1 7 9 44791.008 6 3 2 1 7 9 4[1.00 1.00 0.85 0.74 0.00 1.00 0.83 1.00 0.70]
3690904917.898 6 3 2 1 7 9 45015.078 6 3 2 1 7 9 44917.898 6 3 2 1 7 9 4[1.00 1.00 0.88 0.77 0.00 1.00 0.87 1.00 0.74]
3792455032.608 6 3 2 1 7 9 45125.268 6 3 2 1 7 9 45032.608 6 3 2 1 7 9 4[1.00 1.00 0.91 0.80 0.00 1.00 0.91 1.00 0.78]
4798205511.198 6 3 2 1 7 9 45552.258 6 3 2 1 7 9 45511.198 6 3 2 1 7 9 4[1.00 1.00 1.00 0.90 0.00 1.00 1.00 1.00 0.96]
4898705552.258 6 3 2 1 7 9 4 55582.178 6 3 2 1 7 9 45552.258 6 3 2 1 7 9 4 5[1.00 1.00 0.87 0.76 0.67 1.00 0.86 1.00 0.73]
6387804703.108 6 3 2 1 7 9 44804.118 6 3 2 1 7 9 44703.108 6 3 2 1 7 9 4[1.00 1.00 0.83 0.73 0.00 1.00 0.80 1.00 0.66]
6485554511.808 6 3 2 1 7 94543.818 6 3 2 1 7 94511.808 6 3 2 1 7 9[1.00 1.00 0.96 0.00 0.00 1.00 1.00 1.00 0.88]
6777403907.968 6 3 2 1 7 93985.588 6 3 2 1 7 93907.968 6 3 2 1 7 9[1.00 1.00 0.81 0.00 0.00 1.00 0.76 1.00 0.62]
6874553728.528 6 3 2 1 7 93745.208 6 3 2 1 73728.528 6 3 2 1 7 9[1.00 1.00 0.77 0.00 0.00 0.98 0.71 1.00 0.50]
7264253055.918 6 3 2 1 93115.538 6 3 2 1 73055.918 6 3 2 1 9[1.00 1.00 0.74 0.00 0.00 0.89 0.00 1.00 0.50]
P d a y 296,440.08 kW289,525.25 * kW295,605.88 kW
Table 6. Estimated parameters of the quadratic power consumption model (2) and corresponding  C V % .
Table 6. Estimated parameters of the quadratic power consumption model (2) and corresponding  C V % .
Chiller β 0 , i β 1 , i β 2 , i β 3 , i Q nom
( CV % )( CV % )( CV % )( CV % )[kW]
149.6367124.7681269.57693.19502700
(5.24)(2.50)(0.82)(3.21)
256.2047154.9900278.02112.29972700
(11.97)(7.98)(2.95)(8.32)
3−10.9883520.3125−45.1983−0.33092700
(22.49)(0.85)(6.98)(24.57)
4−159.0637112.298815.752414.14612700
(6.15)(19.80)(8.21)(0.74)
546.7748461.7221−8.3474−1.18352700
(27.74)(1.84)(79.87)(47.23)
Table 7. Execution times for OCL and OSC benchmarks.
Table 7. Execution times for OCL and OSC benchmarks.
N ° ChillersN ° LoadsX-OCL [s]R-OCS [s]X-OCS [s]GAMS [s]
OCL Hsinchu benchmark650.110--1.991
OCS benchmark 1472-0.3020.368-
OCS benchmark 2972-226.77228.29-
OCS Field data568,110-354.03573.96-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Acerbi, F.; Rampazzo, M.; De Nicolao, G. An Exact Algorithm for the Optimal Chiller Loading Problem and Its Application to the Optimal Chiller Sequencing Problem. Energies 2020, 13, 6372. https://doi.org/10.3390/en13236372

AMA Style

Acerbi F, Rampazzo M, De Nicolao G. An Exact Algorithm for the Optimal Chiller Loading Problem and Its Application to the Optimal Chiller Sequencing Problem. Energies. 2020; 13(23):6372. https://doi.org/10.3390/en13236372

Chicago/Turabian Style

Acerbi, Federica, Mirco Rampazzo, and Giuseppe De Nicolao. 2020. "An Exact Algorithm for the Optimal Chiller Loading Problem and Its Application to the Optimal Chiller Sequencing Problem" Energies 13, no. 23: 6372. https://doi.org/10.3390/en13236372

APA Style

Acerbi, F., Rampazzo, M., & De Nicolao, G. (2020). An Exact Algorithm for the Optimal Chiller Loading Problem and Its Application to the Optimal Chiller Sequencing Problem. Energies, 13(23), 6372. https://doi.org/10.3390/en13236372

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop