CN106294283B - 基于泰勒级数展开的时域积分方程快速方法 - Google Patents

基于泰勒级数展开的时域积分方程快速方法 Download PDF

Info

Publication number
CN106294283B
CN106294283B CN201510266484.9A CN201510266484A CN106294283B CN 106294283 B CN106294283 B CN 106294283B CN 201510266484 A CN201510266484 A CN 201510266484A CN 106294283 B CN106294283 B CN 106294283B
Authority
CN
China
Prior art keywords
basic function
field
taylor series
series expansion
time
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.)
Active
Application number
CN201510266484.9A
Other languages
English (en)
Other versions
CN106294283A (zh
Inventor
陈如山
李威
丁大志
樊振宏
程光尚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201510266484.9A priority Critical patent/CN106294283B/zh
Publication of CN106294283A publication Critical patent/CN106294283A/zh
Application granted granted Critical
Publication of CN106294283B publication Critical patent/CN106294283B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于泰勒级数展开的时域积分方程快速算法。传统的时域积分方程方法由于计算时间和内存消耗的限制无法解决大规模的电磁问题。而发展快速算法是解决实际工程问题的必由之路。本发明是利用自由空间中格林函数的泰勒级数展开重构成聚合、转移、配置的形式进而实现矩阵矢量乘的加速计算,本发明对于求解金属复杂目标散射问题需要更少的计算内存以及计算时间。而且编程相对简单易于实现,具有很强的实际工程应用价值。

Description

基于泰勒级数展开的时域积分方程快速方法
技术领域
本发明涉及金属目标瞬态电磁散射特性数值计算技术,属于目标电磁散射特性的快速计算技术领域,具体是一种基于泰勒级数展开的时域积分方程快速算法。
背景技术
随着计算电磁学领域研究的深入,传统的频域方法已经不能满足需要。以计算机硬件技术的发展为契机,人们逐步具有了直接在时域对具有宽频带特性的瞬变电磁场的计算分析能力,从而实现了对物理量和物理现象更深刻、更直观的理解。同频域方法相比,在时域求解目标的电磁特性不仅可以直观的揭示目标与电磁波相互作用的机理,而且通过少量的计算就可以获得目标的宽频带信息,这在宽带电磁问题、瞬态电磁问题分析中具有明显的优势。
同基于微分方程的时域方法(FDTD、FETD等)相比,基于积分方程的时域方法在求解的未知数数量上具有明显的优势,这是因为积分方程利用格林函数建立源和场的关系,求解区域在边界上,离散后未知数的数量与边界面积成正比。其次,积分方程方法自动满足辐射边界条件,不需要强加吸收边界,而吸收边界是基于微分方程的方法所必需的。
由于计算机水平和硬件的限制,传统的时域积分方程方法无法解决大规模的电磁问题。而发展快速算法是解决实际工程问题的必由之路。时域积分方程快速算法的研究是在近20年开展起来的,其中最著名的两种:一是时域平面波算法 (PWTD),二是时域自适应积分方程(TD-AIM)。对于三维目标模型,TD-AIM的计算量和存储量分别为可见其处理三维目标的效率低下。PWTD算法其原理类似于FMM算法,但是因为时间变量的存在,使得 PWTD算法更加复杂,实现起来更加困难。并且该算法的适应性差,当用以求解介质或者有耗媒质等复杂问题时,还需要对算法进行特殊修改。
发明内容
本发明的目的在于提供一种高效、稳定的分析金属目标的时域电磁散射特性的快速方法。由于对远场部分采用泰勒级数展开重构成聚合、转移、配置的形式进而实现矩阵矢量乘的加速计算,本发明对于求解金属复杂目标散射问题需要更少的计算内存以及计算时间。而且编程相对简单易于实现。
实现本发明目的的技术解决方案为:一种基于泰勒级数展开的时域积分方程快速算法,步骤如下:
第一步,建立时域电磁场积分方程。金属目标在入射电磁波的作用下,在金属散射体的表面上会产生面感应面电流,感应电流的不断变化,会向外辐射出电磁场,由感应电流产生的电磁场就是散射场,利用入射电磁场和散射电磁场在金属表面满足的边界条件建立时域电磁场积分方程。
第二步,将散射体表面上离散得到的子散射体分组。任意两个子散射体间的互耦或自耦根据它们所在组的位置关系而分成近场组对和远场组对。当它们是近场组对时,采用直接数值计算。而当它们为远场组对时,则采用泰勒级数展开成聚合-转移-配置方法计算。
第三步,近场阻抗矩阵计算,将金属表面电流密度用空间基函数和时间基函数展开,并在空间域上进行伽辽金测试,时间域上进行点匹配得到矩阵元素值。
第四步,远场部分采用泰勒级数展开重构成聚合、转移、配置的形式。
第五步,矩阵方程求解以及电磁散射参数的计算。利用时间递推的方式求解每个时刻的电流系数。每个时刻空间某处的散射场贡献来源分为近场组源的贡献和远场组源的贡献。前者直接通过近场矩阵和电流系数相乘得到,后者通过泰勒级数展开成聚合-转移-配置方法快速计算。采用迭代法求解出最终的散射电流系数。并根据互易定理求解出雷达散射截面积。
本发明与现有技术相比,其显著优点为:(1)可以快速分析复杂电大目标时域电磁散射特性;(2)远场计算是基于泰勒级数展开的,避免了PWTD算法需要的复杂的时频域转换、谱域积分等操作,编程相对简单易于实现;(3)适用性高,可以只需作较少的改动便可应用于求解介质、有耗、色散等复杂问题的分析中。
附图说明
图1是本发明的应用于基于泰勒级数展开的时域积分方程快速算法中的多层分组示意图。
图2是本发明的应用于基于泰勒级数展开的时域积分方程快速算法中的近远场关系示意图。
图3是本发明的应用于基于泰勒级数展开的时域积分方程快速算法中的场源基函数位置向量分布示意图。
图4是本发明的应用于基于泰勒级数展开的时域积分方程快速算法中的简易导弹模型不同频点处双站RCS比较图。(a)30MHz;(b)150MHz;(c)270MHz。
具体实施方式
下面结合附图对本发明作进一步详细描述。
本发明一种基于泰勒级数展开的时域积分方程快速算法,步骤如下:
第一步,建立时域电磁场积分方程。
金属目标在入射电磁波{Einc(r,t),Hinc(r,t)}的作用下,在金属散射体的表面上会产生面感应面电流J(r,t),感应电流的不断变化,会向外辐射出电磁场,由感应电流产生的电磁场就是散射场{Esca(r,t),Hsca(r,t)},利用入射电磁场和散射电磁场在金属表面满足的边界条件建立时域电磁场积分方程。
其中是单位外法向矢量,α取值为0~1之间的实数,μ和ε分别是自由空间的磁导率和介电参数。η是自由空间的波阻抗。R=|r-r′|,c是自由空间中的光速,τ=t-R/c是延时。
第二步,将散射体表面上离散得到的子散射体分组。
我们借鉴MLFMM里面对空间基函数进行分组的思想,将基函数的相互作用转换为组的相互作用。假设空间中有一个刚好可以把整个目标物体包围的立方体,把这个立方体等分成8个子立方体,接着在将每个子立方体等分成8个更小的立方体,如图1(a-d)所示。依次类推,直到达到预先设置的门限值,停止划分。定义第一层为划分最细的一层,第一层中的立方体中如果包含基函数,定义为第一层组,对于其它层组,依次类推即可。
任意两个子散射体间的互耦或自耦根据它们所在组的位置关系而分成近场组对和远场组对。对于每一层,我们都会设定一个参数β,4<β<6。首先从最高层Nl的各组着手划分,如果两组中心间的距离大于规定的的门限值β,就把其称之为远场组对;然后在其子层中,组对中心间距大于规定的该层的门限值,同时这两个组都没有划分到Nl层的远场组对,那么就规定这样的组对为该层的远场组对;根据以上所述依次类推,对于每一层我们就可以得到其远场组对的信息,如图2所示。当它们是近场组对时,采用直接数值计算。而当它们为远场组对时,则采用泰勒级数展开成聚合-转移-配置方法计算。
第三步,近场阻抗矩阵计算。
把金属表面电流电流密度J(r,t)用RWG基函数Λn(r)作为空间基函数,三角基函数Tj(t)作时间基函数展开:
其中Ns是散射体包含的RWG基函数的个数,Nt是时间基函数的个数,In,j是第 n个RWG基函数第j个时间步上时间基函数的系数。将式(4)代入(1),在空间上进行Galerkin测试,在时间上进行点匹配,最终得到如下的矩阵方程:
第四步,远场部分采用泰勒级数展开重构成聚合、转移、配置的形式。
现在通过计算源点r′处的源信号Jn(r′,t)对场点r处的辐射贡献来描述该算法的基本原理。假设源点r′所在的空间基函数为Λn(r′),因此r′处的源信号 Jn(r′,t)可以展开如下:
实际操作中我们需要将源信号分解为持续时间更短的分段子信号,保证每段子信号的持续时间满足约束条件。源信号Jn(r′,t)可以被分解为L段连续的子信号 Jn,l(r′,t),每一段子信号的持续时间为Ts=(Mt+1)Δt。源信号可以写成如下形式:
那么源点r′处第l段子信号在场点r处产生的测试场为:
其中
现假设源点rn与场点rm分别位于两个组内,如图3所示,两个组分别称为源组和场组,组中心分别为ri和rj,场源基函数之间的矢量可以表示为:
R=rmi+rij-rnj=Rm-Rn (13)
这里,rij=ri-rj,rmi=rm-ri,rnj=rn-rj Rm=rmi+rij/2,Rn=rnj-rij/2
利用泰勒级数展开可以得到如下表达式:
则式(11)(12)可以写成聚合、转移、配置的形式。
第五步,矩阵方程求解以及电磁散射参数的计算。
利用时间递推的方式求解每个时刻的电流系数。现将任意空间上的散射场按来源分为两部分:一部分由该场点所在组的近场组NFP(α)中的源产生的;另一部分由该场点所在组的远场组FFP(α)中的源产生。基于泰勒级数展开的时域积分方程快速算法的递推公式变为:
近场区域产生的贡献是由经典MOT算法计算得到,主要通过矩阵元素值与电流系数相乘得到。在远场区域,也就是对应组中的远场组对,这些组之间的相互作用,过泰勒级数展开成聚合-转移-配置方法快速计算。
由于时刻iΔt以前的Ij在时刻iΔt都是已知的,j=1,2,3,...,i-1,这样每个时间步求解一次式(17)的矩阵方程,就可以得到每个时刻iΔt上的Ii所以可以递推求出各时刻的电流值。最后可以根据求得的瞬态电流系数计算出我们需要的电磁散射参数。
为了验证本发明的正确性与有效性,下面分析了简易导弹模型的的电磁散射特性。
算例:简易导弹模型,几何尺寸为长8m,半径0.25m。激励源设置为:调制高斯脉冲,中心频率150MHz,带宽300MHz,入射波的方向θ=180°,(弹头方向照射),极化方式VV极化,观察角度0°≤θ≤180°,导弹模型采用0.1λ(λ=1m)剖分得到8034个三角形,总未知量为12051。分组大小为0.3λ。 3个盒子之外为远场区域,近场门限0.3λ*3。使用Δt=333.3ps,一共计算了600 个时间步。CFIE的参数α=0.5,采用GMRES迭代方法计算,收敛精度为1e-9。图4给出了导弹模型在30MHz、150MHz以及270MHz时的双站RCS曲线并与经典MOT算法计算结果进行对比,表格1统计了本专利提出的方法与传统的时域积分方法消耗的时间和内存对比。
表1.计算时间和内存消耗对比
从图4中可以看出本发明方法计算的双站RCS数据和MOT计算的结果十分吻合,证明了本发明方法的精确性。从表1看出,采用本发明方法相比经典 MOT算法需要更好的计算内存和计算时间。进一步证明了本发明方法的有效性。
本发明的一种基于泰勒级数展开的时域积分方程快速算法,有益效果如下:实现过程简单,相较于需要的复杂的时频域转换、谱域积分等操作,本发明方法只需在八叉树分组以及多层近远场划分的基础上,对远场组对之间的距离R进行泰勒展开重构成聚合、转移、投射的操作,编程相对简单易于实现。而且本发明方法适用性高,可以只需作较少的改动便可应用于求解介质、有耗、色散等复杂问题的分析中。

Claims (1)

1.一种基于泰勒级数展开的时域积分方程快速方法,其特征在于步骤如下:
第一步,建立时域电磁场积分方程;利用入射电磁场和散射电磁场在金属表面满足的边界条件建立时域电磁场积分方程;
第二步,将散射体表面上离散得到的子散射体分组;用一个立方体盒子将整个目标物体包围,把这个立方体等分成8个子立方体,接着将每个子立方体等分成8个更小的立方体,依次类推,直到达到预先设置的门限值,停止划分;任意两个子散射体间的互耦或自耦根据它们的位置关系而分成近场作用对和远场作用对;当它们是近场作用对时,计算近场阻抗矩阵,当它们为远场作用对时,采用泰勒级数展开成聚合-转移-配置方法计算源组基函数在场组基函数处产生的场;
第三步,计算近场阻抗矩阵,将金属表面电流密度用空间基函数和时间基函数展开,并在空间域上进行伽辽金测试,时间域上进行点匹配得到矩阵元素值;
第四步,远场作用对之间采用泰勒级数展开重构成聚合、转移、配置的形式计算源组基函数在场组基函数处产生的场;
第五步,矩阵方程求解以及电磁散射参数的计算;利用时间递推的方式求解每个时刻的电流系数,采用迭代法求解出最终的感应电流系数,最后根据求得的瞬态电流系数计算出需要的电磁散射参数;
所述步骤四中,在八叉树分组以及多层近远场划分的基础上,对远场作用对之间进行泰勒级数展开重构成聚合、转移、投射的操作来加速矩阵矢量乘;
设源点r'所在的空间基函数为Λn(r'),r'处的源信号Jn(r',t)展开如下:
其中空间基函数Λn(r)为RWG基函数,时间基函数Tj(t)为三角基函数,Nt是时间基函数的个数;
源信号Jn(r',t)被分解为L段连续的子信号Jn,l(r',t),每一段子信号的持续时间为Ts=(Mt+1)Δt,Mt为每段子信号的长度,源信号写成如下形式:
源点r'处第l段子信号在场点r处产生的测试电磁场为:
其中Λm(r)为场点r处的测试基函数,nm为场点r处的单位外法向矢量,En,l(r,t),Hn,l(r,t)为源点r'处第l段子信号在场点r处产生的电磁和磁场,μ和ε分别是自由空间的磁导率和介电常数;R=|r-r'|,c是自由空间中的光速,τ=t-R/c是延时;lm,ln分别为第m和n条边的边长, 为基函数在其支撑域内的积分;
源点rn与场点rm分别位于两个组内,两个组分别称为源组和场组,组中心分别为ri和rj,场源基函数之间的矢量表示为:
R=rmi+rij-rnj=Rm-Rn (5)
这里,rij=ri-rj,rmi=rm-ri,rnj=rn-rj,Rm=rmi+rij/2,Rn=rnj-rij/2
当任意两个子散射体间的互耦或自耦为远场作用对时,利用泰勒级数展开得到如下表达式:
则将(5-8)代入式(3)(4)将写成聚合、转移、配置的形式来加速矩阵矢量乘。
CN201510266484.9A 2015-05-22 2015-05-22 基于泰勒级数展开的时域积分方程快速方法 Active CN106294283B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510266484.9A CN106294283B (zh) 2015-05-22 2015-05-22 基于泰勒级数展开的时域积分方程快速方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510266484.9A CN106294283B (zh) 2015-05-22 2015-05-22 基于泰勒级数展开的时域积分方程快速方法

Publications (2)

Publication Number Publication Date
CN106294283A CN106294283A (zh) 2017-01-04
CN106294283B true CN106294283B (zh) 2019-05-07

Family

ID=57633456

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510266484.9A Active CN106294283B (zh) 2015-05-22 2015-05-22 基于泰勒级数展开的时域积分方程快速方法

Country Status (1)

Country Link
CN (1) CN106294283B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107679287B (zh) * 2017-09-11 2021-07-13 三峡大学 基于3步4阶隐式泰勒级数法的电磁暂态数值计算方法
CN111159637B (zh) * 2020-01-02 2023-06-20 西北工业大学 一种应用于磁化等离子体计算的电磁波时域精细积分方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Efficient Generation of Method of Moments Matrices Using Equivalent Dipole-Moment Method;Jiade Yuan等;《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》;20091231;第8卷;第716-719页
Efficient Generation of Method of Moments Matrices Using the Characteristic Function Method;Junho Yeo等;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20041231;第52卷(第12期);第3405-3410页
Fast Dipole Method for Electromagnetic Scattering From Perfect Electric Conducting Targets;Xinlei Chen等;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20120229;第60卷(第2期);第1186-1191页:摘要,正文第III小节
应用EDM法和ACA算法快速计算电大开缝导体RCS;姚雨帆等;《合肥工业大学学报(自然科学版)》;20140430;第37卷(第4期);第440-443页
时域平面波算法研究;李金艳;《中国优秀硕士学位论文全文数据库-信息科技辑》;20070915(第03期);第I135-5页:摘要,正文第3.1小节、4.1小节
时间步进时域体面积分方程及其快速方法的研究;秦宇晨;《中国优秀硕士学位论文全文数据库-信息科技辑》;20140715(第07期);I135-15页:摘要,正文第3.2小节

Also Published As

Publication number Publication date
CN106294283A (zh) 2017-01-04

Similar Documents

Publication Publication Date Title
Jandhyala et al. A combined steepest descent-fast multipole algorithm for the fast analysis of three-dimensional scattering by rough surfaces
Guan et al. Electromagnetic scattering from objects above a rough surface using the method of moments with half-space Green's function
CN105930567A (zh) 一种基于子区域自适应积分的电磁散射特性获取方法
CN106096267B (zh) 一种腔体电磁散射特性快速计算方法
CN104573289A (zh) 一维导体粗糙海面与二维导体目标的电磁散射仿真方法
CN103235888A (zh) 一种精确计算双各向异性介质球电磁散射的方法
CN104317984A (zh) 基于分域建模的船舶电磁散射预测方法及系统
CN106446470A (zh) 一种高效并行的非均匀介质频域有限差分方法
CN106294283B (zh) 基于泰勒级数展开的时域积分方程快速方法
CN109765538A (zh) 非均匀介质目标体的散射场确定方法
Cao et al. A hybrid approach using machine learning and genetic algorithm to inverse modeling for single sphere scattering in a Gaussian light sheet
CN110274920A (zh) 体面剖分弹跳射线法分析金属介质目标瞬态特性的方法
Su et al. Domain decomposition scheme with equivalence spheres for the analysis of aircraft arrays in a large-scale range
Niu et al. ARCS: Active radar cross section for multi-radiator problems in complex EM environments
Liao et al. 2D inversion of magnetotelluric data using deep learning technology
Ghorbani et al. Through-the-wall radar under electromagnetic complex wall: A deep learning approach
Ergül Fast and accurate solutions of electromagnetics problems involving lossy dielectric objects with the multilevel fast multipole algorithm
CN105184033B (zh) 基于阶数步进金属目标的宽频带电磁特性快速预估方法
Zhao et al. Modifications on parametric models for distributed scattering centres on surfaces with arbitrary shapes
Tsuji et al. A fast directional algorithm for high-frequency electromagnetic scattering
CN111199097B (zh) 瞬态矢量辐射传输理论箔条云散射处理方法
CN105224780A (zh) 分析导体瞬态电磁散射特性的时域高阶Nystrom方法
Mansukhani et al. Rcs based target classification using deep learning methods
Cheng et al. Comparison of particle swarm optimization and self-adaptive dynamic differential evolution for the imaging of a periodic conductor
CN105808804B (zh) 分析金属介质混合目标瞬态电磁散射特性的时域Nystr*m方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant