Identification of Optimal Parameters for a Small-Scale Compressed-Air Energy Storage System Using Real Coded Genetic Algorithm

Compressed-Air energy storage (CAES) is a well-established technology for storing the excess of electricity produced by and available on the power grid during off-peak hours. A drawback of the existing technique relates to the need to burn some fuel in the discharge phase. Sometimes, the design parameters used for the simulation of the new technique are randomly chosen, making their actual construction difficult or impossible. That is why, in this paper, a small-scale CAES without fossil fuel is proposed, analyzed, and optimized to identify the set of its optimal design parameters maximizing its performances. The performance of the system is investigated by global exergy efficiency obtained from energy and exergy analyses methods and used as an objective function for the optimization process. A modified Real Coded Genetic Algorithm (RCGA) is used to maximize the global exergy efficiency depending on thirteen design parameters. The results of the optimization indicate that corresponding to the optimum operating point, the consumed compressor electric energy is 103 . 83 k W h and the electric energy output is 25 . 82 k W h for the system charging and discharging times of about 8.7 and 2 h, respectively. To this same optimum operating point, a global exergy efficiency of 24.87% is achieved. Moreover, if the heat removed during the compression phase is accounted for in system efficiency evaluation based on the First Law of Thermodynamics, an optimal round-trip efficiency of 79.07% can be achieved. By systematically analyzing the variation of all design parameters during evolution in the optimization process, we conclude that the pneumatic motor mass flow rate can be set as constant and equal to its smallest possible value. Finally, a sensitivity analysis performed with the remaining parameters for the change in the global exergy efficiency shows the impact of each of these parameters.


Introduction
The security and reliability of electricity grid need the introduction of a storage system [1]. To these two main objectives of energy storage one can add the reduction of greenhouse gas emissions. Quite important, due to the increased consumption of fossil fuels, the amount of CO 2 emitted increased up to threefold between 1960 and 2008 reaching today more than 32,000 million tons per year. The climate change observed due to these emissions drove many countries to turn to renewable energy (RE) sources for electricity production in order to retain global warming within a 2% range [2,3]. However, these sources are strongly related to meteorology and are intermittent [4]. One of the solutions developed to overcome the problem of intermittency is to couple them with an electrical storage system [5,6].
The storage system will be able to play two roles namely protection and production. By protection it is meant that the system must be able to quickly restore (the response time on the order of a few minutes [7]) the energy stored during the fluctuation of the resource. By production it is meant that the storage system must be able to produce and sustain independently during a sufficiently long period the demands in the absence of the total source.
In the worldwide industries of electricity both mature technologies are used for large scale electricity storage. They are pump hydro (PH) system and compressed air energy storage (CAES) systems [2,[8][9][10][11][12]. As opposed to the hydroelectric pumping stations, the storage systems with compressed air offer flexibility both in size (smaller volumes) and capacity (ranging from several hundred KW to MW). These advantages give CAES the opportunity to be coupled to the power generation system with renewable sources. That is why many researches propose hybrid wind/CAES systems or photovoltaic plant (PV)/CAES systems [13][14][15][16][17][18][19][20][21]. In periods of low and off-peak energy demand, the CAES system stores electricity in the form of compressed air in a natural or artificial tank. The stored compressed air is released and heated in a combustion chamber burning fossil fuel before being expanded in a turbine connected to a generator for electricity reproduction [11,22,23]. Many studies dealing with the partial or total replacement of combustion chamber exist in the literature [3,[24][25][26][27][28]. The heat generated during compression is stored and used to heat air before expansion; such systems are called adiabatic or advanced adiabatic compressed air energy storage (AA-CAES). Unfortunately, conventional and AA-CAES used natural reservoir (underground caverns, rock formations) for storing compressed air which reduces its penetration potentiality due to the geological restriction [29,30]. Therefore, the attention has been recently focused on the usage of artificial air-tanks. The resulting system is known as micro or small scale compressed air energy storage system (SS-CAES). Such system can be used at isolated sites with renewable energy sources or in the residential sector to store electricity during off-peak hours. Generally, in SS-CAES system, fuel combustion is not needed because the compression heat is collected, stored and re-used to heat the compressed air before being expanded in the turbine or the reciprocating air motor. If the cooling energy in the discharged air is collected, the SS-CAES may act as a tri-generative system, for simultaneous production of cold, heat and electricity [4,9,[31][32][33][34][35][36].
To evaluate the performance of SS-CAES system by means of numerical simulations, numerous thermodynamics models have been developed during recent years. Generally, these models are based on the first law of thermodynamics with imposed design and operating parameters of the analyzed systems [9,32,34]. Unfortunately, energy analysis does not provide the information about the locations of energy degradation in a process and does not quantify the irreversibility in different components of the storage system. Therefore, based on both the first and second laws of thermodynamics, exergy analysis appears to be a powerful tool to overcome the limitations of energy analysis [37].
The purpose of this study is to develop a realistic approach to investigate the performance of SS-CAES system using pressure vessels without fossil fuel. This approach is based on the exergy analysis method. The required equations for modeling different components of the system are presented. These equations are used to build the objective function which is the global exergy efficiency of the storage system. We aim to maximize this objective function depending on thirteen design parameters and seven constraints. These design parameters are respectively: number of compression stage (n), compressor pressure ratio (π), volume of air storage tank (V t ), pressure ratio of high-pressure and low-pressure expansion stages (β HP ,β LP ), inlet temperature of high-pressure and low-pressure expansion stages (T in HP ,T in LP ), isentropic efficiency of compressor (η IsC ), isentropic efficiency of pneumatic motor (η Is,m ), mechanical efficiency of compressor (η mC ), mechanical efficiency of pneumatic motor (η mm ), compressor and pneumatic motor mass flow rate (ṁ C ,ṁ m ). The ranges of each parameter have been defined in the light of available technology. The optimization is performed using a modified Real Coded Genetic Algorithm (RCGA) in which two crossover methods are randomly selected from one generation to another during evolution. In order to improve the search efficiency of the RCGA, the Dynamic Random Mutation (DRM) method was used. The coupled modified RCGA-DRM could effectively determine the set of optimal values of influencing parameters that maximizes the global exergy efficiency of the SS-CAES system without fossil fuel used. MATLAB R software is used for all computations.

System description
The system to optimize is shown in Figure 1. The system operates in two phases: charge and discharge. The first phase is also known as compression phase, which is composed of multistage (CS) reciprocating compressor and intercoolers (HE). In this phase, the compressor is powered by the electrical energy available on grid during off-peak load hours or by the electricity generated by a renewable energy source. Intercoolers are used to recover the compression heat with water as the heat transfer fluid. Hot water is stored in an isolated thermal storage tank (HWt). Cooling of compressed air after each stage has the advantage of reducing the required electrical power of the compressor and increasing the compressed air storage tank efficiency due to the high density of the cold compressed air. High pressure cold air from compression and cooling is stored in compressed air storage tank (CASt). In conventional or diabatic CAES systems, at peak load hours, the stored air is released from the underground cavern, throttled through regulating valve and heated by fossil fuel burning in the combustion chamber (CC). But, in an AA-CAES and SS-CAES, the hot water produced during the compression phase can be used to heat air before each stage of expansion through heat exchangers, thus replacing combustion chamber. Two stages pneumatics motor (HPe and LPe) coupled to electric generator (G) are used to achieve the expansion process and generate electricity. In Figure 1, the points i for i ∈ {1, 2, ..., 16} denote the states of air tranformation during all storage process and i ∈ {17, 18, ..., 26} denote the states of water transformation during all transformation storage process.

Thermodynamic modelling
To carry on the optimization problem, it is necessary to build the objective function which in this study is the global exergy efficiency. The analytical form of the objective function results from the thermodynamic analysis of each component of the system. Due to the complexity of the system some simplifying but basic assumptions are made following previous works on thermodynamic modelling of energy systems with air as working fluid [38][39][40][41][42]: • All system components operate under a steady state condition except the CASt for which a dynamic modeling is performed to find the filling and discharge time together with the mean temperature in the tank during the discharge process.
• Air is assumed dry and modelled as an ideal gas.
• Potential and kinetic energy effects are negligible in the energy and exergy balances.
• Pressure drop in the components of system is neglected.
• The reference environment state conditions are T 0 = 20 • C and P 0 = 1.01bar which is also used as the system boundary for energy and exergy analyses. • Isentropic efficiency is assumed constant for compressor and pneumatic motor.
• The mass flow rate of the cooling water is set as constant in every intercoolers. Despite the different pressure ratio of expansions stages, we set the mass flow rate of the heating water in each heater as constant and equal to Six times that of the cooling water in order to increase the heat transfer rate and for more simplicity in modeling. • All the interpolations in thermodynamics tables to find the thermodynamics properties is done using MATLAB built-in function "interp1", with "spline" method. The thermodynamics tables are those of Moran book [43].

Energy analysis
In this section, the First Law of Thermodynamics is used for all system components in order to estimate the temperature of the working fluid (air and water), enthalpy and pressure corresponding to each stage of storage system. These properties depend on the design parameters and finally allow one to evaluate the electrical power consumed by compressor, the electrical power produced by generator coupled to the pneumatic motor, the filling and discharge times, the heat produced during the compression process and the required heat to increase temperature of the compressed air during expansion process.

Compressor
For each stage of compression, the output enthalpy is evaluated as a function of the isentropic efficiency: In Eq. (1), h i c, n i is the inlet specific enthalpy of the i-th stage of compressor estimated at the outlet temperature of intercoolers assumed as constant and equal to 35 • C except for the first stage where it is estimated at ambient temperature T 0 , h c o , u i t is the outlet specific enthalpy of the i-th stage of compressor if the compression process is isentropic, and η IsC is the isentropic efficiency of compressor.
The outlet isentropic temperature of each compressor stage can be determined from: Here π is the pressure ratio which is the same for all stages of compression. Since T The electrical power consumed by the compressor can then be calculated as Where P elc,c (often written asẆ elc,c [43]) is the electrical power consumed by the compressor, η elc,c is the electric efficiency of the compressor assumed to be constant in this study and equal to 98% (upper value of literature ranges from 90% to 98% [44]), η m,c is mechanical efficiency of compressor,ṁ c is the air mass flow rate of compressor and n is the number of compression stages.

Intercoolers
To reduce the work-input required during the compression and prevent the compressor from reaching high temperatures, counter flow air-to-water heat exchangers (called intercoolers and after cooler) cool the compressed air between the stages and after the last stage of the process. By knowing the cooling water mass flow rateṁ cw as well as its inlet temperature, the specific enthalpy of cooling water at intercooler outlet between (i)-th and (i+1)-th stage of compression is computed by an energy balance: Since the inlet temperature of cooling water T in cw,i is known, its specific enthalpy h in cw,i is obtained from Table A-22 in [43] and then its outlet temperature T out cw,i can be obtained by interpolation. The heat stored in the hot water tank (HWt) should be equal to the heat exchanged in the intercoolers (the heat losses are neglected). The total heat transfer rate between air and cooling water in intercoolers during air compression process can be calculated as: The specific enthalpy of cooling water at the inlet of hot water tank is defined as Knowing h in cw,t the final temperature of hot water at the inlet of HWt can be obtained by interpolation from Table A

Pneumatic motor
The expansion process can be regarded as an opposite thermodynamic process of compression. The output specific enthalpy of high and low power stage is respectively, defined as In Eqs. (8) and (9) the outlet isentropic specific enthalpy of each expansion stage h out m,Is,HP and h out m,Is,LP are obtained by interpolation in Table A-22 in [43] knowing S 0 (T out m,Is,HP ) and S 0 (T out m,Is,LP ) given respectively by: Here β HP and β LP are the pressure ratios of high-pressure and low-pressure expansion stages respectively, T in m,HP and T in m,LP are inlet temperature of air in these expansion stages.

of 34
The inlet pressure P i H n P of the high-pressure expansion stage can be calculated knowing the pressure ratios β HP and β LP as follows: The outlet pressure corresponding to each expansion stage is calculated The electric power generated by the electric generator coupled to the pneumatic motor during production phase is evaluated as follows: Here P elc,G is the electric power produced by the generator, η elc,G is the electric efficiency of generator assumed to be constant in this study and equal to 96% (value of literature range 90% to 98% [44]), η mm is the mechanical efficiency of pneumatic motor, and m˙ m is the air mass flow rate of pneumatic motor which is one of the design parameters.

Heater
To eliminate the use of fossil fuels, to prevent the pneumatic motor from reaching low temperatures and to enhance the power production of pneumatic motor, the combustion chambers usually used in conventional CAES system are replaced by counter flow air-to-water heat exchangers, called heater or air preheater. The hot water produced and stored during compression process is used to warm the air up before each expansion stages. With the assumption that HWt process is adiabatic, the inlet specific enthalpy of hot water h i h n w , is known. Then, the specific enthalpy of water at heaters outlet can be respectively expressed as follows: The thermal power required to warm the air up before it is expanded in the pneumatic motor can be computed as follows:Q

Compressed air storage tank
This step is crucial because it allows one to determine the charge and discharge time of the CASt as well as the mean temperature of air at the inlet of heater before the high power expansion stage. We consider CASt as one thermodynamic control volume with the total geometric volume V t . In this study, we assume that the inlet and outlet mass flow rates are constant. The inlet air temperature in the control volume T in c,n is assumed to be equal to that of air exiting the last cooler (35 • C ) and the inlet pressure of compressed air P out c,n is assumed to be equal to that of air exiting the last compression stage. The exit air pressure P out CASt is set to the inlet air pressure of the air motor P in HP and its outlet temperature T out CASt is assumed to be the minimum value of temperature inside the control volume. All these assumptions are summarized in Figure 2, where m, P and T represent respectively air mass, air pressure and air temperature inside the tank andQ CASt is the thermal power lost through the compressed air storage tank walls. Thus, the variations of mass, temperature and pressure of air in the CASt during filling and discharge process are described by the mass conservation equation, the energy conservation equation and the ideal gas equation of state [43].
For the filling process whereṁ out = 0 andṁ in =ṁ c these laws read: And for the discharge process whereṁ in = 0 andṁ out = −ṁ m these laws read: With; In this study, the reference state for enthalpy is 0K and h(T = 0K) = 0 kJ.kg −1 .K −1 . In Eqs. (21) and (22), c p is the specific heat at constant pressure given by (23) [43] which is evaluated during the numerical resolution at the current time step at the temperature corresponding to the previous step time, A in CASt and A out CASt are respectively inlet and outlet cross section of CASt, H t is the height of the CASt and g is the acceleration of the gravity. The heat exchange through the tank wallsQ CASt , is modeled through a quasi-steady process assuming a cylindrical geometry with steel shell structure of the CASt. It is calculated as follows:Q Where H FP/DP , is the heat transfer coefficient between the CASt wall and the air during filling and discharge process (equal to 40 and 45 W.m 2 .K −1 , respectively [45]), S is the area of heat transfer between the air and the CASt wall, T is the temperature of air inside the CASt and T 0 is the environment temperature.
Equations (21) and (22) are solved using the fourth-order Runge-Kutta method presented by Press et al. [46]. As initial conditions of air inside the CASt during filling process we use the environmental temperature and the minimal pressure set to inlet pressure P in HP = P 0 β HP β LP of the high-power expansion stage. The stopping condition relates to the maximum pressure in CASt set to the outlet pressure of last compression stage P out c,n = P 0 × π n . These initial and stopping conditions of air inside the CASt are reversed for the discharge process. At this level of the analysis, the temperature, pressure and enthalpy of each line of storage system to be optimized are known and for each stage of the process. Therefore, the exergy analysis can be carried out with greater flexibility.

Exergy analysis
Exergy is defined as the maximum theoretical work obtainable from an overall system consisting of a system and the environment as the system that comes into equilibrium with the environment [43,47,48]. It can equally be defined as the maximum work that can be obtained from a given form of energy when the reference environment state is defined by the environmental parameters [49]. This second definition is more appropriate to the approach used in this work. In order to carry out the exergy analysis the exergy rate balance should be applied to each component. In steady state, this exergy rate balance for given component with one inlet and one outlet can be expressed as follows: Here E˙x Q is the time rate of exergy transfer associated with heat transfer, E˙x W is the exergy of the work, E˙x D is the time rate of exergy destruction, and E˙x in and E˙x out are the time rate of exergy transfer at inlet and outlet of the considered component respectively. These parameters are defined as follows: Here e x is the specific flow exergy (also known as physical exergy), T b is the temperature of the boundary where heat transfer ( ) occurs. Generally, to evaluate E˙x Q which is associated to exergy loss for a given component it is necessary to know the heat transfer Q˙ across each segment of the boundary and T b . Although it is sometimes possible to calculate Q˙ , the temperature of the boundary is more difficult to obtain and requires experimental measurements. So, an alternative approach that often suffices for modelling is to suppose that the boundary is the outer surface of each component where the temperature corresponds to the ambient temperature taken as the temperature of the exergy reference environment. Thus, the heat transfer occurs at T 0 (T b = T 0 ) and therefore there exists no exergy loss [50]. In this case, the rate of exergy destruction term of Eq. (25) accounts for the exergy destruction owing to friction and the irreversibility of heat transfer within the considered component [50]. Then, using Eqs. (26) to (29), the rate of exergy destruction of all components could be calculated as shown in Table 1.

Storage system performance criteria
We use the efficiency as a mean performance parameter of proposed system. We define the global energy efficiency also known as the Round-Trip Efficiency (RTE) and overall exergy efficiency (η ex ) as follows. The first one is based on the First Law of Thermodynamics that defines the efficiency of an engineering system by the ratio of energy outputs to inputs. For our proposed system, the energy inputs is the electricity used by the compressor to produce compressed air. The energy output is the sum of the electrical energy produced by generator and the part of heat recovered during the compression process that has not been used to reheat the air before its expansion. The second and perhaps the most relevant performance criteria of our proposed system is the exergy efficiency. It is an efficiency based on the Second Law of Thermodynamics and it is defined as the ratio of total exergy outputs to exergy inputs [42,48]. Exergy efficiency is also known as the ratio of the product exergy to the fuel exergy. The fuel exergy is defined by the electrical energy consumed by compressor and the product exergy is the difference between the fuel exergy and the sum of exergy destruction in all components of the storage system and total exergy loss associated with the overall considered system [41,51]. To conclude this part, these efficiencies can be express as follows: Here t FP and t DP are respectively the charge and discharge time, Q˙ h is the total heat transfer rate between air and cooling water in intercoolers during air compression process (see Eq. (6)) and Q˙ Rh the thermal power required to warm the air up before it is expanded in the pneumatic motor (see Eq. (19)).

Definition of objective function
There is a discussion about the definition of global energy efficiency (or RTE) of CAES system. While some, to conform to the First Law of Thermodynamics, define the global energy efficiency by Eq. (30), others suggest taking into account the fact that the electrical energy and heat are different energy forms. For the latter, it would be necessary to convert the heat power of the hot water into its electrical equivalent. For this purpose, they assume a virtual thermal power plant that would use the thermal power of water as heat source [31]. The energy efficiency of this virtual power plant given by Eq. (32) would allow the deduction of the electrical equivalent of heat from hot water (E elc,eq,hw ) using Eq. (33).
Electrical equivalent of the heat power of the hot water(E elc,eq,hw ) Here, η re f is the thermal efficiency of the virtual power plant generally considered equal to that of a reference natural gas power plant (38.2%) [52].
The global energy efficiency of the SS-CAES system would therefore be written as follows; To avoid controversy on good definition of energy efficiency of SS-CAES system, we decided to use as an objective function, the overall exergy efficiency of the system (Eq. (31)). This objective function is subject to thirteen design parameters and seven inequality constraints. The mathematical formulation of the optimization problem is given by: Identify X which maximizes η ex ( X) subject to seven inequality constraints: Here X stands for the solution vector containing the N par = 13 design parameters X = j x 1 , x 2 , ..., x N par and each of them varies in the range of lower and upper bounds x j L , x U .

Constraints
Seven constraints are considered in this optimization procedure: • The minimum power delivered by the storage system during discharge process may not be smaller than 10 4 W (10 kW) so that they can be implemented in a house of the residential sector.
The proposed system needs to be used to store low cost electricity available during off-peak hours and to store electricity from renewable energy sources. Thus, the charge time must not exceed 43200s (12h) and the discharge time should be greater than 7200s (2h).
• The hot water used to reheat air during discharge process is produced during the compression process. Thus, to eliminate the need for infinitely long heat exchangers, the difference between the hot water temperature and the inlet temperature of air of expansion stages must be larger than 5 K.
• Due to mechanical constraints and safety problems, the maximum pressure in compressed air storage tank (CASt) cannot exceed 300 × 10 5 Pa (300 bar).
• The mass of hot water produced during the compression process must be greater than or equal to that required to warm the air up during the discharge process.

Design parameters
As noted in the introduction section, the design parameters selected for use in this study are: number of compression stages (n), compressor pressure ratio (π), volume of air storage tank (V t ), pressure ratio of high-pressure and low-pressure expansion stages of pneumatic motor (β HP , β LP ), inlet temperature of high-pressure and low-pressure expansion stages (T in m,HP , T in m,LP ), compressor and pneumatic motor isentropic efficiency (η IsC , η Is,m ), compressor and pneumatic motor mechanical efficiency (η mC , η mm ), compressor and pneumatic motor mass flow rate (ṁ c ,ṁ m ). The ranges of each parameter have been specified according to the working specifications of each hardware element and summarized in Table 2.

Modified real coded genetic algorithm
Developed by John Holland [53], the genetic algorithm (GA) is an optimization and search technique based on the principles of genetics and natural selection. To perform the optimization, GA produces some random numbers for each design variables that form a population of individuals called initial population, where an individual consists of values of the design variables is a potential solution which maximises the overall exergy efficiency. Based on an analogy with Darwin's laws of natural selection, GA applies to an initial population, the operators of selection, crossover and mutation to allow it to evolve to a new population that is, the next generation. The type of encoding used to represent these design variables are defined by the type of algorithm, thus when the design variables are continuous (as is the case for this study), it is more logical to represent them by floating-point numbers rather than by binary numbers. This is referred to as Real Coded Genetic Algoritm (RCGA) also known as continuous GA [54]. The algorithm has following steps: • Generate the initial population.
• Evaluate the fitness of each individuals of the considered population.
• Select individuals to form the mating pool.
• Select individuals of the mating pool for mating.
• Apply crossover to generate offspring who is individuals of next generation.
• Maintain the diversity in the population by mutation of selected members of the population.
• Terminate the run if the stopping criteria are fulfilled or go back to step 2.
In order to improve the search efficiency and closer simulate the natural selection which is the fundamental principle of genetic algorithm, one modification is introduced in the crossover step. In this step, during the evolution process, two crossover operators are randomly selected to generate the children according to the probability of crossover (Pc). Thus, for a given generation a random number is generated and compared to the probability of crossover (Pc). If this random number is smaller than the probability of crossover, the Simulated Binary crossover (SBX) proposed by Deb and Agrawal [55] is used. Otherwise the Simplex crossover (SPX) developed by Da Ronco and Benini [56] is used. In order to ensure a good exploration of the search space and avoid convergence towards a local optimum value, a newly developed mutation operator named Dynamic Random Mutation (DRM) proposed by Chuang et al. [57] was used. A recent technique of constraint-handling named Violation Constraint-Handling method (VCH) introduced by Chehouri et al. [58] was used in the selection steps during the evolution process. The algorithm configuration of the modified RCGA used in this work is shown in Figure 3.
In this algorithm, normalizing an inequality constraint consists of transforming it to ensure that it admits 1 for maximum as presented in section 4.2. Thus, the normalized constraints used in step 2 in Figure 3, can be expressed as follows: Chehouri et al. [58] define Constraint violation factor (C.V) and number of violation (N.V) that are evaluated here in step 2 for each individual (chromosome) of considered population of possible solutions as follows: N.V = number of violatated constrints 7 The pair-wise feasibility rules used in step 3 of Figure. 3 separate the population into two families; feasible solutions and unfeasible consisting of individuals that violate at least one of seven constraints. The family of feasible solutions is sorted with respect to their fitness value (exergy efficiency) in descending order. The second family is sorted according to these rules: • If two considered chromosomes are infeasible, the best is the one with the lowest Number of Violations (N.V). • If both chromosomes have the same (N.V), the one with the lowest Constraints Violation (C.V) value is the best.
These sorted unfeasible solutions are placed after the sorted feasible solutions and the resulting sorted population is then used in selection of parents to be crossed. The first NKeep individuals of this sorted population are kept to form the mating pool and the rests are discarded and replaced by offspring of parents selected randomly in this mating pool. Since the mating pool is sorted, the tournament selection approach used in this work consists of selecting randomly two chromosomes from the mating pool, the chromosome with the lowest rank becomes a parent. The tournament repeats for every parent needed (twice for SBX and once for SPX). The first individual of mating pool must not be altered by any evolutionary operator (elitism strategy).

Validation of the thermodynamic model
The storage system we propose in this work has not been constructed in reality. It is inspired from an experimental prototype existing in our lab and schematically illustrated in (Figure. 4), due to the poor efficiency 3.4 % obtained by experimentally measuring the total electrical energy produced during the discharge phase (0.45 kW h) and consumed by the compressor during the charge phase (13.12 kW h). In this experimental prototype, a three-stage compressor (cylinders) is used for the production of compressed air. These compression cylinders (a) are separated by intercoolers (b) and the compressed air produced is stored in six storage tanks (c) having a total volume of 300 L. The maximum pressure in storage tanks is set at 180 bar in order to maintain the compression ratios of the three cylinders constant throughout the storage phase. The minimum pressure in storage tanks is limited at 16 bar to ensure good regulation of mass flow rate of expander (d). A fan (f) driven by the electric driving motor (M) of the compressor, stirs the ambient air to cool the compressed air passing through the intercoolers. The rated air mass flow rate of the compressor is 14.4 kg h −1 for the rated power of driven-motor of 4 kW. The pressure and temperature of air at the inlet and outlet of each compression cylinders are acquired. A pneumatic motor (e) coupled to electric generator (G) is used to achieve the expansion process and generate electricity. The inlet pressure and mass flow rate of pneumatic motor are adjusted by expander (d), its rated air mass flow rate is 97.2 kg h −1 and inlet air pressure is 8 bar for a rated output power of 1.2 kW. Because the fan also serves as a flywheel for the drive of the compressor, the ambient air that it stirs also cools the compression cylinders. The measured temperatures cannot fit the model results because of this compression cylinders cooling. That is why, we just use this prototype to validate the model for the filling and the discharge of the tank. The pressure sensor is used for the acquisition of pression inside the storage tanks during charge and discharge process. Figure. 5 shows the comparison of the experimental results and the simulation results. The plot in the left panel shows the variation of pressure inside the storage tank with the time during the filling process; the right shows the same physical quantity during the discharge process. A good agreement between the experimental and the simulation results is observed. The rest of the model (compression, expansion) is validated using the data published by Liu et al. [41] where they thermodynamically analysed a CAES system through an advanced exergetic analysis. We aim to compare the values of the temperature to the different states of their system with those obtained by our model and then to do the same for the specific flow exergy. As shows in Table 3, the differences are quite low (below 2.5 %); one may conclude that our model accurately describes both the filling and the discharge process.
The comparison of the temperature and specific flow exergy predicted by our model with previous results presented in Ref. [41] is illustrated in Table 4. We note that the discrepancy is relatively smaller for the compression process than for the expansion process. This can be explained by the fact that we have considered the combustion gases as an ideal gases with the same properties as the air, therefore it could be concluded that the model is accurate.

Optimization results
In this section, the results of optimization and thermodynamic analysis of optimized system are presented. The simulation code was implemented under Matlab R . In order to understand the impact of RCGA parameters to the optimization procedure, we have distinguished six cases with respect to the population size change and change of the number of maximum generation as shown in Table 5; Here PopSize is the size of population, MaxGens is the maximum number of generations, Pc is the probability of crossover and Pm is the probability of mutation. In all these cases, the probability of crossover and that of mutation are the same and equal respectively to 60 % and 0.5 %. The population size in Case 1 is 50 and the maximum number of generation is 100. For more diversity in initiale population, the population size is doubled in Case 2 for the same maximum number of generation. In the Cases 1 through 5, the population size is held constant at 100 and the evolution time is increased gradually by raising the maximum number of generation from 200 for the Case 3 to 1000 for the Case 5. Finally in Case 6, for the purpose of observing the impact of diversity in the initial population on finding an optimal solution, the population size is doubled with the same maximum number of generation as in the Case 4 (500). The evolutions of the maximum value of the objective function (η ex−Max ) in each generation for these six cases have been plotted as show in Figure 6 where NumGens is the number of generation.
It can been seen in the Case 1 of Figure 6 that, a stable value of 13.85 % is reached for minimum generation number of 70, but this value remains stable for at least 30 generations after the plateau is reached. So, we are not sure whether this stable value is the global optimum of the storage system. One possible explanation of this late convergence may be the lack of diversity of the initial population. Therefore, in Case 2, the population size is doubled with the same maximum number of generation. The representative curve of this second case shows that; the population is gradually converging to the stable objective value of 18.02 % which is the maximum value of the exergy efficiency reached during all the evolution process. Also, this value is reached for minimum generation number of 72 just 28 generation before the maximum number of predefined generation. Therefore, it is not certain that this maximum exergy efficiency of 18.02 % is the optimum solution because the evolution process may not be completed. So, in the Case 3, the maximum generation number is doubled with the same population size as in the Case 2 in order to maintain good diversity in the population and a sufficient evolution time to expect a convergence towards the overall optimum. The result of this case shows that the population is converging to a stable value of exergy efficiency. But it is not sure whether this maximum stable value of 19.88 % is the global optimum since, as can be seen, by doubling the evolution times (NumGens = 200), the maximum value of exergy efficiency is increased by two percentage points (from 18.02 % to 19.88 %). That is why, in the fourth case, the evolution time has increased fivefold (NumGens = 500) while keeping the same diversity in the initial population as in the third case (PopSize = 100).The result of this fourth case shows that the population is converging to a stable value of maximum exergy efficiency that increases by five percentage points to 24.87% in case 4, compared with 19.88% in case 3. This maximum value of 24.87 % remains constant over the last 210 generations during evolution process. It is likely that the algorithm has converged to the global optimum solution.
To make sure of that, we wanted to see whether increasing the evolution time could be able to have an impact on the maximum value of the exergy efficiency. Therefore, in the fifth case, the maximum number of generation is doubled (NumGens = 1000) with the same population size as in the previous case. As shown in the representative curve of the evolution process of this fifth case, the population converges to a stable value but this stable value of 23.40 % is unfortunatly lower than that obtened in Case 4. Since an evolution time beyond 500 generations does not have a relevant impact on the maximum value of exergy efficiency, in the last case, the population size used in the fourth case is doubled (for more diversity in the initial population) with the same maximum number of generation. The stable value of maximum exergy efficiency of 24.81 % obtained at the end of evolution process was indeed close to that reached in the Case 4. Nevertheless, this value remains lower than 24.87 % of Case 4 which is certainly the global optimum of this optimization problem.
To return to Case 4, the observation of the final population at the end of the evolution process shows that, all the chromosomes (sets of design parameters) are identical. This means that almost all the chromosomes of the population have converged to the optimum solution. We have then deduced the set of optimal values of influencing parameters that maximizes the global exergy efficiency of the SS-CAES system without fossil fuel used. These optimal parameters are given in Table 6.
Using these optimal parameters, the thermodynamic properties of each point in the optimized system ( Figure. 7) are shown in Table 7. The performance indicator of the optimized system is calculated and shown in Table 8. As it can be seen in Table 8, under optimal design conditions, the compressors take 8.65 h to fill the CASt while the discharge time of the CASt is equal to 2.02 h. The electric power produced by the generator is 12.77 kW which is well above the electrical energy requirements of an average household (of about 4 individuals). The optimized system also produces 4.64 t of hot water at a temperature of about 30 • C during each cycle of charge/discharge which can be used for heating purposes in a house. The RTE of the optimized system is equal to 79.07 % and its exergy efficiency is 24.87 %. By calculating the contribution of thermal energy to RTE (CTEtoRTE, Eq. (44)), we can see that it represents about 69 % of RTE. Unfortunately, the low-temperature hot water containing this thermal energy does not necessarily have the same value as the electrical energy produced during the discharge process. This is one more reason to use the overall exergy efficiency as a performance evaluation criterion of SS-CAES system proposed in this work.
In order to identify the locations and magnitudes of storage process inefficiencies, the total exergy destruction of the charge/discharge cycle for each component as well as their exergy destruction ratio are listed in Table 9.
As shown in Table 9, the compression stage has the largest exergy destruction, followed by the intercoolers. This can be explained by the fact that, the temperature difference between the inlet and the outlet of each stage is large enough so that the heat transfer to the compression stages walls is no longer negligible compared to the enthalpy change. Another explanation could be given by the advanced exergetic analysis [59,60]. It would certainly indicate that the unavoidable exogenous part of the exergy destruction within these components is important. Indeed, as shown in Table 6, the design of compressor is supposed to be perfect because of the higher values of efficiency (isentropic and mehanical). Thus, the unavoidable endogenous part of the ecergy destruction within the compressor is low compared with the unavoidable exogenous part.

Distribution of design parameters
The behavior between lower and upper bounds (dotted lines in Figure 8 and Figure 9) of each optimal design parameters during the evolution is analysed in this section. The results of this are shown in Figure 8 and Figure 9).
According to this distribution, it can be seen in Figure 9-g that the pneumatic motor mass flow rate reaches its minimum value. This means that the decreased value of this design variable could improve the global exergy efficiency of system. In order to reduce the number of variables and thus to speed up the optimization script, the lower born value of the pneumatic motor mass flow rate can be selected and set as a constant.
As also showns in Figure 9-b to Figure 9-e concerning respectively the distribution in the population during the evolution process, of the values taken by the isentropic efficiency of compressor, by the isentropic efficiency of pneumatic motor, by the mechanical efficiency of compressor and by the mechanical efficiency of pneumatic motor; during the evolution process, these four design parameters only take almost exclusively their maximum value. This means that, the increased value of these designs variables could improve the global exergy efficiency of system. Similarly, the upper born value of these four designs parameters can be selected and set as a constant.
Other design variables have the scattered distribution which mean that, their variations could have some significant impact on the system performance. To have an idea of these impacts on both global exergy efficiency and number of violated constraints at the optimal point, the sensitivity analyses have been made.

Effect of variation of the design variables value on the system efficiency (global exergy efficiency) and on the number of violated constraints
In this section we study the effect of variation of the design variables value on the system efficiency and on the number of violated constraints. For this purpose, we vary the value of each design variable in its allowable range specified in Table 2. 6.3.1. Number of compression stages and the compression ratio values Figure 10-a shows the decrease of the global exergy efficiency with the increasing number of compression stages at fixed optimal compression ratio. The same trend is observed in Figure 10-b for the compression ratio at fixed optimal number of compression stages but, unfortunately, the compression ratio values which provide a maximum exergy efficiency violate some imposed constraints. That is why, as shown in Figure 10-c&d, the optimum values of these two design parameters do not violate any constraint.
Thus, for a SS-CAES system, it is not necessary to use more than three compression stages. Furthermore, the compression ratio value of these compression stages should not exceed four. As can be seen in Figure 11-a, the increase in volume of the air storage tank results in an increase of the global exergy efficiency. But, before the optimal value is reached, all other values violate one constraint as shown in Figure 11-c. This means that if it had been possible to have a larger tank, the system efficiency would have been greater than the optimal value obtained. However, for a given SS-CAES system, a trade-off has to be found between the constraints of space, of charging time, of discharge time and even of the cost of purchasing the storage tanks. Therefore, the optimal value of the air storage tank volume is closely linked to the imposed constraints of the system. 6.3.3. Pressure ratio values of high-pressure and the low-pressure expansion stages of pneumatic motor As it is shown in Figure 11-b&e, the exergy efficiency increases almost linearly with both the pressure ratio of high-pressure and the low-pressure expansion stages of pneumatic motor. But, when the value of pressure ratio of high-pressure expansion stage is greater than 9, the global exergy efficiency begins to decrease (Figure 11-b). This can be explained by the low temperature taken by the air at the exit of such expansion stage. In fact, the temperature of hot water produced during compression process is not high enough. Thus, the pre-heating of air between the expansion stages will no longer be enough to improve the expansion work of the low-pressure stage. Finally, imposed constraints (the minimum power delivered by the storage system during discharge process and the minimum difference between the hot water temperature and the inlet temperature of air of expansion stages) justify the optimal values returned by the optimization algorithm (Figure 11-d&f).
For a SS-CAES system using pneumatic motor as expansion system, it is possible to increase its efficiency by using the pressure ratio values of 9 and 6 for the high-pressure and the low-pressure expansion stages, respectively. However, it should be considered in this case to use a water heating system (like solar water heaters) outside of the SS-CAES system. This water heating system would be used to increase the temperature of the hot water produced during the compression process. As shown in Figure 12-a&b, the exergy efficiency increases linearly with both the inlet temperature of high-pressure and low-pressure expansion stages of pneumatic motor. Moreover, the optimal values of these two design parameters returned by the genetic algorithm are the last ones which do not violate any of the imposed constraints (Figure 12-c&d). If the temperature of the hot water had been higher, the optimal values of these two variables would have been higher too.
For a given SS-CAES, only thermomechanical stresses should limit the supply temperature of the expansion stages. Thus,when designing a SS-CAES using pneumatic air motor, consideration should be given if necessary to an external water heating systems. In that case, during the preheating of the supplying air of the expansion stages during the discharge process, it is possible to move closer to the maximum permissible temperature at the inlet of these stages. 6.3.5. Compressor and pneumatic motor isentropic efficiency values, of the compressor and pneumatic motor mechanical efficiency values As it is shown in Figure 13-a&b and in Figure 14-a&b respectively, the increase in the isentropic and mechanical efficiency of compressor and pneumatic motor leads to a linear increase of the exergy efficiency. Furthermore, it is important to note that this increase is done without any violation of constraints ( Figure 13-c&d and Figure 14-c&d). Therefore, the maximum values of these four design parameters can be selected and set as constants in the cost function for optimization.
We should remember, however, that these results were predictable since for an energy system having turbomachines, the higher is the isentropic efficiency (or mechanical efficiency), the higher will be the system efficiency. Nevertheless, by including these four parameters among the optimization parameters, we are assessing the smooth functioning of the modified RCGA proposed in this work.   Figure 15-a shows that the global exergy efficiency increases with the compressor mass flow rate. In addition, the single value that does not violate any constraint is that returned by the optimization algorithm (Figure 15-c). Unfortunately, its scattered distribution (Figure 9-f) does not allow one to set it as constant during the optimization procedure. Nevertheless, when designing a SS-CAES system, it would be appropriate to give preference to a compressor with a high mass flow rate.
For the pneumatic motor mass flow rate, its increase leads to a decrease in global exergy efficiency (Figure 15-b) and as in the case of compressor mass flow rate, the optimum value of this parameter returned by optimization algorithm is the only one that does not violate any constraint (Figure 15-d).
In contrast to the scattered distribution of the compressor mass flow rate, a homogeneous convergence towards the minimum value can be noted in Figure 9-g. Therefore, this minimum value of the pneumatic motor mass flow rate can be selected and set as a constant in the optimization algorithm. When designing a SS-CAES system, it would be appropriate to give preference to a pneumatic motor with a low mass flow rate.

Conclusion
This paper has presented the energy and exergy analyses of a SS-CAES system without fossil fuel used. The resulting thermodynamic model is fed to a modified RCGA to identify the optimal values of thirteen design parameter of the proposed storage system. The modified RCGA that has been clearly presented and tested to verify its stability and robustness. The results of the optimization indicate that for a maximum efficiency of a SS-CAES system using multistage reciprocating compressor and two stages pneumatic air motor: • The number of compression stage should be less than three.
• The compressor pression ratio of each compression stages should be less than four.
• The maximum value of air storage tank volume allowed by the constraints of spaces, of cost, of charge time and of discharge time should be preferred. • The pressure ratio together with the temperature of the supplying air of the expansion stages are highly dependent on the temperature of the hot water produced during compression process. However, the use of the maximum technologically acceptable values of these four parameters improve the efficiency of the storage systems. • A low mass flow of the pneumatic air motor coupled with a high mass flow rate of the compressor improves the efficiency of the storage system.
The presented modified RCGA can be used for optimizing all scale of CAES system under all kinds of constraints if the optimization problem is well formulated.