A fluid-thermal coupled analysis based on FEM is conducted. The inner structure of the coils is built with consideration of both the structural details and the simplicity; thus, the detailed heat conduction process is coupled with the computational fluid dynamics in the thermal computation of air-core reactors. According to the simulation results, 2D temperature distribution results are given and proved by the thermal test results of a prototype. Then the temperature results are used to calculate the heat flux to predict the detailed heat transfer process in the packages of the reactors. The study in this paper may be useful in the design optimization in air-core reactors.

With the development of power system, the usage amount of power reactors continues to rise more and more large capacity air-core reactors appear. Air-core reactors with large capacities have more packages and greater heights, which may cause uneven distribution of temperature rise. Thus, the accurate temperature field computational methods are in need.

Numerical simulation has been used to analyze the thermal condition of power equipment for a long time. Two categories of numerical approaches were developed as “Network Modeling” [

Air-core reactors have complex structures with paralleled cylinder packages and vertical air ducts. The heat transfer in them is a hybrid process with heat conduction, convection, and radiation. Thus FEM is obviously better suited than the other computational methods in the thermal computation of air-core reactors. In [

In this paper, a fluid-thermal coupled analysis based on FEM is conducted. The inner structure of the coils is built with consideration of both the structural details and the simplicity; thus, the detailed heat conduction process is coupled with the computational fluid dynamics in the thermal computation of air-core reactors. After the analysis, 2D temperature distribution results are given and proved by the test results of a prototype. Then the temperature results are used to calculate the heat flux to predict the detailed heat transfer process in the packages of the reactors. The study in this paper may be useful in the design optimization in air-core reactors.

The typical air-core power reactor contains some cylindrical packages connected in parallel and the vertical air ducts between them, as it is shown in Figure

Structure sketch of air-core reactor.

In this paper, a prototype of air-core power reactor with five packages is taken as an example, the parameters of which are given in Table

The parameter of an air-core power reactor prototype.

Package number | Diameter | Current | Ampere density | Height |
---|---|---|---|---|

(m) | (A) | (A/m^{2}) |
(m) | |

1 | 0.619 | 36.662 | 1.329 × 10^{6} |
0.325 |

2 | 0.684 | 31.081 | 1.329 × 10^{6} |
0.325 |

3 | 0.747 | 39.476 | 1.329 × 10^{6} |
0.325 |

4 | 0.811 | 44.049 | 1.329 × 10^{6} |
0.325 |

5 | 0.880 | 58.732 | 1.329 × 10^{6} |
0.325 |

The heat sources in the packages include the resistive current loss and extra eddy loss [

In most cases, heat transfer process contains heat conduction, convection, and radiation.

Partial differential equation of heat conduction is given as

Heat convection process is closely related to the fluid states, which can be described by the continuity equation, the momentum equation, and the energy equation as shown in (

The heat radiation can be described by

In air-core reactors, most of the packages’ surfaces face the surface of another package, on which the temperature is nearly the same. Thus, for air-core reactors, the radiating heat transfer can be neglected according to (

According to the mathematical modal mentioned above, a fluid-thermal coupled analysis should be adopted. A 2D axial symmetry finite element model in Flotran is built to conduct the analysis.

The simulation model is comprised of five packages and four vertical air ducts. Rectangles are chosen to replace the real cross sections of metal conductors and insulating materials to simplify the modeling. But, in order to study the detailed distribution in the packages, the metal conductor is divided into three parts in the radial direction and 15 parts in the axial direction, as shown in Figure

(a) Metallic conductors. (b) Insulating material.

In the rectangles that indicate the metal conductors, equivalent surface heat densities computed according to (

The insulating materials are also plotted as a combination of rectangles as shown in Figure

The boundary conditions of the analysis are set as follows.

Rotational symmetry conditions are defined on the left boundary line, the radial velocity of which is defined as 0.

The radial and axial velocities of all the fluid-solid interface in the analysis model are set 0. The roughness coefficient of the packages’ surfaces is set 0.2.

In the right boundary, the temperature is set 20°C, the relative pressure is set 1, and the velocities in all directions are set 0.

In the lower boundary, the temperature is set 20°C, the velocities in all directions are set 0, and the pressure is defined variable.

In the lower boundary, the temperature, the pressure, and the velocities in all directions are defined variable.

The turbulence model is open and the SST turbulence model in Flotran is chosen to solve the problem.

After the analysis, results of temperature field, fluid velocity and the pressure distribution, and so on are given. According to the abovementioned results the temperature distribution regulatory and the heat transfer processes of air-core power reactors are studied.

Temperature field of the analyzed area is given in Figure

Temperature field of an air-core reactor prototype.

The temperature field shown in Figure

Data of the temperature along the radial center of the packages are read and plotted as the curves shown in Figure

The axial temperature distribution along the midline of the packages.

Data of the temperature along the radial center of the vertical air ducts are read and plotted as the curves shown in Figure

The axial temperature distribution along the midline of the vertical air ducts.

In order to study the temperature distribution and heat transfer in the radial direction, 3 horizontal paths are defined, as the horizontal lines shown in Figure

Radial temperature distribution along the “Bottom,” “Middle,” and “Top” paths defined in Figure

According to Figures

In the axial direction of the packages, temperature rises from the bottom up. Due to the end effects the maximum temperature appears near 80% the package height.

Along the radial direction of the packages, the maximum temperature appears near the center of the package. Along the path from the center to the side surface, temperature declines slowly and nonlinearly, and the temperature gradient in the insulating material is bigger than that in the metal conductors.

In the air ducts, temperature rises linearly from the bottom up. In the radial direction, along the path from the center of the ducts to the packages’ side surfaces, temperature rises nonlinearly. A remarkable temperature gradient appears in the lamina near the fluid-solid interface, which is defined as thermal boundary layer. From Figure

Heat flux reflects the heat flow in the air-core reactor, which is a good tool to study the heat transfer process. Due to the complex construction, the heat flux inside the packages is hard to acquire. Thus, the heat flux passes through the surfaces of the packages are considered instead.

Using the temperature distribution results and conduction equation given in (

(a) The average heat flux and the axial heat flux through the inner and outer surfaces of the 3rd package. (b) The average heat flux and the axial heat flux through the inner and outer surfaces of the 5th package.

Axial heat flux through the bottom and top of the packages.

The results of

From Figures

Most of the heat load generated inside the packages are dissipated through the side surfaces of the packages and are carried away by the flowing air. A small part of the heat load dissipated through the bottom and top of the packages. The heat flux through the bottom surfaces is larger than that through the top surfaces.

For side surfaces that are facing the air ducts (inside surfaces), the envelope of the heat flux curves decline from the bottom up, which indicates that an axial heat flow goes from above to below inside the package. The mentioned axial heat flux is much smaller than the radial heat flux that goes through the side surfaces.

For inner packages, both the side surfaces of which face the air ducts, the heat flux through the two side surfaces equals each other. For outer packages, the side surfaces of which face the large air space (outside surfaces) and the air ducts (inside surfaces) separately, the heat convection in the air duct is stronger than that in the large air space, which indicates a heat flux from the outer part of the outer packages into the inside surface.

Near the entrance of the air ducts, the local temperature difference between the conductors and air is relatively small, but the local heat flux is considerably higher due to a stronger local heat convective effect.

The agreement of the thermal results between the analysis and the prototype testing proves the computational accuracy of the fluid-thermal coupled analysis.

The fluid-thermal coupled analysis gives more details than the prototype tests. The computational results of temperature field, fluid velocity distribution, and so forth can be used to study the field distribution of a design and thus are useful to judge its reliability.

By means of the data obtained in the fluid-thermal coupled analysis, the heat transfer process, including the heat conduction in the packages and the thermal conduction/convection process in the ducts and large air spaces, can be studied distinctly. As a result, the temperature distribution results could be explained. Furthermore, the study results of detailed heat transfer process in this paper could be used to optimize the heat load distribution in the packages, as well as the structural improvement of the packages and air ducts.

The authors wish to thank The Ministry of Science and Technology of the People’s Republic of China for financial support in the National Key Technology R&D Program (2009BAA19B00).