# Introduction he penetration level of green and renewable energy sources/distributed generation units are expected to grow in the near future as there is a probability of rundown conventional fuels for power generation. The distributed generation is classified as renewable and non-renewable. The distributed generation sources such as Fuel cells, Wind and Solar energy are increasing daily due to increase in demand for electrical power [1]. These energy sources are environmental friendly, reduces transmission and distribution losses, peak load shaving, can be used as backup sources and etc. Fuel cell is a promising device as it is efficient, modular and can be placed at any site for improving system efficiency [2] but it has slow startup response. Solar energy is an important renewable energy source [3] but the intermittent nature of this technology is a major issue. The accessibility of energy is driven by climate and cell temperature however not on the loads of the systems. This innovation can be marked as irregular and typically PV array utilizes a Maximum Power Point Tracking (MPPT) (MPPT) strategy to consistently convey the most highest power to the load when there are variations in irradiation and temperature [4].Because of the intermittent nature of PV array it becomes an uncontrollable source In order to overcome the drawbacks with the slow start-up of fuel cells and intermittent nature of PV cell a nerofuzzy controlled grid associated hybrid photovoltaic as well as proton exchange membrane fuel cell (PEMFC) distributed generation system with battery as energy storage is suggested in this paper. # II. # System Modeling Fig. 1 shows the block diagram of the HRPS (Hybrid Renewable Power Sources) proposed in this paper that connected to main grid in Point Common Coupling (PCC). So by above discussions tww mathematical models recitation the dynamic behavior and each of these constituents are given below. ?? ð??"ð??"ð??"ð??" = ?? 0 (?? 0 + ???? 2?? (log ? ?? ?? 2 ?? ?? 2 0.5 ?? ?? 2 ?? )? ? ?? ?????? ??(1) Where: ?? 0 is the number of cells connected in series; ?? ?????? is the internal resistance of fuel cell stack [?] ?? 0 is the voltage related with the reaction free energy; R is the universal gas constant; T is the temperature; I is the current of the fuel cell stack; F is the Faraday's constant. ?? ?? 2 , ?? ?? 2 ?? , ?? ?? 2 are dictated by the accompanying differential equations. ?? ?? 2 = ? 1 ?? ?? 2 (?? ?? 2 + 1 ?? ?? 2 ??? ?? 2 ???? ? 2?? ?? ?? ð??"ð??"ð??"ð??" ?) ?? ?? 2 ?? = ? 1 ?? ?? 2 ?? (?? ?? 2 ?? + 2 ?? ?? 2 ?? ?? ?? ?? ð??"ð??"ð??"ð??" ))(2)?? ?? 2 = ? 1 ?? ?? 2 (?? ?? 2 + 1 ?? ?? 2 ??? ?? 2 ???? ? ?? ?? ?? ð??"ð??"ð??"ð??" ?) Where, ?? ?? 2 # ???? and ?? ?? 2 ???? are the molar flow of hydrogen and oxygen and where the Kr constant is welldefined by the relation between the rate of reactant hydrogen and the fuel cell current: ?? ?? 2 ?? = ?? 0 ?? 2?? = 2?? ?? ??(3) Besides, a straightforward model of reformer that creates hydrogen through methane has been deliberated and the model is second-order transfer function.So therefore itsmathematical form can be written as follows: ?? ?? 2 ?? ?????????????? = ???? ?? 1 ?? 2 ?? 2 + (?? 1 + ?? 2 )?? + 1 Where?? ?????????????? is methane flow rate [kmol/sec]; CV is conversion factor [kmol of hydrogen per kmol of methane]; ?? 1 , ?? 2 are reformer time constants [sec]. b) Photovoltaic Model (PV) The equivalent circuit shown in Fig. 2 is a one diode model of a solar cell which consists of a diode and a current source connected in parallel with a series resistance R s . The current source foods the photocurrent I ph , which is directly proportional to solar irradiance G. By referring manufacturer's data sheet, the two main parameters used to describe a PV cell are open circuit voltage and another is its short circuit current. The mathematical model [3] of PV cell can be expressed as ?? ???? = ?? ??? ? ?? 0 ?exp ?? ? ?? ???? + ?? ???? . ?? ?? ?????? ? ? 1? (4) Since Photocurrent ?? ??? is directly proportional to solar radiation G. I ph (??) = ?? ??ð??"ð??" ?? ?? ????ð??"ð??"(5) The short-circuit current I SC of solar cell depends linearly on cell temperature. ?? ??ð??"ð??" ??? ?? ? = ?? ??ð??"ð??" [1 + Î?"?? ??ð??"ð??" (?? ?? ? ?? ??????ð??"ð??" ] ?? ??? ???, ?? ?? ? = ?? ??ð??"ð??"?? ?? ?? ð??"ð??"ð??"ð??"ð??"ð??" [1 + Î?"?? ??ð??"ð??" ??? ?? ? ?? ?? ????ð??"ð??" ?(6) I 0 also depends on solar irradiation as well as cell temperature and that can be mathematically expressed as follows ?? 0 ???, ?? ?? ? = ?? ??? ???, ?? ?? ? ?? ? ?? ??ð??"ð??" (?? ?? ) ?? ?? (?? ?? ) ? ? 1 (7) In the writing numerous MPPT strategies are accessible, for instance, incremental conductance (INC), consistent voltage (CV), and perturbation and observation (P&O). Therefore P&O strategy has been generally utilized in view of its basic input structure and less measured parameters. The panel voltage is purposely intentionally agitated (expanded or diminished) then the power is contrasted with the power got before to disturbance. In particular, if the power panel is increased because of the unsettling influence, the accompanying disturbance will be made in the same course and if the power diminishes, the new perturbation is made in the opposite direction.But the demerit with P & O is the outputpower is oscillating in nature. Because of this reason weuse the Fuzzy MPPT technique to deliver the maximumpower and to eliminate perturbations in the output power. # c) Fuzzy MPPT Control The inputs to the fuzzy MPPT control can be measured or computed from the voltage and current of solar panel. Thecontrol rules are indicated in [4]with ?Ppv and ?V pv asinputs and ?V pvref as the output. The membership functions of input and output Variables in which membership functions of input variables ?Ppv and ?V pv are triangular and has seven fuzzy subsets. Seven fuzzy subsets are considered for membership functions of the output variable ?V pvref. These input and output variables are expressed in terms of linguistic variables (such as BN (big negative),MN (Medium negative), SN (small negative), Z (zero), SP(small positive), MP (medium positive), and BP (big positive). # d) Battery Modelling The battery is a device which stores energy inelectro chemical form. Battery is used as energy storage device in wide range of applications like hybrid electric vehicles and hybrid power systems. In this paper, the battery energy storage is combined with hybrid PV/PEMFC distributed generation system. The battery model considered in this paper is shown in fig. 3. The battery model used is based on voltage model proposed by Shepherd [4]. As depicted above boost converter is defined by the following two nonlinear state space averaged equations [7]: ???? ?? ???? = ? ?? ?? ?? ? ? 1??? ?? ? ?? 0 + 1 ?? ?? ?? (8)???? ð??"ð??" ???? = ? 1 ? ?? ?? ? ?? ?? + ?? ?? ?? where ''d'' is the duty cycle of the switching device, ''U'' is the input voltage, ''?? ?? '' is the inductor current, ''?? ð??"ð??" '' is the output voltage and "?? 0 " is the output current. # f) DC-AC Converter Model Converter (VSC) is shown in Fig. 5. To reduce harmonics, LCL filter is connected between the converter and also at the grid side [6]. So therefore the dynamic model of the three-phase VSC is represented in ???? 1?? ???? = ? ?? 1 ?? 1 ?? 1?? + 1 ?? 1 (?? 1?? ? ?? ð??"ð??"?? ) ???? 2?? ???? = ? ?? 2 ?? 2 ?? 2?? + 1 ?? 1 (?? ð??"ð??"?? ? ?? ???? )( 9 )?? ð??"ð??" ???? ð??"ð??"?? ???? = ?? 1?? ? ?? 2?? Whre k={a, b, c} Fig. 5 : DC/ AC Three Phase Inverter # III. Power Control Strategies of Hybrid System By above discussions the Power flow control from hybrid power sources to local AC bus and to/from storage devices is required to keep up power balance at all times while fulfilling the the active and reactive power demanded by the load. Eq.(13) gives power balance expressions that should be satisfied together at the DC-link and at the P CC all the time. The rate and magnitude of fuel cell power P FC and rate, sign and magnitude of battery power P Batt depend on the magnitude and how fast the load changes. ?? ???? = ?? ???? + ?? ???? + ?? ???????? ?? ???????? = ?? ???? + ?? ???????? (10) ?? ???????? = ?? ???? + ?? ???????? According to the control strategy proposed in this paper, ?? ???????? and ?? ???????? are made equal to P ref and Q ref so that the hybrid power system output shadows the load demand only under normal loading conditions also?? ???????? and ?? ???????? are seems to be zero. Soif the local load demand surpasses the hybrid power system capacity, then remaining of the power is supplied from the grid side. Fig. 6 shows the overall structure of the control strategy. ?? ??ð??"ð??" ?? ??ð??"ð??" ???? ??ð??"ð??" ???? = ?? ???? + ?? ???? + ?? ???????? ? ?? ????????( 11 ) Moreover, to meet the requirements of power balance in DC link it is significant to consider the dynamic limitations of fuel cell power. In this case, the fuel cell power could not change rapidly and the fuel cell controller with DCDC converter should regulate the operating point of fuel cell. But the amount of power that should be absorbed by battery energy storage in order to balance the power in DC link is significant and also it is greatly influenced by DC link energy, where its energy measurement is supported with the help of the following calculation: ?? ??ð??"ð??" (??) = (1 2 )?? ??ð??"ð??" ?? ??ð??"ð??" 2 (??) In this paper, a power flow control structure has been established for hybrid power sources during voltage sag. It is based on Fuzzy Logic Control (FLC) strategy that determines the battery energy storage power according to the following inputs: ??(??) = ?? ??ð??"ð??" ?????ð??"ð??" (??) ? ?? ??ð??"ð??" (??) ???(??) = ??(??) ? ??(?? ? 1)(13) where?? ??ð??"ð??" ?????ð??"ð??" is the reference dc link energy which is calculated by reference dc link voltage.Subsequently, it is crucial to outline powerful and stable control technique to ensure the stability of the dc link of hybrid system.For this purpose, a fuzzy neuralcontrol startegy is devloped [8]. IV. # Neuro-Fuzzy Control Startegy In this paper a neuro-fuzzy control strategy, for each of the input, four fuzzy subsets have been employed. These are ZE (zero), L (low), M (medium) and H (high). So for all of these fuzzy sets, a gaussian membership function has been used. As each of the two inputs has four subsets, there are altogether 16 control rules in the neuro-fuzzy logic controller. Fig. 7 : The neuro-fuzzy scheme The neuro-fuzzy calculation utilizes membership functions of gaussian kind. With Gaussian fuzzy sets, the algorithm is fit for using all data contained in the preparation set to calculate each rule conclusion, which is distinctive when utilizing triangular allotments. Fig. 7 represents the neuro-fuzzy scheme for an illustration with two variables (x 1 , x 2 ) and one output variable (y). In the principal phase of the neuro-fuzzy scheme, the two inputs are categorized into philological values by the set of Gaussian membership functions recognized to every variable. The second stage computes every tenet R (l) its separate enactment degree. Last, the derivation system weights every guideline conclusion ? (l) , instated by the group based algorithm, utilizing the enactment degree computed in the second stage. As mentioned the error signal among the model inferred value Yand the particular measured value (or teaching value) y', is employed by the gradient descent scheme to regulate each rule conclusion. Also he algorithm adjusts the values of ? (l) to diminish an objective function E typically expressed by the mean quadratic error (12). In this equation, the value y'(k) is the preferred output value correlated with the condition vector x'(k) = (x 1 ', x 2 ',?xm'). The element Y(x'(k)) is the conditional response to the same condition vector x'(k) and calculated by Eq. ( 14). ?? = 1 2 [??(??(??) ? ??(??)]] 2 ??????(??)? = ? ?? ?? ?? ?? =1 (??) ?? ?? =1 ???? ?? (??)?? ð??"ð??" (??) ð??"ð??" ??=1 (??) ? ?? ?? ?? ?? =1 (1) ?? ?? =1 ???? ?? (??)?? ð??"ð??" ??=1(14) Eq. ( 15) establishes adjustment for every conclusion ? (l) with the aid of gradient-descent method. Hare the symbol ? is the learning rate parameter, and also t designates the number of learning iterations that are executed by the algorithm. ð??"ð??" (??) (?? + 1) = ð??"ð??" (??) (??)?? ???? ??ð??"ð??" (??) (15) # Global Journal of Researches in Engineering ( ) Volum The inference function Eq. ( 14) depends on ? (l) onlythrough its numerator. The expression composing the numerator is now represented by 'a' and is shown in Eq.( 16). ?? = ? ?? ?? ?? ?? =1 (??) ?? ?? =1 ??? ? ?? (??)?? ð??"ð??" (??) ð??"ð??" ??=1 (??) The denominator of function Eq. ( 14) is dependent on a term d (l) , defined in Eq. ( 17), and denoted by b in Eq. (18). ?? (??) = ? ?? ?? ?? =1 (??) ?? ?? =1 ??? ? ?? (??)?(17)?? = ? (?? (??) ) ð??"ð??" ??=1 (18) In order to compute the adjustment of each conclusion value ? (l) , it is necessary to compute the variation of the objective function E, ?E , in relative to the disparity that occurred in ? (l) in the anterior instant, ? ? (l) . So by, using the chain rule to calculate ?E/ ?? (l) results in Eq(19). ?E ??ð??"ð??" (??) ) = ?E ?Y ?Y ?a ?a ??ð??"ð??" (??) ) The utilization of chain rule searches for the term contained in E that is straight forwardly reliant on the quality to be balanced, i.e., the conclusion value ? (l) . Therefore, we can verify by chain Eq. ( 19) that it starts with E dependent of Y value, and it basically depends on ? term also finally, expression a is a function of ? (l) . Now after few moment, the alteration can be done for ? (l) and also be interpreted that is propotional to the error that is between the neuro-fuzzy model response and the supervising value, but it can be weighted by the influence of rule (1), indicated by d (l) , to the final neurofuzzy inference. ð??"ð??" ?? (?? + 1) = ð??"ð??" (??) (??) ? ?? (??(?? ? (??)??? ? (??))?? (??) ? (?? (??) ð??"ð??" ??=1 ) (20) Next, a convergence theorem has been developed to assurance the stability of particular learning algorithm employed for the above-mentioned FNN [10]. A Lyapanov energy function is defined as follows: ?? ?? = ?? ?? = 1 2 ?? ?? 2(21) From Eq. ( 19), we can get ??? = ?? ??+1 ? ?? ?? = 1 2 (?? ??+1 2 ? ?? ?? 2 )(22) The error difference, Î?"E k , can be defined as ??? ?? = ?? ??+1 ? ?? ?? = ???? ?? ??ð??"ð??" ?ð??"ð??" (23) ?ð??"ð??" = ð??"ð??" ??+1 ? ð??"ð??" ?? = ????? ?? ???? ?? ??ð??"ð??"(24) Using Eq. ( 22), we can get ??? = 1 2 (?? ??+1 ? ?? ?? )(?? ??+1 + ?? ?? ) = 1 2 (??? ?? )(2?? ?? + ??? ?? )(25) Substituting eq. ( 24) into eq. ( 23), we have: ??? = 1 2 ???? ?? ??ð??"ð??" ???? ?? ???? ?? ??ð??"ð??" (?2?? ?? + ???? ?? ??ð??"ð??" ???? ?? ???? ?? ??ð??"ð??" ) = 1 2 (?? ?? ???? ?? ?? ð??"ð??" ) 2 [? ???? ?? ?? ð??"ð??" ) 2 ?? 2 ? 2???(26) If ??? < 0, the convergence of the algorithm described in eq.( 26) can be guaranteed. Therefore, we have: ( ???? ?? ??ð??"ð??" ) 2 ?? 2 ? 2?? < 0(27) From eq. ( ), we can obtain: 0 < ?? < 2 ( ?? ?? ?? ??ð??"ð??" ) 2(28) V . # Simulation Results To evaluate the viability of the proposed control strategy, the system is simulated in SIMULINK/ SIMPOWER for a period of 100sec of real and reactive load profiles. The decision of the DC-bus voltage depends on the output voltage of the inverter required which must provide the grid voltage. The association among the DC link voltage Vdc and the line-to-line RMS grid voltage ?? ????,???? , where ?? ?? is the modulation index in the linear region, is specified in Eq. ( 29) [9]. ?? ??ð??"ð??" ? 1.633 ?? ?? ?? ????,???? + ????????????????????????(29) The proposed control strategies are realized using MATLAB/SIMULINK environment by means of the parameters given in Table . 1 [5]. In order to show the response of the power control strategy during the unbalanced voltage condition, another simulation results have been extracted. In this case, The proposed control strategy has been inspected in case of unbalanced voltage conditions. An unbalanced voltage, resulting from unbalanced load, is applied at the grid side. The unbalanced voltage starts at 1.2 sec for duration of 2 sec. The grid voltage during unbalanced voltage has been shown in Fig. 8. The DC link voltage is shown in Figs.13. During the unbalanced voltage, there is an increase on DC linkvoltage but it is not much more than 10% of nominal value. In these conditions, to stabilize the dc-linkpower, the neuro-fuzzy controller manages the power flow between power sources. # Conclusion This paper presents modeling, control and power control in in a grid connected PV/Fuel Cell/Battery hybrid power generation system in a microgrid. Here SIMULINK/SIMPOWER was utilized to to model the system and simulate a power flow control strategy. PV, fuel cell and and battery subsystems with power electronic converters are modeled. In addition, to disseminate the power between power sources, the neuro-fuzzy power controller has been produced to settle the DC-link power. Our simulation results are shown to exhibit the viability and ability of proposed control strategy amid various operating conditions in utility grid. 1![Fig-1 : Block Diagram of Grid connected Hybrid system a) Proton Exchange Membrane Fuel Cell Model A fuel cell functions like a battery by transforming the chemical energy into electrical energy, but it differs from a battery in that as long as the hydrogen and oxygen is supplied it will produce DC electricity continuously. Fuel cells play a vital role in distributed generation because of their advantages such as high efficiency, no contamination gasses and particular structure adaptability. The Nernst's mathematical statement and Ohm's law oversee the normal voltage size of the fuel cell stack. The subsequent mathematical statement representations the voltage of the fuel cell stack:](image-2.png "Fig- 1 :") 2![Fig.2 : Equivalent solar cell model with Rs](image-3.png "Fig. 2 :") 3![Fig. 3 : Battery Model e) DC/DC Boost Converter Model The role of boost DC/DC converters is to deliver power to the user in appropriate form at high efficiency. Generally the Power electronic converters are desirable in PV and fuel cell systems to transform DC voltage to the prerequisite values. Fig. 4 shows the DC/DC converter model.](image-4.png "Fig. 3 :") 4![Fig. 4 : Boost DC/DC Converter Model.](image-5.png "Fig. 4 :") 6![Fig. 6 : Overall system control structure. Hence the control strategy also retains the DClink/battery voltage within a band about the nominal DClink voltage to have the inverter in synchronism with the grid. The following differential equation for DC link power balance is given:](image-6.png "Fig. 6 :") 8![Fig.8 : The unbalanced grid voltages Fig. 9 shows the 3-phase voltages; which are well regulated underthe unbalanced disturbance.](image-7.png "Fig. 8 :") 9![Regulated grid voltageFrom Fig.10and Fig.11the gid voltage and cuurents under the neuro-fuzzy controller have great improvement in case of Total harmonic distortion THD that was depicted in Fig.12.](image-8.png "Fi g. 9 :") 10![Fig. 10](image-9.png "Fig. 10") 11![Fig. 11From Fig.12THD is only 3% in case of neurofuzzy controller and 15.27% in case of without neurofuzzy controller.](image-10.png "Fig. 11 From") 12![Fig. 12](image-11.png "Fig. 12") 13![Fig.13 : DC-link voltage VI.](image-12.png "Fig. 13 :") .Control Strategies of a Neuro-Fuzzy Controlled Grid Connected Hybrid PV/PEMFC/BatterySOC,%70DC/AC Converter ParametersNominal AC Voltage400 VNominal Phase Current125 ANominal DC Voltage (16) R s720 V 0.9 m?Ls0.01mHFs (Hz)50 HzDC/DC Converter ParametersRated Voltage ( V)200/650VResistance (R) Rated Power Capacitance (C)2.3? 50KW 1.5 mfYear 2016Inductor (L)415 µH9(19)XVI Issue II Version IPEMFC Parameters Faraday's Constant Hydrogen time constant(t H2 ) Hydrogen valve molar constant(K H2 )96487000 C/kmol 26.1 sec 8.43x10 -4 ValuesGlobal Journal of Researches in Engineering F ( ) Vol umeKr Constant=N0/4F9.9497x10 -7No load voltage (E 0 )0.6VNumber of cells (N 0 )384Oxygen time constant (t O2 )2.91 secBattery Model parametersValuesMaximum allowable terminal voltage730 VMinimum allowable terminal voltage710 VOperating Terminal Voltage725 V © 2016 Global Journals Inc. (US) * Integration issues of distributed generation in distribution grids" proceedinds of IEEE EdwardJCoster JohannaM AMyrzik BasKruimer LWil Kling January 2011 * Optimization of Photovoltaic Power Systems DjamilaRekioua ErnestMatagne Springers 2012 * Modelling and Control of Fuel Cells MHashem Nehrir CWang 2009 Wiley * Design of primary and secondary cells II. Anequation describing battery discharge CMShepherd J. Electrochem. Soc 112 657 1965 * Intelligent control of hybrid photo voltaic/fuel cell/energy storage power generation system AminHajizadeh SamsonGTesfahunegn ToreMUndeland 10.1063/1.3618743 J. Renewable Sustainable Energy 3 43112 2011 * Maximum photovoltaic power tracking: an algorithm for rapidly changing atmospheric conditions KHussein IMuta THoshino MOsakada IEE Journal of Generation, Transmission and Distribution 142 59 1995 * Sliding mode control of boostand buck-boost power converters using method ofstable system centre YuriBShtessel YBZinober AShkolnikov A Automatica 39 2005 * Intelligent Robust Control of Hybrid Distributed Generation System under Voltage Sag AHajizadeh MAGolkar Journal of Expert Systems with Applications 37 2010 * NMohan TMUndeland WPRobbins 2003 Wiley * Operation of Stand Alone PV Generating System for Supplying Unbalanced AC Loads RNoroozian MAbedi GBGharehpetian SHHosseini Iranian Journal of Electrical & Electronic Engineering (IJEEE) 3 2007