Energy Management Strategy for a Bioethanol Isolated Hybrid System : Simulations and Experiments

Renewable energy sources have significant advantages both from the environmental and the economic point of view. Additionally, renewable energy sources can contribute significantly to the development of isolated areas that currently have no connection to the electricity supply network. In order to make efficient use of these energy sources, it is necessary to develop appropriate energy management strategies. This work presents an energy management strategy for an isolated hybrid renewable energy system with hydrogen production from bioethanol reforming. The system is based on wind-solar energy, batteries and a bioethanol reformer, which produces hydrogen to feed a fuel cell system. Bioethanol can contribute to the development of isolated areas with surplus agricultural production, which can be used to produce bioethanol. The energy management strategy takes the form of a state machine and tries to maximize autonomy time while minimizing recharging time. The proposed rule-based strategy has been validated both by simulation and experimentally in a scale laboratory station. Both tests have shown the viability of the proposed strategy complying with the specifications imposed and a good agreement between experimental and simulation results.


Introduction
Nowadays, energy production from renewable sources, like solar or wind energy and biofuels, is receiving worldwide attention.Fossil fuels have low efficiency with a high level of emissions that contribute to global warming.In addition, the exploitation of hydrocarbons involves serious social and environmental consequences.Renewable energy sources (RESs) are technologies that are contributing significantly to the diversification of the energy matrix due to their sustainability characteristics and low emissions.However, a common drawback is the unpredictable nature and climate dependence of the RESs, which is more serious in isolated systems that do not have a connection to the electricity grid.However, this problem can be overcome by combining different energy sources and using energy storage units.Such systems are called hybrid renewable energy systems (HRESs) [1,2].
The combination of photovoltaic (PV) and wind energy generation has been reported in several works for isolated applications [3,4].However, hybrid systems need energy storage elements to ensure the power balance in the face of sudden changes in power demand/generation.For short-term energy storage, batteries are the most common storage devices due to the high round-trip efficiency and relatively fast time response, which allow mitigating the harmful effect of RESs fluctuations [5].However, during long periods with low generation, an isolated HRES requires a complementary energy source to guarantee the reliability of the system.In our system, this role is fulfilled by polymer electrolyte membrane fuel cells (PEMFCs).PEMFCs generate electricity using hydrogen as the raw material with a higher efficiency than other technologies (e.g., internal combustion engines).Hybridization between fuel cells and wind-solar sources has been reported in several works [6,7].
Regarding the production of hydrogen needed to operate the PEMFC, electrolysis is a widely-used method [8].However, an interesting alternative is bioethanol reforming.The advantages of bioethanol, compared with other fuels, are high hydrogen content, availability, low toxicity and safe storage and handling.Bioethanol reforming for hydrogen production to feed a PEMFC is reported in several works [9][10][11].In addition, bioethanol can be produced by fermentation of biomass sources such as agricultural residues, forestry, livestock and urban organic material [12,13].Furthermore, bioethanol production in agricultural countries like Argentina is significant and is continuously growing as a result of environmental policies that encourage the change of the energy matrix, making it less dependent on fossil fuels [14].Therefore, bioethanol reforming is a promising alternative to produce the hydrogen required by the PEMFC.Henceforth, the reformer that produces hydrogen from bioethanol, with additional cleaning stages, will be called the bioethanol reformer subsystem (BRS).
The energy management strategy (EMS) is the algorithm in charge of the operation of the HRES subsystems.The EMS establishes in real time the power references of each subsystem to comply with the load balance and optimize the performance of HRES.Different approaches of EMS in hybrid systems with multiple renewable energy sources are found in the literature.A complete review can be found in [15].Rule-based strategies are widely used due to the simplicity of their implementation and their low computational cost.The EMS can be represented by a set of possible states in which the power references for each element are defined.The system switches between the states according to heuristic rules [16].More complex strategies based on model predictive control (MPC) and optimization strategies have been widely used in power management applications in recent years [17,18].However, this kind of methodology generally involves high implementation difficulties.Three EMSs based on switching between different modes are presented and compared in [19].The authors applied the proposed EMS to an isolated power system based on a PV array and wind generators that use water electrolysis to produce hydrogen for future use in a PEMFC.In [20], a rule-based EMS implemented through a state machine approach was proposed for a wind/PV/BRS/Fuel Cell system.This strategy is compared in [21] with another that uses the concept of a receding horizon with predictions of the future generation from the renewable sources, the future load and the state of charge in the battery bank.In [17], the authors studied the role of EMSs in microgrids with water electrolysis and presented experimental results.In [22], the authors proposed an energy management strategy based on multiple operating states for a DC microgrid, which is comprised of a PV array, a PEMFC system and a battery bank.This strategy can share the power properly and keep the bus voltage at the desired value under different operating states.However, the design of EMSs for isolated HRESs with bioethanol-based hydrogen production from theoretical and experimental approaches is not addressed in the literature.
This work presents an EMS for an isolated HRES composed of wind-PV generation, a battery bank and a PEMFC, where the hydrogen required by the PEMFC is produced from bioethanol reforming.The proposed EMS is based on a finite state machine.First, the EMS is evaluated by simulation using historical data over long periods of time to study the reliability of the strategy during different scenarios.Then, an experimental validation of the EMS is carried out using the profiles of a typical day.For this purpose, a methodology is proposed to scale the initial scenario so that the tests can be performed in the laboratory.This methodology preserves the most relevant energy-flow characteristics.Finally, the experimental results obtained are compared with those corresponding to the simulations.
The organization of the rest of this paper can be summarized as follows.In Section 2, a description of the HRES and the corresponding mathematical models are addressed.Section 3 describes the proposed energy management strategy.In Section 4.1, a possible real scenario is described.In Section 4.2, the simulations results of the proposed EMS are shown.In Section 4.3, the laboratory test station is described, and in Section 4.4, the scaling methodology is presented.The experimental results are shown and discussed in Section 4.5.Finally, the conclusions and future work lines for further research are stated in Section 5.

System Description
The topology of the HRES under study in this work is shown in Figure 1.The wind and PV sources, as well as the PEMFC are connected to the DC bus via unidirectional DC-DC boost converters.In the case of the battery, a bidirectional boost DC-DC converter must be used to ensure the charging and discharging mode.The power references are imposed by the EMS (explained in detail in Section 3).The size of the system was determined using the sizing methodology presented in [23], where the developed methodology allows one to obtain a sizing for different scenarios and to fulfil the performance specifications.The characteristics of the components of the system in terms of their rated powers or capacities are summarized in Table 1.The power converters in the model are assumed ideal and perfectly controlled.

Bioethanol Reformer Subsystem
Steam reforming is one of the most widely-used methods to obtain hydrogen from bioethanol because it achieves a high concentration of hydrogen, necessary to feed the PEM fuel cell subsystem [24].The main reactor is the ethanol steam reformer (ESR), where most of the conversion of bioethanol into hydrogen takes place.Then, a clean stage composed of three water-gas shift reactors (WGS) removes the carbon monoxide that can damage the catalyst of the fuel cell [25].The endothermic condition of the ESR reaction implies the supply of heat.A solution reported in [26] consists of burning a portion of the bioethanol in a separate reactor to provide the reaction heat from the gases obtained.However, this solution has a negative impact on the efficiency due to the extra bioethanol consumption.To avoid this situation, the burner is replaced by an electric heating subsystem, which is powered from the DC bus [27].
The amount of heat needed in the reaction is ∆H 0 298 = 254.8kJ mol −1 [23].The power required by the heater (P h ) can be obtained from the stoichiometric equation of the reaction: where W H 2 is the rate of hydrogen consumed by the PEMFC.

PEM Fuel Cell Subsystem
In this work, the voltage-current and power-current curves in the data sheet of a 1.2 kW Nexa module are used to build the model [28].The rate of hydrogen consumed in the electrochemical reaction in the PEMFC anode can be approximated by the following expression depending on the current (I f c ): where W H 2 is the rate of hydrogen reacted, η F, f c is the Faraday efficiency (90-95%), n is the number of cells in the FC (a typical value of 26 cells is considered) and F = 96485 C mol −1 is the Faraday number.This hydrogen flow is supplied by the BRS presented in Section 2.1.
The power of the FC must be maintained within certain limits.The upper limit is the nominal power of FC, while the minimum is fixed to guarantee a minimum efficiency.Therefore, the following constraints can be formulated: with a P min f c = 300 W and a P max f c = 1200 W. In addition, the PEMFC power cannot be increased faster than a certain rise rate power (∆P max f c ) to avoid a lack of reactants and the power cannot be decreased faster than a certain fall rate power (∆P min f c ) to prevent overpressure in the PEMFC [29].Therefore, the following constraints are added: where and ∆P min f c = −0.1 P max f c [30].

Wind Power Subsystem
The maximum power that a wind turbine is able to generate can be expressed by: where r is the radius of the blades, ρ a is the air density and v is the wind velocity.The power coefficient (C p ) expresses the efficiency of the turbine with respect to the wind power available.The coefficient C P is a function of λ and the blade pitch angle β, where λ = rω m /v and ω m is the speed angle of the turbine shaft.Then, C p can be expressed by [31]: where C f is a design constant of the blades.
When the wind speed is below the cut-in speed (v w,cin ), no power is generated.Above this value, the generated power responds to Equation (5).When the rated power of the machine (P w,r ) is reached at v w,r , the output power must be operated at constant output power in order to prevent damage.The turbine must be stopped if the wind speed is higher than its cut-out limit (v w,cout ) [32].In this work, a perfect control of the power delivered by the turbine is assumed, working in maximum power point tracking (MPPT) or load power tracking (LPT).The values of the parameters used in the model are presented in Table 2.The value for the power coefficient is C p (λ, β) = 0.59, resulting in a 4.5 kW p wind turbine.

Solar Power Subsystem
The solar power subsystem is composed of a PV array connected to the DC bus through a DC/DC power converter.A PV array is a group of several PV panels, which are electrically connected in series and parallel circuits to generate the required voltage and power.In turn, each PV panel is composed of several PV cells.
PV cells have voltage-current and power-current nonlinear characteristics, strongly dependent on radiation and temperature.A possible expression for the output current I pv of a PV cell is given by [33]: where I ph is the generated current under a reference radiation, I rs is the cell reverse saturation current, V pv is the voltage level on the PV cell terminals, q is the charge of an electron, R s is the intrinsic cell resistance, A c is the cell deviation from the ideal p-n junction characteristic, K is the Boltzmann constant and T is the cell temperature.The photocurrent depends on radiation and temperature according to the following expression: where I sc is the short-circuit cell current at the reference temperature and radiation, K l is the short-circuit current temperature coefficient and ν is the measured radiation expressed in Wm −2 .
The reverse saturation current depends on temperature according to the following expression: where I or is the reverse saturation current at the reference temperature T re f and E go is the band-gap energy of the semiconductor used in the cell.Thus, the available current for a PV panel can be expressed as follows: where V PV is the voltage level in the PV panel terminals, n p is the number of parallel strings and n s is the number of serial connected cells per string.Therefore, the available power generation from a PV panel is: The values of the constants used in the model are given in Table 3 and correspond to a commercial panel of 250 W p .In this work, four PV panels are connected in parallel in a single PV module, obtaining a peak power for each module of 1 kW p .Then, the sizing strategy calculates how many modules to install in the system.The maximum PV available power P pv at each instant is obtained through an MPPT algorithm.

Parameter
Value 60 number of PV panels per module 4

Energy Storage Subsystem
The natural variability of wind-solar power jointly with the isolated application implies the need to include energy storage devices having a reliable power supply.Batteries are the most widely-used energy storage devices in isolated applications due to their high specific energy.
In this work, the energy storage subsystem (ESS) is composed of batteries.The model presented in [34], integrated in the SimPowerSystems toolbox of MATLAB, is used.This model is based on an equivalent circuit with a variable voltage source and a serial internal resistance.The variable voltage source and the internal resistance depend on the battery state of charge (SoC), which represents the remaining capacity still available to be discharged from the battery [35].The SoC can be defined as follows: where C bat is the capacity of the bank of batteries.The SoC must be kept within certain limits to ensure good performance and maximize battery lifetime.The main characteristics of the batteries used are in Table 4.

Energy Management Strategy
The energy management strategy is the algorithm that controls the power references of the local controllers to guarantee the power balance in the DC bus while complying with the constraints of the system.The power balance in the DC bus can be expressed as follows: where P w , P pv , P f c and P bat are the power provided by the wind turbines, the PV array, the fuel cell and the batteries, respectively, while P h and P load are the power consumed by the electric heater and the load, respectively.The power delivered by the batteries is considered positive in the discharge mode and negative when charging.Two important specifications are introduced for the design of the EMS.The first, and the most important in this type of isolated systems, is the autonomy time.It refers to the capacity of the batteries to support the full load demand until it falls below the minimum SoC level.When the load demand exceeds the generation provided by the RESs and the fuel cell, the batteries allow meeting the power balance until the SoC falls below SoC min .When this situation occurs, only the vital loads (P vital load ) are fed.The vital loads are the essential loads that assure the minimum performance of the system.The second specification introduced is the recharge time, which is defined as the time the batteries require to recover the capacity to supply all the regular power demand.
According to these specifications, some basic principles are defined to design the EMS: • The RESs deliver their maximum power (MPPT mode), while the batteries do not exceed the maximum charge level.• When the battery reaches its maximum SoC level, the RESs change to LPT mode, where the power reference follows the load demand.• Minimize the battery cycles of charging and discharging to preserve their lifetime.
• The fuel cell provides energy to the DC bus during low RES generation periods and maintains the SoC of the batteries at a desirable level.
The EMS is implemented using a finite state machine approach.This tool is widely used in decision making algorithms and EMSs [36].A state machine makes a transition between modes (states) according to specific events.In each state, the power references for the local controllers are defined.The algorithm implemented can be represented by a statechart, which is a graphical way to describe a finite state machine, as shown in Figure 2. The description of the events is given in Table 5, while the mathematical description of each state is presented below: • State 1: • State 2: • State 3: P re f bat • State 4: • State 5: where P re f res , P re f f c and P re f bat are the power references to the RESs, the fuel cell and the battery, respectively.Since P h is the power consumed by the heater, P req load is the load power demand and P del load is the power actually delivered to the load.To simplify the notation, the two available powers from the RESs are considered together as P av res = P av pv + P av wind .) is defined between SoC min and SoC max .With an SoC value between SoC L1 and SoC max and the power generated by the RESs higher than the load demand, the EMS is in state S 1 .In this state, the RESs work in MPPT mode; the full load demand is supplied; and the battery is charged with the surplus power.When the SoC reaches the maximum value, defined by Event 4, the EMS changes to state S 4 , and the RESs start to work in LPT mode.On the other hand, if the generation from RESs is lower than the load demand, the battery is discharged to meet the power balance.When the SoC falls below SoC L1 − SoC th , the EMS changes to state S 2 , and the fuel cell is turned on.SoC th is a threshold, which is used to avoid repetitive changes in the states.
In the state S 2 , the fuel cell is in charge of generating the remaining power between P av ren and P req load , plus the power required by the heater and a constant power used to charge the batteries (P ch bat ).The system remains at this state until SoC ≥ SoC L1 .Then, the EMS changes to state S 3 , and P ch bat falls to zero.The fuel cell is turned off when P av ren > P req load .This condition avoids continuous cycles of charging and discharging of the battery.
If the RESs power plus the PEMFC power is not enough to meet the power demand, the batteries are discharged until falling below the SoC min .If this condition occurs (Event 3), the EMS changes to the emergency state S 5 , where only the vital loads are fed.When the SoC reaches again SoC L1 , all the loads are fed again.

SoC Levels Calculation
To simplify future calculations, it is assumed that the battery voltage remains constant at its nominal voltage (V r bat ).Then, the SoC can be expressed in discrete form by Equation (12): for k = 0, 1, . .., with sampling time ∆T.At the beginning of a discharge cycle, the SoC is greater than SoC L1 , where SoC L1 is defined as: Then, from the definition of autonomy time and replacing in Equation ( 34), ∆SoC results: where P max f c is the maximum PEMFC power, P max h is the power required by the heater when the PEMFC works at P max f c and P max load is the maximum load power allowed by the system.Note that both P max f c and P max h are determined by the PEMFC sizing, while P max load is an inherent characteristic of the system.The parameter T auto is a minimum value for this specification.Usually, the SoC min level is imposed to avoid damaging the batteries and is data provided by the manufacturer, while the maximum values P max load , P max f c and P max h are adopted according to the system characteristics.When the EMS reaches State 5, only the vital loads are fed, while the batteries are charged with a constant power equal to P ch bat .As shown in Figure 2, the EMS leaves this state when Event 1 occurs.From Equation (36), and taking into account the definition of recharge time (T ch ) as a system specification, the value of P ch bat is:

Experimental Validation
The performance evaluation of an EMS oriented toward stationary applications, like houses or small buildings, requires the evaluation of the system for at least one full day.In this work, simulations are carried out in order to verify controller reliability through different RES generation and load demand scenarios.Once the EMS was tested through simulations, an experimental test was performed to validate the simulations results.

Scenario Description
Representative profiles are used over a week to show the performance of the EMS under different climatic conditions and, therefore, variable scenarios of generation from RESs.Radiation, temperature and wind speed profiles were collected in the Municipal Observatory of Rosario.They correspond to the first week of January 2011 and are shown in Figure 3a,b.Then, the power profiles generated from the RESs (PV + wind) are obtained through the models presented in Section 2 with the sizing values of Table 1.In addition, the demand profile corresponds to a standard residential building [37].The resulting power profiles are shown in Figure 4.These power values exceed the laboratory station rated powers used to validate the proposed EMS (laboratory station characteristics are presented in Section 4.3).Therefore, a scaling process in terms of power, energy and time was required to properly reproduce a realistic scenario.Before testing the EMS in the laboratory test station, simulations were carried out to verify the behavior of the EMS.

Simulation Results
A complete model of the system was implemented in MATLAB/Simulink according to the description of the models of Section 2. The system specifications, SoC levels, the constant battery charging power, as well as maximum power adopted are shown in Table 6.The profiles corresponding to the generation from RESs and the load demand shown in Figure 4 are used to carry out the simulation of the EMS performance.Figure 5a shows the resulting power distribution on the HRES.As can be seen, the specifications stated in Section 3 are satisfied by applying the proposed EMS.The RESs work in MPPT most of the time, and the batteries are charged if the generation from RESs exceeds the load demand, until SoC max is reached.At this point, the RES goes to work in LPT mode (e.g., between Hour 100 and 116).The PEMFC acts as a secondary source when the RESs produce less than the load demand, while the batteries ensure the power balance during transient periods.Figure 5b shows the behavior of the SoC and the states in which the EMS works.It can be seen that the EMS ensures the power balance, trying to keep the SoC at the desired value (SoC L1 ).This value is between SoC min and SoC max , giving the battery the ability to absorb or deliver energy if there is an excess or lack of power in the DC bus, respectively.Finally, it can be mentioned that the simulation results obtained have demonstrated the ability of the EMS to manage the energy in the system.
It is necessary to evaluate the performance of the EMS and verify the compliance with the system specifications (autonomy time and recharge time) in a critical scenario.Figure 6a shows the power distribution in the simulated scenario, while Figure 6b shows the corresponding SoC evolution and states of the EMS.At time t = 20 h, the generation from the renewable sources is set to zero, and the load demand is set to its rated value.As a result of this, the fuel cell and the battery feed the load demand, and the battery is discharged during the autonomy time from SoC L1 to SoC min .At this moment, the EMS changes from State 2 to State 5, where only the vital loads are fed.Then, batteries start to charge (recharge time) until the SoC L1 is reached, and the EMS leaves State 5.At this point (t = 30 h), the critical scenario is reset, and the normal conditions are restored.The simulations show an adequate fulfillment of the imposed specifications.

Laboratory Station Description
The experimental validation was performed in the Fuel Cell Laboratory of the Institut de Robòtica i Informàtica Industrial (IRI) of the Technical University of Catalonia (CSIC-UPC), in Barcelona, Spain.In this work, a fuel cell test station is used, which consists of a PEMFC stack and a supercapacitor (SC) bank.The SC is used in this experimental setup to replace the batteries as the energy storage.The demand profile is emulated by a programmable load, while the RESs generation is emulated with a programmable current source/sink.A picture of the laboratory station is shown in Figure 7.This station has already been used successfully to validate an EMS applied to hybrid vehicle in [30].
The PEMFC stack is a 1.  7.

Scaling Methodology
In this work, a method is developed to scale a given scenario so that it can be reproduced in the test station.The scaling methodology preserves most of the relevant energy-flow characteristics.Henceforward, the term scaled (superscript s) will be used to refer to the test station, while the term real (superscript r) to the HRES described in Section 2 together with the realistic scenario.
Figure 8 shows the experimental setup used.The EMS is executed in a host computer running a LABView TM platform.The PEMFC stack and the SCs are connected to the DC bus through DC-DC converters.The H&H Source/Sink is used to emulate the RES generation delivered to the bus (P s ps ).The scaled version of RES's available power (P s res ) is obtained using the off-line scaling methodology applied to the real RES generation profile (P r res ), which is formed by the wind turbines' (P r wind ) and PV array's (P r pv ) available powers.The scaled version of the load demand (P s load ) and the power consumed by the BRS heater (P s h ) are emulated using the programmable load (P s pl ).Note that the energy required by the heater is obtained using the power delivered by the Nexa TM .Given a real system with a net power P r during a period T r (Figure 9a), the stored energy results ∆E r = P r T r (Figure 9b).In the same way, in the laboratory test station, there is an equivalent scenario with a power P s and time T s (Figure 9c), resulting in a stored energy ∆E s (Figure 9d).The relative change with respect to the maximum energy capacity E i max (with i = r, s) can be expressed as: The relationship between the real scenario and the scaled scenario can be defined taking into account that the power step must produce the same relative change in the stored energy of both scenarios.From Equation (38), it results: The scaling coefficients can be defined by Equation (39): while the fundamental relationship can be rewritten as: where k e , k p and k t are the energy, power and time scaling coefficients, respectively.Then, the load and RES scaled profiles result: P s res = k p P r wind + P r pv , P s load = k p P r load + P s h . (44) Note that the power consumed by the heater in the scaled system (P s h ) is not multiplied by k p because it is calculated by Equation (1) using the measure of the power delivered by the fuel cell (P s f c ).As described in Section 2.5, the ESS in the real system is formed by a battery bank.Then, the maximum energy storing capacity (E r max ) can be approximated by: Since the ESS in the test station is composed of SCs, the maximum energy storing capacity (E s max ) is calculated as: where C sc is the capacity of the SC and V sc max and V sc min are the maximum and minimum voltage imposed on the SC in the emulation test, respectively.Note that the definition of a minimum and a maximum value in the SC voltage allows us to manipulate its energy capacity, which is one of the scaling process variables.The SoC sc level from the SC can be calculated as: where V sc (k) is the SC voltage.Then, replacing (40), ( 45) and ( 46) in (43), the power scaling coefficient results: To avoid exceeding the laboratory station capacities established in Section 4.3 and taking into account the case study description, some constraints must be imposed on Equation (48).The power scaling coefficient (k p ) must take into account the power range of the real system profiles (RESs, load demand and fuel cell) in such a way that the scaled powers do not exceed power equipment limitations.These constraints impose upper bounds on k p : where P ps max is the rated power of the source/sink and P pl max is the rated power of the programmable load.P r res,max is the maximum power peak of the RES generation profile; P r load,max is maximum load demand; and P r h,max is the power consumed by the electric heater with the fuel cell at its maximum power.P r f c,max is the rated power of the fuel cell in the real scenario, while P r f c,max is the rated power of the Nexa TM fuel cell.
In addition, the maximum current of the SC (I sc max ) also imposes an upper bound on k p .In charging mode, the maximum current occurs with the maximum charging power and the minimum voltage.
With the proposed EMS, the maximum charging power matches the maximum RES generated power (P r res,max ), while the minimum voltage is V sc SoC,min : In discharging mode, the maximum power is when the SC has to provide all the power required by the load: It is expected that the maximum power rate allowed in the laboratory station fuel cell (∆P s f c,max ), once unscaled, would not exceed the real system constraints (∆P r f c,max ): This constraint imposes a lower bound on k p that depends on k t : In Figure 10 are summarized the feasible values of the scaling coefficients (k p , k t ) and the constraints described before.As can be seen, k p can take values in the dashed line (defined by Equation ( 48)) between the lower bound (k p lb ) and the most restrictive upper bound (k p ub,min ).To obtain an experimental test duration as short as possible, a value of k t as small as possible has to be chosen.Then, the adopted scaling coefficient pair ( kp , kt ) is in the intersection defined by Equation (48) and k p ub,min .Finally, Table 8 summarizes the characteristic values of the real scenario, the adopted voltage values for the SCs in the emulation test and the resulting scaling coefficients.The values of the maximum RES generation and the maximum load demand were obtained from the real scenario taking into account a complete year.

Experimental Results
The experimental validation of the EMS was performed using a typical winter day from the annual profiles of the real scenario.Figure 11 shows the evolution of power generation from renewable sources (blue line) and the load power demanded (red line) over the selected day.With regard to the power generated, this power is the sum of the powers corresponding to the two renewable sources: wind and photovoltaic.This power has a peak at 15 o'clockdue to the contribution of sun and wind and a valley during the night because there is no contribution from the Sun.With regard to the curve of the power demanded, it is a profile consistent with the power demand of a home where two peaks are observed, one in the morning near noon and another near dinner time, and two valleys, one in the morning and the other in the afternoon.Then, the power profiles were scaled according to the methodology presented in Section 4.4 to perform the experimental test.Note that, once scaled, the original length of 24 × 3600 = 86,400 s resulted in a duration of the experiment of 86,400 ×k t = 682.5 s.The reduction in experiment duration is an important advantage that results from applying the proposed scaling methodology.

time (hours)
Figure 12a shows the power distribution resulting from the test.As explained in Section 4.3, the RES generation is emulated with the source/sink, while the load demand and the power consumed by the heater are emulated with the programmable load.Note that, while the load demand is scaled from the load profile of Figure 11, the power required by the heater is obtained from Equation (1) using the measured power of the Nexa TM fuel cell.The power delivered by the fuel cell is set by the EMS, and the power delivered/absorbed by the SC arises as a result of the two PI controllers in cascade that attempt to maintain the DC bus voltage at the desired level, while in the simulation, it is computed from a power balance.A constant power delivered by the SC can be seen that corresponds to the necessary power to maintain the DC bus at the desired value.Then, Figure 12b shows the evolution of the state of charge of the batteries and the state in which the EMS works.The commutations between states according to the occurrence of events can be seen.The HRES behavior is simulated using the profiles of Figure 11 and compared with the experiment results.The inverse scaling process is carried out over the profiles presented in Figure 12a.The evolutions of the main variables, both in simulation and in the experimental test, are shown in Figure 13.The RES generation and the load power demand, which are emulated with the source/sink, and the programmable load are shown in Figure 13a,b.In Figure 13c is shown the Nexa TM fuel cell, where the main difference of the experiment evolution with respect to the simulation is a small noise.The major differences observed between experimental and simulation results are found in the SC power of Figure 13d and are due to the low efficiency in the bidirectional DC-DC converter, especially at low powers.This could be solved by incorporating the models of the converters with their efficiency curves.This improvement will be addressed in future works.In summary, the experimental results presented in this section have demonstrated feasibility of implementing the proposed EMS for the HRES under study.In the same way, the proposed scaling methodology provides a useful tool to evaluate the qualitative behavior of the main HRES variables in a real scenario.

Conclusions
This paper describes an energy management system for an isolated hybrid renewable energy system.In this sense, a comprehensive and innovative approach was presented and validated.The EMS proposed is based on a finite state machine designed to obtain good performance in the autonomy and recharge time.The performance of the EMS has been evaluated using historic profiles representing different weather and load conditions; firstly, running the simulation using a MATLAB/Simulink complete model and later in a scaled experimental laboratory setup.In both cases, the EMS offered correct results in terms of compliance with the imposed specifications, and also, the experimental results obtained agree with those obtained by the simulation.As a conclusion, it can be established that the proposed EMS can be implemented easily in real time and offers adequate results.Finally, it is worth mentioning that although the results are focused on a specific HRES, the methodologies presented in this document are general enough to cover a wide range of problems related to energy management in HRES.

Figure 1 .
Figure 1.Schematic diagram of the isolated hybrid renewable energy system with hydrogen production from bioethanol.
and ∆T is the sample time.The values used are ∆P max f c = 0.1 P max f c

Figure 2 .
Figure 2. Statechart diagram describing the energy management strategy based on a finite state machine.

Figure 3 .
Figure 3. Weather variables' profile of the first week of January 2011: (a) (dashed line) temperature and (solid line) radiation; (b) wind velocity.

Figure 4 .
Figure 4. PV and wind turbine energy generation profile and load power demand.

Figure 5 .Figure 6 .
Figure 5. Simulation results: (a) HRES power distribution; (b) (dashed line) state of charge of the batteries (SoC); (solid line) state transitions in the EMS.
2 kW Nexa TM model MAN5100078, manufactured by Ballard TM .The output voltage can be operated safely in the linear range of voltages from 26 V to 36 V and currents from 10 A to 46 A. The ESS consists on a 165 F Maxwell TM Supercapacitor model BMOD0165.The nominal voltage is 48 V with values up to a maximum of 52 V.The SC supports a nominal current of 98 in continuous mode and shows a serial resistance of approximately 6 mΩ.The DC-DC boost converters are implemented with two branches of Semikron TM IGBTs.The converters are controlled with a 20 kHz PWM signal, whose switching duty cycle is set by the EMS.The maximum voltage allowed in each switch is 400 V, and the nominal current is 75 A. The inductances in the converters have a value of 35 µHy.The DC bus voltage adopted is 75 V, and the bus capacitance is 2720 µF.The programmable source is an NL Source-Sink (SS) from Hocherl & Hackl GmbH TM .The output voltage of this source can reach a maximum of 80 V with a nominal power of 3.2 kW.The installed programmable load (PL) is a ZL Electronic DC Load from Hocherl & Hackl GmbH TM , which can work with a maximum voltage of 80 V and 3.4 kW.Finally, the control is implemented using a National Instruments TM controller model Compaq Rio 9035, which has a CPU Dual-Core of 1.33 GHz and an FPGA Xilinx Kintex-7 7K70T.The rated values of the fuel cell test station are summarized in Table

Figure 7 .
Figure 7. Fuel cell test station in the laboratory of IRI-CSIC-UPC.

Figure 8 .
Figure 8. Laboratory system control scheme and experiment configuration.The dotted boxes show the elements of the real system to be emulated, and the arrow indicates the laboratory equipment used to emulate this element.

Figure 9 .
Figure 9. (a) Real system, net power P r ; (b) real system, energy variation ∆E r ; (c) scaled system, net power P s ; (d) scaled system, energy variation ∆E s .

Figure 10 .
Figure 10.Solid line: feasible values of the scaling coefficients (k p , k t ) resulting from Equation (48); dashed lines: lower bound (k p lb ) and most restrictive upper bound (k p ub,min ).o points to the adopted pair ( kp , kt ).

Figure 11 .
Figure 11.Power profiles for testing: RES generation (wind + PV) and load power demand.

Figure 12 .
Figure 12.Laboratory station power profiles: (a) source/sink, programmable load, Nexa TM fuel cell, Maxwell TM SC; (b) (dashed line) state of charge of SCs (SoC sc ) and (solid line) states of the EMS.

Table 3 .
Parameters of each PV panel (250 W p ).

Table 4 .
Parameters of each unit of the energy storage subsystem (1352 Wh).

Table 5 .
Description of the events in the set Σ.
ren An intermediate SoC level (SoC L1

Table 6 .
Real system specifications and EMS parameters.

Table 7 .
Fuel cell test station rated values.

Table 8 .
Experimental validation data and parameters of the scaling methodology.