Power transmission line icing analogue simulation method, device and equipment and storage medium

Abstract

The invention provides a power transmission line icing analogue simulation method, device and equipment and a storage medium, and relates to the technical field of freezing meteorological simulation. The method comprises the steps of creating a target three-dimensional model corresponding to a power transmission line in response to an analogue simulation instruction for icing of the power transmission line; performing grid division on the target three-dimensional model to obtain a computational domain which corresponds to the power transmission line and comprises a plurality of infinitesimal elements; simulating an air flow field in a computational domain based on a computational fluid mechanics method by taking a wind speed and a wind direction in the actually measured meteorological data as boundary conditions to obtain first attribute information of the air flow field; acquiring second attribute information of water drops in the air flow field based on the first attribute information by taking precipitation in the actually measured meteorological data as a boundary condition; and by taking the temperature in the actually measured meteorological data as a boundary condition, simulating the cooling and solidification process of the water drop based on the second attribute information, and obtaining the icing thickness and the icing shape after the set time step length. According to the invention, the icing process and result of the power transmission line can be simulated more accurately.

Description

Technical Field

[0001] The present application relates to the technical field of freezing weather simulation, and in particular to a method, device, equipment and storage medium for simulating icing on transmission lines. Background Art

[0002] In recent years, with climate change, the frequency and intensity of abnormal climate fluctuations have gradually increased. Frequent icing is gradually becoming a major hidden danger to the safe operation of power systems. Therefore, accurately simulating icing on transmission lines has become a key issue in power system management.

[0003] Currently, simulations of transmission line icing are typically performed using finite element analysis under fixed, single meteorological conditions. However, this approach often results in significant discrepancies between the simulation results and actual icing results. Summary of the Invention

[0004] The present application provides a method, device, equipment and medium for simulating icing on a transmission line to solve the problem that when simulating icing on a transmission line using the current method, there is a large deviation between the simulation results and the actual icing results.

[0005] In a first aspect, the present application provides a method for simulating icing on a transmission line, comprising:

[0006] In response to a simulation instruction for ice coating on a transmission line, a target three-dimensional model corresponding to the transmission line is created, wherein the target three-dimensional model includes the transmission line and an external flow field around the transmission line;

[0007] Mesh the target three-dimensional model to obtain a computational domain containing multiple microelements;

[0008] Using wind speed and direction from measured meteorological data as boundary conditions, the air flow field within the computational domain is simulated using computational fluid dynamics methods to obtain primary attribute information of the air flow field, including airflow velocity, airflow direction, and turbulence structure around the transmission line.

[0009] Using the precipitation in the measured meteorological data as a boundary condition, the second attribute information of water droplets in the air flow field is obtained based on the first attribute information. The second attribute information includes the movement path and speed of the water droplets under the action of the air flow field, as well as the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets.

[0010] Taking the temperature in the measured meteorological data as the boundary condition, the cooling and solidification process of water droplets is simulated based on the second attribute information to obtain the ice thickness and shape on the surface of the transmission line after a set time step.

[0011] Optionally, the air flow field in the calculation domain is simulated based on a computational fluid dynamics method to obtain first attribute information of the air flow field, including: simulating the air flow field in the calculation domain based on a computational fluid dynamics method, using the Navier-Stokes equations and the continuity equation to obtain the air flow velocity and air flow direction, and obtaining the turbulence structure through a preset turbulence model.

[0012] Optionally, obtaining second attribute information of water droplets in the air flow field based on the first attribute information includes: obtaining the second attribute information according to the following first formula:

[0013]

[0014] Where α represents the volume fraction of the water droplet phase; Indicates air speed; Indicates the water droplet velocity; C D Represents the water drop resistance coefficient; d Reynolds number for the relative motion of water droplets; ρ a represents the density of air; ρ d represents the density of the water droplet; Fr represents the local Froude number; K represents the inertial parameter of the water droplet; V a,∞ represents the free-flow air velocity; D represents the diameter of the transmission line; L ∞ represents characteristic length; μ represents air kinematic viscosity; M represents precipitation; and v represents the falling velocity of water droplets.

[0015] Optionally, the cooling and solidification process of water droplets is simulated based on the second attribute information to obtain the ice thickness and ice shape on the surface of the transmission line after the set time step, including: simulating the cooling and solidification process of water droplets based on the second attribute information to obtain the ice mass and ice density; based on the ice mass and ice density, obtaining the ice volume and ice thickness on the surface of the transmission line after the set time step; based on the ice volume, obtaining the ice shape on the surface of the transmission line after the set time step.

[0016] Optionally, simulating the cooling and solidification process of the water droplet based on the second attribute information to obtain the ice coverage mass includes: obtaining the ice coverage mass according to the following second formula:

[0017]

[0018] Among them, h f Indicates the height of the water film; represents the velocity of the water film; ρ f Indicates the density of the water film; V ∞ represents the free-flow air velocity; a represents the volume fraction of the water droplet phase; It represents the mass of water droplets evaporated from the microelement; represents the mass of frozen water droplets in the microelement; c frepresents the specific heat capacity of water; c s represents the specific heat capacity of ice; represents the equilibrium temperature of the ice-covered interface; Indicates the free-flow air temperature; represents the average temperature of the ice layer; represents the ambient temperature after taking into account the energy loss caused by friction during convection; represents the local collision velocity of water droplets; L evap represents the latent heat of evaporation of water; L fusion represents the latent heat of melting of ice; σ represents the blackbody radiation constant; ε represents the blackness; A represents the radiated area of the microelement; h c Indicates the convection heat transfer coefficient between the external environment and the transmission line surface; Q anti-icing represents the anti-icing heat flux; β represents the local collision coefficient.

[0019] Optionally, simulating the cooling and solidification process of the water droplet based on the second attribute information to obtain the ice density includes: obtaining the ice density according to the following third formula:

[0020]

[0021] Among them, R M represents the calculation parameter of ice density of transmission line; d represents the diameter of supercooled water droplet; represents the collision velocity of supercooled water droplets; Indicates the ice surface temperature.

[0022] Optionally, meshing the target three-dimensional model to obtain a computational domain including multiple micro-elements includes: meshing the target three-dimensional model in 2.5 dimensions to obtain a computational domain including multiple micro-elements.

[0023] Optionally, after obtaining the ice thickness and ice shape on the surface of the transmission line after the set time step, the transmission line ice simulation method also includes: reconstructing the ice boundary corresponding to the transmission line to obtain a three-dimensional model of the reconstructed ice boundary; using the reconstructed three-dimensional model of the ice boundary as a new target three-dimensional model, performing meshing on the target three-dimensional model, and obtaining a calculation domain containing multiple microelements to obtain the ice thickness and ice shape on the surface of the transmission line after the next set time step.

[0024] In a second aspect, the present application provides a transmission line icing simulation device, comprising:

[0025] A creation module is configured to create a target three-dimensional model corresponding to the transmission line in response to a simulation instruction for icing of the transmission line, wherein the target three-dimensional model includes the transmission line and an external flow field around the transmission line;

[0026] A partitioning module is used to perform mesh partitioning on the target three-dimensional model to obtain a computational domain containing multiple microelements;

[0027] a first acquisition module for simulating the air flow field in the computational domain based on a computational fluid dynamics method using wind speed and wind direction in the measured meteorological data as boundary conditions, thereby obtaining first attribute information of the air flow field, the first attribute information including air flow velocity, air flow direction, and turbulence structure around the transmission line;

[0028] a second acquisition module for acquiring second attribute information of water droplets in the air flow field based on the first attribute information and using the precipitation in the measured meteorological data as a boundary condition, wherein the second attribute information includes the movement path and speed of the water droplets under the action of the air flow field, as well as the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets;

[0029] The third acquisition module is used to simulate the cooling and solidification process of water droplets based on the second attribute information using the temperature in the measured meteorological data as the boundary condition, and obtain the ice thickness and ice shape on the surface of the transmission line after a set time step.

[0030] Optionally, the first acquisition module is specifically used to: simulate the air flow field in the calculation domain based on the computational fluid dynamics method, obtain the airflow velocity and direction using the Navier-Stokes equations and the continuity equation, and obtain the turbulence structure through a preset turbulence model.

[0031] Optionally, the second acquisition module is specifically configured to: acquire the second attribute information according to the following first formula:

[0032]

[0033] Where α represents the volume fraction of the water droplet phase; Indicates air speed; Indicates the water droplet velocity; C D Represents the water drop resistance coefficient; d Reynolds number for the relative motion of water droplets; ρ a represents the density of air; ρ d represents the density of the water droplet; Fr represents the local Froude number; K represents the inertial parameter of the water droplet; V a,∞ represents the free-flow air velocity; D represents the diameter of the transmission line; L ∞ represents characteristic length; μ represents air kinematic viscosity; M represents precipitation; and v represents the falling velocity of water droplets.

[0034] Optionally, the third acquisition module is specifically used to: simulate the cooling and solidification process of water droplets based on the second attribute information to obtain the ice mass and ice density; based on the ice mass and ice density, obtain the ice volume and ice thickness on the surface of the transmission line after a set time step; based on the ice volume, obtain the ice shape on the surface of the transmission line after a set time step.

[0035] Optionally, when the third acquisition module is used to simulate the cooling and solidification process of the water droplet based on the second attribute information to obtain the ice coverage mass, it is specifically used to: obtain the ice coverage mass according to the following second formula:

[0036]

[0037] Among them, h f Indicates the height of the water film; represents the velocity of the water film; ρ f Indicates the density of the water film; V ∞ represents the free-flow air velocity; a represents the volume fraction of the water droplet phase; It represents the mass of water droplets evaporated from the microelement; represents the mass of frozen water droplets in the microelement; c f represents the specific heat capacity of water; c s represents the specific heat capacity of ice; represents the equilibrium temperature of the ice-covered interface; Indicates the free-flow air temperature; represents the average temperature of the ice layer; represents the ambient temperature after taking into account the energy loss caused by friction during convection; represents the local collision velocity of water droplets; L evap represents the latent heat of evaporation of water; L fusion represents the latent heat of melting of ice; σ represents the blackbody radiation constant; ε represents the blackness; A represents the radiated area of the microelement; h c Indicates the convection heat transfer coefficient between the external environment and the transmission line surface; Q anti-icing represents the anti-icing heat flux; β represents the local collision coefficient.

[0038] Optionally, when the third acquisition module is used to simulate the cooling and solidification process of water droplets based on the second attribute information to obtain the ice density, it is specifically used to: obtain the ice density according to the following third formula:

[0039]

[0040] Among them, R M represents the calculation parameter of ice density of transmission line; d represents the diameter of supercooled water droplet; represents the collision velocity of supercooled water droplets; Indicates the ice surface temperature.

[0041] Optionally, the partitioning module is specifically used to: perform 2.5-dimensional grid partitioning on the target three-dimensional model to obtain a calculation domain containing multiple microelements.

[0042] Optionally, the transmission line icing simulation device also includes a processing module, which is used to: after obtaining the ice thickness and ice shape on the transmission line surface after a set time step, reconstruct the ice boundary corresponding to the transmission line to obtain a three-dimensional model of the reconstructed ice boundary; use the three-dimensional model of the reconstructed ice boundary as a new target three-dimensional model, perform meshing on the target three-dimensional model, and obtain a calculation domain containing multiple infinitesimals to obtain the ice thickness and ice shape on the transmission line surface after the next set time step.

[0043] In a third aspect, the present application provides an electronic device, comprising: a processor, and a memory communicatively connected to the processor;

[0044] Memory stores computer-executable instructions;

[0045] The processor executes the computer-executable instructions stored in the memory to implement the transmission line icing simulation method as described in the first aspect of the present application.

[0046] In a fourth aspect, the present application provides a computer-readable storage medium, in which computer program instructions are stored. When the computer program instructions are executed, the transmission line icing simulation method as described in the first aspect of the present application is implemented.

[0047] In a fifth aspect, the present application provides a computer program product, including a computer program, which, when executed, implements the transmission line icing simulation method as described in the first aspect of the present application.

[0048] The present application provides a method, apparatus, equipment and storage medium for simulating icing on a transmission line. In response to a simulation instruction for icing on a transmission line, a target three-dimensional model corresponding to the transmission line is created, the target three-dimensional model including the transmission line and the external flow field around the transmission line; the target three-dimensional model is meshed to obtain a computational domain containing multiple microelements; the wind speed and wind direction in the measured meteorological data are used as boundary conditions, the air flow field in the computational domain is simulated based on a computational fluid dynamics method, and first attribute information of the air flow field is obtained, the first attribute information including the airflow velocity, airflow direction and turbulence structure around the transmission line; the precipitation in the measured meteorological data is used as a boundary condition, and second attribute information of water droplets in the air flow field is obtained based on the first attribute information, the second attribute information including the movement path and speed of the water droplets under the action of the air flow field, and the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets; the temperature in the measured meteorological data is used as a boundary condition, and the cooling and solidification process of the water droplets is simulated based on the second attribute information to obtain the thickness and shape of the ice on the surface of the transmission line after a set time step. This application simulates the ice growth process of transmission lines in a dynamic time-varying environment based on measured meteorological data and finite element analysis methods, rather than simulating the icing of transmission lines based on finite element analysis methods under fixed single meteorological conditions. Therefore, this application can more accurately simulate the icing process and results of transmission lines in actual environments, and effectively solve the problem of large deviations between the simulation results and actual icing results of transmission lines under complex time-varying meteorological conditions. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0050] Figure 1 This is a flow chart of a method for simulating icing on a transmission line provided in one embodiment of the present application;

[0051] Figure 2 A flowchart of a method for simulating icing on a transmission line provided in another embodiment of the present application;

[0052] Figure 3 A schematic diagram of the xy plane meshing result of the target three-dimensional model provided in one embodiment of the present application;

[0053] Figure 4 A schematic diagram of a transmission line wall provided in one embodiment of the present application;

[0054] Figure 5 A flowchart of a method for simulating icing on a transmission line according to another embodiment of the present application;

[0055] Figure 6A schematic diagram of the structure of a transmission line icing simulation device provided in one embodiment of the present application;

[0056] Figure 7 A schematic diagram of the structure of an electronic device provided in one embodiment of the present application.

[0057] The above drawings illustrate specific embodiments of the present application, which will be described in more detail below. These drawings and the textual description are not intended to limit the scope of the present application in any way, but rather to illustrate the concepts of the present application to those skilled in the art by reference to specific embodiments. DETAILED DESCRIPTION

[0058] To make the purpose, technical solutions, and advantages of the embodiments of this application more clear, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0059] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant laws, regulations and standards, and provide corresponding operation entrances for users to choose to authorize or refuse.

[0060] In recent years, with climate change, the frequency and intensity of abnormal climate fluctuations have gradually increased, and the frequent occurrence of icing is gradually becoming a major hidden danger to the safe operation of power systems. Therefore, accurately simulating the icing of transmission lines has become a key issue in power system management. For the study of the icing process and mechanism of transmission line icing under time-varying meteorological factors, the main research methods include: on-site icing observation and measurement, artificial icing experiments, and numerical simulation. Compared with on-site icing observation and artificial icing experiments, numerical simulation based on fluid dynamics and icing growth mechanisms provides a more convenient way to study the thickness and shape of icing on transmission lines.

[0061] Currently, simulations of transmission line icing are typically performed using finite element analysis (FEM) under fixed, single meteorological conditions. These methods fail to consider the impact of real-time micrometeorological factors on transmission line icing. Furthermore, the settings for meteorological parameters such as wind speed, direction, temperature, droplet volume fraction, and droplet diameter are not precise enough. Consequently, while these methods perform well in hypothetical, steady-state meteorological environments, they can produce significant simulation errors when dealing with the complex, time-varying icing conditions found in real-world projects.

[0062] Based on the above problems, the present application provides a method, device, equipment and storage medium for simulating icing on transmission lines, which performs simulation of icing on transmission lines under measured meteorological data based on the icing mechanism and multiphase flow model of transmission lines; wherein, within each set time step, the simulation of air flow field, water droplet movement and collision collection simulation, water droplet freezing and icing growth calculation are carried out in sequence, and simulation results of icing thickness and icing shape that are more in line with reality can be obtained.

[0063] It should be noted that the transmission line icing simulation method provided in the embodiment of the present application can be applied in a server, which can be an independent server or a service cluster.

[0064] The following specific embodiments describe in detail the technical solution of the present application and how the technical solution of the present application solves the above-mentioned technical problems. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments. The embodiments of the present application will be described below in conjunction with the accompanying drawings.

[0065] Figure 1 This is a flow chart of a method for simulating ice coating on a transmission line provided in one embodiment of the present application. Figure 1 As shown, the transmission line icing simulation method of the embodiment of the present application includes:

[0066] S101. In response to a simulation instruction for icing on a transmission line, create a target three-dimensional model corresponding to the transmission line, where the target three-dimensional model includes the transmission line and an external flow field around the transmission line.

[0067] In an embodiment of the present application, the wall of the transmission line and the air around the transmission line constitute the calculation domain corresponding to the transmission line. For example, the condition parameters required for the simulation can be obtained based on the measured meteorological data such as wind speed, wind direction, temperature, humidity and precipitation near the transmission line in the area to be simulated, and then a target three-dimensional model corresponding to the transmission line is created. The target three-dimensional model includes the transmission line and the external flow field around the transmission line. Specifically, the transmission line and the external flow field around the transmission line can be three-dimensionally modeled using engineering simulation software to obtain a target three-dimensional model. For example, the SpaceClaim module (a modeling tool) of ANSYS2024R2 (an engineering simulation software) can be used to create a target three-dimensional model corresponding to the transmission line. Among them, the length of the calculation domain corresponding to the transmission line is, for example, 2 meters (m), and the width of the calculation domain is, for example, 2m, to meet the full development of the wake and vortex shedding. The model of the simulated transmission line is, for example, JLHA2 / G3A-250, and the outer diameter is, for example, 23.8 millimeters (mm).

[0068] S102: Meshing the target three-dimensional model to obtain a calculation domain containing multiple microelements.

[0069] For example, the ICEM CFD module (a meshing tool for computational fluid dynamics simulation) of ANSYS 2024R2 can be used to mesh the target 3D model to obtain a computational domain containing multiple microelements. Specifically, for example, a 2.5-dimensional (2.5D) mesh can be performed on the target 3D model to obtain a computational domain containing multiple microelements.

[0070] S103. Using the wind speed and wind direction in the measured meteorological data as boundary conditions, the air flow field in the calculation domain is simulated based on the computational fluid dynamics method to obtain first attribute information of the air flow field. The first attribute information includes the air flow velocity, air flow direction and turbulence structure around the transmission line.

[0071] It can be understood that based on the principles of computational fluid dynamics, the air flow containing water droplets can be regarded as a gas-liquid two-phase flow. Since the water droplet particle size is small, the influence of the water droplet movement on the air flow field can generally be ignored. Therefore, the air flow field around the transmission line can be first solved through step S103, and then the motion trajectory of the water droplets in the air flow field obtained in step S103 can be analyzed through step S104. For example, with the wind speed and wind direction in the measured meteorological data as boundary conditions, the air flow field in the calculation domain is simulated based on the computational fluid dynamics method, including simulating key parameters such as the air flow velocity, direction and turbulent structure around the transmission line. Accordingly, the first attribute information of the air flow field can be obtained. The first attribute information includes the air flow velocity, air flow direction and turbulent structure around the transmission line, thereby providing accurate air flow field data for the subsequent water droplet collision trajectory simulation. For specific details on how to simulate the air flow field in the calculation domain based on the computational fluid dynamics method and obtain the first attribute information of the air flow field, please refer to the subsequent embodiments and will not be repeated here.

[0072] S104. Using the precipitation in the measured meteorological data as a boundary condition, obtain the second attribute information of the water droplets in the air flow field based on the first attribute information. The second attribute information includes the movement path and speed of the water droplets under the action of the air flow field, as well as the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets.

[0073] In this step, after obtaining the first attribute information of the air flow field, taking the precipitation in the measured meteorological data as the boundary condition, the movement path, speed and collision behavior of water droplets in the air flow field and the surface of the transmission line can be obtained based on the first attribute information, and considering physical phenomena such as evaporation, breakage and adhesion, the water droplet collection amount and spatial distribution of water droplets on the surface of the transmission line can be obtained.

[0074] Further, optionally, obtaining second attribute information of water droplets in the air flow field based on the first attribute information may include: obtaining the second attribute information according to the following formula 1 (i.e., the first formula):

[0075]

[0076] Where α represents the volume fraction of the water droplet phase; Indicates air speed; Indicates the water droplet velocity; C D Represents the water drop resistance coefficient; d Reynolds number for the relative motion of water droplets; ρ a represents the density of air; ρ d represents the density of the water droplet; Fr represents the local Froude number; K represents the inertial parameter of the water droplet; V a,∞ represents the free-flow air velocity; D represents the diameter of the transmission line; L ∞ represents characteristic length; μ represents air kinematic viscosity; M represents precipitation; and v represents the falling velocity of water droplets.

[0077] For example, the Fluent Icing finite element analysis module (a module for simulating and analyzing icing phenomena) of ANSYS 2024R2 can be used to obtain the collection process of water droplets. According to the Fluent Icing user manual, the relevant calculation equation is as shown in Formula 1 above, where α is the volume fraction of the water droplet phase, which represents the ratio of the water droplet volume to the total volume, and the unit is, for example, kilograms per cubic meter (kg / m 3 ), the water droplet phase volume fraction α determines the collection efficiency of the colliding water droplets on the transmission line surface. In the embodiment of the present application, the water droplet phase volume fraction α can be calculated from the precipitation in the weather forecast, which helps to obtain more accurate ice thickness simulation results; L ∞ is a characteristic length, such as the radius of the transmission line. Formula 1 above can be used to determine the direction and velocity of water droplets around the transmission line, the amount of water droplets collected on the surface of the transmission line, and the spatial distribution of water droplets at the impact points on the transmission line, enabling a finite element numerical solution.

[0078] S105. Using the temperature in the measured meteorological data as a boundary condition, simulate the cooling and solidification process of the water droplets based on the second attribute information to obtain the ice thickness and ice shape on the surface of the transmission line after a set time step.

[0079] Exemplarily, the set time step can be determined based on the transmission line icing rate and meteorological conditions. The set time step is, for example, 20 seconds, which is not limited in the embodiments of the present application. It can be understood that the measured meteorological data within the set time step maintains its mean value unchanged. After obtaining the second attribute information of water droplets in the air flow field, considering key factors such as conservation of energy, conservation of mass, and latent heat of evaporation and melting, with the temperature in the measured meteorological data as the boundary condition, the freezing process, ice thickness and shape of the water droplets collected on the transmission line surface are simulated based on the amount of water droplets collected on the transmission line surface and the spatial distribution of water droplets, thereby obtaining the ice thickness and ice shape on the transmission line surface after the set time step. For example, ANSYS2024R2 software can be used to simulate the growth of ice on transmission lines with multiple set time steps. Based on the above steps, the air flow field simulation, water droplet movement and collision collection simulation, water droplet freezing and ice growth calculation are carried out in sequence within each set time step; the ice boundary based on the re-division of the grid is automatically updated, and the ice thickness and ice shape on the transmission line surface after each set time step are iteratively obtained until the preset ice duration is reached. The preset ice duration is, for example, 24 hours, so that the ice thickness and ice shape on the transmission line surface after the preset ice duration can be output.

[0080] The transmission line icing simulation method provided in the embodiment of the present application creates a target three-dimensional model corresponding to the transmission line in response to a simulation instruction for transmission line icing, the target three-dimensional model including the transmission line and the external flow field around the transmission line; meshes the target three-dimensional model to obtain a calculation domain containing multiple microelements; uses the wind speed and wind direction in the measured meteorological data as boundary conditions, simulates the air flow field in the calculation domain based on the computational fluid dynamics method, and obtains first attribute information of the air flow field, the first attribute information including the airflow velocity, airflow direction and turbulence structure around the transmission line; uses the precipitation in the measured meteorological data as a boundary condition, obtains second attribute information of water droplets in the air flow field based on the first attribute information, the second attribute information including the movement path and speed of the water droplets under the action of the air flow field, and the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets; uses the temperature in the measured meteorological data as a boundary condition, simulates the cooling and solidification process of the water droplets based on the second attribute information, and obtains the ice thickness and ice shape on the surface of the transmission line after a set time step. The embodiment of the present application simulates the ice growth process of the transmission line in a dynamic time-varying environment based on measured meteorological data and the finite element analysis method, rather than simulating the icing of the transmission line based on the finite element analysis method under fixed single meteorological conditions. Therefore, the embodiment of the present application can more accurately simulate the icing process and results of the transmission line in the actual environment, and effectively solve the problem of large deviation between the simulation results and the actual icing results of the transmission line under complex time-varying meteorological conditions.

[0081] Figure 2This is a flow chart of a method for simulating ice coating on a transmission line provided by another embodiment of the present application. Based on the above embodiment, this embodiment of the present application further illustrates the method for simulating ice coating on a transmission line. Figure 2 As shown, the transmission line icing simulation method of the embodiment of the present application may include:

[0082] S201. In response to a simulation instruction for icing on a transmission line, create a target three-dimensional model corresponding to the transmission line, where the target three-dimensional model includes the transmission line and an external flow field around the transmission line.

[0083] The detailed description of this step can be found in Figure 1 The relevant description of S101 in the illustrated embodiment will not be repeated here.

[0084] In the embodiment of this application, Figure 1 The step S102 may further include the following step S202:

[0085] S202 , performing 2.5-dimensional grid division on the target three-dimensional model to obtain a calculation domain containing multiple microelements.

[0086] For example, considering the complexity of ice coating on transmission lines and the need for computational efficiency, the ICEM CFD module of ANSYS 2024R2 can be used to perform 2.5-dimensional meshing on the target three-dimensional model to obtain a computational domain containing multiple microelements. Figure 3 A schematic diagram of the xy plane mesh division result of the target three-dimensional model provided in an embodiment of the present application is shown as follows: Figure 3 As shown, a calculation domain 301 containing multiple microelements is obtained after the target three-dimensional model is divided into 2.5-dimensional grids, and the white area in the middle of the calculation domain is the transmission line. Figure 4 A schematic diagram of a transmission line wall provided in one embodiment of the present application, such as Figure 4 The figure shows the transmission line wall. For the area near the transmission line wall, in order to more accurately capture the parameter changes near the transmission line wall, the near-wall grid can be encrypted (i.e., the closer to the transmission line, the denser the grid), with the first layer grid height being, for example, 4e-6 meters (m).

[0087] In the embodiment of this application, Figure 1 The step S103 may further include the following step S203:

[0088] S203. Taking the wind speed and wind direction in the measured meteorological data as boundary conditions, the air flow field in the calculation domain is simulated based on the computational fluid dynamics method, and the Navier-Stokes equations and the continuity equation are used to obtain the airflow velocity and airflow direction in the first attribute information of the air flow field, and the turbulence structure in the first attribute information is obtained through a preset turbulence model.

[0089] It can be understood that based on the principles of computational fluid dynamics, the air flow containing water droplets can be regarded as a gas-liquid two-phase flow. Since the water droplet particle size is small, the influence of the water droplet movement on the air flow field can generally be ignored. Therefore, the air flow field around the transmission line can be solved first in this step, and then the motion trajectory of the water droplets in the air flow field obtained in step S203 can be analyzed in step S204. For example, the air flow velocity and air flow direction of the air flow field in the calculation domain can be obtained using the Navier-Stokes (NS) equation shown in the following formula 2 and the continuity equation shown in the following formula 3:

[0090]

[0091] Among them, u i represents the velocity component of the infinitesimal element in the i-th direction; u j represents the velocity component of the infinitesimal element in the jth direction; t represents time; x i represents the spatial coordinate of the infinitesimal element in the i-th direction; x j represents the spatial coordinate of the infinitesimal element in the jth direction; F i It represents the volume force on unit fluid; ρ represents the air density; p represents the air pressure; μ is the kinematic viscosity of air.

[0092] Considering that the NS equation is a complex second-order linear partial differential equation, it is difficult to solve it directly. In order to obtain the airflow velocity in the calculation domain corresponding to the transmission line, the embodiment of the present application can discretize the NS equation by the finite volume method of FLUENT software (a computational fluid dynamics software). In addition, considering the influence of turbulence on the air flow field during simulation, in order to obtain more accurate air flow field simulation results, the embodiment of the present application uses the k-ωSST (a model for simulating turbulent flow) turbulence model as the preset turbulence model, and obtains the turbulence structure in the first attribute information through the k-ωSST turbulence model.

[0093] S204. Using the precipitation in the measured meteorological data as a boundary condition, obtain the second attribute information of the water droplets in the air flow field based on the first attribute information. The second attribute information includes the movement path and speed of the water droplets under the action of the air flow field, as well as the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets.

[0094] The detailed description of this step can be found in Figure 1 The description related to S104 in the illustrated embodiment will not be repeated here.

[0095] In the embodiment of this application, Figure 1 The step S105 may further include the following three steps S205 to S207:

[0096] S205. Using the temperature in the measured meteorological data as a boundary condition, simulate the cooling and solidification process of the water droplet based on the second attribute information to obtain the ice mass and ice density.

[0097] For example, after obtaining the second property information of water droplets in the air flow field, using the temperature from the measured meteorological data as the boundary condition, the freezing process of water droplets collected on the transmission line surface can be simulated based on the amount of water droplets collected and the spatial distribution of water droplets. The freezing process of water droplets on the transmission line surface is essentially a thermodynamic equilibrium process. This process can be simulated using the mass conservation equation and energy conservation equation shown in Equation 4 below.

[0098] Furthermore, optionally, simulating the cooling and solidification process of the water droplet based on the second attribute information to obtain the ice mass includes: obtaining the ice mass according to the following formula 4 (i.e., the second formula):

[0099]

[0100] Among them, h f Indicates the height of the water film; represents the velocity of the water film; ρ f Indicates the density of the water film; V ∞ represents the free-flow air velocity; a represents the volume fraction of the water droplet phase; It represents the mass of water droplets evaporated from the microelement; represents the mass of frozen water droplets in the microelement; c f represents the specific heat capacity of water; c s represents the specific heat capacity of ice; represents the equilibrium temperature of the ice-covered interface; Indicates the free-flow air temperature; represents the average temperature of the ice layer; represents the ambient temperature after taking into account the energy loss caused by friction during convection; represents the local collision velocity of water droplets; L evap represents the latent heat of evaporation of water; L fusion represents the latent heat of melting of ice; σ represents the blackbody radiation constant; ε represents the blackness; A represents the radiated area of the microelement; h c Indicates the convection heat transfer coefficient between the external environment and the transmission line surface; Q anti-icing represents the anti-icing heat flux; β represents the local collision coefficient.

[0101] For example, β is the local collision coefficient, which represents the ratio of the actual collision rate of water droplets to the maximum possible collision rate. It is an important parameter for obtaining the amount of water droplets collected on the transmission line surface and the spatial distribution of water droplets. The definition of β satisfies the following formula 5:

[0102]

[0103] Where dy represents the difference in coordinates of adjacent incident points in the y-axis direction (vertical direction); dl represents the length of the arc between adjacent collision points on the surface of the transmission line.

[0104] Further, optionally, the cooling and solidification process of the water droplet is simulated based on the second attribute information to obtain the ice density, including: obtaining the ice density ρ according to the following formula 6 (i.e., the third formula): c :

[0105]

[0106] Among them, R M represents the calculation parameter of ice density of transmission line; d represents the diameter of supercooled water droplet; represents the collision velocity of supercooled water droplets; Indicates the ice surface temperature.

[0107] Ice density ρ c It is an important parameter in the prediction of icing on transmission lines. In the embodiment of the present application, for example, the Makkonen-Stallabrass density formula can be used, which takes into account more comprehensive icing influencing parameters.

[0108] S206. Based on the ice mass and ice density, obtain the ice volume and ice thickness on the surface of the transmission line after a set time step.

[0109] In this step, after obtaining the ice mass and ice density, the ice volume of the transmission line surface corresponding to the microelement under the set time step can be obtained based on the ice mass and ice density, and then the ice thickness can be obtained based on the ice volume.

[0110] S207. Based on the ice-covered volume, obtain the ice-covered shape of the transmission line surface after a set time step.

[0111] For example, assuming that the ice growth of each microelement is along the normal direction of the ice surface, the ice shape of the entire transmission line surface after a set time step can be obtained by combining the ice volume growth of all microelements.

[0112] The transmission line icing simulation method provided in the embodiment of the present application creates a target three-dimensional model corresponding to the transmission line in response to a simulation instruction for transmission line icing, the target three-dimensional model including the transmission line and the external flow field around the transmission line; performs 2.5-dimensional grid division on the target three-dimensional model to obtain a calculation domain containing multiple microelements; uses the wind speed and wind direction in the measured meteorological data as boundary conditions, simulates the air flow field in the calculation domain based on the computational fluid dynamics method, uses the Navier-Stokes equation and the continuity equation to obtain the air flow velocity and air flow direction in the first attribute information of the air flow field, and obtains the first attribute information by a preset turbulence model. The turbulent structure in the information is obtained; with the precipitation in the measured meteorological data as the boundary condition, the second attribute information of the water droplets in the air flow field is obtained based on the first attribute information. The second attribute information includes the movement path and speed of the water droplets under the action of the air flow field, as well as the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets; with the temperature in the measured meteorological data as the boundary condition, the cooling and solidification process of the water droplets is simulated based on the second attribute information to obtain the ice mass and ice density; based on the ice mass and ice density, the ice volume and ice thickness on the surface of the transmission line after the set time step are obtained; based on the ice volume, the ice shape on the surface of the transmission line after the set time step is obtained. The embodiment of the present application simulates the ice growth process of the transmission line in a dynamic time-varying environment based on measured meteorological data and the finite element analysis method, rather than simulating the icing of the transmission line based on the finite element analysis method under fixed single meteorological conditions. Therefore, the embodiment of the present application can more accurately simulate the icing process and results of the transmission line in the actual environment, and effectively solve the problem of large deviation between the simulation results and the actual icing results of the transmission line icing under complex time-varying meteorological conditions; and the use of a 2.5D grid generated by stretching a two-dimensional grid for the transmission line icing simulation can effectively improve the computing efficiency compared to directly using a three-dimensional grid for the transmission line icing simulation.

[0113] Based on the above embodiments, Figure 5 This is a flow chart of a method for simulating ice coating on a transmission line provided in another embodiment of the present application. Figure 5 As shown, the transmission line icing simulation method of the embodiment of the present application may include:

[0114] S501. In response to a simulation instruction for icing on a transmission line, create a target three-dimensional model corresponding to the transmission line, where the target three-dimensional model includes the transmission line and an external flow field around the transmission line.

[0115] The detailed description of this step can be found in Figure 1 Relevant description of S101 in the illustrated embodiment.

[0116] S502: Mesh the target three-dimensional model to obtain a calculation domain containing multiple microelements.

[0117] The detailed description of this step can be found in Figure 1 The description related to S102 in the illustrated embodiment will not be repeated here.

[0118] S503. Using the wind speed and wind direction in the measured meteorological data as boundary conditions, the air flow field in the calculation domain is simulated based on the computational fluid dynamics method to obtain the first attribute information of the air flow field. The first attribute information includes the air flow velocity, air flow direction and turbulence structure around the transmission line.

[0119] The detailed description of this step can be found in Figure 1 Relevant description of S103 in the illustrated embodiment.

[0120] S504. Using the precipitation in the measured meteorological data as a boundary condition, obtain the second attribute information of the water droplets in the air flow field based on the first attribute information. The second attribute information includes the movement path and speed of the water droplets under the action of the air flow field, as well as the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets.

[0121] The detailed description of this step can be found in Figure 1 The description related to S104 in the illustrated embodiment will not be repeated here.

[0122] S505. Using the temperature in the measured meteorological data as a boundary condition, the cooling and solidification process of the water droplets is simulated based on the second attribute information to obtain the ice thickness and ice shape on the surface of the transmission line after a set time step.

[0123] The detailed description of this step can be found in Figure 1 The description related to S105 in the illustrated embodiment will not be repeated here.

[0124] S506: Determine whether the preset ice covering time has been reached.

[0125] In this step, the preset ice covering time is, for example, 24 hours. If the preset ice covering time has not been reached, step S507 is executed; if the preset ice covering time has been reached, step S508 is executed to output the ice thickness and ice shape on the transmission line surface after the preset ice covering time.

[0126] S507. Reconstruct the ice boundary corresponding to the transmission line to obtain a reconstructed three-dimensional model of the ice boundary; use the reconstructed three-dimensional model of the ice boundary as a new target three-dimensional model, and execute steps S502 to S506 to obtain the ice thickness and ice shape on the surface of the transmission line after the next set time step.

[0127] It is understood that after each set time step, the ice thickness and shape of the transmission line surface will change. Therefore, it is necessary to reconstruct the ice boundary corresponding to the transmission line to obtain a reconstructed three-dimensional model of the ice boundary. The reconstructed three-dimensional model of the ice boundary is then used as the new target three-dimensional model, and steps S502 to S505 are executed to obtain the ice thickness and shape of the transmission line surface after the next set time step. During the preset ice duration, by iteratively executing steps S502 to S505 to perform a finite element numerical solution, the ice thickness and shape of the transmission line surface after each set time step can be obtained. Thus, after reaching the preset ice duration, the ice thickness and shape of the transmission line surface after the preset ice duration are obtained.

[0128] The transmission line icing simulation method provided in the embodiment of the present application successively carries out the simulation of the air flow field, the simulation of the water droplet movement and collision collection, and the calculation of the water droplet freezing and icing growth in each set time step; the icing boundary based on the grid redivision is automatically updated, and the ice thickness and ice shape on the transmission line surface after each set time step are iteratively obtained until the preset icing duration is reached, so that the ice thickness and ice shape on the transmission line surface after the preset icing duration can be output, and the icing growth process of the transmission line in a dynamic time-varying environment is simulated based on the measured meteorological data and the finite element analysis method, so that the simulation results of the ice thickness and ice shape that are more in line with the actual situation can be obtained.

[0129] The following are device embodiments of the present application, which can be used to implement the method embodiments of the present application. For details not disclosed in the device embodiments of the present application, please refer to the method embodiments of the present application.

[0130] Figure 6 This is a schematic diagram of the structure of a transmission line ice coating simulation device provided in one embodiment of the present application. Figure 6 As shown, the transmission line icing simulation device 600 of the embodiment of the present application includes: a creation module 601, a division module 602, a first acquisition module 603, a second acquisition module 604 and a third acquisition module 605. Among them:

[0131] The creation module 601 is used to create a target three-dimensional model corresponding to the transmission line in response to a simulation instruction for transmission line icing, where the target three-dimensional model includes the transmission line and an external flow field around the transmission line.

[0132] The division module 602 is used to perform mesh division on the target three-dimensional model to obtain a calculation domain containing multiple microelements.

[0133] The first acquisition module 603 is used to simulate the air flow field in the calculation domain based on the computational fluid dynamics method using the wind speed and wind direction in the measured meteorological data as boundary conditions to obtain the first attribute information of the air flow field. The first attribute information includes the air flow velocity, air flow direction and turbulence structure around the transmission line.

[0134] The second acquisition module 604 is used to obtain second attribute information of water droplets in the air flow field based on the first attribute information, using the precipitation in the measured meteorological data as a boundary condition. The second attribute information includes the movement path and speed of the water droplets under the action of the air flow field, as well as the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets.

[0135] The third acquisition module 605 is used to simulate the cooling and solidification process of water droplets based on the second attribute information using the temperature in the measured meteorological data as the boundary condition, and obtain the ice thickness and ice shape on the surface of the transmission line after a set time step.

[0136] In some embodiments, the first acquisition module 603 can be specifically used to: simulate the air flow field in the calculation domain based on the computational fluid dynamics method, obtain the airflow velocity and direction using the Navier-Stokes equations and the continuity equation, and obtain the turbulence structure through a preset turbulence model.

[0137] In some embodiments, the second acquisition module 604 may be specifically configured to acquire the second attribute information according to the following first formula:

[0138]

[0139] Where α represents the volume fraction of the water droplet phase; Indicates air speed; Indicates the water droplet velocity; C D Represents the water drop resistance coefficient; d Reynolds number for the relative motion of water droplets; ρ a represents the density of air; ρ d represents the density of the water droplet; Fr represents the local Froude number; K represents the inertial parameter of the water droplet; V a,∞ represents the free-flow air velocity; D represents the diameter of the transmission line; L ∞ represents characteristic length; μ represents air kinematic viscosity; M represents precipitation; and v represents the falling velocity of water droplets.

[0140] In some embodiments, the third acquisition module 605 can be specifically used to: simulate the cooling and solidification process of water droplets based on the second attribute information to obtain the ice mass and ice density; based on the ice mass and ice density, obtain the ice volume and ice thickness on the surface of the transmission line after a set time step; based on the ice volume, obtain the ice shape on the surface of the transmission line after a set time step.

[0141] Optionally, when the third acquisition module 605 is used to simulate the cooling and solidification process of the water droplet based on the second attribute information to obtain the ice coating mass, it can be specifically used to obtain the ice coating mass according to the following second formula:

[0142]

[0143] Among them, h f Indicates the height of the water film; represents the velocity of the water film; ρ f Indicates the density of the water film; V ∞ represents the free-flow air velocity; a represents the volume fraction of the water droplet phase; It represents the mass of water droplets evaporated from the microelement; represents the mass of frozen water droplets in the microelement; c f represents the specific heat capacity of water; c s represents the specific heat capacity of ice; represents the equilibrium temperature of the ice-covered interface; Indicates the free-flow air temperature; represents the average temperature of the ice layer; represents the ambient temperature after taking into account the energy loss caused by friction during convection; represents the local collision velocity of water droplets; L evap represents the latent heat of evaporation of water; L fusion represents the latent heat of melting of ice; σ represents the blackbody radiation constant; ε represents the blackness; A represents the radiated area of the microelement; h c Indicates the convection heat transfer coefficient between the external environment and the transmission line surface; Q anti-icing represents the anti-icing heat flux; β represents the local collision coefficient.

[0144] Optionally, when the third acquisition module 605 is used to simulate the cooling and solidification process of the water droplet based on the second attribute information to obtain the ice density, it can be specifically used to obtain the ice density according to the following third formula:

[0145]

[0146] Among them, R M represents the calculation parameter of ice density of transmission line; d represents the diameter of supercooled water droplet; represents the collision velocity of supercooled water droplets; Indicates the ice surface temperature.

[0147] In some embodiments, the partitioning module 602 may be specifically configured to perform 2.5-dimensional mesh partitioning on the target three-dimensional model to obtain a computational domain containing a plurality of microelements.

[0148] Optionally, the transmission line icing simulation device 600 may further include a processing module ( Figure 6 (not shown) is used to: after obtaining the ice thickness and ice shape on the transmission line surface after a set time step, reconstruct the ice boundary corresponding to the transmission line to obtain a three-dimensional model of the reconstructed ice boundary; use the three-dimensional model of the reconstructed ice boundary as a new target three-dimensional model, perform meshing on the target three-dimensional model, and obtain a calculation domain containing multiple microelements to obtain the ice thickness and ice shape on the transmission line surface after the next set time step.

[0149] The device of the embodiment of the present application can be used to execute the technical solution of any of the above-mentioned method embodiments. Its implementation principles and technical effects are similar and will not be repeated here.

[0150] Figure 7 This is a schematic diagram of the structure of an electronic device provided in one embodiment of the present application. Figure 7 As shown, the electronic device 700 may include: at least one processor 701 and a memory 702 .

[0151] The memory 702 is used to store programs. Specifically, the programs may include program codes, and the program codes include computer-executable instructions.

[0152] The memory 702 may include a high-speed random access memory (RAM), and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0153] Processor 701 is configured to execute computer-executable instructions stored in memory 702 to implement the transmission line icing simulation method described in the aforementioned method embodiment. Processor 701 may be a central processing unit (CPU), an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application. Specifically, when implementing the transmission line icing simulation method described in the aforementioned method embodiment, the electronic device may be, for example, a server or other electronic device with processing capabilities.

[0154] Optionally, the electronic device 700 may further include a communication interface 703. In a specific implementation, if the communication interface 703, the memory 702, and the processor 701 are implemented independently, the communication interface 703, the memory 702, and the processor 701 may be interconnected via a bus and communicate with each other. The bus may be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus. Buses can be divided into address buses, data buses, control buses, etc., but this does not mean that there is only one bus or only one type of bus.

[0155] Optionally, in a specific implementation, if the communication interface 703, the memory 702 and the processor 701 are integrated on a chip, the communication interface 703, the memory 702 and the processor 701 can complete communication through an internal interface.

[0156] The present application also provides a computer-readable storage medium, in which computer program instructions are stored. When a processor executes the computer program instructions, the above-mentioned method for simulating icing on a transmission line is implemented.

[0157] The present application also provides a computer program product, including a computer program, which, when executed, implements the above-mentioned method for simulating icing on a transmission line.

[0158] The computer-readable storage medium may be implemented by any type of volatile or non-volatile memory device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk. The computer-readable storage medium may be any available medium that can be accessed by a general-purpose or special-purpose computer.

[0159] An exemplary readable storage medium is coupled to the processor, enabling the processor to read information from and write information to the readable storage medium. Of course, the readable storage medium may also be an integral part of the processor. The processor and the readable storage medium may be located in an application-specific integrated circuit. Of course, the processor and the readable storage medium may also exist as discrete components in the transmission line icing simulation device.

[0160] Those skilled in the art will appreciate that all or part of the steps in the above-described method embodiments can be implemented using hardware associated with program instructions. The aforementioned program can be stored in a computer-readable storage medium. When executed, the program performs the steps of the above-described method embodiments. The aforementioned storage medium includes various media capable of storing program code, such as ROM, RAM, magnetic disks, or optical disks.

[0161] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some or all of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for simulating ice coating on a transmission line, characterized in that: include: In response to a simulation instruction for icing on a transmission line, creating a target three-dimensional model corresponding to the transmission line, the target three-dimensional model including the transmission line and an external flow field around the transmission line; Meshing the target three-dimensional model to obtain a computational domain containing a plurality of microelements; The air flow field in the computational domain is simulated based on a computational fluid dynamics method using wind speed and wind direction in measured meteorological data as boundary conditions to obtain first attribute information of the air flow field, the first attribute information including air flow velocity, air flow direction, and turbulence structure around the transmission line; Using the precipitation in the measured meteorological data as a boundary condition, obtaining second attribute information of water droplets in the air flow field based on the first attribute information, the second attribute information including the movement path and speed of the water droplets under the action of the air flow field, as well as the amount of water droplets collected on the surface of the transmission line and the spatial distribution of the water droplets; Taking the temperature in the measured meteorological data as a boundary condition, the cooling and solidification process of the water droplets is simulated based on the second attribute information to obtain the ice thickness and ice shape on the surface of the transmission line after a set time step.

2. The method for simulating ice coating on a transmission line according to claim 1, characterized in that: The simulating the air flow field in the calculation domain based on the computational fluid dynamics method to obtain first attribute information of the air flow field includes: The air flow field in the calculation domain is simulated based on a computational fluid dynamics method, the air flow velocity and the air flow direction are obtained using the Navier-Stokes equation and the continuity equation, and the turbulence structure is obtained through a preset turbulence model.

3. The method for simulating ice coating on a transmission line according to claim 2, characterized in that: The acquiring second attribute information of the water droplets in the air flow field based on the first attribute information includes: The second attribute information is obtained according to the following first formula: Where α represents the volume fraction of the water droplet phase; Indicates air speed; Indicates the water droplet velocity; C D Represents the water drop resistance coefficient; d Reynolds number for the relative motion of water droplets; ρ a represents the density of air; ρ d represents the density of the water droplet; Fr represents the local Froude number; K represents the inertial parameter of the water droplet; V a,∞ represents the free air velocity; D represents the diameter of the transmission line; L ∞ represents characteristic length; μ represents air kinematic viscosity; M represents precipitation; and v represents the falling velocity of water droplets.

4. The method for simulating ice coating on a transmission line according to claim 1, wherein: The simulating the cooling and solidification process of the water droplet based on the second attribute information to obtain the ice thickness and ice shape on the surface of the transmission line after a set time step includes: Simulating the cooling and solidification process of the water droplets based on the second attribute information to obtain the ice mass and ice density; Based on the ice mass and the ice density, obtaining the ice volume and ice thickness on the surface of the transmission line after the set time step; Based on the ice-covered volume, the ice-covered shape of the surface of the transmission line after the set time step is obtained.

5. The method for simulating ice coating on a transmission line according to claim 4, characterized in that: The simulating the cooling and solidification process of the water droplets based on the second attribute information to obtain the ice coating mass includes: The ice coverage mass is obtained according to the following second formula: Among them, h f Indicates the height of the water film; represents the velocity of the water film; ρ f Indicates the density of the water film; V ∞ represents the free-flow air velocity; a represents the volume fraction of the water droplet phase; It represents the mass of water droplets evaporated from the microelement; represents the mass of frozen water droplets in the microelement; c f represents the specific heat capacity of water; c s represents the specific heat capacity of ice; represents the equilibrium temperature of the ice-covered interface; Indicates the free-flow air temperature; represents the average temperature of the ice layer; represents the ambient temperature after taking into account the energy loss caused by friction during convection; represents the local collision velocity of water droplets; L evap represents the latent heat of evaporation of water; L fusion represents the latent heat of melting of ice; σ represents the blackbody radiation constant; ε represents the blackness; A represents the radiated area of the microelement; h c Indicates the convection heat transfer coefficient between the external environment and the transmission line surface; Q anti-icing represents the anti-icing heat flux; β represents the local collision coefficient.

6. The method for simulating ice coating on a transmission line according to claim 4, characterized in that: The cooling and solidification process of the water droplet is simulated based on the second attribute information to obtain the ice density, including: The ice density is obtained according to the following third formula: Among them, R M represents the calculation parameter of ice density of transmission line; d represents the diameter of supercooled water droplet; represents the collision velocity of supercooled water droplets; Indicates the ice-covered surface temperature.

7. The method for simulating ice coating on a transmission line according to any one of claims 1 to 6, characterized in that: The target three-dimensional model is meshed to obtain a computational domain containing a plurality of microelements, including: The target three-dimensional model is divided into 2.5-dimensional grids to obtain a calculation domain containing multiple microelements.

8. The method for simulating ice coating on a transmission line according to any one of claims 1 to 6, characterized in that: After obtaining the ice thickness and ice shape of the transmission line surface after the set time step, the method further includes: Reconstructing the ice boundary corresponding to the transmission line to obtain a reconstructed three-dimensional model of the ice boundary; The three-dimensional model after reconstructing the ice boundary is used as a new target three-dimensional model, and the step of meshing the target three-dimensional model is performed to obtain a calculation domain containing multiple microelements, so as to obtain the ice thickness and ice shape on the surface of the transmission line after the next set time step.

9. A transmission line icing simulation device, characterized in that: include: a creation module, configured to create, in response to a simulation instruction for icing of a transmission line, a target three-dimensional model corresponding to the transmission line, the target three-dimensional model including the transmission line and an external flow field around the transmission line; A partitioning module, configured to perform mesh partitioning on the target three-dimensional model to obtain a computational domain containing a plurality of microelements; a first acquisition module, configured to simulate the air flow field in the calculation domain based on a computational fluid dynamics method using wind speed and wind direction in measured meteorological data as boundary conditions, and obtain first attribute information of the air flow field, wherein the first attribute information includes air flow velocity, air flow direction, and turbulence structure around the transmission line; a second acquisition module, configured to acquire, based on the first attribute information and taking the precipitation in the measured meteorological data as a boundary condition, second attribute information of water droplets in the air flow field, the second attribute information including a movement path and speed of the water droplets under the action of the air flow field, and an amount of water droplets collected on a surface of the transmission line and a spatial distribution of the water droplets; The third acquisition module is used to simulate the cooling and solidification process of water droplets based on the second attribute information using the temperature in the measured meteorological data as a boundary condition, and obtain the ice thickness and ice shape on the surface of the transmission line after a set time step.

10. An electronic device, characterized in that: include: a processor, and a memory communicatively connected to the processor; The memory stores computer-executable instructions; The processor executes the computer-executable instructions stored in the memory to implement the transmission line icing simulation method according to any one of claims 1 to 8.

11. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer program instructions, and when the computer program instructions are executed, the transmission line icing simulation method according to any one of claims 1 to 8 is implemented.

12. A computer program product comprising a computer program, characterized in that When the computer program is executed, the transmission line icing simulation method according to any one of claims 1 to 8 is implemented.

Images