许 勇,黄 勇,余永刚

(中国空气动力研究与发展中心 计算空气动力学研究所,四川 绵阳 621000)

电大尺寸目标电磁散射的并行FVTD计算

许 勇*,黄 勇,余永刚

(中国空气动力研究与发展中心 计算空气动力学研究所,四川 绵阳 621000)

实际目标的高频电磁散射模拟属于电大尺寸目标电磁计算问题，其高精度数值计算通常伴随着大存储量和大计算量的沉重负担。为保证高精度数值计算，构建了直接求解电磁学麦克斯韦方程组，采用贴体多块结构网格的时域有限体积法(FVTD)解算器。为解决大规模网格带来的大计算量问题，则采用MPI并行编程，进行网格多进程分割、负载平衡以及进程通信设置和程序并行化处理，成功地对电大尺寸飞翼外形进行了L波段双站电磁散射场和雷达截面(RCS)计算。结果表明pmbRCS3d这一并行高精度电磁散射模拟软件具有稳定和鲁棒特性，适合进一步应用于目标更高频段电磁散射计算。

麦克斯韦方程组；时域有限体积法(FVTD)；并行计算；L波段；雷达截面(RCS)

0 引 言

飞行器气动/隐身优化设计需要发展快速高效、高精度电磁数值方法，但是复杂外形目标的电磁散射包含了镜面散射、边缘绕射、爬行波等多种复杂现象，准确模拟非常困难。

高频渐进等传统工程方法虽然快捷，能捕获主要散射特征，但误差较大，不能满足隐身设计的精细要求。随着计算机能力的提高，国内外高精度数值方法获得巨大发展，主要有计算微分方程的时域有限差分FDTD[1-3]方法，和计算积分方程的多层快多极子方法[4-7]。其中直接求解电磁学麦克斯韦方程组的时域微分方法，最早来源于Yee K S[1]提出的FDTD，它基于蛙跳格式，电场和磁场在网格交错配置以弥补格式黏性，有易编程和整个散射空间可统一计算的优点，但由于使用笛卡儿直角网格，模拟物面存在阶梯效应。20世纪90年代以来，基于流体力学的欧拉方程组和电磁学的麦克斯韦方程组具有共同的双曲型数学特征，人们开始在计算电磁学(CEM)引入计算流体力学(CFD)技术[8-15]，比较典型的是采用贴体坐标系的时域有限体积法(FVTD)。

FVTD解算器有以下特点：直接求解麦克斯韦方程组；在宽频入射波的一个计算过程中可得到多个频率的结果；能达到接近矩量法(MOM)的数值计算精度；存储仅与未知量数目同数量级，能不需特殊处理地模拟散射、腔激励等复杂现象。但因涉及波在空间长时间传播计算，以及抑制耗散和色散要有一定的网格密度，网格规模会随频率增加。

实际飞机等目标在高频区为电大尺寸，例如B2隐形轰炸机，长度和翼展达数十米，对应高频区(X波段)波长最大到1000量级，每个波长取10～20个网格点，则整个规模高达10亿量级。超大计算规模成为制约气动/隐身优化的瓶颈。目前解决途径是在高性能计算机集群上，采用大规模并行计算技术。

为解决电大目标高频电磁散射大规模计算，本文通过空间网格块分割、负载平衡、多进程MPI通信等，建立了多块结构网格FVTD并行程序pmbRCS3d，在与金属球RCS解析解验证对比基础上，计算了电大尺寸飞翼外形L波段双站电磁散射场和雷达截面(RCS)，结果表明了并行计算程序的正确性和鲁棒性，利于向更高频段(X波段)扩展。

1 数值方法

1.1 控制方程

麦克斯韦方程组的两个时变旋度微分方程，法拉第电磁感应定律和安培环路定理，在自由空间、直角坐标系下可表示为守恒形式：

式中:B是磁感应强度矢量，H是磁场强度，D是电位移矢量，E是电场强度矢量。

对于复杂外形物体，经坐标变换，

得到曲线坐标系下麦克斯韦方程组守恒形：

式中，V是坐标变换的雅可比矩阵行列式值。

1.2 积分方法

采用散射场变量既可保持守恒形式不变，又能避免传播入射场的数值计算误差，入射场解析给出。有限体积法能保持整个网格空间的通量守恒，对式(4)在网格单元作时、空积分，得到数值离散形式：

有限体积法的关键是精确模拟原变量Q在单元分界面的状态变量，以得到相应精确的分界面流通量F，本文F采用STEGER-WARMING分裂：

S、S-为相似矩阵，Λ+、Λ-分别为正负特征值构成的对角矩阵，QL、QR代表分界面处左右状态变量，采用MUSCL格式插值获得，时间推进则采用4步龙格-库塔法。

1.3 入射电磁波波形

入射电磁波选取简谐波(A0是波幅，c为波速)：

1.4 边界条件

完全导电壁处的边界条件：

截断的网格空间外边界相容条件：

1.5 近场到远场的转换

稳定的时变场在频域收敛后，采用格林定理：

R是观察点，R′是虚拟积分面上点矢量，r是观察方向单位矢量，A是虚拟积分面，k是自由空间波数，j是虚数单位。

1.6 雷达散射截面计算

1.7 并行计算

电大尺寸电磁散射伴随大规模三维网格，须用多进程并行计算，包括：(1) 多进程分割和负载平衡；(2) 程序并行化处理，采用MPI接口进行通信。

2 计算结果与分析

2.1 并行化多块FVTD程序pmbRCS3d验证

多块网格FVTD程序mbRCS3d，并行化扩充为pmbRCS3d。首先选择带解析解的金属球(ka=10，其中k为波数，a是半径)为验证目标。

图1是网格块经8个进程分割后表面网格。图2是表面散射场等值线云图。图3是表面诱导电流分布，结果显示前向散射最强。图4是双站RCS数值与精确解析级数解比较，吻合很好，说明整个并行FVTD算法具有很好的精度。

2.2 电大尺寸实际外形RCS计算

针对飞翼外形(类似X-47)的L波段(f=1.2 GHz)电磁散射问题，长度11.64 m，翼展18.9 m，长度对应电尺寸为47个波长，翼展对应76个波长。采用36块网格离散整个计算空间，网格点16 426 816个。网格密度根据电磁波频率选取，要求每个坐标方向上一个波长(250 mm)内至少布15～20个网格点，散射体表面网格加密，远离物面区域适当疏松，以便减少计算量，图5是重新分块后的物面网格边界。在集群上使用16个进程计算，耗时420 h。

该大网格量电磁场并行FVTD计算，得到稳定振荡的简谐场和逐步收敛的频域散射场，从图6和图7散射场云图能看出散射波耦合强烈区域。图8的结果(其中φ为与+X轴所成方位角)表明：前向RCS较大，后向RCS较小在-35 dB附近，有4个较强散射方位(分别对应4个垂直于机身前缘和机翼前缘的方位)。

3 结 论

通过散射电磁场的FVTD解算器和基于MPI的并行编程，一定程度解决了电大尺寸目标高频散射大计算量难题。计算和验证表明：该系统能提供高精度的电磁场和RCS分布，为气动/隐身优化提供精细工具,并有利于向更高频段扩展。

[1]Yee K S.Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media[J].IEEE Trans.Antennas and Propagation,1966,14(3):302-307.

[2]Holland R.Finite-difference solution of Maxwell’s equations in generalized coordinates[J].IEEE Trans.Nuclear Science,2007,30(6):4589-4591.

[3]Fusco M.FDTD algorithm in curvilinear coordinates[J].IEEE Transactions on Antennas and Propagation,1990,38(1):76-89.

[4]Song J M,Chew W C.Fast multipole method solution of combined field integral equation[J].IEEE Trans.on Antennas and Propagation,1994,7(16):760-765.

[5]Nachman A.A brief perspective on computational electromagnetic[J].Journal of Computational Physics,1996,125:237-239.

[6]V elamparambil S,Chew W C,Song J M.10 million unknowns:is it that big?[J].IEEE Antennas and Propagation Magazine,2003,45(2):43-58.

[7]Pan X M,Sheng X Q.A sophisticated parallel MLFMA for scattering by extremely large targets[J].IEEE Antennas and Propagation Magazine,2008,50(3):129-138.

[8]Shankar V,Hill W,Mohammadian A H.A CFD-based finite-volume procedure for computational electromagnetics inter-disiplinary applications of CFD methods[R].A1AA CP 89-1987.

[9]Shankar V,Hill W,Mohammadian A H.A time-domain differential solver for electromagnetic scattering problems[J].Proc.IEEE,1989,77(5):709-721.

[10]Ahuja V,Long L N.A parallel finite-volume Runge-Kutta algorithm for electromagnetic scattering[J].Journal of Computational Physics,1997,137:299-320.

[11]Kabakian A V.A three-dimensional spectral collocation time- domain solver for electromagnetic wave scattering[R].AIAA 98-0980.

[12]Shang J S,Gaitonded.Characteristic-based time-dependent Maxwell equations solver on a general curvilinear frame[R].AIAA 93-3178,1993.

[13]许勇,乐嘉陵.基于CFD的电磁散射对数值模拟[J].空气动力学学报,2004,22(2):185-189.

[14]Shang J S.Time-domain electromagnetic scattering simulations on multicomputers[J].Journal of Computational Physics,1996,128:381-390.

[15]Shang J S,Gaitonde D.Characteristic-based time-dependent Maxwell equations solvers on a general curvilinear frame[J].AIAA Journal,1995,33(3):491-498.

ParellelFVTDcomputationforelectromagneticscatteringofelectricallylargeobjects

XU Yong*,HUANG Yong,YU Yonggang

(ComputationalAerodynamicsInstituteofChinaAerodynamicsResearchandDevelopmentCenter,Mianyang621000,China)

Numerical simulation of real military objects′ electromagnetic scattering at high frequency bands is ordinary corresponding to the electromagnetic field calculation for electrically large objects.This simulation can be computationally expensive with large storage and calculation works when high precision numerical algorithm is used.To ensure high computational precision,we programed a finite volume time domain (FVTD) solver pmbRCS3d based on body-fitted multi-blocks grids,which computes the Maxwell equations of electromagnetism directly.In order to handle the massive calculation work due to huge number of grids,a parallel MPI programing technique was used together with careful design of the grid partition,load balance,data communication,and parallelization.The pmbRCS3d software was applied successfully to calculate the electromagnetic scattering field and the bi-station RCS of an L-band electrically large flying-wing shape.The results show that the developed pmbRCS3d software has stable and robust characteristic and can be used for the RCS calculation of higher frequency band objects.

Maxwell equations; finite volume time domain method; parallel computing; L frequency band; radar cross section

0258-1825(2017)06-0797-04

O441.4

A

10.7638/kqdlxxb-2015.0071

2015-06-03;

2015-07-24

装备预先研究项目(51313010303)

许勇*(1971-),男，四川彭州人，副研究员，研究方向：计算电磁学.E-mail：stephen000@sina.com

许勇,黄勇,余永刚.电大尺寸目标电磁散射的并行FVTD计算[J].空气动力学学报,2017,35(6):797-800.

10.7638/kqdlxxb-2015.0071 XU Y,HUANG Y,YU Y G.Parellel FVTD computation for electromagnetic scattering of electrically large objects[J].Acta Aerodynamica Sinica,2017,35(6):797-800.