Comparative Study of Simulation of Temperature Rise in Ring Main Unit
Summary
This paper addresses the temperature rise problem of the ring main unit by establishing a simplified model and conducting a comparative analysis in two temperature field-solving modules.
Abstract
The Ring Main Unit (RMU) is a critical device in power distribution systems used for connecting and distributing electricity. However, due to its compact internal structure and high current load, heat dissipation issues are particularly prominent. To address this problem, this study innovatively proposes a simplified RMU model, employing finite element simulation methods to accurately solve for the ohmic losses of conductors under actual operating conditions and obtain ohmic loss data for various components. This is the first in-depth investigation of the RMU’s temperature rise problem using such a comprehensive approach. Subsequently, the temperature field was solved using two different temperature field analysis modules, with a detailed comparison and analysis of the simulation results to identify similarities, differences, and trends in temperature distribution. The results indicate that the temperature field solution model, which considers convective heat transfer, is more accurate and aligns with actual operating conditions. This research provides an innovative approach and practical solutions for the design and optimization of RMUs. Future research can further explore multiphysics coupling analysis methods to address structural design and mandatory validation issues for high and ultra-high voltage RMUs and other electrical equipment, thereby providing important insights for engineering design.
Introduction
The ring main unit is a group of high-voltage switchgear mounted in a steel metal cabinet or made of assembled spaced ring network power supply unit of electrical equipment. The overall structure of the load switch and conductive circuit consists of the conductive circuit, which includes a number of components comprising the main core of the ring unit. However, due to its compact internal structure, the ring main unit faces challenges in heat dissipation. This can lead to thermal deformation and aging when operating for extended periods in high-temperature environments. These issues not only affect the service life of the unit but also impact its insulating properties, posing safety risks. In particular, equipment damage and electrical accidents become more likely, posing significant safety hazards.
Within different research areas, scholars have conducted a series of studies on the temperature rise of overhead line switchgear and analyzed various factors affecting the temperature distribution1. In Polykrati et al.2, a mathematical model for the estimation of the temperature rise of components installed on the distribution network during a short-circuit fault is presented. The model was applied to the common disconnecting switches of the network, and the characteristics of the results were plotted according to the different forms of the asymmetrical part of the short-circuit current waveform and the initial value of the short-circuit DC current component. Guan et al., on the other hand, have taken into account the contact resistance and electromagnetic repulsion by building an equivalent contact bridge to simulate the contact interface and further analyzed the electromagnetic-thermal coupling field and temperature rise experiment3. In addition, the researchers investigated the temperature field and thermal stress distribution of the dynamic and static contacts inside the ring main unit by finite element simulation, which provided a basis for the study of circuit breaker life4. Finally, Mueller et al. have focused on the geometrical characteristics of heat sinks and evaluated the effects of material selection, total surface area, temperature uniformity, and maximum surface temperature on thermal performance5. These studies provide valuable insights and methods to improve switchgear performance and reliability, reduce temperature rise, and extend equipment life. Wang et al. proposed a MiNET Deep Learning Model (MDLM) in the UPIOT environment with the purpose of detecting fault diagnosis of electrical ring cabinets, which was validated to have an identification accuracy of 99.1%, which is significantly higher than that of other methods6. Lei et al. studied the thermal performance of a GIS busbar in a steady state using the magneto-fluid-thermal coupling analysis method, thereby optimizing the conductor and tank diameter based on the temperature rise simulation results7. Ouerdani et al. used the RMU temperature rise simulation model to determine the temperature rise at critical locations inside it, thereby fixing the duration of the maximum overload for the components inside the RMU accordingly8. Zheng et al. described a conventional rectangular busbar in a model of high-current switchgear by building a two-dimensional model and applying the finite element method (FEM) for electromagnetic field calculations. It enabled them to obtain the distribution of bus conductor current density and power loss. An irregular busbar was designed after considering the effects of proximity effect and skin effect. This irregular busbar design improved the performance of conventional rectangular busbar9.
As for the aspect of using the icepak simulation, Wang et al. carried out a temperature rise simulation through vortex field, airflow field, and temperature field theories and found that the temperature rise of the ring main unit was more serious under natural convection. They successfully reduced the temperature rise level by adding forced air cooling and making improvements to the internal contact structure10. Zhu et al.11 used the icepak to simulate a thermal model in order to compare the effect of the presence of thermal vias on the PCB and the presence of heat sinks on the temperature of the power devices. Finally, the theoretical analysis is compared with the simulation results to verify the correctness of the theoretical analysis. Mao et al.12 studied the temperature and internal airflow distribution under summer operating conditions by thermal simulation based on the CAE software in the icepak simulation. The problem of how to improve the cooling efficiency and control the temperature rise of multiple silver-plated contacts is given, and the temperature and internal airflow contours captured in the simulation will lay the foundation for the design of the cooling scheme for the six silver-plated contacts mounted in the sealing unit. Conversely, in the use of a steady-state thermal module, Zhang13 Modeling methods are discussed for solving the thermal network of a high-pressure bushing using an alternative transient procedure. Test and simulation results are in good agreement with the thermal steady state and transient states of the bushing. The transient results are then used to evaluate the bushing overload capacity. Vaimann et al.14 developed and analyzed an analytical thermal model of a synchronous reluctance motor for predicting the temperature of its different components and the set total parameter thermal network.
With the continuous advancement of research on electrical equipment such as ring main units, conventional temperature rise tests, and production methods are relatively inefficient. Therefore, by utilizing finite element technology combined with offline tests, not only the design cost issues are addressed, but adjustments and optimizations can be promptly made to real-world problems based on simulations. Based on the research progress mentioned above, the use of ANSYS Icepak and Steady-state thermal coupling for comparative analysis is rarely mentioned. Therefore, the protocol describes the mechanism research of finite elements, uses numerical and morphological combinations to establish a finite element temperature rise simulation model for the enclosure, and discusses the finite element temperature rise simulation model based on the results of the two analytical modules by comparing the results of the two simulation modules. Through the comparison between the two simulation modules, we will get the characteristics of the temperature rise trend of the ring main unit and find the most applicable method so as to provide the necessary basis and research ideas for a strategy to mitigate the temperature rise of the ring main unit.
Protocol
1. Model
NOTE: Due to the complex structure of the ring main unit (Figure 1A), an online design software was chosen to simplify the operation of the ring main unit.
- Modelling simplification
- Partially simplify the model, preserving the air box section of the RMU while removing or simplifying other components such as insulating shafts, fastening bolts, nuts, sealing components, and pressure support brackets. The simplified version is shown in (Figure 1B).
- In the process of simplifying the 630A type ring main unit, remove the insulated shaft connecting the circuit breaker room with the instrument box and many fixed bolts and nuts. Take out the sealing parts and pressure-retaining bracket and connect the static contacts of the isolated static beam with the lower branch busbar under the premise of ensuring that the whole setup has the same conductive current and only the vacuum circuit breaker, the circuit breaker fixing plate, and the static contacts and the vacuum circuit breaker are retained.
- Retain only the vacuum circuit breaker, circuit breaker fixing plate, static contact, and vacuum circuit breaker blocking plate. Overall, remove bolts and gaskets from the model, fill the holes after removing the bolts with solids, reduce the number of mesh parts, and optimize the irregular shapes of the parts. Remove instruments for panel operation, mounting plates, brackets, and other operating parts, such as instrument boxes, that have no effect on the temperature rise simulation process.
- Remove the insulated housings of some components can be disregarded under simulation as they have little effect on the simulation results. In addition, grounding switches that have no effect on the use of the equipment during normal operation remove them, and retain the circuit breaker room for simulation.
- To delete any section, simply select it and click on the Delete option.
- Partially simplify the model, preserving the air box section of the RMU while removing or simplifying other components such as insulating shafts, fastening bolts, nuts, sealing components, and pressure support brackets. The simplified version is shown in (Figure 1B).
2. Eddy field solution
- Pre-processing settings
NOTE: The Eddy current field emulation is the basis for performing the temperature field solution, which requires the subsequent analysis of the solved heat source as a load on the temperature field.- Refer to the equipment documentation for the ring main unit and the relevant manuals to gather information on the physical properties and parameters of each component of the ring main unit. Set the physical attributes and parameters of the ring main unit components in Maxwell based on the information obtained, as detailed in Table 1.
- Set the eddy current field load current at 630 A with a frequency of 50 Hz. In the Maxwell software, select One Side of the Upper and Lower Outgoing Arms, enter the excitation module, and set the current magnitude to 630 A. In the solution settings section, choose a frequency of 50 Hz.
NOTE: In a ring main unit's conductive circuit, the pathway formed by all the components from the upper outlet arm to the lower outlet arm is known as the phase sequence. Therefore, in this paper, the phases A, B, and C are arranged from left to right. - The material parameters of the components of the ring main unit are shown in Table 2.
- Direct the current through the outgoing line arms, flexible connections, busbars, circuit breakers, static contact support busbars, and branch busbars for each phase. The aim is to realize a current path that allows the components to complete the load.
- Utilize Maxwell's adaptive meshing to complete grid control for the model. Use the Maxwell adaptive mesh partitioning method for larger components and the local mesh refinement for smaller internal components.
NOTE: Maxwell can continuously enhance grid precision during the solving process, eliminating the need to click Mesh Operations for additional mesh partitioning. - Set the solution step size. Click on Analysis in the model tree, open the Solve Step settings, and set the Maximum Number of Passes to 10. Keep other settings at their default values without making any changes.
- Principle of Eddy current field calculation15,16.
- Use Maxwell's first equation, which describes the action of charge on the generation of an electric field17.
(1)
where ρ represents the charge density; ε0 represents the vacuum dielectric constant. - Use Maxwell's second equation, which describes the relationship between a changing magnetic field and an electric field and the effect of a magnetic field on the motion of a charge.
(2)
where
represents the magnetic field strength. This equation describes that a varying magnetic field produces a vortex electric field, i.e., the spin of the vortex electric field is equal to the negative of the rate of change of the magnetic field with time. - Use Maxwell's third equation, which describes the effect of magnetic charge on the production of a magnetic field.
(3)
This equation describes the magnetic field produced by a magnetic charge as being passive, i.e., there are no monopoles in the magnetic field. - Use Maxwell's fourth equation, which describes the relationship between a varying electric field and a magnetic field and the effect of an electric current on a magnetic field.
(4)
where
represents the current density and μ0 represents the vacuum permeability. This equation describes that a varying electric field produces a vortex magnetic field, i.e., the spin of the vortex magnetic field is equal to the sum of the current density and the rate of change of the electric field with time. - Based on the above equations, use the Maxwell 3D using Eddy current solver module to solve the ohmic loss generated by the conductive circuit in the RMU, which provides a heat source for the subsequent thermal simulation analysis. Its mathematical expression is given as18
(5)
where σ denotes the conductivity of the conducting loop material; J is the current density in the loop.
- Use Maxwell's first equation, which describes the action of charge on the generation of an electric field17.
- Calculation results
- Click on the Maxwell 3D option in the interface and open the validation check to review all settings for errors. If there are no errors, proceed to click Analyze All to initiate the solving process.
- Utilize Maxwell's post-processing Calculator to compute and plot the ohmic losses in the Eddy current field of the ring main unit, as shown in Table 3.
3. Temperature field solution
NOTE: For comparative purposes, divide the temperature field into Icepak and steady state thermal. Set up and solve each separately to achieve a comparative analysis.
- Icepak model setup
- Set the material properties as follows: designate all circuit solid materials as Cu-Pure, with surfaces using Cu-polished-surface. For the panel components, select Aluminum6061-T6 material, with a surface coating of Paint-AL surface with an emissivity of 0.35. See Table 4 for details. Right-click on the Selected Component, click Edit and then go to Properties to set the material for both the surface and solid materials.
- Select the model and click on Set in the Edit menu, then choose Multilevel Meshing Level to adjust the mesh settings. Set the external cabinet to a mesh level of 2, and all boundaries to a mesh level of 2. For all other components, set the mesh level to 3. Finally, open Mesh control and click Generate to create the mesh.
- To ensure the accuracy and efficiency of the simulation regardless of the grid size, validation of grid independence is necessary. Import the geometric model of the temperature field enclosure, established using the design software for meshing.
- As depicted in Figure 2, the polarization curves of the four grid sets are well-aligned. At a working voltage of 0.5 V, the current densities for the four grid sets are 2.357 A/cm2, 2.358 A/cm2, 2.356 A/cm2, and 2.454 A/cm2, respectively, with the error between the maximum and current densities being less than 1%. To balance efficiency and accuracy, determine the grid size which comes to be 987924.
- Solution setup
- Set the directions of the solution domain Cabinet to Opening.
- In the software, select Problem Step. Under Basic Parameters, check the Surface-to-Surface Radiation Model, choose Zero Equation for Turbulent Flow Regime, select the Gravity option for Natural Convection, and set the ambient temperature to 20 °C.
- In the File settings, choose Volumetric Heat Losses for EM Mapping and select All Objects Shown to complete the loss settings.
- Temperature field calculation
- In icepak, apply three main conservation equations for energy: mass conservation equation, momentum conservation equation, and energy conservation equation. Specifically, use the momentum conservation equation, which is as follows19:
(6)
Energy Conservation Equation:
(7)
Mass conservation equations:
(8)
Energy transfer equation for heat transfer from a solid heat source:
(9)
ρ represents the density of the fluid; v represents the flow velocity vector; T represents the temperature; p is the pressure; τ is the viscous force on the surface of the micro metabolite; κ is the heat transfer coefficient; Sh is the body heat source; h is the specific enthalpy of the fluid and F is the body force of the micro metabolite.
NOTE: The results of the temperature field calculations are shown in Figure 3A and Figure 4A.
- In icepak, apply three main conservation equations for energy: mass conservation equation, momentum conservation equation, and energy conservation equation. Specifically, use the momentum conservation equation, which is as follows19:
- Steady-state thermal model setup
- Maintain the material properties as per Table 3 in the material settings. Generate the ohmic losses resulting from the eddy current field simulation analysis in the Steady State Thermal module by clicking on Thermal Load generation.
- Click on the Convective Temperature Value and set it to 20 °C, with a convective coefficient of 5 (W/m²°C) applied to the inner walls of the cabinet, components, and external cabinet. Apply the settings and generate. Set the output to solve for temperature by clicking on Solve > Output Results.
NOTE: The main temperature field governing equations in the steady-state thermal temperature field calculation20,21,22 principle is usually derived from the law of heat conduction (Fourier's law of heat conduction). In the one-dimensional case, the temperature field heat transfer equation can be expressed as20:
(10)
In this equation, T represents the temperature inside the object, t is time, x is spatial coordinates, and α is the thermal diffusivity. This equation describes the variation of temperature with respect to time and space, where the right side expresses the relationship between heat conduction rate and temperature gradient. In a more general three-dimensional scenario, the heat conduction equation for the temperature field can be expressed in the following form:
(11)
ρ represents the density of the object, c is the specific heat capacity, K is the thermal conductivity and Q is the heat source term within the volume. This equation describes the variation of the temperature field, influenced by heat conduction, heat sources, and thermal capacitance. - The temperature field calculation results are shown in Figure 3. Compare the temperature values summarized in Table 5 and Table 6.
Representative Results
Based on the data in Table 3, the following conclusions can be drawn: The overall losses for Phases A, B, and C are relatively similar. Specifically, the total losses for Phase A are 16.063 W/m³, Phase B is 16.12 W/m³, and Phase C is 19.57 W/m³. The locations with higher losses may be at the connections of various components. This is mainly because contact resistance and conductor resistance typically exist at these connection points. When current passes through these connections, significant heat is generated, leading to an increase in temperature and higher losses in those areas.
The upper and lower outgoing arms of the ring main unit bear some losses, especially when carrying the main load. The losses in the outgoing arms of the three phases are approximately the same. This is because this part of the current is relatively concentrated, and due to the regular shape, the resistance value is small. Therefore, the losses in these parts are approximately equal.
The losses in the branch busbar are relatively high, primarily due to its numerous bends and the presence of angled sections. Most of the current is concentrated in the bend area and near the corners. Apart from the vacuum circuit breaker section, the losses in the copper tube part of the busbar are also relatively high, totaling 22.32 W/m³, accounting for 40% of the overall ohmic losses of the three phases. The losses in the static contact part are relatively small compared to the overall losses.
Other components, such as the circuit breaker fixed plate, circuit breaker baffle, and the outer shell of the ring main unit, have smaller losses. Since they do not directly participate in the load, their losses are mainly generated by the internal components through conduction, resulting in smaller loss values. In post-processing calculations, the losses of these components are not particularly detailed. In summary, the main characteristics of the loss distribution have been outlined, providing valuable insights for further design optimization.
Combining the results from Figure 3 and Figure 4, along with Table 5 and Table 6, a thorough comparison of the temperature field solution results between the Icepak and Steady-state Thermal modules was conducted. In the analysis process, differences and similarities were observed between the two modules.
In the temperature range of 20-31.24 °C, as shown in Figure 2A, the higher temperature areas are concentrated in the regions where the RMU enclosure contacts the internal components. The main reason is that these contact areas typically serve as critical paths for heat conduction. During operation, the internal components generate heat, which is conducted to the enclosure through the contact surfaces with good thermal conductivity, thereby causing an increase in temperature in these regions. These areas include the contact points of the upper outgoing line arm and the portions of the lower outgoing line arm in contact with its nearby shell, with a temperature range of approximately 24-27 °C. Compared with Figure 3B, the warming range of the shells in both cases is roughly the same, and the temperature changes propagate from the hotspots to the cooler areas. Examining Figure 4A,B, the overall temperature trends in both cases are also approximately similar. The primary source of temperature rise is the heat conduction through the lower outgoing line arm, which passes through the internal conductors and dissipates from the upper outgoing line arm. The warming trend is concentrated in areas close to the lower outgoing line arm, including the branch busbars of each phase and the connected supporting busbars. Furthermore, from the observation of the temperature distribution in the internal circuits in Figure 4, whether using Icepak or the Steady-state Thermal module for solving, the overall temperature of Phase B is consistently higher than the other two phases. This mainly suggests that Phase B, during the load process, not only bears the heat generated by the current of that phase but also, due to the overly compact internal structure of the RMU, the heat produced by each phase cannot dissipate promptly. Due to the exchange of heat, the temperature of Phase B remains higher than the other two phases regardless of the solver used. The overall temperature trends also exhibit a high degree of similarity, with the main temperature rise stemming from the heat conduction through the lower outgoing line arm and the dissipation of heat from the upper outgoing line arm. This trend appears highly consistent between the two modules. Moreover, the positions of the hotspots in the temperature distribution maps of both modules are highly consistent. Particularly, in the branch busbar section of Phase C, both solvers show the same highest temperature and the differences are almost negligible. This implies that regardless of the solver used, accurate identification of the hottest temperature positions is crucial for the thermal management of the RMU.
Secondly, it can be seen from the table that the same current produces the same loss, but the temperature value of each phase component in Table 6 is generally lower than that in Table 5; for example, the temperature value of the branch bus in the hot spot area of the entire RMU is the highest temperature in Table 6 at 30.91 °C, and on the other hand, the highest temperature of the branch bus in Table 5 is 31.24 °C. The reason is the solution logic of Icepak: the RMU as a heat source will continue to exchange heat and dissipate heat, and when a fluid is set as a medium for heat transfer, the temperature in the solution domain will gradually be emitted through the medium, which results in a decrease in temperature. The Steady-state Thermal module places less emphasis on heat convection with the surrounding air and instead focuses on a heat conduction-based solution. Compared to models that account for convection, this approach results in a less comprehensive temperature field solution. Consequently, in areas where the temperature is higher, hot spots become more pronounced. In the process of the actual temperature rise test, in addition to considering the characteristics of the model itself, heat transfer at the same time also considers the influence of heat dissipation medium, such as air, etc., on the overall temperature. Combined with the practical needs and experimental operation considerations, the RMU as a heat source continues to carry out heat convection and heat dissipation; the temperature simulation of icepak is more in line with the actual needs. In the case where the ring main unit is used as a heat source with continuous heat convection and dissipation, the temperature simulation of icepak is more in line with the practical needs. On the contrary, the steady-state thermal module mainly focuses on heat conduction, which may not be able to meet the actual needs in some cases.
Figure 1: Ring main unit model. (A) Overall model of ring main unit (B) Simplified model of ring main unit. Please click here to view a larger version of this figure.
Figure 2: Polarization curves for four sets of meshes. Maxwell adaptive meshing is used for larger parts, and local mesh refinement is used for smaller parts. In design software, a simplified 3D model of the ring main cabinet is created and imported into Maxwell, and the model is meshed using its Mesh module. The figure represents the verification of the model grids. The current density of the four groups of grids is 2.357 A/cm2, 2.358 A/cm2, 2.356 A/cm2 and 2.354 A/cm2 under the voltage loading, and the relative error between the maximum current density and the minimum current density is less than 1%, and in order to take into account the efficiency and accuracy of the calculations, the final number of the grids is determined to be 1169091. Please click here to view a larger version of this figure.
Figure 3: Outer cabinet temperature field analysis model. (A) Shell temperature distribution with Icepak solution. (B) Shell temperature distribution with Steady-state thermal solution. Please click here to view a larger version of this figure.
Figure 4: Internal circuit temperature field analysis model. (A) Temperature distribution of the conducting loop with Icepak solution (B) Temperature distribution of the conducting loop with Steady-state thermal solution. Please click here to view a larger version of this figure.
| Material | Air | Galvanized steel sheet | Copper | Aluminum | ||
| Specific heat capacity (J/(kg·K)) | 1007 | 500 | 500 | 897 | ||
| Density (kg/m³) | 1.1614 | 8030 | 8900 | 2689 | ||
| Relative permeability | 1 | 0.3 | 0.9999991 | / | ||
| Relative magnetic permeability | 1 | 2500 | 1 | 1 | ||
| Thermal conductivity (W/(m·K)) | 0.026 | 16 | 386 | 237 | ||
| Emissivity | / | 0.65 | 0.3 | 0.1 | ||
| Electrical conductivity (S/m) | / | 0.8 | 5.80E+07 | / | ||
Table 1: Physical parameters of some materials in the ring main unit.
| Components | Relative permeability | Conductivity/s | Material |
| Current-carrying circuits (busbar, branch busbar, etc.) | 0 | 3,00,00,000 | Aluminium alloy |
| Vacuum circuit breaker | 0.9999 | 2,00,00,000 | Copper alloys |
| Body | 200 | 11,00,000 | Structural steel |
Table 2: List of materials for each component.
| Ohmic loss W/m³ | A | B | C |
| Out let arm | 3.78 | 3.72 | 3.73 |
| Branch busbar | 2.1 | 2.09 | 2.1 |
| Flexible connection | 1.3 | 1.3 | 1.3 |
| Vacuum circuit breaker | 1.023 | 0.95 | 0.98 |
| Static contacts | 0.36 | 0.36 | 0.36 |
| feeder busbar | 1.33 | 1.35 | 1.32 |
| Copper tubes for branch busbars | 6.19 | 6.35 | 9.78 |
Table 3: Ohmic loss values for phase A, B, and C components.
| Number | 1 | 2 |
| Materia | Aluminum6061-T6 | Cu-Pure |
| Thermal conductivity (W/m·K) | 167 | 387.6 |
| Density (kg/m³) | 8933 | 2700 |
| Specific heat capacity/(kg·K) | / | 896 |
| Surface material | Paint-AL surface | Cu-Polished-surface |
| Emissivity | 0.35 | 0.052 |
Table 4: Material parameter settings.
| Temperature monitoring point/temperature (°C) | A | B | C |
| Upper outlet arm | 25.48 | 25.79 | 25.63 |
| Main busbar | 25.93 | 26.28 | 26.13 |
| Branch busbar | 29 | 29.18 | 30.01 |
| Branch busbar copper tube | 31.04 | 31.18 | 31.24 |
| Lower outlet arm | 26.5 | 26.98 | 26.92 |
Table 5: Temperature values of each phase of the ring main unit under the Steady-state thermal temperature field solving module.
| Temperature monitoring point/temperature (°C) | A | B | C |
| Upper outlet arm | 23.73 | 23.82 | 23.81 |
| Main busbar | 25.15 | 25.17 | 25.35 |
| Branch busbar | 27.76 | 28.04 | 29.07 |
| Branch busbar copper tube | 28.42 | 29.31 | 30.91 |
| Lower outlet arm | 24.95 | 24.85 | 26.33 |
Table 6: Temperature values at monitoring points of each phase of the ring cabinet under Icepak temperature field solution module
Discussion
This paper is a comparative simulation analysis of the temperature rise of the ring cabinet based on engineering modeling software and finite element software, and the most suitable solution for the actual temperature rise situation is analyzed by two finite element temperature field solution modules. Thermal management is also described in Icoz23 as a critical and essential component in maintaining the high efficiency and reliability of electronic components. The significance of conducting a comparative analysis is summarized by drawing on the work in Steiner24: a comparative analysis has been carried out using COMSOL and ANSYS Mechanical. Therefore, in the process of manufacturing the actual RMU, its temperature distribution can be analyzed by finite simulation software, which can greatly save manpower and production costs.
The three-dimensional model of the RMU can be created using SolidWorks, as shown in Figure 1A. In Step 1.1, model simplification is emphasized as a crucial step in finite element analysis23. Given the numerous components of the RMU, which may impact solution accuracy and computational results, unnecessary parts are removed after considering only the conductive components. The main working parts are retained, as illustrated in Figure 1B.
In the preprocessing stage of the Eddy current field, a critical aspect involves the addition and verification of load excitation. As described in Step 2.1, it is essential to accurately add the load current to the upper and lower outgoing arms, ensuring that the entire conductive circuit forms a path for the current. This facilitates the calculation of the distribution cloud map of the eddy current field. The difficulty lies in the need to check the circuit before the solution calculation, as an incomplete current path in the circuit can lead to non-convergence during solving or situations where Maxwell cannot compute. Subsequent temperature field analysis relies entirely on the losses obtained from solving this eddy current field, making the inspection of the current path a necessary step.
In the temperature field solution, Step 3 presents two different modules for solving the model. However, ensuring identical initial conditions during the solution process poses a challenge. Since the modules have different emphases in their solutions, it is necessary to approximate or set identical environmental temperatures, convective coefficients, and other solving conditions to ensure consistency under the unique heat source provided by Maxwell. Furthermore, controlling the intervals for temperature field solution in both models to be consistent allows for a direct comparison of temperature displays for the same components in the same temperature range, highlighting differences in the solutions and facilitating intuitive conclusions.
This paper presents an innovative method for accurately measuring and analyzing the temperature distribution of electrical equipment during operation. Addressing the challenge of precisely measuring ohmic losses with traditional methods, this study utilizes the Maxwell eddy current solver for accurate calculations, which then serves as the basis for solving the temperature field. Subsequently, the temperature field-solving module is used to visually display the temperature distribution of each component, significantly enhancing the efficiency and accuracy of engineering tests. The novelty of this paper lies not only in providing an efficient method for solving the temperature field but also in demonstrating how to use these results for structural optimization and heat dissipation analysis. This offers new technical means for the design and optimization of electrical equipment. However, there were limitations in the research process, and in future work, further research will be conducted in the following areas. Firstly, model refinement and multi-field coupling. The model will be further refined to accurately reflect the thermal characteristics of the ring network cabinet during actual operation. Secondly, studying transient temperature fields. In-depth research will be carried out on the temperature field changes of the ring network cabinet under varying load conditions to provide additional support for the stable operation of the ring network cabinet in complex operating conditions.
Disclosures
The authors have nothing to disclose.
Acknowledgements
The authors thank Mr. Wu, MS Sun, Mr. Wang, Mr. Mu, and Mr. Li for their help. This study was supported by the China Postdoctoral Science Foundation (2022M721604) and the Wenzhou Key Science and Technology Tackling Programmer (ZG2023015).
Materials
| Air | / | / | Conventional gases |
| Aluminum | / | / | Alloy Materials |
| Copper | / | / | Alloy Materials |
| Icepak | ANSYS company | ANSYS 2021R1 | A CFD thermal simulation software |
| PC hosting | / | 12th Generation Intel(R) Core(TM) i5-13500F CPU | Host computer equipment |
| SolidWorks | Subsidiary of Dassault Systemes | SolidWorks2021 | An engineering software drawing tool |
| Steady-state thermal | ANSYS company | ANSYS 2021R1 | A thermal simulation solution tool |
References
- Xia, H., et al. Temperature rise test and analysis of high current switchgear in distribution system. J Engg. , 754-757 (2019).
- Polykrati, A. D., Karagiannopoulos, C. G., Bourkas, P. D. Thermal effect on electric power network components under short-circuit currents. Electric Power Syst Res. 72 (3), 261-267 (2004).
- Guan, X., Shu, N., Kang, B., Zou, M. Multiphysics analysis of plug-in connector under steady and short circuit conditions. IEEE Trans Comp Packag Manu Technol. 5 (3), 320-327 (2015).
- Wang, L., Wang, R., Li, X., Jia, S. Simulation analysis on the impact of different filling gases on the temperature rise of C-GIS. IEEE Trans Comp Packag Manu Technol. 9 (10), 2055-2065 (2019).
- Mueller, A., et al. Numerical design and optimization of a novel heatsink using ANSYS steady-state thermal analysis. 2020 27th International Workshop on Electric Drives: MPEI Department of Electric Drives 90th Anniversary (IWED). , 1-5 (2020).
- Wang, Y., Yan, J., Yang, Z., Zhao, Y., Liu, T. Optimizing GIS partial discharge pattern recognition in the ubiquitous power internet of things context: A MiNET deep learning model. Int J Electrical Power Energy Sys. 125, 106484 (2021).
- Lei, J., et al. A 3-D steady-state analysis of thermal behavior in EHV GIS Busbar. J Electr Engg Tech. 11 (3), 781-789 (2016).
- Ouerdani, Y., et al. Temperature rise simulation model of RMU with switchfuse combinations for future load profiles. , 360-364 (2021).
- Zheng, W., Jia, X., Zhou, Z., Yang, J., Wang, Q. Multi-physical field coupling simulation and thermal design of 10 kV-KYN28A high-current switchgear. Thermal Sci Engg Prog. 43, 101954 (2021).
- Wang, L., et al. Electromagnetic-thermal-flow field coupling simulation of 12-kV medium-voltage switchgear. IEEE Trans Comp Packag Manufact Technol. 6 (8), 1208-1220 (2016).
- Zhu, Y., et al. Thermal analysis and design of GaN device of energy storage converter based on Icepak. , 762-767 (2022).
- Mao, Y. e. Thermal simulation of high-current switch cabinet based on Icepak. Electr Ener Mgmt Technol. , 1-7 (2018).
- Zhang, S. Evaluation of thermal transient and overload capability of high-voltage bushings with ATP. IEEE Trans Power Delivery. 24 (3), 1295-1301 (2009).
- Ghahfarokhi, P. S., et al. Steady-state thermal model of a synchronous reluctance motor. , 1-5 (2018).
- Şeker, E. A., Çelik, B., Yildirim, D., Sakaci, E. A., Deniz, A. Temperature field and power loss calculation with coupled simulations for a medium-voltage simplified switchgear. Electrica. 23 (1), 107-120 (2021).
- Ruibo, Y., et al. Research and application of temperature load of switchgear. J Physics: Conf Series. 2378 (2022), 012019 (2022).
- Sheikholeslami, M., Khalili, Z. Simulation for impact of nanofluid spectral splitter on efficiency of concentrated solar photovoltaic thermal system. Sust Cities Soc. 101, 105139 (2024).
- Sheikholeslami, M., Khalili, Z., Scardi, P., Ataollahi, N. Environmental and energy assessment of photovoltaic-thermal system combined with a reflector supported by nanofluid filter and a sustainable thermoelectric generator. J Cleaner Prod. 438, 140659 (2024).
- Sheikholeslami, M., Khalili, Z. Solar photovoltaic-thermal system with novel design of tube containing eco-friendly nanofluid. Renewable Ener. 222, 119862 (2024).
- Sheikholeslami, M., Khalili, Z. Environmental and energy analysis for photovoltaic-thermoelectric solar unit in existence of nanofluid cooling reporting CO2 emission reduction. J Taiwan Inst Chem Eng. 156, 105341 (2024).
- Zhao, L., et al. Research on the temperature rise characteristics of medium-voltage switchgear under different operation conditions. IEEJ Trans Elect Electr Engg. 17 (5), 654-664 (2022).
- Fjeld, E., Rondeel, W., Vaagsaether, K., Attar, E. Influence of heat source location on air temperatures in sealed MV switchgear. , 1-5 (2017).
- Icoz, T., Arik, M. Light weight high performance thermal management with advanced heat sinks and extended surfaces. IEEE Trans Comp Pack Technol. 33 (1), 161-166 (2010).
- Steiner, T. R. High temperature steady-state experiment for computational radiative heat transfer validation using COMSOL and ANSYS. Results Engg. 13, 100354 (2022).
