CN104915501A - 一种强电磁脉冲环境电磁分布仿真方法 - Google Patents
一种强电磁脉冲环境电磁分布仿真方法 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
- time
- electromagnetic
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种仿真强电磁脉冲环境下的电磁分布情况仿真方法。该方法使用双曲余弦函数模拟时域下的电磁脉冲波形,使用时域有限差分算法对环境的电磁场变化进行仿真模拟,使用GPU实现时域有限差分方法的加速,得到时域上的电磁分布情况,最后通过快速傅里叶变换,得到频率域上的电磁分布情况,解决了电磁脉冲武器攻击下超电大尺寸对象的电磁环境分布特性计算问题。
Description
技术领域
本发明涉及计算机仿真中的强电磁脉冲模拟以及并行时域有限差分电磁场仿真技术,解决了电磁脉冲武器攻击下超电大尺寸对象的电磁环境分布特性计算问题。
背景技术
电磁脉冲具有宽频谱和高能量特性,能通过线缆、缝隙等耦合到电子系统中,导致电子系统性能降级甚至损坏。因此,研究复杂电子系统所处的电磁脉冲环境,对提高复杂电子系统对电磁脉冲的防护能力具有重要的意义。
时域有限差分方法是把Maxwell方程式在时间和空间领域上进行差分化。利用蛙跳式方式在空间领域内的电场和磁场进行交替计算,通过时间领域上更新来模仿电磁场的变化,达到数值计算的目的。其优点是能够直接模拟场的分布,精度比较高,是目前使用比较多的数值模拟的方法之一。
考虑到时域有限差分本身算法的限制,仿真范围越大、电磁波频率越高,需要的计算量以及内存使用量越大。而近些年出现的GPU加速方式,可以充分发掘时域有限差分本身的并行特质,使计算速度获得成百倍的提升。
发明内容
本发明使用时域有限差分算法,实现强电磁脉冲环境下的电磁情况仿真,并通过GPU加速仿真过程,实现超电大尺寸目标的电磁分布仿真。
本发明所采用的技术方案:
一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用时域有限差分方法作为电磁场仿真算法进行仿真,建立目标区域时频空多维分析模型。
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用GPU实现时域有限差分方法的并行加速。
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:包含电磁脉冲生成函数、基于GPU的时域有限差分方法、快速傅里叶变换三部分。
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用双曲余弦函数,模拟核电磁脉冲时域波形,作为仿真的信号激励。
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用时域有限差分方法,进行空间中的电磁波模拟,得到空间中电磁场的分布的时域情况。
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用GPU实现对时域有限差分算法的并行加速,减少仿真超电大尺寸目标及高频脉冲时的时间。
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:,使用快速傅里叶变换,从电磁场分布的时域情况得到空间电磁场分布的频域情况。
根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:时域有限差分方法所使用的吸收边界为卷积完美匹配层(CPML)。
本发明相对于现有技术具有如下的优点及特性:
(1)从Maxwell方程直接推导出的时域有限差分算法,仿真精度高。
(2)使用GPU进行时域有限差分方法的并行实现,仿真速度快。
(4)使用快速傅里叶变换,可以从时域结果中直接得到频域的结果。
附图说明
图1是本发明所述强电磁脉冲环境电磁分布仿真方法流程图;
图2是本发明所述强电磁脉冲环境电磁分布仿真方法计算空间划分示意图;
具体实施方式
近年来,强电磁脉冲对电磁环境的影响成为科研研究的重要内容。本发明使用GPU加速的时域有限差分方法,实现强电磁脉冲环境下的目标物体电磁分布情况的仿真。
下面结合附图对本发明的具体结构作进一步的描述。
图1是这个方法的流程图。首先使用双曲余弦函数,实现电磁脉冲时域信号的产生。然后使用时域有限差分方法,对空间进行计算区域划分,通过GPU加速,实现电磁环境的仿真,得到时域的电磁分布情况。最后,使用快速傅里叶变换,从时域信息得到频域的电磁分布结果。
图2是时域有限差分方法的计算空间划分示意图。(1)为电磁脉冲产生的位置,可根据需要自行设置;(2)为目标物体与环境物体等;(3)为仿真所需的CPML吸收边界,(4)为有效的计算区域。
Claims (8)
1.一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用时域有限差分方法作为电磁场仿真算法进行仿真,建立目标区域时频空多维分析模型。
2.根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用GPU实现时域有限差分方法的并行加速。
3.根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:包含电磁脉冲生成函数、基于GPU的时域有限差分方法、快速傅里叶变换三部分。
4.根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用双曲余弦函数,模拟核电磁脉冲时域波形,作为仿真的信号激励。
5.根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用时域有限差分方法,进行空间中的电磁波模拟,得到空间中电磁场的分布的时域情况。
6.根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:使用GPU实现对时域有限差分算法的并行加速,减少仿真电大尺寸目标及高频脉冲时的时间。
7.根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:,使用快速傅里叶变换,从电磁场分布的时域情况得到空间电磁场分布的频域情况。
8.根据权利要求1所述的一种强电磁脉冲环境电磁分布仿真方法,其特征是:时域有限差分方法所使用的吸收边界为卷积完美匹配层(CPML)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510329815.9A CN104915501A (zh) | 2015-06-15 | 2015-06-15 | 一种强电磁脉冲环境电磁分布仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510329815.9A CN104915501A (zh) | 2015-06-15 | 2015-06-15 | 一种强电磁脉冲环境电磁分布仿真方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104915501A true CN104915501A (zh) | 2015-09-16 |
Family
ID=54084564
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510329815.9A Pending CN104915501A (zh) | 2015-06-15 | 2015-06-15 | 一种强电磁脉冲环境电磁分布仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104915501A (zh) |
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 (zh) * | 2014-05-23 | 2014-09-10 | 国家电网公司 | 包含非线性负载的单根架空线路电磁脉冲响应检测方法 |
CN104573240A (zh) * | 2015-01-12 | 2015-04-29 | 西安电子科技大学 | 周期性非均匀介质波导特征模分析的七点频域有限差分方法 |
-
2015
- 2015-06-15 CN CN201510329815.9A patent/CN104915501A/zh 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 (zh) * | 2014-05-23 | 2014-09-10 | 国家电网公司 | 包含非线性负载的单根架空线路电磁脉冲响应检测方法 |
CN104573240A (zh) * | 2015-01-12 | 2015-04-29 | 西安电子科技大学 | 周期性非均匀介质波导特征模分析的七点频域有限差分方法 |
Non-Patent Citations (1)
Title |
---|
毛玉蓉: "时间域电磁响应三维正演计算及GPU实现", 《中国博士学位论文全文数据库 基础科学辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sumithra et al. | Review on computational electromagnetics | |
Liu et al. | Impact of deck and jet blast deflector on the flow and acoustic properties of an imperfectly expanded supersonic jet | |
CN105653747B (zh) | 超高速飞行器共形亚网格电磁散射特性分析的仿真方法 | |
CN113933905A (zh) | 一种圆锥型场源瞬变电磁反演方法 | |
CN106446470A (zh) | 一种高效并行的非均匀介质频域有限差分方法 | |
Izadi et al. | An analytical second-FDTD method for evaluation of electric and magnetic fields at intermediate distances from lightning channel | |
CN106055837A (zh) | 外场激励下有损大地上架空线缆电路模型建立方法和系统 | |
CN104915501A (zh) | 一种强电磁脉冲环境电磁分布仿真方法 | |
CN102332055B (zh) | 一种极低频电磁波的仿真计算方法 | |
CN114925539B (zh) | 声爆波反演方法、装置、设备及存储介质 | |
CN104951606A (zh) | 一种基于频-时变换的强电磁脉冲环境仿真方法 | |
Oikawa et al. | Analysis of lightning electromagnetic field on large‐scale terrain model using three‐dimensional MW‐FDTD parallel computation | |
CN103398832A (zh) | 一种电液加速度伺服系统波形再现控制方法 | |
CN105760596A (zh) | 基于辅助微分方程的二维真空Crank-Nicolson完全匹配层实现算法 | |
Pingenot et al. | Full wave analysis of RF signal attenuation in a lossy rough surface cave using a high order time domain vector finite element method | |
CN105589678A (zh) | 一种用数字信号处理技术实现的时域有限差分方法 | |
Xie et al. | Simulation and analysis of radiated electromagnetic environment from cable in cabin | |
Hou et al. | A complex equivalent source method for scattering effect of aircraft noise | |
Lei et al. | FDTD method for three-dimensional numerical modeling of ESD radiation field | |
Hong-xia et al. | Numerical simulation and experiments on muzzle blast overpressure in large-caliber weapons | |
Lei et al. | High-efficiency Numerical Simulation of Electromagnetic Response of Complex Geoelectric Structures Based on GPU Parallel Symplectic Algorithm | |
Yu et al. | A novel mesh generation method for FDTD without ray-tracing | |
Malinga et al. | Modeling lightning strike behavior in the near field of elevated systems | |
Wang et al. | Simulation of Near-Surface Propagation of Powerful Electromagnetic Pulse Based on Spherical coordinates | |
Liu et al. | High-order dgtd for solving em scattering from hypersonic aircraft with plasma sheath |
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 |
Application publication date: 20150916 |
|
WD01 | Invention patent application deemed withdrawn after publication |