CN104915501A - Method for simulating electromagnetic distribution of strong electromagnetic pulse environment - Google Patents
Method for simulating electromagnetic distribution of strong electromagnetic pulse environment Download PDFInfo
- Publication number
- CN104915501A CN104915501A CN201510329815.9A CN201510329815A CN104915501A CN 104915501 A CN104915501 A CN 104915501A CN 201510329815 A CN201510329815 A CN 201510329815A CN 104915501 A CN104915501 A CN 104915501A
- Authority
- CN
- China
- Prior art keywords
- electromagnetic pulse
- electromagnetic
- time
- domain
- finite difference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 10
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000004088 simulation Methods 0.000 claims description 28
- 238000004458 analytical method Methods 0.000 claims description 2
- 230000005284 excitation Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 8
- 230000001133 acceleration Effects 0.000 abstract description 5
- 230000006870 function Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种仿真强电磁脉冲环境下的电磁分布情况仿真方法。该方法使用双曲余弦函数模拟时域下的电磁脉冲波形，使用时域有限差分算法对环境的电磁场变化进行仿真模拟，使用GPU实现时域有限差分方法的加速，得到时域上的电磁分布情况，最后通过快速傅里叶变换，得到频率域上的电磁分布情况，解决了电磁脉冲武器攻击下超电大尺寸对象的电磁环境分布特性计算问题。
The invention proposes a method for simulating electromagnetic distribution in a strong electromagnetic pulse environment. This method uses the hyperbolic cosine function to simulate the electromagnetic pulse waveform in the time domain, uses the finite difference time domain algorithm to simulate the electromagnetic field changes in the environment, uses the GPU to realize the acceleration of the time domain finite difference method, and obtains the electromagnetic distribution in the time domain , and finally through the fast Fourier transform, the electromagnetic distribution in the frequency domain is obtained, which solves the calculation problem of the electromagnetic environment distribution characteristics of the ultra-electric large-scale objects under the attack of electromagnetic pulse weapons.
Description
技术领域technical field
本发明涉及计算机仿真中的强电磁脉冲模拟以及并行时域有限差分电磁场仿真技术,解决了电磁脉冲武器攻击下超电大尺寸对象的电磁环境分布特性计算问题。The invention relates to strong electromagnetic pulse simulation and parallel time-domain finite difference electromagnetic field simulation technology in computer simulation, and solves the problem of calculating the electromagnetic environment distribution characteristics of superelectric large-scale objects under the attack of electromagnetic pulse weapons.
背景技术Background technique
电磁脉冲具有宽频谱和高能量特性，能通过线缆、缝隙等耦合到电子系统中，导致电子系统性能降级甚至损坏。因此，研究复杂电子系统所处的电磁脉冲环境，对提高复杂电子系统对电磁脉冲的防护能力具有重要的意义。Electromagnetic pulses have wide spectrum and high energy characteristics, and can be coupled into electronic systems through cables, gaps, etc., causing performance degradation or even damage to electronic systems. Therefore, it is of great significance to study the electromagnetic pulse environment in which complex electronic systems are located to improve the protection ability of complex electronic systems against electromagnetic pulses.
时域有限差分方法是把Maxwell方程式在时间和空间领域上进行差分化。利用蛙跳式方式在空间领域内的电场和磁场进行交替计算，通过时间领域上更新来模仿电磁场的变化，达到数值计算的目的。其优点是能够直接模拟场的分布，精度比较高，是目前使用比较多的数值模拟的方法之一。The finite difference time domain method is to differentiate the Maxwell equation in the time and space fields. Alternately calculate the electric field and magnetic field in the space domain by using leapfrog method, and simulate the change of electromagnetic field by updating in the time domain to achieve the purpose of numerical calculation. Its advantage is that it can directly simulate the distribution of the field, and its accuracy is relatively high. It is one of the most widely used numerical simulation methods at present.
考虑到时域有限差分本身算法的限制，仿真范围越大、电磁波频率越高，需要的计算量以及内存使用量越大。而近些年出现的GPU加速方式，可以充分发掘时域有限差分本身的并行特质，使计算速度获得成百倍的提升。Considering the limitations of the finite difference time domain algorithm itself, the larger the simulation range and the higher the frequency of electromagnetic waves, the greater the amount of calculation and memory usage required. The GPU acceleration method that has appeared in recent years can fully exploit the parallel characteristics of the finite difference in time domain itself, and the calculation speed has been increased by hundreds of times.
发明内容Contents of the invention
本发明使用时域有限差分算法，实现强电磁脉冲环境下的电磁情况仿真，并通过GPU加速仿真过程，实现超电大尺寸目标的电磁分布仿真。The invention uses a time-domain finite difference algorithm to realize the simulation of electromagnetic conditions in a strong electromagnetic pulse environment, and accelerates the simulation process through a GPU to realize the simulation of electromagnetic distribution of super-electric large-scale targets.
本发明所采用的技术方案：The technical scheme adopted in the present invention:
一种强电磁脉冲环境电磁分布仿真方法,其特征是：使用时域有限差分方法作为电磁场仿真算法进行仿真,建立目标区域时频空多维分析模型。A method for simulating electromagnetic distribution in a strong electromagnetic pulse environment is characterized by: using the time-domain finite difference method as an electromagnetic field simulation algorithm for simulation, and establishing a time-frequency-space multidimensional analysis model of a target area.
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法，其特征是：使用GPU实现时域有限差分方法的并行加速。The method for simulating electromagnetic distribution in a strong electromagnetic pulse environment according to claim 1, characterized in that: GPU is used to realize the parallel acceleration of the time-domain finite difference method.
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法，其特征是：包含电磁脉冲生成函数、基于GPU的时域有限差分方法、快速傅里叶变换三部分。The electromagnetic distribution simulation method in a strong electromagnetic pulse environment according to claim 1, characterized in that it includes three parts: an electromagnetic pulse generation function, a GPU-based finite difference time domain method, and a fast Fourier transform.
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法，其特征是：使用双曲余弦函数，模拟核电磁脉冲时域波形，作为仿真的信号激励。A method for simulating electromagnetic distribution in a strong electromagnetic pulse environment according to claim 1, characterized in that: hyperbolic cosine function is used to simulate the nuclear electromagnetic pulse time-domain waveform as the signal excitation for simulation.
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法，其特征是：使用时域有限差分方法，进行空间中的电磁波模拟，得到空间中电磁场的分布的时域情况。A method for simulating electromagnetic distribution in a strong electromagnetic pulse environment according to claim 1, characterized in that: the method of finite difference in time domain is used to simulate electromagnetic waves in space to obtain the time domain situation of the distribution of electromagnetic fields in space.
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法，其特征是：使用GPU实现对时域有限差分算法的并行加速，减少仿真超电大尺寸目标及高频脉冲时的时间。The electromagnetic distribution simulation method in a strong electromagnetic pulse environment according to claim 1, characterized in that: GPU is used to realize the parallel acceleration of the time-domain finite difference algorithm, and the time for simulating super-electric large-scale targets and high-frequency pulses is reduced.
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法，其特征是：，使用快速傅里叶变换，从电磁场分布的时域情况得到空间电磁场分布的频域情况。A method for simulating electromagnetic distribution in a strong electromagnetic pulse environment according to claim 1, characterized in that: use fast Fourier transform to obtain the frequency domain situation of the spatial electromagnetic field distribution from the time domain situation of the electromagnetic field distribution.
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法，其特征是：时域有限差分方法所使用的吸收边界为卷积完美匹配层(CPML)。A method for simulating electromagnetic distribution in a strong electromagnetic pulse environment according to claim 1, characterized in that: the absorption boundary used by the finite difference time domain method is a convolutional perfectly matched layer (CPML).
本发明相对于现有技术具有如下的优点及特性：Compared with the prior art, the present invention has the following advantages and characteristics:
(1)从Maxwell方程直接推导出的时域有限差分算法，仿真精度高。(1) The time-domain finite difference algorithm directly derived from the Maxwell equation has high simulation accuracy.
(2)使用GPU进行时域有限差分方法的并行实现，仿真速度快。(2) Parallel implementation of the time-domain finite-difference method using GPU, and the simulation speed is fast.
(4)使用快速傅里叶变换，可以从时域结果中直接得到频域的结果。(4) Using Fast Fourier Transform, the frequency domain result can be obtained directly from the time domain result.
附图说明Description of drawings
图1是本发明所述强电磁脉冲环境电磁分布仿真方法流程图；Fig. 1 is the flow chart of electromagnetic distribution simulation method in strong electromagnetic pulse environment of the present invention;
图2是本发明所述强电磁脉冲环境电磁分布仿真方法计算空间划分示意图；Fig. 2 is a schematic diagram of the calculation space division of the electromagnetic distribution simulation method for the strong electromagnetic pulse environment of the present invention;
具体实施方式Detailed ways
近年来，强电磁脉冲对电磁环境的影响成为科研研究的重要内容。本发明使用GPU加速的时域有限差分方法，实现强电磁脉冲环境下的目标物体电磁分布情况的仿真。In recent years, the impact of strong electromagnetic pulses on the electromagnetic environment has become an important content of scientific research. The invention uses a time-domain finite difference method accelerated by GPU to realize the simulation of the electromagnetic distribution of the target object under the strong electromagnetic pulse environment.
下面结合附图对本发明的具体结构作进一步的描述。The specific structure of the present invention will be further described below in conjunction with the accompanying drawings.
图1是这个方法的流程图。首先使用双曲余弦函数，实现电磁脉冲时域信号的产生。然后使用时域有限差分方法，对空间进行计算区域划分，通过GPU加速，实现电磁环境的仿真，得到时域的电磁分布情况。最后，使用快速傅里叶变换，从时域信息得到频域的电磁分布结果。Figure 1 is a flowchart of this method. Firstly, the hyperbolic cosine function is used to realize the generation of electromagnetic pulse time domain signal. Then use the time domain finite difference method to divide the calculation area of the space, and realize the simulation of the electromagnetic environment through GPU acceleration, and obtain the electromagnetic distribution in the time domain. Finally, using fast Fourier transform, the electromagnetic distribution results in the frequency domain are obtained from the time domain information.
图2是时域有限差分方法的计算空间划分示意图。(1)为电磁脉冲产生的位置，可根据需要自行设置；(2)为目标物体与环境物体等；(3)为仿真所需的CPML吸收边界，(4)为有效的计算区域。Fig. 2 is a schematic diagram of the calculation space division of the finite difference time domain method. (1) is the position where the electromagnetic pulse is generated, which can be set according to the needs; (2) is the target object and the environmental object; (3) is the CPML absorption boundary required for simulation, and (4) is the effective calculation area.
Claims (8)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510329815.9A CN104915501A (en) | 2015-06-15 | 2015-06-15 | Method for simulating electromagnetic distribution of strong electromagnetic pulse environment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510329815.9A CN104915501A (en) | 2015-06-15 | 2015-06-15 | Method for simulating electromagnetic distribution of strong electromagnetic pulse environment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN104915501A true CN104915501A (en) | 2015-09-16 |
Family
ID=54084564
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510329815.9A Pending CN104915501A (en) | 2015-06-15 | 2015-06-15 | Method for simulating electromagnetic distribution of strong electromagnetic pulse environment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104915501A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8021194B2 (en) * | 2005-04-25 | 2011-09-20 | Nvidia Corporation | Controlled impedance display adapter |
| CN104034976A (en) * | 2014-05-23 | 2014-09-10 | 国家电网公司 | Method for electromagnetic pulse response detection of single overhead line comprising nonlinear load |
| CN104573240A (en) * | 2015-01-12 | 2015-04-29 | 西安电子科技大学 | Seven point frequency domain finite difference method for analyzing periodic inhomogeneous dielectric waveguide characteristic modes |
-
2015
- 2015-06-15 CN CN201510329815.9A patent/CN104915501A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8021194B2 (en) * | 2005-04-25 | 2011-09-20 | Nvidia Corporation | Controlled impedance display adapter |
| CN104034976A (en) * | 2014-05-23 | 2014-09-10 | 国家电网公司 | Method for electromagnetic pulse response detection of single overhead line comprising nonlinear load |
| CN104573240A (en) * | 2015-01-12 | 2015-04-29 | 西安电子科技大学 | Seven point frequency domain finite difference method for analyzing periodic inhomogeneous dielectric waveguide characteristic modes |
Non-Patent Citations (1)
| Title |
|---|
| 毛玉蓉: "时间域电磁响应三维正演计算及GPU实现", 《中国博士学位论文全文数据库 基础科学辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wu | Inflow turbulence generation methods | |
| CN102121847B (en) | Method for reestablishing transient sound field | |
| MX356326B (en) | Two stage seismic velocity model generation. | |
| CN105653747B (en) | A Simulation Method for the Analysis of Electromagnetic Scattering Characteristics of Conformal Subgrid for Ultra-High-Speed Vehicles | |
| CN106908835A (en) | Band limit Green's function filters multiple dimensioned full waveform inversion method | |
| CN103593510A (en) | Rough surface and target composite electromagnetic scattering simulation method based on reciprocity principle | |
| CN103279600A (en) | Medium rough surface finite element electromagnetic simulation method based on integral boundary | |
| CN104573289A (en) | Electromagnetic scattering simulation method of one-dimensional conductor rough sea surface and two-dimensional conductor target | |
| CN117332658A (en) | Method and system for determining electromagnetic characteristics of anisotropic time-varying plasma | |
| CN104809286A (en) | Method for realizing perfectly matched absorbing boundary of expanding coordinates in plasma | |
| CN104915465A (en) | Metal target transient electromagnetic scattering analysis method based on delay laguerre polynomials | |
| CN105844019A (en) | Ground point selection method based on vehicle ground current distribution regularities | |
| CN104915501A (en) | Method for simulating electromagnetic distribution of strong electromagnetic pulse environment | |
| CN104951606A (en) | Method for strong electromagnetic pulse environment simulation based on frequency-time transformation | |
| Zheng et al. | Accelerating Computation of Electromagnetic Properties for Low-Frequency Wave Propagating Across the Sea–Air Interface | |
| Chalikov et al. | Three-dimensional periodic fully nonlinear potential waves | |
| Wang et al. | Computation of rotor noise generation in grid turbulence using large-eddy simulation | |
| CN111191392B (en) | A Fast Method for Computing Electromagnetic Scattering Problems of Electrically Large Targets | |
| Hou et al. | A complex equivalent source method for scattering effect of aircraft noise | |
| Yu et al. | Quadrupole noise in turbulent wake interaction problems | |
| Duan et al. | Application of MW-FDTD to simulate the EM wave propagation over ocean with OpenCL | |
| CN104239642B (en) | Magnetosonic couples the vector method for solving of direct problem under a kind of sinusoidal excitation | |
| CN111339673B (en) | Multi-scale noise simulation analysis method | |
| Wang et al. | A One-Dimensional Electronic Fluid Model Based on the Finite-Difference Time-Domain Method | |
| Liu et al. | Sound Propagation Considering Atmospheric Inhomogeneity and Ground Effects |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150916 |