Integrating EESS with a wind power system will reduce the power output fluctuation. In order to reduce both short-term and long-term fluctuations, this paper proposes the combination of the NB and EDLC energy storage system to satisfy both high power and energy capacity requirements.
The typical wind power system as illustrated in
Figure 1 is used for simulation studies in this research. In this system, an induction machine-based wind turbine is integrated with the corresponding HESS. A Vestas V82 wind turbine model is provided in the RSCAD library (Winnipeg, MB, Canada), but the simulation models of the NB and EDLC are not given in the RSCAD. Therefore, the RSCAD simulation models of NB and EDLC are developed based on general mathematical models.
Figure 1.
Wind-hybrid energy storage system (HESS) model.
2.1. Development of Electric Double-Layer Capacitor Simulation Model in RSCAD
EDLC is an electrochemical capacitor that consists of two electrodes to allow a potential to be applied across the cell; two double layers are therefore present, one at each electrode/electrolyte interface. Since EDLC has a number of advantages such as high power density, fast charge/discharge process, long lifetime and low temperature performance, it is considered to be a very promising energy storage device used to mitigate short-term fluctuation of wind power output [
17]. For the simulation studies, EDLC is modeled as an ideal voltage source in series with an equivalent resistance, where the value of the voltage source is calculated from the dynamic model, as shown in
Figure 2.
Figure 2.
Equivalent circuit of EDLC simulation model in RSCAD.
Figure 2.
Equivalent circuit of EDLC simulation model in RSCAD.
The dynamic model represents the electrochemical behavior of EDLC. In this paper, a mathematical model consisting of Stern and Tafel equations is used that gives the relationship between voltage, current and the available amount of charge during charging/discharging processes [
18]. The input of this model is the current obtained from the terminal of EDLC, and the outputs are the voltage and the
SOC of EDLC. The detailed block diagram of the dynamic model is shown in
Figure 3.
Figure 3.
Block diagram of the detailed dynamic model.
Figure 3.
Block diagram of the detailed dynamic model.
The amount of electric charge in EDLC (Q) is calculated as Equation (1), where Qinit is the initial amount of charge, i is the measured current, and isd is the self-discharge current:
The saturation limit function is used to restrict the calculated amount of charge within a specific range to prevent the full discharge and overcharge of EDLC. The Tafel equation given by Equation (2) is used to calculate the leakage current isd that represents the self-discharge phenomenon of EDLC:
Finally, the voltage value of the controlled voltage source is obtained from the Stern equation as follows:
where:
Nc number of layers of electrodes;
Ns number of series EDLC cells;
Np number of parallel EDLC cells;
Q electric charge (C);
r molecular radius (m);
ε permittivity of material;
S interfacial area between electrodes and electrolyte (m2);
R ideal gas constant;
T operating temperature (F);
α charge transfer coefficient;
c molar concentration (mol·m−3);
i0 exchange current density;
F faraday constant;
Vmax surge voltage (V);
ΔV over-potential (V).
The additional details of parameter declaration and calculation are given in [
19]. The other required parameters are computed from basic electrochemical equations and some known constants such as Avogadro, ideal gas and Stefan-Boltzmann constants,
etc.In order to represent different initial SOC level of EDLC, the initial amount of charge of EDLC (Qinit) should be calculated before the simulation starts. Since the capacitance of EDLC is nonlinear and dependent on the open circuit voltage, Qinit can be represented by a function of the initial open circuit voltage (Vinit), governed by:
For the purpose of simplifying the simulation model, the above nonlinear voltage-capacitance dependence is rewritten so that the capacitance C is kept constant while Qinit is approximated as follows:
The function
f(
Vinit) is obtained for a specific EDLC module under test, which consists of 3000 F Maxwell BootsCap Devices (San Diego, CA, USA). For this, the simulation is implemented with an EDLC model provided by Simscape module in MATLAB Simulink [
19]. The
Qinit data is captured by changing the initial voltage from 0 V to 2.8 V (equivalent to
SOC level from 0% to 100%) while the EDLC is either fully charged or fully discharged with a constant current. From the simulation data, the relationship between
Vinit and
Qinit is established by curved fitting as follows:
Based on the above mathematic model, a user-defined model (UDM) of EDLC is made by using the component builder in the RSCAD. From the configuration interface of the developed UDM, the parameters of EDLC can be modified to satisfy the requirements of different simulation cases, as shown in
Figure 4.
Figure 4.
EDLC model parameter configuration menu.
Figure 4.
EDLC model parameter configuration menu.
2.2. Development of Ni-MH Battery Simulation Model in RSCAD
The first NB was introduced in 1990 by Sanyo Electric Company in Osaka, Japan. This battery type has a number of advantages such as having high energy density and high efficiency, and containing no toxic metals. Since its response is slower than EDLC, it is used to reduce long-term wind power fluctuation.
NB is also modeled as a controllable voltage source, where the value of the voltage source is calculated from the dynamic model. In this study, the modified Shepherd-based mathematic model of NB presented in [
20] is used. This model includes an equation to describe the electrochemical behavior of a battery directly in terms of terminal voltage, open circuit voltage, internal resistance, discharge current, and
SOC. This model is based on a number of specific assumptions and limitations, including the internal resistance is supposed to be constant and temperature effect is not considered.
The value of the voltage source during the discharge and charge process can be calculated by using Equations (7) and (8), respectively:
where:
Vbatt battery voltage (V);
E0 battery constant voltage (V);
K polarization constant (V/A·h) or polarization resistance (Ω);
Q battery capacity (A·h);
it actual battery charge (A·h);
A exponential zone amplitude (V);
B exponential zone time constant inverse (A·h)−1;
R internal resistance (Ω);
i battery current (A);
i* filtered current (A).
Based on Equations (7) and (8), the UDM of NB is also created by using the component builder module in the RSCAD environment.
Figure 5 shows the main parameter configuration menu of the NB model. The parameters and stack configuration of the battery can be modified through this interface.
Figure 5.
Battery model parameter configuration menu.
Figure 5.
Battery model parameter configuration menu.