CN115169267A - Numerical simulation method for overflow-free icing growth on surface of power transmission line insulator - Google Patents

Numerical simulation method for overflow-free icing growth on surface of power transmission line insulator Download PDF

Info

Publication number
CN115169267A
CN115169267A CN202211098383.1A CN202211098383A CN115169267A CN 115169267 A CN115169267 A CN 115169267A CN 202211098383 A CN202211098383 A CN 202211098383A CN 115169267 A CN115169267 A CN 115169267A
Authority
CN
China
Prior art keywords
insulator
ice
icing
transmission line
growth
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211098383.1A
Other languages
Chinese (zh)
Other versions
CN115169267B (en
Inventor
胡玉耀
刘宗源
赵冲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN202211098383.1A priority Critical patent/CN115169267B/en
Publication of CN115169267A publication Critical patent/CN115169267A/en
Application granted granted Critical
Publication of CN115169267B publication Critical patent/CN115169267B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Geometry (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Databases & Information Systems (AREA)
  • Computer Graphics (AREA)
  • Operations Research (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Insulators (AREA)
  • Wind Motors (AREA)

Abstract

A numerical simulation method for the increase of ice coating on the surface of an insulator of a power transmission line without overflow belongs to the technical field of prediction of ice coating on the insulator of the power transmission line. S1, establishing an equal-proportion insulator three-dimensional geometric model and constructing an external flow field area; s2, determining the motion tracks of the insulator icing continuous phase air flow field and the dispersed phase water drop; s3, obtaining a local collision coefficient of the surface of the insulator by a triangular area projection method; s4, judging the icing growth type according to the freezing coefficient value; s5, reconstructing an icing boundary by a point-line-surface-body modeling method to obtain a three-dimensional model of the insulator icing growth form; and S6, taking the three-dimensional model as an initial condition of an air flow field in the next time step, and repeating the steps to perform the circulation iteration of ice coating growth until the ice coating form in the required time is obtained. The method can predict the icing growth form, the icing growth thickness at any position on the surface of the insulator and the icing quality, and can be used for constructing an icing early-warning mechanism of the power transmission line in an extreme environment.

Description

Method for simulating numerical value of ice coating growth on surface of insulator of power transmission line without overflow
Technical Field
A numerical simulation method for the increase of ice coating on the surface of an insulator of a power transmission line without overflow belongs to the technical field of prediction of ice coating on the insulator of the power transmission line.
Background
With the rapid development of national economy, in order to meet the increasing power demand and realize the large-scale optimized configuration of power resources in China, a national power grid company builds a plurality of extra-high voltage alternating current and direct current transmission lines. The ultra-high voltage transmission line has obvious advantages in the aspects of transmission capacity, electric energy loss, economy and the like. However, as the transmission distance is long, the insulator inevitably passes through regions with complex climatic environments, and the electrical characteristics of the insulator as an important device in an electric power system play a decisive role in the safe and stable operation of the electric power system. Under normal environment, the performance of the insulator can meet the requirement of line operation. However, in low-temperature weather, the electrical and mechanical strength of the insulator is seriously reduced by the ice and snow on the surface of the insulator, so that accidents such as insulator flashover, pole falling and tower falling seriously affect the safe and stable operation of a power grid occur, and therefore, the design standard and the disaster prevention level of the power transmission line insulator need to be improved aiming at the problems.
At present, most researches on an icing prediction model of a power transmission line are concentrated on the aspect of power transmission conductors, and due to the fact that an insulator is complex in structure, a mature insulator icing prediction model is not available at present.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method overcomes the defects of the prior art, and provides a numerical simulation method for the ice-coating growth of the surface of the insulator of the power transmission line without overflow, which is used for establishing an anti-icing and ice-reducing early warning mechanism of the power transmission line in an extreme environment.
The technical scheme adopted by the invention for solving the technical problems is as follows: the numerical simulation method for the ice-coating growth without overflow on the surface of the insulator of the power transmission line is characterized by comprising the following steps of: the method comprises the following steps:
s1, establishing an equal-proportion insulator three-dimensional geometric model and constructing an outer flow field area;
s2, setting boundary conditions of a flow field according to actual meteorological data to obtain a continuous phase air flow field of the ice-coated insulator and a motion trail of dispersed phase water drops;
s3, extracting three-dimensional coordinates, speed and direction of the water drops in the whole process from releasing to being captured by the insulator, and obtaining a local collision coefficient of the surface of the insulator through a triangular area projection method;
s4, constructing a thermal equilibrium equation according to an energy conservation principle to determine freezing coefficients of all positions on the surface of the insulator, and judging the icing growth type according to the freezing coefficient values;
s5, determining the icing growth within a time step, and reconstructing an icing boundary by a point-line-surface-body modeling method to obtain a three-dimensional model of the icing growth form of the insulator;
and S6, taking the three-dimensional model as an initial condition of the air flow field in the next time step, and repeating the steps to perform the cycle iteration of ice coating growth until the ice coating form in the required time is obtained.
Preferably, the boundary conditions in S2 are:
Figure 261537DEST_PATH_IMAGE001
Figure 352990DEST_PATH_IMAGE002
Figure 593478DEST_PATH_IMAGE003
；
Figure 137592DEST_PATH_IMAGE004
；
wherein,
Figure 837343DEST_PATH_IMAGE005
for the inlet velocity of the turbulent flow field,
Figure 658668DEST_PATH_IMAGE006
is the ambient wind speed and is,
Figure 194692DEST_PATH_IMAGE007
in order to obtain the outlet pressure of the turbulent flow field,
Figure 101468DEST_PATH_IMAGE008
in order to be the intensity of the turbulent flow,
Figure 381140DEST_PATH_IMAGE009
in order to be of a turbulent flow scale,
Figure 56971DEST_PATH_IMAGE010
in order to calculate the hydraulic diameter of the domain,
Figure 763896DEST_PATH_IMAGE011
is the reynolds number.
Preferably, the method further comprises the following steps of regarding the air external flow field of the ice-coated insulator as normal-temperature, low-speed and incompressible turbulent flow, and establishing a turbulent flow model:
Figure 157969DEST_PATH_IMAGE012
；
wherein,
Figure 975752DEST_PATH_IMAGE013
is the velocity vector in the flow field;
Figure 771670DEST_PATH_IMAGE014
is the kinetic viscosity of air;
Figure 383916DEST_PATH_IMAGE015
due to the dynamic viscosity additionally generated by air turbulence;
Figure 530864DEST_PATH_IMAGE016
is the air density;
Figure 152338DEST_PATH_IMAGE017
as kinetic energy of turbulent flow；
Figure 802762DEST_PATH_IMAGE018
Is the turbulent dissipation ratio;
Figure 851490DEST_PATH_IMAGE019
is the principal stress tensor;
Figure 954575DEST_PATH_IMAGE020
is a volume force;
Figure 117091DEST_PATH_IMAGE021
、
Figure 887601DEST_PATH_IMAGE022
are respectively as
Figure 107229DEST_PATH_IMAGE017
And
Figure 697611DEST_PATH_IMAGE018
the inverse of the effective prandtl constant of (a);
Figure 660887DEST_PATH_IMAGE023
、
Figure 285904DEST_PATH_IMAGE024
、
Figure 676434DEST_PATH_IMAGE025
is a turbulence model parameter;
Figure 754111DEST_PATH_IMAGE026
are turbulent flow energy terms.
Preferably, the method further comprises the step of ignoring smaller Saffman lifting force, additional mass force, pressure difference force and the like, and considering that the supercooled water drops are only acted by airflow drag force and gravity, wherein the Lagrange motion control equation of a single water drop is as follows:
Figure 255500DEST_PATH_IMAGE027
wherein the first term on the right side of the equation is the residual gravity of the water droplet, i.e., the residual part of the water droplet gravity minus the air buoyancy;
Figure 735023DEST_PATH_IMAGE028
and
Figure 296454DEST_PATH_IMAGE029
weight and gravitational acceleration of the water droplets, respectively;
Figure 861428DEST_PATH_IMAGE030
is the water drop density; the second term on the right side of the equation is the gas drag experienced by the water droplet;
Figure 166507DEST_PATH_IMAGE031
and
Figure 500536DEST_PATH_IMAGE032
the instantaneous velocities of the air fluid phase and the water droplet phase respectively;
Figure 232869DEST_PATH_IMAGE033
is the droplet diameter;
Figure 285139DEST_PATH_IMAGE034
is the kinetic viscosity of the air fluid.
Preferably, the method for obtaining the local impact coefficient of the surface of the insulator by the triangular area projection method in S3 comprises the following steps:
Figure 594242DEST_PATH_IMAGE035
；
wherein,
Figure 907411DEST_PATH_IMAGE036
is the local collision coefficient;
Figure 686011DEST_PATH_IMAGE037
、
Figure 22315DEST_PATH_IMAGE038
、
Figure 872459DEST_PATH_IMAGE039
the velocities of the three supercooled water droplets when impacting the surface of the insulator are respectively;
Figure 305715DEST_PATH_IMAGE040
the initial speeds of the three water drops during release are consistent;
Figure 255216DEST_PATH_IMAGE041
is the area of a triangle formed by the initial release positions of the three water drops;
Figure 406712DEST_PATH_IMAGE042
the area of a triangle formed by three water drops after impacting the surface of the insulator.
Preferably, the method further comprises the step of deducing the freezing coefficient based on a heat balance equation, wherein the heat balance equation in the insulator icing process is as follows:
Figure 467072DEST_PATH_IMAGE043
；
Figure 754833DEST_PATH_IMAGE044
；
Figure 140815DEST_PATH_IMAGE045
；
Figure 779607DEST_PATH_IMAGE046
；
Figure 643658DEST_PATH_IMAGE047
；
Figure 51506DEST_PATH_IMAGE048
；
Figure 342810DEST_PATH_IMAGE049
；
Figure 206248DEST_PATH_IMAGE050
；
Figure 608410DEST_PATH_IMAGE051
；
Figure 870764DEST_PATH_IMAGE052
；
Figure 332970DEST_PATH_IMAGE053
；
Figure 680775DEST_PATH_IMAGE054
。
wherein,
Figure 886628DEST_PATH_IMAGE055
the latent heat released in the process that a part of a water drop is frozen into ice at 0 ℃ from water at 0 ℃ for a certain control unit on the surface of the insulator to collide and capture;
Figure 3488DEST_PATH_IMAGE056
the area of the control unit is taken;
Figure 902174DEST_PATH_IMAGE057
、
Figure 737275DEST_PATH_IMAGE058
、
Figure 746819DEST_PATH_IMAGE059
respectively is a collision coefficient, a capture coefficient and a freezing coefficient, and the capture coefficient is constant at 1;
Figure 718187DEST_PATH_IMAGE060
is liquid water content;
Figure 787774DEST_PATH_IMAGE061
is the ambient wind speed;
Figure 110171DEST_PATH_IMAGE062
the latent heat of fusion of ice;
Figure 657827DEST_PATH_IMAGE063
friction heating for air flow;
Figure 749279DEST_PATH_IMAGE064
is the collision kinetic energy of water drops;
Figure 989768DEST_PATH_IMAGE065
freezing at 0 deg.C to the stable temperature of the surface of the ice-coated insulator
Figure 530952DEST_PATH_IMAGE066
Heat released at the time;
Figure 882299DEST_PATH_IMAGE067
is the specific heat of ice;
Figure 31521DEST_PATH_IMAGE068
energy obtained for short wave radiation;
Figure 567544DEST_PATH_IMAGE069
heat loss by convection;
Figure 474320DEST_PATH_IMAGE070
the convection heat transfer coefficient of the ice-coated surface;
Figure 488413DEST_PATH_IMAGE071
is ambient temperature;
Figure 695403DEST_PATH_IMAGE072
for icing the surface in dynamic equilibrium(ii) a temperature;
Figure 136749DEST_PATH_IMAGE073
heat carried away by evaporation of liquid water or sublimation of ice;
Figure 796400DEST_PATH_IMAGE074
is the evaporation or sublimation coefficient;
Figure 614183DEST_PATH_IMAGE075
is at a temperature of
Figure 144522DEST_PATH_IMAGE076
The water surface or the saturated vapor pressure of the ice surface of the ice coating;
Figure 22348DEST_PATH_IMAGE077
at a temperature of
Figure 903716DEST_PATH_IMAGE076
Latent heat of vaporization or sublimation in time;
Figure 525191DEST_PATH_IMAGE078
is the specific heat of air;
Figure 175615DEST_PATH_IMAGE079
is air pressure;
Figure 164955DEST_PATH_IMAGE080
the supercooled state is determined when supercooled water drops in the atmosphere freeze on the surface of the insulator
Figure 658253DEST_PATH_IMAGE071
The heat absorbed by the supercooled water drops is quickly increased to 0 ℃;
Figure 693205DEST_PATH_IMAGE081
specific heat of liquid water;
Figure 588349DEST_PATH_IMAGE082
heat lost as long wave radiation;
Figure 417765DEST_PATH_IMAGE083
in order to obtain the emissivity of the ice surface,
Figure 398359DEST_PATH_IMAGE084
Stefan-Boltzman constant;
Figure 237002DEST_PATH_IMAGE085
to conduct heat loss;
Figure 721073DEST_PATH_IMAGE086
a temperature gradient in the normal direction of thermal conduction;
Figure 986969DEST_PATH_IMAGE087
the heat removed for the unfrozen portion of the supercooled water droplets to leave the ice surface.
Preferably, the method further comprises freezing the coefficient
Figure 189281DEST_PATH_IMAGE088
Comprises the following steps:
Figure 831614DEST_PATH_IMAGE089
。
preferably, the method further comprises, when the ice coating is subjected to wet growth, the ice coating grows along the normal direction of the surface of the insulator; dry growth ice coating grows along the direction of droplet impingement.
Preferably, the method further comprises the ice accretion rate being:
Figure 435771DEST_PATH_IMAGE090
；
wherein,
Figure 872569DEST_PATH_IMAGE060
is liquid water content;
Figure 562176DEST_PATH_IMAGE061
is the ambient wind speed;
Figure 742622DEST_PATH_IMAGE091
the calculation method for the ice coating density is as follows:
Figure 198355DEST_PATH_IMAGE092
；
wherein,
Figure 806054DEST_PATH_IMAGE093
，
Figure 982957DEST_PATH_IMAGE094
is the radius of the water drop,
Figure 701515DEST_PATH_IMAGE095
in order to obtain the surface temperature of the ice coating,
Figure 280264DEST_PATH_IMAGE096
is the wind speed.
Preferably, the method further comprises the step of reconstructing the ice coating morphology by a modeling method from a low dimension to a high dimension through a point-line-plane-body.
Preferably, the method further comprises calculating initial conditions by taking the shape as an air flow field in the next time period after the ice shape reconstruction is completed, and repeating iterative calculation until the ice coating shape and the ice coating amount in the required time are obtained.
Compared with the prior art, the invention has the beneficial effects that:
the numerical simulation method for the ice coating growth without overflow on the surface of the insulator of the power transmission line can not only predict the ice coating growth form, but also predict the ice coating growth thickness and the ice coating quality at any position on the surface of the insulator, realizes stronger accuracy compared with the existing method through the calculation method of the iterative air flow field, can be used for constructing an ice coating early warning mechanism of the power transmission line in an extreme environment, and has good application prospect.
Drawings
FIG. 1 is a flow chart of a numerical simulation method for the overflow-free icing growth on the surface of an insulator of a power transmission line;
FIG. 2 is a schematic view of a triangular area projection method;
FIG. 3 is an exemplary illustration of insulator rod diameter icing reconstruction model;
FIG. 4 is a graphical illustration of ice coating morphology prediction and comparison of tests within one hour of ice coating;
FIG. 5 is a comparison of an insulator icing length test example with a simulation;
fig. 6 is a comparison of insulator icing quality test examples and simulations.
Detailed Description
The present invention is further described with reference to the following detailed description, however, it should be understood by those skilled in the art that the detailed description given herein with respect to the accompanying drawings is for better explanation and that the present invention is not necessarily limited to the specific embodiments, but rather, for equivalent alternatives or common approaches, may be omitted from the detailed description, while still remaining within the scope of the present application.
FIGS. 1 to 6 show preferred embodiments of the present invention, and the present invention will be further described with reference to FIGS. 1 to 6.
As shown in fig. 1: a numerical simulation method for the overflow-free icing growth on the surface of an insulator of a power transmission line comprises the following steps:
s1, establishing an equal-proportion insulator three-dimensional geometric model and constructing an outer flow field area.
And S2, setting boundary conditions of a flow field according to actual meteorological data to obtain a continuous phase air flow field and a motion trail of dispersed phase water drops of the ice-coated insulator.
Specifically, boundary conditions of a flow field are set according to actual meteorological data, and the motion tracks of the insulator ice-coated continuous phase air flow field and the dispersed phase water drops are calculated. The required meteorological data comprise ambient temperature, atmospheric pressure, wind speed, wind direction, liquid water content and the median diameter of the supercooled water drops. The boundary of the airflow inlet is set as a speed inlet with the size of wind speed
Figure 58864DEST_PATH_IMAGE097
And the direction is vertical to the inlet section; selection of turbulence intensity
Figure 988643DEST_PATH_IMAGE008
Turbulence scale
Figure 245311DEST_PATH_IMAGE009
To characterize the turbulence parameters of the turbulent flow field, respectively by empirical formula
Figure 678567DEST_PATH_IMAGE003
And
Figure 628068DEST_PATH_IMAGE004
is determined in which
Figure 779564DEST_PATH_IMAGE010
To calculate the hydraulic diameter of the domain; the outlet is a pressure outlet, and static pressure is set to be 0; for discrete phase supercooled water drops, the water drops are uniformly incident from the entrance of the calculation domain in a surface incidence mode, and the initial speed of the water drops is equal to the speed of free incoming flow. The ice-coated insulator air outflow field can be regarded as normal temperature, low speed, incompressible turbulent flow. Turbulence model adoption
Figure 105503DEST_PATH_IMAGE098
Constructing a model, wherein a calculation formula is as follows:
Figure 393265DEST_PATH_IMAGE012
； （1）
wherein,
Figure 513668DEST_PATH_IMAGE013
is the velocity vector in the flow field;
Figure 152459DEST_PATH_IMAGE014
is the kinetic viscosity of air;
Figure 750931DEST_PATH_IMAGE015
due to the dynamic viscosity additionally generated by air turbulence;
Figure 161708DEST_PATH_IMAGE016
is the air density;
Figure 46488DEST_PATH_IMAGE017
is turbulent kinetic energy;
Figure 782363DEST_PATH_IMAGE018
is the turbulent dissipation ratio;
Figure 574738DEST_PATH_IMAGE019
is the principal stress tensor;
Figure 446879DEST_PATH_IMAGE020
is a volume force;
Figure 564877DEST_PATH_IMAGE021
、
Figure 788048DEST_PATH_IMAGE022
are respectively as
Figure 118535DEST_PATH_IMAGE017
And
Figure 110762DEST_PATH_IMAGE018
the inverse of the effective prandtl constant of (a);
Figure 134081DEST_PATH_IMAGE023
、
Figure 578969DEST_PATH_IMAGE024
、
Figure 713147DEST_PATH_IMAGE025
is a turbulence model parameter;
Figure 825460DEST_PATH_IMAGE026
is a turbulent flow energy term.
When the supercooled water drops move around the insulator, smaller Saffman lifting force, additional mass force, differential pressure force and the like are ignored, the supercooled water drops can be considered to be only under the action of airflow drag force and gravity, and the Lagrangian motion control equation of a single water drop is as follows:
Figure 754101DEST_PATH_IMAGE027
； （2）
where the first term on the right side of the equation is the remaining gravity of the water droplet, i.e., the remaining portion of the water droplet gravity minus the buoyancy of the air.
Figure 300268DEST_PATH_IMAGE028
And
Figure 847924DEST_PATH_IMAGE029
the weight and the gravitational acceleration of the water drop, respectively;
Figure 939377DEST_PATH_IMAGE030
is the water drop density; the second term on the right side of the equation is the gas drag experienced by the water droplets.
Figure 914286DEST_PATH_IMAGE031
And
Figure 723980DEST_PATH_IMAGE032
the instantaneous velocities of the air fluid phase and the water droplet phase respectively;
Figure 75326DEST_PATH_IMAGE099
is the droplet diameter;
Figure 286865DEST_PATH_IMAGE100
is the kinetic viscosity of the air fluid.
And (3) integrating the differential equation of the formula (2) to obtain the velocity distribution of each point of the motion track of the water drop, and then integrating the velocity of the water drop in each time step to obtain the motion track of the water drop, thereby obtaining the physical parameters of any position of the water drop.
S3, extracting three-dimensional coordinates, speed and direction of the whole process from releasing of water drops to capturing of the water drops by the insulator, and obtaining the local collision coefficient of the surface of the insulator through a triangular area projection method
Figure 698255DEST_PATH_IMAGE101
。
Specifically, the triangular area projection method is shown in fig. 2, and the calculation formula is as follows:
Figure 729665DEST_PATH_IMAGE035
； （3）
wherein,
Figure 619123DEST_PATH_IMAGE036
is the local collision coefficient;
Figure 685168DEST_PATH_IMAGE037
、
Figure 1880DEST_PATH_IMAGE038
、
Figure 520586DEST_PATH_IMAGE039
the velocities of the three supercooled water droplets when impacting the surface of the insulator are respectively;
Figure 213736DEST_PATH_IMAGE040
the initial speeds of the three water drops during release are consistent;
Figure 134287DEST_PATH_IMAGE041
is the area of a triangle formed by the initial release positions of the three water drops;
Figure 887479DEST_PATH_IMAGE042
the area of a triangle formed by three water drops after impacting the surface of the insulator.
S4, establishing a thermal equilibrium equation according to the energy conservation principle to determine the freezing coefficient of each position on the surface of the insulator
Figure 896411DEST_PATH_IMAGE059
And judging the icing growth type according to the value of the freezing coefficient value.
Specifically, the heat balance equation in the insulator icing process is as follows:
Figure 127672DEST_PATH_IMAGE043
； （4）
wherein,
Figure 105993DEST_PATH_IMAGE055
the part of the insulator surface where some control unit collides and captures water droplets is released from the latent heat during freezing of 0 ℃ water to 0 ℃ ice:
Figure 154720DEST_PATH_IMAGE044
； （5）
wherein,
Figure 523385DEST_PATH_IMAGE056
the area of the control unit is taken;
Figure 682971DEST_PATH_IMAGE057
、
Figure 187901DEST_PATH_IMAGE058
、
Figure 407530DEST_PATH_IMAGE059
respectively is a collision coefficient, a capture coefficient and a freezing coefficient, and the capture coefficient is constant at 1;
Figure 263490DEST_PATH_IMAGE060
is liquid water content;
Figure 961188DEST_PATH_IMAGE061
is the ambient wind speed;
Figure 586204DEST_PATH_IMAGE062
the latent heat of fusion of ice.
Figure 976735DEST_PATH_IMAGE063
For the friction heating of the air flow, the heating of the ice-coating insulator by the air is generated by the friction of the air flow on the ice surface,since the flow rate of air is not large, this term can be ignored, namely:
Figure 54412DEST_PATH_IMAGE045
。 （6）
Figure 821380DEST_PATH_IMAGE064
collision kinetic energy of water drops:
Figure 300903DEST_PATH_IMAGE046
。 （7）
Figure 859404DEST_PATH_IMAGE065
freezing at 0 deg.C to the stable temperature of the surface of the ice-coated insulator
Figure 424378DEST_PATH_IMAGE066
Heat released in the process:
Figure 463878DEST_PATH_IMAGE047
； （8）
wherein,
Figure 63487DEST_PATH_IMAGE067
is the specific heat of the ice and is,
Figure 795819DEST_PATH_IMAGE102
。
Figure 848089DEST_PATH_IMAGE068
the energy obtained for short-wave radiation is generally ignored because icing generally occurs in fog, rain or cloudy days, without direct sunlight, i.e.:
Figure 956859DEST_PATH_IMAGE048
。 （9）
Figure 410974DEST_PATH_IMAGE069
for convective heat losses:
Figure 48629DEST_PATH_IMAGE049
； （10）
wherein,
Figure 853774DEST_PATH_IMAGE070
the convection heat transfer coefficient of the ice-coated surface;
Figure 500656DEST_PATH_IMAGE071
is ambient temperature;
Figure 809277DEST_PATH_IMAGE072
the surface temperature of the ice coating at the time of dynamic equilibrium of the ice coating.
Figure 883413DEST_PATH_IMAGE073
Heat removed for liquid water evaporation or ice sublimation:
Figure 644695DEST_PATH_IMAGE050
； （11）
wherein,
Figure 95268DEST_PATH_IMAGE074
is the evaporation or sublimation coefficient;
Figure 258396DEST_PATH_IMAGE075
is at a temperature of
Figure 506363DEST_PATH_IMAGE076
The water surface or the saturated vapor pressure of the ice surface of the ice coating;
Figure 20521DEST_PATH_IMAGE103
is at a temperature of
Figure 9205DEST_PATH_IMAGE076
Latent heat of vaporization or sublimation in time;
Figure 26840DEST_PATH_IMAGE078
is the specific heat of air;
Figure 708357DEST_PATH_IMAGE079
is air pressure.
Figure 444232DEST_PATH_IMAGE080
The supercooled state is determined when supercooled water drops in the atmosphere freeze on the surface of the insulator
Figure 971028DEST_PATH_IMAGE071
Heat absorbed by the supercooled water droplets during rapid rise to 0 ℃:
Figure 108748DEST_PATH_IMAGE051
； （12）
wherein,
Figure 961167DEST_PATH_IMAGE081
is the specific heat of liquid water.
Figure 246654DEST_PATH_IMAGE082
Heat lost for long wave radiation:
Figure 718087DEST_PATH_IMAGE052
； （13）
wherein,
Figure 834948DEST_PATH_IMAGE083
in order to obtain the emissivity of the ice surface,
Figure 264792DEST_PATH_IMAGE104
，
Figure 303155DEST_PATH_IMAGE084
is the Stefan-Boltzman constant,
Figure 312699DEST_PATH_IMAGE105
。
Figure 956170DEST_PATH_IMAGE085
to conduct heat loss:
Figure 350724DEST_PATH_IMAGE053
； （14）
wherein,
Figure 407542DEST_PATH_IMAGE017
is the thermal conductivity of the medium and,
Figure 220777DEST_PATH_IMAGE086
is the temperature gradient in the normal direction of the heat conduction.
Figure 312230DEST_PATH_IMAGE087
The heat removed for the unfrozen portion of the supercooled water droplets to leave the ice surface:
Figure 287139DEST_PATH_IMAGE054
。 （15）
freezing coefficient
Figure 96832DEST_PATH_IMAGE088
Comprises the following steps:
Figure 448179DEST_PATH_IMAGE089
。 （16）
the icing type under the condition of no overflow on the surface of the insulator can be judged by the existence of a water film in a freezing area, and when the freezing coefficient is
Figure 394138DEST_PATH_IMAGE106
Ice coating for wet growth, and when
Figure 539949DEST_PATH_IMAGE107
The ice coating is dry and long. When the ice coating is in wet growth, the ice coating grows along the normal direction of the surface of the insulator, and when the ice coating is in dry growth, the ice coating grows along the collision direction of liquid drops.
And S5, determining the icing growth within a time step, and reconstructing the icing boundary by a point-line-surface-body modeling method to obtain a three-dimensional model of the insulator icing growth form.
Specifically, the icing growth thickness within a time step is calculated according to an icing growth rate calculation formula, wherein the icing growth rate formula is as follows:
Figure 571359DEST_PATH_IMAGE090
； （17）
wherein,
Figure 991976DEST_PATH_IMAGE060
is liquid water content;
Figure 792441DEST_PATH_IMAGE061
is the ambient wind speed;
Figure 374732DEST_PATH_IMAGE091
the calculation method for the ice coating density is as follows:
Figure 893438DEST_PATH_IMAGE092
； （18）
wherein,
Figure 586588DEST_PATH_IMAGE093
，
Figure 510069DEST_PATH_IMAGE094
is the radius of the water drop,
Figure 591158DEST_PATH_IMAGE095
for icingThe temperature of the surface of the steel sheet is measured,
Figure 206947DEST_PATH_IMAGE096
is the wind speed.
As shown in fig. 3, two curves are formed by connecting adjacent points of the collision point and the newly formed point after the growth, and the end points of the two curves are connected together to form a two-dimensional curved surface. And constructing a three-dimensional entity by connecting the boundaries of each group of curved surfaces, thereby constructing a new ice coating boundary.
And S6, taking the new icing form as an initial condition of the air flow field in the next time step, and repeating the steps to perform the circulation iteration of icing increase until the icing form in the required time is obtained.
Specifically, when the ambient temperature is-10 ℃, the wind speed is 10m/s, and the median diameter of the liquid drop is 50
Figure 828421DEST_PATH_IMAGE108
The liquid water content in the air is
Figure 478845DEST_PATH_IMAGE109
When the temperature of the water is higher than the set temperature,
Figure 527573DEST_PATH_IMAGE110
the prediction and experimental comparison of ice coating morphology within one hour of the insulator is shown in fig. 4. Fig. 5 is a comparison of the ice coating growth length at the edge of the shed and the rod diameter of the composite insulator and the test, and fig. 6 is a comparison of the ice coating quality of the composite insulator and the test.
The method for simulating the ice-coating growth numerical value without overflow on the surface of the insulator of the power transmission line has the advantages of high prediction precision, simplicity, effectiveness and small error. The technical scheme of the invention can be applied to the insulator of the type disclosed by the patent and can be expanded to insulators of any type.
The foregoing is directed to preferred embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow. However, any simple modification, equivalent change and modification of the above embodiments according to the technical essence of the present invention are within the protection scope of the technical solution of the present invention.

Claims (10)

1. A numerical simulation method for the increase of ice coating without overflow on the surface of an insulator of a power transmission line is characterized by comprising the following steps: the method comprises the following steps:
s1, establishing an equal-proportion insulator three-dimensional geometric model and constructing an outer flow field area;
s2, setting boundary conditions of a flow field according to actual meteorological data to obtain a continuous phase air flow field of the ice-coated insulator and a motion trail of dispersed phase water drops;
s3, extracting three-dimensional coordinates, speed and direction of the water drops in the whole process from releasing to being captured by the insulator, and obtaining a local collision coefficient of the surface of the insulator through a triangular area projection method;
s4, constructing a thermal balance equation according to an energy conservation principle to determine the freezing coefficient of each position on the surface of the insulator, and judging the icing growth type according to the freezing coefficient value;
s5, determining the icing growth within a time step, and reconstructing an icing boundary by a point-line-surface-body modeling method to obtain a three-dimensional model of the icing growth form of the insulator;
and S6, taking the three-dimensional model as an initial condition of the air flow field in the next time step, and repeating the steps to perform the cycle iteration of ice coating growth until the ice coating form in the required time is obtained.
2. The method for simulating the overflow-free icing growth numerical value on the surface of the insulator of the power transmission line according to claim 1, wherein the method comprises the following steps: the boundary conditions in S2 are:
Figure DEST_PATH_IMAGE001
Figure 236065DEST_PATH_IMAGE002
Figure DEST_PATH_IMAGE003
；
Figure 683358DEST_PATH_IMAGE004
；
wherein,
Figure 187151DEST_PATH_IMAGE005
is the inlet velocity of the turbulent flow field,
Figure 466252DEST_PATH_IMAGE006
is the ambient wind speed and is,
Figure DEST_PATH_IMAGE007
in order to obtain the outlet pressure of the turbulent flow field,
Figure 623564DEST_PATH_IMAGE008
in order to be the intensity of the turbulent flow,
Figure DEST_PATH_IMAGE009
in order to be of a turbulent flow scale,
Figure 89311DEST_PATH_IMAGE010
in order to calculate the hydraulic diameter of the domain,
Figure 131217DEST_PATH_IMAGE011
is the reynolds number.
3. The method for simulating the numerical value of the ice accretion on the surface of the insulator of the electric transmission line according to claim 1, wherein: the method also comprises the following steps that the ice-coated insulator air external flow field is regarded as normal-temperature, low-speed and incompressible turbulent flow, and a turbulent flow model is established:
Figure 975807DEST_PATH_IMAGE012
；
wherein,
Figure 444966DEST_PATH_IMAGE013
is the velocity vector in the flow field;
Figure 584960DEST_PATH_IMAGE014
is the kinetic viscosity of air;
Figure 758452DEST_PATH_IMAGE015
due to the dynamic viscosity additionally generated by air turbulence;
Figure 582183DEST_PATH_IMAGE016
is the air density;
Figure 487822DEST_PATH_IMAGE017
is turbulent kinetic energy;
Figure 443009DEST_PATH_IMAGE018
is the turbulent dissipation ratio;
Figure 311870DEST_PATH_IMAGE019
is the principal stress tensor;
Figure 177058DEST_PATH_IMAGE020
is a volume force;
Figure 237286DEST_PATH_IMAGE021
、
Figure 181234DEST_PATH_IMAGE022
are respectively as
Figure 368633DEST_PATH_IMAGE017
And
Figure 540857DEST_PATH_IMAGE018
the inverse of the effective prandtl constant of (a);
Figure 319457DEST_PATH_IMAGE023
、
Figure 813018DEST_PATH_IMAGE024
、
Figure 725479DEST_PATH_IMAGE025
is a turbulence model parameter;
Figure 502942DEST_PATH_IMAGE026
are turbulent flow energy terms.
4. The method for simulating the overflow-free icing growth numerical value on the surface of the insulator of the power transmission line according to claim 1, wherein the method comprises the following steps: the method further comprises the following step that the Lagrangian motion control equation of the single water drop is as follows:
Figure 203176DEST_PATH_IMAGE027
wherein,
Figure 636562DEST_PATH_IMAGE028
and
Figure 477348DEST_PATH_IMAGE029
the weight and the gravitational acceleration of the water drop, respectively;
Figure 843739DEST_PATH_IMAGE030
is the water drop density;
Figure DEST_PATH_IMAGE031
and
Figure 777191DEST_PATH_IMAGE032
the instantaneous velocities of the air fluid phase and the water droplet phase respectively;
Figure DEST_PATH_IMAGE033
is the droplet diameter;
Figure 979764DEST_PATH_IMAGE034
is the dynamic viscosity of the air fluid.
5. The method for simulating the overflow-free icing growth numerical value on the surface of the insulator of the power transmission line according to claim 1, wherein the method comprises the following steps: the method for obtaining the local collision coefficient of the surface of the insulator through the triangular area projection method in the S3 comprises the following steps:
Figure DEST_PATH_IMAGE035
；
wherein,
Figure 171711DEST_PATH_IMAGE036
is the local collision coefficient;
Figure DEST_PATH_IMAGE037
、
Figure 346603DEST_PATH_IMAGE038
、
Figure DEST_PATH_IMAGE039
the velocities of the three supercooled water droplets when impacting the surface of the insulator are respectively;
Figure 293699DEST_PATH_IMAGE040
the initial speeds of the three water drops during release are consistent;
Figure DEST_PATH_IMAGE041
is the area of a triangle formed by the initial release positions of the three water drops;
Figure 717990DEST_PATH_IMAGE042
the area of a triangle formed by three water drops after impacting the surface of the insulator.
6. The method for simulating the numerical value of the ice accretion on the surface of the insulator of the electric transmission line according to claim 1, wherein: the method further comprises the step of deducing the freezing coefficient based on a heat balance equation, wherein the heat balance equation in the insulator icing process is as follows:
Figure DEST_PATH_IMAGE043
；
Figure 775944DEST_PATH_IMAGE044
；
Figure DEST_PATH_IMAGE045
；
Figure 726714DEST_PATH_IMAGE046
；
Figure 720078DEST_PATH_IMAGE047
；
Figure 897243DEST_PATH_IMAGE048
；
Figure 837517DEST_PATH_IMAGE049
；
Figure 282274DEST_PATH_IMAGE050
；
Figure 384222DEST_PATH_IMAGE051
；
Figure 986367DEST_PATH_IMAGE052
；
Figure 714021DEST_PATH_IMAGE053
；
Figure 560754DEST_PATH_IMAGE054
；
wherein,
Figure 584336DEST_PATH_IMAGE055
the method is characterized in that a part of a control unit on the surface of an insulator, which collides with and captures water drops, is released from latent heat in the process of freezing water at 0 ℃ into ice at 0 ℃;
Figure 985361DEST_PATH_IMAGE056
the area of the control unit is taken;
Figure 923230DEST_PATH_IMAGE057
、
Figure 86186DEST_PATH_IMAGE058
、
Figure 326674DEST_PATH_IMAGE059
respectively is a collision coefficient, a capture coefficient and a freezing coefficient, and the capture coefficient is constant at 1;
Figure 729843DEST_PATH_IMAGE060
is liquid water content;
Figure 284452DEST_PATH_IMAGE061
is the ambient wind speed;
Figure 325351DEST_PATH_IMAGE062
the latent heat of fusion of ice;
Figure 940004DEST_PATH_IMAGE063
friction heating for air flow;
Figure 971413DEST_PATH_IMAGE064
is the collision kinetic energy of water drops;
Figure 346025DEST_PATH_IMAGE065
freezing at 0 deg.C to the stable temperature of the surface of the ice-coated insulator
Figure 490699DEST_PATH_IMAGE066
Heat released at the time;
Figure 791099DEST_PATH_IMAGE067
is the specific heat of ice;
Figure 450750DEST_PATH_IMAGE068
energy obtained for short wave radiation;
Figure 363474DEST_PATH_IMAGE069
heat loss by convection;
Figure 628233DEST_PATH_IMAGE070
the convection heat transfer coefficient of the ice-coated surface;
Figure 99535DEST_PATH_IMAGE071
is ambient temperature;
Figure 184165DEST_PATH_IMAGE072
the surface temperature of the ice coating during the dynamic equilibrium of the ice coating;
Figure 635000DEST_PATH_IMAGE073
heat removed for liquid water evaporation or ice sublimation;
Figure 754266DEST_PATH_IMAGE074
is the evaporation or sublimation coefficient;
Figure 396469DEST_PATH_IMAGE075
is at a temperature of
Figure 391232DEST_PATH_IMAGE076
The water surface or the saturated vapor pressure of the ice surface of the ice coating;
Figure 691763DEST_PATH_IMAGE077
is at a temperature of
Figure 914803DEST_PATH_IMAGE076
Latent heat of vaporization or sublimation in time;
Figure 213060DEST_PATH_IMAGE078
is the specific heat of air;
Figure 23016DEST_PATH_IMAGE079
is the air pressure;
Figure 251872DEST_PATH_IMAGE080
the supercooled state being a state in which supercooled water droplets in the atmosphere are frozen on the surface of the insulator
Figure 80151DEST_PATH_IMAGE071
The heat absorbed by the supercooled water drops is rapidly increased to 0 ℃;
Figure 831200DEST_PATH_IMAGE081
specific heat of liquid water;
Figure 112140DEST_PATH_IMAGE082
heat lost to long wave radiation;
Figure 207004DEST_PATH_IMAGE083
in order to obtain the emissivity of the ice surface,
Figure 952106DEST_PATH_IMAGE084
Stefan-Boltzman constant;
Figure 342898DEST_PATH_IMAGE085
to conduct heat loss;
Figure 298085DEST_PATH_IMAGE086
a temperature gradient in the normal direction of heat conduction;
Figure 947372DEST_PATH_IMAGE087
the heat removed for the unfrozen portion of the supercooled water droplets to leave the ice surface.
7. The method for simulating the overflow-free icing growth numerical value on the surface of the insulator of the power transmission line according to claim 6, wherein the method comprises the following steps: the method further comprises freezing the coefficients
Figure 766554DEST_PATH_IMAGE088
Comprises the following steps:
Figure DEST_PATH_IMAGE089
。
8. the method for simulating the overflow-free icing growth numerical value on the surface of the insulator of the power transmission line according to claim 1, wherein the method comprises the following steps: the method further comprises, during wet growth icing, the icing grows normally along the surface of the insulator; dry growth ice coating grows along the direction of droplet impingement.
9. The method for simulating the overflow-free icing growth numerical value on the surface of the insulator of the power transmission line according to claim 1, wherein the method comprises the following steps: the method further comprises the following steps of:
Figure 764466DEST_PATH_IMAGE090
；
wherein,
Figure 816736DEST_PATH_IMAGE060
is liquid water content;
Figure 20446DEST_PATH_IMAGE061
is the ambient wind speed;
Figure 677824DEST_PATH_IMAGE091
the ice density is calculated as follows:
Figure 174533DEST_PATH_IMAGE092
；
wherein,
Figure 340197DEST_PATH_IMAGE093
，
Figure 331287DEST_PATH_IMAGE094
is the radius of the water drop,
Figure 358018DEST_PATH_IMAGE095
in order to obtain the surface temperature of the ice coating,
Figure 776361DEST_PATH_IMAGE096
is the wind speed.
10. The method for simulating the overflow-free icing growth numerical value on the surface of the insulator of the power transmission line according to claim 1, wherein the method comprises the following steps: the method also comprises that the reconstruction of the ice coating form is realized by a modeling method from low dimension to high dimension through point-line-plane-body.
CN202211098383.1A 2022-09-09 2022-09-09 Numerical simulation method for overflow-free icing growth on surface of power transmission line insulator Active CN115169267B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211098383.1A CN115169267B (en) 2022-09-09 2022-09-09 Numerical simulation method for overflow-free icing growth on surface of power transmission line insulator

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211098383.1A CN115169267B (en) 2022-09-09 2022-09-09 Numerical simulation method for overflow-free icing growth on surface of power transmission line insulator

Publications (2)

Publication Number Publication Date
CN115169267A true CN115169267A (en) 2022-10-11
CN115169267B CN115169267B (en) 2022-11-22

Family

ID=83482493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211098383.1A Active CN115169267B (en) 2022-09-09 2022-09-09 Numerical simulation method for overflow-free icing growth on surface of power transmission line insulator

Country Status (1)

Country Link
CN (1) CN115169267B (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101710057A (en) * 2009-11-11 2010-05-19 东北电力大学 Ice-coated test system for transmission line structure
CN103604419A (en) * 2013-11-21 2014-02-26 中国科学院西安光学精密机械研究所 Three-dimensional reconstruction-based wire icing monitoring method and system
US20140067271A1 (en) * 2012-08-30 2014-03-06 International Business Machines Corporation Predicting ice coating status on transmission lines
CN104655030A (en) * 2015-02-16 2015-05-27 国网安徽省电力公司铜陵供电公司 Power transmission line icing detecting and early-warning device
CN104913900A (en) * 2015-06-03 2015-09-16 西安工程大学 Measurement method of power transmission line iced conductor galloping aerodynamic parameters
US20160169761A1 (en) * 2014-12-10 2016-06-16 Uchicago Argonne, Llc Method and system for icing condition detection
CN106680603A (en) * 2017-01-17 2017-05-17 清华大学 Method for calculation of characteristic lines of three-dimensional ionized field of direct current transmission line
CN109242159A (en) * 2018-08-17 2019-01-18 华北电力大学 A kind of novel power transmission transmission line icing prediction technique
CN111967201A (en) * 2020-07-30 2020-11-20 华南理工大学 Method for analyzing critical icing type based on numerical simulation model
CN112651188A (en) * 2020-12-02 2021-04-13 国网重庆市电力公司电力科学研究院 CFD-based power transmission line icing state simulation method
CN113591342A (en) * 2021-06-28 2021-11-02 国网天津市电力公司电力科学研究院 10kV overhead insulated line galloping characteristic analysis method under different wind load and ice coating conditions

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101710057A (en) * 2009-11-11 2010-05-19 东北电力大学 Ice-coated test system for transmission line structure
US20140067271A1 (en) * 2012-08-30 2014-03-06 International Business Machines Corporation Predicting ice coating status on transmission lines
CN103604419A (en) * 2013-11-21 2014-02-26 中国科学院西安光学精密机械研究所 Three-dimensional reconstruction-based wire icing monitoring method and system
US20160169761A1 (en) * 2014-12-10 2016-06-16 Uchicago Argonne, Llc Method and system for icing condition detection
CN104655030A (en) * 2015-02-16 2015-05-27 国网安徽省电力公司铜陵供电公司 Power transmission line icing detecting and early-warning device
CN104913900A (en) * 2015-06-03 2015-09-16 西安工程大学 Measurement method of power transmission line iced conductor galloping aerodynamic parameters
CN106680603A (en) * 2017-01-17 2017-05-17 清华大学 Method for calculation of characteristic lines of three-dimensional ionized field of direct current transmission line
CN109242159A (en) * 2018-08-17 2019-01-18 华北电力大学 A kind of novel power transmission transmission line icing prediction technique
CN111967201A (en) * 2020-07-30 2020-11-20 华南理工大学 Method for analyzing critical icing type based on numerical simulation model
CN112651188A (en) * 2020-12-02 2021-04-13 国网重庆市电力公司电力科学研究院 CFD-based power transmission line icing state simulation method
CN113591342A (en) * 2021-06-28 2021-11-02 国网天津市电力公司电力科学研究院 10kV overhead insulated line galloping characteristic analysis method under different wind load and ice coating conditions

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BO LI: "Real-time On-line Anti-icing and Ice-melting Power Calculation Method for Transmission Lines Based on Self-ice-melting Conductors", 《2019 IEEE 3RD CONFERENCE ON ENERGY INTERNET AND ENERGY SYSTEM INTEGRATION (EI2)》 *
杨铭等: "输电线路覆冰原因分析及优化研究", 《电气传动自动化》 *
蒋兴良等: "绝缘子湿增长动态覆冰模型研究", 《中国电机工程学报》 *

Also Published As

Publication number Publication date
CN115169267B (en) 2022-11-22

Similar Documents

Publication Publication Date Title
Makkonen Models for the growth of rime, glaze, icicles and wet snow on structures
Ibrahim et al. Effects of blade design on ice accretion for horizontal axis wind turbines
Kraj et al. Measurement method and results of ice adhesion force on the curved surface of a wind turbine blade
Jin et al. Study of ice accretion and icing effects on aerodynamic characteristics of DU96 wind turbine blade profile
CN111967201B (en) Method for analyzing critical icing type based on numerical simulation model
EP3845465B1 (en) Aerostat icing characteristic numerical simulation and experimental verification system
Reid et al. FENSAP-ICE simulation of icing on wind turbine blades, part 1: performance degradation
CN111396269A (en) Multi-time-step unsteady icing calculation method and system and storage medium
CN114139393B (en) Part electric heating three-dimensional anti-icing numerical simulation method considering water film flow heat transfer
CA2430823A1 (en) Morphogenetic modelling of in-flight icing
CN111967147B (en) Method for calculating key parameters of icing growth model of overhead transmission line
CN109598033B (en) Icing parameter calculation method based on rotating multiple conductors
CN114169077A (en) Strong-coupling three-dimensional numerical simulation method for hot gas anti-icing of aircraft engine inlet part
CN108108533B (en) Method for calculating icing weight of power transmission conductor
Homola et al. Turbine size and temperature dependence of icing on wind turbine blades
CN117010204A (en) Wing icing calculation method and system, storage medium and electronic equipment
Virk et al. Relation between angle of attack and atmospheric ice accretion on large wind turbine's blade
Xu et al. Effects of ambient temperature and wind speed on icing characteristics and anti-icing energy demand of a blade airfoil for wind turbine
CN115169267B (en) Numerical simulation method for overflow-free icing growth on surface of power transmission line insulator
Hu et al. An experimental investigation on the unsteady heat transfer process over an ice accreting NACA 0012 airfoil
Lébatto et al. Conductor icing: Comparison of a glaze icing model with experiments under severe laboratory conditions with moderate wind speed
Virk et al. Multiphysics based numerical study of atmospheric ice accretion on a full scale horizontal axis wind turbine blade
Chitransh et al. Investigation of Single Shot Ice Accretion on Aerofoil of Wind Turbine Blade using ANSYS
Makkonen et al. Numerical modelling of icing on power network equipment
Virk Atmospheric icing of transmission line circular conductor bundles in triplex configuration

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant