Numerical simulation method, device, equipment and medium for electromagnetic field induced by charged particles

Abstract

The present application discloses a numerical simulation method, device, equipment and medium of the induced electromagnetic field of charged particles, which relates to the field of numerical simulation technology, including: obtaining the target fluid velocity vector of the position of the charged particles in the target turbulent field currently output by the preset fluid solver; solving the current motion information of the particles at the current time step based on the target fluid velocity vector and the historical motion information and historical electromagnetic field information of the particles updated at the previous time step; wherein the motion information includes the velocity vector and the displacement vector, and the electromagnetic field information includes the electric field vector and the magnetic field vector; using the current motion information of the particles to update the electric field vector and the current vector of the current time step, and based on the current vector to update the magnetic field vector of the current time step. The present application can more accurately perform numerical simulation of the induced electromagnetic field of charged particles.

Description

Technical Field

[0001] The present invention relates to the field of numerical simulation technology, and in particular to a numerical simulation method, device, equipment and medium for induced electromagnetic fields of charged particles. Background Art

[0002] The generation of supercooled clouds, the combustion of fuel droplets in aircraft engines, or particles in sandstorms all have the phenomenon of charged particles inducing electromagnetic fields. The source of particle charge is the collision and friction between particle groups. Charged particles move in turbulent fields, inducing electromagnetic fields, which in turn deflect the movement trajectory of charged particles. The two are coupled in a closed loop and influence each other.

[0003] However, the traditional particle two-phase flow solver does not consider the particle charging effect, and usually only simulates the electric field force without considering the magnetic field force. Therefore, it cannot accurately reflect the actual situation and does not conform to engineering practice. For aircraft engine design, it is crucial to accurately simulate various physical phenomena in the process of fuel droplet combustion, which can provide a more reliable basis for the optimization design of the engine; in the field of supercooled fog research, numerical simulation can provide a deep understanding of the mechanism of charged particles in the process of fog formation, and provide more accurate theoretical support for weather forecasting, aviation safety and other aspects.

[0004] In summary, how to more accurately perform numerical simulation of the induced electromagnetic field of charged particles is a problem that needs to be solved. Summary of the invention

[0005] In view of this, the purpose of the present invention is to provide a method, device, equipment and medium for numerical simulation of the induced electromagnetic field of charged particles, which can more accurately perform numerical simulation of the induced electromagnetic field of charged particles. The specific scheme is as follows:

[0006] In a first aspect, the present application discloses a numerical simulation method of an induced electromagnetic field of charged particles, comprising:

[0007] Obtain the target fluid velocity vector at the location of the charged particles in the target turbulent field currently output by the preset fluid solver;

[0008] Based on the target fluid velocity vector and according to the particle historical motion information and historical electromagnetic field information updated in the previous time step, the current motion information of the particle in the current time step is solved; wherein the motion information includes a velocity vector and a displacement vector, and the electromagnetic field information includes an electric field vector and a magnetic field vector;

[0009] The electric field vector and the current vector of the current time step are updated using the current movement information of the particle, and the magnetic field vector of the current time step is updated based on the current vector.

[0010] Optionally, the numerical simulation method of the charged particle induced electromagnetic field further includes:

[0011] The fluid velocity vector at the center position of the charged particle is determined as the target fluid velocity vector at the position where the charged particle is located.

[0012] Optionally, solving the current motion information of the particles at the current time step based on the target fluid velocity vector and according to the particle historical motion information and historical electromagnetic field information updated at the previous time step includes:

[0013] Determine a charged particle dynamics equation constructed based on an electric field vector, a magnetic field vector, and a velocity vector and a displacement vector between adjacent time steps; wherein known parameters in the charged particle dynamics equation include particle historical motion information and historical electromagnetic field information updated in the previous time step;

[0014] The target fluid velocity vector and the characteristic time of particle motion are substituted into the charged particle dynamics equation to solve the current velocity vector and the current displacement vector of the charged particle in the current time step.

[0015] Optionally, the process of determining the characteristic time of the particle motion includes:

[0016] The characteristic time of motion of the particles is calculated based on the density, diameter and fluid dynamic viscosity of the charged particles.

[0017] Optionally, updating the electric field vector of the current time step using the current motion information of the particle includes:

[0018] Determining a first charge amount of the charged particle at any target position based on a current displacement vector of the charged particle, and determining a second charge amount on different grids after the target turbulent field is gridded;

[0019] Determine the displacement vector diameter between the position where the first charge quantity is located and the position where the second charge quantity is located, so as to update the electric field vector of the current time step by using the first charge quantity, the second charge quantity and the displacement vector diameter.

[0020] Optionally, updating the magnetic field vector of the current time step based on the current vector includes:

[0021] Determine the preset curl operator;

[0022] Based on the curl operator, a curl operation is performed on the current vector and the displacement vector to update the magnetic field vector of the current time step.

[0023] Optionally, updating the current vector of the current time step using the current motion information of the particle includes:

[0024] The velocity change rate of the charged particle in different directions is determined based on the current velocity vector of the charged particle, and the current vector of the current time step is updated using the charge amount of the charged particle and the velocity change rate.

[0025] In a second aspect, the present application discloses a numerical simulation device for an induced electromagnetic field of charged particles, comprising:

[0026] A fluid velocity vector acquisition module is used to acquire a target fluid velocity vector at a location where charged particles are located in a target turbulent field currently output by a preset fluid solver;

[0027] A motion information updating module, used to solve the current motion information of the particle at the current time step based on the target fluid velocity vector and the particle historical motion information and historical electromagnetic field information updated at the previous time step; wherein the motion information includes a velocity vector and a displacement vector, and the electromagnetic field information includes an electric field vector and a magnetic field vector;

[0028] The electromagnetic field information updating module is used to update the electric field vector and the current vector of the current time step by using the current movement information of the particle, and calculate the magnetic field vector of the current time step based on the current vector.

[0029] In a third aspect, the present application discloses an electronic device, comprising:

[0030] Memory, used to store computer programs;

[0031] The processor is used to execute the computer program to implement the steps of the numerical simulation method of the charged particle induced electromagnetic field disclosed above.

[0032] In a fourth aspect, the present application discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, the steps of the numerical simulation method of charged particle induced electromagnetic field disclosed above are implemented.

[0033] It can be seen that the present application obtains the target fluid velocity vector of the position of the charged particles in the target turbulent field currently output by the preset fluid solver; based on the target fluid velocity vector and according to the historical motion information and historical electromagnetic field information of the particles updated in the previous time step, the current motion information of the particles in the current time step is solved; wherein the motion information includes a velocity vector and a displacement vector, and the electromagnetic field information includes an electric field vector and a magnetic field vector; the electric field vector and current vector of the current time step are updated using the current motion information of the particles, and the magnetic field vector of the current time step is updated based on the current vector.

[0034] It can be seen that after obtaining the target fluid velocity vector of the charged particle position in the target turbulent field currently output by the preset fluid solver, the present application can solve the current motion information of the particle in the current time step based on the target fluid velocity vector and the particle historical motion information and historical electromagnetic field information updated in the previous time step; wherein the motion information specifically includes the velocity vector and the displacement vector, and the electromagnetic field information specifically includes the electric field vector and the magnetic field vector. Furthermore, since the position of the particle is updated, and the position change between different particles will affect the electric field force, it is necessary to update the electric field force information it is subjected to. In addition, the particle velocity changes, causing its current to change, which will also affect the magnetic field force. Therefore, the present application needs to use the particle current motion information to update the electric field vector and current vector of the current time step, and then update the magnetic field vector of the current time step based on the current vector. Through the above scheme, the induced electromagnetic field of charged particles can be numerically simulated more accurately. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0036] Figure 1 A flow chart of a numerical simulation method for an induced electromagnetic field of charged particles disclosed in the present application;

[0037] Figure 2 A flowchart of a specific numerical simulation method for the induced electromagnetic field of charged particles disclosed in this application;

[0038] Figure 3 This is a schematic structural diagram of a numerical simulation device for charged particle induced electromagnetic field disclosed in this application;

[0039] Figure 4 This is a structural diagram of an electronic device disclosed in this application. DETAILED DESCRIPTION

[0040] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0041] Traditional particle two-phase flow solvers do not consider the particle charging effect, and usually only simulate the electric field force without considering the magnetic field force. Therefore, they cannot accurately reflect the actual situation and are not in line with engineering practice. For aircraft engine design, it is crucial to accurately simulate various physical phenomena in the process of fuel droplet combustion, which can provide a more reliable basis for the optimization design of the engine; in the field of supercooled fog research, numerical simulation can provide a deep understanding of the mechanism of charged particles in the process of fog formation, and provide more accurate theoretical support for weather forecasting, aviation safety and other aspects.

[0042] To this end, the embodiments of the present application disclose a method, device, equipment and medium for numerical simulation of the induced electromagnetic field of charged particles, which can more accurately perform numerical simulation of the induced electromagnetic field of charged particles.

[0043] See also Figure 1 and Figure 2 As shown, the embodiment of the present application discloses a numerical simulation method of an induced electromagnetic field of charged particles, the method comprising:

[0044] Step S11: Obtain a target fluid velocity vector at a location where charged particles are located in a target turbulent field currently output by a preset fluid solver.

[0045] In this embodiment, the target turbulent field can be a turbulent field formed by particles in supercooled clouds, aircraft engines or sandstorms during turbulent motion. Charged particles moving in a turbulent field will induce an electromagnetic field, and the induced electromagnetic field will in turn deflect the motion trajectory of the charged particles. The two are closed-loop coupled and influence each other. Therefore, in order to numerically simulate the electromagnetic field in these scenarios, the application first needs to obtain the target fluid velocity vector at the location of the charged particles in the target turbulent field currently output by the preset fluid solver.

[0046] Specifically, the above method also includes: determining the fluid velocity vector at the center position of the charged particle as the target fluid velocity vector at the position of the charged particle. It is understandable that during the simulation process, the fluid velocity at the position of the particle is considered to be uniform and can be approximately represented by the fluid velocity at the center position of the particle. This processing method is common in numerical simulation, especially when using the Lagrangian method to track particle motion, each particle is regarded as an independent entity, and its motion is affected by the surrounding fluid velocity. In order to simplify the calculation, it is often assumed that the fluid velocity around the particle is uniform, so that the overall motion of the particle can be determined by calculating the fluid velocity at the center position of the particle.

[0047] Step S12: Based on the target fluid velocity vector and according to the particle historical motion information and historical electromagnetic field information updated in the previous time step, solve the particle current motion information in the current time step; wherein the motion information includes the velocity vector and the displacement vector, and the electromagnetic field information includes the electric field vector and the magnetic field vector.

[0048] In this embodiment, after obtaining the target fluid velocity vector at the location of the charged particles in the target turbulent field, the current motion information of the particles in the current time step can be solved based on the target fluid velocity vector and the historical motion information and historical electromagnetic field information of the particles updated in the previous time step; wherein the motion information specifically includes the velocity vector and the displacement vector, and the electromagnetic field information specifically includes the electric field vector and the magnetic field vector.

[0049] In a specific embodiment, the above-mentioned method of solving the current motion information of the particles in the current time step based on the target fluid velocity vector and according to the particle historical motion information and historical electromagnetic field information updated in the previous time step includes: determining a charged particle kinetic equation constructed based on the electric field vector, the magnetic field vector, and the velocity vector and the displacement vector between adjacent time steps; wherein the known parameters in the charged particle kinetic equation include the particle historical motion information and historical electromagnetic field information updated in the previous time step; substituting the target fluid velocity vector and the characteristic time of the particle motion into the charged particle kinetic equation to solve the current velocity vector and the current displacement vector of the charged particle in the current time step.

[0050] The dynamic equation for charged particles is as follows:

[0051] ;

[0052] In the formula, all vectors are expressed in the form of Einstein index, such as represents the velocity vector of a particle at any position, represents the particle displacement vector, Indicates the target fluid velocity vector. The values ​​of i, j, and k in the subscripts are 1, 2, and 3, indicating three different directions. The superscript n indicates time n, and n+1 indicates time n+1. represents the discrete time step size; represents the force exerted by the electric field on the particle, represents the force exerted by the magnetic field on the particle, represents the charge value of the charged particle, represents the electric field vector, represents the magnetic field vector, represents the curl operator, Represents the characteristic time of particle motion.

[0053] The process of determining the characteristic time of particle motion includes: calculating the characteristic time of particle motion based on the density, diameter and fluid dynamic viscosity of the charged particles. The specific calculation formula is as follows:

[0054] ;

[0055] In the formula, represents the density of charged particles, represents the diameter of the charged particle, Represents the dynamic viscosity of the fluid.

[0056] From the above content, we can know that the target fluid velocity vector is the fluid velocity vector at the center of the charged particle, and the expression is as follows:

[0057] ;

[0058] This formula represents the target fluid velocity vector at the location of the charged particle in the nth time step. The center of the charged particle The fluid velocity vector .

[0059] Therefore, based on the above solution of the charged particle dynamics equation, the current velocity vector of the charged particle at any point in space can be updated: and the current displacement vector .

[0060] It should also be pointed out that the dynamic equation of charged particles follows the Lagrangian method. When analyzing a single particle, the force it is subjected to can be expressed as the interphase resistance in the first half and the Lorentz force in the second half. In Cartesian coordinates, vector cross products can be expressed using the Levi-Civita notation as follows:

[0061] ;

[0062] is the Levi-Civita symbol, which is defined as follows:

[0063] (1) When When it is an even permutation of 1, 2, 3, that is, an even arrangement, ;

[0064] (2) When When it is an odd permutation of 1, 2, 3, that is, an odd arrangement, ;

[0065] (3) In other cases, ;

[0066] Therefore, the above formula can be written as:

[0067] ;

[0068] In summary, the velocities of the particles in three different directions can be obtained, and similarly, the displacements of the particles in three different directions can also be obtained.

[0069] It should be noted that in a Cartesian coordinate system, the origin is usually an arbitrarily selected point, usually marked as O or simply the origin. In most applications, the origin is set at a convenient place for measurement and calculation.

[0070] Step S13: using the current movement information of the particle to update the electric field vector and the current vector of the current time step, and updating the magnetic field vector of the current time step based on the current vector.

[0071] In this embodiment, since the position of the particles is updated, and the position change between different particles will affect the electric field force, it is necessary to update the electric field force information. In addition, the speed of the particles changes, causing their current to change, which will also affect the magnetic field force. Therefore, this application needs to use the current movement information of the particles to update the electric field vector and current vector of the current time step, and then update the magnetic field vector of the current time step based on the current vector. In this way, through the above scheme, the induced electromagnetic field of charged particles can be numerically simulated more accurately.

[0072] In a specific embodiment, the electric field vector of the current time step is updated using the current motion information of the particle, including: determining the first charge of the charged particle at any target position based on the current displacement vector of the charged particle, and determining the second charge on different grids after the target turbulent field is gridded; determining the displacement vector diameter between the position where the first charge is located and the position where the second charge is located, so as to update the electric field vector of the current time step using the first charge, the second charge and the displacement vector diameter.

[0073] Specifically, the embodiment of the present application discloses the induced electric field equation:

[0074] ;

[0075] In the formula, Use the subscript object to represent any target position, that is, any coordinate point in the coordinate system The charge of the particle at , that is, the first charge of the charged particle, The subscript IJK represents the charge of particles on different grids after the target turbulent field is meshed, that is, the second charge of charged particles. represents the number of numerical grid divisions in three directions, It means along The index of the direction; express and The displacement vector between Represents the radial magnitude of the displacement vector.

[0076] That is, this application divides the entire calculation area into For each part, the electromagnetic field it is subjected to comes from the force exerted on it by all other parts around it. Therefore, by solving the above induced electromagnetic field equation, the electric field vector of the charged particle at the current time step can be obtained.

[0077] In a specific embodiment, the current vector of the current time step is updated using the current motion information of the particle, including: determining the velocity change rate of the charged particle in different directions based on the current velocity vector of the charged particle, and updating the current vector of the current time step using the charge of the charged particle and the velocity change rate.

[0078] Specifically, the present application discloses a current vector equation:

[0079] ;

[0080] The above equation shows that for each block, the current vector in the corresponding direction is equal to the product of the rate of change of the velocity of the charged particle and its charge q.

[0081] Specifically, the above-mentioned updating of the magnetic field vector of the current time step based on the current vector includes: determining a preset curl operator; and performing a curl operation on the current vector and the displacement vector based on the curl operator to update the magnetic field vector of the current time step.

[0082] Specifically, the embodiment of the present application discloses the magnetic field vector equation:

[0083] ;

[0084] That is, the present application uses the current vector and radial vector The curl can be used to obtain the magnetic field vector .

[0085] It can be seen that after obtaining the target fluid velocity vector of the charged particle position in the target turbulent field currently output by the preset fluid solver, the present application can solve the current motion information of the particle in the current time step based on the target fluid velocity vector and the particle historical motion information and historical electromagnetic field information updated in the previous time step; wherein the motion information specifically includes the velocity vector and the displacement vector, and the electromagnetic field information specifically includes the electric field vector and the magnetic field vector. Furthermore, since the position of the particle is updated, and the position change between different particles will affect the electric field force, it is necessary to update the electric field force information it is subjected to. In addition, the particle velocity changes, causing its current to change, which will also affect the magnetic field force. Therefore, the present application needs to use the particle current motion information to update the electric field vector and current vector of the current time step, and then update the magnetic field vector of the current time step based on the current vector. Through the above scheme, the induced electromagnetic field of charged particles can be numerically simulated more accurately.

[0086] In this way, the position and velocity of charged particles and the size of the electromagnetic field are finally obtained, so as to analyze the dispersion of particles or droplets in aircraft engines (closely related to combustion efficiency) and the start-up and development process of particles in sandstorms, as well as the maximum electromagnetic field strength that can be achieved. The above process is completed by code. As long as each step of the particle solver is coupled into the fluid solver, and then the communication of relevant information is completed, the final information obtained will also be output through the code, which is similar to a computing software.

[0087] See also Figure 3 As shown, the embodiment of the present application discloses a numerical simulation device for an induced electromagnetic field of charged particles, the device comprising:

[0088] A fluid velocity vector acquisition module 11 is used to acquire a target fluid velocity vector at a location where charged particles are located in a target turbulent field currently output by a preset fluid solver;

[0089] A motion information updating module 12 is used to solve the current motion information of the particle at the current time step based on the target fluid velocity vector and the particle historical motion information and historical electromagnetic field information updated at the previous time step; wherein the motion information includes a velocity vector and a displacement vector, and the electromagnetic field information includes an electric field vector and a magnetic field vector;

[0090] The electromagnetic field information updating module 13 is used to update the electric field vector and the current vector of the current time step using the current movement information of the particle, and calculate the magnetic field vector of the current time step based on the current vector.

[0091] It can be seen that after obtaining the target fluid velocity vector of the charged particle position in the target turbulent field currently output by the preset fluid solver, the present application can solve the current motion information of the particle in the current time step based on the target fluid velocity vector and the particle historical motion information and historical electromagnetic field information updated in the previous time step; wherein the motion information specifically includes the velocity vector and the displacement vector, and the electromagnetic field information specifically includes the electric field vector and the magnetic field vector. Furthermore, since the position of the particle is updated, and the position change between different particles will affect the electric field force, it is necessary to update the electric field force information it is subjected to. In addition, the particle velocity changes, causing its current to change, which will also affect the magnetic field force. Therefore, the present application needs to use the particle current motion information to update the electric field vector and current vector of the current time step, and then update the magnetic field vector of the current time step based on the current vector. Through the above scheme, the induced electromagnetic field of charged particles can be numerically simulated more accurately.

[0092] Figure 4 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application. Specifically, it may include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the numerical simulation method of the induced electromagnetic field of charged particles performed by the electronic device disclosed in any of the aforementioned embodiments.

[0093] In this embodiment, the power supply 23 is used to provide working voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device, and the communication protocol it follows is any communication protocol that can be applied to the technical solution of the present application, and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0094] Among them, the processor 21 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 21 can be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). The processor 21 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 21 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 21 may also include an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.

[0095] In addition, the memory 22, as a carrier for storing resources, can be a read-only memory, a random access memory, a disk or an optical disk, etc. The resources stored thereon include an operating system 221, a computer program 222 and data 223, etc. The storage method can be temporary storage or permanent storage.

[0096] Among them, the operating system 221 is used to manage and control the hardware devices and computer programs 222 on the electronic device 20, so as to realize the operation and processing of the massive data 223 in the memory 22 by the processor 21, which can be Windows, Unix, Linux, etc. In addition to including a computer program that can be used to complete the numerical simulation method of the charged particle induced electromagnetic field performed by the electronic device 20 disclosed in any of the aforementioned embodiments, the computer program 222 can further include a computer program that can be used to complete other specific tasks. In addition to data received by the electronic device and transmitted from an external device, the data 223 can also include data collected by its own input and output interface 25, etc.

[0097] Furthermore, an embodiment of the present application also discloses a computer-readable storage medium, in which a computer program is stored. When the computer program is loaded and executed by a processor, the steps of the numerical simulation method of the induced electromagnetic field of charged particles disclosed in any of the aforementioned embodiments are implemented.

[0098] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0099] Those skilled in the art may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented with electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art may use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of this application.

[0100] The steps of the method or algorithm described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a compact disc read-only memory (CD-ROM), or any other form of storage medium known in the art.

[0101] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0102] The above is a detailed introduction to the numerical simulation method, device, equipment and storage medium for the induced electromagnetic field of charged particles provided by the present invention. Specific examples are used in this article to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea; at the same time, for general technical personnel in this field, according to the idea of ​​the present invention, there will be changes in the specific implementation method and application scope. In summary, the content of this specification should not be understood as a limitation on the present invention.

Claims

1. A numerical simulation method for the induced electromagnetic field of charged particles, characterized in that: include: Obtain the target fluid velocity vector at the location of the charged particles in the target turbulent field currently output by the preset fluid solver; Based on the target fluid velocity vector and according to the particle historical motion information and historical electromagnetic field information updated in the previous time step, the current motion information of the particle in the current time step is solved; wherein the motion information includes a velocity vector and a displacement vector, and the electromagnetic field information includes an electric field vector and a magnetic field vector; Using the current movement information of the particle, the electric field vector and the current vector of the current time step are updated, and based on the current vector, the magnetic field vector of the current time step is updated; Wherein, updating the electric field vector of the current time step by using the current motion information of the particle includes: Determining a first charge amount of the charged particle at any target position based on a current displacement vector of the charged particle, and determining a second charge amount on different grids after the target turbulent field is gridded; Determine a displacement vector diameter between a position where the first charge quantity is located and a position where the second charge quantity is located, so as to update the electric field vector of the current time step by using the first charge quantity, the second charge quantity and the displacement vector diameter; The updating of the magnetic field vector of the current time step based on the current vector comprises: Determine the preset curl operator; Based on the curl operator, a curl operation is performed on the current vector and the displacement vector to update the magnetic field vector of the current time step; Updating the current vector of the current time step using the current movement information of the particle includes: The velocity change rate of the charged particle in different directions is determined based on the current velocity vector of the charged particle, and the current vector of the current time step is updated using the charge amount of the charged particle and the velocity change rate.

2. The numerical simulation method of charged particle induced electromagnetic field according to claim 1, characterized in that: Also includes: The fluid velocity vector at the center position of the charged particle is determined as the target fluid velocity vector at the position where the charged particle is located.

3. The numerical simulation method of charged particle induced electromagnetic field according to claim 1, characterized in that: The step of solving the particle current motion information at the current time step based on the target fluid velocity vector and according to the particle historical motion information and historical electromagnetic field information updated at the previous time step includes: Determine a charged particle dynamics equation constructed based on an electric field vector, a magnetic field vector, and a velocity vector and a displacement vector between adjacent time steps; wherein known parameters in the charged particle dynamics equation include particle historical motion information and historical electromagnetic field information updated in the previous time step; The target fluid velocity vector and the characteristic time of particle motion are substituted into the charged particle dynamics equation to solve the current velocity vector and the current displacement vector of the charged particle in the current time step.

4. The numerical simulation method of charged particle induced electromagnetic field according to claim 3, characterized in that: The process of determining the characteristic time of the particle motion includes: The characteristic time of motion of the particles is calculated based on the density, diameter and fluid dynamic viscosity of the charged particles.

5. A numerical simulation device for electromagnetic field induced by charged particles, characterized in that: include: A fluid velocity vector acquisition module is used to acquire a target fluid velocity vector at a location where charged particles are located in a target turbulent field currently output by a preset fluid solver; A motion information updating module, used to solve the current motion information of the particle at the current time step based on the target fluid velocity vector and the particle historical motion information and historical electromagnetic field information updated at the previous time step; wherein the motion information includes a velocity vector and a displacement vector, and the electromagnetic field information includes an electric field vector and a magnetic field vector; An electromagnetic field information updating module, used to update the electric field vector and current vector of the current time step using the current motion information of the particle, and calculate the magnetic field vector of the current time step based on the current vector; Wherein, the electromagnetic field information updating module is specifically used to determine the first charge of the charged particle at any target position based on the current displacement vector of the charged particle, and determine the second charge on different grids after the target turbulent field is gridded; determine the displacement vector radial between the position where the first charge is located and the position where the second charge is located, so as to update the electric field vector of the current time step using the first charge, the second charge and the displacement vector radial; The electromagnetic field information updating module is specifically used to determine a preset curl operator; based on the curl operator, perform a curl operation on the current vector and the displacement vector to update the magnetic field vector of the current time step; The electromagnetic field information updating module is specifically used to update the current vector of the current time step using the current motion information of the particle, including: determining the velocity change rate of the charged particle in different directions based on the current velocity vector of the charged particle, and updating the current vector of the current time step using the charge of the charged particle and the velocity change rate.

6. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor is used to execute the computer program to implement the steps of the method for numerical simulation of the induced electromagnetic field of charged particles as claimed in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: Used to store computer programs; wherein, when the computer program is executed by a processor, the steps of the numerical simulation method of the charged particle induced electromagnetic field as claimed in any one of claims 1 to 4 are implemented.

Images