CN112837417B - Three-dimensional space lightning lead simulation method based on improved random model - Google Patents
Three-dimensional space lightning lead simulation method based on improved random model Download PDFInfo
- Publication number
- CN112837417B CN112837417B CN202011430085.9A CN202011430085A CN112837417B CN 112837417 B CN112837417 B CN 112837417B CN 202011430085 A CN202011430085 A CN 202011430085A CN 112837417 B CN112837417 B CN 112837417B
- Authority
- CN
- China
- Prior art keywords
- lightning
- model
- simulation
- client
- comsol
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
The invention discloses a three-dimensional space lightning lead simulation method based on an improved random model, which comprises the following steps of: s1, importing a geometric model of an initial environment, and adjusting the shape and the size of an object and the starting point of lightning lead according to simulation requirements; s2, establishing connection between the Comsol client and the server, and automatically uploading static data of the model to the server by the Comsol client; s3, obtaining model data through codes; s4, the Comsol client obtains an electric field module value of the pilot surface by using a Comsol interface, and the Matlab client interpolates by using a unit shape function to determine a development point of the next iteration pilot; s5, judging whether the electric field module at the boundary of the building surface meets the condition of uplink pilot excitation, and if so, ending the simulation; s6, controlling the Comsol to generate a lightning downlink pilot at the point through an interface to obtain an electric field module value of the whole area, and then performing step S4.
Description
Technical Field
The invention relates to the field of lightning lead simulation in a three-dimensional space, in particular to a three-dimensional space lightning lead simulation method based on an improved random model.
Background
Lightning is a high intensity, long distance transient discharge phenomenon in the atmosphere. Scientific-based lightning protection measures require a thorough understanding of the physical processes and mechanisms of lightning-based lightning. One effective way to study the physical mechanisms of lightning is through modeling of the electrification and discharge processes based on observations. Currently, there are two main types of lightning simulation models, namely physical models and stochastic models.
The shape of the lightning leader tip, the channel, and the induced charge in the corona region are taken into account in the physical model. However, these properties are difficult to obtain. In contrast, the stochastic model requires fewer factors to consider and is easier to implement. To some extent, the lightning leader it develops is more consistent with the branching and tortuosity characteristics that often occur with real lightning.
In addition, the method can be divided into a downlink negative-step pilot model, an uplink head-on positive pilot model and a complete pilot model simultaneously comprising an uplink pilot and a downlink pilot according to the characteristics of the pilot model. With the continuous deep understanding of the pilot physical process and the improvement of the observation technology, a plurality of different pilot models appear to solve the lightning protection problem of ground buildings and various electronic and electric equipment (such as weather radar, high-voltage transmission lines, wind driven generators and the like). The above models may be collectively referred to as a pilot propagation model. In various propagation models, the simulation results are different due to the fact that different assumed conditions or physical quantities are adopted.
Zhang Xin the three-dimensional numerical simulation of lightning on the upper side of a building [ D ]. Nanjing information engineering university, 2017. The three-dimensional numerical simulation of the document is mainly programmed by a Fortran program to realize the solution of electrostatic fields, and MATLAB software is adopted for the output of the graph. Although the generation of side-impact lightning on a building in a three-dimensional space is simulated, in the three-dimensional space, the simulated calculation amount is increased along with the increase of the number of lattice points, and only the side-impact lightning simulation analysis of a single building can be discussed under the limitation of the calculation amount, but the influence of building groups on the side-impact lightning is not deeply discussed.
Ren Xiao simulation study of the interaction of lightning leader with ground object [ D ]. Beijing: the two-dimensional pilot propagation model in the Chinese meteorological science institute 2010 adopts a square grid of 0.1m multiplied by 0.1m, is limited by grid subdivision and computer memory, and only calculates the space within 200m near the ground. However, the two-dimensional pilot propagation model does not reflect the interaction of the pilot and the three-dimensional structure building. Meanwhile, the adoption of a finite difference method is inconvenient for meshing the irregular complex structure target.
Mo Haojiang, wei Guanghui, chen Jiang, etc. three-dimensional numerical simulation of lightning pilot discharge and application of [ J ]. High voltage technology, 2013. Pilot propagation model in this document focuses on random propagation, pursuing similarity of geometry, weakens description of lightning physical process. Due to the purely Matlab-based implementation, no description of the lightning physics process can be added and modeling in complex scenarios is also difficult.
Xu, yazhong & Chen, mini (2013). A3-D Self-Organized Leader Propagation Model and Its Engineering Approximation for Lightning Protection analysis.power Delivery, IEEE Transactions on.28.2342-2355.10.1109/tpwrd.2013.2263846. The three-dimensional Self-organizing downstream leader model in this document weakens the effect of the upstream leader, only considers the interaction of the downstream leader with the ground, and fails to consider the effect of the space charge layer formed by corona ions generated at the tip of the lightning rod or the like. In the numerical simulation, the contradiction between grid sectioning and computer memory distribution of different scale study objects (such as a pilot channel, a lightning rod, a ground building, a simulation space and the like) is also required to be considered. With the development of economy, super-high buildings are more and more, and the current lightning protection standards and standards at home and abroad mainly aim at buildings below 60m, and the research range is narrower.
In summary, the current technical situation is mainly limited by the large calculation amount problem generated under the fine resolution, the mesh division problem under the fine granularity, and the difficulty in describing the shape of the building and its potential distribution through the programming language under the complex scene.
For this reason, a rapid and flexible lightning lead simulation method in a three-dimensional space is required.
Disclosure of Invention
The invention aims to solve the problems of ineffective calculation and simulation in complex scenes caused by long-time and uniform grid division in the current lightning lead simulation process in a three-dimensional space, and by the method, the overall simulation efficiency can be greatly improved, the geometric and image functions of finite element simulation software Comsol can be fully utilized, the development profile of lightning in the simulation process can be fully displayed, and an algorithm capable of rapidly carrying out lightning lead simulation for various different building scenes in the three-dimensional space is provided.
The invention is realized at least by one of the following technical schemes.
A three-dimensional space lightning lead simulation method based on an improved random model comprises the following steps:
s1, importing a three-dimensional model file of an initial environment by a Comsol client, and adjusting the shape and the size of a three-dimensional object in the file, an initial potential and a starting point of lightning lead according to simulation requirements;
s2, establishing connection between the Comsol client and the server, and automatically uploading the three-dimensional object model data to the server by the Comsol client;
s3, establishing connection between the Matlab client and the server, and acquiring three-dimensional object model data by the Matlab client through a Comsol interface;
s4, the Matlab client obtains an electric field module value of the pilot surface in the step S3 by using a Comsol interface, the Matlab client interpolates by using a unit shape function, and a development point of the next iteration pilot is determined according to an improved random lightning pilot development model;
s5, obtaining an electric field module value of the surface boundary of the three-dimensional object, interpolating by utilizing a unit shape function, judging whether the electric field module of the surface boundary of the three-dimensional object meets the condition of uplink pilot excitation, if yes, finishing simulation, otherwise, performing step S6;
s6, controlling the Comsol client to send lightning downlink pilot at the development point through a Comsol interface, adjusting geometric modeling, resetting potential boundaries and grounding boundaries, solving an electrostatic field based on the adjusted geometric model to obtain an electric field model value of a whole area, and then performing step S4.
Preferably, the server is mainly used for transmitting code control instructions and various data to the client, acquiring intermediate states of simulation processes in real time, and for realizing mutual exclusive access of multiple clients to simulation data in the server, closing other clients by utilizing a function of temporarily storing dynamic data of the server, only keeping independent simulation of a Maltab control component, and starting the Comsol client to asynchronously update model data after the simulation is finished.
Preferably, the improved random lightning lead development model is a lightning lead development model based on an improved random model, wherein eta is introduced into the improved random model, each iteration of lightning lead simulation is developed at only one point, and the power distribution model depends on eta and an electric field model E on a point meeting a condition, namely the power distribution probability of each point is as follows:
wherein E is x,y,z η represents the power of η of the electric field modes of points in three-dimensional space at coordinates x, y, z Σ x,y,z η represents the electric field mode of a point in three-dimensional space meeting the pilot development condition in a certain simulation iterationAnd the sum of power eta.
Preferably, the electric field modulus value in the step S4 to the step S6 is obtained by an electrostatic field solving algorithm built in the Comsol software.
Preferably, the interpolation by using the unit shape function is to use local accumulation to form a full area for interpolation, specifically, the coordinates of each iteration development point are recorded, and the interpolation area selection principle is as follows: taking the recorded point as the center, taking the radius of the sphere as an interpolation area, respectively interpolating the spheres at the center point, and selecting the point which meets the condition of generating the downlink pilot.
Preferably, for the simulation area, the improved random lightning lead development model regards the bottom of the thunderstorm cloud and the ground as parallel plate capacitors, the upper boundary as the bottom of the live thundercloud, the lower boundary as the ground surface, and the left and right boundaries far enough from the simulated lightning channel to meet the Newman boundary condition, and is regarded as an infinite open space.
Preferably, only the tip of the upstream pilot can develop to one of the surrounding undeveloped points, the ratio of the downstream pilot to the upstream pilot is set, and the upstream pilot iterates one round each time the downstream pilot iterates several rounds.
Preferably, the condition for judging whether the electric field mode of the boundary of the three-dimensional object surface meets the uplink pilot excitation is as follows: the average electric field intensity between the head of the downward leading of the thunder and lightning and the head of the upward leading or the ground object is larger than the final jump electric field intensity E f Representing that the excitation condition is satisfied, if the excitation condition is satisfied, the simulation is ended, otherwise, the step S6 is performed.
Preferably, the electric field intensity E of a point adjacent to the lightning leading channel in the simulation space x,y,z Greater than critical field strength E s And when the lightning leading channel is in the next development point, the point is the next development point of the lightning leading channel.
Preferably, the Comsol client, matalb client and server communicate in C/S mode
The simulation efficiency is improved mainly by the following modes:
(1) A Comsol built-in electrostatic field solving algorithm is used for realizing optimal solving time and occupying minimum memory space;
(2) The self-adaptive grid dividing method adopted by the Comsol is utilized to intelligently select the thickness of grids in a three-dimensional space, fine grids are used in important areas such as a guide area, a building area and the like, and coarse grids are used in unimportant areas, so that the purposes of optimizing the memory and improving the simulation speed are achieved;
(3) Acquiring an electric field model after solving a Comsol electrostatic field through a Matalb interface opened by the Comsol, and carrying out interpolation through a unit shape function on the Matlab to obtain an electric field model value of uniform grid division in the lightning lead development area;
(4) According to the interpolation result, according to a lightning lead development model based on an improved random model, the development position of the next iteration lead is obtained and transmitted back to the Comsol, the development of the lightning lead is completed by utilizing the geometric drawing function of the Comsol, and the lightning lead is completed by abstracting a cylinder and a sphere;
(5) Realizing automatic simulation by utilizing Matlab codes until the building generates uplink pilot, and finishing the simulation;
in the invention, the Comsol client mainly performs the following two tasks: (1) And (3) completing the geometric drawing of the building and the initial lead according to the simulation task, and completing the development of the downlink lead and the uplink lead in the simulation process. (2) And calculating an electric field model of the electrostatic field in the new scene according to the geographical position of the building and the development condition of the lightning lead.
In the invention, the Matlab client mainly executes the following three tasks:
(1) And finishing the conversion from the non-uniform electric field mode grid to the uniform electric field mode grid through unit shape function interpolation.
(2) And completing the interface call of all the geometric drawing.
(3) And realizing automatic simulation.
Compared with the prior art, the invention has the beneficial effects that:
1. according to the invention, the solution of the electrostatic field in the three-dimensional space is a solution of an electrostatic field solution algorithm under the three-dimensional space built in the Comsol, so that the calculation of an electric field mode in the electrostatic field can be efficiently and stably completed, the abnormal influence caused by improper boundary condition setting, numerical value overflow, external interference and other reasons is reduced, and an effective simulation result is ensured.
2. The method uses the unit shape function to interpolate to obtain the uniform grid, and compensates the non-uniform grid caused by the self-adaptive grid. The conversion from the non-uniform grid to the uniform grid can complete the development of lightning lead by using an improved random model under a regular electric field model, compared with the traditional lightning development model, the granularity of the electric field model in a local area is increased, the characteristics of more discrimination are increased, and more real lightning development behaviors can be simulated relatively accurately. The regular electric field mode is a necessary condition for flexibly designing a lightning lead development model, otherwise, simulation is performed under an error background.
3. By introducing the Comsol and by means of the geometric drawing capability, the invention can simulate the lightning development condition in a complex scene, and improves the universality of the invention.
Drawings
FIG. 1 is a block diagram of a three-dimensional space lightning lead simulation method based on an improved random model;
FIG. 2 is a flow chart of a process for automated simulation in accordance with the present invention;
FIG. 3 is a flow chart of the present invention based on an improved random lead development model;
FIG. 4 is a diagram of the automated simulation initiation conditions of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
As shown in fig. 2, the three-dimensional space lightning lead simulation method based on the improved random model of the embodiment comprises the following steps:
s1, importing mph files of an initial environment by a Comsol client, or describing CAD files of a simulation area, and adjusting the shape and the size of a three-dimensional object in the files, initial potential and a starting point of lightning lead according to simulation requirements;
s2, establishing connection between the Comsol client and the server, and automatically uploading static data of the three-dimensional object model to the server by the Comsol client;
s3, establishing a Matlab client and a server, wherein the Matlab client acquires three-dimensional object model data through a Comsol interface;
the interface code is: model = Model. Model ('Model 2'), model data is obtained;
s4, the Matlab client obtains an electric field module value of the pilot surface in the step S3 by using a Comsol interface, the Matlab client interpolates by using a unit shape function, and the development point of the next iteration pilot is determined according to the improved random lightning pilot development model.
FIG. 3 is a flow chart of the invention based on an improved random lead development model based primarily on a three-dimensional dielectric breakdown model, taking into account both deterministic and stochastic effects in the lightning lead channel development process. For rectangular simulation areas, since the bottom of the thunderstorm cloud and the ground can be regarded as parallel plate capacitors, the upper boundary can be regarded as a charged Lei Yunde part (set constant potential Φ=50000000v), the lower boundary can be regarded as the surface of the earth (set constant potential Φ=0v), and the left and right boundaries can be regarded as infinite open space because they are far enough from the simulated lightning channel to satisfy the newman boundary condition. In the simulation, it is assumed that the pilot channel has been fully ionized and that the potential of the lightning down-lead does not change with the development of the channel. The lightning down leader assumes the same potential as the lightning bottom, the up leader of the super-high building up development assumes the same potential as the ground potential, and no other charge exists in the free space between the cloud and the ground. The descending leader may create a number of branches while the ascending leader generally has no significant branching structure near the ground area. Thus, in the modified random lead development model, each deployment point of the downstream lead can be deployed to a nearby undeployed point. In order to avoid improving branches of the upstream pilot in the random pilot development model, only the tip of the upstream pilot can develop to one of the surrounding undeveloped points. Setting the ratio of the downlink pilot to the uplink pilot as 4:1, that is, each time the downstream pilot iterates four rounds, the upstream pilot iterates one round.
The potential distribution and electric field intensity distribution of the simulation space can be built-in static electricity by using the Comsol softwareAnd obtaining a field solving algorithm. When the electric field intensity E of a point adjacent to the existing lightning leading channel in the simulation space x,y,z Greater than critical field strength E s When each lightning lead simulation iteration is performed, only one point is developed, depending on eta and an electric field mode E on a point meeting a condition, but only one point is developed in each step, and the probability of each point is as follows:
wherein, E x,y,z η represents the power of η of the electric field modes of points in three-dimensional space at coordinates x, y, z Σ x,y,z η represents the sum of the power of η of the electric field modes of points in the three-dimensional space that meet the pilot development conditions in a certain simulation iteration. η is a development probability weight index, the larger η is the size of randomness of the downward pilot development of the table lightning, the stronger the correlation between the downward pilot development of the lightning and the electric field intensity distribution is, the stronger the extensibility of a lightning pilot channel is, and the less the channel bifurcation phenomenon is; and vice versa. When the average electric field intensity between the head of the lightning descending leader and the head of the ascending leader or the ground object is larger than the final jump electric field intensity E f And when the lightning leading connection process and the back striking are considered to be about to happen, the development simulation process of the lightning descending leading is ended.
The implementation of the lightning lead development model based on the improved random model in the three-dimensional space is based on a uniform electric field model grid. The uniform electric field mode grid is interpolated from the cell shapes. Meanwhile, the interpolation area adopts a mode of forming a full area by local accumulation, so that the interpolation of an invalid area is avoided, and memory resources are occupied.
S5, obtaining an electric field module value of the building surface boundary, interpolating by utilizing a unit shape function, and judging whether the electric field module of the building surface boundary meets the uplink pilot (meeting the condition that the average electric field intensity between the head of the lightning downlink pilot and the uplink pilot head or the ground object is larger than the final jump electric field intensity E) f Represents meeting the excitation condition, E f Typically 500,000V/m), if the excitation condition is met, the junction is simulatedAnd if not, performing step S6.
Considering the memory limitation and the granularity of interpolation, the invention abandons a cut full-area interpolation method, but records the coordinates of each iteration development point, and the interpolation area selection principle is as follows: a sphere with a radius of 10m is used as an interpolation area with the recorded point as the center. The spheres at the points are respectively interpolated, and the electric field intensity E which meets the requirement of generating downlink pilot (meeting the condition that a certain point in the simulation space is adjacent to the lightning pilot channel is selected x,y,z Greater than critical field strength E s (E s Generally taking 500000V/m), consider that the point may be the next development point of the lightning pilot channel, and there may be multiple points in the three-dimensional space that meet the pilot development condition, but only one point of each iteration of the lightning pilot simulation is developed, where the probability of each point is as follows:
wherein η represents a development probability index, generally taking 3. And (3) calculating a probability value, and taking the point development with the maximum probability, namely the point with the maximum probability meeting the condition. Compared with the full-area interpolation, the interpolation method using the partial accumulation to form the full-area interpolation can greatly reduce the required memory while ensuring the simulation effect. Because the interpolation is carried out by using a local accumulation mode, not only can the aim of covering all reachable areas of lightning development be achieved, but also unreachable areas can be adaptively screened out, and the feasible solution range of electrostatic field calculation is reduced.
S6, controlling the Comsol software to emit lightning downlink pilot at the point through a Comsol software application program interface, adjusting geometric modeling, resetting potential boundaries and grounding boundaries, solving an electrostatic field based on the adjusted geometric model to obtain an electric field module value of a whole area, and then performing step S4. The interface code is as follows:
% is the current lightning leader designation
gemo_name＝strcat(num2str(i))；
% New lightning development lead
model.geom('geom1').create(gemo_name,'Cylinder')；
model.geom('geom1').feature(gemo_name).set('r',{'1'})；
% calculation of the angle and the length of downward development of lightning leader
base＝newMap([num2str(m(1)),num2str(m(2)),num2str(m(3)+10)])；
base = base'; % transposition
xm＝m(1)-base(1)；
ym＝m(2)-base(2)；
zm＝m(3)-base(3)；
xyz＝[xm ym zm]；
cylen＝sqrt(xm*xm+ym*ym+zm*zm)；
model.geom('geom1').feature(gemo_name).set('ax3',xyz)；
model.geom('geom1').feature(gemo_name).set('h',cylen)；
model.geom('geom1').feature(gemo_name).set('pos',base)；
model.geom('geom1').run；
The calculation of the electric field mode in the electrostatic field is efficiently and stably completed through the built-in electrostatic field solving algorithm of the Comsol software, the abnormal influence caused by improper setting of boundary conditions, numerical overflow, external interference and other reasons is reduced, and an effective simulation result is ensured. Fig. 4 is an automated simulation initial condition diagram of the present embodiment.
The unit shape function of the invention carries out interpolation to obtain uniform grids, and compensates non-uniform grids caused by self-adaptive grids. The conversion from the non-uniform grid to the uniform grid can complete the development of lightning lead by using an improved random model under a regular electric field model, compared with the traditional lightning development model, the granularity of the electric field model in a local area is increased, the characteristics of more discrimination are increased, and more real lightning development behaviors can be simulated relatively accurately.
By introducing the Comsol software and by means of the geometric drawing capability, the invention can simulate the lightning development condition in a complex scene, and improves the universality of the invention. According to the invention, the simulation of the lightning lead is completed by linking the Comsol software and the Matlab software and respectively playing the characteristics of the Comsol software and the Matlab software, so that the simulation speed can be effectively reduced. By utilizing the geometric capability of the Comsol software, lightning simulation research can be carried out in various building scenes, and a foundation is laid for lightning protection research.
As shown in FIG. 1, the system of the three-dimensional space lightning lead simulation method based on the improved random model comprises a Comsol client, a Matalb client and a COMSOL Multiphysics server server, wherein the three components are communicated in a C/S mode, and dynamic data resources and static data resources are stored in the COMSOL Multiphysics server server for each client to use and update. And the Matlab is used as an automatic simulation client component and a lightning lead model development control component by using the Comsol as a client component for solving electrostatic fields and geometrically drawing in a three-dimensional space.
The server side is mainly used for transmitting code control instructions and various data to the client side, acquiring intermediate states of simulation processes in real time, closing other client sides by utilizing the function of temporarily storing dynamic data of the server in order to realize mutual exclusive access of the multiple client sides to the simulation data in the server side, only keeping independent simulation of a Maltab control component, and restarting the Comsol client side to asynchronously update model data after the simulation is finished, so that final consistency of the model data is realized.
The above examples are preferred embodiments of the present invention, but the embodiments of the present invention are not limited to the above examples, and any other changes, modifications, substitutions, combinations, and simplifications that do not depart from the spirit and principle of the present invention should be made in the equivalent manner, and the embodiments are included in the protection scope of the present invention.
Claims (8)
1. The three-dimensional space lightning lead simulation method based on the improved random model is characterized by comprising the following steps of:
s1, importing a three-dimensional object model file of an initial environment by a Comsol client, and adjusting the shape and the size of a three-dimensional object in the file, an initial potential and a starting point of a lightning guide according to simulation requirements;
s2, establishing connection between the Comsol client and the server, and automatically uploading the three-dimensional object model file to the server by the Comsol client;
s3, establishing connection between the Matlab client and the server, and acquiring three-dimensional object model data by the Matlab client through a Comsol interface;
s4, the Matlab client obtains an electric field module value of the pilot surface in the step S3 by using a Comsol interface, the Matlab client interpolates by using a unit shape function, and a development point of the next iteration pilot is determined according to an improved random lightning pilot development model;
the improved random lightning lead development model is based on a three-dimensional dielectric breakdown model, for a simulation area, the bottom of thunderstorm cloud and the ground are regarded as parallel plate capacitors, the upper boundary is regarded as the bottom of live thundercloud, the lower boundary is regarded as the ground surface, the left and right boundaries are far enough from a simulated lightning channel, the Newman boundary condition is met, and the space is regarded as infinite open space;
s5, obtaining an electric field module value of the surface boundary of the three-dimensional object, interpolating by utilizing a unit shape function, judging whether the electric field module value of the surface boundary of the three-dimensional object meets the condition of uplink pilot excitation, if yes, finishing simulation, otherwise, performing step S6;
s6, controlling the Comsol client to emit lightning downlink pilot at the development point through a Comsol interface, adjusting a geometric model, resetting a potential boundary and a grounding boundary, solving an electrostatic field based on the adjusted geometric model to obtain an electric field model value of a whole area, and then performing step S4;
the interpolation by using the unit shape function is to use local accumulation to form a full area for interpolation, specifically, the coordinates of each iteration development point are recorded, and the interpolation area selection principle is as follows: taking the recorded point as the center, taking the radius of the sphere as an interpolation area, respectively interpolating the spheres at the center point, and selecting the point which meets the condition of generating the downlink pilot.
2. The three-dimensional space lightning lead simulation method based on the improved random model according to claim 1, wherein the server side is used for transmitting code control instructions and various data to the client side, acquiring intermediate states of simulation processes in real time, and for realizing mutual exclusive access of the simulation data in the server side by multiple client sides, utilizing the function of temporarily storing dynamic data of the server, closing other client sides, only keeping independent simulation of a Maltab control component, and starting the Comsol client side to asynchronously update model data after the simulation is finished.
3. The three-dimensional space lightning lead simulation method based on the improved random model according to claim 2, wherein the improved random lightning lead development model is a lightning lead development model based on the improved random model, η is introduced into the improved random model, only one point of each iteration of the lightning lead simulation is developed, and the magnitude of the electric field modulus E on the points depending on η and the meeting condition is as follows:
wherein E is x,y,z η Eta power of electric field modulus value representing point with x, y and z coordinates in three-dimensional space, sigma E x,y,z η Representing the sum of power eta of electric field modulus values of points which meet the leading development conditions in a three-dimensional space in a certain simulation iteration; η is a probability of development weight index, representing the magnitude of the randomness of the lightning downstream leader development.
4. The three-dimensional space lightning lead simulation method based on the improved random model according to claim 3, wherein the electric field modulus value in the steps S4 to S6 is obtained through a built-in electrostatic field solving algorithm of Comsol software.
5. The method for three-dimensional space lightning lead simulation based on an improved random model according to claim 4, wherein only the tip of the uplink lead can be developed to surrounding undeveloped points, the proportion of the downlink lead to the uplink lead is set, and each time the downlink lead iterates for a plurality of rounds, the uplink lead iterates for one round.
6. The improved random model-based three-dimensional space lightning lead simulation method of claim 5, wherein the condition for judging whether the electric field modulus of the boundary of the three-dimensional object surface meets the uplink lead excitation is: the average electric field intensity between the head of the downward leading of the thunder and lightning and the head of the upward leading or the ground object is larger than the final jump electric field intensity E f Representing that the excitation condition is satisfied, if the excitation condition is satisfied, the simulation is ended, otherwise, the step S6 is performed.
7. The three-dimensional space lightning lead simulation method based on the improved random model according to claim 6, wherein the electric field intensity E of a point adjacent to the lightning lead channel in the simulation space x,y,z Greater than critical field strength E s And when the lightning leading channel is in the next development point, the point is the next development point of the lightning leading channel.
8. The three-dimensional space lightning lead simulation method based on the improved random model of claim 7, wherein the coomsol client, the Matalb client and the server communicate in a C/S mode.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011430085.9A CN112837417B (en) | 2021-03-02 | 2021-03-02 | Three-dimensional space lightning lead simulation method based on improved random model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011430085.9A CN112837417B (en) | 2021-03-02 | 2021-03-02 | Three-dimensional space lightning lead simulation method based on improved random model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112837417A CN112837417A (en) | 2021-05-25 |
| CN112837417B true CN112837417B (en) | 2023-06-16 |
Family
ID=75923482
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011430085.9A Active CN112837417B (en) | 2021-03-02 | 2021-03-02 | Three-dimensional space lightning lead simulation method based on improved random model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112837417B (en) |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| ID24540A (en) * | 1998-07-27 | 2000-07-27 | Erico Int Corp | AIR LIGHTNING TERMINALS AND DESIGN METHODS AND ITS USE |
| CN108052702B (en) * | 2017-11-23 | 2021-07-20 | 上海无线电设备研究所 | Lightning partitioning method based on lightning shortest path space geometric topology principle |
-
2021
- 2021-03-02 CN CN202011430085.9A patent/CN112837417B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN112837417A (en) | 2021-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bhumkittipich et al. | Optimal placement and sizing of distributed generation for power loss reduction using particle swarm optimization | |
| Yammani et al. | A Multi-objective Shuffled Bat algorithm for optimal placement and sizing of multi distributed generations with different load models | |
| Lopes et al. | Impact of the combined integration of wind generation and small hydropower plants on the system reliability | |
| Martinez et al. | A parallel Monte Carlo method for optimum allocation of distributed generation | |
| D'alessandro et al. | A “Collection Volume Method” for the placement of air terminals for the protection of structures against lightning | |
| Ugranlı et al. | Optimal wind turbine sizing to minimize energy loss | |
| Alhazmi et al. | Power grid optimal topology control considering correlations of system uncertainties | |
| Padhi et al. | Solving dynamic economic emission dispatch problem with uncertainty of wind and load using whale optimization algorithm | |
| Alik et al. | Minimization of grounding system cost using PSO, GAO, and HPSGAO techniques | |
| Nasiraghdam et al. | Load model effect assessment on optimal distributed generation (DG) sizing and allocation using improved harmony search algorithm | |
| CN112837417B (en) | Three-dimensional space lightning lead simulation method based on improved random model | |
| Cai et al. | Cause-effect modeling and spatial-temporal simulation of power distribution fault events | |
| CN113806947A (en) | Offshore wind farm layout processing method, device and equipment | |
| Sabiha et al. | Assessment of grounding grid for enhancing wind turbine service sustainability | |
| Lückehe et al. | On evolutionary approaches to wind turbine placement with geo-constraints | |
| Fakhri et al. | Overall tidal farm optimal design–Application to the alderney race and the fromveur strait (France) | |
| CN110428164A (en) | A kind of operation of power networks state monitoring method, device, equipment and readable storage medium storing program for executing | |
| CN116362194A (en) | Wiring resource pre-allocation method, device, computing equipment and storage medium | |
| CN111146815A (en) | Intelligent power distribution network distributed power generation planning configuration method | |
| Ali et al. | Reconfiguration with simultaneous DG installation to improve the Voltage Profile in Distribution Network using Harmony Search Algorithm | |
| Ioannidis et al. | Accelerating stochastic lightning attachment simulations for the estimation of lightning incidence to overhead lines | |
| Phuangpornpitak et al. | Research Article Principle Optimal Placement and Sizing of Single Distributed Generation for Power Loss Reduction using Particle Swarm Optimization | |
| Akpınar et al. | Statistical analysis of wind speed distribution of Turkey as regional | |
| Thakkar et al. | Multi-criteria Valuation for Sustainable Autonomous Microgrid Planning: a Comparative Analysis of Technology Mix with Different MADM Techniques | |
| CN108549985B (en) | Improved Monte Carlo method for solving interval direct current power flow model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |