A method, device and medium for calculating water droplet capture rate of wire icing considering electric field effect

Abstract

The present application relates to the technical field of power system transmission line anti-icing, and discloses a kind of wire icing water droplet capture rate calculation method, equipment and medium considering electric field effect, comprising, by obtaining wire, environment and electrical parameter, water droplet charge amount model is established;Water droplet comprehensive force model is constructed;Based on water droplet comprehensive force model, water droplet motion equation is established, and the motion trajectory of water droplet in different initial position is solved using numerical method;By establishing icing wire geometric model, the collision situation of each water droplet trajectory and icing surface is judged and position is recorded;According to the collision statistical result, the local capture rate and overall capture rate of each region of icing wire are calculated.The present application quantifies the influence of electric field on water droplet motion, significantly improves the prediction accuracy of charged wire icing growth, and can provide reliable basis for the anti-icing design and safe operation of transmission line.

Description

Technical Field

[0001] This invention relates to the field of anti-icing technology for power transmission lines, and in particular to a method, equipment, and medium for calculating the ice droplet capture rate of conductors considering the electric field effect. Background Technology

[0002] Icing on transmission lines is a significant factor threatening the safe operation of power grids. In freezing rain and snow, supercooled water droplets impact and freeze on the surface of conductors, forming icing. In severe cases, this can lead to conductor galloping, breakage, and collapse, causing widespread power outages and substantial economic losses. The water droplet capture rate is a key parameter describing the icing process on conductors. Accurate calculation of the water droplet capture rate is crucial for predicting the icing growth rate and morphology, as well as assessing the risk of icing disasters. Traditional methods for calculating the water droplet capture rate are applicable to icing calculations under conditions without an electric field, such as aircraft wing icing or wind turbine blade icing. However, actual high-voltage transmission lines operate in strong electric field environments, and the influence of the electric field on water droplet movement cannot be ignored.

[0003] Existing technologies lack a comprehensive method for calculating the water droplet capture rate that takes into account the aforementioned electric field effects. This results in insufficient accuracy in predicting icing on energized conductors and fails to accurately reflect the icing characteristics under actual energized operating conditions. Therefore, there is an urgent need to develop a method for calculating the water droplet capture rate of conductors that considers the electric field effect, in order to improve the accuracy of icing prediction and provide a scientific basis for the anti-icing design and operation and maintenance decisions of transmission lines. Summary of the Invention

[0004] In view of the above-mentioned existing problems, the present invention provides a method, device and medium for calculating the ice droplet capture rate of wires considering the electric field effect.

[0005] Therefore, the technical problem solved by this invention is: by establishing a force model of water droplets that comprehensively considers gravity, air viscous drag, electric field polarization force and corona wind force, and combining the water droplet charge characteristics, the motion trajectory of water droplets around high-voltage conductors can be solved, and then the capture rate of water droplets on the ice-covered surface of the conductor and various parts of the ice floe can be calculated, providing a more accurate calculation method for predicting ice accumulation on charged conductors.

[0006] To address the aforementioned technical problems, this invention provides the following technical solution: a method for calculating the capture rate of ice-covered water droplets on a conductor considering the electric field effect, comprising: acquiring conductor parameters, environmental parameters, and electrical parameters; determining the charge of the water droplets based on the droplet diameter and electric field strength, and establishing a water droplet charge model; establishing a water droplet force model and calculating the net force acting on the water droplets during their motion; establishing the water droplet motion equation based on the water droplet force model and the water droplet charge model, and using numerical methods to solve for the motion trajectory of water droplets at different initial positions; establishing a geometric model of the ice-covered conductor, determining the collision situation between each water droplet's motion trajectory and the ice-covered surface, and recording the collision positions; and calculating the local capture rate and overall capture rate of each region of the ice-covered conductor based on the collision statistics.

[0007] As a preferred embodiment of the method for calculating the capture rate of ice-covered water droplets on a conductor considering the electric field effect described in this invention, the step of establishing the water droplet charge model includes: experimentally measuring the charge of water droplets of different diameters under different electric field intensities, establishing a relationship model between the water droplet charge and the water droplet diameter and electric field intensity, and determining the model parameters using a fitting method.

[0008] As a preferred embodiment of the method for calculating the capture rate of iced water droplets on a conductor considering the electric field effect described in this invention, the method for establishing the force model of the water droplet includes obtaining the resultant force on the water droplet during its motion by vector superimposing gravity, air viscous drag, electric field polarization force, corona wind force and Coulomb force, wherein the air viscous drag coefficient is determined piecewise based on the water droplet Reynolds number.

[0009] As a preferred embodiment of the method for calculating the capture rate of ice-covered water droplets on a conductor considering the electric field effect described in this invention, the method for solving the water droplet trajectory includes setting a water droplet release surface upstream of the conductor, with the release surface perpendicular to the direction of the incoming flow. Water droplet release points are evenly arranged on the release surface at preset intervals; Based on the median diameter and particle size distribution of water droplets, water droplets of different diameters are allocated to each release point; Calculate the trajectory of the water droplets released at each release point under the action of the resultant force.

[0010] As a preferred embodiment of the method for calculating the ice-covered water droplet capture rate of a conductor considering the electric field effect described in this invention, the step of establishing the geometric model of the ice-covered conductor includes establishing a geometric model of the cylindrical surface of the conductor based on the conductor diameter. Based on the type and thickness of icing, a geometric model of the icing layer is established on the surface of the conductor; When it is determined to be rime ice, an ice geometric model is established below the ice layer. The ice geometric model includes parameters such as ice length, ice diameter, and ice spacing.

[0011] As a preferred embodiment of the method for calculating the water droplet capture rate of an ice-covered conductor considering the electric field effect described in this invention, the method for determining the collision between the trajectory of each water droplet and the ice-covered surface includes determining whether the position of the water droplet enters the boundary of the geometric model of the ice-covered conductor within each calculation time step. When a water droplet enters the boundary of the ice-covered surface, it is determined that the current water droplet has been captured, and the collision position coordinates are recorded. If the trajectory of a water droplet bypasses the icing wire or reaches the boundary of the computational domain without intersecting the icing surface, it is determined that the current water droplet has not been captured.

[0012] As a preferred embodiment of the method for calculating the water droplet capture rate of an ice-covered conductor considering the electric field effect described in this invention, the calculation of the local capture rate and the overall capture rate of each region of the ice-covered conductor includes dividing the surface of the ice-covered conductor into several regions, including the windward side region of the conductor, the leeward side region of the conductor, the surface region of the ice layer, the tip region of the ice floe, and the side region of the ice floe. Local capture rate β in each region i The calculation formula is: In the formula, N ci N represents the number of water droplets captured in the i-th region. 0i The total number of water droplets within the projection range of the i-th region in the incoming flow; The formula for calculating the overall capture rate β is: In the formula, N0 is the total number of water droplets in the incoming flow corresponding to the entire projection range of the icing guide, and n is the total number of regions divided.

[0013] As a preferred embodiment of the method for calculating the capture rate of ice-covered water droplets on a conductor considering the electric field effect described in this invention, the method further includes comparing and verifying the calculated local capture rate and overall capture rate with the results of artificial climate chamber experiments; and correcting the parameters in the water droplet force model or water droplet charge model based on the comparison results to improve the calculation accuracy.

[0014] The present invention provides a computer device, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps of a method for calculating the capture rate of ice-covered water droplets on a wire considering the electric field effect.

[0015] The present invention provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of a method for calculating the capture rate of ice-covered water droplets on a wire considering the electric field effect.

[0016] The beneficial effects of this invention are as follows: By establishing a force model for water droplets that includes electric field polarization force, corona wind force, and Coulomb force, and combining this with numerical simulation of the water droplet's charge characteristics and trajectory, the icing process of high-voltage transmission lines under energized operating conditions can be predicted more accurately. This overcomes the shortcomings of traditional methods that do not consider the influence of the electric field, significantly improves the accuracy of icing growth prediction, and provides a reliable theoretical basis and calculation tool for the anti-icing design and safe operation and maintenance of transmission lines. Attached Figure Description

[0017] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the following description of the embodiments will be briefly introduced. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0018] Figure 1 This is a schematic flowchart illustrating a method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect, as provided in one embodiment of the present invention.

[0019] Figure 2 This is a schematic diagram of the force analysis of water droplets in a method for calculating the water droplet capture rate of an ice-covered conductor considering the electric field effect, provided as an embodiment of the present invention.

[0020] Figure 3 This is a schematic diagram of the water droplet trajectory calculation results provided by a method for calculating the capture rate of ice-covered water droplets on a conductor considering the electric field effect, according to an embodiment of the present invention.

[0021] Figure 4 This is a schematic diagram of the local capture rate distribution curve along the ice-covered surface, provided by an embodiment of the present invention for a method of calculating the capture rate of ice-covered water droplets on a conductor considering the electric field effect. Detailed Implementation

[0022] To make the above-mentioned objects, features, and advantages of the present invention more apparent and understandable, specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of them. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the protection scope of the present invention.

[0023] Example 1, referring to Figure 1 This is the first embodiment of the present invention, which provides a method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect, including: S1: Obtain conductor parameters, environmental parameters, and electrical parameters; determine the charge of the water droplet based on the droplet diameter and electric field strength; and establish a water droplet charge model.

[0024] S2: Establish a force model of the water droplet and calculate the resultant force on the water droplet during its motion.

[0025] S3: Based on the force model and charge model of the water droplet, establish the motion equation of the water droplet, and use numerical methods to solve the motion trajectory of the water droplet at different initial positions.

[0026] S4: Establish a geometric model of the ice-covered guide wire, determine the collision between the trajectory of each water droplet and the ice-covered surface, and record the collision position.

[0027] S5: Based on the collision statistics, calculate the local capture rate and overall capture rate of each region of the icy conductor.

[0028] It should be noted that, compared with traditional methods for calculating the water droplet capture rate that do not consider the influence of the electric field, this method establishes a water droplet force model that includes electric field polarization force, corona wind force and Coulomb force, and couples the water droplet charge characteristics with the motion trajectory analysis. For the first time, it has achieved an accurate quantitative description of the icing process of high-voltage transmission lines under energized operation conditions, thus solving the key problem that the existing technology's ignoring of the electric field effect leads to discrepancies between the icing prediction results and the actual operating conditions.

[0029] Example 2, refer to Figure 2 - Figure 4 As an embodiment of the present invention, based on the above embodiment, a method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect is provided.

[0030] Furthermore, in this embodiment, step S1 involves acquiring conductor parameters, environmental parameters, and electrical parameters, determining the water droplet's charge based on the droplet diameter and electric field strength, and establishing a water droplet charge model. Specific steps include: The conductor parameters include: conductor diameter D, in meters; conductor surface roughness; icing geometry, including icing thickness δ and icing type (rime, glaze, or mixed glaze); if it is glaze icing, it also includes parameters such as icicle length, icicle root radius r0, and icicle spacing.

[0031] Environmental parameters include: ambient temperature T (°C); incoming air velocity U∞ (m / s); median water droplet diameter dMVD (μm); liquid water content W (g / m³); air density ρa (kg / m³); and aerodynamic viscosity μa (Pa·s).

[0032] The electrical parameters include: conductor operating voltage V, in kV; electric field intensity distribution E(x,y) around the conductor, which can be obtained through finite element simulation or analytical formula calculation; corona initiation voltage and corona wind speed distribution vc(x,y), which can be obtained through corona discharge simulation or experimental measurement.

[0033] Water droplets acquire charge in a high-voltage electric field through various mechanisms, including: inductive charging (charge separation caused by electric field polarization), collisional charging (charge transfer when colliding with other particles or surfaces), and corona charging (capturing ions generated by corona discharge).

[0034] The charge of a water droplet is related to factors such as its diameter and the electric field strength. This invention establishes a model relating the water droplet charge q to its diameter d and the electric field strength E by constructing an experimental platform for measuring the charge of water droplets of different diameters under different electric field strengths. In the formula, α, β, and γ are fitting coefficients, which are determined based on experimental data using the least squares method or other optimization algorithms.

[0035] In an alternative embodiment, the water droplet charge model can also be established using a polynomial form or by establishing a neural network model to fit the experimental data, thereby obtaining a more accurate charge prediction.

[0036] In another alternative embodiment, the water droplet charge model can be constructed by using experimentally measured water droplet diameter and electric field strength as input features and charge as the output label to build a dataset. Classic machine learning algorithms such as random forest regression are then used for training to establish a black-box prediction model from the input parameters to the charge.

[0037] Furthermore, in this embodiment of the application, step S2 establishes a force model of the water droplet and calculates the resultant force on the water droplet during its motion. Specific steps include: As a water droplet moves around a conductor, it is subjected to various forces. This invention establishes a force model that comprehensively considers gravity, air viscous resistance, electric field polarization force, and corona wind force.

[0038] The net force F acting on the water droplet is expressed as: In the formula, F g For gravity; F d For air viscous drag; F dep This is the electric field polarization force (dielectrophoresis force); F c The force of the corona wind; F q This refers to Coulomb force. All units mentioned above are in N (units unspecified). (1) Gravity Fg Gravity is the gravitational force exerted on a water droplet by the Earth, directed vertically downwards, and is represented as: In the formula, Let r be the mass of the water droplet, kg; and r be the radius of the water droplet, m. Let g be the density of water, taken as 998 kg / m³; g is the gravitational acceleration vector, g = 9.81 m / s², directed vertically downwards.

[0039] (2) Air viscous resistance Fd Air viscous drag is the resistance experienced by a water droplet relative to the air. Its direction is opposite to the direction of the droplet's velocity relative to the air, and it is expressed as: In the formula, C d The drag coefficient is dimensionless. air density, kg / m³; Let be the velocity of the air relative to the water droplet, in m / s; The local air velocity is in m / s; Let be the velocity of the water droplet, in m / s.

[0040] Drag coefficient C d It is the Reynolds number The function is determined according to the Schiller-Naumann formula: In the formula, Let be the aerodynamic viscosity, Pa·s. The drag coefficient is taken as follows, depending on the Reynolds number range: when <1 o'clock (Stokes District): C d =24 / When 1≤ <1000 hours (transition zone): when When the value is ≥1000 (Newtonian zone): Cd=0.44 For the common water droplet diameter (20-100μm) and wind speed (3-15m / s) conditions of icing on transmission lines, the Reynolds number is usually in the range of 1-100, which is in the transition zone.

[0041] In an optional embodiment, the drag coefficient Cd is calculated using the Oseen approximation formula for the Reynolds number in the transition region (Rep is approximately in the range of 1-100).

[0042] In an alternative embodiment, the air viscous drag coefficient can also be determined using a single continuous function expression applicable to a wider range of Reynolds numbers. The drag coefficient Cd can be solved directly by substituting the calculated droplet Reynolds number (Rep) into the correlation, thus avoiding the use of piecewise functions.

[0043] (3) Electric field polarization force F dep In a non-uniform electric field, even electrically neutral polar molecules (such as water) will be polarized, generating an induced dipole moment, and thus experiencing a dielectric force pointing towards the region of high field strength. This is one of the key electric field effects considered in this invention.

[0044] The electric field polarization force is expressed as: In the formula, The vacuum permittivity, F / m; K is the Clausius-Mossotti factor, dimensionless; E is the electric field intensity vector, V / m; ∇|E|² is the gradient of the square of the electric field intensity, V² / m³.

[0045] The Clausius-Mosotti factor reflects the difference in dielectric constant between the water droplet and the surrounding medium (air), and is calculated using the following formula: In the formula, is the relative permittivity of water, which is approximately 80 at room temperature; Let K be the relative permittivity of air, which is approximately 1. Substituting this into the calculation, we get K ≈ 0.988, which is close to 1.

[0046] The electric polarization force is proportional to the square of the electric field strength and points in the direction of increasing electric field strength. For power transmission lines, the surface field strength is the highest, so the polarization force points towards the surface of the line, increasing the probability of water droplets colliding with the line.

[0047] (4) Corona wind force F c When the electric field strength on the surface of a conductor exceeds the corona initiation field strength of air (approximately 30 kV / cm peak value), corona discharge occurs. Corona discharge generates a large number of ions, which move outward under the influence of the electric field, frequently colliding with air molecules and transferring momentum to them, forming an ion wind blowing outward from the surface of the conductor, i.e., corona wind.

[0048] The corona wind alters the airflow distribution around the conductor, exerting an additional aerodynamic force on the water droplets. The corona wind force can be represented in a form similar to air resistance. In the formula, v c The corona wind velocity is denoted in m / s, and its distribution can be obtained through numerical simulation or experimental measurement of corona discharge.

[0049] The velocity distribution of corona wind can be approximated by an exponential decay model: In the formula, The maximum corona wind velocity on the surface of the conductor is given in m / s. The equivalent radius of the icy conductor is in meters. The characteristic length of corona wind attenuation is in meters (m), typically 5-15 mm.

[0050] (5) Coulomb force Fq After a water droplet gains an electric charge in a high-voltage electric field, the charged water droplet will experience a Coulomb force in the electric field: In the formula, q is the charge of the water droplet, C, which is determined by the charge model; E is the electric field strength, V / m.

[0051] The direction of the Coulomb force depends on the sign of the water droplet's charge and the direction of the electric field. For a power transmission line, the electric field radiates outward from the line. If the water droplet is positively charged, the Coulomb force is directed away from the line; if the water droplet is negatively charged, the Coulomb force is directed towards the line.

[0052] Furthermore, in this embodiment, step S3 establishes the water droplet motion equation based on the aforementioned water droplet force model and water droplet charge model, and uses numerical methods to solve for the motion trajectory of the water droplet at different initial positions. Specific steps include: Based on Newton's second law, the equation of motion of the water droplet is established: In the formula, m p Let v be the mass of the water droplet, in kg; p x is the velocity vector of the water droplet, in m / s; p Let be the position vector of the water droplet, m; F be the net force acting on the water droplet, N; and t be the time, s.

[0053] Expanding the above vector equation into two-dimensional coordinate component form (taking the xy plane as an example): In the formula, (u p , v p (x) represents the components of the water droplet velocity in the x and y directions, in m / s; (x) p , y p (F) represents the position coordinates of the water droplet, in meters; x , F y () represents the components of the resultant force in the x and y directions, in N.

[0054] The above system of ordinary differential equations is solved numerically using the fourth-order Runge-Kutta method. Let the state vector be... The iteration format is then: in: In the formula, Δt is the time step, in seconds; f(·) is the right-hand side function of the state equation: The selection of the time step must meet the requirements of computational stability and accuracy. For typical transmission line icing calculation conditions, the recommended time step is... ,s.

[0055] In an optional embodiment, the system of ordinary differential equations governing the motion of the water droplet can also be solved using the improved Euler method, which iteratively solves the equations step by step. Within each time step, an explicit Euler prediction is first performed, followed by a correction calculation using the prediction results.

[0056] In another alternative embodiment, the system of ordinary differential equations governing the motion of the water droplet can also be solved numerically using the forward Euler method. The state at the next time step is updated using the state and derivative of the current time step through simple linear extrapolation.

[0057] Water droplet release point settings: A water droplet release surface is installed at a certain distance L0 upstream of the guide wire, perpendicular to the incoming flow direction. The height of the release surface ranges from [-H, H], and it must cover the projected area of ​​the icy guide wire with a certain margin. Water droplets are placed on the release surface at preset intervals. N release points are evenly distributed.

[0058] The initial conditions for the water droplet are set as follows: Initial position: ( ) = (-L0, ), i = 1, 2, ..., N Initial velocity: (u0, v0) = (U∞, 0) That is, the initial velocity of the water droplet is equal to the incoming wind speed, and the direction is the same as the wind direction.

[0059] Furthermore, in this embodiment, step S4 establishes a geometric model of the ice-covered guide wire, determines the collision between the trajectory of each water droplet and the ice-covered surface, and records the collision position. Specific steps include: (1) Establish the geometric model of the icing conductor Establish a rectangular coordinate system with the axis of the conductor as the origin, with the x-axis along the direction of the incoming flow and the y-axis vertically upward.

[0060] The conductor has a circular cross-section, and the equation is: An ice layer covers the surface of the conductor. Assuming the ice thickness is uniform, the equation of the outer boundary after icing is: In the formula, δ is the equivalent radius of the icing conductor, in meters; δ is the ice thickness, in meters.

[0061] For rime ice accumulation, icicles will form below the conductor. These icicles are simplified to a cone shape, with the following geometric parameters: Location of the icicle root: (0, -r c ) The tip of the icicle: The radius of the ice crystal varies linearly along its length: In the formula, Let be the length of the icicle, in meters (m). Let be the radius of the icicle's root, in meters (m).

[0062] (2) Collision determination method After each calculation time step, determine the current position (x) of the water droplet. p ,y p Whether it has entered the geometric boundary of the icy guide or ice floes.

[0063] Criteria for determining collision with icing conductors: The criteria for determining collision with icicles (for the icicle region) ): If any of the above conditions are met, the water droplet is determined to be captured, and the collision position coordinates (x, y) are recorded. c ,y c And determine the area (windward, side, leeward or ice) based on the location of the collision.

[0064] If the trajectory of the water droplet bypasses the icy conductor (x>0 and |y|>r) c If the water droplet reaches the boundary of the computational domain (x>xend) but still does not intersect with the icy surface, it is determined that the water droplet has not been captured.

[0065] In an optional embodiment, the collision determination method may also consider the movement line segment of the water droplet from the starting point to the ending point over the entire time step ∆t. Collision determination is performed by calculating whether this movement line segment intersects with the circular boundary of the icing guide or the conical surface of the ice floe. If an intersection occurs, the intersection point is recorded as the collision location.

[0066] In another optional embodiment, the collision determination method can also predefine several key "detection surfaces" around the geometric model of the icy guide wire and icicle. Only when the trajectory of the water droplet passes through these preset detection surfaces will a fine judgment be triggered on the current position and geometric boundary of the water droplet to determine whether a collision has occurred.

[0067] (3) Regional division To analyze the differences in water droplet capture ability at different parts of the iced conductor, the iced surface was divided into the following regions: Region 1 – Windward side: Angle θ∈[-45°,45°], θ=atan2(y,x); Region 2 – Lateral region: Angle θ∈(45°,135°) or θ∈(-135°,-45°); Region 3 – Leeward side: Angle θ∈[135°,180°] or θ∈[-180°,-135°]; Area 4 – Ice Zone: <-r c .

[0068] Furthermore, in this embodiment of the application, step S5 calculates the local capture rate and overall capture rate of each region of the icing conductor based on the collision statistics results. The specific steps include: (1) Local capture rate The local capture rate reflects the ability of different regions on the surface of an icing conductor to capture water droplets, and is defined as the ratio of the number of water droplets captured in that region to the number of water droplets in the incoming flow within the corresponding projected area. In the formula, N ci N represents the number of water droplets captured in the i-th region. 0i The total number of water droplets within the projection range of the i-th region in the incoming flow; The local capture rate can also be defined in a continuous form, based on the ratio of the initial spacing between adjacent trajectories to the spacing between collision points: In the formula, Δy0 is the initial distance between two adjacent trajectories captured by the i-th region on the release surface, in meters (m). Let m be the arc length of the point of collision between the two trajectories on the icy surface.

[0069] (2) Overall capture rate The overall capture rate β reflects the total capture capacity of the icing guide for incoming water droplets, and is defined as the ratio of the total number of water droplets captured to the total number of water droplets in the incoming flow that pass through the projected area of ​​the icing guide: In the formula, N0 is the total number of water droplets in the incoming flow corresponding to the entire projection range of the icing guide; n is the total number of regions divided.

[0070] Quality-weighted calculation of capture rate If the droplet size distribution in the incoming stream is not uniform, the capture rate can be calculated using a mass-weighted method. In the formula, M represents the total number of water droplets released; m j Let χ be the mass of the j-th water droplet, in kg; j Let χ be the collision indicator function; if the j-th droplet is captured, then... j =1, otherwise χ j =0.

[0071] Example 3 is the third embodiment of the present invention, which differs from the previous two embodiments in that: This embodiment also provides an electronic device, including: a memory and a processor; the memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions to implement the method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect as proposed in the above embodiment.

[0072] This embodiment also provides a storage medium storing a computer program that, when executed by a processor, implements a method for calculating the capture rate of ice-covered water droplets on a wire, considering the electric field effect, as proposed in the above embodiment.

[0073] The storage medium proposed in this embodiment and the method for calculating the capture rate of iced water droplets on wires considering the electric field effect proposed in the above embodiment belong to the same inventive concept. Technical details not described in detail in this embodiment can be found in the above embodiments, and this embodiment has the same beneficial effects as the above embodiments.

[0074] Based on the above description of the implementation methods, those skilled in the art can clearly understand that the present invention can be implemented using software and necessary general-purpose hardware, and of course, it can also be implemented using hardware. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as a computer floppy disk, read-only memory (ROM), random access memory (RAM), flash memory, hard disk, or optical disk, etc., including several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods of the various embodiments of the present invention.

[0075] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A method for calculating the water droplet capture rate of an ice-covered wire considering the electric field effect, characterized in that: include, Obtain conductor parameters, environmental parameters, and electrical parameters; determine the charge of the water droplet based on the droplet diameter and electric field strength; and establish a water droplet charge model. Establish a force model of the water droplet and calculate the net force acting on the water droplet during its motion; Based on the aforementioned force model and charge model of the water droplet, the motion equation of the water droplet is established, and the motion trajectory of the water droplet at different initial positions is solved using numerical methods. Establish a geometric model of the ice-covered guide wire, determine the collision between the trajectory of each water droplet and the ice-covered surface, and record the collision positions; Based on the collision statistics, the local capture rate and overall capture rate of each region of the icy conductor are calculated.

2. The method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect as described in claim 1, characterized in that: The establishment of the water droplet charge model includes: experimentally measuring the charge of water droplets of different diameters under different electric field intensities, establishing a model relating the water droplet charge to the droplet diameter and electric field intensity, and using a fitting method to determine the model parameters.

3. The method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect as described in claim 2, characterized in that: The establishment of the force model of the water droplet includes vector superposition of gravity, air viscous drag, electric field polarization force, corona wind force and Coulomb force to obtain the resultant force on the water droplet during its motion, wherein the air viscous drag coefficient is determined piecewise based on the water droplet Reynolds number.

4. The method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect as described in claim 3, characterized in that: The method for solving the water droplet trajectory includes setting a water droplet release surface upstream of the guide wire, with the release surface perpendicular to the incoming flow direction; Water droplet release points are evenly arranged on the release surface at preset intervals; Based on the median diameter and particle size distribution of water droplets, water droplets of different diameters are allocated to each release point; Calculate the trajectory of the water droplets released at each release point under the action of the resultant force.

5. The method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect as described in claim 4, characterized in that: The establishment of the geometric model of the icy conductor includes establishing a geometric model of the cylindrical surface of the conductor based on the conductor diameter; Based on the type and thickness of icing, a geometric model of the icing layer is established on the surface of the conductor; When it is determined to be rime ice, an ice geometric model is established below the ice layer. The ice geometric model includes parameters such as ice length, ice diameter, and ice spacing.

6. The method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect as described in claim 5, characterized in that: The determination of the collision between the trajectory of each water droplet and the ice-covered surface includes determining whether the position of the water droplet enters the boundary of the geometric model of the ice-covered conductor within each calculation time step. When a water droplet enters the boundary of the ice-covered surface, it is determined that the current water droplet has been captured, and the collision position coordinates are recorded. If the trajectory of a water droplet bypasses the icing wire or reaches the boundary of the computational domain without intersecting the icing surface, it is determined that the current water droplet has not been captured.

7. The method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect as described in claim 6, characterized in that: The calculation of the local capture rate and overall capture rate of each region of the icy conductor includes dividing the surface of the icy conductor into several regions, including the windward region of the conductor, the leeward region of the conductor, the surface region of the icing layer, the tip region of the ice floe, and the side region of the ice floe. Local capture rate β in each region i The calculation formula is: In the formula, N ci N represents the number of water droplets captured in the i-th region. 0i The total number of water droplets within the projection range of the i-th region in the incoming flow; The formula for calculating the overall capture rate β is: In the formula, N0 is the total number of water droplets in the incoming flow corresponding to the entire projection range of the icing guide, and n is the total number of regions divided.

8. The method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect as described in claim 7, characterized in that, The calculated local and overall capture rates were compared and verified with the results of artificial climate chamber experiments. Based on the comparison results, the parameters in the water droplet force model or water droplet charge model were corrected to improve the calculation accuracy.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect, as described in any one of claims 1 to 8.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for calculating the ice-covered water droplet capture rate of a wire considering the electric field effect, as described in any one of claims 1 to 8.

Images