The Optimal Pumping Power under Different Ice Slurry Concentrations Using Evolutionary Strategy Algorithms

: A suitable ice slurry ﬂuid with a suitable ice concentration ratio can save operational costs. The design of the optimal ice slurry concentration focuses on ﬁnding an evolution strategy, which can further minimize the power consumption of the pump. A theoretical model was established to simulate the effect of different ice concentrations and ﬂow rates on the performance of the pump. The data obtained were ﬁtted by curve-ﬁtting function. The process was modeled in the MATLAB evolutionary strategy algorithm to obtain the conﬁguration scheme of the ice concentration and ﬂow under different refrigeration capacities. The simulation results showed that when the required cooling capacity was 13.889 kWh, ice concentration was set to 19.68%, and ﬂow rate was set to 2.1075 × 10 − 4 m 3 /s, the power consumption could be reduced by 23%.


Introduction
Owing to the intensification of global warming effects, the international community has paid more and more attention to the reduction of greenhouse gas emissions, and the role of secondary refrigerant in industrial and commercial refrigeration has also drew more and more attention [1]. The awareness of the environmental impact of fluorinated gases (F-gases) used in refrigeration is instigating the development of technologies to recover and recycle them [2]. There are extensive opportunities to reduce emissions using existing technology and alternative substances with low global warming potential [3]. The increasing awareness of the release of fluorinated gases (F-gases) into the atmosphere is also instigating the development of techniques to secondary refrigerant [4].
Ice slurry, as an excellent cold carrier, is extremely important in refrigeration applications [5]. Ice crystal particles of fluidized ice will undergo phase change during heat transfer, which can quickly cool down and meet the requirements of high cooling load. It has the advantages of high energy storage density, good heat transfer performance, and large latent heat of phase change [6]. In practical application, the diameter of the conveying pipe, the energy consumption of the pump, and the structural size of the heat exchanger can be greatly reduced. Under the condition of the same size of the pipe, it can carry more cooling capacity than that of the liquid in the same temperature range [7]. In the application of ice slurry systems, the energy consumption accounts for a considerable part of the operating cost, and the optimal design and control of ice slurry systems is a promising solution to reduce the operating costs and the peak power consumption. Vetterli [8] paid attention to the thermal modeling of the building using mixed-integer Linear programming to solve the optimization design problem, using this model to derive optimal designs under various electric tariff schemes. In [9], Henze compared four different control strategies using optimal control strategies as a benchmark. In [10], Lee focused on the optimal design of an ice storage system using particle swarm algorithms.
Ice slurry has received increasing attention for air conditioning systems in buildings due to the widespread concerns over shifting peak load, environmental friendliness, and minimizing tank storage. CAPCOM in Japan is the one of successful installations of the ice slurry systems, similarly to the Crystal Tower of Osaka [11].
This study took the headquarters building of CAPCOM in Japan as an example. For an ice slurry refrigeration system, its energy consumption accounts for 4% of the total energy consumption of the building [12], among which, due to the high floor, the transportation loss is huge. Therefore, studying the pump power under the ice slurry with different ice mass concentrations, so as to improve the pumping efficiency of the ice slurry system can significantly save the transportation cost of the ice slurry. Many researchers conducted this kind of study previously. Frei [13] found that the performance of the pump decreased with the increase of the ice concentration. Nørgaard [14] measured the pressure difference between the inlet and outlet of the test pump to obtain roughly the actual head of the pump in operation, and then measured the pump flow rate to obtain pump performance curves with different ice concentrations. The influence of ice slurry on the pump power is not a simple linear problem [15]. The optimal design of an ice slurry system should analyze and support a design process characterized by low energy consumption [16] through the optimal selection of power, e.g., pumps [17][18][19], drive unit [20][21][22], or control algorithm [23][24][25], depending on the conditions of use. Currently, most studies on ice slurry systems are based on experiments, and the microscopic characteristics of ice particles passing through different channels (vertical Slit channel [26], Rectangular Channel [27], etc.) are studied, but there are few articles about the application of algorithms to improve the characteristics. In order to solve ice slurry flow problem, this paper adopted an evolutionary strategy algorithm to study the optimal pump power under different flow ice concentrations as a black box optimization problem. This paper studied a typical ice slurry cooling system from the angle of economic benefit and cost, sought its optimal pump power, and provided suggestions for the operation of the ice slurry system.

Ice Slurry Refrigeration System
The CAPCOM building [28] is located in the central district of Osaka, Japan, and is the R&D department of CAPCOM. It has 20 floors and a total floor area of 16,784 m 2 . The company's offices are located on floors 2 to 16, each with an area of 570 m 2 . It has been calculated that the heat load in the office is 151 W/m 2 , which is relatively higher than that of the normal office building, due to the heat generated by the computers in the office. Thus, a gas-driven absorption refrigeration system was installed to provide cooling for the office area.
The cooling of the internal area of the office space of the building is provided by an ice slurry cooling system, as shown in Figure 1. The system consists of two ice slurry units, each with a cooling capacity of 272 kW. The ice slurry machine on the 17th floor produces ice stored in two tanks. The ice slurry from the storage tank is pumped to the air handling units on each floor through the ice slurry separator, which controls the ice mass fraction of the ice slurry up to 20%. Each of the 15 office floors has two air handling units with flow rates of 2.5 m 3 /s, which in turn provide six variable air volume (VAV) terminals. The low temperature of the ice slurry led to the choice of a cold air distribution temperature of 12 • C rather than 15 • C, which is a normal design practice in Japan. This reduced air flow requirements from 41 m 3 /h to 32 m 3 /h, resulting in smaller equipment sizes and lower power requirements for air distribution, and savings in capital and operating costs. cooling energy battery of the building. Thus, the energy cost and electricity bills are decreased by using this approach. Through the energy analysis of the system, the 16 layers that carry the ice slurry pipeline cost the main energy consumption. It is of great significance to study how to design the ice concentration and flow rate to reasonably reduce the pump power of the ice slurry system at different cooling capacity rates.

Simulation Models/Methods
According to their properties, fluids are usually classified as Newtonian or non-Newtonian. At present, the research on the flow characteristics of ice slurry, especially the non-Newtonian fluid characteristics, usually adopts two methods: fitting experimental data to obtain relevant empirical or semi-empirical formulae or using the existing non-Newtonian flow model to carry out theoretical calculation and analysis. Due to the limitation of experimental equipment, selection of working conditions, and slurry variability, the empirical or semi-empirical formulas obtained in the study have a certain applicable scope, and it is very difficult to determine the flow model of ice slurry accurately. Therefore, the work of this paper was not to study the pumping situation of ice slurry in the pipeline, but to The working process of the ice slurry cooling system is shown in Figure 1 [29]. During the night time, the ice slurry was produced and stored in the tank as the electricity cost is relatively low during night time. In the daytime, owing to the increasing needs of the cooling demands, the ice slurry was pumped to the VAV terminals in each floor to compensate for the rising temperature in the offices. The ice slurry tank can be regarded as a cooling energy battery of the building. Thus, the energy cost and electricity bills are decreased by using this approach. Through the energy analysis of the system, the 16 layers that carry the ice slurry pipeline cost the main energy consumption. It is of great significance to study how to design the ice concentration and flow rate to reasonably reduce the pump power of the ice slurry system at different cooling capacity rates.

Simulation Models/Methods
According to their properties, fluids are usually classified as Newtonian or non-Newtonian. At present, the research on the flow characteristics of ice slurry, especially the non-Newtonian fluid characteristics, usually adopts two methods: fitting experimental data to obtain relevant empirical or semi-empirical formulae or using the existing non-Newtonian flow model to carry out theoretical calculation and analysis. Due to the limitation of experimental equipment, selection of working conditions, and slurry variability, the empirical or semi-empirical formulas obtained in the study have a certain applicable scope, and it is very difficult to determine the flow model of ice slurry accurately. Therefore, the work of this paper was not to study the pumping situation of ice slurry in the pipeline, but to determine the optimal concentration and flow of the pump inlet so that it can be transported to the pipeline through the internal pressure of the pump so as to reduce the consumption of pumping power and save energy.
The model of a single pump in the ice slurry system can be simplified as Figure 2, a simple ice slurry conveying pipeline composed of a storage tank, a transport pump, a test pump, and a control valve. The storage tank is the location of the ice slurry. The transport pump refers to all the pumps except the pump under study, and the test pump is the pump under study. The two flow meters and two pressure detectors in the figure are for data collecting. In this paper, the pumping model of ice slurry was studied by analyzing the experimental data of Nørgaard. The pump used in the experiment was a GrundfosCr2-50 centrifugal pump, in which the liquid solution was 16% propanediol, the diameter of the ice particles was between 0.1 mm and 0.3 mm, and the pump speed was 2800 r/min. Their experimental figure was selected as experimental data, and the data are shown in Tables 1 and 2. According to the data in Tables 1 and 2, MATLAB was used for data fitting, and the empirical formula related to ice concentration could be obtained by the pump efficiency calculation formula: and the empirical formula related to ice concentration could be obtained by the pump efficiency calculation formula: The ice concentration is defined by: where: -the ice concentration in ice slurry (-), -the mass of ice particles (kg), -the mass of ice slurry fluid (kg). Pump efficiency calculation formula: where: -the pump efficiency (%), -the pump power output (W), -the shaft power of the pump (W). The pump power output can also be written as: where: -the pump power output (W), -the fluid density (kg/m 3 ), -the volume flow rate (m 3 /h), -the gravitational acceleration (m/s 2 ), -the total head of the pump (m).      The ice concentration is defined by: where: C-the ice concentration in ice slurry (−), m ice -the mass of ice particles (kg), m-the mass of ice slurry fluid (kg). Pump efficiency calculation formula: where: η-the pump efficiency (%), P e -the pump power output (W), P-the shaft power of the pump (W). The pump power output can also be written as: where: P e -the pump power output (W), ρ-the fluid density (kg/m 3 ), Q-the volume flow rate (m 3 /h), g-the gravitational acceleration (m/s 2 ), H-the total head of the pump (m).
In ice slurry two-phase flow, the density of ice crystal and solution is different to some extent, and the density of ice slurry can be calculated by linear weighting of the two phase states: where: C-the ice concentration (−), ρ-the fluid density (kg/m 3 ), ρ s -the density of the solution (kg/m 3 ), ρ i -the density of ice crystals (kg/m 3 ), According to the ASHRAE 2020 manual, the density of 16% propylene glycol solution is about 1020 kg/m 3 , and the density of ice is about 917 kg/m 3 .
Substituting Equation (3) into Equation (2), we can acquire [30]: where: η-the pump efficiency (%), Q-the volume flow rate (m 3 /h), P-the power consumption (i.e., the pump power) (W), H-the total head of the pump (m), ρ-the fluid density (kg/m 3 ), g-the gravitational acceleration (m/s 2 ). By using the curve fit tool of MATLAB to fit the data in Table 2, we can acquire Figure 3:  The R-square is 0.9117. R-square is generally used to evaluate the quality of the model in regression. The closer its value is to 1, the better the model performance is.
The fitting results are as follows:  (6) where: -the volume flow rate (m 3 /h), -the total head of the pump (m), -the ice concentration (-). Such empirical formulas can only be applied in the range of 0% to 30% ice concentration.
By mathematically fitting the data of ice concentration and efficiency in Table 1 in the same way, we can acquire Figure 4 and Equation (7):  The R-square is 0.9117. R-square is generally used to evaluate the quality of the model in regression. The closer its value is to 1, the better the model performance is.
The fitting results are as follows: where: Q-the volume flow rate (m 3 /h), H-the total head of the pump (m), C-the ice concentration (−). Such empirical formulas can only be applied in the range of 0% to 30% ice concentration. By mathematically fitting the data of ice concentration and efficiency in Table 1 in the same way, we can acquire Figure 4 and Equation (7):  The R-square is 0.9117. R-square is generally used to evaluate the quality of the model in regression. The closer its value is to 1, the better the model performance is.
The fitting results are as follows:  (6) where: -the volume flow rate (m 3 /h), -the total head of the pump (m), -the ice concentration (-). Such empirical formulas can only be applied in the range of 0% to 30% ice concentration.
By mathematically fitting the data of ice concentration and efficiency in Table 1 in the same way, we can acquire Figure 4 and Equation (7):  The R-square is 0.984. This shows that our model performs well.
Such an empirical formula of efficiency is applicable in a certain range and meaningless in another range. Meaningful ranges mainly include: the efficiency range being between 0 and 1, and therefore, the efficiency beyond this range is meaningless. Due to the performance of the pump, it does not make sense to exceed the flow rate and ice concentration that the pump can transport.
Substituting the above empirical head and efficiency Equations (5) and (6) into the theoretical pump efficiency Equation (4) mentioned above, the pump power can be obtained as follows: where: Q-the volume flow rate (m 3 /h), ρ-the fluid density (kg/m 3 ), C-the ice concentration (−). We use (µ,λ) evolutionary strategy algorithm [31] for optimization. Here are the main steps of the algorithm [32]: First, the problem to be solved by the algorithm must be determined. The problem of solving the optimal concentration into a mathematical problem should be transformed: solving the minimum value of Equation (7). The algorithm flow chart is shown in Figure 5. as to carry out the next evolution. Eventually, generation after generation, the population must converge to the optimal solution we need. Essentially, this is an evolutionary strategy algorithm, a probabilistic algorithm that takes a random value and then calculates a comparative fitness [35]. The advantage of this algorithm is that the optimal concentration of ice slurry is a practical problem, and its mathematical model is necessarily very complicated, which is difficult to analyze by ordinary methods. By using the evolutionary strategy algorithm, the function properties can be ignored, the optimal value can be obtained easily [36], and the problem of local optimization can be avoided. However, with the increase of iteration times, the optimal value becomes more accurate, but the running time is also greatly increased.

Empirical Function Diagram
MATLAB was used to draw Equations (6) and (7). The drawing range is 0~30% of ice  Set up an initial population, which is the first generation of parents for evolution (which can be randomly generated). The parent should be a matrix composed of concentration and flow, the first behavior is concentration, and the second behavior is flow, one of which represents a set of solutions to the function P. Considering the running time and the accuracy of the results, the initial population size as 10 was set. That is, there are 10 solutions of P functions in the initial population matrix.
Next, the evolution of the population should be started [33].
Step one, cross out parents. For example, two groups of random solutions (concentration 0.1 flow 2) and (concentration 0.2 flow 1) in the parent, then the concentration or flow of the two groups of solutions needed to be interchanged, which becomes (concentration 0.1 flow 1) and (concentration 0.2 flow 2). This is to make the range of data solution more extensive, to avoid only small concentration against small flow, large concentration against large flow situations.
Step two, the program has to start mutating, which is basically adding a variable to each of the two quantities of each set of solutions. This variable satisfies a Gaussian distribution of some zero mean and some variance. The variance is called the degree of variation in the evolutionary strategy algorithm. So, if you add this variable, you are adding a degree of variation to each of the solutions, which is an indication of how much the solutions have changed. When the population evolves to convergence, the degree of variation will also gradually decrease, making the whole population converge easier. In general, crossover and variation are designed to make the values of concentration and flow more "random" within the allowable range. In this way, the group of solutions which makes the pump power minimum has a higher probability of appearing.
The data completed by mutation, namely, the offspring, were selected. Generally speaking, the ratio of parents to offspring in (µ, λ) selection strategy was 1:7 [34]. Since the parent was set as 10 above, the offspring should be 70. The pump power matrix corresponding to the subgeneration solution can be obtained by calculating 70 subgenerations into the pump power formula. The pump power matrices of 70 groups of solutions were compared, and the 10 groups with the least pump power were finally selected. These 10 groups of solutions were taken as new parents, namely the second generation parents, so as to carry out the next evolution. Eventually, generation after generation, the population must converge to the optimal solution we need.
Essentially, this is an evolutionary strategy algorithm, a probabilistic algorithm that takes a random value and then calculates a comparative fitness [35]. The advantage of this algorithm is that the optimal concentration of ice slurry is a practical problem, and its mathematical model is necessarily very complicated, which is difficult to analyze by ordinary methods. By using the evolutionary strategy algorithm, the function properties can be ignored, the optimal value can be obtained easily [36], and the problem of local optimization can be avoided. However, with the increase of iteration times, the optimal value becomes more accurate, but the running time is also greatly increased.

Empirical Function Diagram
MATLAB was used to draw Equations (6) and (7). The drawing range is 0~30% of ice concentration and 0~5 m 3 /h of flow, and Figure 6a,b can be obtained: As can be seen from Figure 6a, when the concentration increased, the head of the pump did not change much at the same flow rate, but the flow that the pump can transport was even smaller. When the ice concentration reached the system limit of 30%, the pump flow was only about 0.5 m 3 /h. The reason for this situation is that the power of the pump has a limit. Under the condition of a certain speed, with the increase of concentration, the energy consumption of conveying ice slurry per unit flow rate is higher. In order to make the ice slurry reach the required head, only a small part of the ice slurry can be transported at the expense of the flow. It also shows that the properties of ice slurry are different from an ordinary water solution. Obviously, the energy consumption of conveying water is much lower than that of conveying ice slurry. Due to the Q-H characteristics of the pump itself, the flow rate must decrease with the rise of the head, but at high ice concentrations, the flow rate decreases more rapidly. Thus, in general, as the concentration increases, the pump can deliver less flow. As can be seen from Figure 6a, when the concentration increased, the head of the pump did not change much at the same flow rate, but the flow that the pump can transport was even smaller. When the ice concentration reached the system limit of 30%, the pump flow was only about 0.5 m 3 /h. The reason for this situation is that the power of the pump has a limit. Under the condition of a certain speed, with the increase of concentration, the energy consumption of conveying ice slurry per unit flow rate is higher. In order to make the ice slurry reach the required head, only a small part of the ice slurry can be transported at the expense of the flow. It also shows that the properties of ice slurry are different from an ordinary water solution. Obviously, the energy consumption of conveying water is much lower than that of conveying ice slurry. Due to the Q-H characteristics of the pump itself, the flow rate must decrease with the rise of the head, but at high ice concentrations, the flow rate decreases more rapidly. Thus, in general, as the concentration increases, the pump can deliver less flow.
As can be seen from Figure 6b, when the concentration increased, the efficiency of the pump decreased continuously. That is to say, the higher the concentration, the more difficult the pump is to pump, which is consistent with the pump characteristics and the rise of the flow, so that the pump efficiency rose to the maximum value and then decreased. The flow rose to the level that pump can withstand the extreme value; at this time the pump power is infinite and meaningless, and the efficiency is zero.
According to the analysis of the above two figures, it is consistent with the influence of two-phase flow on the pump: when the concentration of conveying ice particles increased, the head decreased less at a small flow rate, but the drop was larger at a large flow rate. As the concentration increased, the maximum flow rate of the pump gradually decreased, and in order to achieve the required head, the power will increase and the efficiency will decrease. This also shows that such an empirical formula is valid in a certain range, and this analysis method has a certain practical significance.
Then, the pump power was analyzed in the same way, and Equation (8) was plotted using MATLAB. The drawing range was 0~30% ice concentration, 0~5 m 3 /h flow, and 0~750 W pump power, and Figure 7 could be obtained. As can be seen from Figure 6b, when the concentration increased, the efficiency of the pump decreased continuously. That is to say, the higher the concentration, the more difficult the pump is to pump, which is consistent with the pump characteristics and the rise of the flow, so that the pump efficiency rose to the maximum value and then decreased. The flow rose to the level that pump can withstand the extreme value; at this time the pump power is infinite and meaningless, and the efficiency is zero.
According to the analysis of the above two figures, it is consistent with the influence of two-phase flow on the pump: when the concentration of conveying ice particles increased, the head decreased less at a small flow rate, but the drop was larger at a large flow rate. As the concentration increased, the maximum flow rate of the pump gradually decreased, and in order to achieve the required head, the power will increase and the efficiency will decrease. This also shows that such an empirical formula is valid in a certain range, and this analysis method has a certain practical significance.
Then, the pump power was analyzed in the same way, and Equation (8) was plotted using MATLAB. The drawing range was 0~30% ice concentration, 0~5 m 3 /h flow, and 0~750 W pump power, and Figure 7 could be obtained. There are abnormal data in the upper left corner of Figure 7. This is because our els are based on experimental practice and only consider their mathematical chara tics. In fact, when the ice slurry transportation rate exceeds the pump's capac  There are abnormal data in the upper left corner of Figure 7. This is because our models are based on experimental practice and only consider their mathematical characteristics. In fact, when the ice slurry transportation rate exceeds the pump's capacity to transport the limit, the pump cannot work, which is shown in the mathematical model that the pump power is infinite. Obviously, too much flow and ice concentration can cause infinite pump power.

The Operation Range of Pump
Combining with the influence of flow rate and ice concentration, the flow rate of high ice concentration ice slurry transported by this system should not exceed 3.5 m 3 /h. The operating range of the pump is approximately shown in Figure 8. There are abnormal data in the upper left corner of Figure 7. This is because ou els are based on experimental practice and only consider their mathematical chara tics. In fact, when the ice slurry transportation rate exceeds the pump's capa transport the limit, the pump cannot work, which is shown in the mathematical that the pump power is infinite. Obviously, too much flow and ice concentration can infinite pump power.

The Operation Range of Pump
Combining with the influence of flow rate and ice concentration, the flow rate ice concentration ice slurry transported by this system should not exceed 3.5 m 3 operating range of the pump is approximately shown in Figure 8.  The actual operating range of the pump is larger than the theoretical operation, and the actual range of the flow is larger. The reason is that, at a certain flow rate and at the ice concentration (where the flow does not exceed the maximum capacity of the pump), the pump efficiency is not zero (the pump is still running at rated power or over-rated power), but the it is already very small, and the pump power required to reach the pump required flow, although not infinite, exceeds the pump's limit power, in which case the pump is prone to blockage, damage, and system failure. Therefore, in the design of the algorithm, we avoided this situation, so in theory, the flow rate of the pump runs in a smaller range [37].

Apply the Optimal Algorithm
To apply the optimal concentration algorithm, first we should define what is "optimal", which will be the key to the whole algorithm. The criterion of "optimal" was treated as the fitness in the evolutionary strategy algorithm. The process of finding the best fitness by the algorithm is the process of optimizing the concentration.
From the engineering practice, the optimal concentration is the concentration with the highest economic benefit. Therefore, it is defined that the concentration that minimizes the work of the conveying pump is the optimal concentration under the condition of satisfying an equal cooling capacity. In the algorithm, the pump power is the fitness. The smaller the pump power is, the better the fitness is.
The cooling capacity [38] of energy is set to meet the requirements. According to the energy absorbed by ice melting phase change and the conservation law of heat transfer, the following equation can be obtained: where: h-the enthalpy change (J), Q-the volume flow rate (m 3 /h), ρ-the fluid density (kg/m 3 ), t-the time (s), h ls -the latent heat of ice (kJ /kg), C-the ice concentration (−). The greater the change in enthalpy, the more energy is absorbed during the phase transition. Therefore, the larger is, the better its cold storage performance is. That is, when: where: h-the enthalpy change (J), Q-the volume flow rate (m 3 /h), ρ-the fluid density (kg/m 3 ), t-the time (s). C-the ice concentration (−). That is, it can be judged to satisfy the supply demand. This expression can be used as a selection condition of the algorithm. The optimal concentration was obtained by using the algorithm.
Taking h = 13.889 kWh (50,000 kJ) and the maximum operating power of the pump as 750 W, we can obtain Figure 9 and C = 0.1968 (the ice concentration) and Q = 0.7587 (the volume flow rate).
Algorithm running results shows that when the required supply h = 13.889 kWh is met, the ice concentration is 19.68% and the volume flow rate is 0.7587 m 3 /h, that is, 2.1075 × 10 −4 m 3 /s, making the pump power a minimum of 417.7 W.

Comparison under Different Cooling Capacity
Adjust the cooling capacity and put it into the algorithm. After multiple calculations, the data shown in Table 3 were obtained.
It can be found from Table 3 that the optimal concentration is stable at about 20% with little change, indicating that the optimal concentration does not change with the cooling capacity, while the optimal flow rate changes greatly. For that, we take concentration for a constant 20%, will be cooling capacity and flow rate, pump power, and data fitting, and acquiring the graphics as shown in Figure 10; as can be seen, the relationship is linear, and with the rise of the cooling capacity, the best flow and pump power will rise to meet the needs of cooling capacity, but after the flow rapidly rises, it reached the peak. Pump power is slowly rising due to the limitation of its own pump power, but also slowly reaches the peak, no longer rising. Algorithm running results shows that when the required sup met, the ice concentration is 19.68% and the volume flow rate is 0.75 * 10 −4 m 3 /s, making the pump power a minimum of 417.7 W.

Comparison under Different Cooling Capacity
Adjust the cooling capacity and put it into the algorithm. After the data shown in Table 3 were obtained.
It can be found from Table 3 that the optimal concentration with little change, indicating that the optimal concentration does no ing capacity, while the optimal flow rate changes greatly. For that, for a constant 20%, will be cooling capacity and flow rate, pump p and acquiring the graphics as shown in Figure 10; as can be seen, th and with the rise of the cooling capacity, the best flow and pump p the needs of cooling capacity, but after the flow rapidly rises, it re power is slowly rising due to the limitation of its own pump p reaches the peak, no longer rising.    It should be noted that there is a limit to the cooling capacity h that the system supply for a period of time due to the upper limit of pump power. When applying algorithm, be careful not to set the h value of cooling capacity beyond the limit, otherw there will be operational error.

Conclusions
In this paper, a single pump in the ice slurry system was taken as the research ob and a set of optimal fluid pumping layout schemes were given by using the evolution strategy algorithm through the constructed pumping model. When the rated spee Grundfos CR2-50 centrifugal pump or other centrifugal pumps with similar perform was 2800 r/min, to meet the cooling capacity of 13.889 kWh, the ice concentration at 19. and the ice slurry fluid flow of 2.1075 * 10 −4 m 3 /s could minimize the pump operation save transportation cost. In order to facilitate comparison, when the required cooling pacity was constant, the required concentration and flow were calculated by the a rithm. Only 417.7 W power is needed to meet the requirements. The calculation show 23% reduction in power consumption. In addition, the optimum concentration, flow and pump power under different cooling capacities were analyzed, and the linear tionship was obtained. The ideas provided in this paper for the construction of ice sl system are as follows: When the cooling demand needs to be calculated, we can optim it from the algorithm level to get the best flow, concentration, and power so as to energy and avoid wasting resources.   It should be noted that there is a limit to the cooling capacity h that the system can supply for a period of time due to the upper limit of pump power. When applying the algorithm, be careful not to set the h value of cooling capacity beyond the limit, otherwise there will be operational error.

Conclusions
In this paper, a single pump in the ice slurry system was taken as the research object, and a set of optimal fluid pumping layout schemes were given by using the evolutionary strategy algorithm through the constructed pumping model. When the rated speed of Grundfos CR2-50 centrifugal pump or other centrifugal pumps with similar performance was 2800 r/min, to meet the cooling capacity of 13.889 kWh, the ice concentration at 19.68% and the ice slurry fluid flow of 2.1075 × 10 −4 m 3 /s could minimize the pump operation and save transportation cost. In order to facilitate comparison, when the required cooling capacity was constant, the required concentration and flow were calculated by the algorithm. Only 417.7 W power is needed to meet the requirements. The calculation showed a 23% reduction in power consumption. In addition, the optimum concentration, flow rate, and pump power under different cooling capacities were analyzed, and the linear relationship was obtained. The ideas provided in this paper for the construction of ice slurry system are as follows: When the cooling demand needs to be calculated, we can optimize it from the algorithm level to get the best flow, concentration, and power so as to save energy and avoid wasting resources.