翼型外形高气动效率／低可探测性的优化.pdfVIP

翼型外形高气动效率／低可探测性的优化。 朱自强ZubairIslam朱一琨李海明 (北京航空航天大学流体力学研究所，北京，100083) 摘要采用VanLeer矢通量分裂格式求解EuMr方程的方法计算了绕翼型的气动特性；采用矢通量分裂方法计算了绕翼型的时 域电磁散射场特性及雷达散射截面积(RCS)；采用一种简单而有效的数值优化方法对流场解和电磁场解进行了翼型外形高 气动效率，低可探测性的优化计算．算例结果表明，本方法提供了一种对嚣型既可作气动优化设计亦可进行多学科综合设计 的有效工具． 0引言 飞行器的隐形性能已成为飞行器不可缺少的一种重要性能．现代战争中雷达是探测飞行器最有效的 方法。因此减弱飞行嚣的雷达反射信号是隐形技术中最关键、晟重要的部份。隐形与气动的综合设计耳标 是妥善处理两者相互矛盾的要求，获得既有低可探测性又有先进飞行性能的外形．通过流场解，电磁场解 的数值优化计算可以实现气动与隐形的综合设计。 计算流体动力学(CFD)的先进计算方法可扩展应用于计算科学的其它领域中实现与其它学科的交 方程)，并采用数值优化方法作了气动／隐形的综台优化计算。 1流场计算 本文在流场计算中基于无粘流的模型求解了Euler方程．无量纲可压缩二维Euler方程的守恒型向量 形式可写作‘1}： 3 盟+一dF+一8G；0(1 dt 毋{ 口口 采用有限体积法将主控方程(1)作空间离散，其中无粘通量向量的离散格式取VanLeer矢通量分裂格式 f2J。用Lu分解方法对方程(】)进行时间积分。为加速收敛，汁算中采用了局部时间步长和多重网格技 术。壁面边界条件采用通常的无穿透条件。无穷远处采用无反射边界条件。用代数网格生成方法形成128 ×32的O型计算网格，共五层。 ■ 2电磁场计算 Maxwell方程是描述电磁场的基本方程，求解时域Maxwell方程完全可描述非定常波的传播和散射现 象【3】。此方程组是双曲型的，与气动计算的Euler方程性质相同，可借鉴CFD中的计算方法和经验。 2．1主控方程和边界条件 在一般曲线座标系中二维Maxwell方程的TM模态可写作pI： Q+疋+G。=o (2) 或 gl+一绞+Bg=0 (3) 将岛和硒无因次化为1。按矩阵A和B的正负特征值可将通量F、G分别分裂为r，F．和G+，G．-并分 别采用后向和前向离散格式，可得最终离散方程f”。对理想导体在物面边界应满足总电场完全反射的条 件。在外边界采用Mur给出的无反射边界条件吼 2．2雷达散射截面积(RCS)的计算 ● 求解最终的离散方程和上述边界条件，可得近区电磁散射场H5和E5。RCS的定义是基于远区散射 ’航空基金资助项目 438 场的特性，故需按等价理论将计算所得的近场值转换成远场值㈣。对远场值作离散付里哀变换，取得频域 值，再按 RCS=磐z华剖2 (4) 求得RCS。式中F为入射电磁场，E邵为散射场远区值：p为散射体中央到观察点的距离． 本文中入射波取为 E‘‘=Eosin[k(x COS∥+ysin∥一cjf)】 s—infl (5) _Vx一：—Ezj ，H，．．：墨：：!!!生 胁c0 心c0 其中Eo是振幅， k=2万／五是波数，／3是入射角．