131 results about "Charge density" patented technology

The invention provides a tunnel blasting scheme determining method and device, and relates to the technical field of tunnel engineering data processing.The method compensates stress wave attenuation caused by large-aperture blast hole air gaps through radial explosive density gradient design, so that the energy utilization rate is higher; the axial charge density is dynamically adjusted on the basis of a three-dimensional clamping coefficient, the problem that energy of a hole bottom high constraint area is insufficient is solved, the explosive density and the hole net topology are collaboratively optimized, the contour forming quality is greatly improved, back break is reduced, and therefore the cost is reduced, and the time is shortened. Due to the fact that explosive and hole net parameters are optimized, energy distribution is more reasonable, blasting disturbance is reduced, stability of surrounding rock is enhanced, safety in a hole is greatly improved, and casualties are reduced.

The invention belongs to the field of electromagnetic particle simulation, and relates to a vacuum confluence region electromagnetic particle simulation method, system, equipment and medium, and the method comprises the following steps: S1, constructing a particle simulation domain of an electric field and a magnetic field of a vacuum confluence region, injecting particles into the particle simulation domain, and carrying out the regional decomposition of the particle simulation domain to obtain a plurality of sub-domains; s2, acquiring charge density and current density of grid points in each sub-domain; discretizing a Maxwell equation set into a form of a Cx = D equation, substituting the charge density and the current density of the grid points in each sub-domain into the Cx = D equation, and solving the Cx = D equation by adopting a GMRES algorithm to obtain potential and vector magnetic potential distribution in each sub-domain; s3, the electric field intensity and the magnetic induction intensity of the grid points in each sub-domain are solved, the electric field and the magnetic field borne by the particles are obtained, and therefore the positions of the particles in the next time step in each sub-domain are determined; and S4, updating the electric field and the magnetic field borne by the particles in each sub-domain in the next time step.

The invention relates to a direct current GIL / GIS surface charge and surface discharge on-line monitoring method, which is used for a direct current GIL / GIS, and comprises the following steps: S1, carrying out the corresponding division of the surface area of an insulator based on a number; s2, establishing a GIL / GIS model, and performing corresponding numbering and division; s3, the simulated induced voltage is calculated, and a contribution coefficient matrix is calculated based on the simulated induced voltage; s4, on the basis of the induction current, actual induction voltage is obtained through calculation in combination with the measurement resistance, and actual surface charge density is obtained through calculation; and S5, judging the position area of the creeping discharge point according to the amplitude of the pulse current on each segmented electrode. Compared with the prior art, the direct-current GIL / GIS basin-type insulator surface charge online measurement device has the advantages of realizing online measurement of the direct-current GIL / GIS basin-type insulator surface charge and the

The invention discloses a self-adaptive Monte Carlo macro particle simulation method and system based on parallel computing, and relates to the technical field of particle simulation, and the method comprises the steps: initializing parameters, if the current step number is not the maximum step number, storing node data, balancing particle loading, and calculating charge density and battery density; using macro particles to represent a group of real particles, updating the positions and speeds of all the particles in parallel, and calculating new charge distribution; setting adaptive grid configuration, and dynamically adjusting grid density according to electric field intensity and local density of charge distribution; and based on the updated particle position and speed, evaluating and processing a collision event, storing collision data, re-calculating surface charge density data, adding 1 to the current step number, and circularly executing to finish simulation when the maximum step number is reached. According to the method, the accuracy of a simulation result is improved through the improved collision model, electrical discharge behaviors under different conditions are better predicted, and more reliable data support is provided for design and safety evaluation of high-voltage equipment.

The invention discloses a VASP-based metal surface adsorption energy differential charge analysis method and system, and relates to the technical field of metal surface adsorption energy differential charge analysis, and the method comprises the steps: carrying out the modeling and standardization processing of a metal surface-adsorbate system, and obtaining a standard data set; training the standard data set by adopting a deep neural network model to obtain a prediction model, inputting the standard data set into the prediction model, and outputting a prediction result; dynamically adjusting VASP calculation parameters by adopting an environmental effect correction method, and generating a VASP input file based on a prediction result and an actual operation condition; performing high-precision energy and electron density calculation on the input file by adopting a first principle calculation tool VASP to obtain electron density distribution; performing differential operation on the electron density distribution to obtain a differential charge density map; and performing three-dimensional rendering on the differential charge density map by adopting a visual analysis tool to obtain a visual analysis result.

The invention discloses a high-performance layered nitride thermoelectric material screening method and system. The method comprises the following steps: acquiring a series of crystal structures of initial layered nitride; carrying out magnetic test on the material by utilizing first principle calculation, and selecting a non-magnetic material; carrying out structure optimization on the non-magnetic material, carrying out static calculation and electron energy band structure calculation after calculation convergence, and selecting a material with a band gap range of 0.1-3 eV and high energy band structure degeneracy; calculating partial charge density of the material, and judging whether two-dimensional electron gas exists or not, namely, alternate occurrence of a conductive layer and an insulating layer; stability judgment: carrying out phonon spectrum and molecular dynamics simulation calculation on the material, and selecting the material with phonon spectrum without imaginary frequency and stable dynamics; the heat transport performance of the material is calculated, the lattice heat conductivity and Green Ehson parameters are obtained, and the material with high non-harmony is found; and finally, carrying out electric transport calculation to obtain the ZT value of the material.

The invention belongs to the technical field of electrical performance research of insulating materials, and relates to a creeping discharge surface charge high-precision measurement method, system and equipment and a medium. According to the method, basic data are provided for subsequent analysis by obtaining the creeping discharge path image, the potential scanning result and the dust diagram distribution characteristic of the insulating material; performing inversion calculation on a potential scanning result according to the shape of the insulating material in combination with a corresponding inversion algorithm to obtain a surface charge distribution characteristic; according to the surface discharge path image, the surface charge distribution characteristic and the dust diagram distribution characteristic, the characteristics of the stream-pilot channel in the surface discharge process are obtained, the charge distribution and space distribution characteristics of the high-precision stream-pilot discharge channel of the stream can be obtained, and the charge density in the stream-pilot channel can be quantitatively distinguished. In addition, the microcosmic physical process in the discharge process can be obtained more finely, and the microcosmic physical process comprises reverse discharge, charge neutralization in a plasma channel and an electron collapse generation process around a streamer channel.

The low-pressure ionization chamber comprises a closed chamber, a low-pressure environment lower than the standard atmospheric pressure is formed in the closed chamber, and an electrode plate assembly is arranged in the closed chamber to collect electron-ion pairs in the closed chamber in the low-pressure environment; by arranging the negative-pressure closed chamber, the air density is reduced under the condition of low air pressure, so that the number of gas molecules in unit volume is reduced, the total number of charges generated by radiation ionization of a gas medium is reduced, and then the charge density and the ion recombination probability are reduced, so that the saturation characteristic of the ionization chamber is influenced.

A method and corresponding devices to reduce the presence of unwanted ions within the pseudo-potential well of the RF fields for inductively coupled plasma mass spectrometry is disclosed. This method reduces space-charge effects and leaves more room available for the ions of interest and, hence, leads to increased transmission of the desired ions into the mass analyzer. In the present invention, the space charge is reduced by adding a low-band resolving DC potential, or a fast sweeping resolving DC potential with a notch, or a fast constant resolving DC potential with a notch, or by adding auxiliary quadrupolar or dipolar RF potential to the quadrupole ion guide, such that all the unwanted ions will be outside the stable region and therefore ejected. This reduces the charge density in the ion beam and ameliorate space charge effects.

The invention belongs to the field of nano energy, and provides a surface charge density prediction method of a friction nano generator based on a modified chitosan material, and the method comprises the following steps: S1, preparing modified chitosan; s2, carrying out triboelectric polarity test and quantitative characterization on the modified chitosan; S3, obtaining surface potential energy and saturated charge density of the modified chitosan and the pairing material; and S4, constructing a surface charge density prediction model of the friction nano-generator based on the modified chitosan material. According to the modification degree of the modified chitosan material, the surface charge density of the friction nano-generator based on the modified chitosan material can be directly obtained, and a solid theoretical basis is provided for subsequent deep exploration of application of chitosan in the triboelectricity field.

The invention relates to the technical field of partial discharge processing, and discloses a partial discharge processing method and system based on electric field analysis, and the method comprises the steps: synchronously collecting the multi-point electric field intensity and charge density of a device surface, recognizing a charge abnormal region, and calculating an electric field non-uniformity index; inputting the information into a finite element analysis model, simulating an electric field distortion process, and positioning potential discharge induction points; predicting a charge accumulation trend of each potential point in combination with historical operation data, and constructing an optimization function taking minimization of global electric field distortion and charge accumulation risk as a target according to the charge accumulation trend; and with the function as a convergence target, an optimal electric field compensation parameter is automatically solved through iterative calculation. According to the partial discharge processing method and system based on electric field analysis, predictive recognition and counteracting are carried out on discharge core inducements, and active and accurate suppression of partial discharge of equipment is achieved.

The invention relates to the technical field of monitoring modeling, in particular to an electrostatic precipitation smoke emission dynamic modeling method based on multi-parameter coupling analysis, which comprises the following steps: calculating electric field rate and acceleration through sliding window difference, judging electric field abrupt change trigger parameter adjustment, normalizing particle size and charge density, and mapping disturbance intensity index. And calculating response delay and a gradient field to generate a disturbance coefficient, transferring a model weight when the disturbance exceeds a threshold value, performing fitting reconstruction, and performing regression weight reverse normalization to output a steady-state parameter group. According to the method, electric field dynamic characteristics are captured through equal-interval sampling in combination with a sliding window, an acceleration super-threshold model is constructed through second-order difference, a dynamic adjustment mechanism is triggered, particle size and charge density are normalized, disturbance intensity is mapped through linear interpolation, speed and density abrupt change are quantized through gradient weighting, and weight transfer and a least square method are introduced to reconstruct parameters. The transient analysis capability and robustness of the model are improved, and high-precision synchronous mapping of the emission trend and the equipment state is realized.

The low-pressure ionization chamber comprises a closed chamber, a low-pressure environment lower than the standard atmospheric pressure is formed in the closed chamber, and an electrode plate assembly is arranged in the closed chamber to collect electron-ion pairs in the closed chamber in the low-pressure environment; by arranging the negative-pressure closed chamber, the air density is reduced under the condition of low air pressure, so that the number of gas molecules in unit volume is reduced, the total number of charges generated by radiation ionization of a gas medium is reduced, and then the charge density and the ion recombination probability are reduced, so that the saturation characteristic of the ionization chamber is influenced.

The invention discloses a side-to-boundary blasting method for sedimentary metamorphic surface iron ore, and relates to the technical field of iron ore blasting construction, which comprises the following steps of: drilling: selecting a proper drilling position according to a rock mass exploration result; charging: uniformly distributing ammonium nitrate fuel oil explosives in the drill hole by adopting a segmented charging mode, and installing a primer detonator; blocking: after charging is completed, cement paste and gravel are mixed to form a blocking material to tightly block the main blasting hole; detonating: transmitting a detonating signal, detonating the detonator, and enabling the explosive to explode in the main blasting hole; according to the method, the charge density of the presplitting line is adjusted to 0.5 kg / m from 0.7 kg / m, the presplitting line is segmented into 10 m sections, the single-shot charge is controlled to be 60 kg or below, the use amount of explosives is reduced, the propagation range and strength of shock waves and stress waves generated during blasting can be reduced, direct impact and damage to the side slope can be reduced, and the mode is beneficial to saving cost and generating economic benefits.

The invention relates to the technical field of nano materials, and provides a preparation method and application of a friction nano generator with high surface charge density. The invention aims to overcome the defects of low surface charge density, environmental influence, easy aging and the like in high-frequency mechanical energy collection of the existing friction nano-generator, and the application of the friction nano-generator in environmental mechanical energy collection is improved. The main scheme is as follows: 4, the preparation method comprises the following steps: dissolving 4, 4 '-diaminodiphenyl ether in N, N'-dimethylacetamide, adding diphenyl ether dianhydride into the solution to synthesize a polyamide acid solution, carrying out electrospinning and high-temperature treatment on the solution to obtain an imidized nanofiber membrane, treating the surface of the nanofiber membrane through a plasma sputtering instrument, selecting a fluorinated ethylene-propylene copolymer as an electric negative electrode material, and carrying out electrospinning and high-temperature treatment on the surface of the nanofiber membrane. And pasting an electrode and packaging to construct the friction nano generator. The invention is used in the fields of friction power generation, mechanical sensing and the like, and has the advantages of simple process, high preparation speed, high yield, low cost, good robustness, no environmental influence and high output performance.

The utility model belongs to the field of ionization electrodes for separating dispersed particles from gas by using an electrostatic effect, and particularly relates to a discharge electrode for an electrostatic trap and the electrostatic trap. The discharge electrode for the electrostatic trap comprises an electrode body, one end of the electrode body is provided with a fixing structure used for being fixed to an electrode fixing pipe of the electrostatic trap, the electrode body comprises a cylindrical discharge section used for discharging, the periphery of the cylindrical discharge section is provided with external threads, and threads of the external threads form a continuous discharge tip. The cylindrical structure can increase the wall thickness of the discharge electrode, so that the discharge electrode is not easy to wear or deform, the service life is prolonged, and meanwhile, the thread of the external thread can form a discharge tip, so that discharge occurs at the spiral end, the charge density is improved, and the electric field distribution is more uniform.

The application discloses a GIS insulator partial discharge charge on-line detection method and system, and relates to the technical field of charge detection.The method comprises the following steps: according to a GIS simulation model, simulating the electric field distribution of charge generation in the partial discharge process, determining the performance index of the required electric field sensing probe, and preparing the electric field sensing probe; according to an epoxy resin model, obtaining the charge density, the charge distribution range and the electric field signal collected by the electric field sensing probe under different voltage pulses, and constructing a test data set; according to the GIS simulation model, obtaining the electric field signal under different charge densities and different charge distribution radii, and constructing a training data set; training a neural network model according to the training data set, and correcting according to the test data set; and collecting the partial discharge electric field signal by using the electric field sensing probe array arranged on the outer wall of the GIS insulator, and obtaining the charge density and the charge distribution range on the surface of the insulator by using the neural network. Non-invasive and non-contact on-line measurement is realized.

The application provides a method for calculating interface polarization charge and current distribution based on reflection coefficient, and relates to the technical field of computational electromagnetics. The method comprises the following steps: establishing a reflection coefficient-dielectric parameter mapping, directly correlating the reflection characteristics of electromagnetic waves with the intrinsic parameters of the medium through a reflection coefficient formula, and laying a physical foundation for polarization analysis; interface electromagnetic field reconstruction, using boundary conditions and constitutive relations to determine the electric field distribution on both sides of the interface and the corresponding polarization intensity. The distribution rule of the electromagnetic field in the medium is determined, the macroscopic electric field is converted into the spatial description of the polarization intensity, and an intermediate variable is provided for charge and current calculation; polarization dynamic response quantification, through the spatial and temporal derivatives of the polarization intensity, the interface polarization charge density and the polarization current density are derived. The theoretical model is converted into actual physical quantities, and the dynamic response mechanism of charges and currents when electromagnetic waves interact with the medium is revealed.

A method and an apparatus for measuring an electrical property of a sample material are provided. In the method, an alternating voltage is applied to a piezoelectric ceramic so that a probe fixed on the piezoelectric ceramic vibrates above a surface of the sample material, a target contact potential difference between the probe and the sample material is measured, and a target amplitude of an induced alternating current generated due to the vibration of the probe above the surface of the sample material, a target work function of the surface of the sample material is determined according to the target contact potential difference and a determined work function of the probe, and a target charge density of the surface of the sample material is determined according to the target amplitude, the target contact potential difference, and a pre-stored target amplitude determination function for the induced alternating current.

The invention relates to a method for measuring the surface charge of a double-sided asymmetric basin-type insulator. The method comprises the following steps: S1, designing a coaxial electrode system matched with the basin-type insulator; s2, connecting the basin-type insulator and the coaxial electrode system with the experiment cavity; s3, charging the basin-type insulator; s4, acquiring potential data of concave and convex surfaces of the basin-type insulator; and S5, measuring the surface charge of the basin-type insulator. Based on an active probe theory, a three-axis displacement platform is built to control three-dimensional operation of a Kelvin probe, a rotating platform is designed to be used for constant-speed coaxial rotation of the basin-type insulator, and a three-dimensional motion track of the Kelvin probe suitable for surface scanning of the basin-type insulator of a double-sided asymmetric structure is designed. The surface potential data is converted into charge density data through inversion calculation, so that the surface charge of the basin-type insulator with the asymmetric surface can be quickly and accurately measured.

The invention belongs to the field of electrochemical behaviors of seawater battery anodes, and particularly relates to a method for predicting constant-current discharge performance of a magnesium anode by using a first principle. Comprising the following steps: (1) establishing a matrix model according to atomic configurations and arrangement modes of different anodes; cutting the sample surface, and arranging a vacuum layer and a solvent layer; (2) carrying out geometric optimization on the structure formed in the step (1); (3) in combination with a joint density functional theory, putting the geometrically optimized structure into a solvent environment for solvation, and calculating the total energy, interface electron density, charge density and electrode potential of the structure; and (4) predicting a potential-time relationship of anode constant current discharge within the set acting time by using a calculation result in the step (3) and combining a discharge dynamic law of a chronopotentiometry. According to the method, the magnesium anode surface electric double-layer model in seawater can be obtained, and the electrochemical properties of the magnesium anode surface electric double-layer model can be quickly and effectively predicted based on the joint density functional theory.

The present invention relates to the fields of chemistry and materials science, and discloses a visualization system for large-scale quantum chemistry, which calculates the energy information of each structure in the calculation system, and finds the minimum-energy structure in the system as the optimized structure by means of gradient descent, and displays the structures before and after optimization in the form of an atomic coordinate matrix for easy visualization; at the same time, it can also obtain information such as charge density, differential charge density tensor, difference tensor, etc. according to the energy information of the structure, and draw visual images by using these information.

The invention relates to a superlattice ferroelectric enhancement structure in a hafnium oxide-based film as well as a preparation method and application thereof, and aims to realize remarkable improvement of ferroelectric properties of a hafnium oxide-doped film by stabilizing a ferroelectric orthogonal phase. The superlattice structure of the doped hafnium oxide thin film comprises bottom electrode layers and doped hafnium oxide thin film layers, wherein the bottom electrode layers are sequentially arranged in the direction perpendicular to the substrate, and the doped hafnium oxide thin film layers are alternately stacked. The superlattice interface engineering prepared through the atomic layer deposition or magnetron sputtering technology can effectively inhibit generation of non-ferroelectric phases (such as a monoclinic phase and a tetragonal phase), so that ferroelectric remanent polarization is remarkably enhanced. According to the doped hafnium oxide thin film device based on the superlattice structure, the voltage difference in the polarization inversion period can be enhanced, the state distinguishing in the reading operation period can be improved, the reading error rate can be reduced, the charge density stored in the unit area can be increased, the device size can be reduced, and the device performance can be improved. Therefore, the integration density and the storage density of the ferroelectric-based memory are improved.

The invention relates to a physical simulation model of a ferroelectric neuron circuit and a design method thereof, and the method comprises the steps: building the physical simulation model of the ferroelectric neuron circuit according to the corresponding relation between the polarization intensity of a ferroelectric layer in a ferroelectric field effect transistor and the voltage of the ferroelectric layer, the surface charge density of the ferroelectric layer, and the source voltage of the ferroelectric field effect transistor, and predicting the working state of the ferroelectric neuron circuit according to the physical simulation model. Meanwhile, a corresponding experimental design method is further developed according to the provided model. The physical simulation model combining the ferroelectric polarization flipping and the dynamic characteristics of the ferroelectric neuron circuit has very high physical reliability; the design method provided by the invention can guide the specific design of the experiment, significantly reduces the trial and error cost of the experiment manufacturing of the ferroelectric neuron circuit, and provides a physical basis for the further development of the ferroelectric neuron circuit in electronic design automation software.

The invention discloses an electrode structure, a treatment head and a radio frequency therapeutic apparatus, the electrode structure comprises an electrode body, a first electrode and a second electrode, the first electrode and the second electrode are distributed on the surface of the electrode body, the total area of the first electrode is consistent with the total area of the second electrode, the first electrode and the second electrode are opposite in electrical property, and the first electrode and the second electrode are opposite in electrical property. The first electrode and the second electrode are adjacent to each other and are arranged at an interval; wherein the width of the electrode located on the outer periphery is smaller than that of the electrode located on the inner periphery. According to the technical scheme, the width of the peripheral electrode is set to be smaller than that of the inner peripheral electrode, so that the cross sectional area of the peripheral electrode is increased, the difference value of the charge density of the two ends of the peripheral electrode is reduced, charges are evenly distributed on the electrodes, and the phenomenon of local overheating is avoided.

The invention discloses a method for determining blasting parameters of a pre-splitting hole in pre-splitting blasting. The method comprises the steps that the blast hole diameter, the charge non-coupling coefficient, the charge density, the blast hole density coefficient and the blast hole length of the pre-splitting hole are determined; according to the blast hole diameter of the pre-splitting hole and the charge non-coupling coefficient, the charge diameter of the pre-splitting hole is obtained; according to the charging diameter and the charging density of the pre-split hole, the linear charging density of the pre-split hole is obtained; according to the linear charge density of the pre-splitting hole and the length of the blast hole, the single-hole charge amount of the pre-splitting hole is obtained; according to the single-hole explosive charge / linear charge density of the pre-split holes, the pitch of the pre-split holes is obtained based on the blast hole density coefficient and the pre-split blasting unit consumption; and the filling length and the charging structure of the pre-split holes are obtained according to the pitch of the pre-split holes. By means of the calculation method, the design threshold of presplitting hole parameters in presplitting blasting is lowered, and the blasting success rate of primary design (or one-time engineering) is remarkably increased.

The invention discloses a spatial gravitational wave detection inertial sensor plaque field effect simulation and related physical quantity calculation method, a sensitive structure of an inertial sensor comprises an electrode cage and a test mass, and the plaque field effect comprises that bumps appear on the outer surface of the test mass and / or the inner surface of the electrode cage; chemical pollution is formed to change the local potential and influence the electric field; the method comprises the following steps: respectively dividing the outer surface of a test mass without plaques and the inner surface of an electrode cage into a plurality of triangular units; moving partial dissection nodes which are not on the edges to enable the dissection nodes to leave the original plane to simulate bumps caused by plaques; based on a set condition that the potential at the centroid of the triangular unit is the same as the potential of the area where the triangular unit is located, simulating the influence of plaques on the potential; based on setting of uniform distribution of charges on the surface of the same triangular unit, the charge density of the surface of each triangular unit is solved; and according to the solved charge density, calculating the electric field force borne by the test mass and the rigidity of the electric field force.

Disclosed in the present invention are a particle accommodating device, and a particle contactless processing apparatus and a light sensing structure thereof. The particle accommodating device is used for accommodating a liquid specimen that contains a plurality of particles. The accommodating device comprises a light sensing structure and a mating structure spaced apart from the light sensing structure. The light sensing structure comprises a substrate, an electrode layer and a photoelectric layer which are formed on the substrate, and an insulator formed on the photoelectric layer. The insulator comprises a charge-type track protruding from the photoelectric layer. The photoelectric layer is operable to generate a plurality of charges that are concentrated on the charge-type track, so that the charge density of the charge-type track is greater than that of the outer surface of the photoelectric layer, and thus the charge-type track can position at least one of the particles thereon by electrostatic adsorption. Thus, the particles are moved and positioned by means of the electrostatic adsorption, thereby achieving precise positioning of the particles.

The application discloses a composite defect signal quantification detection method, relates to the technical field of industrial nondestructive testing, and is characterized in that: a defect space leakage magnetic field is established, a three-dimensional rectangular coordinate system is established based on a magnetic charge model, magnetic charge density corresponding to a magnetic charge unit segmented on a defect end face is established, the relationship between a ferromagnetic material and stress is determined, magnetic induction intensity and unit area corresponding to the magnetic charge unit segmented on the defect end face are established, total magnetic charge density, magnetic induction intensity and area of the defect end face are established, and then axial and radial components of a defect magnetic signal are solved. The application fully considers the influence of stress on defects of the ferromagnetic material, simulates and calculates signal characteristics of the magnetic signal of the ferromagnetic material under the action of different defect stresses, and provides a more comprehensive, efficient and stable detection method for defect evaluation of the ferromagnetic component by the leakage magnetic detection technology.

The invention discloses an implicit particle simulation analysis method for atmospheric pressure micro-scale arc discharge kinetic characteristics, and the method specifically comprises the steps: 1, carrying out the uniform grid division of a calculation region, and placing a plurality of charged simulation particles with initial positions and speeds in each grid; step 2, calculating charge density and current density on each grid node; 3, electric potential and electric field distribution on each grid node is calculated; 4, calculating the stress of the charged simulation particles and updating the coordinates and speed of the charged simulation particles; 5, calculating a new speed after collision based on a collision theory; and step 6, taking the updated coordinates and speed as the coordinates and speed in the step 2, and repeatedly executing the steps 2-6 in a time step until a set number of iterations is reached. According to the method, the problem that simulation time consumption is remarkably increased due to the fact that the space-time step length must be small enough to distinguish the Debye length and the oscillation period of electrons in the existing simulation technology is solved.