超高速飞行器共形亚网格电磁散射特性分析的仿真方法技术

本发明专利技术公开了一种超高速飞行器共形亚网格电磁散射特性分析的仿真方法。步骤如下：根据飞行器的气动外形和飞行的速度高度进行气动热力学仿真，由仿真信息确定各部分等离子体碰撞频率和等离子体震荡频率；采用四面体对飞行器模型连同等离子鞘套一起进行剖分，得到飞行器模型的结构信息；将流场点信息映射到时域有限差分法计算网格的棱线上，根据棱线上的碰撞频率和震荡频率确定相对介电常数大于6的区域进行亚网格处理；时域有限差分法计算网格的棱线与金属面交界的地方采用共性时域有限差分法进行处理；等离子体区域采用等离子体的迭代公式进行计算，最终确定超高速飞行器的雷达散射截面积。本发明专利技术具有很好的适应性和更高的计算效率。

【技术实现步骤摘要】

本专利技术属于超高速飞行器电磁散射特性
，特别是一种超高速飞行器共形亚网格电磁散射特性分析的仿真方法。
技术介绍
在现代化战争中，武器系统的雷达散射截面积(RCS)是雷达系统对目标“可观测性”的一个重要指标。开展临近空间高超声速目标电磁特性研究，可为有效探测临近空间高超声速飞行器及武器突防提供技术支持。传统的RCS计算软件大多是基于飞行器本体目标开发的，等离子体鞘套与电磁波相互作用机理十分复杂，传统的光学近似分析方法显得较难分析。所以需要采用别的数值方法对临高超声速涂覆隐身目标电磁特性进行计算仿真，从高超声速飞行器流场分布出发，实现临近空间高超声速隐身目标宽频带、宽角度、多站，全极化的电磁特性高效数值分析，为临近空间高超声速飞行器雷达目标特性变化规律的研究以及突防评估与设计提供技术基础。现有的利用时域有限差分法FDTD方法分析高超声速飞行器还存在两个问题：(1)时域有限差分法FDTD方法建模不够灵活，由于高超速飞行器高速飞行产生的等离子体鞘套形状各异，飞行器的气动外形各异，普通时域有限差分法FDTD无法模拟高超速飞行器的气动外形；(2)由于等离子鞘套各部分电磁参数各异，而计算网格的剖分尺寸是按照最大的介电参数确定的，而介电参数较大的部分可能是很小的一部分，这样势必会为了那一小部分相对介电参数较大的部分整体采用细网格，这样就牺牲了整体的计算时间。
技术实现思路
本专利技术的目的在于提供一种适应性好、计算效率高的超高速飞行器共形亚网格电磁散射特性分析的仿真方法。实现本专利技术目的的技术解决方案为：一种超高速飞行器共形亚网格电磁散射特性分析的仿真方法，步骤如下：步骤1，根据飞行器的气动外形和飞行的速度高度进行气动热力学仿真，得到目标周围的温度、压强、离子浓度，由仿真信息确定各部分等离子体碰撞频率υ和等离子体震荡频率ωp；步骤2，采用四面体对飞行器模型连同等离子鞘套一起进行剖分，得到飞行器模型的结构信息，所述飞行器模型的结构信息包括四面体的个数和每个四面体上点的坐标；步骤3，将流场点信息映射到时域有限差分法FDTD计算网格的棱线上，根据棱线上的碰撞频率和震荡频率确定相对介电常数εr大于6的区域进行亚网格处理；步骤4，时域有限差分法FDTD计算网格的棱线与金属面交界的地方采用共性时域有限差分法CFDFD进行处理；步骤5，等离子体区域采用等离子体的迭代公式进行计算，最终确定超高速飞行器的雷达散射截面积。本专利技术与现有技术相比，其显著优点为：(1)将物体进行四面体剖分，可以很好的拟合复杂物体的外形，提高计算精度；(2)采用亚网格技术，根据介电常数自适应找到需要进行细网格处理的区域，其余部分仍然采用正常网格，可大大的减小了计算时间；(3)细分的区域可拟合相对介电常数较大部分的形状，具有很好的适应性和更高的计算效率。下面结合附图对本专利技术作进一步详细描述。附图说明图1是本专利技术超高速飞行器与等离子体鞘套剖分示意图。图2是本专利技术超高速飞行器亚网格处理的示意图。图3是本专利技术金属共形示意图。具体实施方式对于超高速飞行器，在高速飞行过程中飞行器周围会形成等离子鞘套，等离子体鞘套与电磁波相互作用机理十分复杂，传统的高频方法较难分析，本专利技术超高速飞行器共形亚网格电磁散射特性分析的仿真方法，在时域有限差分法(FDTD)的基础上对超高速飞行器的电磁特性进行建模分析，加入金属曲面共形技术，提高了计算精度，具体步骤如下：步骤1，根据飞行器的气动外形和飞行的速度高度利用FD-FASTRAN软件进行气动热力学仿真，得到所需要的流场信息即目标周围的温度、压强、离子浓度，由仿真信息确定各部分等离子体碰撞频率υ和等离子体震荡频率ωp。FD-FASTRAN软件特别适合应用于航空、航天工业中，尤其在高超声速CFD计算领域以其高效、稳健、合理的计算结果而获得业内人士的广泛认可。步骤2，采用四面体对飞行器模型连同等离子鞘套一起利用商用软件进行剖分，如图1所示，得到飞行器模型的结构信息，所述飞行器模型的结构信息包括四面体的个数和每个四面体上点的坐标。步骤3，将流场点信息映射到时域有限差分法FDTD计算网格的棱线上，由映射找到这些棱线中相对介电参数较大的部分，根据棱线上的碰撞频率和震荡频率确定相对介电常数εr大于6的区域进行亚网格处理；如图2所示的是一个等离子体鞘套示意图，颜色较深的区域(两条细带之间)，相对介电参数较大，示意图中采用了细网格处理，用细网格拟合这部分形状，如图2中的阴影部分所示。所述将流场点信息映射到时域有限差分法FDTD计算网格的棱线上，即取棱线的三等分点，并依次确定离每个等分点最近的流场点的等离子体碰撞频率υ和等离子体震荡频率ωp，将确定的三个等离子体碰撞频率υ取平均得到该条棱线的等离子体碰撞频率将三个等离子体震荡频率ωp取平均得到该条棱线的等离子体震荡频率根据棱线上的碰撞频率和震荡频率结合公式(1)确定相对介电常数εr大于6的区域，ω是入射波的角频率，则：ϵr=1-ω‾p2ω2+υ‾2-jυ‾ωω‾p2ω2+υ‾2---(1)]]>相对介电常数εr大于6的区域进行亚网格处理。步骤4，时域有限差分法FDTD计算网格的棱线与金属面交界的地方采用共性时域有限差分法CFDFD进行处理。图3为理想导体的共形网格在xoy平面内的截面示意图，阴影部分为理想导体区域，空白处为自由空间区域，普通的磁场迭代公式如(2)所示：Hz(i,j,k)n+1/2=Hz(i,j,k)n-1/2+Δtu·[Ex(i,j+1,k)n-Ex(i,j,k)nΔy-Ey(i+1,j,k)n-Ey(i,j,k)nΔx]---(2)]]>在计算网格被飞行器外表面所截时候采用共形处理，将磁场迭代公式(2)进行修正如(3)所示，而电场迭代公式保持不变Hz(i,j,k)n+1/2=Hz(i,j,k)n-1/2+Δtu·Sz(i,j,k)·Ex(i,j+1,k)n·lx(i,j+1,k)-Ex(i,j,k)n·lx(i,j,k)-Ey(i+1,j,k)n·ly(i+1,j,k)+Ey(i,j,k)n·ly(i,j,k本文档来自技高网...

【技术保护点】
一种超高速飞行器共形亚网格电磁散射特性分析的仿真方法，其特征在于，步骤如下：步骤1，根据飞行器的气动外形和飞行的速度高度进行气动热力学仿真，得到目标周围的温度、压强、离子浓度，由仿真信息确定各部分等离子体碰撞频率υ和等离子体震荡频率ωp；步骤2，采用四面体对飞行器模型连同等离子鞘套一起进行剖分，得到飞行器模型的结构信息，所述飞行器模型的结构信息包括四面体的个数和每个四面体上点的坐标；步骤3，将流场点信息映射到时域有限差分法FDTD计算网格的棱线上，根据棱线上的碰撞频率和震荡频率确定相对介电常数εr大于6的区域进行亚网格处理；步骤4，时域有限差分法FDTD计算网格的棱线与金属面交界的地方采用共性时域有限差分法CFDFD进行处理；步骤5，等离子体区域采用等离子体的迭代公式进行计算，最终确定超高速飞行器的雷达散射截面积。

【技术特征摘要】
1.一种超高速飞行器共形亚网格电磁散射特性分析的仿真方法，其特征在于，步
骤如下：
步骤1，根据飞行器的气动外形和飞行的速度高度进行气动热力学仿真，得到目标
周围的温度、压强、离子浓度，由仿真信息确定各部分等离子体碰撞频率υ和等离子体
震荡频率ωp；
步骤2，采用四面体对飞行器模型连同等离子鞘套一起进行剖分，得到飞行器模型
的结构信息，所述飞行器模型的结构信息包括四面体的个数和每个四面体上点的坐标；
步骤3，将流场点信息映射到时域有限差分法FDTD计算网格的棱线上，根据棱线
上的碰撞频率和震荡频率确定相对介电常数εr大于6的区域进行亚网格处理；
步骤4，时域有限差分法FDTD计算网格的棱线与金属面交界的地方采用共性时域
有限差分法CFDFD进行处理；
步骤5，等离子体区域采用等离子体的迭代公式进行计算，最终确定超高速飞行器
的雷达散射截面积。
2.根据权利要求1所述的超高速飞行器共形亚网格电磁散射特性分析的仿真方法，
其特征在于，步骤3所述将流场点信息映射到时域有限差分法FDTD计算网格的棱线上，
即取棱线的三等分点，并依次确定离每个等分点最近的流场点的等离子体碰撞频率υ和
等离子体震荡频率ωp，将确定的三个等离子体碰撞频率υ取平均得到该条棱线的等离子
体碰撞频率将三个等离子体震荡频率ωp取平均得到该条棱线的等离子体震荡频率
根据棱线上的碰撞频率和震荡频率结合公式(1)确定相对介电常数εr大于6
的区域，ω是入射波的角频率，则：
ϵr=1-ω‾p2ω2+υ‾2-jυ‾ωω‾p2ω2+υ‾2---(1)]]>相对介电常数εr大于6的区域进行亚网格处理。
3.根据权利要求1所述的超高速飞行器共形亚网格电磁散射特性分析的仿真方法，
其特征在于，步骤4所述共性时域有限差分法CFDFD具体如下：
普通的磁场迭代公式如(2)所示：
Hz(i,j,k)n+1/2=Hz(i,j,k)n-1/2...

【专利技术属性】
技术研发人员：陈如山，丁大志，樊振宏，张超，黄静，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏;32