High Frequency Transformers for Solid-State Transformer Applications

: This paper focuses on the study of the high frequency transformer incorporated in solid-state transformers, speciﬁcally on the development of the steps that enable the design of an optimized high frequency transformer and its equivalent model based on the desired characteristics. The impact of operating a transformer at high frequency and the respective solutions that allow this impact to be reduced are analyzed, alongside the numerous advantages that the utilization of these transformers has over traditional 50/60 Hz transformers. Furthermore, the power scheme of the solid-state transformer is outlined, focusing on the power converters, which are immediately before and after the high frequency transformer (HFT). We also investigate a control technique that allows for correct operation and the existence of power bidirectionality. In a novel approach, this paper demonstrates the systematic steps for designing an HFT according to the desired speciﬁcations of each given project, helping students and engineers achieve their objectives in power-electronic applications. Moreover, this paper aims at increasing the knowledge of this area of power electronics and facilitating the development of new topologies with high power density, which are very important to the integration of renewable power sources and other applications. Finally, a simulation is presented to validate a high frequency transformer and its control technique.


Introduction
The solid-state transformer (SST) is an emerging technology that can influence power transmission/distribution networks, their integration with renewable energy sources, and traction systems, among others.In recent years, the integration of renewable energy sources has increased significantly, leading to numerous investigations to integrate this principle with electrical distribution networks and their optimization [1,2].Some advantages of using SSTs instead of conventional transformers, which operate at a low frequency of 50/60 Hz, include reactive power compensation, voltage regulation, power flow control, and the possibility of easy interconnection with equipment, allowing the improvement of power quality.One of the main advantages of using SSTs instead of traditional transformers is the reduction of volume and weight.For example, according to [1], a three-phase STT would accomplish a volume reduction of 80% compared to a conventional 50/60 Hz transformer.This significant difference has a great impact as it facilitates portability, decreases installation costs, and significantly reduces the consumption of raw materials [1].
Some relevant works about the design of HFTs based on a two-level voltage converter include [3], which presents and discusses the effect and influence of HFT parasitic capacitances on the operation and design of an SST with two-level converts.The optimization issues of HFT are addressed in [4], via a physics-based design optimization algorithm combined with multi-objective genetic algorithms.In [5], HFT design in relation to Multi Port SST (MPSST) applications is discussed.Different winding placement and distribution on the same core are suggested, and the inductance and coupling coefficient matrices are calculated using ANSYS Maxwell 3D simulation software.In [6], an optimum design of HFT based on winding spacing, minimum weight, and efficiency constraints is presented.This transformer design takes into account series resonant converter applications.Several other papers with HFT for resonant converters can be found in the literature.In [7], the design of HFT for LLC resonant converters with SiC-based semiconductors and high efficiency is proposed.A similar design solution using a GaN power semiconductor is also proposed in [8].An extended comparison among different HFTs for resonant power converters can be found in [9].Several papers about multilevel inverters for SSTs have also been published.In [10], an HFT with a high-voltage insulation structure for SST for medium-voltage multilevel converters is presented.A comparative analysis of a multilevel high frequency link and multilevel DClink DC-DC HFT based on Modular Multilevel Converters (MMC) and a Dual-Active Bridge (DAB) for MVDC application can be found in [11].Several other studies about HFT and SST discuss different types of transformers.For example, [12], proposes a comparative study of two types of HFT (one made of planar cores and windings in the form of thin copper foil and the other built with E-type cores and Litz wire) in a 20-kW phase-shifted full-bridge converter application.An analysis of the electromagnetic characteristics of HFT based on C-type nanocrystalline magnetic cores is proposed in [13].Another study presents the design and the main advantages of dry-type high-power HFTs based on triangular closed cores [14].Considerations about the design, modeling, and analysis of compact planar transformers can be found in [15].The high-frequency modeling of foil-type transformers for shipboard power-electronic distribution systems is studied and analyzed in [16].
However, to achieve the desired operational objectives of SST, many challenges must still be overcome, especially in the dimensions of the components of the High Frequency Transformer.These include determining the optimal flux density and minimizing power losses to maximize efficiency.To overcome these challenges, it is necessary to know the impact of high frequency operation on the transformer and the options that exist to minimize this impact, such as utilizing different topologies and materials for the magnetic material structure, windings, and conductors [17].
Another challenge to overcome is the constitution of the power converters, namely the power semiconductors that allow the operation of the SST at high power and frequency.Moreover, the development of a control solution allows for maximizing the system efficiency, but this last challenge is not addressed in this article.
The main research objective of this paper is to study the operation of an HFT by considering different types of magnetic cores to achieve higher power densities, and furthermore to investigate different equivalent circuits of the transformer, such as the impact of stray parasitic capacitances.The novelty of this paper lies in demonstrating the systematic steps to design an HFT according to the desired specifications of a given project, helping students and engineers to achieve these objectives in power electronic applications.Moreover, we aim at widening the knowledge of this area of power electronics and improving the development of new topologies with high power density, an area important to the integration of renewable power sources and other applications.
The paper is organized as follows.Section 2 presents the modulation of the converter system, showing the application of SST, its characteristics and design issues.Furthermore, we discuss modeling aspects of converters with SST transformers.Section 3 presents some control aspects of SST and proposes a control solution.Section 4 presents the results of simulations for a fast charger of an electric vehicle with power bidirectionality.Finally, Section 5 is dedicated to conclusions.Table 1 shows the variables and parameters used throughout the paper.The SST can assume several topologies; some examples are given in Figure 1.The main difference between the various topologies is the number of constituent stages, the availability of direct current (DC) buses, the allowed functions, and the control complexity [1,18].
The three-stage topology, Type D, has the greatest popularity due to its features such as reactive power compensation, voltage drop compensation, integration of renewable energy sources, and integration with energy storage systems, among others.It is designed to provide two DC buses, a low voltage (LV), and a high voltage (HV), which allow the integration of other components to improve the power quality or to supply direct current equipment.The main disadvantage of this topology is the large number of power semiconductors required, which generally translates into lower efficiency and lower reliability [18,19].
In all the topologies depicted in Figure 1, power bidirectionality can be ensured by using power semiconductors that support operation in all four quadrants.The topology of the transformer is as important as the components in meeting the requirements of power density and efficiency.As SST applications, we highlight the solenoidal structures, namely the core type, shell type, and matrix type in Figure 2 [2,18].The three-stage topology, Type D, has the greatest popularity due to its features such as reactive power compensation, voltage drop compensation, integration of renewable energy sources, and integration with energy storage systems, among others.It is designed to provide two DC buses, a low voltage (LV), and a high voltage (HV), which allow the integration of other components to improve the power quality or to supply direct current equipment.The main disadvantage of this topology is the large number of power semiconductors required, which generally translates into lower efficiency and lower reliability [18,19].
In all the topologies depicted in Figure 1, power bidirectionality can be ensured by using power semiconductors that support operation in all four quadrants.The topology of the transformer is as important as the components in meeting the requirements of power density and efficiency.As SST applications, we highlight the solenoidal structures, namely the core type, shell type, and matrix type in Figure 2 [2,18].The main advantage of the core type is that the windings are around a single magnetic core, which results in a weight reduction of the transformer [2].However, a large part of the magnetic core is surrounded by windings, which, associated with the poor thermal conductivity of the magnetic core, become a disadvantage of the core type and consequently an advantage of the shell type from a thermal point of view.Apparently, the matrix type has the same advantages as the shell type, but it has two major disadvantages that restrict its use: first, the number of magnetic cores, which results in the increased volume and weight of the transformer and, consequently, increased losses and cost; second, the primary and secondary are separated at opposite ends of the magnetic core, leading to higher leakage inductance, an unfavorable quality in large power converters [2].The three-stage topology, Type D, has the greatest popularity due to its features such as reactive power compensation, voltage drop compensation, integration of renewable energy sources, and integration with energy storage systems, among others.It is designed to provide two DC buses, a low voltage (LV), and a high voltage (HV), which allow the integration of other components to improve the power quality or to supply direct current equipment.The main disadvantage of this topology is the large number of power semiconductors required, which generally translates into lower efficiency and lower reliability [18,19].
In all the topologies depicted in Figure 1, power bidirectionality can be ensured by using power semiconductors that support operation in all four quadrants.The topology of the transformer is as important as the components in meeting the requirements of power density and efficiency.As SST applications, we highlight the solenoidal structures, namely the core type, shell type, and matrix type in Figure 2 [2,18].The main advantage of the core type is that the windings are around a single magnetic core, which results in a weight reduction of the transformer [2].However, a large part of the magnetic core is surrounded by windings, which, associated with the poor thermal conductivity of the magnetic core, become a disadvantage of the core type and consequently an advantage of the shell type from a thermal point of view.Apparently, the matrix type has the same advantages as the shell type, but it has two major disadvantages that restrict its use: first, the number of magnetic cores, which results in the increased volume and weight of the transformer and, consequently, increased losses and cost; second, the primary and secondary are separated at opposite ends of the magnetic core, leading to higher leakage inductance, an unfavorable quality in large power converters [2].The main advantage of the core type is that the windings are around a single magnetic core, which results in a weight reduction of the transformer [2].However, a large part of the magnetic core is surrounded by windings, which, associated with the poor thermal conductivity of the magnetic core, become a disadvantage of the core type and consequently an advantage of the shell type from a thermal point of view.Apparently, the matrix type has the same advantages as the shell type, but it has two major disadvantages that restrict its use: first, the number of magnetic cores, which results in the increased volume and weight of the transformer and, consequently, increased losses and cost; second, the primary and secondary are separated at opposite ends of the magnetic core, leading to higher leakage inductance, an unfavorable quality in large power converters [2].
The main challenge in the design of conductors for high frequency operation is the effect of eddy currents in the windings and their contributions to system losses, namely, the losses due to skin effect and proximity effect, outlined as follows:

•
Skin effect: the passage of an alternating current through a conductor causes a timevarying magnetic field, which, according to Faraday's law, causes an electric field.The electric field induces eddy currents in the conductor, which oppose the variation in current density.The combination of the magnetic field and the induced electric field, both time-varying, creates electromagnetic induction in the conductor.This forces most of the electric current to flow at the surface of the conductor, instead of flowing uniformly.The skin effect is aggravated with increased operating frequency [20,21]; • Proximity effect: when two conductors close to each other are crossed by an alternating current, the current distribution in each conductor is affected by the changing magnetic field of the other conductor.Due to the induced eddy currents, the current will flow in undesired patterns, and this effect is also aggravated with increasing frequency [20,21].
Ultimately, these effects can also be interpreted as increases in the resistivity of the conductor as the frequency increases.
Of the mechanisms developed so far, the most efficient for reduction is the Litz wire.This consists of a set of single conductors isolated from each other, which can also be braided together in a uniform pattern, as illustrated in Figure 3 [22].
time-varying, creates electromagnetic induction in the conductor.This forces most of the electric current to flow at the surface of the conductor, instead of flowing uniformly.The skin effect is aggravated with increased operating frequency [20,21]; • Proximity effect: when two conductors close to each other are crossed by an alternating current, the current distribution in each conductor is affected by the changing magnetic field of the other conductor.Due to the induced eddy currents, the current will flow in undesired patterns, and this effect is also aggravated with increasing frequency [20,21].
Ultimately, these effects can also be interpreted as increases in the resistivity of the conductor as the frequency increases.
Of the mechanisms developed so far, the most efficient for reduction is the Litz wire.This consists of a set of single conductors isolated from each other, which can also be braided together in a uniform pattern, as illustrated in Figure 3 [22].For a given application, there is not just one ideal conductor topology.The best solution will depend on the objectives, e.g., minimization of volume, minimization of losses, or minimization of costs.Determining the best conductor topology for a given application requires calculation methods, as in [20].

Transformer Modeling
Over the last decades, the trend has been to increase the power and frequency (typically from some kHz to dozens of MHz) of magnetic circuits.Nevertheless, this leads to a challenging design of magnetic circuits where parasitic effects play a key role.According to [23,24], most of the proposed designs are based on optimization routines, depending on the temperature rise and assuming uniform heat dissipation.Some authors consider the desired efficiency as a design input for the HFT [24]; otherwise, the optimization is focused on the Ap area [25,26].In [27], the magnetizing inductance is considered the first step for the design process of HFT.In [28], The leakage inductance is first investigated, and then the remaining parameters of HFT are given.
This paper proposes a step-by-step methodology for designing medium and highfrequency transformers and their optimum electrical and magnetic specifications to reduce volume and increase efficiency.The proposed design methodology for the HFT first finds the optimum magnetic-flux density and then arrives at the other parameters according to the desired criteria such as leakage inductance and losses, as shown in Figure 4.For a given application, there is not just one ideal conductor topology.The best solution will depend on the objectives, e.g., minimization of volume, minimization of losses, or minimization of costs.Determining the best conductor topology for a given application requires calculation methods, as in [20].

Transformer Modeling
Over the last decades, the trend has been to increase the power and frequency (typically from some kHz to dozens of MHz) of magnetic circuits.Nevertheless, this leads to a challenging design of magnetic circuits where parasitic effects play a key role.According to [23,24], most of the proposed designs are based on optimization routines, depending on the temperature rise and assuming uniform heat dissipation.Some authors consider the desired efficiency as a design input for the HFT [24]; otherwise, the optimization is focused on the A p area [25,26].In [27], the magnetizing inductance is considered the first step for the design process of HFT.In [28], The leakage inductance is first investigated, and then the remaining parameters of HFT are given.
This paper proposes a step-by-step methodology for designing medium and highfrequency transformers and their optimum electrical and magnetic specifications to reduce volume and increase efficiency.The proposed design methodology for the HFT first finds the optimum magnetic-flux density and then arrives at the other parameters according to the desired criteria such as leakage inductance and losses, as shown in Figure 4.

1.
System and Topology Specifications: Definition of the main system specifications, namely: apparent power (S), primary voltage (V P ), secondary voltage (V S ), operating frequency (f ), ambient temperature (T a ), admissible temperature variation (∆T), wave factor (K v ), window utilization factor (k u ) and stacking factor (k f ); 2.
Magnetic Core Properties: Selection of the material to be used for the magnetic core to find the saturation flux density (B sat ) and parameters of the constants α, β and Kc of the respective material, essential parameters for the following steps; most common cores, such as amorphous, ferrite, and nanocrystalline magnetic materials are compared in [29] for core losses at 10, 20, 50, and 100 kHz, respectively.

3.
Calculate Optimum Flux Density: The optimal flux density corresponds to the flux density value that minimizes the sum of copper losses with iron losses.Knowing the expression that describes the losses in copper and iron, as a function of the frequency and flux density, and after numerical simplification, it is possible to obtain the expression of the optimal flux density by Equation ( 1), as laid out in [30].
where k a , k c , k w are dimensionless constants, which usually assume the values 40, 5.6 and 10 respectively [30]; ∑ VA represents the apparent power of the transformer, h c is the heat transfer coefficient, and ρ w is the electrical resistance of the conductor that composes the windings.

4.
Magnetic Core Sizing: For the optimal flux density, determined in the previous point, the dimensions of the magnetic core of the transformer are selected based on the product of the window area (W a ) by the cross-sectional area (A c ); this product is given the nomenclature A p , usually expressed in cm 4 .For the analytical determination of this parameter, it is necessary to perform a comparison between the value of the optimal flux density (B opt ) and the saturation flux density value (B sat ).Based on the exposition in [30]: , the value of the parameter A p is given by Equation (2). (2) , the value of the parameter A p is given by Equation (3).
Based on the determined analytical value, A p , a magnetic core with an equal or greater value should be chosen.

5.
Calculate the Number of Turns: The number of turns of each winding can be described by Equation (7).
where the parameter A m represents the effective cross-sectional area of the core, which is given by the product between the physical cross-sectional area of the core (A c ) and the staking factor (k f ) [17,30].

6.
Conductor Sizing: For the dimensioning of the conductors to be used it is necessary to determine the current density (J o ).To determine this parameter, it is necessary to compare the value of the optimum flux density (B opt ) and the value of the saturation flux density (B sat ).Based on [30]: , the value of J o is given by Equation (8).
• If B opt > B sat , the value of J o is given by Equation (9).
For parameter A p , the value of the previously chosen magnetic core should be used.Where the parameters V c , V w represent the volume of the core and the volume of the windings respectively, these parameters can be calculated from Equations ( 10) and (11), respectively.Where l c is the average length of the magnetic part of the core, MLT is the average length of a coil, and W a is the area of the magnetic core window.
The l c and MLT length can be obtained using Equations ( 12) and ( 13) respectively, where dimensions a, b, d, f and u are obtained from the core dimensions, as in Figure 5.
To minimize the skin effect and the proximity effect, the choice of the Litz wire should consider the previously determined current density.However, the sizing of the internal structure of the Litz wire is highly complex.In [20,31], simplified methods for sizing Litz wire are presented.In this study, we used the data provided by the manufacturers, namely the relationship between section/number of single conductors and the type of construction.

7.
Estimation of Leakage Inductance: The estimation of leakage inductance can be obtained through Dowell's model, but this model was designed to calculate the leakage inductance in a transformer composed of foil conductor windings.Considering that the model under study uses Litz wire windings, it will first be necessary to perform the equivalence between Litz wire and foil conductors and subsequently estimate the leakage inductance value through Dowell's model.
The thickness of a sheet conductor, d eq , may be given by Equation (14), where d represents the diameter of each conductor that makes up the Litz wire.For the height of the foil winding, it will be considered that it occupies the entire height of the window, a situation that will allow simplification of the calculation.The equivalent section can be given by Equation (15), where h w and d w represent, respectively, the height and width of the winding composed of Litz conductors.Thus, the number of sheet conductors, m, will be given by Equation (16), where N represents the number of turns of each winding, and n represents the number of single conductors that make up the Litz wire [32].
Using Dowell's model, the estimate of the leakage inductance, referenced to the primary, is given by Equation (17).Where d g is the distance between windings and d i is the distance between consecutive sheet conductors.The indices "1" refer to the primary and the indices "2" refer to the secondary of the HFT [30][31][32][33].
The inductance factor, F L , is given by Equation (18).The parameter ∆ is given by the ratio between the thickness of the sheet conductor, d eq , and the skin effect, δ, using Equations ( 19) and (20).
Core and Winding Losses and Efficiency: The power losses in the iron can be translated by (21), and the value obtained by this equation will be expressed in W/m 3 [30,34].
where the parameter k i is given by Equation ( 22), the parameter T represents the period (inverse of the operating frequency), the parameter D represents the cycle factor, and the parameter ∆B represents the flux density variation.
α+1.354 (22) Particularizing for the case where a Dual Active Bridge, a DAB converter with a Dual Phase Shift (DPS) modulation technique is used, the variation of the flux density (∆B), could be translated by twice the maximum flux density (B max ), where the maximum flux density is given by Equation ( 23) [30,34].
The power losses in each winding, considering the harmonic components of the current, can be translated by Equation (24) [30].
The parameter I DC represents the value of the direct current, and the parameter I n represents the value of the alternating current for harmonic n.It is also necessary to know the value of the DC resistance (R DC ) and the value of the AC resistance (R AC ); the latter also depends on the frequency of the harmonic.The value of the DC resistance can be given by Equation (25), where ρ is the electrical resistance of the conductor.
The value of the AC resistance, for harmonic n, based on Dowell's equation is given by Equation (26).It is designed to foil type windings, so the power losses in the windings will be calculated for the equivalent foil winding determined in Topic 7 [35].
The term A str is described by Equation ( 27), which results from the ratio between sheet conductor thickness, d eq , and the skin effect, δ, for the operating frequency [35].
The efficiency of HFT can be determined by (28).
In HFTs, the accounting of parasitic parameters, namely parasitic capacitances, is fundamental for the correct modeling of the system.These capacitances are responsible for unwanted resonance frequencies and oscillations of the temporal evolution of the current and consequently a decrease in the reliability and efficiency of the system.
The parasitic capacitances in a winding are essentially between conductors of the same winding, of different windings, and between the conductors and the magnetic core.The three factors with the greatest contribution to the value of the parasitic capacitance are the dielectric constant of the insulating materials used, conductor geometry (diameter and size of the conductor, space between conductors and their distance from the magnetic core) and the characteristics of the winding (organization of conductors and winding direction) [36,37].
The equation for determining the C stray is given by Equation ( 29), based on the equivalent circuit of six capacitors and the simplification to a single equivalent capacitance [38].
The equivalent resistance of the transformer's magnetic core, R c , can be described by Equation (30), in which l m represents the average length that the magnetic flux describes in the magnetic core [40].
The equivalent inductance of the magnetic core of the transformer, L m , can be by Equation (31) [17].
Figure 6 shows the simplified equivalent scheme of an HFT, referenced to the transformer's primary [41].With the transition of a parameter from the secondary to the primary, as a rule, the same nomenclature is maintained, but a pellicle is added.1. System and Topology Specifications: Definition of the main system specifications, namely: apparent power (), primary voltage ( ), secondary voltage ( ), operating frequency (), ambient temperature ( ), admissible temperature variation (Δ), wave factor ( ), window utilization factor ( ) and stacking factor ( ); 2. Magnetic Core Properties: Selection of the material to be used for the magnetic core to find the saturation flux density ( ) and parameters of the constants ,  and  of the respective material, essential parameters for the following steps; most common cores, such as amorphous, ferrite, and nanocrystalline magnetic materials are compared in [29] for core losses at 10, 20, 50, and 100 kHz, respectively.3. Calculate Optimum Flux Density: The optimal flux density corresponds to the flux density value that minimizes the sum of copper losses with iron losses.Knowing the expression that describes the losses in copper and iron, as a function of the frequency and flux density, and after numerical simplification, it is possible to obtain the expression of the optimal flux density by Equation ( 1), as laid out in [30].
where  ,  e  are dimensionless constants, which usually assume the values 40 , 5.6 and 10 respectively [30]; ∑  represents the apparent power of the respectively.Where  is the average length of the magnetic part of the core,  is the average length of a coil, and  is the area of the magnetic core window.
The lc and MLT length can be obtained using Equations ( 12) and ( 13) respectively, where dimensions a, b, d, f and u are obtained from the core dimensions, as in Figure 5.To minimize the skin effect and the proximity effect, the choice of the Litz wire should consider the previously determined current density.However, the sizing of the internal structure of the Litz wire is highly complex.In [20,31], simplified methods for sizing Litz wire are presented.In this study, we used the data provided by the manufacturers, namely the relationship between section/number of single conductors and the type of construction.
7. Estimation of Leakage Inductance: The estimation of leakage inductance can be obtained through Dowell's model, but this model was designed to calculate the leakage inductance in a transformer composed of foil conductor windings.
Considering that the model under study uses Litz wire windings, it will first be The equivalent inductance of the magnetic core of the transformer,  , can be by Equation ( 31) [17].
Figure 6 shows the simplified equivalent scheme of an HFT, referenced to the transformer's primary [41].With the transition of a parameter from the secondary to the primary, as a rule, the same nomenclature is maintained, but a pellicle is added.

Converter Modeling
Based on the HFT modeling exposed above, it is possible to model the converter under study.Figure 7 represents the simplified power and control scheme, consisting of a Type D SST (see Figure 1), the only one between the two DC buses, that are the target of this study.

Converter Modeling
Based on the HFT modeling exposed above, it is possible to model the converter under study.Figure 7 represents the simplified power and control scheme, consisting of a Type D SST (see Figure 1), the only one between the two DC buses, that are the target of this study.The schematic depicted in Figure 7 also corresponds to the structure of a DAB, consisting of two single phase bridges, connected to each other via an HFT.To allow power bidirectionality, the semiconductors should be fully controllable, e.g., IGBT or MOSFET.Several control techniques can e used, which are described in more detail in [42,43].For the present study, rectangular modeling was used and is described below [42].

Control System
The implemented output voltage control establishes a phase shift angle control between the primary and secondary voltage waveforms, which is only applied to the second converter of the DAB.The first converter of the DAB assumes a null angle and is constant over time [42].
The purpose of the phase shift angle control is to regulate the power delivered to the load so that the voltage at the load terminals assumes a constant and predefined value, regardless of the load power, considering the power limitations supported by the system.Thus, the phase shift angle is determined by measuring the voltage at the load terminals.
Figure 8 presents the block diagram of the controller.The power converter is modulated with a gain,  , a delay time,  , and the load dynamics.Subtracting the reference voltage,  , from the measured voltage at the output of the converter,  , provides the value of the voltage error.The gain  depends on the configuration of the instrumentation used to perform the output voltage measurement, and for simulation purposes this parameter is assumed to be a unitary.The schematic depicted in Figure 7 also corresponds to the structure of a DAB, consisting of two single phase bridges, connected to each other via an HFT.To allow power bidirectionality, the semiconductors should be fully controllable, e.g., IGBT or MOSFET.Several control techniques can e used, which are described in more detail in [42,43].For the present study, rectangular modeling was used and is described below [42].

Control System
The implemented output voltage control establishes a phase shift angle control between the primary and secondary voltage waveforms, which is only applied to the second converter of the DAB.The first converter of the DAB assumes a null angle and is constant over time [42].
The purpose of the phase shift angle control is to regulate the power delivered to the load so that the voltage at the load terminals assumes a constant and predefined value, regardless of the load power, considering the power limitations supported by the system.Thus, the phase shift angle is determined by measuring the voltage at the load terminals.
Figure 8 presents the block diagram of the controller.The power converter is modulated with a gain, K d , a delay time, T d, and the load dynamics.Subtracting the reference voltage, V 0 Re f , from the measured voltage at the output of the converter, V 0 , provides the value of the voltage error.The gain α depends on the configuration of the instrumentation used to perform the output voltage measurement, and for simulation purposes this parameter is assumed to be a unitary., from the measured voltage at the output of the converter,  , provides the value of the voltage error.The gain  depends on the configuration of the instrumentation used to perform the output voltage measurement, and for simulation purposes this parameter is assumed to be a unitary.The controller   , has as output the phase shift angle, given by sin  , allows a zero error to be obtained, so that the output voltage,  , assumes the value of the reference voltage,  .The controller C(s), has as output the phase shift angle, given by sin(δ), allows a zero error to be obtained, so that the output voltage, V 0 , assumes the value of the reference voltage, V 0 Re f .
Based on the simplified equivalent circuit of the HFT, represented in Figure 9, where − V P and − V S represent the voltage at the primary and secondary terminals, respectively, performing the mesh analysis and after algebraic simplification, it is possible to obtain Equation (32).Based on the simplified equivalent circuit of the HFT, represented in Figure 9, where  and  represent the voltage at the primary and secondary terminals, respectively, performing the mesh analysis and after algebraic simplification, it is possible to obtain Equation (32).To design the controller C(s), it is necessary to establish the relationship between the input power of the DAB converter,  , and the output power of the DAB converter,  , which corresponds to the power delivered to the load.Assuming that the converter is ideal, to simplify the controller design, i.e., that the input power is equal to the output power, it is possible to derive the expressions exposed in Equation (33).
Since the controller (s) outputs sin  , and based on the equivalent controller scheme depicted in Figure 8, it is possible to conclude that the gain  will be given by Equation (34).
Considering the use of a proportional integral (PI) controller, starting from the diagram represented in Figure 8, consider  = 1 and cancelling the dominant pole with the controller zero, it is possible to obtain the closed chain transfer function, represented in Equation (35).
By simplifying Equation (35) and comparing it with the canonical second-order equation, it is possible to determine the expressions that define the proportional and integral gain, represented in Equations ( 36) and (37) [44], respectively.
Figure 9. Simplified HFT equivalent circuit, with the primary and secondary voltages represented.
To design the controller C(s), it is necessary to establish the relationship between the input power of the DAB converter, P in , and the output power of the DAB converter, P 0 , which corresponds to the power delivered to the load.Assuming that the converter is ideal, to simplify the controller design, i.e., that the input power is equal to the output power, it is possible to derive the expressions exposed in Equation (33).
Since the controller C(s) outputs sin (δ), and based on the equivalent controller scheme depicted in Figure 8, it is possible to conclude that the gain K d will be given by Equation (34).
Considering the use of a proportional integral (PI) controller, starting from the diagram represented in Figure 8, consider α = 1 and cancelling the dominant pole with the controller zero, it is possible to obtain the closed chain transfer function, represented in Equation (35).
By simplifying Equation (35) and comparing it with the canonical second-order equation, it is possible to determine the expressions that define the proportional and integral gain, represented in Equations ( 36) and (37) [44], respectively.
assumed that there is power bidirectionality.Figure 10 shows the generic block diagram that summarizes the various elements inherent to the application.The present study focuses only on the "Module", which consists of the two converters and the HFT, as represented generically in Figure 7.The main considerations assumed are as follows: • "Three-phase bridge rectifiers", fed from the conventional three-phase grid (400 V 50 Hz), composed of fully controlled semiconductors.For simulation purposes, consider an ideal DC voltage source of equal effective value; To simplify the analysis, the power requested by the "Charger/Battery" will be considered as a pure resistive load with equal power value.Consider that the charger should be powered with 1000 V DC.
Based on the flowchart presented in Figure 4 and the above considerations, Tables 2  and 3 show the main dimensioned parameters.The present study focuses only on the "Module", which consists of the two converters and the HFT, as represented generically in Figure 7.The main considerations assumed are as follows: • "Three-phase bridge rectifiers", fed from the conventional three-phase grid (400 V 50 Hz), composed of fully controlled semiconductors.For simulation purposes, consider an ideal DC voltage source of equal effective value;

•
To simplify the analysis, the power requested by the "Charger/Battery" will be considered as a pure resistive load with equal power value.Consider that the charger should be powered with 1000 V DC.
Based on the flowchart presented in Figure 4 and the above considerations, Tables 2 and 3 show the main dimensioned parameters.For the constitution of the windings, we used the characteristics of the conductors provided by the manufacturer New England Wire Technologies.The size of the windings must be compatible with the size of the magnetic core chosen; Table 3 summarizes this.
For the previous data, we obtain a value of theoretical power losses in the iron of 98.5W and in the copper of 526.5W, for calculation purposes until the 13th harmonic.This results in a theoretical efficiency of 98.8%.
The parameters of the HFT, which respond to the equivalent scheme presented 6, are shown in Table 4.The need to add a capacitor in parallel with the load and an inductance in series with the transformer were verified to improve the behavior/performance of the controller and consequently of the system to be modeled.

Numerical Simulation Results
To verify the performance of the scaled transformer model and the ability of the control to adapt, three tests were employed:

•
Test 1-Constant Load: Consider that the system feeds a constant, purely resistive load in the direct power direction and that the load has a power of 50 kVA.
By analyzing the waveforms of the voltage represented in Figure 11, it is possible to verify that the control system worked as expected, acting on the phase shift angle so that the voltage at the load terminals assumes the previously established reference value.The value of the phase shift angle depends essentially on the equivalent inductance of the model under study and the value of the current that the load requests, having as a maximum value +Ω/2.To evaluate the performance of the developed HFT, its performance can be analyzed.Figure 11 represents the average instantaneous power at the primary and secondary terminals.By using the measurement tools integrated in the software, Matlab/Simulink, it is possible to determine the obtained efficiency, having verified a value of 99.5% .Analyzing the obtained efficiency, we conclude that it does not consider the winding resistance variation as a function of frequency and does not account for the skin effect.When the theoretical efficiency is determined without accounting for the skin effect, a theoretical efficiency of 99.4% is obtained.
The analysis of the voltage waveform at the transformer terminals will be performed in Test 3.

•
Test 2-Variable load: Consider that the system is feeding a variable load, purely resistive, in the direct power direction.The initial load presents a power of 25 kVA and, at instant  = 1 s, the load doubles, resulting in a power of 50 kVA.All the analysis performed for Test 1 is valid for Test 2. The focus of the analysis described below is the impact of varying the value of the load on the system under study.
Analyzing Figure 12 it is possible to verify the occurrence of a voltage drop at the load terminals and a consequent controller action to guarantee the reference voltage at the load terminals.The voltage drop is significant, but it should be considered that the instantaneous doubling of the load is an extremely demanding transient for the system.Nevertheless, it is an excellent case to test the good performance of the controller.The shift angle after load variation and stabilization assumes approximately twice the value because the load has also doubled from 25 kVA to 50 kVA.To evaluate the performance of the developed HFT, its performance can be analyzed.Figure 11 represents the average instantaneous power at the primary and secondary terminals.By using the measurement tools integrated in the software, Matlab/Simulink, it is possible to determine the obtained efficiency, having verified a value of 99.5%.Analyzing the obtained efficiency, we conclude that it does not consider the winding resistance variation as a function of frequency and does not account for the skin effect.When the theoretical efficiency is determined without accounting for the skin effect, a theoretical efficiency of 99.4% is obtained.
The analysis of the voltage waveform at the transformer terminals will be performed in Test 3.

•
Test 2-Variable load: Consider that the system is feeding a variable load, purely resistive, in the direct power direction.The initial load presents a power of 25 kVA and, at instant t = 1 s the load doubles, resulting in a power of 50 kVA.All the analysis performed for Test 1 is valid for Test 2. The focus of the analysis described below is the impact of varying the value of the load on the system under study.
Analyzing Figure 12 it is possible to verify the occurrence of a voltage drop at the load terminals and a consequent controller action to guarantee the reference voltage at the load terminals.The voltage drop is significant, but it should be considered that the instantaneous doubling of the load is an extremely demanding transient for the system.Nevertheless, it is an excellent case to test the good performance of the controller.The shift angle after load variation and stabilization assumes approximately twice the value because the load has also doubled from 25 kVA to 50 kVA.
Figure 12 shows the variation of the average instantaneous power at the HFT terminals.Although this is a different test, after the stabilization of the load of 50 kVA, the performance values verified are the same as those in Test 1.When the load assumes a value of 25 kVA, an efficiency in simulation of 99.7% is verified.

•
Test 3-Power Bidirectionality: Consider that the system is supplying a purely resistive load with a power of 50 kVA in the direct power direction, and at instant t = 1 s, the power flow changes from the direct direction to the reverse direction, at which point the load is supplying the grid with a power of 50 kVA.
The shift angle is only applied to the bridge on the secondary side of the HFT, operating the bridge on the primary side of the HFT with a zero angle.When the power flow switches from the forward to the reverse direction, the phase shift angle will have to take on a negative value, as verified in Figure 13. Figure 12 shows the variation of the average instantaneous power at the HFT terminals.Although this is a different test, after the stabilization of the load of 50 kVA, the performance values verified are the same as those in Test 1.When the load assumes a value of 25 kVA, an efficiency in simulation of 99.7% is verified.

•
Test 3-Power Bidirectionality: Consider that the system is supplying a purely resistive load with a power of 50 kVA in the direct power direction, and at instant  = 1 s, the power flow changes from the direct direction to the reverse direction, at which point the load is supplying the grid with a power of 50 kVA.
The shift angle is only applied to the bridge on the secondary side of the HFT, operating the bridge on the primary side of the HFT with a zero angle.When the power flow switches from the forward to the reverse direction, the phase shift angle will have to take on a negative value, as verified in Figure 13.   Figure 12 shows the variation of the average instantaneous power at the HFT terminals.Although this is a different test, after the stabilization of the load of 50 kVA, the performance values verified are the same as those in Test 1.When the load assumes a value of 25 kVA, an efficiency in simulation of 99.7% is verified.

•
Test 3-Power Bidirectionality: Consider that the system is supplying a purely resistive load with a power of 50 kVA in the direct power direction, and at instant  = 1 s, the power flow changes from the direct direction to the reverse direction, at which point the load is supplying the grid with a power of 50 kVA.
The shift angle is only applied to the bridge on the secondary side of the HFT, operating the bridge on the primary side of the HFT with a zero angle.When the power flow switches from the forward to the reverse direction, the phase shift angle will have to take on a negative value, as verified in Figure 13.  Figure 12 also shows the variation of the average instantaneous power at the transformer terminals, where it is possible to see the power bidirectionality.The efficiency verified when the power flow is in the direct direction is equal to the efficiency described in Test 1.When the power flow is in the reverse direction, the same efficiency value is verified, 99.5%, and in this instance, the value of the average instantaneous power in the secondary is higher than the value of the average power in the primary.This allows the conclusion that regardless of the power flow direction, for the same power value, the losses in the transformer will be approximately equal.
The temporal evolutions of the voltage at the load terminals present the same temporal progression as described in Tests 1 and 2.
The impact of the variation of the shift phase angle can also be seen in Figure 14, where the temporal evolution of the voltage at the HFT terminals is represented.When the power flow is in the forward direction, the secondary voltage leads the primary voltage.On the contrary, when the power flow is in the opposite direction, it is verified that the secondary voltage is laggs the primary voltage.
in Test 1.When the power flow is in the reverse direction, the same efficiency value is verified, 99.5%, and in this instance, the value of the average instantaneous power in the secondary is higher than the value of the average power in the primary.This allows the conclusion that regardless of the power flow direction, for the same power value, the losses in the transformer will be approximately equal.
The temporal evolutions of the voltage at the load terminals present the same temporal progression as described in Tests 1 and 2.
The impact of the variation of the shift phase angle can also be seen in Figure 14, where the temporal evolution of the voltage at the HFT terminals is represented.When the power flow is in the forward direction, the secondary voltage leads the primary voltage.On the contrary, when the power flow is in the opposite direction, it is verified that the secondary voltage is laggs the primary voltage.

Conclusions
SST's are increasingly used to replace traditional transformers (50/60 Hz) because of the many advantages they provide, such as smaller dimensions, higher power density, ease of integration with alternative systems, supply of loads in direct current, and integration of systems for the improvement of power quality, among others.However, utilizing HFTs instead of low frequency transformers may present higher costs because numerous power components and auxiliary control circuits are required, which can be relatively complex and expensive.
Several simulation tests were performed in Matlab/Simulink on the application of the SST to a fast charger of an electric vehicle supplied by an electrical network.In these simulation tests, an efficiency of around 99.5% was obtained for a 50 kVA constant resistive load operating at 1kHz.Another test was also performed, in which the load was changed from 25 kVA to 50 kVA.Despite the extreme load variation, the controller presented good performance and fast stabilization.The maximum efficiency was 99.7% obtained at 25 kVA.The bidirectionality power flow of the solution was also tested and

Conclusions
SST's are increasingly used to replace traditional transformers (50/60 Hz) because of the many advantages they provide, such as smaller dimensions, higher power density, ease of integration with alternative systems, supply of loads in direct current, and integration of systems for the improvement of power quality, among others.However, utilizing HFTs instead of low frequency transformers may present higher costs because numerous power components and auxiliary control circuits are required, which can be relatively complex and expensive.
Several simulation tests were performed in Matlab/Simulink on the application of the SST to a fast charger of an electric vehicle supplied by an electrical network.In these simulation tests, an efficiency of around 99.5% was obtained for a 50 kVA constant resistive load operating at 1 kHz.Another test was also performed, in which the load was changed from 25 kVA to 50 kVA.Despite the extreme load variation, the controller presented good performance and fast stabilization.The maximum efficiency was 99.7% obtained at 25 kVA.The bidirectionality power flow of the solution was also tested and confirmed, showing a gradual transition from battery charging (Grid-to-Vehicle-G2V) to battery discharging (Vehicle-to-Grid-V2G) while the phase shift angle changed from positive to negative values to regulate the power delivered to the load.
According to the simulation results obtained, all our objectives were verified, such as the validation of the theoretical model of the HFT, the validation of the control scheme developed for SST control, and the verification of the ability to support power bidirectionality.Furthermore, temporal evolutions and efficiency values consistent with the theoretical concepts were obtained.

Figure 4 .
Figure 4. Methodology of the transformer design.

Figure 4 .
Figure 4. Methodology of the transformer design.

Figure 8
Figure8presents the block diagram of the controller.The power converter is modulated with a gain,  , a delay time,  , and the load dynamics.Subtracting the reference voltage,  , from the measured voltage at the output of the converter,  , provides the value of the voltage error.The gain  depends on the configuration of the instrumentation used to perform the output voltage measurement, and for simulation purposes this parameter is assumed to be a unitary.

Figure 9 .
Figure 9. Simplified HFT equivalent circuit, with the primary and secondary voltages represented.

Figure 11 .
Figure 11.Test 1: Temporal evolution of voltage, phase shift and power.

Figure 11 .
Figure 11.Test 1: Temporal evolution of voltage, phase shift and power.

Figure 13 .
Figure 13.Test 3: Temporal evolution of phase shift and power.

Figure 12 .
Figure 12.Test 2: Temporal evolution of voltage, phase shift and power.

Figure 12 .
Figure 12.Test 2: Temporal evolution of voltage, phase shift and power.

Figure 13 .
Figure 13.Test 3: Temporal evolution of phase shift and power.Figure 13.Test 3: Temporal evolution of phase shift and power.

Figure 13 .
Figure 13.Test 3: Temporal evolution of phase shift and power.Figure 13.Test 3: Temporal evolution of phase shift and power.

Figure 14 .
Figure 14.Temporal voltage evolution at of the terminal's transformer.

Figure 14 .
Figure 14.Temporal voltage evolution at of the terminal's transformer.

Table 1 .
List of variables and parameters.

. Modulation of the Converter System 2
.1.Solid-State Transformer

Table 3 .
Constitution of the windings.

Table 3 .
Constitution of the windings.

Table 4 .
The HFT equivalent parameters.