CN106407541A - Numerical simulation method for atmospheric electric field distortion effect of buildings under different parameters - Google Patents
Numerical simulation method for atmospheric electric field distortion effect of buildings under different parameters Download PDFInfo
- Publication number
- CN106407541A CN106407541A CN201610810981.5A CN201610810981A CN106407541A CN 106407541 A CN106407541 A CN 106407541A CN 201610810981 A CN201610810981 A CN 201610810981A CN 106407541 A CN106407541 A CN 106407541A
- Authority
- CN
- China
- Prior art keywords
- electric field
- building
- atmospheric electric
- equation
- value
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Atmospheric Sciences (AREA)
- Environmental Sciences (AREA)
- Ecology (AREA)
- Biodiversity & Conservation Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a numerical simulation method for an atmospheric electric field distortion effect of buildings under different parameters. The numerical simulation method comprises the following steps of (1) solving an atmospheric electric field distortion coefficient of the buildings through a finite difference method, wherein the finite difference method is as follows specifically: performing discretization on a researched spatial simulation region, using a finite difference equation to replace a differential equation, performing uniform division after grid discretization for data processing, and solving electric potential values of respective grid points in a discretized space range; and (2) solving the finite difference equation, performing sub-division on grids in the spatial simulation region, solving electromagnetic problems in the spatial simulation region, selecting an iteration factor, and then solving the atmospheric electric field distortion coefficient of the buildings. By performing numerical simulation on the atmospheric electric field distortion effect of the buildings, the actual measurement value of an atmospheric electric field instrument can be revised, so that the accuracy of thunder and lightning monitoring and early warning can be improved.
Description
Technical field
The present invention relates to a kind of method for numerical simulation to atmospheric electric field distortion effect for building different parameters, belong to thunder and lightning
Monitoring and warning technical field.
Background technology
Thunderstorm refers to the drastic discharge phenomenon lightninged accompanied by peals of thunder occurring in cumulonimbus, existing as a kind of local strong convective weather
As often linking together with diastrous weathers such as heavy rain, strong wind, hail, tornados, in recent years, with regard to thunder and lightning to people's life
Property safety constitutes serious threat, and therefore thunderstorm monitoring early warning work is particularly important.Atmospheric electric field detector is that thunderstorm monitoring is pre-
The important component part of alarm system, with the presence of building around the position installed due to atmospheric electric field detector, building is to atmospheric electricity
Field has certain distortion effect it is therefore desirable to study to the distortion effect of building atmospheric electric field, thus to atmospheric electricity
The installation site of field instrument carries out place revision, it is possible to increase the accuracy of lightning monitoring early warning.
It is just at present, less to the correlational study of the impact to atmospheric electric field distortion effect numerical value for the building parameter differing needles,
And not perfect, system, it is badly in need of related content is studied, the conclusion improving, can be effectively to atmospheric electric field
The revision of instrument and the expansion offer scientific guidance that thunderstorm is detected with early warning work.
Content of the invention
In view of the shortcomings of the prior art, it is an object of the present invention to provide a kind of building different parameters are abnormal to atmospheric electric field
The method for numerical simulation of change effect, by numerical simulation is carried out to the distortion effect of building atmospheric electric field, so as to big
The measurement data of gas electric field instrument is revised, it is possible to increase the accuracy of lightning monitoring early warning.
To achieve these goals, the present invention is to realize by the following technical solutions：
A kind of method for numerical simulation to atmospheric electric field distortion effect for the building different parameters of the present invention, walks including following
Suddenly：
(1) obtain the distortion factor of building atmospheric electric field using finite difference method, for the reality to atmospheric electric field detector
Border measured value is revised；
Described finite difference method concrete grammar is as follows：
Sliding-model control is carried out to the spatial mode near-field studied, meets the two dimensional electric field basic representation of Poisson's equation such as
Under：
In formula, ρ is the free charge density value of institute's survey region, and ε is dielectric constant,Potential for institute's survey region；
If under the influence of there is no free charge in the spatial mode near-field studied, taking ρ=0 in above-mentioned formula, then above-mentioned Poisson's equation
With regard to equivalent one-tenth Laplace's equation：
Then by above-mentioned equivalent after Laplace's equation write as the expression-form in rectangular coordinate system：
Then adopt the form of five points difference, sliding-model control is carried out to the spatial mode near-field studied it is assumed that in space
The step-length of the meshing on x and y direction on yardstick is all equal, and step value is h, the difference in the field domain of space that the factor is simulated
Point equation expression formula is：
Wherein, i is the lattice point number in x direction, and j is the lattice point number in y direction；If i=1 is the 1st lattice point in x direction, j=2
I.e. the 2nd lattice point in y direction；
Then, equation (2) and (3) formula are brought in differentiated equation (4) and obtain：
Last it is considered to a series of boundary conditions in spatial simulation region of being located of building, equation (5) is calculated so that it may
To obtain the potential value of each lattice point in the spatial dimension after discretization；
(2) potential value of each lattice point is calculated, during calculating, i, j are former according to from small to large
Then calculated, be then iterated, iterate to setting locus when, calculate newEquation (5) is adopted
With over-relaxation iteration algorithm, and equation is rewritten into：
Wherein, n is iterationses, and ω is overrelaxation iteration factor, using equation (6) just it can be calculated that whole space
The potential value on each lattice point in analog domain, then calculates the electric field intensity on lattice point by calculated potential value,
Relational expression between electric field intensity and current potential is：
I.e. electric field intensity is equal to the gradient of negative current potential；
By the ratio of the electric field intensity value on lattice point and the background electric field value of the spatial mode near-field studied, i.e. distortion factor.
In step (1), the coboundary of analog domain and right boundary meet Neumann boundary condition, ground and building
Thing surface meets Dirichlet boundary condition.
The value of above-mentioned overrelaxation iteration factor ω is 1<ω<2.
Compared with prior art, the invention has the advantages that：Single building is studied respectively by finite difference method
Build the height of height, width and multiple building, the numerical simulation to atmospheric electric field distortion effect for the width of thing, and obtain abnormal
Variable coefficient, for revising to the actual measured value of atmospheric electric field detector, it is possible to increase the accuracy of lightning monitoring early warning.
Brief description
Fig. 1 is the structured flowchart of the present invention；
Fig. 2 is analog domain space structure figure；
Fig. 3 (a) is to atmospheric electric field distortion effect analog result during 20m for depth of building；
Fig. 3 (b) is to atmospheric electric field distortion effect analog result during 25m for depth of building；
Fig. 3 (c) is to atmospheric electric field distortion effect analog result during 30m for depth of building；
Fig. 3 (d) is to atmospheric electric field distortion effect analog result during 35m for depth of building；
Fig. 4 is distortion factor and linear fit result under building differing heights；
Fig. 5 (a) is building width value is to atmospheric electric field distortion effect analog result during 5m；
Fig. 5 (b) is building width value is to atmospheric electric field distortion effect analog result during 7m；
Fig. 5 (c) is building width value is to atmospheric electric field distortion effect analog result during 9m；
Fig. 5 (d) is building width value is to atmospheric electric field distortion effect analog result during 11m；
Fig. 6 is distortion factor and linear fit result under building different in width；
Fig. 7 (a) is positioned at for atmospheric electric field detector and builds value apart from building is to atmospheric electric field distortion effect mould during 2m
Intend result；
Fig. 7 (b) is positioned at for atmospheric electric field detector and builds value apart from building is to atmospheric electric field distortion effect mould during 4m
Intend result；
Fig. 7 (c) is positioned at for atmospheric electric field detector and builds value apart from building is to atmospheric electric field distortion effect mould during 8m
Intend result；
Fig. 7 (d) is positioned at for atmospheric electric field detector and builds value apart from building is to atmospheric electric field distortion effect mould during 16m
Intend result；
Distortion factor under different distance between Fig. 8 ground electric field instrument installation site and building and linear fit knot
Really；
Fig. 9 (a) atmospheric electric field detector is arranged on height and is respectively 35m, 45m, and width is respectively two buildings of 10m, 10m
Atmospheric electric field distortion effect analog result when between building；
Fig. 9 (b) atmospheric electric field detector is arranged on height and is respectively 35m, 70m, and width is respectively two buildings of 10m, 10m
Atmospheric electric field distortion effect analog result when between building；
Fig. 9 (c) atmospheric electric field detector is arranged on height and is respectively 35m, 45m, and width is respectively two buildings of 10m, 6m
Atmospheric electric field distortion effect analog result when between building；
Fig. 9 (d) atmospheric electric field detector is arranged on height and is respectively 35m, 70m, and width is respectively two buildings of 10m, 6m
Atmospheric electric field distortion effect analog result when between building；
Atmospheric electric field distortion effect when Figure 10 (a) is arranged on the low asymmetric buildings thing side of complexity for atmospheric electric field detector
Analog result；
Atmospheric electric field distortion effect when Figure 10 (b) is arranged on the high asymmetric buildings thing side of complexity for atmospheric electric field detector
Analog result.
Specific embodiment
Technological means, creation characteristic, reached purpose and effect for making the present invention realize are easy to understand, with reference to
Specific embodiment, is expanded on further the present invention.
The present invention, when studying the distortion effect research that building causes to atmospheric electric field, needs to obtain its distortion factor table
The ratio of the electric field intensity value being shown as on the lattice point obtained and the background electric field value of spatial mode near-field studied.One kind is provided to build
Build the numerical value to atmospheric electric field distortion effect for the thing different parameters, by Numerical-Mode is carried out to the distortion effect of building atmospheric electric field
Intend, so as to revising to the actual measured results of atmospheric electric field detector.
Referring to Fig. 1, the topmost thought of method of finite difference is exactly that the spatial mode studied near-field is carried out discretization
Process, the present invention can be solved to such electromagnetism problem using finite difference.Meet the two dimensional electric field base of Poisson's equation
This expression formula is as follows：
In formula, ρ is the free charge density value of institute's survey region, and ε is dielectric constant,Potential for institute's survey region.
If under the influence of there is no free charge in the spatial mode near-field studied, taking ρ=0 in above-mentioned formula, then above-mentioned Poisson's equation
With regard to equivalent one-tenth Laplace's equation：
Then by above-mentioned equivalent after Laplace's equation write as the expression-form in rectangular coordinate system：
For two-dimensional problems, then adopt the form of five points difference herein, the spatial mode near-field studied is carried out discrete
Change is processed, and assumes that the step-length of mesh generation on x and y direction on space scale is all equal herein, step value is h, the factor
The difference equation expression formula in the field domain of space simulated is：
Then equation (2) and (3) formula are brought in differentiated equation (4) and obtain：
Last considering some row boundary conditions in spatial simulation region that building is located, equation (5) is calculated so that it may
To obtain the potential value of each lattice point in the spatial dimension after discretization.
The present invention, in the computational engineering to atmospheric electric field distortion effect, has carried out discrete stroke to the space field domain of simulation
Point, then the electric field value of each lattice point is calculated, during calculating, i, j are to enter according to principle from small to large
Row calculates, and is then iterated, when iterating to some locus, calculates newOverrelaxation is adopted to equation (5)
Iterative algorithm, and equation is rewritten into：
In order to ensure the convergence of iteration in calculating process, it is 1 generally to overrelaxation iteration factor value<ω<2.Using
Equation (6) just it can be calculated that potential value on each lattice point in whole spatial mode near-field, then passes through calculated electricity
Place value calculates the electric field intensity on lattice point, and the relational expression between electric field intensity and current potential is：
I.e. electric field intensity is equal to the gradient of negative current potential.The distortion effect in research building, atmospheric electric field being caused herein
During research, need to obtain its distortion factor and can be expressed as the electric field intensity value on the lattice point obtained and the spatial simulation studied
The ratio of the background electric field value in domain.
The present invention, in the research work to the distortion effect of atmospheric electric field, establishes two-dimensional space analog domain, spatial domain
Size be 300m × 100m, be the average of fine day atmospheric electric field to the background electric field value of spatial mode near-field herein, i.e. 130V/
M, fine day atmospheric electric field direction is perpendicularly oriented to the impact not having other free charges in ground, and this research spatial mode near-field, and
Assume that the surface of building and ground are in that equipotential is distributed, that is, potential value is 0.Fig. 2 is the space structure figure of analog domain, in figure
W is the width of building, and H is depth of building.The boundary condition of spatial mode near-field is set to, the coboundary of analog domain and
Right boundary meets Neumann boundary condition；Ground and building surface meet Dirichlet boundary condition.Studied area
The spatial resolution of the spatial mode near-field in domain is that the grid sum that after that is, discretization is taken in space, institute's subdivision goes out is (300/h) *
(100/h).
The present invention studies in building width one timing first, the height of the different buildings distortion effect to atmospheric electric field
Impact.The width of selected building be 20m, then respectively choose building height be 20m, 25m, 30m, 35m this
Four kinds of situations, the size of spatial mode near-field is 300m × 100m.Atmospheric electric field detector is arranged on the top of building, electric field instrument
Height value is 0.5m, and width value is no other free charges in 0.08m, and simulated domain.Using matlab software, write
Go out described above to meet the program under certain boundary conditions, Poisson's equation being solved.Fig. 3 (a), (b), (c), (d)
It is respectively when 20m, 25m, 30m, 35m to atmospheric electric field distortion effect analog result for depth of building, it can be seen that
If simulated domain does not have building, atmospheric electric field should be in an isopotential line, when there is building in simulated domain, building
Thing can cause certain distortion effect to atmospheric electric field, and the more intensive place of contour, and building causes to atmospheric electric field
Distortion effect is more notable.When spatial mode near-field is built-in build thing height be 20m when, the atmospheric electric field distortion factor obtained is
2.0091；As a height of 25m of building, the atmospheric electric field distortion factor obtained is 2.2610；As a height of 30m of building, ask
The atmospheric electric field distortion factor going out is 2.5125；As a height of 30m of building, the atmospheric electric field distortion factor obtained is
2.7638 it can be seen that when being gradually increased with the height built, the atmospheric electric field distortion factor of electric field instrument top end is also therewith
It is gradually increased.
Atmospheric electric field distortion factor under table 1 building differing heights
Table 1 counted building height span between 20m-90m when atmospheric electric field distortion factor statistical result,
It can be seen that under differing heights, atmospheric electric field distortion factor is all higher than 1, further illustrating building all can be to atmospheric electric field
Cause certain distortion effect.Wherein still further it can be seen that, exist certain between atmospheric electric field distortion factor and depth of building
Positive correlation relation, that is, building is higher, and atmospheric electric field distortion factor value is also bigger.To the distortion factor calculating in table 1
As a result, carry out curve fitting, fitting result is as shown in figure 4, it can be seen that atmospheric electric field distortion factor and building
There is significant unary linear relation, the unitary linear equation simulating is y=0.0502*x+1.0058, equation between height
In x be building height, y be atmospheric electric field distortion factor, the increase with x value for the atmospheric electric field distortion factor can be seen
Present the variation tendency of linear increment.Therefore, if it is desired to obtain the air outside research altitude range (20m-90m) herein
During electric field distortion coefficient value, as long as depth of building value is brought into above-mentioned matching equation out.
The present invention secondly research in depth of building one timing, the distortion to atmospheric electric field of the width of the different buildings of research
The impact of effect.The height of selected building be 35m, then respectively choose building width be 5m, 7m, 9m, 11m this
Four kinds of situations, the size of spatial mode near-field is 300m × 100m.Atmospheric electric field detector is arranged on the top centre position of building
Place, the height value of electric field instrument is 0.5m, and width value is to appoint so no other free charges in 0.08m, and simulated domain.Using
Matlab software, writes out and described above meets the program under certain boundary conditions, Poisson's equation being solved.Fig. 5
A (), (b), (c), (d) are depth of building one timing, value is when 5m, 7m, 9m, 11m, atmospheric electric field to be distorted to width respectively
Effect analog result, can intuitively find out from figure, and when the width value of building is different, what atmospheric electric field was caused is abnormal
Change effect also differs, and wherein building width value gets over hour, and the distortion effect that atmospheric electric field is caused is more notable.Work as space
When in analog domain, the width value of building is 5m, the atmospheric electric field distortion factor obtained is 4.9525；When building width takes
When being worth for 7m, the atmospheric electric field distortion factor obtained is 4.2000；When building width value is 9m, the atmospheric electric field obtained
Distortion factor is 3.7511；When building width value is 11m, the atmospheric electric field distortion factor obtained is 3.4521, permissible
Find out, when being gradually increased with the width built, the distortion factor of atmospheric electric field is gradually reduced therewith, on top of building etc.
The more intensive place of gesture line, the distortion effect causing is more notable, the size of electric field value, equal to the gradient of the potential at this position.
Get over the isopotential line of eminence apart from ground, more parallel to ground.
Atmospheric electric field distortion factor under table 2 building different in width
Table 2 counted building width span between 5m-33m when atmospheric electric field distortion factor statistical result,
It can be seen that under different building width values, atmospheric electric field distortion factor is all higher than 1, further illustrates building width
Value all can cause certain distortion effect to atmospheric electric field.Wherein still further it can be seen that, atmospheric electric field distortion factor and building
There is certain negative correlation relation, that is, building is wider, atmospheric electric field distortion factor value is less between width.To in table 2
The distortion factor result calculating, carries out curve fitting, and fitting result is as shown in fig. 6, it can be seen that atmospheric electric field is abnormal
It is nonlinear relation between variable coefficient and the width of building, present the relation of conic section, two Yuan songs simulating
Line equation is y=0.0007*x2- 0.1030*x+6.4071, the x in equation are the width of building, and y distorts for atmospheric electric field
Coefficient, can see that atmospheric electric field distortion factor assumes the variation tendency of decreases in non-linear with the increase of x value.And in building
Width is between 5m-7m, and with the increase of building width, atmospheric electric field distortion factor is in the decline trend of significance, from building
Build thing width increase to 7m after, with building width continuation increase, the amplitude that atmospheric electric field distortion factor successively decreases is relatively
Little.In addition if it is desired to when obtaining the atmospheric electric field distortion factor value outside research width range (5m-33m) herein, as long as will
Building width value brings above-mentioned matching binary curvilinear equation out into.
The present invention continues research in building size one timing, and research atmospheric electric field detector installation site is different apart from building
The distortion effect of atmospheric electric field when at distance.The height of selected building is 35m, and width is 20m.Atmospheric electric field detector is pacified
On the ground, the height value of electric field instrument is 0.5m to dress, and width value is 0.08m, then chooses atmospheric electric field detector respectively and installs
Position is 2m with the distance between building difference value, these four situations of 4m, 8m, 16m, and the size of spatial mode near-field is
300m×100m.Fig. 7 (a), (b), (c), (d) for atmospheric electric field detector be positioned at apart from building build respectively value 2m, 4m,
To atmospheric electric field distortion effect analog result when 8m, 16m, it can be seen that around the position that atmospheric electric field detector is installed
In the presence of having building, building can form a kind of effect of shielding to electric field instrument, and the atmospheric electric field detector calculating is arranged on to be built
The atmospheric electric field distortion factor on electric field instrument top when building near thing is respectively less than 1, and this is mainly due to the shielding action of building
Cause, due to the impact of shielding action, so that the actual measurement measured value of atmospheric electric field detector is less than normal.When ground electric field instrument is installed
Positional distance building apart from value for 2m when, the atmospheric electric field distortion factor on the electric field instrument top obtained is 0.0450；When
Ground electric field instrument installation site apart from building apart from value for 4m when, the electric field instrument top obtained atmospheric electric field distortion system
Number is 0.1794；When ground electric field instrument installation site apart from building apart from value for 8m when, the electric field instrument top obtained
Atmospheric electric field distortion factor is 0.4204；When ground electric field instrument installation site apart from building apart from value for 16m when, obtain
Electric field instrument top atmospheric electric field distortion factor be 0.8282.
Atmospheric electric field distortion factor when table 3 electric field instrument is at building different distance
Table 3 has counted between ground electric field instrument installation site and building horizontal range span between 2m-16m
When electric field instrument top end atmospheric electric field distortion factor statistical result it can be seen that at the position different apart from building,
The atmospheric electric field distortion factor that the surface air electric field instrument top calculating goes out is respectively less than 1 it is generally the case that electric field instrument installation position
Put around no other building when, the distortion factor on the electric field instrument top calculating is 1, and the statistical data from form can enter
One step illustrates, building can form a kind of shielding action to the electric field instrument being arranged on about, so that the reality of electric field instrument
Survey measured value and there is very big error, in the place the closer to building, atmospheric electric field distortion factor value is less, further away from building
Build the place of thing, distortion factor value is bigger, mainly building is less to screen effect remotely.
To the distortion factor result calculating in table 3, carry out curve fitting, fitting result is as shown in figure 8, Fig. 8 is ground
Distortion factor under different distance between electric field instrument installation site and building and linear fit result, can from figure
Go out, ground electric field instrument has positive correlation apart between the horizontal range of building and the atmospheric electric field distortion factor on electric field instrument top
The relation of property, that is, electric field instrument apart from building more remote when, distortion factor value is bigger.It can be seen that ground electric field instrument distance building
It is linear relationship between the atmospheric electric field distortion factor on the horizontal range of thing and electric field instrument top, the unitary linear equation simulating
For y=0.0559*x-0.0397, the x in equation is the horizontal range between ground electric field instrument and building, and y is electric field instrument top
End atmospheric electric field distortion factor, can visually see from equation, atmospheric electric field distortion factor assumes line with the increase of x value
Property be incremented by variation tendency.In addition if it is desired to obtain the atmospheric electric field distortion outside research width range (2m-16m) herein
During coefficient value, as long as the horizontal range value between ground electric field instrument and building is brought into above-mentioned matching unitary straight line side out
Journey.
Multiple buildings to atmospheric electric field distortion effect, with the presence of two or more buildings around the electric field instrument of ground
When, multiple buildings will necessarily to electric field instrument cause a kind of shielding action it is therefore necessary to study multiple buildings to ground
The distortion effect of atmospheric electric field.The size of the spatial mode near-field that the present invention chooses is 300m × 100m, and appoints in simulated domain
So no other free charges.Atmospheric electric field detector is arranged on level ground, and respectively has a solitary building floor in electric field instrument both sides,
Simulate the atmospheric electric field distortion effect under 4 groups of different story heights respectively.First group of situation be：The building of the ground electric field instrument left side
Building thing height is 35m, and building width is 10m, and building is 25m apart from electric field instrument horizontal range, on the right of ground electric field instrument
Depth of building is 45m, and building width is 10m, and building is 20m apart from electric field instrument horizontal range.Second group of situation be：Ground
The building size on the face electric field instrument left side is 70m with the situation of first group of setting, the depth of building on the right of ground electric field instrument, builds
Building thing width is 10m, and building is 20m apart from electric field instrument horizontal range.3rd group of situation be：The building on the ground electric field instrument left side
With the situation of first group of setting, the depth of building on the right of ground electric field instrument is 45m to thing size, and building width is 6m, building
Thing is 22m apart from electric field instrument horizontal range.4th group of situation be：The building size on the ground electric field instrument left side is with first group of setting
Situation, the depth of building on the right of ground electric field instrument be 70m, building width be 6m, building apart from electric field instrument level away from
From for 22m.
Fig. 9 (a), (b), (c), (d) for two buildings that atmospheric electric field detector is arranged under different size build between when big
Pneumoelectric field distortion effect analog result, it can be seen that the effect that different building sizes distorts to atmospheric electric field is
Differ, but on the whole, the more intensive place of contour, the distortion effect that building causes to atmospheric electric field is more notable,
Atmospheric electric field isopotential line in the place more distant place away from building for the vertical dimension is more parallel to ground.For the first situation,
The distortion factor to the electric field instrument top being arranged between both size buildings that we go out with matlab computed in software
It is worth for 0.8303, distortion factor is less than 1, the screen effect mainly due to building causes；Second group of situation, only changes
The height of right side building, the atmospheric electric field distortion factor value on the electric field instrument top calculating is also 0.8303, and this illustrates
In the case of the building size constancy of left side, only change the height of right side building, the electric field distortion value calculating is identical,
This may be nearer apart from left side building mainly due to electric field instrument, and electric field instrument is within the shield ranges of left side building；
For the 3rd group of situation, the electric field distortion coefficient of the electric field instrument top end of calculating is 0.8644, and this distortion result is greater than the
First, two groups of situation, mainly the 3rd group experiment reduces the width of right side building, and narrower building is made to atmospheric electric field
The distortion effect becoming is more notable；The electric field distortion coefficient of the electric field instrument top end calculating in the case of the 4th group is also 0.8644.
It is abnormal to atmospheric electric field that the present invention continues have studied the building when atmospheric electric field detector is arranged on asymmetric buildings side
The impact of change effect, and the screen effect to ground electric field instrument.Figure 10 (a), (b) are arranged on for atmospheric electric field detector and irregularly build
Atmospheric electric field distortion effect analog result when building thing side, it can be seen that abnormal in building more most advanced and sophisticated place atmospheric electric field
Change effect is more notable, and in figure simulates two kinds of complex building situations, and the complex building calculating is to ground atmospheric electric field detector
Distortion factor be respectively 0.4611,0.4453.
Ultimate principle and principal character and the advantages of the present invention of the present invention have been shown and described above.The technology of the industry
, it should be appreciated that the present invention is not restricted to the described embodiments, the simply explanation described in above-described embodiment and description is originally for personnel
The principle of invention, without departing from the spirit and scope of the present invention, the present invention also has various changes and modifications, these changes
Change and improvement both falls within scope of the claimed invention.Claimed scope by appending claims and its
Equivalent thereof.
Claims (3)
1. a kind of building different parameters to the method for numerical simulation of atmospheric electric field distortion effect it is characterised in that including following
Step：
(1) obtain the distortion factor of building atmospheric electric field using finite difference method, for the actual survey to atmospheric electric field detector
Value is revised；
Described finite difference method concrete grammar is as follows：
Sliding-model control is carried out to the spatial mode near-field studied, the two dimensional electric field basic representation meeting Poisson's equation is as follows：
In formula, ρ is the free charge density value of institute's survey region, and ε is dielectric constant,Potential for institute's survey region；If institute
Under the influence of there is no free charge in the spatial mode near-field of research, in above-mentioned formula, take ρ=0, then above-mentioned Poisson's equation is just equivalent
Become Laplace's equation：
Then by above-mentioned equivalent after Laplace's equation write as the expression-form in rectangular coordinate system：
Then adopt the form of five points difference, sliding-model control is carried out to the spatial mode near-field studied it is assumed that in space scale
The step-length of the meshing on upper x and y direction is all equal, and step value is h, the difference side in the field domain of space that the factor is simulated
Journey expression formula is：
Wherein, i is the lattice point number in x direction, and j is the lattice point number in y direction；
Then, equation (2) and (3) formula are brought in differentiated equation (4) and obtain：
Finally it is considered to the boundary condition in the spatial simulation region at building place, equation (5) is calculated it is possible to obtain discrete
The potential value of each lattice point in spatial dimension after change；
(2) potential value of each lattice point is calculated, during calculating, i, j are to enter according to principle from small to large
Row calculate, be then iterated, iterate to setting locus when, calculate newTo equation (5) using super
Relaxative iteration algorithm, and equation is rewritten into：
Wherein, n is iterationses, and ω is overrelaxation iteration factor, using equation (6) just it can be calculated that whole spatial simulation
The potential value on each lattice point in domain, then calculates the electric field intensity on lattice point, electric field by calculated potential value
Relational expression between intensity and current potential is：
I.e. electric field intensity is equal to the gradient of negative current potential；
By the ratio of the electric field intensity value on lattice point and the background electric field value of the spatial mode near-field studied, i.e. distortion factor.
2. the method for numerical simulation to atmospheric electric field distortion effect for the building different parameters according to claim 1, it is special
Levy and be, in step (1), the coboundary of analog domain and right boundary meet Neumann boundary condition, ground and building
Surface meets Dirichlet boundary condition.
3. the method for numerical simulation to atmospheric electric field distortion effect for the building different parameters according to claim 1, it is special
Levy and be, the value of described overrelaxation iteration factor ω is 1<ω<2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610810981.5A CN106407541A (en) | 2016-09-08 | 2016-09-08 | Numerical simulation method for atmospheric electric field distortion effect of buildings under different parameters |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610810981.5A CN106407541A (en) | 2016-09-08 | 2016-09-08 | Numerical simulation method for atmospheric electric field distortion effect of buildings under different parameters |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106407541A true CN106407541A (en) | 2017-02-15 |
Family
ID=57999415
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610810981.5A Pending CN106407541A (en) | 2016-09-08 | 2016-09-08 | Numerical simulation method for atmospheric electric field distortion effect of buildings under different parameters |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106407541A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112067908A (en) * | 2020-08-20 | 2020-12-11 | 国网山东省电力公司电力科学研究院 | Fitting method and system for distortion electric field when transformer substation robot measures power frequency electric field |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103926475A (en) * | 2014-04-08 | 2014-07-16 | 沈阳工业大学 | Electric field measuring method and measuring system thereof |
| JP5652568B1 (en) * | 2014-07-23 | 2015-01-14 | Ｊｆｅスチール株式会社 | Manufacturing method of ferritic stainless steel foil for solar cell substrate |
| CN104880650A (en) * | 2015-05-26 | 2015-09-02 | 武汉大学 | Air gap breakdown voltage prediction method |
| CN104913847A (en) * | 2015-06-04 | 2015-09-16 | 深圳大学 | Device and method for measuring high-intensity terahertz time-domain spectroscopy |
| CN105203856A (en) * | 2015-09-28 | 2015-12-30 | 国网重庆市电力公司电力科学研究院 | Electric field intensity calculating method and device of oil paper insulation system |
-
2016
- 2016-09-08 CN CN201610810981.5A patent/CN106407541A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103926475A (en) * | 2014-04-08 | 2014-07-16 | 沈阳工业大学 | Electric field measuring method and measuring system thereof |
| JP5652568B1 (en) * | 2014-07-23 | 2015-01-14 | Ｊｆｅスチール株式会社 | Manufacturing method of ferritic stainless steel foil for solar cell substrate |
| CN104880650A (en) * | 2015-05-26 | 2015-09-02 | 武汉大学 | Air gap breakdown voltage prediction method |
| CN104913847A (en) * | 2015-06-04 | 2015-09-16 | 深圳大学 | Device and method for measuring high-intensity terahertz time-domain spectroscopy |
| CN105203856A (en) * | 2015-09-28 | 2015-12-30 | 国网重庆市电力公司电力科学研究院 | Electric field intensity calculating method and device of oil paper insulation system |
Non-Patent Citations (2)
| Title |
|---|
| 吴坚铃: "建筑物对大气电场畸变效应的数值模拟", 《广东气象》 * |
| 郭秀峰: "建筑物尖端对大气电场畸变影响的数值计算", 《应用气象学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112067908A (en) * | 2020-08-20 | 2020-12-11 | 国网山东省电力公司电力科学研究院 | Fitting method and system for distortion electric field when transformer substation robot measures power frequency electric field |
| CN112067908B (en) * | 2020-08-20 | 2023-06-16 | 国网山东省电力公司电力科学研究院 | Method and system for fitting distorted electric field during power frequency electric field measurement by substation robot |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103020342B (en) | Method for extracting contour and corner of building from ground LiDAR data | |
| CN107704641A (en) | Fine simulation of wind method based on outdoor scene vegetation spatial distribution roughness | |
| KR101332559B1 (en) | Numerical simulation system and method for atmospheric wind flow by computational fluid dynamics | |
| CN104463700A (en) | Power transmission line tower lightning strike risk evaluation method based on data mining technology | |
| CN110598290A (en) | Method and system for predicting future hydropower generation capacity of basin considering climate change | |
| CN105576650A (en) | Wind power prediction method and system based on assimilation technology and same wind band theory | |
| CN112381358A (en) | Near-real-time assessment method for wind disaster damage risk of greening trees facing urban area | |
| CN108509718B (en) | Far-field wake two-dimensional analytic model based on mass conservation | |
| Hågbo et al. | Pedestrian wind comfort assessment using computational fluid dynamics simulations with varying number of wind directions | |
| CN106407541A (en) | Numerical simulation method for atmospheric electric field distortion effect of buildings under different parameters | |
| CN104951999A (en) | Mountain land photovoltaic power station wind pressure computing method based on landform fitting and CFD (Computational Fluid Dynamics) | |
| CN102236746B (en) | Wind resource simulated estimation method for region without wind measurement records | |
| CN109344356B (en) | Method for determining direct lightning protection efficiency of high-magnitude lightning protection system | |
| CN114139391A (en) | Wind energy resource assessment method and device based on separation vortex turbulence model and storage medium | |
| CN112966404B (en) | Method for generating three-dimensional lightning leader development path | |
| CN109670229B (en) | Method for estimating ground lightning falling density near overhead bridge overhead contact system of high-speed railway | |
| Guo et al. | Improvement of corona discharge model and its application on simulating corona discharge in the presence of wind | |
| CN107886573A (en) | A kind of complex geological condition slope three-dimensional finite element mesh generation method | |
| CN103399980B (en) | A kind of optimized Selection collocation method of power transmission and transformation insulator | |
| CN108764602B (en) | Building group natural ventilation potential evaluation method and system | |
| Perera et al. | Fractal nature of simulated lightning channels | |
| JP4043779B2 (en) | Basic wind speed map creation method by wind direction and basic wind speed estimation method by wind direction | |
| CN107315855A (en) | Wind power plant turbulent flow optimization method and system | |
| CN114021303A (en) | Urban ventilation corridor intelligent excavation method based on high-precision oblique photographic image | |
| Adekitan et al. | A further look at dynamic electro-geometrical model: Its fundamentals and implementation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170215 |
|
| RJ01 | Rejection of invention patent application after publication |