Optimal Operation of Combined Photovoltaic Electrolyzer Systems

Abstract: In this study, the design and simulation of a combination of a photovoltaic (PV) array with an alkaline electrolyzer is performed to maximize the production of hydrogen as a reliable power resource. Detailed electrical model of PV system, as long as thermal and electrochemical model of electrolyzer is used. Since an electrolyzer is a non-linear load, its coupling with PV systems to get the best power transfer is very important. Solar irradiation calculations were done for the region of Miami (FL, USA), giving an optimal surface slope of 25.7 ̋ for the PV array. The size of the PV array is optimized, considering maximum hydrogen production and minimum excess power production in a diurnal operation of a system using the imperialistic competitive algorithm (ICA). The results show that for a 10 kW alkaline electrolyzer, a PV array with a nominal power of 12.3 kW The results show that 12.3 kW photvoltaic system can be utilized for supplying a 10 kW electrolyzer. Hydrogen production and Faraday efficiency of the system are 697.21 mol and 0.3905 mol, respectively.


Introduction
The demand for environmentally-friendly renewable energy resources has been growing recently because of a decreasing demand for fossil fuels, and an increasing demand for power [1,2].Among the various kinds of renewable power resources, wind and solar power are the best-known.The power sources of these systems are abundant, free, and environmentally friendly.Photovoltaic (PV) systems can be installed to generate electricity anywhere there is a suitable amount of solar irradiation, however, the irradiation is only available during the day.Additionally, this daily irradiation fluctuates with non-linear behavior because of unstable climatic conditions.Therefore, an energy storage system is needed to increase the reliability of solar systems.Batteries are good choices, however, their use for bulk power production during extended time periods has some issues.First, during operation, they have significant amounts of hourly power leakage, which makes them useless for long-term applications.Batteries are also not capable of storing large amounts of energy in commercial applications.There are other novel techniques, such as superconductivity, which have been developed in order to store and deliver power instantaneously [3][4][5][6][7][8].Hydrogen is another suitable option because of its ability to work as a reliable fuel for almost every application, especially in transportation devices.Furthermore, hydrogen can be converted to electricity, and is able to heat more efficiently than fossil fuels [9][10][11][12][13].Combining hydrogen fuel with other sources of energy makes the system more reliable, secure, and flexible with respect to different energy management techniques.One of the most environmentally-friendly and technically-mature procedures to produce hydrogen from electricity is by using an electrolyzer, which is able to operate in a large range of capacity with high efficiency, combined with various energy production sources [14,15].The connection of a PV system Energies 2016, 9, 332 2 of 12 to an electrolyzer is possible, either directly or by powered electronic devices [16].However, with an optimal sizing of the system, the effects of lack of power of electronic devices can be minimized, which is beneficial in an economic point of view, making hydrogen fuel cost competitive with respect to conventional, commonly-used fuels.
Various studies have been done regarding the commercial, technological, techno-economical, and environmental aspects of this issue.Moreover, studies of different conditions of energy management procedures, with various scenarios about how to combine them to get the most efficient operation, have been performed [17].
In this study, generation of hydrogen with a stand-alone PV power system is evaluated.The model of components of the PV array, electrolyzer, and storage tank are described in details.The optimal sizing of the PV array is simulated and discussed.As the cost of the components of this system varies over time, another index for optimization is used, which is independent of time, but can also serve to describe the cost.The introduced index is based on maximum hydrogen production with a minimum excess of power production.

Analysis and Modeling
Detailed models of the PV system and electrolyzer are given, and the individual characteristics of the various components are discussed.The direct coupling of the components must be optimized in order to get the most of the produced power to the electrolyzer.A PV array generates the power needed for the electrolysis process in an alkaline electrolyzer.The PV system should be optimally sized in order to get as much hydrogen as possible from a 10-kW electrolyzer.

Photovoltaic System
The circuit diagram of a one-diode model of a PV cell is shown in Figure 1.The model can be described based on Kirchhoff's current law [18,19].
Energies 2016, 9, 332 2 of 11 devices can be minimized, which is beneficial in an economic point of view, making hydrogen fuel cost competitive with respect to conventional, commonly-used fuels.
Various studies have been done regarding the commercial, technological, techno-economical, and environmental aspects of this issue.Moreover, studies of different conditions of energy management procedures, with various scenarios about how to combine them to get the most efficient operation, have been performed [17].
In this study, generation of hydrogen with a stand-alone PV power system is evaluated.The model of components of the PV array, electrolyzer, and storage tank are described in details.The optimal sizing of the PV array is simulated and discussed.As the cost of the components of this system varies over time, another index for optimization is used, which is independent of time, but can also serve to describe the cost.The introduced index is based on maximum hydrogen production with a minimum excess of power production.

Analysis and Modeling
Detailed models of the PV system and electrolyzer are given, and the individual characteristics of the various components are discussed.The direct coupling of the components must be optimized in order to get the most of the produced power to the electrolyzer.A PV array generates the power needed for the electrolysis process in an alkaline electrolyzer.The PV system should be optimally sized in order to get as much hydrogen as possible from a 10-kW electrolyzer.

Photovoltaic System
The circuit diagram of a one-diode model of a PV cell is shown in Figure 1.The model can be described based on Kirchhoff's current law [18,19].As can be seen, the open-circuit voltages pattern is non-linear and the short-circuit current is directly proportional to the irradiance.During darkness, assuming that the reverse saturation current is negligible, the solar cell works as a diode, which is off, and produces neither voltage nor current.The determination of I-V characteristics of the cell by diode is given by Equation (1): The characteristics of the PV system are tabulated in Table 1.
As derived from Equation (1), the difference between the generated current IL, from the shunt leakage current ISH, and diode current ID, gives the output current.ID is defined by: As can be seen, the open-circuit voltages pattern is non-linear and the short-circuit current is directly proportional to the irradiance.During darkness, assuming that the reverse saturation current is negligible, the solar cell works as a diode, which is off, and produces neither voltage nor current.The determination of I-V characteristics of the cell by diode is given by Equation (1): The characteristics of the PV system are tabulated in Table 1.As derived from Equation (1), the difference between the generated current I L , from the shunt leakage current I SH , and diode current I D , gives the output current.I D is defined by: where I 0 is the saturation current of a diode, q is the electronic charge (c), V is the solar cell terminal voltage, R S is the cell series resistance, A is the ideal factor of PV cell, k is the Boltzmann constant.T c is the cells' working temperature (K).The shunt leakage current can be obtained by: Shunt resistance is generated because of manufacturing defects causing power losses and therefore, the I-V curve of the PV is affected by the shunt resistance.It should be mentioned that at low voltages the effect of this resistance is intensified.

Electrolyzer
An electrolyzer is a device that can split hydrogen and oxygen, and form the water using electricity injected into it.The electrolyzer used in this study is an alkaline electrolyzer.The electrolyzer model is based on a combination of fundamental thermodynamics, heat transfer theory and empirical electrochemical relationships [20][21][22][23].

Electrochemical Model
In the operation of an electrolyser, the voltage of the cell is dependent on the three reversible, activation and ohmic overvoltages, as given by [20]: The reversible voltage is the minimum electric voltage that must be delivered to the electrolyzer because of its chemical characteristics, based on the Gibbs equation as follows: where ∆G is the Gibbs energy which is 237.2 kJ/kmol, z is the number of electrons involved in the reaction which is 2, and F is the Faraday constant equal to 96,485 C/mol.Therefore, the reversible voltage for the alkaline electrolyzer is 1.229 V. Also, the activation and ohmic overvoltages can be determined by: where r i is the ohmic resistance of the electrolyte (Ωm 2 ), s and t i are electrode response overvoltage coefficients which are temperature dependent (V), (A ´1m 2˝C ) [20].The basic form of the U-I characteristic for a given temperature used in this study is expressed by: where r is the ohmic resistance of electrolyte (Ωm 2 ), s and t are overvoltage coefficients which are temperature dependent (V), (A ´1m 2˝C ).The parameters used in this study are given in Table 2.
Figure 2 shows the U-I characteristics of a typical alkaline electrolyzer at temperatures of 20 ˝C and 80 ˝C.The difference between the curves is caused because of temperature dependence of the overvoltages, so when the temperature rises, the current increases drastically.However, the parasitic bulbs of the can affect the current in very high temperatures [24][25][26].
Energies 2016, 9, 332 4 of 11 Figure 2 shows the U-I characteristics of a typical alkaline electrolyzer at temperatures of 20 °C and 80 °C.The difference between the curves is caused because of temperature dependence of the overvoltages, so when the temperature rises, the current increases drastically.However, the parasitic bulbs of the can affect the current in very high temperatures [24][25][26].The actual amount of hydrogen production is less than the theoretical one due to Faraday efficiency which is the ratio of actual hydrogen production to theoretical amount.An empirical expression to depict the Faraday efficiency is given by: where f1 (mA 2 cm −4 ), and f2 are parameters in the Faraday efficiency calculation.The Faraday efficiency parameters at 60 °C are given in Table 3.
The hydrogen production rate for several cells connected in series is expressed by: where nc is the number of cells in the series.

Thermal Model
The operating temperature of an electrolyzer affects its performance.The variation of the temperature is caused by the difference in generated heat, less heat loss and cooling.In other words, as experimentally shown in [27,28], the cell performance increases with a corresponding increase in temperature as a result of reduced reversible potential, kinetic loss and ohmic loss.The overall thermal energy balance is expressed by: The actual amount of hydrogen production is less than the theoretical one due to Faraday efficiency which is the ratio of actual hydrogen production to theoretical amount.An empirical expression to depict the Faraday efficiency is given by: where f 1 (mA 2 cm ´4), and f 2 are parameters in the Faraday efficiency calculation.The Faraday efficiency parameters at 60 ˝C are given in Table 3.

´4
The hydrogen production rate for several cells connected in series is expressed by: .
where n c is the number of cells in the series.

Thermal Model
The operating temperature of an electrolyzer affects its performance.The variation of the temperature is caused by the difference in generated heat, less heat loss and cooling.In other words, Energies 2016, 9, 332 5 of 12 as experimentally shown in [27,28], the cell performance increases with a corresponding increase in temperature as a result of reduced reversible potential, kinetic loss and ohmic loss.The overall thermal energy balance is expressed by:

Storage
The storage of hydrogen is greatly affected by the choice of storage conditions, which include but are not limited to, storage pressure, the attainable volumetric density under that pressure, and the amount of thermal losses that accompany the mode of storage [29].Physical hydrogen storage is a common technique to store hydrogen using tanks.The model of a storage system which calculates tank pressure using the flow rate of produced hydrogen is expressed by [30,31]: where M H2 is the molar mass of hydrogen (kg¨kmol ´1), P b is the pressure of tank (Pa), P bi is the initial pressure of the storage tank (Pa), R is the universal gas constant (J¨kmol ´1 K ´1), T b is the operating temperature (K), V b is the volume of the tank (m 3 ), z is the compressibility factor as a function of the temperature of the storage tank.

Imperialist Competitive Algorithm
The primary objective of this paper to estimate the optimal size of the PV subsystem, taking into account the efficiency of the electrolyzer as well, which are most directly determined by maximal hydrogen production and minimal excess power generation.Most recently, meta-heuristic-based optimization algorithms inspired by biological phenomena and natural occurrences have been proven to perform much better in comparison with other widely available algorithms [32][33][34].One immensely useful addition to this arsenal is the imperialistic competitive algorithm (ICA), which does not require the function gradient (or slope) during its optimization process.While GA views the problem as a biological evolution problem, the ICA views the problem as a socio-economic evolution specific to human civilization, which was witnessed at the dawn of industrial revolution and consequential imperialistic gains.
Having a good search capacity in addition to a powerful exploitation strategy renders this method very advantageous in identifying the global optimal solution [35]. Figure 3 shows the generic flowchart of the ICA, wherein the optimization commences by generating a candidate set of random solutions in the search space of the defined problem.

Storage
The storage of hydrogen is greatly affected by the choice of storage conditions, which include but are not limited to, storage pressure, the attainable volumetric density under that pressure, and the amount of thermal losses that accompany the mode of storage [29].Physical hydrogen storage is a common technique to store hydrogen using tanks.The model of a storage system which calculates tank pressure using the flow rate of produced hydrogen is expressed by [30,31]: where MH2 is the molar mass of hydrogen (kg•kmol −1 ), Pb is the pressure of tank (Pa), Pbi is the initial pressure of the storage tank (Pa), R is the universal gas constant (J•kmol −1 K −1 ), Tb is the operating temperature (K), Vb is the volume of the tank (m 3 ), z is the compressibility factor as a function of the temperature of the storage tank.

Imperialist Competitive Algorithm
The primary objective of this paper to estimate the optimal size of the PV subsystem, taking into account the efficiency of the electrolyzer as well, which are most directly determined by maximal hydrogen production and minimal excess power generation.Most recently, meta-heuristic-based optimization algorithms inspired by biological phenomena and natural occurrences have been proven to perform much better in comparison with other widely available algorithms [32][33][34].One immensely useful addition to this arsenal is the imperialistic competitive algorithm (ICA), which does not require the function gradient (or slope) during its optimization process.While GA views the problem as a biological evolution problem, the ICA views the problem as a socio-economic evolution specific to human civilization, which was witnessed at the dawn of industrial revolution and consequential imperialistic gains.
Having a good search capacity in addition to a powerful exploitation strategy renders this method very advantageous in identifying the global optimal solution [35]. Figure 3 shows the generic flowchart of the ICA, wherein the optimization commences by generating a candidate set of random solutions in the search space of the defined problem.This set of random points is called "countries", which are essentially analogous to the "chromosomes" in GA and "particles" in PSO [35].Like every optimization problem, this is fundamentally built on a cost function as its primary objective as well.However, here, the cost function is contextually associated with the "imperialistic power of a country".Since minimizing the cost is the goal of the objective function, the optimization works towards the "countries" with lower This set of random points is called "countries", which are essentially analogous to the "chromosomes" in GA and "particles" in PSO [35].Like every optimization problem, this is fundamentally built on a cost function as its primary objective as well.However, here, the cost function is contextually associated with the "imperialistic power of a country".Since minimizing the cost is the goal of the objective function, the optimization works towards the "countries" with lower cost (and hence greater power) "acquiring" or "annexing" the "countries" with higher cost (and hence lower power, called "colonies") into their space, called "empires" [37,38].As can be seen from the flowchart, the ICA has two major operators: assimilation and revolution, both of which are antagonistic yet complementary to each other.While assimilation deals with the process of "annexation" of "colonies" to "countries", thereby progressing the solution closer to the intended objective, revolution introduces sudden randomness in the process by altering the positions of some "countries" in the search space.Hence, the continuous operations of assimilation and revolution might transform a "colony" into a "country", enabling it to acquire an entire "empire".The optimization can also be viewed as a game wherein there is a constant struggle between the different "countries" and "colonies" and within the "countries" as well, to attain the greatest "empire", that is, the least cost.Like other optimization strategies in this area, ICA is also iterative in nature.After a series of assimilation, revolution and games (also called competitions), the algorithm reaches its stop condition when the objective is satisfied.Provided below in Table 4 is a brief pseudocode of the generic ICA [39], following which, in Table 5, the pseudocode for the ICA specific to the proposed method is provided.This forms the setup for further competition to happen between the different imperialists wherein each imperialist has the goal of maximizing the power by adding as many colonies as possible.Thus in the process, those empires that lose all of their colonies to others collapse and are weeded out of the optimization problem.In the final state, the colonies and imperialists will all have the same power.

Solar Calculation
Total hourly radiation on a tilted surface is given by: where I b and I d are the beam and diffuse irradiations, respectively.R b is the geometric factor which is the ratio of beam radiation on a tilted surface to that on a horizontal surface at any time.Further, F is the view factor in which the subscripts "c-s", "c-hz", and "c-g" designate cell to sky, cell to horizon, and cell to ground, respectively.The integration of total hourly radiation for one year can give the optimal surface angle.Figure 4 shows the annual radiation for Miami (FL, USA).As it can be seen, the maximum energy is absorbed when the angle is 25.7 ˝which is equal to latitude of the region.Also, the surface azimuth angle should be zero, making the tilts always directed to south.
Energies 2016, 9, 332 7 of 11 where Ib and Id are the beam and diffuse irradiations, respectively.Rb is the geometric factor which is the ratio of beam radiation on a tilted surface to that on a horizontal surface at any time.Further, F is the view factor in which the subscripts "c-s", "c-hz", and "c-g" designate cell to sky, cell to horizon, and cell to ground, respectively.The integration of total hourly radiation for one year can give the optimal surface angle.Figure 4 shows the annual radiation for Miami (FL, USA).As it can be seen, the maximum energy is absorbed when the angle is 25.7° which is equal to latitude of the region.Also, the surface azimuth angle should be zero, making the tilts always directed to south.

Hybrid System Simulation
To obtain the maximum amount of hydrogen from the source, the PV and the electrolyzer are combined to yield the hybrid system being proposed.In this section, the simulation setup for this hybrid system is explained with appropriate results and discussions in the subsequent section.The irradiance computed using Equation ( 13) represents the starting point for the simulations.Derivation of the electrical working point of this entire system is done by equating the electrical characteristics of the PV with those of the electrolyzer, as elaborated by Equations ( 2) and (8).As described earlier, the electrolyzer's operation is dependent to its temperature, which is depicted by Equation ( 6) and (9).Section 3 details the optimization procedure, wherein the step of resolving the colonies also computes the PV power.Through this, the cost function is generated following which the Simulink model runs to obtain results in the step where total cost of the empires is computed.The remainder of the steps of iteration are as expressed in the aforementioned section.

Results and Discussion
According to the mentioned models, the coupling of a PV array and an electrolyzer is simulated in this study.The simulation is based on diurnal solar calculation for the region of Miami which is shown in Figure 5.The optimal size of the system is achieved through ICA [23].The various costs of the system components are not included in optimization process.The dependency of the cost of the system to time makes this factor inappropriate to be considered as cost function.Another index to optimization is introduced, based on maximum hydrogen production and storage through minimum average excess power production.As minimum excess power production tends to lower the dimensions of the system, its cost is indirectly accounted for in the optimization.In this case, it gives 12.3 kW for the size of PV array for a 10 kW advanced alkaline electrolyzer.The results of the simulation are given in Table 6.This condition produces 697 mol of hydrogen in which the power generation of PV is 2.94 kW.For a PV panel of 85 W, 145 panels are needed.The results are almost the same in other profiles of irradiation despite the simulation being performed for a sunny day.

Hybrid System Simulation
To obtain the maximum amount of hydrogen from the source, the PV and the electrolyzer are combined to yield the hybrid system being proposed.In this section, the simulation setup for this hybrid system is explained with appropriate results and discussions in the subsequent section.The irradiance computed using Equation ( 13) represents the starting point for the simulations.Derivation of the electrical working point of this entire system is done by equating the electrical characteristics of the PV with those of the electrolyzer, as elaborated by Equations ( 2) and (8).As described earlier, the electrolyzer's operation is dependent to its temperature, which is depicted by Equation ( 6) and ( 9).Section 3 details the optimization procedure, wherein the step of resolving the colonies also computes the PV power.Through this, the cost function is generated following which the Simulink model runs to obtain results in the step where total cost of the empires is computed.The remainder of the steps of iteration are as expressed in the aforementioned section.

Results and Discussion
According to the mentioned models, the coupling of a PV array and an electrolyzer is simulated in this study.The simulation is based on diurnal solar calculation for the region of Miami which is shown in Figure 5.The optimal size of the system is achieved through ICA [23].The various costs of the system components are not included in optimization process.The dependency of the cost of the system to time makes this factor inappropriate to be considered as cost function.Another index to optimization is introduced, based on maximum hydrogen production and storage through minimum average excess power production.As minimum excess power production tends to lower the dimensions of the system, its cost is indirectly accounted for in the optimization.In this case, it gives 12.3 kW for the size of PV array for a 10 kW advanced alkaline electrolyzer.The results of the simulation are given in Table 6.This condition produces 697 mol of hydrogen in which the power generation of PV is 2.94 kW.For a PV panel of 85 W, 145 panels are needed.The results are almost the same in other profiles of irradiation despite the simulation being performed for a sunny day.Power loss of the system is shown in Figure 6.As it is shown, excess power is produced at the middle of the day, when the irradiation is maximum.The average daily unused power is 0.029 kW, which spends 2500 kJ of energy for that day.The transient switching power of the system cuts off the power flow from the electrolyzer sub-system, causing spikes in the excess power as can be observed.The optimization based on minimizing excess power, means minimizing the dimensions of the system which indirectly reduces the cost.Instead of evaluating the cost in objective function which is variable during time, we minimize the excess power which is independent of time.It should be mentioned that this factor is different from the transfer loss, which is the difference between the actual power point and the maximum power point.
Faraday efficiency directly affects hydrogen production, which is shown in Figure 7.The average Faraday efficiency is 0.52 below the maximum possible value.However, when the PV starts to generate electricity, the efficiency immediately raises to its maximum value.The effect of radiation to faraday efficiency is negligible, so during the operation of the electrolyzer, the efficiency is almost fixed.Power loss of the system is shown in Figure 6.As it is shown, excess power is produced at the middle of the day, when the irradiation is maximum.The average daily unused power is 0.029 kW, which spends 2500 kJ of energy for that day.The transient switching power of the system cuts off the power flow from the electrolyzer sub-system, causing spikes in the excess power as can be observed.The optimization based on minimizing excess power, means minimizing the dimensions of the system which indirectly reduces the cost.Power loss of the system is shown in Figure 6.As it is shown, excess power is produced at the middle of the day, when the irradiation is maximum.The average daily unused power is 0.029 kW, which spends 2500 kJ of energy for that day.The transient switching power of the system cuts off the power flow from the electrolyzer sub-system, causing spikes in the excess power as can be observed.The optimization based on minimizing excess power, means minimizing the dimensions of the system which indirectly reduces the cost.Instead of evaluating the cost in objective function which is variable during time, we minimize the excess power which is independent of time.It should be mentioned that this factor is different from the transfer loss, which is the difference between the actual power point and the maximum power point.
Faraday efficiency directly affects hydrogen production, which is shown in Figure 7.The average Faraday efficiency is 0.52 below the maximum possible value.However, when the PV starts to generate electricity, the efficiency immediately raises to its maximum value.The effect of radiation to faraday efficiency is negligible, so during the operation of the electrolyzer, the efficiency is almost fixed.Instead of evaluating the cost in objective function which is variable during time, we minimize the excess power which is independent of time.It should be mentioned that this factor is different from the transfer loss, which is the difference between the actual power point and the maximum power point.
Faraday efficiency directly affects hydrogen production, which is shown in Figure 7.The average Faraday efficiency is 0.52 below the maximum possible value.However, when the PV starts to generate electricity, the efficiency immediately raises to its maximum value.The effect of radiation to faraday efficiency is negligible, so during the operation of the electrolyzer, the efficiency is almost fixed.
the excess power which is independent of time.It should be mentioned that this factor is different from the transfer loss, which is the difference between the actual power point and the maximum power point.
Faraday efficiency directly affects hydrogen production, which is shown in Figure 7.The average Faraday efficiency is 0.52 below the maximum possible value.However, when the PV starts to generate electricity, the efficiency immediately raises to its maximum value.The effect of radiation to faraday efficiency is negligible, so during the operation of the electrolyzer, the efficiency is almost fixed.The system voltage is another factor which is shown in Figure 8.The proposed system operates at voltages around 35 V.The peak voltage point for the system is 37 V which can produce the maximum amount of hydrogen.As we have not used power electronics devices in our system to control the voltage and make the system work at the maximum power point of PV output, the voltage is set by the connecting point of the U-I characteristics of PV system and electrolyzer.
Energies 2016, 9, 332 9 of 11 The system voltage is another factor which is shown in Figure 8.The proposed system operates at voltages around 35 V.The peak voltage point for the system is 37 V which can produce the maximum amount of hydrogen.As we have not used power electronics devices in our system to control the voltage and make the system work at the maximum power point of PV output, the voltage is set by the connecting point of the U-I characteristics of PV system and electrolyzer.At the end of the operation of the PV/hydrogen system, the produced hydrogen is stored in storage tanks.The storage pressure of a tank for daily operation is shown in Figure 9. Maximum pressure for the storage tank with the capacity of 0.25 m 3 is 3.62 MPa.Therefore, a storage tank with a capacity higher than the given pressure is needed.The stored hydrogen is portable and can be used as a fuel for the input hydrogen of fuel cells.In this study, it is assumed that the pressurizing power of the storage tanks is derived from another source.

Conclusions
The operation of a combined PV-electrolyzer system to produce hydrogen as a fuel is designed and simulated.The solar radiation for one day is used for solar calculations.The maximum hydrogen production and minimum average excess power are considered as the optimization indices, which are 697 mol and 0.029 kW, respectively.This condition occurs when the tilts are in their appropriate location with an angle of 38°.The ratio of the PV array's average power production to its maximum At the end of the operation of the PV/hydrogen system, the produced hydrogen is stored in storage tanks.The storage pressure of a tank for daily operation is shown in Figure 9. Maximum pressure for the storage tank with the capacity of 0.25 m 3 is 3.62 MPa.Therefore, a storage tank with a capacity higher than the given pressure is needed.The stored hydrogen is portable and can be used as a fuel for the input hydrogen of fuel cells.In this study, it is assumed that the pressurizing power of the storage tanks is derived from another source.
Energies 2016, 9, 332 9 of 11 The system voltage is another factor which is shown in Figure 8.The proposed system operates at voltages around 35 V.The peak voltage point for the system is 37 V which can produce the maximum amount of hydrogen.As we have not used power electronics devices in our system to control the voltage and make the system work at the maximum power point of PV output, the voltage is set by the connecting point of the U-I characteristics of PV system and electrolyzer.At the end of the operation of the PV/hydrogen system, the produced hydrogen is stored in storage tanks.The storage pressure of a tank for daily operation is shown in Figure 9. Maximum pressure for the storage tank with the capacity of 0.25 m 3 is 3.62 MPa.Therefore, a storage tank with a capacity higher than the given pressure is needed.The stored hydrogen is portable and can be used as a fuel for the input hydrogen of fuel cells.In this study, it is assumed that the pressurizing power of the storage tanks is derived from another source.

Conclusions
The operation of a combined PV-electrolyzer system to produce hydrogen as a fuel is designed and simulated.The solar radiation for one day is used for solar calculations.The maximum hydrogen production and minimum average excess power are considered as the optimization indices, which are 697 mol and 0.029 kW, respectively.This condition occurs when the tilts are in their appropriate

Conclusions
The operation of a combined PV-electrolyzer system to produce hydrogen as a fuel is designed and simulated.The solar radiation for one day is used for solar calculations.The maximum hydrogen production and minimum average excess power are considered as the optimization indices, which are 697 mol and 0.029 kW, respectively.This condition occurs when the tilts are in their appropriate location with an angle of 38 ˝.The ratio of the PV array's average power production to its maximum power is about 0.23, which shows how efficient it works.The results can be extrapolated for various other power values for the same rate, although the results have been shown herein only for a 10 kW electrolyzer.For further studies, batteries can be added to system to decrease the size of the electrolyzer.The excess surplus power which is wasted because of the excess over the nominal power of the electrolyser can be stored in batteries to be delivered to the electrolyser at other times.

Figure 2 .
Figure 2. Cell voltage of an electrolyzer in the 20-80 °C temperature range

Figure 2 .
Figure 2. Cell voltage of an electrolyzer in the 20-80 ˝C temperature range.

Figure 4 .
Figure 4. Annual Energy received to tilted surface.

Figure 4 .
Figure 4. Annual Energy received to tilted surface.

Figure 6 .
Figure 6.Power loss of the system.

Figure 5 .
Figure 5. Energy received on a tilted surface.

Figure 6 .
Figure 6.Power loss of the system.

Figure 6 .
Figure 6.Power loss of the system.

Figure 8 .
Figure 8. Operating voltage of the system.

Figure 9 .
Figure 9. Pressure of the hydrogen storage tank.

Figure 8 .
Figure 8. Operating voltage of the system.

Figure 8 .
Figure 8. Operating voltage of the system.

Figure 9 .
Figure 9. Pressure of the hydrogen storage tank.

Figure 9 .
Figure 9. Pressure of the hydrogen storage tank.

Table 1 .
Characteristics of the photovoltaic (PV) modules.

Table 1 .
Characteristics of the photovoltaic (PV) modules.

Table 5 .
Pseudocode for ICA optimization steps.