CC BY
86
UDC 621.314
Ivankov V. F.1, Basova A. V.2
PhD, Chief of bureau of PJSC «ZTR», Ukraine, E-mail: [email protected]
2Leading engineer-designer of PJSC «ZTR», Ukraine
CALCULATION OF CFD-THERMAL MODELS OF OIL-COOLED
TRANSFORMER EQUIPMENT
The purpose ofthe study is to ensure designing of full-function and stable CFD-simulation procedure for integrated thermal models of the transformers and the reactors, as well as to receive the approval of method of quality and abilities using the calculation examples of full-scale models of the equipment, along with autonomous models of coil-type windings having various design versions of heat-transfer intensification.
Research Methods. Computational Fluid Dynamics (CFD) method of mathematical simulation of nonlinear processes as concerns hydrodynamics and heat transfer in the transformer equipment using finite-element analysis is employed.
The results obtained. The paper presents the main elements oftechnique for creation ofmathematical models; it also contains the examples of CFD-calculations as referred to axisymmetrical integrated models of furnace transformer and gapped-core shunt reactor, as well as the models of windings having design approaches of heat exchange intensification owing to «labyrinth» (partitions) and «alternation» (ofnumber and locations) ofaxial cooling ducts.
Scientific novelty. Scientific value ofapplied methodological approach lies in the fact that the developed models are the integrated ones, i.e., they consider geometry, loss, thermal parameters not only of the windings, but also of the main structural elements and cooling system. This ensures the quality and the accuracy of simulation of heat-and-mass transfer processes in complex structure of oil ducts and coils in the windings, enables to avoid erroneous «zigzag» oil flow movement through the groups of coil regular structures (without labyrinth and «alternation» of number and locations of axial ducts under conditions of transformer oil natural cooling as was deduced in the certain studies.
Practical significance. Integrated models ensure calculation of oil temperature distribution within active part, including windingfields, oil temperature field between the tank and the windings, temperatures at oil outletfrom the tank (top) and oil inlet into the tank (bottom). Calculations allow estimation of mean temperature distribution over the cross-section ofwinding coils, mean winding temperatures by means ofaveraging of the temperatures within the coils, detection of location and maximum temperature on the surface of conductors relevant to the most heated coil. The latter is treated as winding hot spot temperature (HST) and used to evaluate the aging of the contacting insulation. Determination of winding hot spot locations and temperatures (HST) is used as support data for installation areas of fiber optic probes for measurement during type testing, as well as in operational monitoring systems of the equipment. The results presented above are practically appliedfor industrial designing and testing of transformers and reactors.
Key words: transformers, heat-and-mass transfer, CFD-simulation.
INTRODUCTION
An important condition for ensuring reliability and duration of operation (load capacity) of the transformer equipment is efficient removal of active energy released as the heat in the magnetic system (MS), the windings, the conductive elements of active part structure, and the tanks. Namely, efficiency of heat removal in the windings determines essentially technical, mass-and-overall and economical characteristics of the equipment. Therefore, within designing availability of reliable calculation methods of winding heating is of great importance [1]. This is especially true in the cases of approaches of heat transfer intensification in the coils due to oil labyrinth flow in the coil groups, provision of axial ducts with alternation of their number and locations along the coil width, etc.
Known methods of oil heating calculations, transformer equipment windings within steady-state temperature mode are conventionally divided in the following groups:
- «empirical method of overheating» [1-5] based on application of heat equation, and empirical heat transfer coefficients averaged over the coils' surfaces, taking into
© Ivankov V. F., Basova A. V., 2016 DOI 10.15588/1607-6761-2016-2-3
account oil temperatures in place of coils arrangement along the winding height and summation of temperature rises of coil conductor surface over the oil, temperature gradient along insulation thickness of elementary conductors and common insulation of the conductors. Empirical methods are widely used in industrial designing [6], however, their application is limited to the range of empirical data by type, size and design of the windings; in many cases these methods do not allow the required accuracy to determine hot-spot locations (HST) of the windings and their temperature;
- methods of nonlinear thermal-hydraulic circuits both within the windings [7-9], and for common thermal-hydraulic circuit of the transformer, including external cooling system (CS) [10]. When used, the complex processes in continuous media as described by equation system (1) - (3), are represented simplistically by zero-dimensional thermal and hydraulic elements having the lumped parameters, that do not allow to determine the local values of temperature and velocity;
- methods of field simulation of thermal, hydrostatic and hydrodynamic winding fields [11-14] based on the
numerical solution of the problems (1) - (3) using various additional simplifications as to their formulation; modern advanced CFD methods (Computational Fluid Dynamics) of mathematical modeling of hydrodynamics and heat transfer processes in the transformers, in particular, using the system of finite element analysis (FEA) [9, 15-23].
A similar group of methods is also emphasized in the summary paper [24], where the first of these methods has been defined as the empirical method of «correlation»; there are also large-scale and detailed modern studies as regards to CFD-modeling, investigation of transfer thermal processes, measurements of windings HST, including the method using fiber optic sensors.
Methodological approaches as to the CFD-modeling of power transformer thermal modes of using the complex models of transformer equipment composed of the MS, the conductive structural elements, the windings and the external cooling systems are presented in [16-21].
At consideration of mentioned works as regards to CFD-modeling, the following notes should be made.
In contrast to [16-21], [9, 15, 22-24], the complex CFD-modes of the transformers are omitted, and there are no studies as regards to electrical reactors.
In [16], one of methodological approaches is used for consideration of development and investigation of the systems of consistent field macro- and micro-models of the devices, the windings, the coil groups. However, experience of practical calculations has shown that this approach has the restricted application. Computation based on micromodels of winding equivalent parameters and their further application in the complex thermal macro-models has been problematic.
At the same time, simulation of the windings, the coil groups by means of consideration of their autonomous models having specified boundary conditions is appropriate [9, 22-24].
In a sense, the most stable thermal-hydraulic processes of oil motion (velocity fields) and the distribution of heating (temperature fields) are demonstrated in all of these studies using the windings with «arranged» direction of cooling oil. Typically, this is achieved by means of so-called labyrinth movement of oil in groups of horizontal ducts along winding height separated by the guiding spacers. In the number of studies [19, 20, 22, 23], it has been turned up that oil velocity vectors can change oil direction in the case of natural convection of oil in the windings with regular structure of horizontal and vertical ducts. So, in the paper [19, 20] «Analysis of flow pattern in the system of interconnected ducts allows detecting of spontaneous zigzag flowing of oil in groups, that directly affects the coil temperature». It should be noted that in the certain experimental studies of such nature of the velocity field is reported. The numerical studies conducted by the authors of this paper have shown that such results can be obtained by faulty settings of computational process.
The above mentioned problem is to be focused on creation of adequate mathematical model of oil flow in the
system of interconnected horizontal and vertical ducts of the winding by means of invariant software. This issue was especially considered in the study [17, 18]. The studies were carried out using conventional winding models with limited coil number (10, 3) placed in the thank with external SW. On assumption as for laminar nature of oil flow it was determined that in horizontal ducts oil flow makes a loop-like movement with gradual increase of oil penetration depth to the coil center as it flows from the bottom upwards along the vertical ducts. Oil flow in the duct prompts the flow during formation of Reyleigh-Bernard convection cells [25]. In [25] the effect of destruction and detachment of oil boundary layer at the area nearby the vertical surfaces of the coils is also noted, which does not mean the beginning of turbulent (eddy) motion of oil, but is caused by mixing of oil particles having different temperatures. In both of mentioned studies [17, 18], there were no results as referred to straight-through oil flow in the horizontal ducts. It is emphasized that numerical studies were carried out using the winding fragments consisting of limited number of coils.
But further in the studies [19, 20], during investigation of complete integrated model of power transformer 140 MVA, it was found in HV winding consisting of 156 coils that «in horizontal and vertical ducts of the winding the steady-state structures in the form of zigzag and/or loop-like oil flow within the groups of horizontal ducts are created.» The resulting numerical result is considered to be erroneous.
The relevance of the study is determined by the need to creation of high-performance transformer equipment with increased specific electromagnetic and thermal loads, reliability of which is achieved also by means of heat transfer intensification in the windings.
The purpose of the study is to ensure designing of full-function and stable CFD-simulation procedure for integrated thermal models of the transformers and the reactors, as well as to receive of the approval of method quality and abilities using the calculation examples of full-scale models of the equipment, as well as autonomous models of coil-type windings having various design versions of heat-transfer intensification.
KEY ELEMENTS OF PROCEDURE
For the purpose of formation of mathematical models and CFD-studies of oil-cooled transformer equipment the following method basic elements (including the papers [1621], are presented as thesis).
1 Thermal processes in oil-cooled electrical devices in the most complete formulation should be considered using known Navier-Stokes equations of motion and continuity of cooling liquid.
p(V -V)v = -V-(PI) + V-X + pg ,
X =
(vv + VvT )-
(VV + WT )--V- VI
(1)
3
V-(pv) = 0, (2)
to be supplemented by energy-conservation equation
v-(-xve)=qv -pC (v -ve), (3)
where the following is designated: e, P and V - fields of temperatures, pressures and velocities of cooling medium,
Qv (e) - apparent densities of heat sources, p(e), |(e),
C(e) and X(e) - oil density, its dynamic viscosity, heat
capacity and thermal conductivity, g - gravity vector, matrix
I determines the direction of unit vector of selected coordinates system.
2 Complex three-dimensional structure of core-type transformer equipment is reduced to complex axisymmetrical models in view of the fact that MS cores, windings, insulating cylinders nearby the windings and structure of yoke insulation have the pronounced cylindrical symmetry.
3 Dimensions of cylindrical tank models are determined based on the number of cooling oil in the tank of the device under consideration. Heat dissipation of the tank into ambient air is modeled by specified factors of heat transfer from the horizontal and vertical surface, taking into account the radiation and the convective heat exchange with air. In cases of multi-phase designs the single-core thermal structure is formed (with proportional release of active losses).
4 The winding models do not consider the discrete structure of vertical rods (or supports in the form of «bridges») and horizontal spacers between the coils.
5 For qualitative (comparative) studies, for example, of possible ways of intensification of heat transfer in the windings, the models of autonomous models of the windings or the coil groups with appropriate boundary conditions are employed.
6 It is advisable to build up the special grid-scale models based on the objectives of investigations and assumptions of laminar or turbulent nature of oil flow.
7 To consider nonlinear properties of transformer oil, the experimental data are used. [25]
8 For complex CFD-models of the transformer equipment, two approaches to consider CS are approbated. The first approach [16-20] lies in the fact that external CS is modeled as the cylindrical element with appropriate volume of oil that determines its dimensions as to design model. The characteristic of heat removal known from engineering calculations is specified for CS element [26]. The points of CS connection to the tank correspond to actual design. The width of the gap of CS design model branch pipes is determined based on condition of equity of total flow passage of the branch pipes in the model gap. As to the second approach [21], the boundary conditions in the form of velocity and temperature of inlet oil are specified. In both approaches, the thermal parameters of CS and inlet oil are determined by the appropriate calculation using the empirical methods [6, 26].
9 At setting of heat sources, the losses in MS, structural elements, windings and tank are assumed to be known (constant), and could be estimated using the complex of recalculation [6] or best estimate calculations using FEA [27].
10 Both at calculation of equipment integrated models and winding autonomous models, one of the key issues is equivalenting of the parameters of complex discrete structures of the winding coils due to the fact that number of coils over winding height may be several hundred, and number of elementary wires in each coil, for example, as referred to CTC, could reach up to hundred wires.
The coils of transformer and reactor windings are wound in radial direction using different types of conductors: conventional conductors (CC), divided (DT) or transposed (CTC) conductors. Elementary conductors can have enamel or paper insulation, DT and CTC conductors have external turn insulation. Two horizontal rows of elementary CTC conductors are divided from each other by paper, cardboard or fiberglass spacers. During numerical studies, the mentioned complex discrete structure is presented as homogeneous bodies with equivalent thermal conductivity in horizontal and vertical directions.
Effective thermal conductivity \e of the material (corresponds to the term «coefficient of thermal conductivity») is estimated as the ratio of thickness L of the material sample to its thermal resistance R and the area in direction of heat flow ££ [28]:
Xeff
L
AS • R
eff
(4)
The laws for estimation of total thermal resistance in case of parallel and series connection of individual thermal resistances similar for linear electric circuits are employed. If heat flow is directed along the layers, Xef should be calculated in the same way as electrical conductivity of the circuit with parallel resistance, in case of flow directed perpendicularly to the layers - according to series circuit.
Resistance of the material based on estimation of effective thermal conductivity:
Reff -
L
AS •X
eff
(5)
For schematic contact between two bodies (Fig. 1), calculation of effective thermal conductivity in the axial and radial directions is analyzed.
In the axial direction z there is series connection of two samples having thicknesses L,, L2 in direction of heat flow, with the same section AB, when thermal conductivity of
each sample X1, X2is:
(
Tfser _
Reff -
A • B
L1 + V X1 X 2
(6)
1
For schematic contact between two bodies (Fig. 1), calculation of effective thermal conductivity in the axial and radial directions is analyzed.
In the axial direction z there is series connection of two samples having thicknesses L,, L2 in direction of heat flow, with the same section AB, when thermal conductivity of
each sample X1, X2is:
nser _
Reff =
A - B
L1 + L2_ ^ 2 y
(6)
Then, the effective thermal conductivity of the sample having total thickness L+L2 and cross-section AB is:
eff
(L1 + L2 2
L{k 2 + L2^1
(7)
Fig. 1 - The model of thermal contact between two bodies
In the radial direction r there is parallel connection of two samples having thickness B in the direction of heat flow, with cross-sections proportional to lengths L, L2:
Rf =
B/A
L1^1 + L2^2
(8)
Then, the effective thermal conductivity of the sample having the thickness B and total cross-sectionA^L+LJ is:
^ par = L1^1 + L2^ 2
eff = Li + L2
(9)
Obtained relationships allow estimation of the parameters of anisotropic thermal conductivity depending on the specific parameters of the coils and used conductors.
As an example, Fig. 2 shows the calculation results of autonomous coils made of CTC, which ensure verification of the present method of equivalent estimation of discrete structure of coil conductors using anisotropic thermal conductivity.
11 The experience of calculations has shown the significance of quality settings of computational processes which should ensure stable solutions of nonlinear systems (1) - (3) and elimination of problem as to «spontaneous zigzag oil flowing through the ducts groups» (without labyrinth and alteration of number and locations of axial cooling ducts) under conditions of natural cooling with transformer oil. [19, 20, 22, 23]. The above mentioned problem is to be focused on creation of adequate mathematical model of oil flow in the system of interconnected horizontal and vertical ducts of the winding by means of invariant software.
In the judgment of the authors the main method of solution of the problem under consideration is application
1
Fig. 2 22
- Distribution of temperature field of DC coil in oil: at the top - with discrete conductor specification at the bottom - with
equivalent anisotropic conductivity
of improved oil flow models if compared with standard models. In particular, use of recommendations [29] related to simulation of heat transfer in the ducts with intensifiers by means of application of k-e, k-w models of turbulence and wall functions separating the purely turbulent motion zone which comply with logarithmic law, and the zone of viscous sub-layers close to wall surface was found to be reasonably effective.
After design experiments as to use of turbulence models, for example, k-e with different wall functions, the desired result was obtained for the scalable wall function (SWF). The main idea of SWF is to restrict relative size of wall viscous layer. This statement has showed impact-decreasing effect of oil rate in the horizontal duct in regard to coil heating. Effect of oil boundary layers detachment nearby the vertical surface of the coils is observed [25]. Presence of insignificant areas of stagnant oil in the ducts does not cause coil overheating; zigzag oil flow along the groups of winding ducts is eliminated. Moreover, distribution of maximum and average temperature rise of the coils over the height becomes almost linear, that is in line with experimental studies [2].
CALCULATION EXAMPLES OF INTEGRATED CFD-MODELS OF TRANSFORMERS AND REACTORS
As an example, Fig. 3 shows the sketch of the design model of one of the cores (phase) of of furnace transformer EOCNK - 45000 / 34.5.
Core losses in the magnetic system, the coils of internal EW winding, the external LV2 winding and in the tank have been estimated using appropriate calculations for the transformer under consideration. The special element which simulates CS heat transfer is included into the model. The positions of oil inlet and outlet from the tank into the CS represent the actual condition; at the area of inlet into the tank an additional oil head is specified to simulate forced oil circulation in the system due to installed pumps.
The geometry of the transformer model and calculation results are shown in Fig. 4.
As the results of calculation during steady-state thermal mode the following is determined: distribution of oil temperature in the active part, including the winding fields, temperature field between the tank and the windings, the temperature at the points of oil outlet from the tank (top) and oil inlet into the tank (bottom); distribution of average temperatures over the section of the winding coils; the average temperatures of the windings as averaging of the temperatures in the coils, maximum temperature on the surface of the conductors of the most heated coil treated as HST during estimation of aging of insulation in contact.
Application of temperature fields of coils and oil in the winding allows compliance with the requirements of IEC 60076-7 standard [30] as for estimation of HST temperature rise above oil temperature in the windings directly, but not in simplified way [1-6] - by calculation or based on measured temperatures in central and top parts of the tank.
Measurement results obtained during testing including those of fiber optic sensors at the area of upper winding coils, correspond to calculation data with appropriate accuracy.
Fig. 3 - Sketch of CFD-model of the transformer: 1 - MS core; 2 - tank, 3 - CS; 4 - common pressing ring made of electric cardboard, 5 - top yoke insulation, 6 - cylinders of the main insulation, 7 - bottom yoke insulation, 8 - lower support segment from of electric
cardboard
Under similar assumptions, CFD-model of the reactor type ROM-80000/765 is considered (system of natural air cooling ONAN) - Fig. 5. In contrast to the CFD-model of the transformer, the reactor model has external CS with boundary conditions at area of oil inlet into the tank.
The model takes into account the core of discrete discs, the tie rod, the vertical cardboard cylinders of the main insulation, the elements of upper and lower yoke insulation, the common pressing ring made of electric cardboard. The reactor winding is represented by four concenters separated by axial ducts with «bridge» -type spacers. The geometry of the reactor and the temperature field model are shown in Fig. 5a.
The mentioned Figure should be provided with further clarifications. The tie rod located inside the core has a complex, but not cylindrical shape, which does not allow
If - —
= = =
a) b)
Fig. 4 - CFD-model of the transformer EOCNK - 45000/34.5: a) - temperature fields, b) - velocity fields
proper heating based on specified losses, so it is omitted in the geometry of design model, and oil duct between the rod and the internal surface of the core corresponds to the actual condition.
Oil velocity field in the tank and the active part of the reactor is shown in Fig. 5b.
The enlarged fragment of the reactor upper is shown in Fig. 6.
Detailed analysis of temperature fields and rates of cooling oil over entire height of EW internal winding, in LV2 external winding of the transformer, as well as in winding of the reactor, including the fragment of its upper part, is the evidence of achieving the quality simulation of heat transfer in both integrated models of the equipment and also as referred to the effectiveness of applied mathematical CFD-simulation techniques. Oil movement in the windings has no numerical effects of cross-through flow of oil.
Type testing of the rector has also confirmed the appropriate compliance of calculation with measurements of oil temperature in the tank, average winding temperature and HST of winding
a)
b)
Fig. 5 - CFD-model of the reactor R0M-80000/765: a) - temperature field, b) - velocity field
CALCULATION EXAMPLES CFD-MODELS OF AUTONOMOUS WINDING
Calculation of temperature rises of common winding (CW) with labyrinth (with the partitions) and of series winding (SW) with «alternation» of number and locations of axial cooling ducts of autotransformer rated for capacity 83.3 MVA is investigated as an examples of CFD-models of autonomous windings. Both methods of heat transfer improvement are well-known, employed for production of power transformers and reactors, and represented in the industrial procedures [1-6] with application of «overheating empirical method.» However, empirical methods do not allow estimation of HST of the coils with required accuracy, and this is the task of CFD-calculations.
CW winding has been simulated as autonomous group of upper coils. At the right and the left (along the coil ends) the model is limited by vertical insulation cylinders, on the top - by pressing rings, at the bottom the section is made at the center of height of lower coil model. Dimensions of the coils, vertical and horizontal cooling ducts correspond to design sizes; the geometry of the coils has simplified (rectangular) shape; filling the coil with the conductors is considered in specified thermal properties of the coils. Upper top coils generate the upper passage of labyrinth oil movement. Labyrinth oil flow is formed by means of horizontal partitions arranged in external vertical cooling duct above coil 1 and below coil 8.
Oil temperature at inlet to the bottom planes of vertical cooling ducts of the model is determined based on method [26]. Oil velocity in the lower ducts of the model is assumed from known assumption that oil consumption is proportional to the winding losses.
Calculation results are shown in Fig. 7. It has been determined that the maximum HST was in the fourth top
a)
Fig. 6 - Fragment of CFD-model of the reactor ROM-!
coil. Its absolute temperature makes up 94,4°C, which is appropriate at ambient temperature 20°C due to permissible value 78 K of HST temperature rise above ambient temperature.
Obtained result is explained as follows. To top three coils of upper path of the labyrinth the higher horizontal ducts are adjacent than those adjacent to the coils 4-8. Increased oil rates in these ducts are characterized by reduced hydraulic resistance of oil, as well as by forced acceleration of oil rate due to oil labyrinth flow - Fig. 7b. The factors mentioned, as well as uneven losses in the coils caused by eddy currents in the coil conductors due to axial and radial components of magnetic induction have formed the conditions under which the coil 4 is the most heated one. Labyrinth oil flow has an effect both on both the uneven temperature field of oil, conductors of the coils - Fig. 7a, and also on oil rate field in vertical and horizontal ducts - Fig. 7b.
SW winding with alteration of number and locations of axial cooling ducts was studied using the same approach -Fig. 8. The coil pairs with small horizontal duct between them are arranged along the winding height. In the top coil of each pair there are two internal axial ducts; in the bottom coil - one duct. Such arrangement of axial ducts provide the conditions similar to labyrinth oil flow. This causes increase of oil rate both in said internal axial ducts and in adjacent areas of horizontal ducts - Fig. 8b, due to unevenness of additional losses in the coils, both over coil height (maximum losses occur in the conductors of the first coil) and in different groups of the conductors along radial dimension of each coil. These factors create the conditions of uneven temperature field in the coils and within the surrounding oil
- Fig. 8a. These factors create the conditions of uneven temperature field in the coils and within the surrounding oil
- Fig. 8a. HST of upper coils makes up 94,7°C.
b)
: a) temperature field, °C; b) oil velocity field, m/s
T, °C
R
I
§p
93 1 91 3 9®
sg:i
i'8'7 G §P
85.2
Si.B Ik eg "m®
77:5 75.9 74.6
73.3 72.0 70 6 ■fe.3 63.0
a)
c -
—«■a*
L m «mm*
1 'J
I
V, m/s
0.068 0.0S5 0 0B1 0.058 0 055 0 051 0 048 0.044 0.041 0 038 0.034 0.031 0.027 0 024 0 020 0.017
0 014
1 lli'l 0.007 0 003 0 000
b)
Fig. 7 - Model of upper part of CW winding with «labyrinth»: a) temperature field, °C; b) oil velocity field, m/s
T, °C
V, m/s
a) b)
Fig. 8 - SW having «alternate arrangement» of vertical ducts: a) temperature field, °C; b) oil velocity field, m/s
Medium voltage windings (MVW) having nine coils of upper passage of «labyrinth», and low voltage windings (LV) with additional vertical duct at upper four coils of the transformer also rated for 83.3 MVA have been estimated. The results are shown in Fig. 9 and 10.
It has been determined that in MV winding with «labyrinth», the most heated coils are the first, the third and the fourth (from the top) having the same HST temperature 93,7°C.
Upper coils in LV winding have rather low heating, it makes up 94,9°C in the sixth coil, and 96,7°C in the seventh coil.
Labyrinth oil flow in MV winding and presence of additional axial ducts at upper four coils of LV winding create uneven temperature and rate (velocity) fields similarly to winding structures shown in Fig. 7, 8, and this determines specific location of windings' HST and their temperature values. It should be noted that it would be problematic to determine with required accuracy the locations and HST temperatures using presented calculation examples of CFD-models of autonomous windings by means of empirical procedures [1-6].
The obtained results have confirmed the reduction of HST heating and average winding temperature rises at the expense of adopted methods of heat transfer improvement,
and are used to justify the locations of installation of fiber optic sensors for the measurement both during transformer type testing and within operation.
CONCLUSION
Complex CFD-models of oil-cooled transformers equipment using magnetic systems, main structural elements, windings, tanks and external cooling systems and external cooling systems have been developed, the results of which are employed in practice of industrial designing and testing of the equipment.
The study ensures quality simulation of heat and mass transfer in integrated models of the transformer equipment that is the evidence of the effectiveness of applied mathematical methods as regard to CFD-simulation.
For the conditions of natural oil flow, recommendation is given as to application of turbulence models, for example, k-e with scalable wall function (SWF). This statement has showed impact-decreasing effect of oil rate in horizontal duct in regard to coil heating. Effect of oil boundary layers detachment nearby the vertical surface of the coils is observed. Presence of insignificant areas of stagnant oil in the ducts does not cause coil overheating; zigzag oil flow along the groups of winding ducts is eliminated.
E
Fig. 9 - Model of MV winding with «labyrinth»: a) temperature field, °C; b) oil velocity field, m/s
Calculation examples using autonomous CFD-models of windings with different method of heat transfer improvement by means of specific design solutions: «labyrinth» oil flow, alternate arrangement and introduction of additional vertical ducts along the coil width, represents the possibility of improved estimation of hot spot temperatures of the windings, and validation of the points to be used for measurement of local temperatures using fiber optic sensors.
REFERENCES
1 Бики М. А. Проектирование силовых трансформаторов. Расчеты основных параметров / М. А. Бики.
- М. : Знак, 2013. - 612 с.
2 Теплоотдача внутренних катушечных обмоток трансформаторов / [А. Б. Швидлер, Ю. А. Михайловский, Г. Б. Чередниченко, Л. А. Клименко] // Электротехника. - 1980. - N 7. - C. 19-21.
3 Методы расчета превышений температуры обмоток силовых трансформаторов / [И. Д. Воеводин, Ю. А. Михайловский, В. М. Черноготский, и др.] // Трансформаторы : Перенапряж. и координация изоляции. Пер. докл. междунар. конф. по больш. электр. системам СИГРЭ-84. - М. : Энергоатомиздат, 1986.
- С. 190-198.
4 Интенсификация теплообмена в катушечных обмотках трансформаторов / [А. Б. Швидлер , Ю. А. Михайловский, Г. Б. Чередниченко и др.] // Электротехника. - 1987. - № 5. - С. 8-10.
a)
Fig. 10 - Model of LV winding with additional upper duct at the
5 РД16 303-86. Руководящий документ. Трансформаторы силовые масляные. Тепловой расчет обмоток с различными способами интенсификации теплообмена / [Разработчики : О.И. Сисуненко, Ю. С. Ги-тин, Ю. А. Михайловский, А. Б. Швидлер, Е. А. Сереб-рицкая]. - Введен с 01.01.87. Запорожье : ВИТ. - 42 с.
6 1ванков В. Ф. Розрахункова шдсистема автоматизо -ваного проектування трансформаторiв i реакторiв /
B. Ф. 1ванков , Ю. Н. Шафiр // Пращ 1н-ту електроди-нашки НАН Украши. Збiрник наукових праць. - К. : 1н-т електродинашки НАН Украши., 2008. - N 18. -
C. 123-131.
7 Oliver A. J. Estimation of transformer winding temperatures and coolant flows using a general network method / A. J. Oliver // IEE Proc., Vol. 127, Pt. C, No. 6, November 1980. - P. 395-405.
8 Transformer design principles: with application to core-form transformers / [R. M. Del Vecchio, B. Poulin, P. T. Feghali and etc] // Gordon and Breach Science Publishers, 2001. - 599 p.
9 Codd J. Assessment of a Hydraulic Network Model for Zig-zag Cooled Power Transformer Windings / Joris Codd,, Wim Van der Veken, Martine Baelman // Preprint submitted to Applied Thermal Engineering January 21, 2015. - Access mode : http://www.sciencedirect.com/ science/article/pii/S 1359431115000885.
10 Pat. US 6,609,079 B1, G01K 17/00, 702/13. Method and arrangement for ascertaining state variables / Seitlinger
b)
four coils: a) temperature field, °C; b) oil velocity field, m/s
28
ISSN 1607-6761. Електротехшка та елекгроенергетика. 2016. N° 2
W. ; Assignee: VA TECH Elin Transformatoren GmbH. -№ US006609079B1; filed 06.05.1999 ; pub. 25.11.1999, № WO99/60682. Access mode : https:// patentimages.storage.googleapis.com.
11 Mufuta J.-M. Modelling of the mixed convection in the winding of a disc-type power transformer / J. -M. Mufuta, E. Bulck // Applied Thermal Engineering 20 (2000). -P. 417-437. - Access mode : http://www.elsevier.com.
12 Кучерявая И. Н. Численное исследование влияния различных факторов процесса теплообмена на эффективность охлаждения трансформаторно-выпря-мительных модулей / И. Н. Кучерявая // Техн. елект-родинашка. - 2000. - № 3. - С. 56-61.
13 Козлов В. В. Щдвищення ефекгивносп примусового охолодження високовольтних масляних трансформа-торiв: автореф. дис. на здобуття наук. ступеня канд. техн. наук: спец. 05.09.01 «Електричт машини й апа-рати» / В. В. Козлов; Нац. техн. ун-т «Харк. полггехн. ш-т». - Харкiв, 2004. - 19 с.
14 Рассальський О. М. Аналiз i розрахунок теплових режимiв силового електроустаткування [Текст] : навч. поаб. для студ. вищ. навч. закл. / О. М. Рассальський ; ЗНТУ - Запорiжжя : [б.в.], 2005. - 146 с
15 Thermal performance of power transformers: thermal calculation tools focused on new operating requirements / Tanguy A, Patelli J. P., Taisne J. P. [and etc] // CIGRE, Paris Session 2004, Report A2-105. -Access mode : http://www.cigre.org.
16 Методические подходы к CFD-моделированию тепловых режимов силовых масляных трансформаторов / П. Г. Круковский, В. А. Яцевский, Л. Н. Конто-рович, В. Ф. Иванков, Д. Д. Юрченко // Промышленная теплотехника. - 2008. - Т. 30. - № 6. - С. 57 - 66.
17 CFD-моделирование теплообмена в катушечных обмотках трансформаторов при естественной конвекции охлаждающего масла / П. Г. Круковский, В. А. Яцевский, Л. Н. Конторович, В. Ф. Иванков, Д. Д. Юрченко // Промышленная теплотехника. - 2009. -Т. 31. - № 4. - С. 17-26.
18 CFD-анализ теплового состояния обмоток масляных трансформаторов при критических размерах горизонтальных охлаждающих каналов / [П. Г. Кру-ковский, В. А. Яцевский, Л. Н. Конторович и др.] // Промышленная теплотехника. - 2009. - Т. 31, № 6. -С. 35-45.
19 Яцевский В. А. Особенности гидродинамики и теплообмена при течении масла во взаимосвязанных каналах катушечных обмоток силовых трансформа-
торов / В. А. Яцевский, П. Г. Круковский // Промышленная теплотехника. - 2011. - №2. - С. 24-34.
20 Yatsevsky V. A. Hydrodynamics and heat transfer in cooling channels of oil-filled power transformers with multicoil windings / V A. Yatsevsky // Applied Thermal Engineering 63 (1) (2014) 347 - 353. - Access mode : http://www.sciencedirect.com/science/article/pii/ S1359431113007746.
21 .Ivankov V F. Verification methods for electrothermal calculations of electric reactors without steel / V F. Ivankov, A. V Basova, N. V Shulga // Електротех-тка та електроенергетика. - 2015. - № 5. - С. 26-34.
22 JiaoY. CFD Study On The Thermal Performance of Transformer Disc Windings Without Oil Guides / Yuhe Jiao // M. SС Thesis. EGI 2012: 089MSC EKV915. Stockholm. - Access mode : http://www.diva-portal.se/ smash/get/diva2:555330/FULLTEXT01.pdf.
23 Wittmaack R. Thermal Design of Power Transformers via CFD / Ralf Wittmaack // Journal of Energy and Power Engineering 9 (2015) 102-107 doi: 10.17265/1934-8975/ 2015.01.012. - Access mode : http:// www.davidpublisher.org/Public/uploads/Contribute/ 550682c1907a5.pdf.
24 Transformer Thermal Modeling. Working Group A2.38, CIGRE. June 2016, ISBN: 978-2-85873-362-0. - 188 p. -Access mode : http://www.cigre.org.
25 Киш Л. Нагрев и охлаждение трансформаторов / Л. Киш ; пер. с венгерского М. А. Бики ; [под ред. Г. Е. Тарле]. - М. : Энергия, 1980. - 208 с., ил. - (Трансформаторы; Вып. 36).
26 Басова А. В. Инженерный расчет комбинированной масляной системы охлаждения трансформаторов / А. В.Басова, В. Ф. Иванков // Пром. теплотехника. -2013. - Т. 35. - № 4. - С. 9-17.
27 Иванков В. Ф. Электротепловые расчетные модели элементов конструкции трансформаторного оборудования / В. Ф. Иванков, А. В. Басова, Н. В. Шульга // Електротехшка та електроенергетика. - 2014. -№ 2. - С. 41-53.
28 Михеев М. А. Основы теплопередачи / М. А. Михеев, И. М. Михеева. - Изд. 2-е, стереотип. - М. : Энергия, 1977. - 344 с.
29 Юрченко Д. / Моделирование теплообмена в каналах с интенсификаторами / Д. Юрченко // ANSYS Advantage. - 2009. - № 10. - C. 32-34. Режим доступа : http://www.ansyssolutions.ru.
30 IEC 60076-7. Power transformers - Part 7: Loading guide for oil-immersed transformers. - 2005. - 114 p. : www.iec.ch.(Eng.)
1ванков В.Ф.1, Басова А.В.2
'Канд. техн. наук, начальник бюро ПАТ «Зап^жтрансформатор», Украша 2Провiдний шженер-конструктор ПАТ «Зап^жтрансформатор», Украша
РОЗРАХУНОК ТЕПЛОВИХ CFD-МОДЕЛЕЙ ТРАНСФОРМАТОРНОГО УСТАТКУВАННЯ З МАСЛЯНИМ ОХОЛОДЖЕННЯМ
Метароботи полягае в забезпечент розрахункового проектування повнофункцюнальною i стткою в об-числювальному процеЫ CFD-методикою моделювання комплексных теплових моделей трансформаторiв iреак-торiв, а також в апробацп и якостi i можливостей на прикладахрозрахунюв натурных зразтвустаткування,
а також автономних моделей котушкових обмоток з pi3HUMU конструктивними способами ттенсифжацп теплообмту.
Методи дослгджень. Використаний CFD-метод (Computational Fluid Dynamics) математичного моделю-вання нелiнiйних процеЫв гiдродинамiки i теплообмту в трансформаторному обладнанн з застосуванням систем сктчено-елементного аналiзу.
Отримат результаты. Представлен основш елементи методики формування математичних моделей i приклади CFD-розрахункiв осесиметричних комплекснихмоделей перетворювального трансформатора i шун-тувального реактора з промiжками в стрижн магнтног системи, а також моделей обмоток з конструктивними способами ттенсифжацп теплообмту за рахунок «лабiринту» (перегородок) та «почережностi» (числа та мкць) осьових охолодних каналiв.
Наукова новизна. Наукова цттсть використаного методологiчного тдходу полягае в тому, щорозроблен моделi е комплексними, тобто враховують геометрiю, втрати, тепловi параметри не тыьки обмоток, а й основних елементiв конструкци i системи охолодження. Науковою новизною роботи е досягнення авторами певним яюсним налаштування обчислювального процесу сталого ршення засобами iнварiантноi системи чи-сельного моделювання нелтшних рiвнянь Нав 'е-Стокса. Це забезпечило яюсть i точтсть моделювання процеЫв тепломасопереносу в складнш структурi масляних каналiв i котушок в обмотках, дозволило уникнути отриманого в деяких до^дженнях помилкового «зиГзаГоподiбного» руху масла по групахрегулярних структур котушок (без лабiринту i без «почережностi» числа та мкць осьових каналiв) за умов природного охолодження трансформаторним маслом.
Практична значимгсть. Комплексн моделi забезпечують розрахунок розподыу температур масла в активны частини, включаючи областi обмоток, поле температур масла мiж баком i обмотками, температури в мiсцях виходу масла з бака (верх) i входу в бак (низ). Розрахунки дозволяють визначити розподы середтх температур по перетину котушок обмоток, середн температури обмоток шляхом усереднення температур в котушках, визначення мiсця i максимальну температуру на поверхн провiдникiв найбыьш нагртог котушки. Остант, яю тлумачаться як температури найбыьш нагртог точки (ННТ) обмоток, використовуються при оцтщ стартня доторкноi iзоляцií. Визначення мiсць i температур ННТ обмоток використовуються для обгрунтування мкцяустановки оптоволоконних датчиюв для вимiрiв при типових випробуваннях i в експлуата-цшних системах монторингу обладнання. Представлен результати застосован в практиц промислового проектування i випробування трансформаторiв i реакторiв.
Ключовi слова: трансформатори, тепломасоперенос, CFD-моделювання.
Иванков В. Ф.1, Басова А. В.2
1Канд. техн. наук, начальник бюро ПАО «Запорожтрансформатор», Украина
2Ведущий инженер-конструктор ПАО «Запорожтрансформатор», Украина
РАСЧЕТ ТЕПЛОВЫХ CFD-МОДЕЛЕЙ ТРАНСФОРМАТОРНОГО ОБОРУДОВАНИЯ С МАСЛЯНЫМ ОХЛАЖДЕНИЕМ
Цель работы заключается в обеспечении расчетного проектирования полнофункциональной и устойчивой в вычислительном процессе CFD-методикой моделирования комплексных тепловых моделей трансформаторов и реакторов, а также в апробации ее качества и возможностей на примерах расчетов натурных образцов оборудования, а также автономных моделей катушечных обмоток с различными конструктивными способами интенсификации теплообмена.
Методы исследований. Использован CFD-метод (Computational Fluid Dynamics) математического моделирования нелинейных процессов гидродинамики и теплообмена в трансформаторном оборудовании с применением систем конечно-элементного анализа.
Полученные результаты. Представлены основные элементы методики формирования математических моделей и примеры CFD-расчетов осесимметричных комплексных моделей печного трансформатора и шунтирующего реактора с зазорами в стержне магнитной системы, а также моделей обмоток с конструктивными способами интенсификации теплообмена за счет «лабиринта» (перегородок) и с «чередованием» (числа и мест) осевых охлаждающей каналов.
Научная новизна. Научная ценность использованного методологического подхода заключается в том, что разработанные модели являются комплексными, то есть учитывают геометрию, потери, тепловые параметры не только обмоток, но и основных элементов конструкции и системы охлаждения. Научной новизной работы является достижение авторами определенными качественными настройками вычислительного процесса устойчивого решения средствами инвариантной системы численного моделирования нелинейных уравнений Навье-Стокса. Это обеспечило качество и точность моделирования процессов тепломассопереноса в сложной структуре масляных каналов и катушек в обмотках, позволило избежать полученного в некоторых исследованиях ложного «зигзагообразного» движения масла по группам регулярных структур катушек (без лабиринта и без «чередования» числа и мест осевых каналов) при условиях естественного охлаждения трансформаторным маслом.
Практическая значимость. Комплексные модели обеспечивают расчет распределения температур масла в активной части, включая области обмоток, поле температур масла между баком и обмотками, температуры в местах выхода масла из бака (верх) и входа в бак (низ). Расчеты позволяют определить распределение
30
ISSN 1607-6761. Електротехшка та електроенергетика. 2016. № 2
средних температур по сечению катушек обмоток, средние температуры обмоток путем усреднения температур в катушках, определение места и максимальной температуру на поверхности проводников наиболее нагретой катушки. Последние, трактуемые как температуры наиболее нагретой точки (ННТ) обмоток, используются при оценке старения соприкасающейся изоляции. Определение мест и температур ННТ обмоток используются для обоснования места установки оптоволоконных датчиков для измерений при типовых испытаниях и в эксплуатационных системах мониторинга оборудования. Представленные результаты применены в практике промышленного проектирования и испытания трансформаторов и реакторов.
Ключевые слова: трансформаторы, тепломассоперенос, CFD-моделирование.
REFERENCES
1 Biki M.A. Proektirovanie silovykh transformatorov. Raschety osnovnykh parametrov. Moscow, Znak, 2013, 612 p.
2 Shvidler A. B., Mikhailovskii Yu. A., Cherednichenko H. B., Klimenko L. A. Teplootdacha vnutrennikh katushechnykh obmotok transformatorov, Elektrotekhnika. 1980, No 7. pp. 19-21.
3 Voevodin I. D., Mikhailovskii Yu. A., Chernohotskii V M., Shvidler A. B. Metody rascheta prevyshenii temperatury obmotok silovykh transformatorov, Transformatory, Perenapriazh. i koordinatsiia izoliatsii, Per. dokl. mezhdunar. konf. po bolsh. elektr. sistemam SIHRE-84, Moscow, Enerhoatomizdat, 1986, S. 190-198.
4 Shvidler A. B., Mikhailovskii Yu. A., Cherednichenko H. B. Intensifikatsiia teploobmena v katushechnykh obmotkakh transformatorov, Elektrotekhnika, 1987, No 5, pp. 8-10.
5 Sisunenko O. I., Hitin Yu. S., Mikhailovskii Yu. A., Shvidler A. B., Serebritskaia E. A. RD16 303-86. Rukovodiashchii dokument. Transformatory silovye maslianye. Teplovoi raschet obmotok s razlichnymi sposobami intensifikatsii teploobmena, Vveden s 01.01.87. Zaporozhe, VIT, 42 s.
6 Ivankov V F., Shafir Yu. N. Rozrakhunkova pidsistema avtomatizovanoho proektuvannia transformatoriv i reaktoriv. Pratsi In-tu elektrodinamiki NAN Ukraini. Zbirnik naukovikh prats, Kiev, In-t elektrodinamiki NAN Ukraini, 2008, No 18, S. 123-131.
7 Oliver A. J. Estimation of transformer winding temperatures and coolant flows using a general network method. IEE Proc., Vol. 127, No. 6, November 1980, P. 395-405.
8 Del Vecchio R. M., Poulin B., Feghali P. T. Transformer design principles: with application to core-form transformers, Gordon and Breach Science Publishers, 2001, 599 p.
9 Codd J. Assessment of a Hydraulic Network Model for Zig-zag Cooled Power Transformer Windings / Joris Codd Wim Van der Veken, Martine Baelman, Preprint submitted to Applied Thermal Engineering January 21, 2015. - Access mode : http://www.sciencedirect.com/ science/article/pii/S 1359431115000885.
10 Pat. US 6,609,079 B1, G01K 17/00, 702/13. Method and arrangement for ascertaining state variables, Seitlinger W. ; Assignee: VA TECH Elin Transformatoren GmbH, No US006609079B1 ; filed 06.05.1999 ; pub. 25.11.1999, № W099/60682. Access mode : https:// patentimages.storage.googleapis.com.
11 Mufuta J.-M., Bulck E. Modelling of the mixed convection in the winding of a disc-type power transformer, Applied Thermal Engineering 20 (2000), P. 417-437. Access mode : http://www.elsevier.com.
12 Kucheriavaia I. N. Chislennoe issledovanie vliianiia razlichnykh faktorov protsessa teploobmena na effektivnost okhlazhdeniia transformatorno-vypriamitelnykh modulei. Tekhn. Elektrodinamika, 2000, No 3, pp. 56-61.
13 Kozlov V V Pidvishchennia efektivnosti primusovoho okholodzhennia visokovoltnikh maslianikh transformatoriv: avtoref. dis. na zdobuttia nauk. stupenia kand. tekhn. nauk: spets. 05.09.01 «Elektrichni mashini i aparati», Nats. tekhn. un-t «Khark. politekhn. in-t», Kharkiv, 2004, 19 s.
14 Rassalskii O. M. Analiz i rozrakhunok teplovikh rezhimiv silovoho elektroustatkuvannia, Zaporizkii natsionalnii tekhnichnii un-t, Zaporizhzhia, 2005, 146 s.
15 Tanguy A, Patelli J. P., Taisne J. P. Thermal performance of power transformers: thermal calculation tools focused on new operating requirements, CIGRE, Paris Session, 2004, Report A2-105
16 Krukovskii P. H., Yatsevskii V. A., Kontorovich L. N., Ivankov V F., Yurchenko D. D. Metodicheskie podkhody k CFD-modelirovaniiu teplovykh rezhimov silovykh maslianykh transformatorov. Promyshlennaia teplotekhnika 2008, Vol. 30, No 6, S. 57-66.
17 Krukovskii P. H., Yatsevskii V. A., Kontorovich L. N., Ivankov V. F., Yurchenko D. D. CFD-modelirovanie teploobmena v katushechnykh obmotkakh transformatorov pri este-stvennoi konvektsii okhlazhdaiushcheho masla. Promyshlennaia teplotekhnika, 2009, Vol. 31, No 4, pp. 17-26.
18 Krukovskii P. H., Yatsevskii V. A., Kontorovich L. N., Ivankov V F., Khutornyi V M. CFD-analiz teplovoho sostoianiia obmotok maslianykh transformatorov pri kritiche-skikh razmerakh horizontalnykh
okhlazhdaiushchikh kanalov. Promyshlennaia teplotekhnika, 2009, Vol. 31, No 6, pp. 35 - 45.
19 Yatsevskii V. A., Krukovskii P. H. Osobennosti hidrodinamiki i teploobmena pri techenii masla vo vzaimosviazannykh kanalakh katushechnykh obmotok silovykh transformatorov. Promyshlennaia teplotekhnika, 2011, No 2, pp. 24-34.
20 Yatsevsky V. A. Hydrodynamics and heat transfer in cooling channels of oil-filled power transformers with multicoil windings. Applied Thermal Engineering, 2014, pp. 347-353.
21 Ivankov V. F., Basova A. V, Shulga N. V. Verification methods for electrothermal calculations of electric reactors without steel. Elektrotekhnika ta elektroenerhetika, 2015, No 5, pp. 26-34.
22 JiaoY. CFD Study On The Thermal Performance of Transformer Disc Windings Without Oil Guides / Yuhe Jiao // M. SC Thesis. EGI 2012: 089MSC EKV915. Stockholm.
23 Wittmaack R. Thermal Design of Power Transformers via CFD, Journal of Energy and Power Engineering, 2015, No 9, pp. 102-107.
24 Transformer Thermal Modeling. Working Group A2.38, CIGRE. June 2016, ISBN: 978-2-85873-362-0, 188 p.
25 Kish L. Nahrev i okhlazhdenie transformatorov. Moscow, Enerhiia, 1980, 208 s.,
26 Basova A. V., Ivankov V. F. Inzhenernyi raschet kombinirovannoi maslianoi sistemy okhlazhdeniia transformatorov. Prom. teplotekhnika, 2013, Vol. 35, No 4, pp. 9-17.
27 Ivankov V F., Basova A. V, Shulha N. V Elektroteplovye raschetnye modeli elementov konstruktsii transformatornoho oborudovaniia. Elektrotekhnika ta elektroenerhetika, 2014, No .№ 2, pp. 41-53.
28 Mikheev M. A., Mikheeva I. M. Osnovy teploperedachi. Moscow, Enerhiia, 1977, 344 s.
29 Yurchenko D. Modelirovanie teploobmena v kanalakh s intensifikatorami. ANSYS Advantage, 2009, No 10, pp. 32-34.
30 IEC 60076-7. Power transformers Part 7: Loading guide for oil-immersed transformers, 2005, 114 p.
