Ice Storage Air-conditioning System Simulation with Dynamic Electricity Pricing: a Demand Response Study

This paper presents an optimal dispatch model of an ice storage airconditioning system for participants to quickly and accurately perform energy saving and demand response, and to avoid the over contact with electricity price peak. The schedule planning for an ice storage airconditioning system of demand response is mainly to transfer energy consumption from the peak load to the partial-peak or off-peak load. Least Squares Regression (LSR) is used to obtain the polynomial function for the cooling capacity and the cost of power consumption with a real ice storage airconditioning system. Based on the dynamic electricity pricing, the requirements of cooling loads, and all technical constraints, the dispatch model of the ice-storage airconditioning system is formulated to minimize the operation cost. The Improved Ripple Bee Swarm Optimization (IRBSO) algorithm is proposed to solve the dispatch model of the ice storage airconditioning system in a daily schedule on summer. Simulation results indicate that reasonable solutions provide a practical and flexible framework allowing the demand response of ice storage airconditioning systems to demonstrate the optimization of its energy savings and operational efficiency and offering greater energy efficiency.


Introduction
In recent years, due to rapid industrial development, the electricity demand growth has been substantial.In the summer, the peak load continues to rise year by year, but international oil costs and production lack stability.Therefore, there is an urgent need to seek and develop alternative sources of energy, as well as to carry out comprehensive reviews of the efficiency of using these energy sources.Thus, improving the efficiency of the chiller units to reduce electricity costs is very important.Currently, public buildings, such as large hospitals, office buildings and shopping malls, use central air-conditioning systems.In the current design of large-scale central air-conditioning systems, the main structure includes fans, ice water pumps, chiller units, cooling water pumps and cooling fans.By using the hot and cold exchange system, the indoor temperature heat load is transferred to the outdoors and the cooling load is about 40% of the total power load.The system peak loads may not be met by the air-conditioning system, so the ice storage technology has begun to be used for reducing Energies 2016, 9, 113 2 of 16 energy consumption.An ice storage system can reduce the off-peak demand load so the power charges are at the off-peak rate.The central ice storage air-conditioning system [1,2], which integrates chillers and an ice storage tank, has been effectively applied to the management of demand response.
The ice storage system runs the refrigerating compressor with off-peak loads during the night.The water is frozen to store plenty of latent heat by phase change, and then it is discharged to reduce the load profile of the air conditioning on peak load or partial-peak load during the day [3].
The economical and operational benefits offered by the chiller dispatch under time-of-use (TOU) [4] in Taiwan.This not only reduces the on-peak load for demand response, but also lowers the electricity price as the preferential rates are shifted from the peak load profile to off-peak periods.The ice storage tank shifts the peak load to off-peak hours in order to reduce the problem of overloads during peak hours, thereby improving the efficiency of off-peak electricity use and balancing electricity consumption.Therefore, the demand response strategies in an ice storage air-conditioning system are very important for saving energy [5].The ice storage air-conditioning system is a combination of different chiller units and an ice storage tank, each with their own characteristic curves and different operating limitations.In related studies, the cubic equation was used to construct the characteristic curve to achieve higher accuracy; the accuracy of the efficiency curve of each chiller unit and the ice storage tank was improved, but the difficulty of the problem increased.
The problem is scheduling of the ice storage air-conditioning system for the industry, service industry, and large commercial office centers.Sun and Reddy [6] used a Taylor-Series expansion to solve the economic dispatch of chiller plants problem in order to find the optimal cooling load share for chiller units with different efficiencies.In reference [7] a particle swarm algorithm was used to facilitate optimization of ice-storage air-conditioning systems and to develop optimal operating strategies.In [3] the use of ice thermal energy storage to reduce and shift the electricity consumption of air-conditioning systems from on-peak hours to off-peak hours was proposed.In [8] ice storage systems were used to shift the peak cooling demand to off-peak periods in office buildings and several other studies have discussed the optimal design and control of cooling systems with ice storage [9][10][11].
A high-dimensional ice storage air-conditioning system presents a complex problem.In recent years, numerous experts and scholars have adopted the labor division and cooperation procedures of biological communities to develop many heuristic algorithms.Several studies have been carried out on the effectiveness of energy storage systems and the related literature includes the combination of sequential quadratic programming (SQP) [12,13], evolutionary programming (EP) [14], or the particle swarm algorithm (PSO) [15,16].Also, optimization methods for initial solutions have often been adopted [17,18].In particular, the particle swarm algorithm containing time-varying acceleration coefficients [19,20] was proposed for application in a variety of optimization problems.By using time-varying weight factors, the search space and accuracy of the algorithm have been improved.These hybrid and weight factor adjustment optimization methods are worth further study and discussion.
The enhanced bee swarm optimization algorithm (EBSO) [21] is a kind of optimization technique based on collective intelligence that originates from the foraging behavior of bees.Bee swarms communicate and collaborate with each other through pheromones to exhibit intelligent behavior.This behavior has been used in the development of several related algorithms [22] and applied to actual optimization problems.In this paper, the improved ripple bee swarm optimization (IRBSO) uses nonlinear ripple weight factors to improve the behavior pattern of each bee swarm and increase its search efficiency and accuracy in high dimensions.The cognitive and social parameter factors use nonlinear ripple weight to improve the search ability.In this paper, the IRBSO is proposed to solve the system's dynamic electricity price for demand response requires.In our results show that the proposed method is feasible, robust, and more effective than many previously developed algorithms.

System Structure
The main structure of a large-scale central ice storage air-conditioning system includes fans, ice water pumps, chillers, an ice storage tank, ice storage pump, cooling water pumps and cooling fans, as shown in Figure 1.To meet the cooling load and provide large buildings with a stable supply, multiple chiller units are operated in parallel and connected to a common supply system.By using the hot and cold exchange system, the indoor temperature heat load is transferred to the outdoors.This architecture can effectively make the chiller and ice storage operations more flexible and reliable, as well as provide a standby capacity, while requiring lower unit maintenance costs.The ice storage air-conditioning system of demand response works by operating the ice storage tank to take advantage of the off-peak electricity price.This architecture can effectively make the chiller and ice storage operations more flexible and reliable, as well as provide a standby capacity, while requiring lower unit maintenance costs.The ice storage air-conditioning system of demand response works by operating the ice storage tank to take advantage of the off-peak electricity price.The cold energy is stored in the form of sensible heat, and the ice is melted into water to release cold energy in order to provide the required cooling load.In this way, air-conditioning demands are met, electricity use during peak hours in the daytime is reduced, the stored ice is released to transfer the peak electricity demands and power costs are reduced.This paper is to derive the best single-day schedule planning for ice storage air-conditioning systems.The mathematical formulas of the chillers and ice storage tank are introduced as follows: . . .

Figure 1.
Ice storage air-conditioning system.

Chiller Capacity Cooling Load
The cooling load capacity is calculated based on the return water temperature and supply water temperature of the chillers, as well as the flow rate of the chilled water.The calculations of the cooling load capacity for the chillers are expressed in Equations ( 1) and ( 2 The cold energy is stored in the form of sensible heat, and the ice is melted into water to release cold energy in order to provide the required cooling load.In this way, air-conditioning demands are met, electricity use during peak hours in the daytime is reduced, the stored ice is released to transfer the peak electricity demands and power costs are reduced.This paper is to derive the best single-day schedule planning for ice storage air-conditioning systems.The mathematical formulas of the chillers and ice storage tank are introduced as follows:

Chiller Capacity Cooling Load
The cooling load capacity is calculated based on the return water temperature and supply water temperature of the chillers, as well as the flow rate of the chilled water.The calculations of the cooling load capacity for the chillers are expressed in Equations ( 1) and (2): Q chiller,i " `LPM chiller,i ˆp∆T chw,i q ˆρw ˆCpw ˘ˆ60 min (1) Energies 2016, 9, 113 4 of 16 ∆T chw,i : temperature difference of chilled water ( ˝C) T chwrt,i : the return temperature of chilled water ( ˝C) T chwst,i : the supply temperature of chilled water ( ˝C) LPM chiller,i : the flow rate of chilled water (L/min) ρ w : the density of chilled water (1 kg/m 3 ) C pw : the specific heat of chilled water (4.186 kJ/kg¨K) Q chiller,i : the cooling load of the ith chiller (kJ/h).
The power consumption of chillers is a convex function of the cooling load capacity and the each chiller represents the operational control in Equations ( 3) and (4): where P chiller,i is the power consumption (kW) of chiller i; and a i , b i , c i and d i are the regression coefficients of the operation curve of the ith chiller.

Ice Storage Cooling System
The ice storage tank is stored in the form of sensible heat, and the ice is melted into water to release cold energy in order to provide the required cooling load.The ice storage operation during off-peaks of electricity mainly involves freezing water into ice in order to store the cooling capability; the cooling load capability of the melting ice is calculated using the water temperature differences and the water flow valve.The cooling load capacity of an ice storage tank is calculated based on the return water temperature, supply water temperature and the flow rate of the chilled water.The calculation of the cooling load capacity is represented in Equations ( 5) and (6): ∆T isw " T iswr ´Tisws ∆T isw : the temperature difference of ice storage water ( ˝C) T iswr : the return temperature of ice storage water ( ˝C) T isws : the supply temperature of ice storage water ( ˝C) LPM: the flow rate of ice storage water (L/min) Q ice : cooling load capacity of ice storage tank (kJ/h).
The operation modes of the ice storage tank can be divided into a "charge process" and a "discharge process".The ice storage operation in the charging process mainly involves freezing water into ice in order to store the cooling capability, while the ice melting in the discharging process supplies the cooling capability.The ice storage power consumption in the charging process is briefly formulated in Equation (7) and the discharging process in Equation (8): where a cp , b cp , c cp and d cp are the regression coefficients of Q ice,cp and P ice,cp is the ice storage charge process power (kW).The cooling load capacity of the ice melting is calculated based on the amount of melting ice, where a dp , b dp , c dp and d dp are the regression coefficients of Q ice,dp and P ice,dp is the ice storage pump power (kW).

Ice Storage Air-Conditioning System
The ice storage air-conditioning system's dynamic electricity price for demand response requires minimizing the total operating cost.This paper also considers the connection with a utility company regarding the TOU rate [4].The total operating cost includes the charging and discharging process of the power from an ice storage tank, which is formulated as Equation ( 9): where h is the scheduling hour; j is the total number of chiller units; U t i is the ith chiller's operation status (on/off) at hour t; SU and SD are the start-on and shut-off power of the chiller units, respectively; price t and price ice are the dynamic electricity price at hour t and the favorable price for ice storage.The ice storage charge process time is from t = 1 to k 1 = 10 (22:00 p.m. to 07:00 a.m.) have favorable price and discharge processing times is from t = 11 to k 2 = 24 (8:00 a.m. to 21:00 p.m.) used summer power prices.The cooling load of the system can be used according to Equation ( 10): where CL t is the system total cooling load at hour t; hour t i is the time duration for which unit i has been on or off at time t, The chiller unit operation minimum on-time, T i,on , and minimum off-time, T i,o f f constraints are shown in Equations (11) and (12).The l loss is transmission loss for the ice storage system and l loss is set to 0.05.The ice storage tank for the charge process and discharge process for a day can be used according to Equation (13).

Proposed Methodology
This study proposes an IRBSO in order to achieve the optimal demand response of an ice storage air-conditioning system.The objective is to minimize cost in a daily schedule while satisfying all constraints including cooling loading.The bee swarm algorithm significantly changes the operation mode of bee swarms, which can be divided into forager bees, onlooker bees and scout bees.To further strengthen the global search capability, each bee type uses a different approach.The input data include system total cooling load, dynamic electricity price, operation status of each chiller and cooling load, temperature of chilled water, ice storage charge and discharge processes and the temperature of the ice storage water.The IRBSO parameters include the forager bees, onlooker bees and scout bees, with the total population of bees as X, the number of as iter, the nonlinear ripple weight factor as H and the interference factor as pr.
) be an individual, x = 1, 2, ..., X and X is set to 100 in this paper.The population sizes in the IRBSO algorithm are forager bees 40%, onlooker bees 40% and scout bees 20%.All individuals are set between the lower and upper limits with a uniform distribution, as shown in Equation ( 14 where rand is a random variable between 0 and 1.The fitness score of each B t k is obtained by calculating the objective function.

Forager Bees
The messages of the releaser pheromones and primer pheromones spread in the message exchange area are retrieved by experienced forager bees.The forager bees remember the best location in their own searches, Pbest t x 1 , and x 1 is the number of forager bees.This method can effectively strengthen the search range and accuracy, resulting in a higher probability of finding the global optimal solution.The forager bees follow the location of the current optimal solution to search Gbest t , as shown in Equation ( 15): where r b and r g are two random numbers between 0 and 1.Some bees flying over some parts of the search space may acquire profitable information from other bee swarms.The bees need to avoid premature convergence and to enlarge the search area to cover places that might have been neglected.The w b and w g determine the importance of the social and cognitive information for each iteration.
In this study, w b,max w b,min , w g,max , w g,min and itermax are set to 1.5, 0.8, 0.8, 1.5 and 200, respectively.
In this study, we used the nonlinear sine function and temporal variation as givens in Equations ( 16) and (17).Based on the resonance of weight factors w br f and w gr f , staggered search results can be achieved in which the nonlinear ripple weight factor is more effective at making the algorithm search better both locally and globally: where H is set to 0.05, 0.1 and 1; the non-linear ripple weight factors can be derived from the w br f and w gr f ripple resonance curves placed into w b and w g , as in Equations ( 18) and ( 19), Figures 2 and 3 where rand is a random variable between 0 and 1.The fitness score of each t k B is obtained by calculating the objective function.

Forager Bees
The messages of the releaser pheromones and primer pheromones spread in the message exchange area are retrieved by experienced forager bees.The forager bees remember the best location in their own searches, Pbest 1 , and x1 is the number of forager bees.This method can effectively strengthen the search range and accuracy, resulting in a higher probability of finding the global optimal solution.The forager bees follow the location of the current optimal solution to search Gbest t , as shown in Equation ( 15): where rb and rg are two random numbers between 0 and 1.Some bees flying over some parts of the search space may acquire profitable information from other bee swarms.The bees need to avoid premature convergence and to enlarge the search area to cover places that might have been neglected.The wb and wg determine the importance of the social and cognitive information for each iteration.In this study, wb,max wb,min, wg,max, wg,min and itermax are set to 1.5, 0.8, 0.8, 1.5 and 200, respectively.
In this study, we used the nonlinear sine function and temporal variation as givens in Equations ( 16) and (17).Based on the resonance of weight factors wbrf and wgrf, staggered search results can be achieved in which the nonlinear ripple weight factor is more effective at making the algorithm search better both locally and globally: (16) brf grf w w   (17) where H is set to 0.05, 0.1 and 1; the non-linear ripple weight factors can be derived from the wbrf and wgrf ripple resonance curves placed into wb and wg, as in Equations ( 18) and ( 19), Figures 2 and 3

Onlooker Bees
The onlooker bees in the improved bee swarm algorithm follow the experienced forager bees to obtain nectar information.Equation ( 20) is used for judging whether or not to follow the forager bees for foraging.The onlooker bees use the probabilistic selection method, as shown in Equation ( 21), to follow the experienced forager bees to gather nectar.In the working mode of the onlooker bees, the interference factor is also included in order to enlarge the search area: x is the number of onlooker bees; and r is a random number between 0 and 1.

Scout Bees
In the bee swarm algorithm, the model of the scout bees is no longer a baseless random search.In [20], the working model of the scout bees is modified to the average value of the global optimum solution and all swarm locations, as in Equations ( 22) and (23): where 3 x is the number of scout bees; and M is the average of all variable solutions in the iteration.

Self-Adaptation Repulsion Factor of Bee Swarm
To further strengthen the global search capability, the experienced forager bees use different approaches.This procedure causes them to fly over some parts of the search space and may include profitable information from the bee swarm.The increasing diversity of the bee swarm is incorporated in order to avoid premature convergence.To enlarge a search area that might have been neglected, a relevant study suggests the concept of the interference factor, sign.When the randomly generated ) 1 , 0 ( rand is larger than the predefined pr, a reverse search, as given in Equation (24), takes place and pr is initially set to 0.8:

Onlooker Bees
The onlooker bees in the improved bee swarm algorithm follow the experienced forager bees to obtain nectar information.Equation ( 20) is used for judging whether or not to follow the forager bees for foraging.The onlooker bees use the probabilistic selection method, as shown in Equation ( 21), to follow the experienced forager bees to gather nectar.In the working mode of the onlooker bees, the interference factor is also included in order to enlarge the search area: B t x 2 ,iter`1 " B t x 2 ,iter `wg ˆr ˆpB t probpeq ´Bt x 2 ,iter q where prob peq is the fitness value of the source food; B t probpeq is the best source; x 2 is the number of onlooker bees; and r is a random number between 0 and 1.

Scout Bees
In the bee swarm algorithm, the model of the scout bees is no longer a baseless random search.In [20], the working model of the scout bees is modified to the average value of the global optimum solution and all swarm locations, as in Equations ( 22) and (23): x 3 ,iter`1 " B t x 3 ,iter `r ˆpGbest t ´l ¨Miter q l " round p1 `randp0, 1qq where x 3 is the number of scout bees; and M is the average of all variable solutions in the iteration.

Self-Adaptation Repulsion Factor of Bee Swarm
To further strengthen the global search capability, the experienced forager bees use different approaches.This procedure causes them to fly over some parts of the search space and may include profitable information from the bee swarm.The increasing diversity of the bee swarm is incorporated in order to avoid premature convergence.To enlarge a search area that might have been neglected, a relevant study suggests the concept of the interference factor, sign.When the randomly generated randp0, 1q is larger than the predefined pr, a reverse search, as given in Equation (24), takes place and pr is initially set to 0.8: Energies 2016, 9, 113 8 of 16 We modified the interference factor, sign, to the self-adaptation repulsion factor, sign, with an initial setting of pr " 0.7 and the random variable randp0, 1q for the total number of bee swarms.The sign values used by the bee swarms were recorded, and the pr value based on self-adaptation repulsion was adjusted based on the fitness value for iteration, as in Figure 4.
if Min Cost comes from sign " 1 pr " pr `0.1 if pr ą pr max , then pr " pr max and ct " ct `1 otherwise ct " 0 if ct " ct max , then pr " 0.7 and ct " 0 otherwise Min Cost comes from sign " ´1 pr " pr ´0.1 if pr ă pr min , then pr " pr min and ct " ct `1 otherwise ct " 0 if ct " ct max , then pr " 0.7 and ct " 0 end In this paper, pr max " 0.9 and pr min " 0.1 to prevent pr from becoming extreme, and ct max " 5.The updating velocities used in this study improved the diversity of the solutions.This behavior is referred to as the self-adaptation repulsion factor.
Energies 2016, 9, 113 8 of 16 We modified the interference factor, sign, to the self-adaptation repulsion factor, sign, with an initial setting of The sign values used by the bee swarms were recorded, and the pr value based on self-adaptation repulsion was adjusted based on the fitness value for iteration, as in Figure 4 The updating velocities used in this study improved the diversity of the solutions.This behavior is referred to as the self-adaptation repulsion factor.The effect of the repulsion factor for the bee swarm can be expressed as Equations ( 25) to (27):

Simulation Results
In this study we investigated the actual practice in the operation planning of the chiller units and ice storage tank for a hospital.IRBSO, TVAC-PSO and PSO were used to calculate the optimal single-day cooling loads of the dynamic electricity price in summer, along with the derived switch status and which include six chillers and one ice storage tank.The TOU rates for the chillers and ice storage tank was calculated based on the summer date listed in Table 1.The effect of the repulsion factor for the bee swarm can be expressed as Equations ( 25) to (27): x 2 ,iter`1 " B t x 2 ,iter `sign ˆwg ˆr ˆpB t probpeq ´Bt x 2 ,iter q (26) B t x 3 ,iter`1 " B t x 3 ,iter `sign ˆr ˆpGbest t ´l ¨Miter q (27)

Simulation Results
In this study we investigated the actual practice in the operation planning of the chiller units and ice storage tank for a hospital.IRBSO, TVAC-PSO and PSO were used to calculate the optimal single-day cooling loads of the dynamic electricity price in summer, along with the derived switch status and which include six chillers and one ice storage tank.The TOU rates for the chillers and ice storage tank was calculated based on the summer date listed in Table 1.All the simulations were carried out with MATLAB 7.6 on a Core i5, 3.0 GHz personal computer with 8 GB of RAM.The relevant information for each chiller and the ice storage included the hourly recorded operating coefficients from March to December in 2013 for the data modeling.The cooling capability for each refrigerating ton (RT) is 12,658.46kJ/h and operating parameters of the 6 chiller units and the ice storage tank of 8000 RT.By using the measurement data, the least squares method (LSR) was used to get the input-output (I/O) operation curves of the chillers and the ice storage tank as shown in Tables 2 and 3. Customers operating ice storage central air conditioning system are eligible for a 40% discount on the energy charge for the system's off-peak kWh consumption in Taiwan.The ice storage tank model included the demand response, the average discharge and the situation without an ice storage tank to determine the dynamic electricity price in the case tests.Table 4 shows the respective simulation results for the objective programming of the case for the total scheduling of a day.That minimized the total cost and dispatch of an ice storage tank in summer.Different models had different results, and all constraints were met for the case study.The results obtained by the IRBSO for this case and other methods, which have been reported in the literature, are summarized in Table 4, which shows the total cost values achieved using the IRBSO algorithm for the demand response, average discharge and no ice storage tank to be better than those of the PSO, TVAC-PSO and BSO methods.From Figure 5, the ice storage tank of demand response supplies the required cooling load during peak periods and all chillers are appropriately dispatched to achieve minimal cost on a summer day for IRBSO.The improvement of the IRBSO over other algorithms was clear.Figure 6 illustrates the convergence characteristics of PSO, TVAC-PSO, BSO and IRBSO in the demand response scenario, which also showed the capacity of IRBSO to explore a more likely global optimum.The demand response modeling total cost of 201142 (NT$) for the IRBSO used in case study.The results, showed the capability of IRBSO to explore a more likely global optimum.Table 5 shows the ice storage air conditioning of average discharge and Tables 6 and 7 show the IRBSO at the dynamic electricity price intervals for demand response.Table 5 shows the ice storage air conditioning of average discharge and Tables 6 and 7 show the IRBSO at the dynamic electricity price intervals for demand response.Table 5 shows the ice storage air conditioning of average discharge and Tables 6 and 7 show the IRBSO at the dynamic electricity price intervals for demand response.Figure 7 shows the respective approximate thermal energy stored in the ice storage tank during the charging and discharging processes, i.e., the building cooling load.As Figures 7 and 8 show, during off-peak hours for the cooling load of chiller units, the ice-storage tank is in a charging process.During peak-load hours, the discharging process of the ice-storage tank provides the building's required cooling load.As the results show, the ice storage air-conditioning system utilizing IRBSO for the demand response evidenced better planning of the charging and discharging procedures.Figure 7 shows the respective approximate thermal energy stored in the ice storage tank during the charging and discharging processes, i.e., the building cooling load.As Figures 7 and 8 show, during off-peak hours for the cooling load of chiller units, the ice-storage tank is in a charging process.During peak-load hours, the discharging process of the ice-storage tank provides the building's required cooling load.As the results show, the ice storage air-conditioning system utilizing IRBSO for the demand response evidenced better planning of the charging and discharging procedures.

Simulation of Demand Response
From Table 8, the total power consumption of a day did not clearly change but Table 9 clearly shows the maximum peak-load reduction in the summer, with the drop in power cost dependent on the dynamic electricity price.Table 9 shows the load range of the demand response to be 449.73kW less than the average discharge of 601.72 kW, and 1282.13 kW with no ice storage tank in the summer.As Tables 7 and 9 and the Figure 9 show, the average peak load of demand response average was 2809.27kW less than the average discharge of 324.11 kW, and that with no ice storage tank as 624.91 kW.The system's maximum partial-peak load for the demand response was higher than the average discharge, but less than with no ice storage tank.There was a clear load shifting of the demand response in summer, as shown in Figure 9.
and the Figure 9 show, the average peak load of demand response average was 2809.27kW less than the average discharge of 324.11 kW, and that with no ice storage tank as 624.91 kW.The system's maximum partial-peak load for the demand response was higher than the average discharge, but less than with no ice storage tank.There was a clear load shifting of the demand response in summer, as shown in Figure 9.

Conclusions
This paper has presented the IRBSO algorithm to solve the ice storage air-conditioning system simulation and dynamic electricity price for the demand response problem.The IRBSO improved

Conclusions
This paper has presented the IRBSO algorithm to solve the ice storage air-conditioning system simulation and dynamic electricity price for the demand response problem.The IRBSO improved moving patterns and repulsion techniques for the diversity of solutions.In the test cases, the operating coefficients and cooling capacity of the chillers and ice storage tank were measured, and LSM was used to construct the curve model of the power consumption and cooling load.Based on the dynamic electricity price and all technical constraints, the dispatch model of the ice storage air-conditioning system was formulated by considering the ice storage discharging scheduling of demand response.This method could also improve the operating efficiency of ice storage and air-conditioning equipment.The planning of the charge and discharge processes of the ice storage tank for demand response study enhanced the operating efficiency of the chillers.This paper also provide greater energy efficiency in dispatching chillers and ice storage, thus reducing a user's electricity bill.The actual cases were used to verify the effectiveness of the proposed method.
number of bee swarms.

Figure 5 .
Figure 5.The dynamic electricity price on a summer day (16 July 2013) for IRBSO.

Figure 5 .
Figure 5.The dynamic electricity price on a summer day (16 July 2013) for IRBSO.

Figure 7 .
Figure 7. Cooling capability provided by each chiller and the ice storage tank on a summer day (16 July 2013).

Figure 7 .
Figure 7. Cooling capability provided by each chiller and the ice storage tank on a summer day (16 July 2013).

Figure 7 .
Figure 7. Cooling capability provided by each chiller and the ice storage tank on a summer day (16 July 2013).

Figure 8 .
Figure 8. Operating types of ice storage tank in summer.

Figure 8 .
Figure 8. Operating types of ice storage tank in summer.

Figure 9 .
Figure 9. Simulation of dynamic electricity price.

Figure 9 .
Figure 9. Simulation of dynamic electricity price.
show the behavior of wb and wg with different values of H: .

Table 3 .
Chillers 1 to 6 and ice storage tank.

Table 4 .
The ice storage air conditioning with dynamic electricity price on a summer day (16 July 2013).
* Average Discharge: The cooling load of discharge is identical.

Table 5 .
The IRBSO for ice storage air conditioning of average discharge.The dynamic electricity price on a summer day (16 July 2013) for IRBSO.

Table 5 .
The IRBSO for ice storage air conditioning of average discharge.

Table 5 .
The IRBSO for ice storage air conditioning of average discharge.

Table 6 .
The IRBSO for ice storage air conditioning of demand response.

Table 7 .
Power loads of demand response for IRBSO.

Table 8 .
Power load for dynamic electricity price in case study.

Table 9 .
Power load and dynamic electricity price for case results.