Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC in a grid/charging station connected hybrid power system

This paper focuses on the indirect adaptive tracking control of renewable energy sources in a grid-connected hybrid power system. The renewable energy systems have low efficiency and intermittent nature due to unpredictable meteorological conditions. The domestic load and the conventional charging stations behave in an uncertain manner. To operate the renewable energy sources efficiently for harvesting maximum power, instantaneous nonlinear dynamics should be captured online. A Chebyshev-wavelet embedded NeuroFuzzy indirect adaptive MPPT (maximum power point tracking) control paradigm is proposed for variable speed wind turbine-permanent synchronous generator (VSWT-PMSG). A Hermite-wavelet incorporated NeuroFuzzy indirect adaptive MPPT control strategy for photovoltaic (PV) system to extract maximum power and indirect adaptive tracking control scheme for Solid Oxide Fuel Cell (SOFC) is developed. A comprehensive simulation test-bed for a grid-connected hybrid power system is developed in Matlab/Simulink. The robustness of the suggested indirect adaptive control paradigms are evaluated through simulation results in a grid-connected hybrid power system test-bed by comparison with conventional and intelligent control techniques. The simulation results validate the effectiveness of the proposed control paradigms.


Introduction
The global warming and environmental deterioration are considered as a factor of supply and consumption of energy based on the fossil fuels. The fastest growing electrical energy demand, continuously diminishing fossil fuel energy sources and pollution crises bring renewed interest in renewable energy. Renewable energy is clean, sustainable and inexhaustible.
Among numerous renewable energy sources, wind energy conversion system (WECS) and PV energy conversion system (PVECS) are the fastest growing energy sources. In WECS, the wind turbine blades are used to harvest the kinetic energy of wind. The kinetic energy is PLOS ONE | https://doi.org/10.1371/journal.pone.0183750 September 6, 2017 1 / 32 a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 converted to electrical power by using an appropriate generator. WECS can be categorized into fixed-speed generation and variable-speed generation. The variable-speed WECS has several advantages as compared to fixed-speed WECS in-terms of operation, efficiency and power quality [1]. The use of a direct-drive permanent magnet synchronous generator (PMSG) with a wind turbine, significantly enhances the reliability of the VSWT-PMSG, because, the absence of magnetizing current in PMSG allows it to operate at high power factor [2]. The PVECS converts the solar energy into electrical energy. The power generation capability of WECS and/or PVECS is heavily dependent on local weather patterns, i.e., wind speed, irradiance and temperature. To improve the output efficiency of WECS and PVECS, it is crucial to operate WECS and PVECS near maximum power point (MPP). For WECS, Optimal torque control (OTC) MPPT strategy is quite popular, as it extracts the optimum wind energy [3], [4]. In OTC MPPT method, the torque of PMSG at the given wind speed is amended on the basis of the maximum power reference torque. Although, the OTC MPPT strategy is simple and fast but the efficiency of this method is low, because, the wind speed is not measured directly. Therefore, the wind fluctuations cannot be captured substantially and abruptly on the reference signal [5]. Tip speed ratio (TSR) is another commonly used MPPT technique for WECS. The optimal TSR remains constant irrespective of wind speed and guarantees that the maximum power is extracted [6]. Although, TSR MPPT scheme is simple as it continuously and directly measures the wind speed but in reality, measuring the accurate wind speed becomes challenging and also, increases the cost of the system. In power signal feedback (PSF) MPPT for WECS, the wind turbine reference optimum power curve is obtained through experimental results. Then, the data of MPPs versus wind speeds are used from a lookup table [7]. However, the PSF MPPT scheme is complex to implement. In perturbation and observation (P&O) or hill climbing (HC) MPPT for WECS, the control variable is perturbed with a small step-size and the subsequent variations are observed in the objective function till the slope becomes zero [8]. The main disadvantage of P&O or HC MPPT scheme is that the indistinct difference between the powers results in incorrect decision in defining the direction for next step. The intelligent MPPT controls for WECS are also reported in the literature. The intelligent MPPT controls include artificial neural network (ANN) [9], fuzzy logic [10] and NeuroFuzzy [11]. Though, the neural network has the strong capability of learning but it needs a tremendous amount of training data and requires a long time to train the network. The fuzzy logic controllers have high convergence speed and acceptance of ill-defined signals but their design depends upon trial-and-error method. The NeuroFuzzy MPPT is successfully implemented to acquire the maximum power from WECS. However, the inherent shortcoming of the NeuroFuzzy system is that it has long computational time and becomes trapped in local minima of the search space [12].
There are two main categorizes of MPPT control schemes for PV named as conventional and soft computing. Conventional MPPT control schemes include P&O, HC and incremental conductance (IC), whereas, the soft computing MPPT include ANN, fuzzy logic and evolutionary algorithms. The efficiency of the P&O MPPT is improved by subsiding the steady-state oscillation and reducing the probability of losing its tracking direction [13]. However, a boundary condition is applied to confine the MPP voltage in the range of the MPP locus. The HC MPPT control scheme is successfully applied for stand-alone parallel PV power generation [14]. However, the dynamic performance of HC suffers the divergence of MPP. The IC MPPT control scheme has a fast and accurate response but the perturbation step and marginal error are computed on the basis of trial and error approach [15]. The ANN based MPPT method is implemented to spontaneously identify the global MPP based on the preselected number of power measurements [16]. Although, the results yield the robustness of the proposed MPPT control scheme but at the cost of long computational time. In [17], two types of fuzzy logic controllers are implemented to track the PV MPP. One is an adaptive fuzzy logic controller and other is the conventional fuzzy logic controller. However, the adaptive fuzzy logic controller still requires fine tuning via trial and error, because, it exhibits the same behavior as that of the conventional fuzzy logic controller. A wavelet-based NeuroFuzzy control scheme is used to track the MPP of PV [18], [19]. Wavelets are introduced in the structure of NeuroFuzzy to improve the performance. The simulation results reveal high efficiency and fast response of wavelet-based NeuroFuzzy control scheme. The most popular evolutionary algorithms used to track the MPP of PV are the genetic algorithm (GA) [20], particle swarm optimization (PSO) [21] and differential evolution (DE) [22]. Although, evolutionary algorithms are the stochastic methods which are quite efficient to solve a real valued, nonlinear, multi-modal optimization problems. But the appropriate selection of control parameters, initial values, solution archive and locality of the search space are still the potential areas of concern.
SOFC is versatile, efficient and alternative energy source which generates electrical power directly from hydrocarbon fuels at 800-1000˚C. The limitation of SOFC is that it has dynamic load following issue. In the case of large power variations, the hydrogen starvation occurs owing to the slow fuel supply which results in voltage drop, anode oxidation, and catalyst corrosion. In literature, the load following issue of SOFC is addressed by using two different control strategies. In one control strategy, the input hydrogen which is directly proportional to the SOFC stack current is controlled to resolve the load following issue [23]. In the other control strategy, the SOFC terminal voltage is maintained at a constant value to get the swift response of SOFC [24]. To address the load following issue of SOFC, different conventional and advanced control strategies are available. A feedforward controller is designed to resolve the fuel starvation problem but the system is not very efficient due to the weaknesses of feedforward controller, i.e., sensitive to external disturbances and the occurrence of steady state error [25]. In [26], a multi-loop feedforward/feedback control scheme is presented. Although, the plant is well stable under a tight linear region of operation but the overall system becomes complex and complicated. A master control PID feedback approach is adopted for load following but large load variations need more effective control scheme [27]. Model predictive control (MPC) is a predictive model and receding horizon optimization based feedback control system. This is an attractive approach for SOFC, because, a wide-range constraints of input/output variables of a nonlinear system are directly handled [28]. However, MPC is complex and computationally slow to get a reliable and accurate first-principles prediction model for a large-scale nonlinear system. Likewise, the advanced control methods are also used for the SOFC system which includes ANN predictive control [29] and fuzzy logic predictive control [30]. The effectiveness and stability of ANN predictive control and fuzzy logic predictive control is satisfactory. However, the resulting controllers are too convoluted and increase the cost of the system for hardware implementation. The hybrid power system (HPS) is an emerging power generation scheme, because, the integration of renewable energy sources along with the storage systems improve the efficiency and energy supply, and its environmental and economic sustainability. A PV and SOFC along with electrolyzer based HPS is developed for supplying the electricity to residential load [31]. However, the stated HPS has small scale application. The dynamic operation and control of PV, wind and SOFC based HPS is presented [32]. A proficient power sharing method for all energy sources is also presented. However, the system is unable to store the surplus power produced by renewable energy sources. The dynamic performance of a stand-alone wind-solarbattery based HPS is investigated [33]. However, the HPS is not connected to the utility grid. A stand-alone PV, fuel cell and ultra-capacitor based HPS for the residential load is presented in [34]. The sizing and designing of PV, fuel cell and ultra-capacitor are also addressed to enhance the performance of the HPS. But the application of stated HPS is stand-alone. The dynamic modeling and control of a grid-connected wind-PV-battery based HPS is designed with versatile power transfer [35].The HPS is stable and has the dispatch-ability to provide the power to the grid. However, a small scale application of HPS is presented. The solar, wind and diesel engine based stand-alone HPS is used to meet the load [36]. However, the quality of power is not addressed. A novel operation and control strategy for a stand-alone HPS with battery storage is presented [37]. Different wind and load conditions are investigated. The optimum power from wind is extracted using machine side converter. A boost converter is used to control fuel cell and a buck converter is used to control electrolyzer. The stated HPS prevents the blackouts. However, the stand alone application of HPS is offered. A PV-FC-WES based grid-integrated HPS is used to meet the load [38]. The P&O algorithm is used to extract the MPP from PV. A pitch angle controller is used to extract the MPP from wind turbine. However, the load power is still unmanageable, therefore, local power shedding is used. The stochastic energy management is applied for a PV power based smart home with plug-in electric vehicle (PEV) to minimize the consumer's energy charges [39]. By using the stochastic dynamic programming (SDP), home demand as well as the PEV charging requirements are satisfied. The modeling of PEV energy storage is reported for the first time. However, only a single home is considered. A PV array, battery storage system and the utility grid are used to provide the power to a building's electrical loads [40]. The PV and battery are coupled to a DC bus. A DC/AC inverter connects the DC bus to AC loads and with the utility grid. A nonlinear predictive energy management strategy is used for optimal power flow. The energy management strategy uses real-time forecasted load, weather conditions and electricity cost. However, wind turbine and SOFC with appropriate models need to integrate in the system.
To address all the afore mentioned hitches, an efficient Chebyshev wavelet embedded Neu-roFuzzy indirect adaptive MPPT control scheme for VSWT-PMSG, an efficient Hermite wavelet embedded NeuroFuzzy indirect adaptive MPPT control scheme for PV system and an effective Hermite wavelet incorporated NeuroFuzzy indirect adaptive control scheme for SOFC system integrated into a grid-connected HPS are proposed. In order to supply uninterrupted and consistent power to the load, a supervisory control policy which consists of nine different modes of operation is also presented. In the stated HPS, tracking the MPP for VSWT-PMSG system, tracking the MPP of PV system and getting the swift response of SOFC are quite difficult, because, this system is highly characterized by nonlinearity. The nonlinearity arises due to the erratic load, variable wind speed, dynamic solar radiation and inconsistent temperature. This work is actually an extension of the work presented in [12], [18], [19]. The rest of the paper is organized into five main sections. Section 2 presents the problem formulation. Section 3 presents the mathematical modeling for VSWT-PMSG, PV and SOFC system. Section 4 gives the details of operation strategy and supervisory control. Simulation results are discussed in section 5. Section 6 concludes the outcomes of this research work.

Problem formulation
The nonlinear MIMO hybrid power system with renewable energy sources is shown in Fig 1. The auto-regression NeuroFuzzy model of the nonlinear dynamic HPS can be described as: Where Ω(k) = y(k-1), . . ., y(k-n), u(k), u(k-1), . . ., u(k-m). ŷ NF-WT (k) gives predictive output at time step k for a SISO VSWT system. ŷ NF-PV (k) gives predictive output at time step k for a SISO PV system. ŷ NF-SOFC (k) gives predictive output at time step k for a SISO SOFC system. The nonlinear dynamic models for VSWT system, PV system and SOFC system is captured online if: Grid stability ) The μTurbine, battery storage system, super capacitor and the charging station are interfaced to either DC bus or AC bus via PID-PWM based converters. The DC and AC buses are connected through the main inverter which are controlled by a PID-hysteresis based PWM technique. All these devices are operated by the supervisory control system. They are made to follow their respective reference power signals computed by the supervisory control system. In a guaranteed stable HPS, the change in all the system variables (states) diminishes asymptotically as the time goes to infinity. For the safe and stable operation of the HPS, it is extremely important to supply quality and reliable power to the load. It will also facilitate the smooth and uninterrupted bidirectional flow of power between the grid and HPS.

VSWT-PMSG MPPT subsystem adaptive control design
For the wind-turbine subsystem control strategy design, the objective is to acquire the MPP of the wind-turbine subsystem through the indirect adaptive Chebyshev wavelet based Neuro-Fuzzy algorithm. The unknown subsystem model, f NF-WT (Ω(K)) is identified online as defined in (1) and depicted in Fig 2. The nonlinear subsystem dynamic model, f NF-WT (Ω(K)) for VSWT-PMSG is identified by using the objective function Where V wind is the wind speed, V ci is the cut-in wind speed, V co is the cut-out wind speed, V R is the reference wind speed, η gears is the gearbox efficiency, η gen is the generator efficiency, N WT is the capacity factor, P R is the rated power, A S is the swept area, ρ is the air density, C p is the power coefficient. C p is the function of TSR λ and the blade pitch angle β. Based on X WT -Ide, the parameters, ξ ij {m ij , σ ij , w ij } of the Chebyshev wavelet based NeuroFuzzy model, f NF-WT (Ω(K)) are optimized adaptively. Where m ij , σ ij and w ij are the mean, variance and weight of the Gaussian membership function for ith input, jth rule. Where i = 1, . . ., n and i = 1, . . ., m. The update equations for all parameters are: Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC Where α WT is the learning rate and e I = I NF-WT (k)-Î NF-WT (k) is the identification error.
WhereG m are Chebyshev polynomials and can be calculated as: The cost function for the controller is: The control law u WT (k) is: The generalized update equation for control law u WT (k) is given as: The e c = I NF-WT (k)-I NF-WT-ref (k) is used to optimize the parameters χ ij {κ ij , ϑ ij , ν ij } of the controller. Where Δe c (k + 1) = e c (k)−e c (k−1), η WT is the learning rate, κ ij , ϑ ij and ν ij are the mean, variance and weight of the Gaussian membership function for ith input, jth rule. The term ℧ WT (k) can be calculated as: Where @℧ WT ðkÞ @w ij ðkÞ can be simplified as: The term @Î NFÀ WT ðkÞ @u WT ðkÞ can be calculated as: Where c i 11 and c i 12 are wavelet coefficients for ith input. The MPPT of VSWT-PMSG nonlinear subsystem is guaranteed by the successful convergence of the variables given below through the implementation of the proposed adaptive control algorithm. Where T e is the electrical torque, T m is the mechanical torque and T ref is the reference torque for VSWT-PMSG. The proposed intelligent adaptive control law, u WT will operate the subsystem via AC/DC/DC/DC converter to generate maximum power based on the available wind speed. It ensures the maximum power acquired for a safe wind speed profile (between cut-in and cut-out speeds).

PV MPPT subsystem adaptive control design
The objective of the photovoltaic subsystem control strategy is to track the maximum operating point through the indirect adaptive Hermite wavelet based NeuroFuzzy algorithm. The unknown subsystem model, f NF-PV (Ω(K)) is identified online as given in (1)  The nonlinear subsystem dynamic model, f NF-PV (Ω(K)) for PV system is identified by using the objective function s.t. (PV subsystem model constraints) Where P peak is the peak power, V oc is the open circuit voltage, I sc is the short circuit current, φ max is the maximum solar irradiance and T max is the maximum temperature. They are specified under standard test conditions (STC). When the change of PV power with respect to the operating voltage is zero, then the MPP is achieved as follows: Where s is the slope of PV power with respect to the operating voltage. At MPP, the dynamic input conductance of the subsystem model is equal to the negative of the input static conductance. The change in duty cycle also alters the impedance seen at the input of a DC/DC boost converter. For boost converter, the load impedance (output impedance) is always greater than the input impedance. Based on X PV -Ide , the parameters, z ij {ƛ ij ,γ ij ,ϖ ij } of the Hermite wavelet based NeuroFuzzy model, f NF-PV (Ω(K)) are optimized adaptively. The Hermite wavelet ψ a,b (x i ) is defined on the interval [0, 1] as follows: Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC Where x i is the ith input, c = 1, 2, . . ., c n is the level of resolution, a = 1, 2, . . ., 2 c-1 , â = 2a-1 is the translation parameter, and b = 1, 2, . . ., B-1 is the order of the polynomial, B>0. Where H b is orthogonal pertaining to the weight function as: The Hermite polynomial H b of order b is defined on the interval [-1,1] and is given as The update equations for all parameters are: Where α PV is the learning rate and e IÀ PV ¼ sðkÞ ÀŝðkÞ is the identification error. ƛ ij , γ ij and ϖ ij are the mean, variance and weight of Gaussian membership function for ith input, jth rule, i.e., The cost function for the controller is: The control law u pv (k) is: The e c-PV = s(k)-r(k) is used to optimize the parameters ℵ ij {} ij , ϛ ij , τ ij } of control law u pv (k) . The generalized update equation is given as follows: Where Δe c -PV (k+1) = Δe c -PV (k)-Δe c -PV (k-1), η PV is the learning rate, } ij , ϛ ij , and τ ij are the mean, variance and weight of the Gaussian membership function for ith input, jth rule. The term ℧ PV is calculated as: Where @℧ PV ðkÞ @@ ij ðkÞ can be simplified using the following equation: @℧ PV ðkÞ @@ ij ðkÞ ¼ e cÀ PV ðkÞ @ŝðkÞ @u PV ðkÞ À Z PV u PV ðkÞ ! @u PV ðkÞ @@ ij ðkÞ ð30Þ The term @ŝðkÞ @u PV ðkÞ can be calculated as follows: Where q i 11 and q i 12 are Hermite wavelet coefficients.
The MPPT of PV subsystem is accomplished by the successful convergence of the variables given below through the implementation of the proposed adaptive control algorithm. The proposed intelligent adaptive control law, u PV will operate the subsystem via DC/DC converter to extract maximum power according to the irradiance and temperature. The maximum power extraction is guaranteed under the unpredictable climatic conditions.

SOFC adaptive control problem
The SOFCs have sluggish dynamic response as compared to the dynamic responses of the power conditioner and the load. The SOFC incapability to alter its electrical output, i.e., current as swiftly as the electrical load variation has significant consequences on the overall HPS. Therefore, the dynamic response of SOFC has substantial importance. For the SOFC swift response, the quantity of input hydrogen flow rate must be controlled. The input hydrogen flow rate is proportional to the SOFC stack current. Therefore, the optimal flow rate of input hydrogen is achieved by controlling the SOFC stack current. The relationship for the SOFC stack current is given by: Where H uti 2 is the optimal hydrogen utilization, and m in H 2 is the molar flow of input hydrogen. H uti 2 has a typical range of 80-90%. For optimal hydrogen utilization, the SOFC current lies within the following limits: Where ' is the constant that gives the amount of hydrogen reacting in the SOFC, and 0:8m in H 2 and 0:9m in H 2 are the minimum and maximum limits of molar flow of hydrogen, respectively. I SOFC-min I SOFC-ref and I SOFC-max are the minimum reference and maximum SOFC currents, respectively. The SOFC power demand is converted into the current as follows: The unknown subsystem model, f NF-SOFC (Ω(K)) is identified online as defined in (1). The SOFC-electrolyzer subsystem along with hydrogen storage tank is depicted by Fig 4. The nonlinear subsystem dynamic model, f NF-SOFC (Ω(K)) for SOFC is identified by using the objective function
The SOFC nonlinear subsystem convergence is guaranteed if the following variables converge successfully using the proposed adaptive control algorithm. The proposed intelligent adaptive control law, u SOFC (k) will utilize the subsystem via DC/ DC converter to deliver hydrogen swiftly and consequently the output current, I NF-SOFC (k). The control law u SOFC (k) is: 3 Mathematical modeling

Variable Speed Wind Energy Conversion System (VS-WECS) modeling
3.1.1 Wind turbine model for maximum power extraction. The aero dynamic power P a caught by the wind turbine is: Where, V co > V wind > V ci . λ is the ratio of the tip speed of the turbine blades to the wind speed is l ¼ Ro r V wind . Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC The rotor power (aerodynamic power) is also defined by: The aerodynamic torque T a is: The T a drives the wind turbine at the ω r (speed). For a specific wind speed, the output power is proportional to the rotor speed, i.e., P a = kω r 2 . Where k ¼ 0:5rAC p R l À Á 3 and the optimum aerodynamic rotor power is captured by controlling the rotor speed ω r . For a certain wind speed, the optimum power is: Where k opt ¼ 0:5rAC p;opt R l opt 3 .
The optimum power acquisition from the wind refers to acquiring the necessary power under fluctuating wind speed condition.

Modeling of photovoltaic subsystem.
A photovoltaic cell converts photon energy directly into electrical energy in the form of direct current which means that a photovoltaic device model must be based on the electrical characteristics. The average photovoltaic array model consists of two controlled current sources I pv-1 and I pv-2 . The current source I pv-1 is controlled by I L .
The controlled current source I pv-2 depends upon I L2 .
Where s = irradiance, I d = diode current, v d = diode voltage, I 0 = diode saturation current, v T = temperature voltage = kT cellÀ K qnIN cell N ser , T cell-K = cell temperature, k = Boltzman constant = 1.3806×10 −23 Jk −1 , q = electron constant = 1.6022×10 −19 C, nI = diode identity factor, N cell = no. of series connected cells per module, N ser = no. of series connected modules per string, N par = no. of parallel connected modules per string. The output voltage of the PV array V PV is given as: Where v d is the voltage across R sh and R sÀ array ¼ R s N ser N par . The output power of the PV array at the PV terminals is given as: Solid Oxide Fuel Cell Energy Conversion System (SOFCECS) modeling Fig 5. The molar flow of hydrogen H 2 through the valve is proportional to its partial pressure, and is given as:

Modeling of SOFC subsystem. The SOFC dynamic model is shown in
Where k is the anode valve constant, m H 2 is the molecular mass of hydrogen, q H 2 is the molar flow of hydrogen, p H 2 is the partial pressure of hydrogen and K_H 2 is the valve molar constant of hydrogen. The molar flow of hydrogen that reacts is calculated as: Where N 0 is the number of series cells in the stack, F is the Faraday's constant, K r is the constant and I SOFC is the stack current. The partial pressure of hydrogen P_H 2 yields the following expression: Where, The v anode represents the volume of the anode, R is the universal gas constant and T is the absolute temperature. A similar operation is made to calculate the partial pressures of oxygen and water. The Nernst's equation used to calculate the stack output voltage is represented as: Where E 0 is the voltage associated with the reaction free energy, p_O 2 and p_H 2 O are the partial pressures of oxygen and water. The output voltage of SOFC is given as: where V SOFC is the SOFC output voltage, V Nernst is the Nernst potential, V A is the activation polarization, V 0 is the ohmic polarization, and V C is the concentration polarization.

General description of the operation strategy
The prime obligation of the supervisory control is to ensure continuous and reliable power supply to the load. The stated HPS has two types of dynamic loads, one is residential load and other is charging station load. The total load is the accumulative of both, i.e., load = residential load + charging station load. There are nine different modes of operation and the supervisory control switches between these modes depending upon the operation strategy. These operation modes are determined by the energy balance between the total generation and the total load. The Modes 1-5 deal with the deficient power, i.e., generation is less than the load, whereas, Modes 6-9 deal with the excess power, i.e., generation is greater than the load.

Supervisory control policy
The following control actions take place in every mode of operation.
Where P D_ref = P Load = P residential + P CS = I L × V L , SOC Bat is the state-of-charge of battery and P Bat-d represents the battery in discharge mode. This mode of operation is characterized by both generation subsystems, i.e., WT and PV set to operate at their maximum energy conversion points but the load is not met by the maximum power of WT and PV. The battery is in discharge mode to meet the load. During this mode, only renewable energy sources and battery system deliver the power to the load.
Where SOC SC is the state-of-charge of SC. In this mode of operation, WT and PV subsystems deliver their maximum power but the load still requires the power. The battery is in discharge mode which means the battery is delivering power to the load. The remaining deficient power is delivered by SC. The SOFC and MT are kept OFF and also, there is no need to take Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC the power from the utility grid, because, the load is met by power taken from renewable sources and storage system.
Where P SOFC is the SOFC power. During this mode of operation, the load is met by the cumulative power taken from WT, PV, battery, SC and SOFC. The battery and SC both are in discharge mode but the load requires more power. The SOFC is set to track a power reference. This reference corresponds to the power required to satisfy the total power demand.

Mode 4:
Where P Grid represents the grid power. In this mode, WT, PV, battery, SC and SOFC deliver the power to the load but the load requires more power which is given by the utility grid.

Mode 5:
Where P MT represents the micro-turbine power. In this mode of operation, the utility grid is unavailable due to peak demand hours. The WT and PV subsystem deliver their maximum powers to the load. The battery and SC both are in the discharge mode, SOFC also delivers the power to the load. The MT supplies the remaining deficient power to the load.

Mode 6: If P D ref P Ren
Where P Elect represents the electrolyzer power. In this mode of operation, there is some excess power P Exss in the system. This excess power is utilized by electrolyzer. The battery remains in the discharge mode.

Mode 7:
Where P SC-r represents the SC is in charge mode. During this mode, the renewable energy sources generate more power than load. The battery is in discharge mode but the SOC of SC is less than 90%, so the excess power is utilized to charge SC and then by the electrolyzer.

Mode 8:
P Bat ¼ P BatÀ d P SCÀ c ¼ P Exss P GridÀ r ¼ P Exss À P SCÀ c P Elect ¼ P Exss À P SCÀ c À P GridÀ r 8 > > > > < > > > > : Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC Where P Grid is the electrolyzer power. In this mode, the WT and PV generate more power than load. The SC is in charge mode. The excess power after charging the SC is given to the utility grid. There is still some excess power which is utilized by electrolyzer to produce hydrogen for SOFC.
In this mode of operation, the battery is disconnected, i.e., neither it will charge nor discharge. The generation is greater than load. The SC is in charge mode. The excess power is given to utility and then to the electrolyzer.

Results and discussion
The performance of the stated HPS and proposed controllers is evaluated in MATLAB/Simulink R2015a. In an 11 kV grid-connected HPS, wind generation of 100 kW, PV of 260 kW, SOFC of 200 kW, electrolyzer of 150 kW, and MT of 200 kVA along with backup sources (200 Ah battery and 165 F Super-Capacitor) are modeled for the accumulative dynamic residential and charging station load. Defense Housing Authority (DHA), Islamabad, Pakistan, is taken as a case study. The hourly basis wind speed (m/s), irradiance (W/m 2 ) and ambient temperature (˚C) levels are recorded by the Pakistan Meteorological Department (PMD). The stated HPS with suggested controllers is tested for the 24-hrs period of 22nd June 2015. In this study, a Chebyshev wavelet incorporated NeuroFuzzy indirect adaptive controller (CWNFC) is used to extract the maximum power from WT subsystem. A Hermite wavelet embedded NeuroFuzzy indirect adaptive controller (HWNFC) is used to acquire the maximum power from PV subsystem. A Hermite wavelet incorporated NeuroFuzzy indirect adaptive controller (HWNFC) is used to get the swift response from SOFC subsystem. To evaluate the performance of CWNFC, an adaptive indirect NeuroFuzzy Takagi-Sugeno-Kang controller (ITSKC), an adaptive direct NeuroFuzzy TSK controller (TSKC) and a hysteresis based PI controller (Hyst-PI) are also used to extract the MPP from WT subsystem. For PV subsystem, the performance of HWNFC is compared with ITSKC, TSKC and incremental conductance based PI controller (InCond-PI). The swift response from SOFC subsystem with HWNFC is also compared with ITSKC, TSKC and conventional PI controller.
The MPP current of the WT subsystem is shown in Fig 6 (A). The variable speed wind-turbine is used. The base wind speed is taken as 12 m/s, whereas, wind speed varies between 11 and 14 m/s as shown in Fig 6 (A). For t = 0-4 hrs, wind speed is 14 m/s. During this time period, the current of WT subsystem also increases to 256.4 A. From 4-8 hrs, wind speed drops to 11 m/s and WT current also starts decreasing and reaches 220A at t = 8 hrs. From 8-12 hrs, wind speed again increases to 12 m/s and keeps increasing to 14 m/s till t = 16 hrs. The WT current also changes according to the wind speed, i.e., starts increasing from 8-12 hrs and reaches 226 A at t = 12 hrs and keeps increasing to 249.5 A at t = 16 hrs. In next time interval, i.e., t = 16-20 hrs, again wind speed decreases from 14 m/s to 11 m/s and the WT current also starts decreasing and reaches to 242.8 A at t = 18.53 hrs. In 20-24 hrs, a small change in wind speed, i.e., from 11 m/s to 12 m/s also alters the WT current accordingly. As the WT dynamics cannot change abruptly with the wind speed, therefore, the effect of the change in current is slow. The CWNF, ITSKC, TSKC and Hyst-PI all extract the power from WT. The power extraction using CWNFC is more precise, because, the percentage error using CWNFC is 3.2%, ITSKC is 6.4%, TSKC is 12.4% and Hyst-PI is 42%. The WT MPPT errors with CWNF, ITSKC, TSKC and Hyst-PI are shown in Fig 6 (B). The CWNFC has the minimum error as compared to ITSKC, TSKC and Hyst-PI.
The WT subsystem has two types of torques, one is T m produced by wind-turbine and the other is T e produced by PMSG. The variation in wind speed alters the T m of wind-turbine. The T e of PMSG changes according to the T m of wind-turbine. Fig 7 (A) shows the change in T e with respect to T m . So, T m becomes the reference for T e . Fig 7 (B) shows the change in T m according to reference torque. All the four controllers, i.e., CWNFC, ITSKC, TSKC and Hyst-PI capture both the torques but the CWNFC acquires the torques with the minimum steady error.
In The CWNFC can track with a more stable output WT power as compared to other controllers.
The output power of the PV subsystem depends upon the ambient temperature and irradiance level which are shown in Fig 9. For t = 0-6 hrs, when the sun power is not available, the ambient temperature also decreases and reaches 27˚C at t = 6 hrs. During 5.5-12 hrs, the irradiance level keeps increasing and reaches to 1058 W/m 2 . After t = 12 hrs, the irradiance level starts decreasing and become zero at t = 19 hrs. The ambient temperature reaches a maximum value of 41˚C at t = 15 hrs, but after 15 hrs, the temperature starts decreasing again.
The maximum powers extracted by PV subsystem using HWNFC, ITSKC, TSKC and InCond-PI are shown Fig 10 (A). The HWNFC extracts the PV subsystem output power with steady-state error = 0.4 kW, % overshoot = 0 and % undershoot = -2. The ITSKC extracts the PV subsystem output power with steady-state error = 20.03 kW, % overshoot = 6 and % undershoot = -13. The TSKC extracts the PV subsystem output power with steady-state error = 105 kW, % overshoot = 61 and % undershoot = -92. The InCond-PI extracts the PV subsystem output power with steady-state error = 135 kW, % overshoot = 2 and % undershoot = -8. The HWNFC has minimum steady-state error, % overshoot and % undershoot as shown in The power drawn from the SOFC depends upon the molar flow of input hydrogen. The molar flow of input hydrogen is controlled to get the swift response of SOFC. The HWNFC, ITSKC, TSKC and PI are used to get the swift response from SOFC. In the case of load variations, after a short transient period, the HWNFC quickly achieves the stable condition compared to the other controllers as shown in Fig 12 (A). The ITSKC, TSKC and PI take time and fluctuate more for load variations. The HWNFC has steady-state error = 410 Watts, % overshoot = 0 and % undershoot = -7. The ITSK has steady-state error = 16400 Watts, % overshoot = 21 and % undershoot = -12. The TSK has steady-state error = 49880 Watts, % overshoot = 368 and % undershoot = -4. The PI has steady-state error = 61890 Watts, % overshoot = 12.75 and % undershoot = 0. The HWNFC provides a better control than ITSKC, TSKC and PI as shown in Fig 12 (B).
The performance indexes ITAE, IAE, ISE and ITSE are computed for SOFC. These indexes are computed using e SOFC (t) = P SOFC-ref (t)-P SOFC (t) for HWNFC, ITSKC, TSKC and PI. The HWNFC indexes have least and flattest profile as shown in Fig 13. The HWNFC has a more stable output SOFC power as compared to other controllers.
The reference power is shown in Fig 14 is the load required power. The load is actually the cumulative load of residential load and charging station load, i.e., P Load = P residenitial +P CS . At each hour, the load required power is essentially satisfied by the power extracted from the generating sources used in the HPS. In Fig 14, the P WNFC is the load power when the WT has CWNFC, PV has HWNFC and SOFC has HWNFC, whereas, the P TSKC is the load power when the WT has ITSKC, PV has ITSKC and SOFC has ITSKC, the P TSKC is the load power when the WT has TSKC, PV has TSKC and SOFC has TSKC and, the P PI is the load power when the WT has Hyst-PI, PV has InCond-PI and SOFC has PI. The active power of the load provided by P WNFC has 800 watts steady-state error, whereas, P ITSKC has 10 kWatts, P TSKC has 35 kWatts and P PI has 50 kWatts. The reactive power of the load provided by Q WNFC has 1000 watts steady-state error, whereas, Q ITSKC has 7.2 kWatts, Q TSKC has 25.3 kWatts and Q PI has 36.5 kWatts. Moreover, the active and reactive powers of the load provided by P WNFC and Q WNFC have less fluctuations as compared to other controllers. operation, the load power is met by the power taken from renewable sources and battery storage system. Mode 2 is activated during 0-0.5 hrs, 2-6 hrs and then again 18.5-24 hrs. In this mode, renewable sources, battery and SC deliver the power to the load. In mode 3, renewable sources, battery, SC and SOFC deliver the power to the load. In mode 4, the load power is met by the power taken from renewable sources, battery, SC, SOFC and utility grid. The mode 5 remains active during 0-0.5 hrs, 2.08-6 hrs and 22.25-24 hrs. During this mode of operation, the grid is unavailable due to having peak load hours, so, renewable sources, battery, SC, SOFC and MT deliver the power to the load. The mode 6 is activated during 0.5-2.015 hrs, 2.63 hrs, 2.84 hrs and 6-18 hrs. Mode 6, the HPS has excess power, because, the renewable sources generate more power than the load. This excess power is utilized by electrolyzer to produce In order to ensure the stability of HPS, the net power on DC as well as AC bus must be zero. Fig 16 shows the net power on the DC bus with all four controllers. The ΔP WNFC is the net DC bus power when the WT has CWNFC, PV has HWNFC and SOFC has HWNFC, whereas, the ΔP ITSKC is the net DC bus power when the WT has ITSKC, PV has ITSKC and SOFC has ITSKC, the ΔP TSKC is the net DC bus power when the WT has TSKC, PV has TSKC and SOFC has TSKC and the ΔP PI is the net DC bus power when the WT has Hyst-PI, PV has InCond-PI and SOFC has PI. The net power on DC bus is close to zero in ΔP WNFC case.
Similarly, the net active and reactive powers on AC bus are shown in Fig 17. Among all the four cases, the net active and reactive powers on AC bus are approximately zero in ΔP WNFC and ΔQ WNFC case.
Power quality is the most important aspect of the stated HPS. Power quality is particularly addressed in HPS, which is clearly shown in Fig 18. Total harmonic distortion (THD) for load voltage and load RMS voltage are all in their acceptable limits in all four cases, which ensures that the system is quite stable. The maximum deviations of voltage THD and RMS voltage are quite less in the case of THD WNFC and V RMS-WNFC as compared to other two cases.
The power conversion efficiency for WT, PV and SOFC is the ratio of its output power to the input power, over a specified operating period. The power conversion efficiency for WT, PV and SOFC can be either static or dynamic. The static efficiency refers to the stable conditions of wind speed, irradiance level and load, whereas, the dynamic efficiency is associated with the varying conditions of wind speed, irradiance levels and load. The dynamic efficiency Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC of adaptive MPPT-VSWT-PMSG, MPPT-PV and SOFC tracking control system can be computed using the following equation: Where η 100 is the efficiency and P is P WT-CWNFC , P PV-HWNFC and P SOFC-HWNFC in the case of VSWT-PMSG, PV and SOFC, respectively. The dynamic efficiency for MPPT-VSWT-PMSG, MPPT-PV and SOFC tracking control with all four controllers is shown in Fig 19.
Where P is P WT-CWNFC , P PV-HWNFC and P SOFC-HWNFC in the case of VSWT-PMSG, PV and SOFC, respectively. The rate of change is minimum in VSWT-PMSG power using CWNFC, PV power using HWNFC and SOFC power using HWNFC. The rate of change of the powers extracted by the proposed control paradigms remains minimum during simulation. It clearly Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC indicates that the converters switches will be subjected to minimum stress in case of the proposed adaptive control schemes. Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC

Conclusion
This paper presents the integration of VSWT-PMSG, PV, SOFC, electrolyzer, MT along with storage systems (battery and SC) in a grid-connected HPS to meet the cumulative load of residential load and charging station load. The stated HPS is a large-scaled nonlinear model of various subsystems. In order to improve the performance of the HPS, three different wavelets based indirect adaptive NeuroFuzzy control schemes are used which is the main contribution of this research work. A Chebyshev wavelet embedded NeuroFuzzy indirect adaptive MPPT control scheme is used to extract the optimum power from VSWT-PMSG. A Hermite wavelet embedded NeuroFuzzy indirect adaptive control scheme is used to acquire the optimum power from PV system. Moreover, the load following issue of SOFC is solved using a Hermite wavelet incorporated NeuroFuzzy indirect adaptive control scheme. To maintain a balance between the generation and load in the HPS, a supervisory control policy which consists of nine different modes of operation is used. To validate the performance, the proposed control techniques are tested in a Matlab/Simulink environment. The proposed control paradigms have higher precision than ITSKC, TSKC and PI in terms of efficiency, steady-state error and performance indexes. The effectiveness of the proposed schemes is confirmed by the simulation results. Indirect adaptive soft computing based wavelet-embedded control paradigms for WT/PV/SOFC

Future work
The future work is to replace the conventional PID controllers by adaptive PID controllers used to control the converters more effectively.