高超声速飞行体亚密湍流尾迹RCS特性的相似规律研究

再入或临近空间高超声速飞行器在大气层中飞行时,会形成等离子体绕流和尾迹,将严重影响飞行器的电磁散射特性,从而使得临近空间高超声速飞行器表现出与大气层内常规飞行器和大气层外轨道飞行器在电磁散射特性上的明显不同,对飞行器的探测、跟踪带来新的挑战。弹道靶是高超声速飞行器尾迹电磁散射特性研究的重要设备,为了由弹道靶试验结果外推出真实飞行条件下亚密湍流尾迹的RCS,本文从高超声速流场模拟的双缩尺率和亚密湍流尾迹RCS模拟的Born近似出发,推导出真实飞行条件下和地面弹道靶试验之间亚密湍流尾迹RCS模拟的一种相似规律,推导结果显示当D1/D2=η,V∞1=V∞2,ρ∞1/ρ∞2=1/η,ω1/ω2=λ2/λ1=1/η时,雷达散射截面服从σ1/σ2=η4的相似规律。推导出的相似规律可以为超高速飞行器亚密湍流尾迹RCS研究提供技术支持。     

一种快速模拟振荡叶栅非定常流的数值方法

为了提高叶轮机内周期性非定常流数值计算效率,发展了一种计算振荡叶栅非线性非定常流的谐波平衡方法,将周期性非定常流动参数采用傅里叶级数表示,从而将传统的定常流动数值计算方法应用于非定常流问题求解中,大大提高了非定常流计算效率。计算结果表明,发展的谐波平衡方法与线性方法、非线性时域法计算结果吻合较好,且其计算效率比时域法高1个数量级。另外,系统研究了谐波阶次对计算效率与计算精度的影响,认为在气流未出现严重分离时,应采用低阶谐波平衡方法,以减少计算耗时。此外,从压气机振荡叶栅非定常流场计算结果可以看出,叶片吸力面诱导涡的生成频率与振荡频率并不一致,涡生成频率要高于振荡频率。     

旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

基于北斗的车载接收机航向测量技术研究

研究一种BD2双天线的双频基线测量算法,实现了在车载动态条件下利用BD2卫星进行载体航向实时测量,满足军用车辆、无人机、船舶舰艇等载体高精度、快速实时、连续稳定长时间工作、低成本的要求,对算法实现进行了试验验证,取得了理想的成效,具有较强的工程实用价值。     

基于人工免疫算法的载人登月任务地月转移双脉冲中止策略研究彭 坤,果琳丽,向开恒,王 平,杨 雷

载人登月任务中,任务中止策略设计是确保航天员安全返回的重要基础。首先结合"星座"计划飞行方案分析了载人登月任务各飞行阶段的中止策略;其次针对地月转移巡航段进行了双脉冲中止策略设计,以速度增量数值、方位角以及变轨时间间隔为控制变量,加入轨道同向、近地点高度、偏心率以及飞行时间约束,提出双脉冲变轨计算流程;最后采用人工免疫算法对该问题进行了求解和优化。仿真算例表明,双脉冲中止策略存在多组解,其全局分布特性为:飞行时间越短速度增量需求越大;飞行时间相近时,大偏心率中止轨道对应的速度增量小;故障点离地月加速点越近,所需速度增量越小。同时也验证了人工免疫算法求解双脉冲中止策略问题的有效性。     

弯曲网格上的间断有限元湍流数值解法研究

采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。     

多尾鳍仿生航行器推进性能的三维数值研究

采用格子波尔兹曼方法及浸没边界法,发展了一套适用于自由游动的三维数值研究程序。进而以多尾鳍推进仿生航行器为原型,讨论了尾鳍数目、形状、材料刚度等构型参数及振幅、频率等尾鳍摆动参数对推进性能的影响。结果表明对称布置、反对称摆动的双尾鳍能够明显消除侧向力而避免航行器主体的横向晃动,且推力大于两个单独尾鳍的简单加和;在中等Reynolds数下,具有波动性质的柔性尾鳍不论是推进速度还是推进效率均优于刚性尾鳍;当尾鳍单纯摆动推进时,具有完整鳍面的半椭圆形尾鳍的推进性能优于后部有缺口的深叉形尾鳍。     

目标函数误差估算及网格自适应处理

发展了一种基于目标函数误差估算的网格自适应准则,进而通过网格自适应处理提高目标函数计算的效率和准确性。首先描述了目标函数的误差估算及修正的方法,该方法通过伴随方程将原方程的残值误差与目标函数联系起来。然后,基于误差估算建立网格自适应准则,以减少目标函数修正后的剩余误差,提高目标函数计算结果的准确性,并将方法进一步发展至多目标问题。最后将该准则应用于Euler方程控制的NACA0012翼型无粘可压流动的网格自适应模拟。数值计算成功地捕获了与升力、阻力和力矩等目标函数相关的特征流动区域,计算出符合指定精度要求的目标函数值,验证了本文所发展的方法。     

悬停状态倾转旋翼/机翼干扰流场及气动力的CFD计算

基于一套高效通用的多层运动嵌套网格技术,建立了适合悬停状态倾转旋翼机旋翼/机翼气动干扰特性分析的高效混合计算流体力学（CFD）方法。在倾转旋翼/机翼贴体网格区采用可压雷诺平均Navier-Stokes（RANS）方程作为主控方程,过渡/背景网格区选用Euler方程,湍流模型选用Spalart-Allmaras模型。时间推进上采用双时间推进格式进行非定常求解,并在方法中运用了SPMD（Single Program Multiple Data）模式的并行加速技术。在此基础上,首先分别采用UH-60A直升机旋翼及XV-15倾转旋翼机旋翼作为数值算例,验证了CFD方法的有效性。然后着重对倾转旋翼/机翼的非定常干扰流场及气动力分布特征进行了数值研究,模拟得到与实际情况相符的“喷泉效应”干扰现象。计算结果表明,干扰作用使得倾转旋翼相对于孤立旋翼拉力数值减小了3%,但总的拉力系数损失达到了17%,证明悬停状态下气动干扰对飞行器气动性能有重要影响。     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。