A method and device for simulating and calculating spatial distribution of electromagnetic field based on Maxwell equation

Abstract

The application provides a Maxwell equation-based electromagnetic field space distribution simulation calculation method and device, and is realized based on a pre-constructed electromagnetic field space distribution calculation model; the electromagnetic field space distribution calculation model is a mathematical model for simulating high-altitude electromagnetic pulses in an explosion mode and describing electromagnetic field space distribution rules by using Maxwell equations; the application defines explosion parameters, defines an observation point position, defines simulation calculation parameters, and then performs simulation calculation by using the electromagnetic field space distribution calculation model based on the defined various types of parameters, so as to obtain space electric field distribution data at the observation point position. The electromagnetic field space distribution calculation model is constructed by simulating high-altitude electromagnetic pulses in an explosion mode and according to Maxwell equations, and the model can be used to obtain space electric field distribution data at the observation point, and can provide a strong basis for the design and optimization of related protection schemes.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of high-altitude electromagnetic pulse environment electric field distribution calculation, and particularly relates to an electromagnetic field space distribution simulation calculation method and device based on Maxwell equations. BACKGROUND

[0002] High-altitude electromagnetic pulses have characteristics such as large electric field intensity, fast front rising speed, and wide frequency spectrum range, can reach a peak value in a very short time, generate thousands of volts of transient voltage on bare cables or be coupled into equipment, and cause irreversible damage to power electronic equipment. The environment is divided into three stages: early, middle and late. The early stage contains medium frequency to ultra-high frequency band signals, has wide radiation, large intensity and wide spectrum, can invade intelligent electronic devices such as data acquisition and monitoring control systems, and cause insulator flashover. The middle stage is similar to lightning electromagnetic pulses, and can destroy power main equipment after the early stage environment removes protective equipment. The late stage is similar to geomagnetic storms, and threatens low-frequency equipment and long-distance power transmission and distribution lines due to long duration. With the development of smart grids, high-precision electronic devices are widely used, which are sensitive to electromagnetic pulses and are easily damaged. Key electrical connections, gaps and cables in power systems are main coupling points of electromagnetic pulses. Power cables, as the connecting link of equipment, are exposed to electromagnetic pulse environments, which can couple high voltage and current, invade power networks, threaten end devices or damage cables and disrupt system stable operation.

[0003] In order to ensure the stable operation of the power system, it is necessary to study the coupling effect of key equipment in the high-altitude electromagnetic pulse environment to design a protection scheme. However, it is difficult to effectively carry out related research and protection scheme development work in actual production due to the inability to generate a real high-altitude electromagnetic pulse environment. Therefore, an effective simulation technology is urgently needed to simulate the high-altitude electromagnetic pulse environment and its influence on the power system, so as to provide a basis and support for protection scheme design. SUMMARY

[0004] Therefore, the application provides an electromagnetic field space distribution simulation calculation method and device based on Maxwell equations, which aims to simulate the generation and propagation of high-altitude electromagnetic pulse environment by numerical methods, and is used for calculating the electric field distribution of high-altitude electromagnetic pulse environment at a certain point.

[0005] In order to achieve the above-mentioned purpose, the technical scheme provided by the application is as follows:

[0006] In a first aspect, the application provides an electromagnetic field space distribution simulation calculation method based on Maxwell equations, which is realized based on a pre-constructed electromagnetic field space distribution calculation model. The electromagnetic field space distribution calculation model is a mathematical model for simulating high-altitude electromagnetic pulses in an explosive manner and describing the electromagnetic field space distribution law by using Maxwell equations.

[0007] The simulation calculation method comprises the following steps:

[0008] define explosion parameters;

[0009] define observation point position; the observation point position is a point in the high-altitude electromagnetic pulse simulation scene;

[0010] define simulation calculation parameters;

[0011] Based on the defined various types of parameters, the electromagnetic field space distribution calculation model is used for simulation calculation to obtain the spatial electric field distribution data at the observation point position.

[0012] Further, the mathematical expression of the electromagnetic field space distribution calculation model is as follows:

[0013]

[0014]

[0015] In the formula, is a distance variable, is a component of the magnetic field intensity, is a component of the magnetic field intensity, is a component of the current density, is a component of the current density, is an electrical conductivity, is a dielectric constant, is a magnetic permeability, is a time variable. Further, for the electromagnetic field space distribution calculation model, the equal phase method of delay time is used for conversion, and based on the converted electromagnetic field space distribution calculation model, the spatial electric field distribution data at the observation point position is obtained.

[0016] Further, the mathematical expression of the converted electromagnetic field space distribution calculation model is as follows:

[0017]

[0018]

[0019]

[0020]

[0021] In the formula, and are right and left traveling waves respectively, is a distance variable for variable substitution by the equal phase method of delay time, is a time variable for variable substitution by the equal phase method of delay time.

[0022] ​​Further, in solving the converted electromagnetic field spatial distribution calculation model, further comprising:

[0023] Discretize the converted electromagnetic field spatial distribution calculation model into difference equations.

[0024] In a second aspect, the application provides an electromagnetic field spatial distribution simulation calculation device based on Maxwell equation, which is realized based on a pre-constructed electromagnetic field spatial distribution calculation model, the electromagnetic field spatial distribution calculation model is a mathematical model for simulating high-altitude electromagnetic pulse in an explosion mode and describing electromagnetic field spatial distribution law by using Maxwell equation;

[0025] The simulation calculation device comprises:

[0026] A parameter definition module, configured to define explosion parameters, define an observation point position, the observation point position being a point in a high-altitude electromagnetic pulse simulation scene, and define simulation calculation parameters.

[0027] A simulation calculation module, configured to perform simulation calculation by using the electromagnetic field spatial distribution calculation model based on the defined parameters of each type, and obtain spatial electric field distribution data at the observation point position.

[0028] Further, the mathematical expression of the electromagnetic field spatial distribution calculation model is as follows:

[0029]

[0030]

[0031] In the formula, is a distance variable, is a component of magnetic field intensity, is a component of magnetic field intensity, is a component of current density, is electric conductivity, is dielectric constant, is magnetic permeability, is a time variable.

[0032] Further, for the electromagnetic field spatial distribution calculation model, the equal phase method of time delay is used for conversion, and the converted electromagnetic field spatial distribution calculation model is used for solving, so as to obtain the spatial electric field distribution data at the observation point position.

[0033] Further, the mathematical expression of the converted electromagnetic field spatial distribution calculation model is as follows:

[0034]

[0035]

[0036]

[0037] In the formula, and are right and left traveling waves respectively, is a distance variable for variable substitution by using the equal phase method of delay time, is a time variable for variable substitution by using the equal phase method of delay time.

[0038] Further, in solving the converted electromagnetic field spatial distribution calculation model, the method further comprises:

[0039] discretizing the converted electromagnetic field spatial distribution calculation model into difference equations.

[0040] In summary, the application provides an electromagnetic field spatial distribution simulation calculation method and device based on Maxwell equation, which is realized based on a pre-constructed electromagnetic field spatial distribution calculation model, the electromagnetic field spatial distribution calculation model is a mathematical model for simulating high-altitude electromagnetic pulse in an explosion mode and describing electromagnetic field spatial distribution law by using Maxwell equation; the application defines explosion parameters, defines an observation point position, the observation point position is a point in a high-altitude electromagnetic pulse simulation scene, defines simulation calculation parameters, and then performs simulation calculation by using the electromagnetic field spatial distribution calculation model based on the defined various types of parameters to obtain spatial electric field distribution data at the observation point position. The electromagnetic field spatial distribution calculation model of the application simulates high-altitude electromagnetic pulse in an explosion mode and is constructed according to Maxwell equation, the model can obtain spatial electric field distribution data at the observation point, and can provide a strong basis for design and optimization of relevant protection schemes. BRIEF DESCRIPTION OF DRAWINGS

[0041] In order to more clearly illustrate the technical solutions of the embodiments of the application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or prior art description. Obviously, the drawings in the following description are only some embodiments of the application, and for those skilled in the art, other drawings can also be obtained without creative labor based on these drawings.

[0042] Figure 1 A flowchart of an electromagnetic field spatial distribution simulation calculation method based on Maxwell equation provided by the embodiment of the application is shown in the figure.

[0043] Figure 2 A schematic diagram of an observation point provided by the embodiment of the application is shown in the figure.

[0044] Figure 3 A spatial electric field distribution diagram of an observation point provided by the embodiment of the application is shown in the figure.

[0045] Figure 4 A composition block diagram of an electromagnetic field space distribution simulation calculation device based on Maxwell equations is provided for an embodiment of the present application.

[0046] Figure 5 A composition block diagram of a computer device is provided for an embodiment of the present application. DETAILED DESCRIPTION

[0047] In order to make the purpose, features and advantages of the present application more obvious and easy to understand, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the embodiments described below are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0048] Please refer to Figure 1 The present embodiment provides a simulation calculation method of electromagnetic field space distribution based on Maxwell equations, which is realized based on a pre-constructed electromagnetic field space distribution calculation model. The electromagnetic field space distribution calculation model is a mathematical model that simulates high-altitude electromagnetic pulses in an explosion manner and describes the mathematical model of electromagnetic field space distribution law by using Maxwell equations.

[0049] The electromagnetic field space distribution calculation model is a mathematical model that simulates high-altitude electromagnetic pulses in an explosion manner and describes the mathematical model of electromagnetic field space distribution law according to Maxwell equations. Maxwell equations are the basis of classical electromagnetism theory, which comprehensively describes the basic laws of electromagnetic field generation, propagation, interaction, etc. The generation process of high-altitude electromagnetic pulses is simulated in an explosion manner. When the explosion occurs, a huge amount of energy is released instantaneously, causing physical changes such as ionization of the surrounding medium (such as air), thereby generating electromagnetic pulses that propagate outward. By setting the explosion equivalent, the energy intensity of the simulated electromagnetic pulse source can be controlled, and by determining the explosion position, the starting propagation point of the electromagnetic pulse in space can be determined.

[0050] The simulation based on this kind of way can more realistically reflect the source situation of high-altitude electromagnetic pulses, and provides initial conditions that conform to actual physical scenarios for subsequent accurate calculation of electromagnetic field space distribution. Different explosion equivalent and position settings will result in different electromagnetic pulse characteristics (such as intensity, propagation direction, etc.), and by reasonably setting these parameters, the corresponding electromagnetic pulse environment can be simulated according to specific research needs.

[0051] The simulation calculation method includes the following steps:

[0052] Step one: define the explosion parameters.

[0053] It should be noted that in this step, the defined explosion parameters mainly refer to parameters used to describe the relevant characteristics of the explosion event, such as the explosion equivalent, explosion location, etc. The explosion equivalent is an indicator to measure the size of the energy released by the explosion, usually in units of kiloton (kt) or the like; the explosion location clearly indicates the specific location of the explosion in space.

[0054] Step two: define the observation point position; the observation point position is a point in the high-altitude electromagnetic pulse simulation scenario.

[0055] It should be noted that when studying the spatial distribution of electromagnetic fields, it is usually not possible to comprehensively and meticulously monitor the electromagnetic field situation in the entire space in real time, so it is necessary to select representative specific points for in-depth analysis. This observation point position is like a monitoring station set up in the actual environment, by focusing on this point, the influence of the electromagnetic pulse on the electromagnetic field when it propagates to this point can be understood.

[0056] The observation point position is a specific point selected in the simulated high-altitude electromagnetic pulse scenario, used to focus on and analyze the electromagnetic field-related characteristics at this point, such as spatial electric field distribution, etc.

[0057] Step three: define the simulation calculation parameters.

[0058] It should be noted that the simulation calculation parameters cover other related parameters in addition to the explosion parameters and the observation point position in the process of simulating the spatial distribution of electromagnetic fields, which may include the setting of calculation duration, calculation accuracy requirements, calculation methods (such as whether to use a certain numerical calculation method, whether to perform approximation, etc.) and other aspects.

[0059] Different simulation calculation parameter settings will affect the process and results of the calculation. For example, the calculation duration determines the time range of the simulation of electromagnetic pulse propagation and electromagnetic field changes, a longer calculation duration can more comprehensively capture the influence of electromagnetic pulses at different stages, but it may also increase the calculation cost and time consumption; the calculation accuracy requirement is related to the accuracy of the calculation results, higher accuracy requirements usually require more refined calculation methods and more calculation resources.

[0060] By reasonably defining these parameters, a balance between calculation efficiency and calculation result accuracy can be found according to the research purpose and actual calculation resources, etc., to ensure that accurate electromagnetic field spatial distribution calculation results that meet the research requirements can be obtained within an acceptable time.

[0061] Step four: based on the defined various types of parameters, use the electromagnetic field spatial distribution calculation model to perform simulation calculation, and obtain the spatial electric field distribution data at the observation point position.

[0062] It should be noted that, based on the previously defined explosion parameters, observation point locations, and simulation calculation parameters, this information about the actual physical scenario and computational requirements is input into the electromagnetic field spatial distribution calculation model. This model utilizes Maxwell's equations and, through a series of mathematical operations (such as variable substitution of electromagnetic field-related physical quantities, derivation of isophase equations, and discretization of continuous equations into difference equations), simulates the propagation of electromagnetic pulses in space and the changing process of the electromagnetic field. In this process, based on the relationship between electric and magnetic fields in Maxwell's equations and the spatiotemporal variation law of the electromagnetic field, the electromagnetic field state at each location (especially the observation point location) at different time points is calculated step by step.

[0063] Through the above simulation calculations, accurate spatial electric field distribution data at the observation point can be obtained. This data can intuitively reflect the spatial distribution of the electric field when the electromagnetic pulse propagates to the observation point, such as the magnitude and direction of the electric field strength. This is of great significance for studying the coupling effect of electromagnetic pulses on equipment and assessing the safety of equipment in electromagnetic pulse environments.

[0064] This embodiment provides a simulation method for the spatial distribution of electromagnetic fields based on Maxwell's equations. Based on the fundamental laws of electromagnetic fields described by the classical Maxwell's equations, and combined with a simulation of the generation process of high-altitude electromagnetic pulses via explosion, it integrates the actual physical scenario and computational requirements into a constructed electromagnetic field spatial distribution calculation model by reasonably defining various parameters related to the generation, propagation, and calculation of electromagnetic pulses (explosion parameters, observation point location, simulation calculation parameters, etc.). Then, a series of mathematical operations are performed using this model to simulate the propagation of electromagnetic pulses in space and the changes in the electromagnetic field, thereby calculating the spatial electric field distribution data at the observation point location, thus enabling simulation studies of the spatial distribution of electromagnetic fields under specific scenarios.

[0065] In one embodiment, the mathematical expression for the electromagnetic field spatial distribution calculation model is as follows:

[0066] (1)

[0067] (2)

[0068] In the formula, For distance variables, For magnetic field strength Quantity, For magnetic field strength Quantity, For current density Quantity, For electrical conductivity, Where is the dielectric constant. Permeability, is a time variable.

[0069] In this embodiment, the distribution of the electric field intensity in space can be inferred from the changes in the magnetic field by using the derivatives of the magnetic field intensity in space, the derivatives of the electric field intensity in space and time, and other physical quantities. The correlation between the rate of change of the electric field intensity in space and the rate of change of the magnetic field over time is also considered. Given boundary conditions such as explosion parameters and observation point positions, the spatial electric field distribution can be obtained by solving the equations.

[0070] In further embodiments, for the electromagnetic field spatial distribution calculation model, a time-delayed equal-phase method is used for conversion, and based on the converted electromagnetic field spatial distribution calculation model, the spatial electric field distribution data at the observation point position is obtained by solving.

[0071] In this embodiment, the time-delayed equal-phase method is an effective method commonly used in the field of electromagnetic field calculation. The core idea is to simplify the calculation process and improve the calculation efficiency through specific mathematical transformation, while ensuring that the required electromagnetic field related data can be accurately obtained.

[0072] In dealing with complex problems such as electromagnetic field spatial distribution calculation, the original electromagnetic field spatial distribution calculation model (such as the equation form mentioned earlier based on Maxwell's equations) may face difficulties such as large amount of calculation and long calculation time in direct calculation, especially when dealing with large-scale spatial regions or complex electromagnetic pulse scenarios (such as simulating high-altitude electromagnetic pulse environments). By using the time-delayed equal-phase method for conversion, the spatial electric field distribution data at the observation point position can be obtained more efficiently while ensuring the accuracy of the calculation results.

[0073] In further embodiments, the process of using the time-delayed equal-phase method for conversion is as follows:

[0074] Substitute variables for equations (1) and (2) Equation (1) becomes:

[0075] (3)

[0076] Equation (2) becomes:

[0077]

[0078] (3) + (4):

[0079] (5)

[0080] (3) - (4):

[0081]

[0082] Let: , . We can get:

[0083]

[0084] where is the right-going wave, is the left-going wave. Substituting and into equations (5) and (6), we get the equal-phase equation of the transverse magnetic wave:

[0085]

[0086] In further embodiments, in solving the converted electromagnetic field spatial distribution calculation model, further comprising: discretizing the converted electromagnetic field spatial distribution calculation model into difference equations.

[0087] In this embodiment, after the equal-phase method with time delay is used to convert the electromagnetic field spatial distribution calculation model in the previous, the form of the equation has been optimized to a certain extent, which is more convenient to reflect the equal-phase characteristics in the propagation process of electromagnetic pulses. On this basis, the operation of discretizing into difference equations can further utilize this optimized equation structure. For example, in the equation converted by the equal-phase method, there may be some variables and terms related to the time delay. In the discretization process, these terms can be better handled, so that the discretized difference equation can more accurately reflect the propagation characteristics of electromagnetic pulses at discrete space and time points.

[0088] By discretizing into difference equations, numerical calculation methods (such as iterative method, direct solution method, etc.) can be used to solve these equations, thereby obtaining the values of electromagnetic field physical quantities at each discrete point. For the observation point position, since it is already one of the discrete space points (or an approximate value of the observation point position can be obtained from the values of the discrete points through interpolation, etc.), the spatial electric field distribution data at the observation point position can be obtained. This method can effectively handle complex electromagnetic field problems, and the accuracy of spatial and temporal discretization can be adjusted as needed to balance the calculation efficiency and the accuracy of the calculation results.

[0089] The simulation calculation method provided by the above embodiments will be introduced below in combination with an example.

[0090] 1. Define the explosion parameters: define the explosion equivalent as 1000 kt, and the explosion position as the 0 point in Figure 2 .

[0091] Figure 2The distance relationship between the observation point and the source point (e.g., the center of an electric dipole) is shown, with the radius r and angle a reflecting the range and direction of the explosion's influence.

[0092] 2. Set the observation point method to "Observation Point Location".

[0093] 3. Define the "observation point location", and select the observation point as... Figure 2 Point A in the diagram.

[0094] 4. Select "Single Observation Point Calculation" as the calculation mode;

[0095] 5. Define the total computation time as 100ns.

[0096] 6. Calculate the spatial electric field distribution at the observation point, such as... Figure 3 As shown.

[0097] Based on the same inventive concept, this application also provides a Maxwell's equations-based electromagnetic field spatial distribution simulation calculation device for implementing the Maxwell's equations-based electromagnetic field spatial distribution simulation calculation method described above. The solution provided by this device is similar to the solution described in the above method; therefore, the specific limitations in the embodiments of the Maxwell's equations-based electromagnetic field spatial distribution simulation calculation device provided below can be found in the limitations of the Maxwell's equations-based electromagnetic field spatial distribution simulation calculation method described above, and will not be repeated here.

[0098] Please see Figure 4 The present invention also provides an electromagnetic field spatial distribution simulation calculation device based on Maxwell's equations, which is based on a pre-constructed electromagnetic field spatial distribution calculation model. The electromagnetic field spatial distribution calculation model is a mathematical model that simulates high-altitude electromagnetic pulses in the form of explosions and uses Maxwell's equations to describe the electromagnetic field spatial distribution law.

[0099] The simulation computing device includes:

[0100] The parameter definition module is used to define explosion parameters; it is also used to define the observation point location, which is a point in the high-altitude electromagnetic pulse simulation scenario; and it is also used to define simulation calculation parameters.

[0101] The simulation calculation module is used to perform simulation calculations based on the defined parameters and the electromagnetic field spatial distribution calculation model to obtain the spatial electric field distribution data at the observation point.

[0102] Furthermore, the mathematical expression for the electromagnetic field spatial distribution calculation model is as follows:

[0103]

[0104]

[0105] wherein, is a distance variable, is a component of magnetic field intensity, is a component of magnetic field intensity, is a component of magnetic field intensity, is a component of magnetic field intensity, is a component of current density, is a component of current density, is a conductivity, is a dielectric constant, is a permeability, is a time variable.

[0106] Further, for the electromagnetic field spatial distribution calculation model, a delayed time equal phase method is used for conversion, and based on the converted electromagnetic field spatial distribution calculation model, a solution is obtained to obtain spatial electric field distribution data at an observation point position.

[0107] Further, the mathematical expression of the converted electromagnetic field spatial distribution calculation model is as follows:

[0108]

[0109]

[0110]

[0111] wherein, and are right and left traveling waves respectively, is a distance variable for variable replacement by the delayed time equal phase method, is a time variable for variable replacement by the delayed time equal phase method.

[0112] Further, in solving the converted electromagnetic field spatial distribution calculation model, it also includes:

[0113] discretizing the converted electromagnetic field spatial distribution calculation model into difference equations.

[0114] Those skilled in the art can clearly understand that, for the convenience and brevity of description, only the division of the above functional units and modules is exemplified, and in actual application, the above functions can be completed by different functional units and modules according to needs, that is, the internal structure of the system is divided into different functional units or modules to complete all or part of the functions described above. Each functional unit or module in the embodiment can be integrated in one processing unit, or each unit can exist physically, or two or more units can be integrated in one unit. The integrated unit can be realized in the form of hardware or in the form of software functional unit. In addition, the specific names of each functional unit or module are only for convenient distinction, and do not limit the protection scope of the present application. The specific working process of the unit or module in the system can refer to the corresponding process in the foregoing method embodiment, which will not be repeated here.

[0115] With reference to Figure 5 The embodiment of the present application also provides a computer device, comprising a memory and a processor and a computer program stored in the memory, when the computer program is executed on the processor, the Maxwell equation based electromagnetic field space distribution simulation calculation method is realized.

[0116] The computer device can be a desktop computer, a notebook computer, a palm computer and a cloud server and the like. The computer device can include, but is not limited to, a processor and a memory. Those skilled in the art can understand that, Figure 5 It is only an example of the computer device, and does not constitute a limitation on the computer device, and can include more or fewer components than the illustration, or combine certain components, or different components, for example, it can also include an input and output device, a network access device and the like.

[0117] The processor can be a central processing unit (CPU), and the processor can also be other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor.

[0118] The memory can be an internal storage unit of the computer device in some embodiments, such as a hard disk or a memory of the computer device. The memory can also be an external storage device of the computer device in other embodiments, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, and the like. Further, the memory can include both an internal storage unit and an external storage device of the computer device. The memory is used to store an operating system, an application program, a boot loader, data, and other programs, such as program codes of the computer program. The memory can also be used to temporarily store data that has been output or is to be output.

[0119] The embodiment of the present application further provides a computer readable storage medium, which stores a computer program. The computer program is run by a processor to implement the electromagnetic field space distribution simulation calculation method based on Maxwell equations according to any one of the above methods.

[0120] In the embodiment, the integrated unit, if implemented in the form of a software functional unit and sold or used as an independent product, can be stored in a computer readable storage medium. Based on such understanding, the present application can implement all or part of the processes of the above-mentioned embodiment methods by a computer program to instruct related hardware to complete, and the computer program can be stored in a computer readable storage medium. The computer program is executed by a processor to implement the steps of each method embodiment. The computer program includes computer program codes, which can be in the form of source code, object code, executable files or some intermediate forms, etc. The computer readable medium at least includes any entity or device capable of carrying the computer program code to a photographing device / terminal device, a recording medium, a computer memory, a read-only memory (ROM), a random access memory (RAM), an electrical carrier signal, a telecommunications signal and a software distribution medium. For example, a U disk, a mobile hard disk, a magnetic disk or an optical disk, etc. In some jurisdictions, according to legislation and patent practice, the computer readable medium can not be an electrical carrier signal and a telecommunications signal.

[0121] In the above embodiments, the description of each embodiment has its own focus, and the parts not described or recorded in detail in a certain embodiment can be referred to the relevant description of other embodiments.

[0122] Those skilled in the art can understand that the units and algorithm steps of each example described in combination with the embodiments disclosed herein can be realized by electronic hardware or a combination of computer software and electronic hardware. Whether the functions are realized in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to realize the described functions for each specific application, but such implementation should not be considered beyond the scope of the present application.

[0123] In the embodiments disclosed in the present application, it should be understood that the disclosed apparatus / terminal device and method can be implemented in other ways. For example, the apparatus / terminal device embodiments described above are merely schematic, for example, the division of the modules or units is merely a logical function division, and actual implementation can have another division manner, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the displayed or discussed each other can be indirect coupling or communication connection through some interfaces, devices or units, and can be electrical, mechanical or other forms.

[0124] The above embodiments are only used to illustrate the technical solutions of the present application, rather than limit them; although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that the technical solutions recorded in the foregoing embodiments can be modified, or some technical features can be replaced by equivalents; and these modifications or replacements do not make the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for simulating Maxwell equation-based electromagnetic field spatial distribution, characterized in that, Based on a pre-constructed electromagnetic field space distribution calculation model, the electromagnetic field space distribution calculation model is a mathematical model for simulating high-altitude electromagnetic pulse in an explosion manner and describing electromagnetic field space distribution law by using Maxwell equation; The simulation calculation method comprises the following steps: Defining explosion parameters; Defining an observation point position, which is a point in a high-altitude electromagnetic pulse simulation scene; Defining simulation calculation parameters; Based on the defined various types of parameters, the electromagnetic field space distribution calculation model is used for simulation calculation to obtain spatial electric field distribution data at the observation point position; For the electromagnetic field space distribution calculation model, a delayed time equal phase method is used for conversion, and the converted electromagnetic field space distribution calculation model is used for solving to obtain the spatial electric field distribution data at the observation point position; The mathematical expression of the converted electromagnetic field space distribution calculation model is as follows: ； ； ； wherein and are right and left traveling waves, respectively, is a distance variable for variable substitution using the equal phase method with retardation time, is a time variable for variable substitution using the equal phase method with retardation time; is a component of magnetic field intensity; is electric conductivity; is dielectric constant; is a component of magnetic field intensity; is a component of current density.

2. The electromagnetic field spatial distribution simulation calculation method based on Maxwell equations according to claim 1, characterized in that, The mathematical expression of the electromagnetic field space distribution calculation model is as follows: ； ； wherein is a distance variable, is a magnetic permeability, is a time variable.

3. The electromagnetic field spatial distribution simulation calculation method based on Maxwell equations according to claim 1, characterized in that, When solving the converted electromagnetic field space distribution calculation model, the following steps are further included: The converted electromagnetic field space distribution calculation model is discretized into a difference equation.

4. A device for simulating a spatial distribution of an electromagnetic field based on Maxwell's equations, characterized by Based on a pre-constructed electromagnetic field space distribution calculation model, the electromagnetic field space distribution calculation model is a mathematical model for simulating high-altitude electromagnetic pulse in an explosion manner and describing electromagnetic field space distribution law by using Maxwell equation; The simulation calculation device comprises: A parameter definition module, configured to define explosion parameters, define an observation point position, which is a point in a high-altitude electromagnetic pulse simulation scene, and define simulation calculation parameters; A simulation calculation module, configured to, based on the defined various types of parameters, use the electromagnetic field space distribution calculation model for simulation calculation to obtain spatial electric field distribution data at the observation point position; For the electromagnetic field space distribution calculation model, a delayed time equal phase method is used for conversion, and the converted electromagnetic field space distribution calculation model is used for solving to obtain the spatial electric field distribution data at the observation point position; The mathematical expression of the converted electromagnetic field space distribution calculation model is as follows: ； ； ； wherein and are right and left traveling waves, respectively, is a distance variable for variable substitution using the equal phase method of retardation time, is a time variable for variable substitution using the equal phase method of retardation time; is a component of magnetic field intensity; is electrical conductivity; is dielectric constant; is a component of magnetic field intensity; is a component of current density.

5. The Maxwell equation based electromagnetic field space distribution emulation computing device of claim 4, wherein, The mathematical expression of the electromagnetic field space distribution calculation model is as follows: ； ； wherein is a distance variable, is a magnetic permeability, is a time variable.

6. The Maxwell equation based electromagnetic field space distribution emulation computing device of claim 4, wherein, When solving the converted electromagnetic field space distribution calculation model, the following steps are further included: The converted electromagnetic field space distribution calculation model is discretized into a difference equation.

Images