发明名称 超高速飞行目标的电磁散射分析方法
摘要 本发明公开了一种超高速飞行目标的电磁散射分析方法。针对包裹在超高速飞行目标周围等离子体的非均匀特性,采用了体面积分方程方法进行分析,对等离子体等效相对介电常数接近于1的部分当作空气进行处理,等效相对介电常数较大的区域采用自适应网格加密处理,以便达到所需的求解精度。相对于采用统一网格体剖分等离子壳套的传统方法,本发明中的方法可以大大的节省计算资源,同时由于体面积分方程中采用的格林函数为真空中的格林函数,多层快速多级子技术被用来进一步加速求解,使得本发明对于求解超高速飞行目标散射问题需要更少的计算内存以及计算时间。
申请公布号 CN103198227A 申请公布日期 2013.07.10
申请号 CN201310136199.6 申请日期 2013.04.18
申请人 南京理工大学 发明人 陈如山;樊振宏;丁大志;陶诗飞;盛亦军;陈峰;呼延龙;沙侃;叶晓东
分类号 G06F19/00(2006.01)I 主分类号 G06F19/00(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱显国
主权项 1.一种超高速飞行目标的电磁散射分析方法,其特征在于步骤如下:第一步,建立超高速飞行目标等离子壳套模型,根据飞行目标的飞行高度、攻角和飞行马赫数参数,对超高速飞行目标进行气动模拟计算,得到目标的电子数密度、温度和压强信息数据,由此得到等离子体特征频率以及碰撞频率,再由以下公式得到等离子壳套各空间位置的等效相对介电常数,<maths num="0001"><![CDATA[<math><mrow><msub><mi>&epsiv;</mi><mi>r</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><msub><msup><mi>&omega;</mi><mn>2</mn></msup><mi>pe</mi></msub><mrow><msup><mi>&omega;</mi><mn>2</mn></msup><mo>+</mo><msup><mi>v</mi><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mi>j</mi><mfrac><mi>v</mi><mi>&omega;</mi></mfrac><mfrac><msub><msup><mi>&omega;</mi><mn>2</mn></msup><mi>pe</mi></msub><mrow><msup><mi>&omega;</mi><mn>2</mn></msup><mo>+</mo><msup><mi>v</mi><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中ω<sub>pe</sub>为等离子体特征频率,ω为电磁波频率,v为等离子体碰撞频率;第二步,对飞行目标金属部分采用三角形剖分,等离子壳套等效相对介电常数接近于1的部分当作空气,不需要进行网格剖分,等离子壳套非空气区域部分采用四面体剖分,并对等效相对介电常数较大的区域采用自适应网格加密处理;第三步,根据包裹等离子壳套高超速飞行目标结构的散射特性,采用矩量法基础理论,得到体面积分方程,其矩阵方程形式为:<maths num="0002"><![CDATA[<math><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>Z</mi><mi>mn</mi><mi>DD</mi></msubsup></mtd><mtd><msubsup><mi>Z</mi><mi>mn</mi><mi>MD</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>Z</mi><mi>mn</mi><mi>DM</mi></msubsup></mtd><mtd><msubsup><mi>Z</mi><mi>mn</mi><mi>MM</mi></msubsup></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>D</mi><mi>n</mi></msub></mtd></mtr><mtr><mtd><msub><mi>I</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>v</mi><mi>m</mi><mi>V</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>v</mi><mi>m</mi><mi>S</mi></msubsup></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>(2)Z<sup>DD</sup>代表介质对介质的作用,Z<sup>DM</sup>表示介质对金属的作用,Z<sup>MD</sup>都表示介质和金属的相互作用部分,Z<sup>MM</sup>表示金属对金属的作用,D<sub>n</sub>和I<sub>n</sub>是待求未知系数,<img file="FDA00003070385900013.GIF" wi="59" he="64" />和<img file="FDA00003070385900014.GIF" wi="58" he="70" />是右边向量激励;第四步,对第三步中的自由空间格林函数<img file="FDA00003070385900015.GIF" wi="162" he="74" />按照加法定理展开,并给出体面积分方程中的体聚合因子、转移因子和配置因子的具体表达式;第五步,求解矩阵方程(2),得到电流系数,再根据互易定理由电流系数计算电磁散射参量。
地址 210094 江苏省南京市孝陵卫200号