求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

随机集的概率假设密度粒子滤波

多目标跟踪问题中,当目标数已知时,可以用概率数据互联(PDA)或联合概率数据互联(JPDA)算法。而当目标数未知或随时间变化时,需要对不同目标数的跟踪进行比较。可以把目标集看作随机集进行讨论,目标数N是随机变量。随机集的跟踪通过有限集统计(FISST)理论来完成。文中讨论了用粒子滤波实现跟踪随机集的方法。实验表明,在杂波环境下,粒子滤波可以稳健跟踪目标状态和目标数。     

直升机桨叶刚柔耦合特性及计算方法分析

为了研究桨叶刚柔耦合特性,采用中等变形梁模型处理桨叶弹性变形,分别采用小转角和有限转角两种方法处理桨叶绕铰的刚性运动。以桨叶扬起下坠碰撞问题为研究对象。计算分析表明:(1)弹性桨叶存在较明显的刚性运动和弹性变形耦合;(2)挥舞角较小时,两种处理方法计算结果一致,均与试验数据吻合较好;(3)挥舞角较大时,刚柔耦合增强,两种处理方法的计算值差别较为明显;(4)桨叶弯曲刚度对桨叶刚柔耦合特性有较明显影响,刚度越小,耦合越强;(5)桨叶如果比较柔软,小转角处理方法计算值相对有限转角处理方法,相位滞后、幅值增加均较为明显,计算桨叶动响应时需引起注意。      

等离子体涡电磁散射特性及隐身性能

 针对非均匀等离子体在飞行器隐身中的应用,采用分段线性电流密度递归卷积时域有限差分（PLJERC-FDTD）方法计算等离子体涡及涡串电磁散射特性,分析等离子体涡对飞行器隐身性能影响。计算表明,等离子体涡在很大频率区间对电磁波吸收效果显著,RCS降低很大,具有明显的隐身效果。等离子体涡表现出一定规律性的极化特性,对L,S和C波段电磁波具有不同的吸收、反射特性。     

超燃冲压发动机燃烧室的准一维计算与分析

在考虑有限速率化学反应的准一维Euler方程基础上，通过增加截面面积变化、壁面摩擦和添质的源项，发展了适用于超燃燃烧室性能分析的准一维计算方法。依次以中国空气动力研究与发展中心（CARDC）和日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为验证算例，分别采用传统的一步反应模型和发展的有限速率反应模型，模拟了燃烧室流场，并基于NAL燃烧室，计算分析了不同当量比和进口压强对燃烧室流动特性的影响。结果表明：两种方法都能得到与实验数据吻合良好的结果；和一步反应模型相比，有限速率反应模型不仅可以更细致地捕捉流场细节，而且能够初步分析化学非平衡效应的影响；对于NAL燃烧室，当量比≥0.6时，压强随当量比的升高而增大，当达到1.0时，反压已推进隔离段，且推进速度随当量比增大而增加；进口压强不大于110.444kPa时，反压随进口压强增大而升高，且当反压不小于82.833kPa时，反压被隔离在等直段燃烧室入口处；过小的当量比和过大的进口压强均会导致燃烧室出口马赫数严重下降，甚至出现亚声速出流状态。     

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

感性耦合夹层等离子体隐身天线罩电磁散射分析

设计了一种石英夹层感性耦合等离子体（ICP）隐身天线罩模型,采用有限元与Z变换时域有限差分（ZT-FDTD）联合仿真的方法,建立了ICP放电的流体模型,得到不同气压及功率下与电磁散射相关的电子密度空间分布,在此基础上建立Z变换时域有限差分法模型,对石英夹层等离子体隐身天线罩的宽频段后向散射进行计算,同时利用微波干涉法及XFDTD软件对算法及程序进行验证。结果表明：感性耦合等离子体可产生较高密度等离子体,能有效实现雷达散射截面（RCS）的减缩,在气压为2 Pa时,碰撞衰减较弱,等离子体密度分布较均匀,衰减带宽集中在等离子体振荡频率附近,功率增加会使衰减峰值向高频方向移动,气压为20 Pa时,碰撞衰减增强,且等离子体密度分布有较大梯度,衰减带宽有效增加,同时RCS曲线的波动特性加强。     

带空间协同的多导弹时间协同制导律

针对多导弹在平面内从期望的弹目视线（LOS）相对方向同时击中固定目标问题,提出了一种带空间协同的多导弹时间协同制导律。基于平面内的导弹-目标相对运动方程,建立了带空间协同的多导弹时间协同制导模型;基于多智能体协同控制理论,在视线方向设计了分布式时间协同制导律,可保证所有导弹的打击时刻在有限时间内达到一致,在视线法向方向设计了分布式空间协同制导律,可保证所有导弹的相对视线角在有限时间内收敛到期望值;最后,通过仿真验证了所设计的协同制导律可使多导弹从期望的弹目视线相对方向同时击中目标。     

动态时频谱分析、探测和跟踪的随机有限集法

动态出现和消失多分量信号的时频分析问题一直是非平稳信号处理的难点之一。为此,提出了一种分析、探测和跟踪多分量信号的随机有限集法。该算法利用时频变换,如短时傅里叶变换或自适应谱估计法,以及多项式预测模型,将多分量信号的时频分析问题归纳成可利用随机有限集进行多目标追踪的问题。分析表明：借助于提出的初始权重赋值算法,以及谱分量幅度和频率的联合似然函数,就可利用高斯混合概率假设密度滤波器来实现对动态时频谱的分析、探测和跟踪。在仿真实验中,所提算法有效提升了动态时频谱的跟踪精度,其对微弱时频谱分量的探测能力,以及对载频差异的分析能力均优于文献报道的算法。     

用Tam-Ganesan湍流模型研究排气系统红外辐射特性

通过三维流场、介质的质量浓度场与辐射场耦合求解方法,计算了1台有中心锥的轴对称喷管的空腔-喷流组合红外辐射特性.利用专门为自由剪切流动发展的Tam-Ganesan湍流模型,有效提高排气系统红外辐射模拟精度.采用有限体积法(FVM)结合窄带模型(NBM)计算了排气系统在3~5μm波段的辐射传输问题,并考虑大气衰减作用.对比了采用标准 k - ε 模型和Tam-Ganesan模型对排气系统核心区长度的影响,以及由此引起的排气系统红外辐射特性的变化.通过与试验数据对比证明了该结果的准确性,结果表明选择合适的湍流模型对计算排气系统红外辐射特性有重要的意义.