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
, denote the cooling power delivered by the
i-th chiller and
the associated power consumption. For a prescribed overall cooling load demand
, the goal of the OCL problem is finding the cooling powers
that the
n chillers have to deliver in order to minimize the system total energy consumption
:
For each chiller, PLR (part load ratio) denotes the cooling load fraction, given by
where
is the maximum power supplied under full capacity operation. The vector of all PLRs is denoted by
When the
i-th chiller is turned on, it should not operate under a minimum
, denoted by
. For the subsequent derivations, it is convenient to introduce a binary variable
that indicates the status of the
i-th chiller and a real-valued variable
, such that:
The power consumption of the i-th chiller is assumed to depend mainly on and the condenser inlet water temperature , that is .
The consumption surface is either obtained from laboratory experiments or field data collected during operation.
Note that the problem of finding the optimal part load ratios
can be formulated as a mixed-integer nonlinear program (MINLP):
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
, 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 of the i-th chiller obeys the following model:where are the model parameters and is a suitable function. Moreover, it is assumed that . When the OCL problem is solved in a given time slot, the condenser inlet water temperature
can be assumed to be known. Then, for a given
, in the interval
, the consumption surface is a quadratic function of
alone:
with
,
, and
.
In view of the quadratic nature of the cost function, the system total energy consumption can be expressed in matrix form as follows:
where
Moreover, the equality constraint (1b) can be rewritten as
where
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
for
is
where each subset
is associated with one of the following operating conditions: switched off (
), minimum part load (
), maximum part load (
), intermediate part load (
). In the sequel,
will denote the state of the
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.,
Since the admissible set
of a single chiller can be partitioned in four subsets, the overall admissible set
can be partitioned in
subsets
, each of which is in a one-to-one correspondence with the
digit multichiller code
formed by the state codes
, of the
n chillers.
To make an example, consider the case of
chillers. Then, the possible
subsets
are associated to the multi-chiller codes as follows:
The
n elements of the set
are in a one-to-one correspondence with the chiller number
. Given the chiller number
i, the notation
will denote the operating condition, either a point or a range, of the
i-th chiller. For instance, within
, whose multi-chiller code is
, the three chillers
will operate at the following conditions:
2.3. Divide and Conquer Strategy
For , let indicate the OCL problem (4) restricted to the subset .
The problems can be partitioned in two subsets, and :
- -
: all the partial load ratios are fixed so that has cardinality one;
- -
, when there is at least one chiller operating at intermediate part load ().
Consider, for example,
, whose multi-chiller code is
: 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’
s are fixed and
has cardinality one. Therefore, only a feasibility check is required: if constraint (8b) is satisfied, then
is the optimal solution of
. Otherwise,
does not admit a solution. It is easy to see that the number of elements of the subset
is
.
As a second example, consider
, whose multichiller code is
: here, the chiller #3 is operating at intermediate part load, i.e.
with
. In this case,
is a constrained quadratic programming problem in the unknown
.
As will be shown later, for , the 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, will denote:
the optimal solution of , if ;
the unique feasible critical point of , if .
Let
denote the set of integers
j s.t.
exists. The associated value of the cost function will be denoted by
Then, the key idea is to exploit the partition of the solution space by a two-step procedure:
In the next subsection, it is shown how to reduce the inequality-constrained problems associated with to equality-constrained quadratic problems (EQP), for which a closed form solution is available.
2.4. Reduction to Equality-Constrained Problems
For a given
, we denote by
the set of chillers operating at intermediate part loads, i.e., between
and 1. We also let
denote the cardinality of
. It is then possible to rewrite
as a reduced-order quadratic problem in the
unknowns
.
above,
is
is the cooling load that must be supplied by the chillers operating at intermediate part loads.
Now, we associate to each the corresponding equality-constrained quadratic problem , that is obtained by removing the inequality constraints (10c).
It is convenient to rewrite in matrix form. For this purpose, we introduce a selection matrix that selects elements out of n.
Note that, being integers, the elements of
admit an obvious ordering. Then,
In short,
will be denoted by
. The reduced-order matrices are thus given by:
Then,
can be restated as:
By applying the first-order Karush–Kuhn–Tucker (KKT) necessary condition to the
EQP problem, the following linear system is obtained:
where
identifies a critical point, either maximum, minimum or saddle, and
is the associated Lagrange multiplier. In order to guarantee the existence of the solution, a technical assumption is introduced.
Assumption 2. .
Note that Assumption 2 is immediately satisfied if
, that is when the quadratic power consumption curves (
3) are all convex, although this is not necessary.
Theorem 1. Under Assumption 1, admits a unique critical point , given by Proof of Theorem 1. In view of Assumption 1, , because . Moreover, Assumption 2 guarantees that . Then, it is immediate to see that the KKT condition (15) admits (16)–(17) as a unique solution. □
Given the critical point
, it is easy to obtain a critical point for
, as well. For this purpose, it is convenient to introduce an auxiliary vector
, such that:
The candidate critical point for
is thus given by:
The keystone of the solution procedure is the connection between the critical points of and , as stated in the following theorem.
Theorem 2. The critical point exists if . In such a case, .
Proof of Theorem 2. Sufficiency. Assume that the critical point for EQP(j) belongs to . 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 . 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 subsets allows one to divide the MINLP problem into sub-problems (8). For each subset , two situations can occur:
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 problem is reduced to the solution of an problem that admits a unique critical point, easily computable in closed form. As stated in Theorem 2, if the critical point of belongs to , it coincides with the critical point of . Otherwise, does not admit a critical point.
Once all the
have been processed, the corresponding set of critical points
include the optimal solution of the overall OCL problem, which can be found just by comparing the associated system power consumptions
. A pseudo-code summary is reported in Algorithm 1.
Algorithm 1: X-OCL. |
- 1:
Input:, , - 2:
Output: - 3:
- 4:
fordo - 5:
, - 6:
if then - 7:
feasibility check} - 8:
all the chillers part load ratios are fixed - 9:
- 10:
if then - 11:
- 12:
- 13:
else - 14:
- 15:
- 16:
end if - 17:
else - 18:
solve the EQP associated to QP} - 19:
- 20:
- 21:
- 22:
- 23:
- 24:
if then - 25:
- 26:
, - 27:
else - 28:
- 29:
- 30:
end if - 31:
end if - 32:
end for - 33:
- 34:
- 35:
- 36:
|
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:
- -
and
from [
7];
- -
, and of the chillers’ maximum capacity ();
- -
.
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 ( [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 ( 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]:
- -
at load 75% and 70%;
- -
at load 75% and 70%.
In fact,
Table 1 reports
[kW] at 70% cooling load, which is inconsistent with the
reported in the same table. Such inconsistency is easily verified by plugging the
’s into the chillers’ power consumption Equation (
3) in order to obtain the individual chiller power consumptions of the
[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
and the corresponding corrected values
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
.
Once the errors have been corrected, IFA and DCSA are no more optimal at
and
loads. At load
rises from
to
[kW] (
), while at
, it rises from
to
[kW] (
). Analogously, for DCSA,
goes from
to
[kW] (
) at
and from
to
[kW] (
) at
. 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. where
and
are the
i-th chiller’s minimum up-time and minimum down-time limits, expressed in number of time steps. The time counters
and
are expressed as:
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 .
At each time step t, two situations can occur:
- -
all the chillers’ states 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 is constrained to be online or offline ( or ) 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
. At each time step
t, the solution space of the greedy OCS is therefore given by:
where
are the sets of chillers which must be on and off, respectively.
The idea is to exploit the partitions of the solution space by the typical two-step procedure of X-OCL:
Herein, , which denotes that the set of integers j s.t. 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:
where
,
and
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:
- -
;
- -
;
- -
.
Recall that means that the i-th chiller must be on at least 3 consecutive steps before being turned off. Likewise, 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: 645,220.08 [kW] instead of 64,883.84 [kW]
- -
Case study 2: 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:
- -
;
- -
(one hour), ;
- -
(half an hour), .
4.6. Chiller Energy Consumption Models
The chillers’ power consumption models were estimated using the evaporator cooling capacity [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
of the model (
2) were estimated via least squares fitting of the training data, discarding data with
. The values of the estimated parameters
are reported in
Table 6, together with the percentage coefficient of variation, defined as
. The cooling capacity
are the nominal values provided by the manufacturer.
In
Figure 5, the fitted surfaces
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
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
.
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 [kWh]. This figure can be compared with the R-OCS lower bound, equal to [kWh]. This means that the potential margin of improvement is not larger than .
When the X-OCS algorithm was applied, the cumulative energy power consumption was
[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
). On average, the power saving is
[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
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 ( s) and, more importantly, definitely lower than the time ( 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: s for OCS benchmark 1 (4 chillers) and 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 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
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
, 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
| Electric power consumption of the i-th chiller [kW] |
| Multi-chiller system total power consumption [kW] |
| Cooling power delivered by the i-th chiller [kW] |
| Maximum cooling power supplied under full capacity by the i-th chiller [kW] |
| Cooling load demand [kW] |
| Part load ratio of the i-th chiller |
| Minimum part load ratio of the i-th chiller |
| Condenser inlet water temperature on the i-th chiller [C] |
| Optimal Chiller Loading problem |
| Optimal Chiller Sequencing problem |
| Exact Optimal Chiller Sequencing algorithm |
| Exact Optimal Chiller Loading algorithm |
| Relaxed Optimal Chiller Sequencing problem |
| 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 [,1] obtaining the pairs ;
Compute the chiller’s power consumptions as
Using the training set made of input-output paired samples 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.
Systems | Chiller | | | | |
---|
Case 1 | 1 | 0.1561 | 3.9023 | −2.5909 | 450 |
2 | 0.9000 | 1.8432 | −1.4188 | 450 |
3 | 0.2932 | 3.0419 | −2.0054 | 1000 |
4 | 0.1415 | 3.6376 | −2.2469 | 1000 |
Case 2 | 1 | 1.5652 | 1.8094 | −0.9803 | 1250 |
2 | 1.0519 | 4.1471 | −2.4173 | 1250 |
3 | 0.5703 | 3.1602 | −2.0912 | 1250 |
4 | 0.3257 | 2.3513 | −1.4265 | 1250 |
5 | 0.5438 | 1.8668 | −1.2361 | 1250 |
6 | 1.5271 | 1.0634 | −0.7238 | 1250 |
7 | 0.7865 | 1.8473 | −1.1633 | 1250 |
8 | 0.8499 | 3.7768 | −2.2859 | 1250 |
9 | 1.1191 | 1.0228 | −0.7542 | 1250 |
Table A2.
Estimated parameters of chillers’ P-PLR curves.
Table A2.
Estimated parameters of chillers’ P-PLR curves.
Systems | Chiller | | | | |
---|
Case 1 | 1 | 243.58 | −398.01 | 504.00 | 450 |
2 | 130.81 | −103.53 | 309.65 | 450 |
3 | 417.51 | −444.57 | 771.99 | 1000 |
4 | 383.79 | −347.84 | 611.64 | 1000 |
Case 2 | 1 | 95.54 | 321.92 | 103.60 | 1250 |
2 | 170.68 | 41.43 | 235.46 | 1250 |
3 | 371.09 | −307.99 | 693.76 | 1250 |
4 | 477.85 | −217.09 | 733.35 | 1250 |
5 | 433.17 | −186–20 | 810.89 | 1250 |
6 | 104.21 | 358.47 | 205.50 | 1250 |
7 | 272.33 | 116.22 | 457.96 | 1250 |
8 | 218.68 | −20–94 | 333.72 | 1250 |
9 | 191.69 | 276.56 | 429.51 | 1250 |
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.
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.
References
- Pérez-Lombard, L.; Ortiz, J.; Pout, C. A review on buildings energy consumption information. Energy Build. 2008, 40, 394–398. [Google Scholar] [CrossRef]
- 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]
- Chang, Y.C. A novel energy conservation method—optimal chiller loading. Electr. Power Syst. Res. 2004, 69, 221–226. [Google Scholar] [CrossRef]
- 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]
- Geem, Z.W. Solution quality improvement in chiller loading optimization. Appl. Therm. Eng. 2011, 31, 1848–1851. [Google Scholar] [CrossRef]
- 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]
- Chang, Y.C. Genetic algorithm based optimal chiller loading for energy conservation. Appl. Therm. Eng. 2005, 25, 2800–2815. [Google Scholar] [CrossRef]
- 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]
- 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]
- Chang, Y.C. An innovative approach for demand side management—optimal chiller loading by simulated annealing. Energy 2006, 31, 1883–1896. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Chang, Y.C. An outstanding method for saving energy-optimal chiller operation. IEEE Trans. Energy Convers. 2006, 21, 527–532. [Google Scholar] [CrossRef]
- 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]
- 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 C.
Figure 1.
Hsinchu benchmark: P-PLR curves at C.
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).
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).
Figure 4.
Cooling load demand time series.
Figure 4.
Cooling load demand time series.
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.
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.
Figure 7.
3-D histogram of the pairs recorded from February 2017 to January 2019.
Figure 7.
3-D histogram of the pairs recorded from February 2017 to January 2019.
Figure 8.
Condenser inlet water temperature time series.
Figure 8.
Condenser inlet water temperature time series.
Figure 9.
P-PLR curves of Field data benchmark at 21.5 C.
Figure 9.
P-PLR curves of Field data benchmark at 21.5 C.
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).
Figure 11.
Simulated X-OCL experiment: average computation time per load.
Figure 11.
Simulated X-OCL experiment: average computation time per load.
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.
[kW] | i | IFA (2013) | DCSA (2014) | GAMS (2015) | TLBO (2017) | EIWO (2018) | DCEDA (2020) |
---|
| [kW] | | [kW] | | [kW] | | [kW] | | [kW] | | [kW] |
---|
| 1 | 0.812774 | - | 0.812726 | - | 0.8127 | 808.9736 | 0.8186 | - | 0.8127 | - | 0.8126 | - |
2 | 0.749527 | - | 0.749619 | - | 0.7496 | 740.7275 | 0.7523 | - | 0.7492 | - | 0.7489 | - |
3 | 1.000000 | - | 1.000000 | - | 1.0000 | 903.3450 | 1.0000 | - | 1.0000 | - | 1.0000 | - |
4 | 1.000000 | - | 1.000000 | - | 1.0000 | 781.4890 | 1.0000 | - | 1.0000 | - | 1.0000 | - |
5 | 1.000000 | - | 1.000000 | - | 1.0000 | 755.2010 | 1.0000 | - | 1.0000 | - | 1.0000 | - |
6 | 0.838603 | - | 0.838559 | - | 0.8386 | 748.8392 | 0.8297 | - | 0.8390 | - | 0.8395 | - |
∑ | | | | | | | | | | | | |
| 1 | 0.727803 | - | 0.727731 | - | 0.7277 | 718.5040 | 0.727731 | - | 0.7275 | - | 0.7280 | - |
2 | 0.656174 | - | 0.656132 | - | 0.6561 | 641.1960 | 0.656132 | - | 0.6563 | - | 0.6564 | - |
3 | 1.000000 | - | 1.000000 | - | 1.0000 | 903.3450 | 1.000000 | - | 1.0000 | - | 1.0000 | - |
4 | 1.000000 | - | 1.000000 | - | 1.0000 | 781.4890 | 1.000000 | - | 1.0000 | - | 1.0000 | - |
5 | 1.000000 | - | 1.000000 | - | 1.0000 | 755.2010 | 1.000000 | - | 1.0000 | - | 1.0000 | - |
6 | 0.716408 | - | 0.716524 | - | 0.7165 | 621.9135 | 0.716524 | - | 0.7166 | - | 0.7160 | - |
∑ | | | | | | | | | | | | |
| 1 | 0.642725 | - | 0.642735 | - | 0.6427 | 639.1411 | 0.6431 | - | 0.6427 | - | 0.6431 | - |
2 | 0.562642 | - | 0.562645 | - | 0.5626 | 553.8955 | 0.5621 | - | 0.5628 | - | 0.5622 | - |
3 | 1.000000 | - | 1.000000 | - | 1.0000 | 903.3450 | 1.000000 | - | 1.0000 | - | 1.0000 | - |
4 | 1.000000 | - | 1.000000 | - | 1.0000 | 781.4890 | 1.000000 | - | 1.0000 | - | 1.0000 | - |
5 | 1.000000 | - | 1.000000 | - | 1.0000 | 755.2010 | 1.000000 | - | 1.0000 | - | 1.0000 | - |
6 | 0.594504 | - | 0.594490 | - | 0.5945 | 510.6347 | 0.5946 | - | 0.5944 | - | 0.5946 | - |
∑ | | | | | | | | | | | | |
| 1 | 0.842218 | - | 0.843697 | - | 0.0000 | 0.0000 | 0.55765 | - | 0.0000 | - | 0.0000 | - |
2 | 0.781365 | - | 0.783794 | - | 0.7150 | 702.4809 | 0.46918 | - | 0.7151 | - | 0.7144 | - |
3 | 0.000002 | - | 0.000001 | - | 1.0000 | 903.3450 | 0.99995 | - | 1.0000 | - | 1.0000 | - |
4 | 0.999995 | - | 1.000000 | - | 1.0000 | 781.4890 | 1.00000 | - | 1.0000 | - | 1.0000 | - |
5 | 1.000000 | - | 1.000000 | - | 1.0000 | 755.2010 | 1.00000 | - | 1.0000 | - | 1.0000 | - |
6 | 0.887053 | - | 0.883049 | - | 0.7934 | 700.0373 | 0.47250 | - | 0.7933 | - | 0.7941 | - |
∑ | | | | | | | | | | | | |
| 1 | 0.759350 | - | 0.749969 | - | 0.0000 | 0.0000 | 0.64179 | - | 0.0000 | - | 0.0000 | - |
2 | 0.691121 | - | 0.682477 | - | 0.5835 | 572.3074 | 0.66219 | - | 0.5834 | - | 0.5831 | - |
3 | 0.000021 | - | 0.000012 | - | 1.0000 | 903.3450 | 0.33009 | - | 1.0000 | - | 1.0000 | - |
4 | 1.000000 | - | 1.000000 | - | 1.0000 | 781.4890 | 0.99059 | - | 1.0000 | - | 1.0000 | - |
5 | 1.000000 | - | 1.000000 | - | 1.0000 | 755.2010 | 0.99900 | - | 1.0000 | - | 1.0000 | - |
6 | 0.757897 | - | 0.776363 | - | 0.6217 | 534.0950 | 0.58047 | - | 0.6218 | - | 0.6221 | - |
∑ | | | | | | | | | | | | |
Table 2.
Published chillers’ power consumptions and the corresponding corrected values.
Table 2.
Published chillers’ power consumptions and the corresponding corrected values.
[kW] | i | IFA | DCSA |
---|
| [kW] | [kW] | | [kW] | [kW] |
---|
| 1 | 0.842218 | 843.0026 | 843.0026 | 0.843697 | 844.7431 | 844.7431 |
2 | 0.781365 | 777.3214 | 777.3214 | 0.783794 | 780.1789 | 780.1789 |
3 | 0.000002 | −120.5019 | 0.0000 | 0.000001 | −120.5035 | 0.0000 |
4 | 0.999995 | 781.4855 | 781.4855 | 1.000000 | 781.4890 | 781.4890 |
5 | 1.000000 | 755.2010 | 755.2010 | 1.000000 | 755.2010 | 755.2010 |
6 | 0.887053 | 803.5604 | 803.5604 | 0.883049 | 798.9460 | 798.9460 |
∑ | | | | | | |
| 1 | 0.759350 | 750.8680 | 750.8680 | 0.749969 | 741.1047 | 741.1047 |
2 | 0.691121 | 677.0164 | 677.0164 | 0.682477 | 668.0074 | 668.0074 |
3 | 0.000021 | −120.4730 | 0.0000 | 0.000012 | −120.4867 | 0.0000 |
4 | 1.000000 | 781.4890 | 781.4890 | 1.000000 | 781.4890 | 781.4890 |
5 | 1.000000 | 755.2010 | 755.2010 | 1.000000 | 755.2010 | 755.2010 |
6 | 0.757897 | 663.1834 | 663.1834 | 0.776363 | 682.1826 | 682.1826 |
∑ | | | | | | |
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.
[kW] | i | IFA (2013) | DCSA (2014) | GAMS (2015) | TLBO (2017) | EIWO (2018) | DCEDA (2020) | X-OCL |
---|
| [kW] | | [kW] | | [kW] | | [kW] | | [kW] | | [kW] | | [kW] |
---|
| 1 | 0.812774 | 809.0541 | 0.812726 | 808.9999 | 0.8127 | 808.9736 | 0.8186 | 815.6654 | 0.8127 | 808.9705 | 0.8126 | 808.8595 | 0.8127 | 809.0002 |
2 | 0.749527 | 740.6313 | 0.749619 | 740.7353 | 0.7496 | 740.7275 | 0.7523 | 743.7705 | 0.7492 | 740.2618 | 0.7489 | 739.9230 | 0.7496 | 740.7352 |
3 | 1.000000 | 903.3450 | 1.000000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 |
4 | 1.000000 | 781.4890 | 1.000000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 |
5 | 1.000000 | 755.2010 | 1.000000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 |
6 | 0.838603 | 748.8520 | 0.838559 | 748.8035 | 0.8386 | 748.8392 | 0.8297 | 739.0666 | 0.8390 | 749.2903 | 0.8395 | 749.8425 | 0.8386 | 748.8030 |
∑ | | 4738.576 | | 4738.575 | | 4738.5753 | | 4738.54 | | 4738.575 | | 4738.66 | | 4738.5733 |
| 1 | 0.727803 | 718.5742 | 0.727731 | 718.5023 | 0.7277 | 718.5040 | 0.727731 | 718.5023 | 0.7275 | 718.2715 | 0.7280 | 718.7731 | 0.7277 | 718.5010 |
2 | 0.656174 | 641.2379 | 0.656132 | 641.1959 | 0.6561 | 641.1960 | 0.656132 | 641.1959 | 0.6563 | 641.3638 | 0.6564 | 641.4637 | 0.6561 | 641.1960 |
3 | 1.000000 | 903.3450 | 1.000000 | 903.3450 | 1.0000 | 903.3450 | 1.000000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 |
4 | 1.000000 | 781.4890 | 1.000000 | 781.4890 | 1.0000 | 781.4890 | 1.000000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 |
5 | 1.000000 | 755.2010 | 1.000000 | 755.2010 | 1.0000 | 755.2010 | 1.000000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 |
6 | 0.716408 | 621.8000 | 0.716524 | 621.9132 | 0.7165 | 621.9135 | 0.716524 | 621.9132 | 0.7166 | 621.9874 | 0.7160 | 621.4020 | 0.7165 | 621.9136 |
∑ | | 4421.649 | | 4421.649 | | 4421.6486 | | 4421.65 | | 4421.649 | | 4421.6739 | | 4421.6466 |
| 1 | 0.642725 | 639.1269 | 0.642735 | 639.1356 | 0.6427 | 639.1411 | 0.6431 | 639.4526 | 0.6427 | 639.1052 | 0.6431 | 639.4546 | 0.6427 | 639.1359 |
2 | 0.562642 | 553.8980 | 0.562645 | 553.9006 | 0.5626 | 553.8955 | 0.5621 | 553.4276 | 0.5628 | 554.0352 | 0.5622 | 553.5144 | 0.5626 | 553.9010 |
3 | 1.000000 | 903.3450 | 1.000000 | 903.3450 | 1.0000 | 903.3450 | 1.000000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 |
4 | 1.000000 | 781.4890 | 1.000000 | 781.4890 | 1.0000 | 781.4890 | 1.000000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 |
5 | 1.000000 | 755.2010 | 1.000000 | 755.2010 | 1.0000 | 755.2010 | 1.000000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 |
6 | 0.594504 | 510.6442 | 0.594490 | 510.6324 | 0.5945 | 510.6347 | 0.5946 | 510.7257 | 0.5944 | 510.5561 | 0.5946 | 510.7257 | 0.5945 | 510.6325 |
∑ | | | | | | | | | | | | | | |
| 1 | 0.842218 | 843.0026 | 0.843697 | 844.7431 | 0.0000 | 0.0000 | 0.55765 | 570.8354 | 0.0000 | 0.0000 | 0.0000 | 0.0000 | 0.0000 | 0.0000 |
2 | 0.781365 | 777.3214 | 0.783794 | 780.1780 | 0.7150 | 702.4809 | 0.46918 | 478.8657 | 0.7151 | 702.5557 | 0.5831 | 701.7987 | 0.7150 | 702.4808 |
3 | 0.000002 | 0.0000 | 0.000001 | 0.0000 | 1.0000 | 903.3450 | 0.99995 | 903.3189 | 1.0000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 |
4 | 0.999995 | 781.4855 | 1.000000 | 781.4890 | 1.0000 | 781.4890 | 1.00000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 |
5 | 1.000000 | 755.2010 | 1.000000 | 755.2010 | 1.0000 | 755.2010 | 1.00000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 |
6 | 0.887053 | 803.5604 | 0.883049 | 798.9460 | 0.7934 | 700.0373 | 0.47250 | 414.9916 | 0.7933 | 699.9229 | 0.7941 | 700.7682 | 0.7934 | 700.0374 |
∑ | | | | | | | | | | | | | | |
| 1 | 0.759350 | 750.8680 | 0.749969 | 741.1047 | 0.0000 | 0.0000 | 0.64179 | 638.3158 | 0.0000 | 0.0000 | 0.000 | 0.0000 | 0.0000 | 0.0000 |
2 | 0.691121 | 677.0164 | 0.682477 | 668.0074 | 0.5835 | 572.3074 | 0.66219 | 647.2751 | 0.5834 | 572.2236 | 0.5831 | 571.9545 | 0.5835 | 572.3068 |
3 | 0.000021 | 0.0000 | 0.000012 | 0.0000 | 1.0000 | 903.3450 | 0.33009 | 328.4962 | 1.0000 | 903.3450 | 1.0000 | 903.3450 | 1.0000 | 903.3450 |
4 | 1.000000 | 781.4890 | 1.000000 | 781.4890 | 1.0000 | 781.4890 | 0.99059 | 774.8702 | 1.0000 | 781.4890 | 1.0000 | 781.4890 | 1.0000 | 781.4890 |
5 | 1.000000 | 781.2010 | 1.000000 | 755.2010 | 1.0000 | 755.2010 | 0.99900 | 754.7026 | 1.0000 | 755.2010 | 1.0000 | 755.2010 | 1.0000 | 755.2010 |
6 | 0.757897 | 663.1834 | 0.776363 | 682.1826 | 0.6217 | 534.0950 | 0.58047 | 489.8474 | 0.6218 | 534.1803 | 0.6221 | 534.4433 | 0.6217 | 534.0957 |
∑ | | | | | | | | | | | | | | |
Table 4.
Case study 1.
| | R-OCS | DP | X-OCS |
---|
| | | | | | |
---|
1 | 700 | 441.13 | 4 | 441.15 | 4 | 441.13 | 4 | [0.00, 0.00, 0.00, 0.70] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
20 | 890 | 556.53 | 4 | 556.55 | 4 | 556.53 | 4 | [0.00, 0.00, 0.00, 0.89] |
21 | 1100 | 695.32 | 4 2 | 695.56 | 4 2 | 695.32 | 4 2 | [0.00, 0.63, 0.00, 0.82] |
22 | 1250 | 801.66 | 4 2 | 801.78 | 4 2 | 801.66 | 4 2 | [0.00, 0.73, 0.00, 0.92] |
23 | 1100 | 695.32 | 4 2 | 695.56 | 4 2 | 695.32 | 4 2 | [0.00, 0.63, 0,00, 0.82] |
24 | 1000 | 634.42 | 4 2 | 636.73 | 4 2 | 634.42 | 4 2 | [0.00, 0.50, 0,00, 0.77] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
28 | 920 | 579.96 | 4 | 579.98 | 4 | 579.96 | 4 | [0.00, 0.00, 0.00, 0.92] |
29 | 1610 | 1048.59 | 4 2 1 | 1048.43 | 4 2 1 | 1048.59 | 4 2 1 | [0.75, 0.75, 0.00, 0.94] |
30 | 1650 | 1081.63 | 4 3 | 1082.07 | 4 2 1 | 1082.07 | 4 2 1 | [0.76, 0.77, 0.00, 0.96] |
31 | 1670 | 1096.43 | 4 3 | 1099.69 | 4 2 1 | 1099.43 | 4 2 1 | [0.77, 0.78, 0.00, 0.97] |
32 | 1700 | 1118.96 | 4 3 | 1122.12 | 4 2 3 | 1118.96 | 4 3 | [0.00, 0.00, 0.78, 0.91] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
36 | 1900 | 1246.74 | 4 3 2 | 1247.61 | 4 2 3 | 1246.74 | 4 3 2 | [0.00, 0.67, 0.74, 0.86] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
43 | 1750 | 1150.68 | 4 3 2 | 1151.60 | 4 2 3 | 1150.68 | 4 2 3 | [0.00, 0.61, 0.69, 0.79] |
44 | 1575 | 1020.82 | 4 2 1 | 1030.56 | 4 3 | 1020.82 | 4 2 1 | [0.74, 0.73, 0.00, 0.92] |
45 | 1900 | 1246.81 | 4 3 2 | 1250.32 | 4 3 1 | 1249.42 | 4 3 1 | [0.70, 0.00, 0.73, 0.85] |
46 | 1790 | 1175.08 | 4 3 2 | 1180.75 | 4 3 1 | 1179.56 | 4 3 1 | [0.67, 0.00, 0.69, 0.79] |
47 | 1725 | 1135.32 | 4 3 2 | 1142.76 | 4 3 1 | 1135.32 | 4 3 2 | [0.00, 0.50, 0.70, 0.80] |
48 | 1775 | 1165.89 | 4 3 2 | 1166.78 | 4 3 2 | 1165.89 | 4 3 2 | [0.00, 0.62, 0.69, 0.80] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
57 | 1900 | 1246.75 | 4 3 2 | 1247.61 | 4 3 2 | 1246.75 | 4 3 2 | [0.00, 0.67, 0.74, 0.86] |
58 | 2100 | 1385.03 | 4 3 2 1 | 1386.08 | 4 3 2 1 | 1385.03 | 4 3 2 1 | [0.68, 0.63, 0.70, 0.81] |
59 | 2000 | 1319.66 | 4 3 1 | 1326.20 | 4 3 2 1 | 1319.66 | 4 3 1 | [0.73, 0.00, 0.77, 0.90] |
60 | 2100 | 1385.03 | 4 3 2 1 | 1386.08 | 4 3 2 1 | 1385.03 | 4 3 2 1 | [0.68, 0.63, 0.70, 0.81] |
61 | 1900 | 1246.75 | 4 3 2 | 1247.61 | 4 3 2 | 1246.75 | 4 3 2 | [0.00, 0.67, 0.74, 0.86] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
67 | 1750 | 1150.68 | 4 3 2 | 1151.60 | 4 3 2 | 1150.68 | 4 3 2 | [0.00, 0.61, 0.69, 0.79] |
68 | 1475 | 948.04 | 4 2 1 | 948.08 | 4 2 1 | 948.04 | 4 2 1 | [0.71, 0.67, 0.00, 0.85] |
69 | 1200 | 763.59 | 4 2 | 765.20 | 4 1 | 765.45 | 4 1 | [0.72, 0.00, 0.00, 0.88] |
70 | 1050 | 664.47 | 4 2 | 671.50 | 4 1 | 670.46 | 4 1 | [0.65, 0.00, 0.00, 0.76] |
71 | 1100 | 695.32 | 4 2 | 700.10 | 4 1 | 695.32 | 4 2 | [0.00, 0.63, 0.00, 0.82] |
72 | 800 | 495.80 | 4 | 495.82 | 4 | 539.23 | 4 2 | [0.00, 0.50, 0.00, 0.57] |
| 64,432.56 kW | 64,883.84 *
kW | 64,495.91 kW |
Table 5.
Case study 3.
| | X-OCL | DP | X-OCS |
---|
| | | | | | |
---|
1 | 6210 | 2899.68 | 8 6 3 2 1 | 2899.69 | 8 6 3 2 1 | 2899.68 | 8 6 3 2 1 | [1.00 1.00 0.97 0.00 0.00 1.00 0.00 1.00 0.00] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
21 | 6150 | 2848.50 | 8 6 3 2 1 | 2852.34 | 8 6 3 2 1 | 2848.50 | 8 6 3 2 1 | [1.00 1.00 0.92 0.00 0.00 1.00 0.00 1.00 0.00] |
22 | 6280 | 2974.84 | 8 6 3 2 1 9 | 3039.86 | 8 6 3 2 1 7 | 2974.84 | 8 6 3 2 1 9 | [1.00 1.00 0.72 0.00 0.00 0.80 0.00 1.00 0.50] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
27 | 7355 | 3658.93 | 8 6 3 2 1 7 | 3672.69 | 8 6 3 2 1 7 | 3658.93 | 8 6 3 2 1 7 | [1.00 1.00 0.93 0.00 0.00 1.00 0.95 1.00 0.00] |
28 | 7520 | 3768.28 | 8 6 3 2 1 7 9 | 3847.81 | 8 6 3 2 1 7 9 | 3768.28 | 8 6 3 2 1 7 9 | [1.00 1.00 0.79 0.00 0.00 1.00 0.73 1.00 0.50] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
34 | 8720 | 4659.67 | 8 6 3 2 1 7 9 | 4659.65 | 8 6 3 2 1 7 9 | 4659.67 | 8 6 3 2 1 7 9 | [1.00 1.00 0.98 0.00 0.00 1.00 1.00 1.00 1.00] |
35 | 8910 | 4791.00 | 8 6 3 2 1 7 9 4 | 4891.13 | 8 6 3 2 1 7 9 4 | 4791.00 | 8 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] |
36 | 9090 | 4917.89 | 8 6 3 2 1 7 9 4 | 5015.07 | 8 6 3 2 1 7 9 4 | 4917.89 | 8 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] |
37 | 9245 | 5032.60 | 8 6 3 2 1 7 9 4 | 5125.26 | 8 6 3 2 1 7 9 4 | 5032.60 | 8 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] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
47 | 9820 | 5511.19 | 8 6 3 2 1 7 9 4 | 5552.25 | 8 6 3 2 1 7 9 4 | 5511.19 | 8 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] |
48 | 9870 | 5552.25 | 8 6 3 2 1 7 9 4 5 | 5582.17 | 8 6 3 2 1 7 9 4 | 5552.25 | 8 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] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
63 | 8780 | 4703.10 | 8 6 3 2 1 7 9 4 | 4804.11 | 8 6 3 2 1 7 9 4 | 4703.10 | 8 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] |
64 | 8555 | 4511.80 | 8 6 3 2 1 7 9 | 4543.81 | 8 6 3 2 1 7 9 | 4511.80 | 8 6 3 2 1 7 9 | [1.00 1.00 0.96 0.00 0.00 1.00 1.00 1.00 0.88] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
67 | 7740 | 3907.96 | 8 6 3 2 1 7 9 | 3985.58 | 8 6 3 2 1 7 9 | 3907.96 | 8 6 3 2 1 7 9 | [1.00 1.00 0.81 0.00 0.00 1.00 0.76 1.00 0.62] |
68 | 7455 | 3728.52 | 8 6 3 2 1 7 9 | 3745.20 | 8 6 3 2 1 7 | 3728.52 | 8 6 3 2 1 7 9 | [1.00 1.00 0.77 0.00 0.00 0.98 0.71 1.00 0.50] |
⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
72 | 6425 | 3055.91 | 8 6 3 2 1 9 | 3115.53 | 8 6 3 2 1 7 | 3055.91 | 8 6 3 2 1 9 | [1.00 1.00 0.74 0.00 0.00 0.89 0.00 1.00 0.50] |
| 296,440.08 kW | 289,525.25 * kW | 295,605.88 kW |
Table 6.
Estimated parameters of the quadratic power consumption model (
2) and corresponding
.
Table 6.
Estimated parameters of the quadratic power consumption model (
2) and corresponding
.
Chiller | | | | | |
---|
() | () | () | () | [kW] |
---|
1 | 49.6367 | 124.7681 | 269.5769 | 3.1950 | 2700 |
(5.24) | (2.50) | (0.82) | (3.21) |
2 | 56.2047 | 154.9900 | 278.0211 | 2.2997 | 2700 |
(11.97) | (7.98) | (2.95) | (8.32) |
3 | −10.9883 | 520.3125 | −45.1983 | −0.3309 | 2700 |
(22.49) | (0.85) | (6.98) | (24.57) |
4 | −159.0637 | 112.2988 | 15.7524 | 14.1461 | 2700 |
(6.15) | (19.80) | (8.21) | (0.74) |
5 | 46.7748 | 461.7221 | −8.3474 | −1.1835 | 2700 |
(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 Chillers | N Loads | X-OCL [s] | R-OCS [s] | X-OCS [s] | GAMS [s] |
---|
OCL Hsinchu benchmark | 6 | 5 | 0.110 | - | - | 1.991 |
OCS benchmark 1 | 4 | 72 | - | 0.302 | 0.368 | - |
OCS benchmark 2 | 9 | 72 | - | 226.77 | 228.29 | - |
OCS Field data | 5 | 68,110 | - | 354.03 | 573.96 | - |
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).