CN110414053A - A kind of time-domain numerical simulation method of quick determining component micro-discharge threshold - Google Patents
A kind of time-domain numerical simulation method of quick determining component micro-discharge threshold Download PDFInfo
- Publication number
- CN110414053A CN110414053A CN201910569890.0A CN201910569890A CN110414053A CN 110414053 A CN110414053 A CN 110414053A CN 201910569890 A CN201910569890 A CN 201910569890A CN 110414053 A CN110414053 A CN 110414053A
- Authority
- CN
- China
- Prior art keywords
- micro
- discharge
- numerical simulation
- macro
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 111
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000002245 particle Substances 0.000 claims abstract description 106
- 239000000463 material Substances 0.000 claims abstract description 35
- 239000003574 free electron Substances 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 30
- 238000006073 displacement reaction Methods 0.000 claims description 15
- 230000005672 electromagnetic field Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 5
- 230000035699 permeability Effects 0.000 claims description 3
- 208000028659 discharge Diseases 0.000 description 117
- 238000013461 design Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000000126 substance Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of time-domain numerical simulation methods of quick determining component micro-discharge threshold, comprising: establishes the 3-D geometric model of description component minimum dimension structure, and the 3-D geometric model is split into multiple hexahedral mesh;Determine the boundary of 3-D geometric model and the secondary electron emission characteristic parameter of material properties and material;Using the free electron in macro particle simulation space, primary power, the quantity of electric charge, quality, initial distribution state and the initial motion direction of defmacro particle;Macro particle chained list is established, and carries out macro particle chained list initialization;Based on 3-D geometric model, the boundary of 3-D geometric model and material properties, the secondary electron emission characteristic parameter of material and macro particle chained list, from micro-discharge threshold scan power initial value P0Start, carries out micro-discharge threshold numerical simulation.The simulation efficiency of micro discharge numerical simulation is improved through the invention.
Description
Technical Field
The invention belongs to the technical field of satellite-borne microwave components, and particularly relates to a time domain numerical simulation method for rapidly determining a micro-discharge threshold of a component.
Background
The high micro-discharge risk of the high-power microwave component of the spacecraft is a key factor influencing the long service life and high reliability of the effective load of the spacecraft, and is also the largest single-point failure link of the satellite under the high-power application. For the micro-discharge design of a high-power microwave component, effective numerical simulation and threshold analysis are key steps for ground verification, and the problems that the development period is too long and even the maximum working power reliability design cannot be obtained due to experiments and repeated design change can be solved.
The micro-discharge three-dimensional numerical simulation technology based on full-wave electromagnetic calculation and particle trajectory propulsion is the most advanced micro-discharge analysis technology in the world at present, is a key technology for verifying the micro-discharge mechanism, is a necessary means for optimizing and overcoming the inexhaustible experimental parameters by the micro-discharge inhibition technology, and is a necessary tool for designing a high-power microwave component of a spacecraft.
For complex microwave components, how to improve the micro-discharge effect numerical simulation speed of a spacecraft high-power microwave device and reduce the actual micro-discharge threshold prediction time of the microwave device is a problem to be solved urgently by technical personnel in the field.
Disclosure of Invention
The technical problem of the invention is solved: the method overcomes the defects of the prior art, and provides a time domain numerical simulation method for rapidly determining the micro-discharge threshold of the component so as to improve the micro-discharge numerical simulation efficiency.
In order to solve the technical problem, the invention discloses a time domain numerical simulation method for rapidly determining a micro-discharge threshold of a component, which comprises the following steps:
establishing a three-dimensional geometric model describing the minimum-size structure of a component, and dividing the three-dimensional geometric model into a plurality of hexahedral meshes;
determining the boundary and material attribute of the three-dimensional geometric model and the secondary electron emission characteristic parameter of the material;
adopting macro particles to simulate free electrons in space, and defining initial energy, charge quantity, mass, initial distribution state and initial movement direction of the macro particles; establishing a macro particle chain table, and initializing the macro particle chain table;
scanning a power starting value P from a micro-discharge threshold value based on a three-dimensional geometric model, a boundary and material property of the three-dimensional geometric model, a secondary electron emission characteristic parameter of a material, and a macro-particle chain table0Initially, a microdischarge threshold value numerical simulation was performed.
In the time domain numerical simulation method for rapidly determining the micro-discharge threshold of the component, the material property includes: dielectric constant εrMagnetic permeability murAnd the electrical conductivity σ.
In the time domain numerical simulation method for rapidly determining the micro-discharge threshold of the component, the secondary electron emission characteristic parameter comprises: secondary electron emission coefficient deltase、δseCorresponding incident electron energy value EmIncident electron energy value E corresponding to a secondary electron emission coefficient of 11And E2。
In the time domain numerical simulation method for quickly determining the micro-discharge threshold of the component, the macro particle sublist records: the motion velocity vector and the displacement vector of each macro-particle.
In the time domain numerical simulation method for quickly determining the micro-discharge threshold of the component, the initialization of the macro-particle chain table includes: determining that initial displacement vectors of macro particles are randomly distributed in a micro-discharge numerical simulation area, and the initial movement speed is zero; determining the total macro particle number in all hexahedron grids as N0。
In the time domain numerical simulation method for quickly determining the micro-discharge threshold of the component, the method further includes: defining: a hollow area in the three-dimensional geometric model is a micro-discharge numerical simulation area; the material of the micro-discharge numerical simulation area is vacuum; the micro-discharge numerical simulation time step is N delta t; the initial value of micro-discharge threshold scanning power is P0The power scanning precision is a, and the micro-discharge numerical simulation time is T1(ii) a Wherein, DeltatAnd N is an integer greater than or equal to 1.
Time domain number of micro-discharge threshold of the rapid determination componentThe value simulation method further comprises the following steps: working center frequency f according to three-dimensional geometric model0And obtaining the geometric characteristic of an electromagnetic wave input port of the three-dimensional geometric model by solving Maxwell equations to obtain a time period T0Electric field distribution E in all hexahedral meshes in the internal three-dimensional geometric modeliAnd magnetic field distribution Bi(ii) a Wherein, T0=1/f0N, n is the number of hexahedral meshes.
In the time domain numerical simulation method for rapidly determining the micro-discharge threshold of the component, the initial power value P is scanned from the micro-discharge threshold based on the three-dimensional geometric model, the boundary and material properties of the three-dimensional geometric model, the secondary electron emission characteristic parameter of the material, and the macro-particle chain table0Initially, a microdischarge threshold value simulation is performed, comprising:
(1) determining a power scaling factor xp(ii) a Wherein the content of the first and second substances,p represents micro-discharge numerical simulation input power;
(2) the one time period T0Electric field distribution E in all hexahedral meshes in the internal three-dimensional geometric modeliAnd magnetic field distribution BiAre respectively multiplied by power scale factors xpCirculating according to periods to obtain electromagnetic field distribution at any time in all hexahedral meshes when the micro-discharge numerical simulation input power P is obtained;
(3) according to the micro-discharge numerical simulation time step N delta t, performing linear interpolation on the electromagnetic field distribution at any time in all hexahedral grids when the micro-discharge numerical simulation input power P is input, and obtaining the electromagnetic field at the position of each macro particle; performing iterative calculation by using a Lorentz equation system for describing the movement of the macro particles to obtain a macro particle movement velocity vector and a displacement vector of each micro-discharge numerical simulation time step; judging whether the macro particles collide with the boundary of the three-dimensional geometric model in the time step according to the coordinates of the macro particles; if the collision is determined, deleting all the collided macro particles from the macro particle sublist, and not counting the calculation of the next time step; judging whether macro particles are emitted from the material or not according to the motion velocity vector of the macro particles and the secondary electron emission characteristic parameter of the material of the three-dimensional geometric model at the collision position; if the macro particles are determined to be emitted from the material, determining the motion velocity vector and the displacement vector of the emitted macro particles, and taking the determined motion velocity vector and the determined displacement vector of the emitted macro particles as the initial state of the macro particles in the next time step; recording the total number Nx of macro particles in all hexahedron grids in the time step, and updating a macro particle chain table;
(4) simulating a time step according to the micro-discharge numerical value to carry out propulsion, recording the change of the total macro-particle number in all hexahedral meshes in all time steps along with the time, and repeating the step (3); when the total number of macro particles in all hexahedron grids is larger than the initial value N of the total number of macro particles in a certain time step010 of3When the micro discharge threshold value is doubled, ending the micro discharge threshold value numerical simulation; or ending the micro-discharge threshold value numerical simulation when the micro-discharge numerical simulation time reaches the set simulation time threshold value.
In the time domain numerical simulation method for quickly determining the micro-discharge threshold of the component, the method further includes:
defining the number of macro particles in the c-th radio frequency period as Nc, and if the micro-discharge numerical simulation time exceeds c radio frequency periods, each time period T till the micro-discharge numerical simulation0The total number of macro particles in all the hexahedral meshes is less than the next time period T0And determining that micro-discharge occurs if the total number of macro-particles in all the hexahedral meshes is equal to or less than the total number of macro-particles in all the hexahedral meshes, and otherwise, determining that micro-discharge does not occur.
In the time domain numerical simulation method for quickly determining the micro-discharge threshold of the component, the method further includes:
when the micro-discharge numerical simulation input power P is the micro-discharge threshold value scanning power initial value P0If micro-discharge is determined to occur, then according to the power scaling factor PS1Scanning the micro-discharge threshold with the initial power value P0Is set to be Ps1P0Repeating steps (1) to (1)4) Carrying out micro-discharge threshold value numerical simulation until determining that micro-discharge does not occur, and recording the non-discharge power P at the momentm-1And, no discharge power Pm-1Last input power P ofm;
When the micro-discharge numerical simulation input power P is the micro-discharge threshold value scanning power initial value P0If it is determined that no micro-discharge has occurred, the power is increased by a power increase factor Pz1Scanning the micro-discharge threshold with the initial power value P0Is set to be Pz1P0Repeating the steps (1) to (4), carrying out micro-discharge threshold value numerical simulation until micro-discharge is determined to occur, and recording the discharge power at the momentAnd, discharge powerLast input power ofIf it isThen, the micro-discharge numerical value is set as the analog input power PRepeating the steps (1) to (4) to carry out micro-discharge threshold value numerical simulation untilWill be provided withDetermined at the timeThe value of (d) is used as the microdischarge threshold.
The invention has the following advantages:
the invention discloses a time domain numerical simulation method for rapidly determining a micro-discharge threshold of a component, aiming at the characteristics that the number of particles at the initial stage of micro-discharge is slowly increased and the influence on an electromagnetic field is small, firstly, the device is subjected to primary electromagnetic field solving, and then, the solved electromagnetic field result is repeatedly utilized in micro-discharge simulation under different powers, so that the time for calculating the micro-discharge threshold is greatly shortened, and the simulation efficiency is improved.
The method automatically judges the micro-discharge threshold value through the change characteristic of the particle number curve, realizes quick and effective micro-discharge numerical simulation and threshold value determination, is widely applied to micro-discharge analysis of various high-power microwave components in the follow-up process, and has wide potential market application prospect and huge economic benefit.
The method is beneficial to improving the micro-discharge effect numerical simulation speed of the high-power microwave device (complex microwave component) of the spacecraft, shortening the micro-discharge threshold prediction time of the actual microwave device, improving the practicability and the usability of MSAT software and providing more powerful support for the engineering design of the micro-discharge effect of the high-power microwave device of the spacecraft.
Drawings
FIG. 1 is a flow chart illustrating steps of a time domain numerical simulation method for rapidly determining a micro-discharge threshold of a component according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a time step of a three-dimensional geometric model according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the embodiments of the present invention will be described in detail with reference to the accompanying drawings.
As shown in fig. 1, in this embodiment, the time domain numerical simulation method for quickly determining the micro-discharge threshold of the component includes:
step 101, establishing a three-dimensional geometric model describing a minimum-size structure of a component, and dividing the three-dimensional geometric model into a plurality of hexahedral meshes.
In the present embodiment, the frequency f of the working center can be determined according to the three-dimensional geometric model0And obtaining an electromagnetic wave input port geometric characteristic of the three-dimensional geometric model by solving Maxwell equationsTime period T0Electric field distribution E in all hexahedral meshes in the internal three-dimensional geometric modeliAnd magnetic field distribution Bi. Wherein, T0=1/f0N, n is the number of hexahedral meshes.
Preferably, a time period T0Electric field distribution E in all hexahedral meshes in the internal three-dimensional geometric modeliAnd magnetic field distribution BiThe solution flow of (c) may be as follows:
working center frequency f according to three-dimensional geometric model0And the geometric characteristics of the electromagnetic wave input port of the three-dimensional geometric model, solving Maxwell equations to obtain the initial value of the electromagnetic field in each hexahedral mesh on the electromagnetic wave input port of the three-dimensional geometric model when the input power is 1W, continuously solving Maxwell equations describing the electromagnetic product evolution in each hexahedral mesh to obtain a time period T0Electric field distribution E in all hexahedral meshes in the internal three-dimensional geometric modeliAnd magnetic field distribution Bi。
And 102, determining the boundary and material properties of the three-dimensional geometric model and secondary electron emission characteristic parameters of the material.
In the present embodiment, the material properties include, but are not limited to: dielectric constant εrMagnetic permeability murAnd the electrical conductivity σ. Secondary electron emission characteristic parameters include, but are not limited to: : secondary electron emission coefficient deltase、δseCorresponding incident electron energy value EmIncident electron energy value E corresponding to a secondary electron emission coefficient of 11And E2。
In this embodiment, the following definitions can be made: a hollow area in the three-dimensional geometric model is a micro-discharge numerical simulation area; the material of the micro-discharge numerical simulation area is vacuum; the micro-discharge numerical simulation time step is N delta t; the initial value of micro-discharge threshold scanning power is P0The power scanning precision is a, and the micro-discharge numerical simulation time is T1. In fig. 2, Δ t represents a time step of the three-dimensional geometric model, and N is an integer greater than or equal to 1.
103, adopting macro particles to simulate free electrons in space, and defining initial energy, charge quantity, mass, initial distribution state and initial movement direction of the macro particles; and establishing a macro particle chain table and initializing the macro particle chain table.
In this embodiment, the macro particle sublist has recorded therein: the motion velocity vector and the displacement vector of each macro-particle. Wherein, the macro particle chain table initialization comprises: determining that initial displacement vectors of macro particles are randomly distributed in a micro-discharge numerical simulation area, and the initial movement speed is zero; determining the total macro particle number in all hexahedron grids as N0。
Step 104, scanning the initial power value P from the micro-discharge threshold value based on the three-dimensional geometric model, the boundary and material attribute of the three-dimensional geometric model, the secondary electron emission characteristic parameter of the material and the macro-particle chain table0Initially, a microdischarge threshold value numerical simulation was performed.
In this embodiment, the step 104 may specifically include:
(1) determining a power scaling factor xp。
Wherein the content of the first and second substances,p represents the micro-discharge numerical analog input power.
(2) The one time period T0Electric field distribution E in all hexahedral meshes in the internal three-dimensional geometric modeliAnd magnetic field distribution BiAre respectively multiplied by power scale factors xpAnd circulating according to periods to obtain the electromagnetic field distribution at any time in all hexahedral meshes when the micro-discharge numerical simulation input power P is obtained.
(3) According to the micro-discharge numerical simulation time step N delta t, performing linear interpolation on the electromagnetic field distribution at any time in all hexahedral grids when the micro-discharge numerical simulation input power P is input, and obtaining the electromagnetic field at the position of each macro particle; performing iterative calculation by using a Lorentz equation system for describing the movement of the macro particles to obtain a macro particle movement velocity vector and a displacement vector of each micro-discharge numerical simulation time step; judging whether the macro particles collide with the boundary of the three-dimensional geometric model in the time step according to the coordinates of the macro particles; if the collision is determined, deleting all the collided macro particles from the macro particle sublist, and not counting the calculation of the next time step; judging whether macro particles are emitted from the material or not according to the motion velocity vector of the macro particles and the secondary electron emission characteristic parameter of the material of the three-dimensional geometric model at the collision position; if the macro particles are determined to be emitted from the material, determining the motion velocity vector and the displacement vector of the emitted macro particles, and taking the determined motion velocity vector and the determined displacement vector of the emitted macro particles as the initial state of the macro particles in the next time step; and recording the total number Nx of macro particles in all hexahedron grids in the time step, and updating the macro particle chain table.
(4) Simulating a time step according to the micro-discharge numerical value to carry out propulsion, recording the change of the total macro-particle number in all hexahedral meshes in all time steps along with the time, and repeating the step (3); when the total number of macro particles in all hexahedron grids is larger than the initial value N of the total number of macro particles in a certain time step010 of3When the micro discharge threshold value is doubled, ending the micro discharge threshold value numerical simulation; or ending the micro-discharge threshold value numerical simulation when the micro-discharge numerical simulation time reaches the set simulation time threshold value.
In a preferred embodiment of the present invention, the time domain numerical simulation method for rapidly determining the micro-discharge threshold of the component may further include:
defining the number of macro particles in the c-th radio frequency period as Nc, and if the micro-discharge numerical simulation time exceeds c radio frequency periods, each time period T till the micro-discharge numerical simulation0The total number of macro particles in all the hexahedral meshes is less than the next time period T0And determining that micro-discharge occurs if the total number of macro-particles in all the hexahedral meshes is equal to or less than the total number of macro-particles in all the hexahedral meshes, and otherwise, determining that micro-discharge does not occur.
In a preferred embodiment of the present invention, the time domain numerical simulation method for rapidly determining the micro-discharge threshold of the component may further include:
when the micro-discharge numerical simulation input power P is the micro-discharge threshold value scanning power initial value P0If micro-discharge is determined to occur, then according to the power scaling factor PS1Scanning the micro-discharge threshold with the initial power value P0Is set to be Ps1P0Repeating the steps (1) to (4), carrying out micro-discharge threshold value numerical simulation until determining that micro-discharge does not occur, and recording the non-discharge power P at the momentm-1And, no discharge power Pm-1Last input power P ofm。
When the micro-discharge numerical simulation input power P is the micro-discharge threshold value scanning power initial value P0If it is determined that no micro-discharge has occurred, the power is increased by a power increase factor Pz1Scanning the micro-discharge threshold with the initial power value P0Is set to be Pz1P0Repeating the steps (1) to (4), carrying out micro-discharge threshold value numerical simulation until micro-discharge is determined to occur, and recording the discharge power at the momentAnd, discharge powerLast input power ofIf it isThen, the micro-discharge numerical value is set as the analog input power PRepeating the steps (1) to (4) to carry out micro-discharge threshold value numerical simulation untilWill be provided withDetermined at the timeThe value of (d) is used as the microdischarge threshold.
In this embodiment, as shown in table 1, the threshold accuracy and the calculation efficiency at different time steps are shown:
TABLE 1
Therefore, the time domain numerical simulation method for rapidly determining the micro-discharge threshold of the component improves the micro-discharge numerical simulation efficiency.
The embodiments in the present description are all described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The above description is only for the best mode of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention.
Those skilled in the art will appreciate that the invention may be practiced without these specific details.
Claims (10)
1. A time domain numerical simulation method for rapidly determining a micro-discharge threshold of a component is characterized by comprising the following steps:
establishing a three-dimensional geometric model describing the minimum-size structure of a component, and dividing the three-dimensional geometric model into a plurality of hexahedral meshes;
determining the boundary and material attribute of the three-dimensional geometric model and the secondary electron emission characteristic parameter of the material;
adopting macro particles to simulate free electrons in space, and defining initial energy, charge quantity, mass, initial distribution state and initial movement direction of the macro particles; establishing a macro particle chain table, and initializing the macro particle chain table;
based on three-dimensional tableWhich model, the boundary and material properties of the three-dimensional geometric model, the secondary electron emission characteristic parameters of the material, and the macro-particle chain table, scanning the power starting value P from the micro-discharge threshold0Initially, a microdischarge threshold value numerical simulation was performed.
2. The time-domain numerical simulation method for rapidly determining a component microdischarge threshold of claim 1, wherein the material property comprises: dielectric constant εrMagnetic permeability murAnd the electrical conductivity σ.
3. The time-domain numerical simulation method for rapidly determining a component microdischarge threshold of claim 2, wherein the secondary electron emission characteristic parameters comprise: secondary electron emission coefficient deltase、δseCorresponding incident electron energy value EmIncident electron energy value E corresponding to a secondary electron emission coefficient of 11 and E2。
4. A time domain numerical simulation method for rapidly determining a component microdischarge threshold as claimed in claim 3, wherein the macro particle chain table has recorded therein: the motion velocity vector and the displacement vector of each macro-particle.
5. A time-domain numerical simulation method for rapidly determining a component microdischarge threshold of claim 4, wherein the macro-particle chain table initialization comprises: determining that initial displacement vectors of macro particles are randomly distributed in a micro-discharge numerical simulation area, and the initial movement speed is zero; determining the total macro particle number in all hexahedron grids as N0。
6. The time-domain numerical simulation method for rapidly determining a component microdischarge threshold of claim 5, further comprising: defining: a hollow area in the three-dimensional geometric model is a micro-discharge numerical simulation area; the material of the micro-discharge numerical simulation area is vacuum; the micro-discharge numerical simulation time step is N delta t; micro-discharge threshold scanningPower starting value is P0The power scanning precision is a, and the micro-discharge numerical simulation time is T1(ii) a Wherein, Δ t represents the time step of the three-dimensional geometric model, and N is an integer greater than or equal to 1.
7. The time-domain numerical simulation method for rapidly determining a component microdischarge threshold of claim 6, further comprising: working center frequency f according to three-dimensional geometric model0And obtaining the geometric characteristic of an electromagnetic wave input port of the three-dimensional geometric model by solving Maxwell equations to obtain a time period T0Electric field distribution E in all hexahedral meshes in the internal three-dimensional geometric modeliAnd magnetic field distribution Bi; wherein ,T0=1/f0And i is 1, …, and n is the number of hexahedral meshes.
8. The time-domain numerical simulation method for rapidly determining a micro-discharge threshold of a component of claim 7, wherein the power starting value P is scanned from the micro-discharge threshold based on a three-dimensional geometric model, boundary and material properties of the three-dimensional geometric model, secondary electron emission characteristic parameters of the material, and a macro-particle sublist0Initially, a microdischarge threshold value simulation is performed, comprising:
(1) determining a power scaling factor xp; wherein ,p represents micro-discharge numerical simulation input power;
(2) the one time period T0Electric field distribution E in all hexahedral meshes in the internal three-dimensional geometric modeliAnd magnetic field distribution BiAre respectively multiplied by power scale factors xpCirculating according to periods to obtain electromagnetic field distribution at any time in all hexahedral meshes when the micro-discharge numerical simulation input power P is obtained;
(3) according to the micro-discharge numerical simulation time step N delta t, performing linear interpolation on the electromagnetic field distribution at any time in all hexahedral grids when the micro-discharge numerical simulation input power P is input, and obtaining the electromagnetic field at the position of each macro particle; performing iterative calculation by using a Lorentz equation system for describing the movement of the macro particles to obtain a macro particle movement velocity vector and a displacement vector of each micro-discharge numerical simulation time step; judging whether the macro particles collide with the boundary of the three-dimensional geometric model in the time step according to the coordinates of the macro particles; if the collision is determined, deleting all the collided macro particles from the macro particle sublist, and not counting the calculation of the next time step; judging whether macro particles are emitted from the material or not according to the motion velocity vector of the macro particles and the secondary electron emission characteristic parameter of the material of the three-dimensional geometric model at the collision position; if the macro particles are determined to be emitted from the material, determining the motion velocity vector and the displacement vector of the emitted macro particles, and taking the determined motion velocity vector and the determined displacement vector of the emitted macro particles as the initial state of the macro particles in the next time step; recording the total number Nx of macro particles in all hexahedron grids in the time step, and updating a macro particle chain table;
(4) simulating a time step according to the micro-discharge numerical value to carry out propulsion, recording the change of the total macro-particle number in all hexahedral meshes in all time steps along with the time, and repeating the step (3); when the total number of macro particles in all hexahedron grids is larger than the initial value N of the total number of macro particles in a certain time step010 of3When the micro discharge threshold value is doubled, ending the micro discharge threshold value numerical simulation; or ending the micro-discharge threshold value numerical simulation when the micro-discharge numerical simulation time reaches the set simulation time threshold value.
9. The time-domain numerical simulation method for rapidly determining a component microdischarge threshold of claim 8, further comprising:
defining the number of macro particles in the c-th radio frequency period as Nc, and if the micro-discharge numerical simulation time exceeds c radio frequency periods, each time period T till the micro-discharge numerical simulation0The total number of macro particles in all the hexahedral meshes is less than the next time period T0All six inAnd determining that the micro-discharge occurs if the total number of macro-particles in the surface body grid is less than the total number of macro-particles in the surface body grid, and otherwise, determining that the micro-discharge does not occur.
10. The time-domain numerical simulation method for rapidly determining a component microdischarge threshold of claim 9, further comprising:
when the micro-discharge numerical simulation input power P is the micro-discharge threshold value scanning power initial value P0If micro-discharge is determined to occur, then according to the power scaling factor PS1Scanning the micro-discharge threshold with the initial power value P0Is set to be Ps1P0Repeating the steps (1) to (4), carrying out micro-discharge threshold value numerical simulation until determining that micro-discharge does not occur, and recording the non-discharge power P at the momentm-1And, no discharge power Pm-1Last input power P ofm;
When the micro-discharge numerical simulation input power P is the micro-discharge threshold value scanning power initial value P0If it is determined that no micro-discharge has occurred, the power is increased by a power increase factor Pz1Scanning the micro-discharge threshold with the initial power value P0Is set to be Pz1P0Repeating the steps (1) to (4), carrying out micro-discharge threshold value numerical simulation until micro-discharge is determined to occur, and recording the discharge power at the momentAnd, discharge powerLast input power ofIf it isThen, the micro-discharge numerical value is set as the analog input power PRepeating the stepsStep (1) to step (4), micro-discharge threshold value numerical simulation is carried out untilWill be provided withDetermined at the timeThe value of (d) is used as the microdischarge threshold.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910569890.0A CN110414053B (en) | 2019-06-27 | 2019-06-27 | Time domain numerical simulation method for rapidly determining micro-discharge threshold of component |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910569890.0A CN110414053B (en) | 2019-06-27 | 2019-06-27 | Time domain numerical simulation method for rapidly determining micro-discharge threshold of component |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110414053A true CN110414053A (en) | 2019-11-05 |
CN110414053B CN110414053B (en) | 2023-06-09 |
Family
ID=68359917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910569890.0A Active CN110414053B (en) | 2019-06-27 | 2019-06-27 | Time domain numerical simulation method for rapidly determining micro-discharge threshold of component |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110414053B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259584A (en) * | 2020-01-14 | 2020-06-09 | 北京航空航天大学 | Time-frequency field combined non-reciprocal medium micro-discharge threshold prediction method |
CN111709179A (en) * | 2020-05-28 | 2020-09-25 | 西安交通大学 | Rapid transition method for micro-discharge development process of microwave component |
CN112836421A (en) * | 2021-02-23 | 2021-05-25 | 西安交通大学 | Multi-scale correlation analysis method for device micro-discharge inhibition |
CN114621778A (en) * | 2020-12-11 | 2022-06-14 | 中国石油化工股份有限公司 | Memory, temperature control method, device and equipment for biomass microwave pyrolysis process |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120012816A1 (en) * | 2009-08-18 | 2012-01-19 | U.S. Government As Represented By The Secretary Of The Army | Photodetectors using resonance and method of making |
CN102930100A (en) * | 2012-10-31 | 2013-02-13 | 西安空间无线电技术研究所 | Method for predicting micro discharge threshold of cavity filter |
WO2017084105A1 (en) * | 2015-11-20 | 2017-05-26 | 田川 | System and method for numerical simulation of plasma discharges |
-
2019
- 2019-06-27 CN CN201910569890.0A patent/CN110414053B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120012816A1 (en) * | 2009-08-18 | 2012-01-19 | U.S. Government As Represented By The Secretary Of The Army | Photodetectors using resonance and method of making |
CN102930100A (en) * | 2012-10-31 | 2013-02-13 | 西安空间无线电技术研究所 | Method for predicting micro discharge threshold of cavity filter |
WO2017084105A1 (en) * | 2015-11-20 | 2017-05-26 | 田川 | System and method for numerical simulation of plasma discharges |
Non-Patent Citations (1)
Title |
---|
王洪广等: "基于频域电磁场的微波器件微放电阈值快速粒子模拟", 《物理学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259584A (en) * | 2020-01-14 | 2020-06-09 | 北京航空航天大学 | Time-frequency field combined non-reciprocal medium micro-discharge threshold prediction method |
CN111709179A (en) * | 2020-05-28 | 2020-09-25 | 西安交通大学 | Rapid transition method for micro-discharge development process of microwave component |
CN114621778A (en) * | 2020-12-11 | 2022-06-14 | 中国石油化工股份有限公司 | Memory, temperature control method, device and equipment for biomass microwave pyrolysis process |
CN114621778B (en) * | 2020-12-11 | 2023-09-01 | 中国石油化工股份有限公司 | Memory, biomass microwave pyrolysis process temperature control method, device and equipment |
CN112836421A (en) * | 2021-02-23 | 2021-05-25 | 西安交通大学 | Multi-scale correlation analysis method for device micro-discharge inhibition |
CN112836421B (en) * | 2021-02-23 | 2022-10-25 | 西安交通大学 | Multi-scale correlation analysis method for micro-discharge inhibition of device |
Also Published As
Publication number | Publication date |
---|---|
CN110414053B (en) | 2023-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110414053A (en) | A kind of time-domain numerical simulation method of quick determining component micro-discharge threshold | |
CN110442919B (en) | Microwave component micro-discharge numerical simulation method based on GPU (graphics processing Unit) architecture | |
CN108920732B (en) | Method for rapidly determining micro-discharge threshold of microwave component loaded by dielectric material | |
CN110276109B (en) | Simulation method for electromagnetic property of plasma sheath of hypersonic aircraft | |
Shoubao et al. | Improvement of current filament method and its application in performance analysis of induction coil gun | |
CN106294894B (en) | Finite element boundary integration method for rapidly analyzing electromagnetic scattering characteristics of non-uniform target | |
Wu et al. | Performance enhanced Crank-Nicolson boundary conditions for EM problems | |
CN112329247A (en) | Simulation method for discharge plasma and self-sputtering of Hall electric thruster | |
Maolin et al. | Particle simulation of grid system for krypton ion thrusters | |
Lu et al. | A fully decoupled iterative method with three-dimensional anisotropic immersed finite elements for Kaufman-type discharge problems | |
CN117350137B (en) | Finite element simulation method for transient characteristics of discharge plasma channel | |
Juanjuan et al. | Investigation of variable aperture on the performance and lifetime of ion thruster | |
CN107944113A (en) | A kind of method for calculating three-dimensional high-speed translation Electromagnetic Scattering of Target field | |
Irishkov et al. | Fully kinetic modeling of low-power hall thrusters | |
CN113517032A (en) | Secondary electron yield calculation method based on first sex principle | |
Berrutti et al. | Multipacting Simulations of SSR2 Cavity at FNAL. | |
Zhang et al. | Preliminary study on the effects of emitter surface roughness on gyrotron electron beam quality | |
CN111563345B (en) | Particle merging method for micro-discharge numerical simulation based on K-D tree data structure | |
Skudnova et al. | On Approach For Resonant Frequency Tuning In Drift Tube Structures On The Designing Stage | |
CN115099074B (en) | Microdischarge threshold prediction method and system based on 1D3V statistical theory modeling | |
Peiyu et al. | Hybrid domain multipactor prediction algorithm and its CUDA parallel implementation | |
Liu et al. | High-order dgtd for solving em scattering from hypersonic aircraft with plasma sheath | |
Wu et al. | Study on the parameters of cavity system-generate electromagnetic pulse (SGEMP) calculation model base on Maxwell-Valsov equation | |
Hao et al. | Study on fluid model of disconnector arc applied in very fast transient overvoltage | |
Ding et al. | A Novel EMC Simulation Analysis Method for Ship-Borne Phased Array Radar System |
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 |