三角翼大迎角风洞试验支架干扰数值模拟研究

现代战争要求战斗机能够在大迎角（AOA）状态下进行过失速飞行,对飞机大迎角绕流流场的研究主要的方法有风洞试验和数值模拟。在大迎角风洞试验中,常用的是尾支撑方法,支架的存在会对模型的试验结果产生一定的影响,本文通过数值模拟来对这个影响进行研究。以开源计算流体力学软件OpenFOAM 2.3为平台,采用PIMPLE算法求解Navier-Stokes（N-S）方程,PIMPLE算法是SIMPLE（Semi-Implicit Method for Pressure-linked Equations）算法和PISO（Pressure Implicit with Splitting of Operator）算法的结合体;采用基于有限体积的空间离散方法和空间二阶精度的线性插值方法,时间离散采用后向差分方法,湍流模型采用SA-DDES（Spalart-Allmaras-Delayed Detached Eddy Simulation）模型。为了验证方法的可靠性,首先对0°、10°、30°、50°、70°以及90°迎角下的有支架三角翼绕流流场进行计算,并将计算结果与试验结果进行对比,两者吻合较好。在此基础上,数值模拟了无支架的三角翼绕流流场,对比有/无支架情况下数值模拟结果,得到支架对三角翼绕流流场、背风面压强分布和气动力的影响。计算结果表明：大迎角情况下,有支架与无支架时相比,支架的存在会影响三角翼附近的流场（但是不会改变涡系等流动结构）、改变翼表面压强分布,从而导致三角翼的法向力系数和俯仰力矩系数发生明显变化。     

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

多段翼构型的积冰数值模拟

基于非结构网格提出一种针对多段翼构型的积冰数值模拟方法.通过求解气体欧拉方程得到多段翼的绕流流场,基于空气绕流流场的气流速度分布,采用迭代法对水滴控制方程进行数值求解,得到多段翼各翼面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对某典型三段翼霜冰和明冰的结冰情况进行了预测研究,根据多段翼的流场特点提出一种明冰冻结模型的求解思路,并分析了积冰对多段翼性能的影响.      

基于结构网格的栅格翼绕流数值模拟

针对栅格翼这种复杂的构造形式,采取有效的多块结构网格生成策略,对栅格翼导弹生成高质量贴体的结构网格。在此基础上,通过求解N-S方程,对栅格翼导弹亚、跨、超声速绕流流场进行了数值模拟。计算结果与风洞试验结果吻合较好,能够准确地捕捉到栅格框内复杂流场的流动特性,从而结构网格生成策略与计算方法得到验证,为栅格翼导弹绕流流场的精确数值模拟提供了技术途径。     

固体火箭发动机喷管燃气舵铰链力矩参变分析

以固体火箭发动机喷管燃气舵绕流流场为研究对象,采用数值模拟技术,通过离散求解流动方程及相关湍流模型方程,对包含舵基、舵片和发动机喷管在内的复杂流场进行了变参数分析,获得并揭示了舵基、舵片附近的流动现象和流动特点.验证了所引入的数值模拟方法的有效性,分析结果对舵基、舵片的结构安装以及控制机构的设计具有借鉴意义.     

简化的积冰热力学模型及其应用

依据飞行器表面液态水凝结过程的热力学理论分析,提出了一种适用于不同环境温度的翼型结冰数值模拟方法.采用中心型有限体积方法求解N-S方程来获得翼型的绕流流场,结合给出的结冰数值模拟方法对两种典型温度下的翼型积冰情况进行了模拟.模拟所得的积冰冰形与文献中的实验数据吻合良好,表明该方法的可行性,并进一步给出了翼型表面的冻结系数分布.     

三角翼匀速上仰跨声速非定常流场数值模拟

用有限差分法求解Ｎ－Ｓ方程，数值模拟了细长三角翼作匀速上仰运动时的距声速绕流流场。给出了非定常流场结构，涡破裂特性和非定常升力，阻力特性，研究了上仰运动对动态失速的影响。     

多段翼型低速流场与气动特性数值计算研究

在结构化网格求解RANS方程的基础上,将嵌套网格技术应用在低速绕流气动特性计算中,并以NLR-7301两段翼型和某运输机5段翼型在起飞构型下对低速绕流流场进行数值模拟,将气动特性及流场特性的数值模拟数据与对接网格计算以及风洞实验数据进行了比较.结果表明,该气动特性计算方法有着较高的精度和可靠性.     

基于阻力气动特性计算的飞艇艇身外形研究

采用SIMPLE算法、S-A湍流模型求解了三维不可压雷诺平均N-S方程,数值模拟了绕飞艇的低速不可压粘性流动,计算了不同艇身外形的阻力特性.本文方法计算结果与实验结果或其他文献的计算结果符合良好,可以用于飞艇的气动特性计算.通过对6个不同外形的相同体积飞艇阻力特性计算,得到了最佳艇身外形,表明本文方法可以用于艇身外形的选型设计.     

基于气动数值模拟的翼型反设计方法

将计算流体动力学（CFD）与反设计技术相结合,通过数值求解欧拉方程,对翼型绕流流场作出数值模拟,再用几何和流动控制方程,反复迭代求得满足给定流场的翼型。以NACA0012为初始翼型,RAE-2822为目标翼型,选取两种工况,都取得了满意的结果。     

自适应变搜索域遗传算法及其在发动机模型中的应用

针对发动机数学模型中非线性方程插值解法的不足,提出非线性方程自适应变搜索域遗传算法解法。论文详细分析了方程求解转化为遗传算法优化的数学描述、依据方程解临域特征的自适应变搜索域机理及算法的具体实现技术。实例计算与发动机仿真结果说明:所设计的变搜索域遗传算法可作为通用的非线性方程解法,相比较于简单遗传算法,能成倍提高计算效率;替代插值解法应用于发动机模型,可有效提高模型的收敛性。     

三维分块非结构化网格上的可压缩流计算方法

提出了分块非结构化网格生成方法,通过子块划分、子块网格划分、整体网格拓扑连接生成可局部加密的分块非结构化网格,既具备了非结构化网格对复杂区域的处理能力,又避免了分块结构化网格块间信息传递的复杂计算过程.通过将密度修正值引入压力修正方程,发展了非结构化网格上可压缩流动的SIMPLE算法.对孤立翼型及透平叶栅的流动进行了数值计算,压力系数、速度、折转角的数值结果与已有的实验值吻合较好,验证了数值方法的精度.      

燃烧室贫油熄火极限数值预测

在任意曲线坐标系下分别对双级轴向涡流器矩形燃烧室及涡流杯回流环形燃烧室三维两相反应流场进行数值模拟,采用RNG(renormalization group)k-ε,SOM-EBU(second order moment-eddy break up model)及六通量辐射等模型考虑紊流流动、化学反应及辐射传热对燃烧过程的影响;采用SI MPLE(semi-implicit method for pressure linked equation)算法求解流场,气液两相耦合采用颗粒源项法处理.采用提出的燃油稳态逐次逼近法对贫油熄火极限进行数值预测,所得计算值与实验数据较为相符,表明此数值预测方法可靠,可用于工程实践.      

敷设黏弹性阻尼层的复合材料结构非平稳随机振动分析

提出了对敷设黏弹性阻尼层的复合材料层合结构非平稳随机振动分析的高效精确算法.考虑了黏弹性阻尼材料对激振频率的依赖性及Adams提出的复合材料层合结构的阻尼模型.为解决依赖频率变化的非经典阻尼的复杂结构,结合实模态降阶直接解法对虚拟激励法做了相应的发展,并采用精细积分法高效精确地求解了随机振动方程.以国产有限元程序系统DDJ为平台开发了相应的程序,并对某飞机的水平尾翼的复合材料安定面结构进行模拟研究,从复杂的计算及合理的结论可以看出,本文中所提出的方法对于这类十分复杂复合材料结构在航空领域的非平稳随机振动分析十分有效.     

稀薄流到连续流的气体运动论统一数值算法初步研究

从非线性模型Ｂｏｌｔｚｍａｎｎ方程出发,引入简化速度分布函数、使用离散速度坐标法对速度空间进行离散、降维,去掉分布函数对速度分量的连续依赖性;采用时间分裂法,将简化速度分布函数松驰变化方程分解为源项碰撞变化方程、对流运动方程,进行耦合计算,应用ＮＮＤ耗散差分方法直接模拟气体分子速度分布函数;发展离散速度数值积分法,通过宏观取矩获取物理空间各点的流动参数,从而建立一套能有效模拟各流域气动问题的简化的     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

基于Chebyshev正交分解的曲线运动轨迹SAR的Chirp Scaling算法

针对具有三维速度和加速度的曲线运动轨迹合成孔径雷达（SAR），传统的斜距模型无法精确描述其运动特性，曲线历程增加了距离走动现象和方位向时间的高次项，使二维耦合现象更为复杂。本文提出了一种考虑载体平台三维速度和加速度的Chirp Scaling算法以解决曲线运动轨迹SAR成像问题。首先根据运动方程建立斜距表达式，然后对其进行Chebyshev近似，并构造其等效双曲方程形式的斜距模型，推导了具有空变性的距离徙动函数，Chirp Scaling因子以及适用于曲线轨迹的Chirp Scaling成像算法。仿真结果证实了此扩展的等效斜距模型和Chirp Scaling算法在大合成孔径时间下的有效性，并给出了三维加速度的边界值。     

基于锅炉炉膛三维温度场检测算法的研究

鉴于大型锅炉内燃烧火焰变化的多样性,为了实时掌握炉内火焰燃烧状况,提出了一种改进的检测炉内三维温度场的算法.该算法通过对炉膛火焰二维累加辐射能的图像处理,并利用辐射传递方程的逆求解,实现了对炉内火焰三维温度场的确定.仿真研究和实验结果表明此算法是可行的.     

复杂前机身—进气道组合体的流场计算

本文提出了一种能够对复杂前机身-进气道组合体进行流场计算与气动分析的方法。通过求解Euler方程来模拟前机身-进气道组合体在亚、跨、超声速,不同攻角和发动机流量条件下的各类流场。