Analysis of Thermal Behavior of High Frequency Transformers Using Finite Element Method ()
1. Introduction
Magnetic components design plays a key role in achieving high efficiency, low volume and reasonable price of power electronic equipment. With the advent of higher switching frequencies and power densities in power electronic circuits, it is important to ensure that magnetic components, such as transformers and inductors, operate within the limits defined by the thermal specifications of the circuit. Temperature rise depends on the power losses through the laws of the heat transfer theory. In order to obtain an accurate value of the maximum temperature of the device during the design process, it is necessary to apply an accurate thermal model. It is necessary of a trade-off between thermal distribution accuracy (directions through heat transfer is considered, constant thermal properties, feedback with the magnetic model, steady state models, etc) and the complexity of the thermal model that is obtained.
Losses in the magnetic components are important design parameters. In many high-frequency designs, magnets are limited by their losses. Thus, it is extremely important for the designer to have a good practical model for estimating the losses under various excitations that are frequently encountered in power electronics design. There are two major components of the losses [1-9] in a magnetic component: the core loss (i.e., the losses in the magnetic material which is used as core) and the winding losses. High-frequency conductor losses for power electronics applications have been considered by a number of authors [8,9]. There are models provided by different manufacturers [10-12], that are synthesized in a simple expression for the calculation of the temperature rise. They provide the average temperature in the outer surface of the device.
Most of thermal models for magnetic components are analytical and assume 1D heat transfer [13-15]. These models are based on thermal networks and commonly consider only steady state, constant thermal properties, concentrated and uniform power losses and no thermal feedback with the electric properties. However, the thermal distribution is 2D/3D for most of cases, even if the magnetic field exhibits a 1D distribution.
The thermal design is usually somewhat neglected as it is often not clear which theory and coefficients should be used and experiments are time consuming. In this article, by designing a typical high-frequency transformer with ferrite core, its thermal behavior is investigated. Many parameters such as flux density, winding loss value, using a fan to cool the winding and core and thermal conductivity have effects on the ferrite temperature which is simulated with ANSYS software.
2. Calculating the Heat Generation Caused by Winding and Investigating Its Effective Parameters
Heat generating sources are transformer windings which generate heat by passing the current and this heat is transmitted to other parts. For a winding with resistance Rj, the power loss can be written as:
(1)
Which the jth winding resistance is as follow:
(2)
Where:
(3)
(4)
In Equation (3), MLT is mean length per turn (cm). By considering Equation (2) to Equation (4), the power loss of jth winding is as follow:
(5)
And the total power dissipation is the sum of the power losses for each winding:
(6)
The loss in winding is optimal providing that [16]:
(7)
According to Equation (7), it can be found out that the power loss or generated heat depends on parameters such as ρ (wire effective resistivity), Ku (winding fill factor), turns, winding current and MLT (mean length per turn). So, these are effective parameters on generation heat source. In a transformer, if the primary voltage is like Figure 1, then the following equations will be written:
(8)
(9)
By substituting Equation (9) in Equation (7), the result can be written as:
(10)
where:
(11)
Equation (10) consists of three parts in which part A is electrical characteristics and B and C include core and magnetic characteristics respectively.
The core loss is obtained from the following equation:
(12)
Kfe is core loss coefficient which is different for various frequencies. Ac is core cross-sectional area and lm is magnetic path length. The typical value of β for ferrite core is 2.6 < β < 2.8. Kfe is drastically increased by increasing the frequency. In addition, Kfe depends on core temperature and Bmax. The dependency of Kfe on frequency, Bmax and temperature can be obtained from core characteristics. In choosing the core when combined with different alloys, there has been always trade-off between saturated magnetic flux density and the core loss. Using materials with high Bsat leads to volume, size and price reduction. But these materials cause more core loss.
3. Transformer Characteristic for Full-Bridge Converter
To study the temperature rise in ferrite core and its effective parameters, a sample transformer for full bridge converter with two secondary windings and one primary winding is chosen which can be seen in Figure 2 [16].
Figure 2. Transformer for full bridge converter with two outputs and one input.
Table 1. Parameters of designed transformer for full bridge converter
The obtained parameters for this transformer are listed in Table 1. The core and winding losses calculated for transformer are equal to 230 mW and 3.89 W respectively. These losses can be used as heat generating source in modeling and simulating the transformer in ANSYS.
4. Transformer Modeling in ANSYS
In this part, transformer with two secondary windings is modeled in ANSYS software which is shown in Figure 3. The distance between each winding is 0.5 mm. the primary winding is 20 mm long and 2 mm wide which this cross-sectional area includes 22 turns primary winding. This area for secondary winding is 20 mm2 (1 mm × 20 mm) which is for one turn. Also, another secondary winding is 20mm long and 0.5mm wide which consists of three turns.
The electrical and thermal characteristics for different parts of transformer are listed in Table 2.
5. Simulation Results
In this simulation, it is assumed that the heat transfer between winding and core of the transformer has been naturally done (natural convection). So, the film coefficient varies from 10 W/m2∙ºC to 25 W/m2∙ºC when there is natural convection. The film coefficient is considered equal to 10 W/m2∙ºC. The total cross-section of area is 70 mm2 and the power dissipation in this area is 3.89W. Therefore, the value of watt per square meter of area is equal to 55570 W/m2 which is used for heat generation source in ANSYS. In Figure 4, temperature distribution
Table 2. Electrical and thermal characteristics for different parts of transformer
Figure 4. Thermal distribution in transformer.
through the transformer is shown. In this figure, the hottest point is primary winding and middle leg of ferrite core which is about 61.962ºC. By paying attention to this figure, it is clear that the core corners have lower temperature in comparison with other parts. It should be noticed that in this case, the secondary windings temperature are lower than that in primary one. The heat distribution in ferrite core is shown in Figure 5.
Thermal flux in this case is shown in Figure 6 and it is maximized in primary winding corners which are about 428.69 W/m2. Regarding this matter and according to Figure 4, it can be pointed out that the corners temperatures are lower than that in other parts.
5.1. Variation of Flux Density
In this part, the effect of flux density variation on temperature distribution in transformer has been investigated. For this purpose, the flux density is decreased from 80 mT to 120 mT. According to Equation (10), the power loss is in inverse proportion with the square of flux density and by this change, the power loss is decreased 2.25 times.
The heat distribution for this condition is shown in Figure 7 in which the hottest point is about 55.309ºC. In
Figure 7. Thermal distribution in transformer (Decreasing flux density).
this case, one of the secondary windings temperatures (neighbor of primary winding) is approximately equal to primary winding temperature (the highest temperature).
5.2. Variation of the Ambient Temperature
In this case, it is supposed that the transformer is used in 40ºC temperature. So, it is observed that the maximum temperature inside the transformer has been slightly increased and has reached to 61.971ºC. In this case, the maximum temperature is in primary winding and core middle leg as well.
By comparing Figure 4 with Figure 8, we can make the conclusion that regarding the proximity of core sides to heat source and warmer environment, the sides’ temperature are higher than that in up and down of core.
5.3. Using a Fan to Cool the Winding
The film coefficient varies from 50 W/m2∙ºC to 120 W/m2∙ºC when there is forced convection. For this case the film coefficient has been changed from 10 W/m2∙ºC to 50 W/m2∙ºC. The temperature distribution is shown in Figure 9.
By using the fan, almost all windings have the same temperature of 52.508ºC. By comparing this figure with Figure 4, it can be understood that using the fan leads to winding temperature reduction from 61.962ºC to 52.508ºC. By this analysis, the appropriate fan can be chosen.
5.4. Using a Fan to Cool the Core
In this case, it is supposed that core is cooled by a fan (the film coefficient is 50 W/m2.ºC). The heat distribution in this condition is plotted in Figure 10. According to this figure, the core ambient temperature is almost
Figure 8. Thermal distribution in transformer (Increasing the ambient temperature).
Figure 9. Thermal distribution in transformer (Using fan to cool the winding).
equal. This figure shows that changing the film coefficient causes suitable increase in core temperature (61.95ºC).
5.5. Variation of Thermal Conductivity
On of the other studied parameters is thermal conductivity that it is changed from 0.004 W/(m∙K) to 0.008 W/(m.K). Thermal distribution for this mode is shown in Figure 11.
This figure shows that the thermal conductivity has a little effect (61.95ºC) on the temperature but the thermaldistribution around the transformer is uniform as a circle from centre of transformer.
6. Conclusions
In this paper, the thermal behavior of different cores has been studied. For this purpose, first for a typical trans-
Figure 10. Thermal distribution in transformer (Using fan to cool the core).
Figure 11. Thermal distribution in transformer (Increasing thermal conductivity).
former, the generation heat is computed then it is totally simulated in ANSYS software. In this simulation, the thermal analysis and heat distribution were studied. Then, the effects of parameters such as flux density, winding loss value, the film coefficient have been investigated.
According to obtained results, in case of using a fan to cool the winding, the maximum temperature of transformer has been decreased. Furthermore, it was shown that by decreasing the flux density, the temperature of the hottest point in transformer has been decreased which is due to the power loss reduction in windings. Moreover, it is explained in the paper that the core temperature can be controlled by increasing the film coefficient (by using a fan). As a result, with the findings of this analysis, an appropriate fan can be chosen in order to cool the whole transformer. In the paper, it is also elaborated that the core corners have the lower temperature in comparison with other parts. This part has the maximum thermal flux as well.
Nomenclature