Next Article in Journal
How to Meet the Green Deal Objectives—Is It Possible to Obtain 100% RES at the Regional Level in the EU?
Next Article in Special Issue
Three-Phase Grid-Connected Inverter Power Control under Unbalanced Grid Conditions Using a Proportional-Resonant Control Method
Previous Article in Journal
Research on Fuel Cell Fault Diagnosis Based on Genetic Algorithm Optimization of Support Vector Machine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dual-Active-Bridge Model and Control for Supporting Fast Synthetic Inertial Action

Department of Electrical, Electronic and Information Engineering, University of Bologna, 40136 Bologna, Italy
*
Author to whom correspondence should be addressed.
Energies 2022, 15(6), 2295; https://doi.org/10.3390/en15062295
Submission received: 31 January 2022 / Revised: 15 March 2022 / Accepted: 19 March 2022 / Published: 21 March 2022
(This article belongs to the Special Issue Control and Stability of Grid-Connected Power Electronic Converters)

Abstract

:
This article proposes a dual-active-bridge control to support the fast synthetic inertial action in DC microgrids. First of all, the selection of the isolated DC/DC converter to link an energy storage system with the DC bus in a microgrid is analyzed and the advantages of the dual-active-bridge converter controlled by a single-phase shift modulation justify its selection. An active front-end can be then adapted to connect the DC bus with an AC grid. Secondly, this paper presents the design of a discrete PI controller for supporting fast synthetic inertial action. In particular, a discrete dual-active-bridge model based on the transferred power between both converter bridges, which overcomes the approximations of the output current linearization model, is proposed. Moreover, the article introduces a novel equation set to directly and dynamically tune discrete PI parameters to fulfill the design frequency specifications based on the inversion formulae method. In this way, during the voltage/power transients on the DC bus, the controller actively responds and recovers those transients within a grid fundamental cycle. Since the developed set of control equations is very simple, it can be easily implemented by a discrete control algorithm, avoiding the use of offline trial and error procedures which may lead to system instability under large load variations. Finally, the proposed control system is evaluated and validated in PLECS simulations and hardware-in-the-loop tests.

1. Introduction

In conventional power plants, synchronous generators (SGs) are commonly used to convert mechanical power into electrical power. Electrical networks characterized by the presence of these generators are typically referred to as generator-dominated systems. In such networks, SGs slightly accelerate or decelerate during load transients to preserve power balance and system voltage and frequency stability [1]. The dynamics with which these changes of angular speed occur depend strongly on the inertia of the rotor [2]. As a result, whenever there is a sudden disturbance or transient in the grid, such as load or generation connection or disconnection, the variation of kinetic energy stored in the rotor temporarily restores the power imbalances until the intervention of the control on the turbines, which steadily clears out power level mismatches, voltage fluctuations, and frequency deviations. Therefore, sufficiently significant inertia acts as temporary energy storage before appropriate control action takes place [3].
In most cases, SGs operate within conventional fossil fuel-based power plants. However, extensive efforts are being made by the international community to reduce the use of fossil fuels in view of their long-term depletion and because of the known environmental and health issues caused by their combustion [4]. A solution to challenges caused by fossil fuels combustion is to move towards sustainable sources [5]. Such solutions are characterized by distributed generation (DG) systems, based on renewable energy sources (RESs), spreading widely throughout the electrical networks. The large penetration of DG systems has led to the widespread development of the microgrid concept [6]. These microgrids can be connected to the centralized power grid or operated in the islanded mode [7]. With a centralized power grid, microgrids are dominated by the SGs, which provide sufficient inertial support during power imbalances. On the contrary, in the case in which the microgrids are predominant, the diffusion of RES-based DGs, as in the extreme case of the islanded mode, the generator-dominated system becomes an inverter-dominated system (IDS) [8]. This IDS network suffers from low inertia since it has few rotating masses. The stability of such systems is supported by mimicking the synchronous generator rotor activity and providing the inertia virtually using energy storage systems (ESSs), such as batteries, flywheels and supercapacitors coupled with dedicated power converters [9].
To improve the stability of IDS, virtual inertia control (VIC) techniques have been studied and implemented in [10,11]. VIC is proposed mainly to reduce the voltage fluctuations and frequency deviation in the AC bus by the inverter in IDS. This process involves power injection/absorption, mainly to/from the batteries connected with a bidirectional DC-DC converter [12]. It stands clear that the DC bus voltage should be kept constant, permitting optimal inverter operability and achieving frequency and voltage stability in the AC bus using VIC. The control of IDS systems is very complex, and its effectiveness is unequivocally related by one side to the selection of the converter and by the other side to the considered control strategy.
In this paper, the bidirectional DC-DC converter model has been analyzed to achieve a fast response controller. The batteries with the proposed control algorithm for a DC-DC converter can quickly stabilize the DC bus voltage during power imbalances, which later connect to the grid-connected inverter that operates VIC. A similar approach can be adopted to improve the robustness of the DC-link in DC microgrids. In this paper, instead of using a non-isolated bidirectional DC-DC converter, as done in previous literature [13,14], an isolated bidirectional DC-DC converter is employed. Specifically, a dual-active-bridge (DAB) converter is used because it provides isolation through the high-frequency transformer, protecting the batteries during huge load transients or power network collapse, which even leads to avoiding the usage of a low-frequency transformer within the power network. Moreover, DABs can be controlled with different modulation techniques, such as the single-phase-shift (SPS), extended-phase-shift (EPS), dual-phase-shift (DPS), triple-phase-shift (TPS) which already exists in the previous literature [15]. Compared to the SPS, the other modulation techniques improve the degree of freedom, leading to optimizing the converter operation. However, with more degrees of freedom, there will be several control variables needed to be controlled, which makes it difficult to achieve sophisticated calculations and cohesive control [16]. Considering these difficulties with other modulation techniques, the SPS technique has been selected in this paper as it provides some major advantages, such as low switching losses, ease of implementation, high reliability and high efficiency (especially if wide bandgap technology is considered) [17,18]. This work demonstrates the possibility of obtaining performances, in terms of control dynamics, comparable or superior to the systems proposed in [12,19], indicating how the SPS-DAB converter’s fast response can ensure the stability of the network. Therefore, the proposed results can apply to inverter-dominated low and medium-voltage microgrids. The general schematic of the microgrid with isolated DC-DC converters is shown in Figure 1.
Broad research on the SPS-DAB converter for widespread diffusion in various applications has been carried out. However, modeling and stability analysis of the digital controlled DAB converter remains a challenge [20]. In particular, [21] demonstrated that the classical state-space averaging technique is unsuitable for capturing the DAB converter’s small-signal dynamic. This is due to the zero steady-state average value of the transformer current [22]. An improved generalized averaging model (GAM) was presented in [23]. However, this method considers only the first harmonic of the inductor current Fourier series. This approximation leads to a significant steady-state error for the large-signal model [24]. Third-order modeling was proposed in [25] to improve the overall accuracy. However, the complexity of the state-space model renders it challenging to find a direct relation in terms of the circuit parameters and small-signal transfer function useful for control design [26]. An alternative approach called the output current linearization (OCL) model is based on the power flow equations from the leading to lagging bridge. A comparative analysis between GAM and OCL models has been carried out in [27]. In particular, the GAM model presents a steady-state error in the large-signal model, as previously described, a magnitude error, and phase oscillations at the converter’s switching frequency. For these reasons, the OCL model has been selected for the design of the regulator in [27]. Generally, the controller design is achieved in the continuous-time domain meeting specifications on the open-loop transfer function, as the gain crossover frequency and phase margin. Alternatively, the design requirements refer to the settling time and the overshoot of the step response. The discrete-time controller is then achieved using discretization methods [27,28].
This paper presents a discrete DAB model based on the transferred power between both converter bridges, which overcomes the approximations introduced by the OCL model, which may affect the controller tuning process. The main idea is to change the generally adopted regulator control variable, i.e., the PMW phase shift, into the averaged output converter current. In this way, an effortless and linear first-order discrete transfer function of the system to be controlled can be obtained. The discrete controller parameters can be set based on the given design specifications. The PMW phase shift is then computed at the end of each sampling period, applying a proposed nonlinear function.
Based on DAB models, several regulators and control procedures can be found in the literature. However, there are still open issues to be tackled concerning the different application fields of the DAB converter. In particular, the essence of nonlinearity of the DAB converter may cause instability of the closed-loop system under a large load variation [29]. For this reason, easy trial and error procedures for classical PI/PID controllers are not suitable. More complex heuristic approaches, like the genetic algorithm technique, are used to minimize the error functions by applying optimization techniques, see [30]. In this case, the fractional PI controller has been considered. The main advantage of this controller over classical PI regulator is that it has a more degree of freedom, providing a broader range of controllability and more robustness to the system [31]. However, any inappropriate selection of the algorithm parameters, such as the population size, the rate of mutation, and crossover will produce inadequate results, see [32].
In this paper, the selected form of the regulator is the well-known discrete PI controller based on the trapezoidal integration method [33]. This choice is due to the simple structure of the PI controller, its wide use in microgrid controls, and its capability to automatically guarantee zero steady-state error on DC voltage step response. Moreover, discrete formulae for obtaining PI parameters as a function of DAB converter parameters are derived and presented. These parameters are derived through the so-called inversion formulae design method to exactly meet the gain crossover frequency and phase margin design specifications [34]. Since a phase margin equal to 75° guarantees good robustness on DAB converter and control system parameter variations, the PI proportional gain and integral time parameters can be analytically computed to exactly meet the given gain crossover frequency. Since the variation of this specification affects the response dynamic, it could be seen as a design requirement degree of freedom on the base of different operating modes of the DAB converter. For example, during the charging process of a battery, the DAB converter can be controlled to realize constant current (CC) mode, constant power (CP) mode, and constant voltage (CV) mode, while, when the battery works as an energy generator, the DAB converter control challenge is to stabilize undesired oscillations due to load variations promptly. When a large load power variation is generated, the control could be improved by changing the gain crossover frequency requirement. The presented discrete formulae for PI parameters can be easily implemented in a control algorithm for real-time parameters variation.
Firstly, this paper highlights the fast response controller, i.e., to stabilize the DC bus voltage, and the rest of the paper is divided into the following sections: Section 2 discusses the DAB converter with the single-phase shift controller. Section 3 discusses the proposed discrete control model and formulae for PI tuning. The simulation results, comparison with the method proposed in [35] for fast inertial action, and comparison of the proposed discrete control model with a classical pole-placement methodology are discussed in Section 4, respectively. Hardware-In-The-Loop (HIL) results and the comparison with simulation results are discussed in Section 5. Finally, the conclusions and perspectives about future developments of the work are reported in Section 6.

2. Dual-Active-Bridge Converter Connected to Energy Storage

The considered DAB converter topology is shown in Figure 2. Here, the converter input and output voltages represent the battery voltage v b a t and the DC Bus voltage v o u t , while i 1 and i 2 denote the DAB input and output currents, respectively. The output capacitor C is characterized by the voltage v C , the current i C and the equivalent series resistance R C . The two symmetric H-bridges are interconnected by a high-frequency transformer. This is characterized by a turn ratio n , primary and secondary voltages v p and v s , and a primary referred equivalent leakage inductance L , and its current i L . Variables Q i are used to refer to the switching functions of the H-bridge switches. In the considered SPS modulation, the voltages of each bridge are square-wave modulated with a phase shift δ π / 2 : π / 2 between them. The switches operate at a 50% duty ratio and period T s = 1 / f s   ( where f s   is the switching frequency) and following operation modes are described in Table 1. The corresponding switching sequence and the main signals are represented in Figure 3. When δ > 0 , v p leads v s , conversely δ < 0   means that v p lags v s . Note: All the abbreviations are well defined in the text and are also reported at the end of the paper.
The considered DAB converter model hinges on the conservation of the average value of the power P δ , t flowing in the DC input section, high-frequency switching bridge section, and DC output section,
P δ , t = v ¯ b a t t i ¯ 1 δ , t = v ¯ p t i ¯ L δ , t = v ¯ o u t t i ¯ 2 δ , t
neglecting the power dissipation in the converter. Notice that the variables in Equation (1) represent the average values over a switching period of the correspondent signals. The power P δ , t can be computed by the average value of inductor current over a half switching period for symmetry conditions,
i ¯ L δ , t = 1 π 0 δ i L θ d θ + δ π i L θ d θ
where θ = 2 π f s t = ω t . From Equations (1) and (2), it follows that:
P δ , t = v ¯ b a t t v ¯ o u t t ω L n δ t 1 δ t π ,
i ¯ 1 δ , t = v ¯ o u t t ω L n δ t 1 δ t π ,
i ¯ 2 δ , t = v ¯ b a t t ω L n δ t 1 δ t π ,
and the correspondent SPS-DAB model can be represented as in Figure 4 [29,36].
Classical compensator design requires perturbing and linearizing Equation (5) to obtain a linear small-signal model of the DAB converter. For this purpose, the output current linearization (OCL) method can be applied, see [27]. In particular, the phase shift can be rewritten as δ t = φ + δ ^ t , where φ is the phase shift required to transfer the desired power P 2 = v o u t 2 / R at the correspondent operating point, and δ ^ t is the small-signal variable. The phase shift at the steady-state condition can be calculated from Equations (1) and (5) as:
φ = π 2 1 1 P 2 8 n f s L v b a t v o u t   .
In the same way, we can write i ¯ 2 = I 2 + i ^ 2 , where I 2 is the average current around the operation point, while i ^ 2 is the correspondent small-signal variable. Assuming that the first-order terms will be significantly more predominant than second-order terms, Equation (5) can be linearized around the operation point φ , I 2 . It follows that:
i ¯ 2 δ , t = i ¯ 2 φ + δ ^ t i ¯ 2 δ δ | δ = φ = I 2 + δ ^ t φ G i 2 d ,
where,
I 2 = i ¯ 2 φ = v b a t ω L n φ 1 φ π ;   G i 2 d = i ^ 2 δ ^ t = v b a t ω L n 1 2 φ π .
Graphical representations of Equations (5) and (7) are shown in blue and orange lines in Figure 5.
The small-signal transfer function relating the phase angle to the output voltage,
G v d = v ^ o u t s δ ^ s = b s + a ,
where,
a = 1 R C ,   b = v b a t ω L C n 1 2 φ π ,
can be obtained as in [27], neglecting R C . The controller C d z can be obtained by the discretization of Equation (9) and meeting design specifications in terms of phase margin and gain crossover frequency of the open-loop transfer function, see Figure 6.

3. The Proposed Discrete Control Method

This paper proposes a digital technique that overcomes the approximations introduced by the inductor current linearization and the consequent use of the small-signal DAB model in the controller tuning process. The main idea is to change the system control variable from the PMW phase shift δ , see Figure 6, to the mean value of the output converter i ¯ 2 z , see Figure 7.
In this way, the discrete-time system to be controlled has the following linear transfer function:
G v i z = v ¯ o u t z i ¯ 2 z = R p z β z α   ,
where,
R p = R R C R + R C   ,   α = e T s C R + R C   ,   β = R + R C α R R C   .
Equations (11) and (12) have been obtained by the continuous-time counterpart G v i s = v ¯ o u t s / i ¯ 2 s and the zero-order hold (ZOH) discretization method. In particular, the transfer function
G v i s = R p s + β ¯ s + α ¯  
where,
α ¯ = 1 C R + R C ,   β ¯ = 1 C R C ,
has been derived by the circuit Kirchhoff’s laws of the circuit, shown in Figure 4, including the equivalent series resistance R C .
The discrete compensator C i z can be designed to calculate the control variable i ¯ 2 z to satisfy the given static and dynamic specifications.
The system is characterized by G v i z and is composed of the subsystems shown in Figure 8. In particular, the nonlinear function block f i ¯ 2 digitally computes the phase shift δ k as function of i ¯ 2 k using the following relation:
δ k = f i ¯ 2 = π 2 1 1 4 i ¯ 2 k ω n L π v ¯ b a t ,
which directly derives from Equation (5). The exact correspondence between i ¯ 2 and δ in Equation (15) avoids the linear approximation introduced in Equation (7). The SPS modulator generates the appropriate switching signals Q i according to the operation modes described in Table 1. The DC Bus voltage v ¯ o u t is then sampled and held by the ADC converter.

PI Controller Model

The considered discrete-time PI controller can be described by the following expression:
C i z = K P 1 + 1 T i z + 1 z 1 ,
and the corresponding frequency response in rectangular form is given in Equation (17):
C i ω , T s = K P 1 j 1 T i tan ω T s 2 .
Equation (17) can also be rewritten in polar form as a function of magnitude M ω , T s and phase φ ω , T s ,
C i ω , T s = M e j φ = M cos φ + j sin φ ,
while the frequency response of the discrete system of Equation (13) is:
G v i ω , T s = R p 1 + α β α + β cos ω T s + j β α sin ω T s 1 + α 2 2 α cos ω T s .
The proposed discrete formulae to exactly meet design specifications on zero position error, given gain crossover frequency ω g , and phase margin Φ m of the loop gain transfer function L z = C i z G v i z , are:
K P = M g cos φ g ,
T i = 1 tan ω g T s 2 tan φ g ,
where,
M g = 1 R p   1 + α 2 2 α cos ω g T s 1 + β 2 2 β cos ω g T s ,
φ g = Φ m π tan 1 β α sin ω g T s 1 + α β α + β cos ω g T s .
Since the DAB converter model is a type zero system, the pole at 1 of the PI controller guarantees zero steady-state error on nominal output voltage step response. Equations (20)–(23) have been obtained imposing that the open-loop response L ω , T s as unity magnitude and phase are equal to Φ m π at frequency ω = ω g . This is when C i ω g , T s G v i ω g , T s = 1 and C i ω g , T s + G v i ω g , T s = Φ m π . From Equation (19), magnitude M ω g = C i ω g , T s and the phase φ ω g = C i ω g , T s that the controller should have at a given frequency ω g , are:
M g = M ω g = 1 / G v i e j ω g T s
φ g = φ ω g = Φ m π G e j ω g T s .
By substituting Equations (17)–(19) in Equations (24) and (25), the following Equations (20)–(23) can be obtained, and this solution is admissible only if K P > 0 and T i > 0 [34].

4. Simulation Results

The PI controller and the SPS-DAB have been designed and carried out in the PLECS simulation software on two different resistive load values to verify the effectiveness of the proposed control in stabilizing the DC bus voltage. The simulation parameters are tabulated in Table 2.

4.1. PI Design

The discrete-time system Equation (11), with the parameters shown in Table 2. The circuit parameters of the DAB converter are:
G v i z = v ¯ o u t z i ¯ 2 z = 0.001 · z + 0.28 z 0.99   .
The PI parameters can be designed to yield zero position error, a phase margin Φ m and a gain crossover frequency ω g of the open-loop frequency response equal to 75° and 1200 rad/s, respectively. They are K P = 0.41 ,   T i = 60.58   s . Using inversion Formulae (20)–(23), the PI controller has the following expression:
C i z = 0.41   z 0.39 z 1 .
The Nyquist plot of the open-loop frequency response L ω , T s is shown in Figure 9.
Notice that the design specifications are precisely met. The green lines represent the Nyquist plot of the open-loop frequency response L ω , T s = C i ω , T s G v i ω , T s under a power load variation from 4 to 20 kW.

4.2. PLECS Simulation Results

The DAB converter and the controller are designed in two separate subsystems that can be used for offline and real-time simulation purposes. In particular, the controller subsystem has been developed, as shown in Figure 10. Nevertheless, the PI parameters have been derived for the discrete PI controller from Equation (16); they can also be used with other discrete-time controller types. Since PLECS discrete controller block implements the following parallel PI form selecting the trapezoidal integration method
C z = K P + K i T s 2   z   + 1 z 1   ,
the integral constant K i can be derived as
K i = K p T i 2 T s   ,
To evaluate the performance of the implemented control, the system is brought to a steady-state condition by connecting a 60 Ω resistor as a load. This resistor value represents an equivalent 6 kW load. At instant t = 10 ms, a 4 kW increase in the load is emulated by lowering the resistor value to 36 Ω. As shown in Figure 11, the controller can react against the load variation and restore the rated output voltage in roughly 11 ms. At any moment of the transient, the voltage remains well confined within the ±5% of the nominal voltage (represented with the gray band).
Figure 12 shows the output power and voltage changes when load abruptly decreases from 10 kW to 6 kW at t = 10 ms. During this load transient, the DC bus voltage reaches a peak value of 614 V, and again the control allows the new steady-state within t = 11 ms. Moreso, in this case, voltage evolution is well limited within the ±5% boundary.

4.3. Simulated Performance Comparison for Fast Inertial Action

To evaluate the performance of the proposed controller, a comparison with the results presented in [35] has been made. This work, in fact, proposes a different control technique applied to the same DAB architecture supporting the fast synthetic inertial action. This work is taken as a reference for the comparison because, to the best of the authors’ knowledge, the technique presented guarantees one of the fastest responses found in the literature. The converter operates on a DC line at the rated voltage of 500 V, and the results summary of [35] are as follows: when the DC load changes from 12 kW to 6 kW, the DC bus voltage reaches 520 V, and when the DC load changes from 6 kW to 12 kW, the DC bus voltage reduces to 490 V with the settling time 0.4 s to 0.5 s. The authors of the paper claim that the dynamic of the control can be accelerated by modifying the virtual inertial time constant and damping coefficient, but with the cost of ripples in output power and output voltage.
Considering the same working conditions of the converter studied in [35], the controller proposed in the present work demonstrated better results in recovering the rated DC bus voltage in a much shorter time and guaranteed no oscillations. The following results are achieved with the phase margin Φ m and the gain crossover frequency ω g of the open-loop frequency response equal to 75° and 700 rad/s, respectively. When the DC load changes from 12 kW to 6 kW, the bus voltage varies by only 3.3 V, and the nominal value is restored in 14 ms. In such a case, the K P   value remains fixed at 0.8006, while the T i value changes from 441.6 s to 488.9 s. When the DC load is increased from 6 kW to 12 kW, the DC bus voltage reduces to 496.9 V and the same settling time of 14 ms is maintained. In all the cases, the voltage does not overpass the ±5% nominal voltage threshold. The results are shown in Figure 13 and Figure 14, respectively.

4.4. Comparison of Proposed PI Controller Model with a Classical Pole-Placement Method

Apart from the previous discussion, to evaluate the effectiveness and to understand the robustness of the proposed discrete model, a comparison with the classical pole-placement method is studied in this section. Considering the same Table 2 circuit design specifications, new PI parameters are calculated using the classical pole-placement method for the second-order system [28,37]. In particular, considering the converter model (9), the parameters of the continuous-time PI controller C s = K P + K i S can be tuned by selecting the damping ratio and the natural frequency of the desired second-order closed-loop polynomial equal to 0.89 and 676 rad/s, respectively. These values lead to the same frequency design specifications considered in the proposed method when the load is at its nominal value. The resulting PI parameters are K P = 0.0051 and K i = 2.06. Bearing in mind that a discrete PI regulator has been used in the proposed model, the continuous-time PI controller in a classical pole-placement method is approximated using the trapezoidal discretization method and implemented in the PLECS simulation. Notice that the discretization process leads to a phase margin degradation from 75° to 73°. However, the PI parameters of a classical method are fixed, therefore the proposed model showed better results during the load transients. Figure 15 and Figure 16 show the simulation results during the load change from 6 to 10 kW and from 10 to 6 kW, respectively. Although in both the methods, the maximum overshoot and undershoot values are comparable and succeed in preserving the nominal voltage within the ±5% band, the settling time of the proposed controller is approximately 35% faster than the classical pole-placement method, maintaining the transient duration well within a grid fundamental period.

5. Hardware-in-the-Loop Results

The simulation results have been validated using the hardware-in-the-loop real-time testing platform RT Box 1 powered by PLECS. This device is fully integrated with the Xilinx Zynq system on a chip consisting of 64 digital and 32 analog I/O ports. The control hardware, i.e., Texas Instrument LaunchPad XL running the TMS320F28069M DSP, is connected to the RT Box for real-time operations using an RT Box LaunchPad Interface (Figure 17) [38]. In this experiment, the RT box is used to emulate the DAB converter, and the DSP runs the control scheme. The PWMs are generated with the switching frequency f s   equal to 20 kHz, whereas the DAB converter operations are emulated with a sampling frequency of 100 kHz. The increase in the number of samples allows calculating the DAB converter’s output voltage accurately over a switching period. Due to the hardware limitations of the DSP, the operating frequency of the controller is limited to 10 kHz. The external setup has been made by connecting the RT Box 1 using an ethernet cable, and the DSP using the USB cable to the human interface device (Laptop/PC), which has PLECS (version 4.6.1) software installed in it. Once this external connection is made, the DAB and controller models are converted into the C code using the PLECS embedded coder. This PLECS embedded coder is a setting available in the PLECS software that enables the connection between the RT Box 1, DSP, and the PC and generates the C code for both the RT Box1 and DSP. More details regarding the setup and the code implementation procedure can be found in [39]. PLECS library has an inbuilt DAB converter block, which is suitable to use in the RT Box 1. Along with this converter block, analog I/O, digital I/O, analog-to-digital converter (ADC) blocks were used for transmitting and receiving the signals in between the RT Box 1 and DSP. Table 2 and Table 3 parameters are considered and executed in the HIL experiment. The output power and output voltage are measured during the real-time operations, similar to the offline simulation.
The results of the real-time experiment shown in Figure 18 and Figure 19 are in good accordance with the results obtained from the offline simulation shown in Figure 11 and Figure 12.
The HIL results show that the DC bus voltage reaches the steady-state within 11 ms and the bus voltage varies by 12 V (i.e., 2%) during the DC load transients. A sub-cycle synthetic inertial action can be effectively supported by the settling time lower than a fundamental period and the voltage transient within the classical ±5% band.
For better understanding and visualizing purposes, Table 4 provides the comparison of obtained simulation results with HIL results. In both cases, the settling time is 11   ms, whereas the voltage maximum overshoot and undershoot values are given below.

6. Conclusions

This paper addresses the discrete DAB converter model and the controller design for fast recovering voltage/power transients on the DC bus with sub-cycle transients compatible with synthetic inertia action. The proposed PI controller satisfies the exact standard design constraints, such as gain crossover frequency, phase margin, and zero position error. Moreover, the design has been developed so that PI parameters are updated in real-time without additional dedicated sensors. To prove the performance of the proposed controller, a comparison with the DAB parameters considered in [35] in the virtual inertial application point of view has been made and achieved the settling time of 14 ms with a 3.3 V voltage variation. Moreover, to understand the controller robustness, another comparison with the classical pole-placement method has been made. Even in this case, the results were promising, as it is proved that the proposed controller is 35% faster than the classical pole-placement method. HIL testing has also been carried out to verify that the controller can be implemented in real systems. Practically, the power/voltage transients compromise the safe and stable operation of DC microgrids or IDS networks, unless they are recovered within tens of milliseconds. Therefore, the proposed controller can be utilized in order to recover those transients and ultimately maintain grid stability. Since the controller design is based on the DAB converter power flow, more complex modulation strategies like dual-, triple-, and extended-phase-shift as well as the DAB in its three-phase structures, might be studied in the future to try and overcome the difficulties mentioned in the introduction, such as reducing the complex calculations and achieve unified control. The possibility to introduce supercapacitors and hybrid battery packs will also be considered in further contributions.

Author Contributions

Conceptualization, S.C., R.M., L.K.P., V.C. and M.R.; methodology, S.C., R.M. and L.K.P.; software, S.C., L.K.P. and V.C.; validation, S.C., R.M. and L.K.P.; formal analysis, S.C., R.M., V.C. and M.R.; investigation, S.C. and R.M.; resources, S.C., R.M. and L.K.P.; data curation, S.C., R.M., L.K.P., V.C. and M.R.; writing—original draft preparation, S.C., R.M., L.K.P. and V.C.; writing—review and editing, S.C., R.M., L.K.P., V.C. and M.R.; visualization, S.C., R.M. and L.K.P.; supervision, R.M., V.C. and M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is available in the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

v b a t , v o u t Converter’s input voltage (battery voltage) and output voltage (DC bus voltage)
i 1 , i 2 Converter’s input current and output current
I 2 ,  i ^ 2 Average current, and corresponding small-signal variable
v C ,   i C ,   R C Converter’s output capacitor voltage, current, equivalent series resistor
v p ,   v s Primary and secondary voltages of a high-frequency transformer
Q i Converter’s switches with index i (from 1 to 8)
R Equivalent resistive load
δ ,  φ ,   δ ^ t   Phase-shift between two bridges, nominal/constant phase-shift, and corresponding small-signal variable
P 1 ,   P 2 ,   P Power of the primary full-bridge, secondary full-bridge and total output rated power of the converter
G i 2 d s , G v d s   Transfer function relating the phase angle to the output current and output voltage in s-domain
G i 2 d z , G v d z Transfer function relating the phase angle to the output current and output voltage in z-domain
G v i ,   G v i ω , T s Transfer function relating the output voltage and output current in both s, z-domain, and corresponding frequency response
C i z , C i ω , T s Discrete PI controller relating to output current in z-domain and corresponding frequency response
C d z Discrete PI controller relating to phase angle in z-domain
ω g ,   Φ m Gain crossover frequency and phase margin
L z ,   L ω , T s   An open-loop gain transfer function in z-domain and corresponding frequency response

References

  1. Long, B.; Liao, Y.; Chong, K.T.; Rodríguez, J.; Guerrero, J.M. MPC-Controlled Virtual Synchronous Generator to Enhance Frequency and Voltage Dynamic Performance in Islanded Microgrids. IEEE Trans. Smart Grid 2021, 12, 953–964. [Google Scholar] [CrossRef]
  2. Ramana, N.V. Power System Analysis; Pearson Education: London, UK, 2011; ISBN 9788131755921. [Google Scholar]
  3. Tamrakar, U.; Shrestha, D.; Maharjan, M.; Bhattarai, B.P.; Hansen, T.M.; Tonkoski, R. Virtual inertia: Current trends and future directions. Appl. Sci. 2017, 7, 654. [Google Scholar] [CrossRef]
  4. Algarni, A.S.; Member, G.S.; Suryanarayanan, S.; Member, S.; Siegel, H.J.; Fellow, L.; Maciejewski, A.A. Combined Impact of Demand Response Aggregators and Carbon Taxation on Emissions Reduction in Electric Power Systems. IEEE Trans. Smart Grid 2021, 12, 1825–1827. [Google Scholar] [CrossRef]
  5. Obi, M.; Bass, R. Trends and challenges of grid-connected photovoltaic systems—A review. Renew. Sustain. Energy Rev. 2016, 58, 1082–1094. [Google Scholar] [CrossRef]
  6. Raza, S.A.; Jiang, J.I.N. A Benchmark Distribution System for Investigation of Residential Microgrids. IEEE Open-Access J. Power Energy 2020, 7, 41–50. [Google Scholar] [CrossRef]
  7. Gomez, J.S.; Rodriguez, J.; Garcia, C.; Tarisciotti, L.; Flores-Bahamonde, F.; Pereda, J.; Nunez, F.; Cipriano, A.Z.; Salas, J.C. An Overview of Microgrids Challenges in the Mining Industry. IEEE Access 2020, 8, 191378–191393. [Google Scholar] [CrossRef]
  8. Makolo, P.; Zamora, R.; Lie, T.T. The role of inertia for grid flexibility under high penetration of variable renewables—A review of challenges and solutions. Renew. Sustain. Energy Rev. 2021, 147, 111223. [Google Scholar] [CrossRef]
  9. Abubakr, H.; Mohamed, T.H.; Hussein, M.M.; Guerrero, J.M.; Agundis-Tinajero, G. Adaptive frequency regulation strategy in multi-area microgrids including renewable energy and electric vehicles supported by virtual inertia. Int. J. Electr. Power Energy Syst. 2021, 129, 106814. [Google Scholar] [CrossRef]
  10. Adu, J.A.; Rios Penaloza, J.D.; Napolitano, F.; Tossani, F. Virtual Inertia in a Microgrid with Renewable Generation and a Battery Energy Storage System in Islanding Transition. In Proceedings of the 2019 1st International Conference on Energy Transition in the Mediterranean Area (SyNERGY MED), Cagliari, Italy, 28–30 May 2019; pp. 1–5. [Google Scholar] [CrossRef]
  11. Mallemaci, V.; Mandrile, F.; Rubino, S.; Mazza, A.; Carpaneto, E.A. comprehensive comparison of Virtual Synchronous Generators with focus on virtual inertia and frequency regulation. Electr. Power Syst. Res. 2021, 201, 107516. [Google Scholar] [CrossRef]
  12. Sarojini, R.K.; Palanisamy, K. Emulated inertia control for the stand-alone microgrid with high penetration of renewable energy sources. Int. J. Renew. Energy Res. 2020, 10, 831–842. [Google Scholar] [CrossRef]
  13. Sarojini, R.K.; Kaliannan, P.; Teekaraman, Y.; Nikolovski, S.; Baghaee, H.R. An Enhanced Emulated Inertia Control for Grid-Connected PV Systems with HESS in a Weak Grid. Energies 2021, 14, 1721. [Google Scholar] [CrossRef]
  14. Bagheri, S.; CheshmehBeigi, H.M. DC Microgrid Voltage Stability through Inertia Enhancement Using a Bidirectional DC-DC Converter. In Proceedings of the 7th Iran Wind Energy Conference (IWEC2021), Shahrood, Iran, 17–18 May 2021; pp. 1–5. [Google Scholar] [CrossRef]
  15. Zhao, B.; Song, Q.; Liu, W.; Sun, Y. Overview of dual-active-bridge isolated bidirectional DC-DC converter for high-frequency-link power-conversion system. IEEE Trans. Power Electron. 2014, 29, 4091–4106. [Google Scholar] [CrossRef]
  16. Chen, X.; Xu, G.; Han, H.; Liu, D.; Sun, Y.; Su, M. Light-Load Efficiency Enhancement of High-frequency Dual-Active-Bridge Converter under SPS Control. IEEE Trans. Ind. Electron. 2021, 68, 12941–12946. [Google Scholar] [CrossRef]
  17. Zhao, B.; Yu, Q.; Sun, W. Extended-phase-shift control of isolated bidirectional DC-DC converter for power distribution in microgrid. IEEE Trans. Power Electron. 2012, 27, 4667–4680. [Google Scholar] [CrossRef]
  18. Inoue, S.; Akagi, H. A Bidirectional DC-DC Converter for an Energy Storage System With Galvanic Isolation. IEEE Trans. Power Electron. 2007, 22, 2299–2306. [Google Scholar] [CrossRef]
  19. Choi, K.; Kim, Y.; Kim, K.; Kim, S. Output voltage tracking controller embedding auto-tuning algorithm for DC/DC boost converters. IET Power Electron. 2019, 12, 3767–3773. [Google Scholar] [CrossRef]
  20. Pal, A.; Kapat, S.; Jha, K.; Tiwari, A. Discrete-time framework for digital control design in a high-frequency dual active bridge converter. In Proceedings of the 2018 IEEE Applied Power Electronics Conference and Exposition (APEC), San Antonio, TX, USA, 4–8 March 2018; pp. 2264–2270. [Google Scholar] [CrossRef]
  21. Middlebrook, R.D.; Cuk, S. A general unified approach to modelling switching-converter power stages. In Proceedings of the 1976 IEEE Power Electronics Specialists Conference, Cleveland, OH, USA, 8–10 June 1976; pp. 18–34. [Google Scholar] [CrossRef]
  22. De Doncker, R.W.; Divan, D.M.; Kheraluwala, M.H. A three-phase soft-switched high power density DC/DC converter for high power applications. In Proceedings of the Conference Record of the 1988 IEEE Industry Applications Society Annual Meeting, Pittsburgh, PA, USA, 2–7 October 1988; Volume 1, pp. 796–805. [Google Scholar] [CrossRef]
  23. Sanders, S.R.; Noworolski, J.M.; Liu, X.Z.; Verghese, G.C. Generalized averaging method for power conversion circuits. In Proceedings of the 21st Annual IEEE Conference on Power Electronics Specialists, San Antonio, TX, USA; 1990; pp. 333–340. [Google Scholar] [CrossRef] [Green Version]
  24. Qin, H.; Kimball, J.W. Generalized average modeling of dual active bridge DC-DC converter. IEEE Trans. Power Electron. 2012, 27, 2078–2084. [Google Scholar] [CrossRef]
  25. Liu, B.; Davari, P.; Blaabjerg, F. An Enhanced Generalized Average Modeling of Dual Active Bridge Converters. In Proceedings of the 2020 IEEE Applied Power Electronics Conference and Exposition (APEC), New Orleans, LA, USA, 15–19 March 2020; pp. 85–90. [Google Scholar] [CrossRef]
  26. Iqbal, M.T.; Member, S.; Maswood, A.I.; Member, S. An Explicit Discrete-Time Large- and Small-Signal Modeling of the Dual Active Bridge DC—DC Converter Based on the Time Scale Methodology. IEEE J. Emerg. Sel. Top. Ind. Electron. 2021, 2, 545–555. [Google Scholar] [CrossRef]
  27. Costa, P.F.S.; Löbler, P.H.B.; Roggia, L.; Schuch, L. Modeling and Control of DAB Converter Applied to Batteries Charging. IEEE Trans. Energy Convers. 2022, 37, 175–184. [Google Scholar] [CrossRef]
  28. Santos da Silva, E.L.; Luís Kirsten, A.; Pagano, D.J. Discrete SPS Control of a DAB converter using partial Feedback Linearization. In Proceedings of the 2019 IEEE 15th Brazilian Power Electronics Conference and 5th IEEE Southern Power Electronics Conference (COBEP/SPEC), Santos, Brazil, 1–4 December 2019; pp. 1–6. [Google Scholar] [CrossRef]
  29. Tong, A.; Hang, L.; Li, G.; Huang, J. Nonlinear characteristics of DAB converter and linearized control method. In Proceedings of the 2018 IEEE Applied Power Electronics Conference and Exposition (APEC), San Antonio, TX, USA, 4–8 March 2018; pp. 331–337. [Google Scholar] [CrossRef]
  30. Tiwari, S.; Hanif, O.; Sarangi, S. Fractional Order PI Control of Dual Active Bridge Converter Using Generalized Average Modelling. In Proceedings of the 2019 20th International Conference on Intelligent System Application to Power Systems (ISAP), New Delhi, India, 10–14 December 2019; pp. 253–259. [Google Scholar] [CrossRef]
  31. Xue, D.; Zhao, C.; Chen, Y. Fractional order PID control of a DC-motor with elastic shaft: A case study. In Proceedings of the 2006 American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; p. 6. [Google Scholar] [CrossRef]
  32. Yang, X. Nature-Inspired Optimization Algorithms; Elsevier: Amsterdam, The Netherlands, 2020; ISBN 0124167454. [Google Scholar]
  33. Korlinchak, C.; Comanescu, M. Discrete time integration of observers with continuous feedback based on Tustin’s method with variable prewarping. In Proceedings of the 6th IET International Conference on Power Electronics, Machines and Drives (PEMD 2012), Bristol, UK, 27–29 March 2012; pp. 1–6. [Google Scholar] [CrossRef]
  34. Cuoghi, S.; Mandrioli, R.; Ntogramatzidis, L.; Gabriele, G. Multileg Interleaved Buck Converter for EV Charging: Discrete-Time Model and Direct Control Design. Energies 2020, 13, 466. [Google Scholar] [CrossRef] [Green Version]
  35. Guo, Y.; Ma, W.; Meng, J.; Wang, Y. A Virtual Inertia Control Strategy for Dual Active Bridge DC-DC Converter. In Proceedings of the 2018 2nd IEEE Conference on Energy Internet and Energy System Integration (EI2), Beijing, China, 20–22 October 2018; pp. 1–5. [Google Scholar] [CrossRef]
  36. Kheraluwala, M.H.; De Doncker, R.W. Single phase unity power factor control for dual active bridge converter. In Proceedings of the Conference Record of the 1993 IEEE Industry Applications Conference Twenty-Eighth IAS Annual Meeting, Toronto, ON, Canada, 2–8 October 1993; Volume 2, pp. 909–916. [Google Scholar] [CrossRef]
  37. Ogata, K. Modern Control. Engineering; Prentice Hall: Hoboken, NJ, USA, 2010; ISBN 9780136156734. [Google Scholar]
  38. RT Box LaunchPad Interface. May 2020. Available online: https://www.plexim.com/sites/default/files/launchpadinterfacemanual.pdf (accessed on 15 March 2022).
  39. RT Box User Manual. January 2021. Available online: https://www.plexim.com/sites/default/files/rtboxmanual.pdf (accessed on 15 March 2022).
Figure 1. General scheme of an inverter-dominated microgrid.
Figure 1. General scheme of an inverter-dominated microgrid.
Energies 15 02295 g001
Figure 2. DAB converter circuit for battery energy storage system (BESS) connected to DC Bus.
Figure 2. DAB converter circuit for battery energy storage system (BESS) connected to DC Bus.
Energies 15 02295 g002
Figure 3. General SPS-DAB switching sequence and related main quantities.
Figure 3. General SPS-DAB switching sequence and related main quantities.
Energies 15 02295 g003
Figure 4. SPS-DAB model developed by averaging the converter output current.
Figure 4. SPS-DAB model developed by averaging the converter output current.
Energies 15 02295 g004
Figure 5. Graphical presentation of Equation (5) (blue) and its linearization around the operation point of Equation (7) (orange).
Figure 5. Graphical presentation of Equation (5) (blue) and its linearization around the operation point of Equation (7) (orange).
Energies 15 02295 g005
Figure 6. Control block scheme using OCL model method.
Figure 6. Control block scheme using OCL model method.
Energies 15 02295 g006
Figure 7. The proposed voltage control scheme.
Figure 7. The proposed voltage control scheme.
Energies 15 02295 g007
Figure 8. Detail view of the proposed voltage control scheme.
Figure 8. Detail view of the proposed voltage control scheme.
Energies 15 02295 g008
Figure 9. Nyquist plot of the open-loop frequency response L ω , T s = C i ω , T s G v i ω , T s with the nominal power load value (red) and under load variation from 4 to 20 kW (green).
Figure 9. Nyquist plot of the open-loop frequency response L ω , T s = C i ω , T s G v i ω , T s with the nominal power load value (red) and under load variation from 4 to 20 kW (green).
Energies 15 02295 g009
Figure 10. PLECS controller subsystem.
Figure 10. PLECS controller subsystem.
Energies 15 02295 g010
Figure 11. Output power (a) and output voltage (b) change on the DC bus when load changes from 6 kW to 10 kW with the proposed control method.
Figure 11. Output power (a) and output voltage (b) change on the DC bus when load changes from 6 kW to 10 kW with the proposed control method.
Energies 15 02295 g011
Figure 12. Output power (a) and output voltage (b) change on the DC bus when load changes from 10 kW to 6 kW with the proposed control method.
Figure 12. Output power (a) and output voltage (b) change on the DC bus when load changes from 10 kW to 6 kW with the proposed control method.
Energies 15 02295 g012
Figure 13. Output power (a) and output voltage (b) change on the DC bus when load changes from 12 kW to 6 kW in comparison with the same parameters of the comparative article [35].
Figure 13. Output power (a) and output voltage (b) change on the DC bus when load changes from 12 kW to 6 kW in comparison with the same parameters of the comparative article [35].
Energies 15 02295 g013
Figure 14. Output power (a) and output voltage (b) change on the DC bus when load changes from 6 kW to 12 kW in comparison with the same parameters of the comparative article [35].
Figure 14. Output power (a) and output voltage (b) change on the DC bus when load changes from 6 kW to 12 kW in comparison with the same parameters of the comparative article [35].
Energies 15 02295 g014
Figure 15. Output power (a) and output voltage (b) change on the DC bus when load changes from 6 kW to 10 kW of both classical pole-placement and proposed control methods.
Figure 15. Output power (a) and output voltage (b) change on the DC bus when load changes from 6 kW to 10 kW of both classical pole-placement and proposed control methods.
Energies 15 02295 g015
Figure 16. Output power (a) and output voltage (b) change on the DC bus when load changes from 10 kW to 6 kW of both classical pole-placement and proposed control methods.
Figure 16. Output power (a) and output voltage (b) change on the DC bus when load changes from 10 kW to 6 kW of both classical pole-placement and proposed control methods.
Energies 15 02295 g016
Figure 17. Hardware-in-the-loop testing setup.
Figure 17. Hardware-in-the-loop testing setup.
Energies 15 02295 g017
Figure 18. Real-Time HIL results (a) output power, (b) output voltage during load variation from 6 kW to 10 kW with the proposed control method.
Figure 18. Real-Time HIL results (a) output power, (b) output voltage during load variation from 6 kW to 10 kW with the proposed control method.
Energies 15 02295 g018
Figure 19. Real-Time HIL results (a) output power, (b) output voltage during load variation from 10 kW to 6 kW with the proposed control method.
Figure 19. Real-Time HIL results (a) output power, (b) output voltage during load variation from 10 kW to 6 kW with the proposed control method.
Energies 15 02295 g019
Table 1. SPS-DAB Operation Modes.
Table 1. SPS-DAB Operation Modes.
ModeConducting SwitchesDuration
1 Q 1   Q 4   Q 6   Q 7 δ π   T s 2
2 Q 1   Q 4   Q 5   Q 8 1 δ π   T s 2
3 Q 2   Q 3   Q 5   Q 8 δ π   T s 2
4 Q 2   Q 3   Q 6   Q 7 1 δ π   T s 2
Table 2. Circuit parameters of the dual-active-bridge converter.
Table 2. Circuit parameters of the dual-active-bridge converter.
ParameterSymbolValueUnits
Battery voltage v b a t 600 V
Rated DC bus voltage v o u t 600 V
Filter capacitance C 350 μ F
Series resistance R c 1 m Ω
Leakage inductance L 53.64 μ H
Rated power P 10 kW
Switching frequency f s 20 kHz
Sampling time T s 0.1 ms
Transformer ratio n 1 -
Table 3. PLECS hardware-in-the-loop parameters.
Table 3. PLECS hardware-in-the-loop parameters.
ParameterValueUnits
Solver step size 1 ms
Discretization step size for converter 10 μ s
Discretization step size for controller 100 μ s
Analog output scaling factor for voltage 3.3 / 800 V
Analog output scaling factor for current 3.3 / 40 A
Analog output offset value (voltage and current) 0 -
ADC voltage scaling factor 800 / 3.3 V
ADC current scaling factor 40 / 3.3 A
ADC offset value 0 -
Table 4. Comparison of offline simulation and hardware-in-the-loop results.
Table 4. Comparison of offline simulation and hardware-in-the-loop results.
ParameterNominal Values6 kW to 10 kW
Load Change
10 kW to 6 kW
Load Change
Simulation ResultHIL ResultSimulation ResultHIL Result
Power (kW) 10 9.58 maximum undershoot value 9.60 maximum undershoot value 6.30 maximum overshoot value 6.25 maximum overshoot value
DC bus voltage (V) 600 588 maximum undershoot value 588 maximum undershoot value 614 maximum overshoot value 614 maximum overshoot value
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cuoghi, S.; Mandrioli, R.; Pittala, L.K.; Cirimele, V.; Ricco, M. Dual-Active-Bridge Model and Control for Supporting Fast Synthetic Inertial Action. Energies 2022, 15, 2295. https://doi.org/10.3390/en15062295

AMA Style

Cuoghi S, Mandrioli R, Pittala LK, Cirimele V, Ricco M. Dual-Active-Bridge Model and Control for Supporting Fast Synthetic Inertial Action. Energies. 2022; 15(6):2295. https://doi.org/10.3390/en15062295

Chicago/Turabian Style

Cuoghi, Stefania, Riccardo Mandrioli, Lohith Kumar Pittala, Vincenzo Cirimele, and Mattia Ricco. 2022. "Dual-Active-Bridge Model and Control for Supporting Fast Synthetic Inertial Action" Energies 15, no. 6: 2295. https://doi.org/10.3390/en15062295

APA Style

Cuoghi, S., Mandrioli, R., Pittala, L. K., Cirimele, V., & Ricco, M. (2022). Dual-Active-Bridge Model and Control for Supporting Fast Synthetic Inertial Action. Energies, 15(6), 2295. https://doi.org/10.3390/en15062295

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop